Dependence of earthquake recurrence times and independence of magnitudes on seismicity history
Introduction
Although seismology has intensively analyzed the characteristics of single earthquakes with great success, a clear description of the properties of seismicity (i.e., the structures formed by all earthquakes in space, time, magnitude, etc.) is still lacking. It is true that a number of empirical laws have been established since the end of the XIX century, but these laws are not enough to provide a coherent phenomenology for the overall picture of earthquake occurrence (Mulargia and Geller, 2003). In general, the approach to this phenomenon has to be statistical, as the Gutenberg–Richter law for the number of occurrences as a function of magnitude exemplifies. Moreover, for some studies it is necessary to assume that an earthquake can be reduced to a point event in space and time; despite this drastic simplification, the resulting point process for seismicity retains a high degree of complexity. In this context, the problem turns out to be very similar to those usually studied in statistical physics. Remarkably, in the last years the so-called ETAS model (which combines the Gutenberg–Richter law, the Omori law, and the law of productivity of aftershocks) has been widely studied (Ogata, 1988, Ogata, 1999, Helmstetter and Sornette, 2002); however, this model is still a first order approximation to real seismicity, useful mainly as a null-hypothesis model.
In contrast to the issues of the size of earthquakes or the temporal variations of the rate of occurrence, the timing of earthquakes (despite its potential interest for forecasting and risk assessment) has been the subject of relatively few empirical statistical studies, with the additional drawback that no unifying law for the time occurrence has emerged, at variance with the general acceptation of the Gutenberg–Richter law and the Omori law (Gardner and Knopoff, 1974, Udías and Rice, 1975, Koyama et al., 1995, Wang and Kuo, 1998, Ellsworth et al., 1999, Helmstetter and Sornette, 2004), see also the citations of Smalley et al. (1987) and Sornette and Knopoff (1997). Rather, results turn out to be in contradiction with each other, from the paradigm of regular cycles and characteristic earthquakes to the view of totally random occurrence (Bakun and Lindh, 1985, Sieh et al., 1989, Stein, 1995, Sieh, 1996, Murray and Segall, 2002, Stein, 2002, Kerr, 2004, Weldon et al., 2004, Weldon et al., 2005). In the opposite side of cycle regularity is the view of earthquake occurrence in the form of clusters, which has received more support from catalog analysis (Kagan and Jackson, 1991, Kagan and Jackson, 1995, Kagan, 1997).
A first step towards the solution of the time-occurrence problem and its insertion in a unifying description of seismicity was taken by Bak et al., who related, by means of a unified scaling law, the distribution of earthquake recurrence times with the Gutenberg–Richter law and the fractal distribution of epicenters (Bak et al., 2002, Christensen et al., 2002, Corral, 2003, Corral, 2004b, Corral and Christensen, 2006, Davidsen and Goltz, 2004). One of the goals of this line of research (apart from the applications in forecasting) should be to build a theory of seismicity that clarifies which kind of dynamical process we are dealing with. Simple behaviors such as periodicity, random occurrence, or even chaos seem too simplistic to account for the complexity of seismicity; in contrast, self-organized criticality fits better with this purpose (Bak, 1996, Jensen, 1998, Turcotte, 1999, Sornette, 2000, Hergarten, 2002), although different alternatives have been proposed more recently, see Rundle et al. (2003) and Sornette and Ouillon (2005).
Nevertheless, the recurrence–time distribution introduced by Bak et al. (2002) was defined in a somewhat complicated way through the mixing of recurrence times coming from (more or less) squared equally-sized regions with disparate seismic rates. A simpler, alternative perspective, in which each (arbitrary) region has its own recurrence–time distribution, has been developed (Corral, 2004a, Corral, 2005b). In this case, a universal scaling law involving the Gutenberg–Richter relation describes these distributions for a wide range of magnitudes and sizes of the regions, as long as seismicity shows a stationary behavior.
In this paper we extend this procedure to study conditional probability distributions. These distributions will provide important information about correlations in seismicity, which are fundamental for the existence of the scaling law for recurrence–time distributions. It turns out that recurrence times depend on previous recurrence times and magnitudes, in such a way that small values of the recurrence time and large magnitudes lead to small recurrence times, and vice versa. In contrast, the magnitude is fairly independent on history (although obviously we cannot reject the hypothesis of a dependence weaker than our uncertainty).
Section snippets
Scaling of recurrence-time distributions and renormalization-group transformations
The main subject of our study is the earthquake recurrence time (also called waiting time, interevent time, interoccurrence time, etc.), defined as the time between consecutive events, once a window in space, time and magnitude has been selected. In this way, the i-th recurrence time is defined aswhere ti and ti−1 denote, respectively, the time of occurrence of the i-th and i − 1-th earthquakes in the selected space–time–magnitude window. Note that, following Bak et al., once
Conditional probability densities and correlations
In the previous section we have seen how the fulfillment of a scaling law for the recurrence-time distributions with a non-exponential scaling function is an indication of the existence of important correlations in seismicity. Once a spatial area and a magnitude range have been selected for study, we can consider the process as taking place in the magnitude–time domain (intentionally disregarding the spatial degrees of freedom, which will be the subject of future research, see also Davidsen and
Conclusions
We have argued how a scaling law for earthquake recurrence-time distributions can be understood as an invariance under a transformation analogous to those of the renormalization group (as in the critical phenomena found in condensed matter physics). Such invariance can only exist due to an interplay between the probability distribution and the correlations of the seismicity process. We have seen how the correlations are reflected on conditional probability densities, whose study for stationary
Acknowledgements
The author is grateful to Y. Y. Kagan for his referral of a previous work and other interesting comments. M. Kukito, S. Santitoh, and M. Delduero are also acknowledged for quick inspiration, as well as the Ramón y Cajal program, the National Earthquake Information Center, the Southern California Earthquake Data Center, and the participation at the research projects BFM2003-06033 (MCyT) and 2001-SGR-00186 (DGRGC).
References (63)
Universal local versus unified global scaling laws in the statistics of seismicity
Physica, A
(2004)Statistical aspects of Parkfield earthquake sequence and Parkfield prediction experiment
Tectonophysics
(1997)- et al.
Recurrence intervals between earthquakes strongly depend on history
Physica, A
(2005) How Nature Works: The Science of Self-Organized Criticality
(1996)- et al.
Unified scaling law for earthquakes
Phys. Rev. Lett.
(2002) - et al.
The Parkfield, California, earthquake prediction experiment
Science
(1985) - et al.
Complexity and Criticality
(2005) - et al.
Unified scaling law for earthquakes
Proc. Natl. Acad. Sci. U. S. A.
(2002) Predecessors of the giant 1960 Chile earthquake
Nature
(2005)Local distributions and rate fluctuations in a unified scaling law for earthquakes
Phys. Rev., E
(2003)
Long-term clustering, scaling, and universality in the temporal occurrence of earthquakes
Phys. Rev. Lett.
Comment on “Do earthquakes exhibit self-organized criticality?”
Phys. Rev. Lett.
Mixing of rescaled data and Bayesian inference for earthquake recurrence times
Nonlinear Process. Geophys.
Renormalization-group transformations and correlations of seismicity
Phys. Rev. Lett.
Time-decreasing hazard and increasing time until the next earthquake
Phys. Rev., E
Comment on “Earthquakes Descaled: On Waiting Time Distributions and Scaling Laws"
Phys. Rev. Lett.
The Statistical Analysis of Series of Events
Analysis of Survival Data
An Introduction to the Theory of Point Processes
Are seismic waiting time distributions universal?
Geophys. Res. Lett.
Analysis of the spatial distribution between successive earthquakes
Phys. Rev. Lett.
The longer it has been since the last earthquake, the longer the expected time till the next?
Bull. Seismol. Soc. Am.
A physically-based earthquake recurrence model for estimation of long-term earthquake probabilities
U.S. Geol. Surv. Open-File Rep.
Is the sequence of earthquakes in Southern California, with aftershocks removed Poissonian?
Bull. Seismol. Soc. Am.
Analysis of Temporal Structures of Seismic Events on Different Scale Levels
Applications of Statistical Mechanics
Estimating background activity based on interevent-time distribution
Bull. Seismol. Soc. Am.
Diffusion of epicenters of earthquake aftershocks, Omori's law, and generalized continuous-time random walk models
Phys. Rev., E
Comment on “Power-law time distribution of large earthquakes”
Phys. Rev. Lett.
Comparison of short-term and time-independent earthquake forecast models for Southern California
Bull. Seismol. Soc. Am.
Self-Organized Criticality in Earth Systems
Cited by (82)
Spatiotemporal analysis of earthquake occurrence in synthetic and worldwide data
2022, Chaos, Solitons and FractalsGenerating a seismogenic source zone model for the Pyrenees: A GIS-assisted triclustering approach
2021, Computers and GeosciencesCitation Excerpt :In this case, a memory process refers to the fact that shocks do not happen independently at any time. By contrast, they remember the occurrence (time and magnitude) of the last earthquake (Corral, 2006). The main strength of using zones is that they are efficient because their associated parameters are somehow easily obtained, and their calculations are fast.
A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media
2021, Journal of Contaminant HydrologyCitation Excerpt :The work presented here focuses exclusively on the development and application of such models in the context of geologic media; however, it must be noted that correlated CTRWs, similar in nature to the SMM presented here, play an important role in other branches of physics also (Montero and Masoliver, 2007; Tejedor and Metzler, 2010). As noted by Magdziarz et al. (Magdziarz et al., 2012), correlations arise when dealing with living systems, including bacterial motion (Rogers et al., 2008), the ecology of animal motion (Nathan et al., 2008) and human mobility (Song et al., 2010), as well as in other dynamic systems such as financial markets (Scalas, 2006; Montero and Masoliver, 2017), seismology (Corral, 2006) and chaotic and turbulent flows (Manneville, 2010). It is assumed that the process is stationary such that f(τ(n+1)| τ(n)) = f(τ2| τ1).
Polysemy and brevity versus frequency in language
2019, Computer Speech and LanguageA spatio-temporal binary grid-based clustering model for seismicity analysis
2024, Pattern Analysis and Applications