ETAS: Epidemic-Type Aftershock Sequence
This code was written by Leila Mizrahi with help from Shyam Nandan for the following article:
Leila Mizrahi, Shyam Nandan, Stefan Wiemer 2021;
The Effect of Declustering on the Size Distribution of Mainshocks.
Seismological Research Letters; doi: https://doi.org/10.1785/0220200231
The option for (space-time-)varying completeness magnitude in the parameter inversion is described in:
Leila Mizrahi, Shyam Nandan, Stefan Wiemer 2021;
Embracing Data Incompleteness for Better Earthquake Forecasting. (Section 3.1)
Journal of Geophysical Research: Solid Earth; doi: https://doi.org/10.1029/2021JB022379
To cite the code, please cite the article(s).
For more documentation on the code, see the (electronic supplement of the) articles.
For Probabilistic, Epidemic-Type Aftershock Incomplenteness, see PETAI.
In case of questions or comments, contact me: [email protected].
Contents:
- inversion.py: Expectation Maximization algorithm to estimate ETAS parameters.
- simulation.py: Catalog simulation given ETAS parameters.
- mc_b_est.py: Coupled estimation of beta and Mc. including example estimation.
- invert_etas.py: run example parameter inversion.
- invert_etas_mc_var.py: run example parameter inversion with varying magnitude of completeness.
- simulate_catalog.py: simulate one example catalog.
- magnitudes.npy: magnitude sample for mc/beta estimation.
- california_shape.npy: polygon coordinates of California.
- synthetic_catalog.csv: synthetic catalog to be inverted by invert_etas.py
- synthetic_catalog_mc_var.csv: synthetic catalog to be inverted by invert_etas_mc_var.py
Previous dependencies on certain package versions should now be gone, but just in case, here is my pip freeze:
- geopandas==0.9.0
- numpy==1.19.1
- pandas==1.1.1
- pymap3d==2.4.3
- pynverse==0.1.4.4
- pyproj==3.0.1
- scikit-learn==0.23.2
- scipy==1.5.2
- Shapely==1.7.1