Subcritical and supercritical regimes in epidemic models of earthquake aftershocks
- 18 October 2002
- journal article
- Published by American Geophysical Union (AGU) in Journal of Geophysical Research
- Vol. 107 (B10) , ESE 10-1-ESE 10-21
- https://doi.org/10.1029/2001jb001580
Abstract
We present an analytical solution and numerical tests of the epidemic‐type aftershock (ETAS) model for aftershocks, which describes foreshocks, aftershocks, and main shocks on the same footing. In this model, each earthquake of magnitude m triggers aftershocks with a rate proportional to 10αm. The occurrence rate of direct aftershocks triggered by a single main shock decreases with the time from the main shock according to the “local” modified Omori law K/(t + c)p with p = 1 + θ. Contrary to the usual definition, the ETAS model does not impose an aftershock to have a magnitude smaller than the main shock. Starting with a main shock at time t = 0 that triggers aftershocks according to the local Omori law, which in turn trigger their own aftershocks and so on, we study the seismicity rate of the global aftershock sequence composed of all the secondary and subsequent aftershock sequences. The effective branching parameter n, defined as the mean aftershock number triggered per event, controls the transition between a subcritical regime n < 1 and a supercritical regime n > 1. A characteristic time t*, function of all the ETAS parameters, marks the transition from the early time behavior to the large time behavior. In the subcritical regime, we recover and document the crossover from an Omori exponent 1 − θ for t < t* to 1 + θ for t > t* found previously in the work of Sornette and Sornette for a special case of the ETAS model. In the supercritical regime n > 1 and θ > 0, we find a novel transition from an Omori decay law with exponent 1 − θ for t < t* to an explosive exponential increase of the seismicity rate for t > t*. The case θ < 0 yields an infinite n‐value. In this case, we find another characteristic time τ controlling the crossover from an Omori law with exponent 1 − |θ| for t < τ, similar to the local law, to an exponential increase at large times. These results can rationalize many of the stylized facts reported for aftershock and foreshock sequences, such as (1) the suggestion that a small p‐value may be a precursor of a large earthquake, (2) the relative seismic quiescence sometimes observed before large aftershocks, (3) the positive correlation between b and p values, (4) the observation that great earthquakes are sometimes preceded by a decrease of b‐value, and (5) the acceleration of the seismicity preceding great earthquakes.Keywords
All Related Versions
This publication has 73 references indexed in Scilit:
- A simple and testable model for earthquake clusteringJournal of Geophysical Research, 2001
- Novel mechanism for discrete scale invariance in sandpile modelsZeitschrift für Physik B Condensed Matter, 2000
- Scaling laws in blocks dynamics and dynamic self-organized criticalityPhysics of the Earth and Planetary Interiors, 1997
- Fox function representation of non-debye relaxation processesJournal of Statistical Physics, 1993
- Temporal variations in seismic event rate and b-values from stress corrosion constitutive lawsTectonophysics, 1992
- Temporal variations in seismicity during quasi-static and dynamic rock failureTectonophysics, 1990
- Statistical model for standard seismicity and detection of anomalies by residual analysisTectonophysics, 1989
- Statistical Models for Earthquake Occurrences and Residual Analysis for Point ProcessesJournal of the American Statistical Association, 1988
- The precursory earthquake swarmPhysics of the Earth and Planetary Interiors, 1977
- Multiplicative population chainsProceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 1962