Multiplicative cascades and seismicity in natural time

The debate on the earthquake magnitude correlations: a meta-analysis

Among the most important questions that await an answer in seismology, perhaps one is whether there is a correlation between the magnitudes of two successive seismic events. The answer to this question is considered of fundamental importance given the potential effect in forecasting models, such as Epidemic Type Aftershock Sequence models. After a meta-analysis of 29 papers, we speculate that given the lack of studies carried out with realistic physical models and given the possible bias due to the lack of events recorded in the experimental seismic catalogs, important improvements are necessary on both fronts to be sure to provide a statistically relevant answer.

Natural Time Analysis of Global Seismicity

Natural time analysis enables the introduction of an order parameter for seismicity, which is just the variance of natural time χ, κ1=⟨χ2⟩−⟨χ⟩2. During the last years, there has been significant progress in the natural time analysis of seismicity. Milestones in this progress are the identification of clearly distiguishable minima of the fluctuations of the order parameter κ1 of seismicity both in the regional and global scale, the emergence of an interrelation between the time correlations of the earthquake (EQ) magnitude time series and these minima, and the introduction by Turcotte, Rundle and coworkers of EQ nowcasting. Here, we apply all these recent advances in the global seismicity by employing the Global Centroid Moment Tensor (GCMT) catalog. We show that the combination of the above three milestones may provide useful precursory information for the time of occurrence and epicenter location of strong EQs with M≥8.5 in GCMT. This can be achieved with high statistical significance (p-values of the order of 10−5), while the epicentral areas lie within a region covering only 4% of that investigated.

Order Parameter and Entropy of Seismicity in Natural Time before Major Earthquakes: Recent Results

A lot of work in geosciences has been completed during the last decade on the analysis in the new concept of time, termed natural time, introduced in 2001. The main advances are presented, including, among others, the following: First, the direct experimental verification of the interconnection between a Seismic Electric Signals (SES) activity and seismicity, i.e., the order parameter fluctuations of seismicity exhibit a clearly detectable minimum when an SES activity starts. These two phenomena are also linked closely in space. Second, the identification of the epicentral area and the occurrence time of an impending major earthquake (EQ) by means of the order parameter of seismicity and the entropy change of seismicity under time reversal as well as the extrema of their fluctuations. An indicative example is the M9 Tohoku EQ in Japan on 11 March 2011. Third, to answer the crucial question—when a magnitude 7 class EQ occurs—whether it is a foreshock or a mainshock. This can be answered by means of the key quantities already mentioned, i.e., the order parameter of seismicity and the entropy change of seismicity under time reversal along with their fluctuations. The explanation of the experimental findings identified before major EQs is given in a unified way on the basis of a physical model already proposed in the 1980s.

The Evolution of Preseismic Patterns Related to the Central Crete (Mw6.0) Strong Earthquake on 27 September 2021 Revealed by Multiresolution Wavelets and Natural Time Analysis

On 27 September 2021, a shallow earthquake with focal depth of 10 km and moment magnitude Mw6.0 occurred onshore in central Crete (Greece). The evolution of possible preseismic patterns in the area of central Crete before the Mw6.0 event was investigated by applying the method of multiresolution wavelet analysis (MRWA), along with that of natural time (NT). The monitoring of preseismic patterns by critical parameters defined by NT analysis, integrated with the results of MRWA as the initiation point for the NT analysis, forms a promising framework that may lead to new universal principles that describe the evolution patterns before strong earthquakes. Initially, we apply MRWA to the interevent time series of the successive regional earthquakes in order to investigate the approach of the regional seismicity towards critical stages and to define the starting point of the natural time domain. Then, using the results of MRWA, we apply the NT analysis, showing that the regional seismicity approached criticality for a prolonged period of ~40 days before the occurrence of the Mw6.0 earthquake, when the κ1 natural time parameter reached the critical value of κ1 = 0.070, as suggested by the NT method.

Estimating the Epicenter of a Future Strong Earthquake in Southern California, Mexico, and Central America by Means of Natural Time Analysis and Earthquake Nowcasting

It has recently been shown in the Eastern Mediterranean that by combining natural time analysis of seismicity with earthquake networks based on similar activity patterns and earthquake nowcasting, an estimate of the epicenter location of a future strong earthquake can be obtained. This is based on the construction of average earthquake potential score maps. Here, we propose a method of obtaining such estimates for a highly seismically active area that includes Southern California, Mexico and part of Central America, i.e., the area N1035W80120. The study includes 28 strong earthquakes of magnitude M ≥7.0 that occurred during the time period from 1989 to 2020. The results indicate that there is a strong correlation between the epicenter of a future strong earthquake and the average earthquake potential score maps. Moreover, the method is also applied to the very recent 7 September 2021 Guerrero, Mexico, M7 earthquake as well as to the 22 September 2021 Jiquilillo, Nicaragua, M6.5 earthquake with successful results. We also show that in 28 out of the 29 strong M ≥7.0 EQs studied, their epicenters lie close to an estimated zone covering only 8.5% of the total area.

Estimating the Epicenter of an Impending Strong Earthquake by Combining the Seismicity Order Parameter Variability Analysis with Earthquake Networks and Nowcasting: Application in the Eastern Mediterranean

The variance κ1 of the natural time analysis of earthquake catalogs was proposed in 2005 as an order parameter for seismicity, whose fluctuations proved, in 2011, to be minimized a few months before the strongest mainshock when studying the earthquakes in a given area. After the introduction of earthquake networks based on similar activity patterns, in 2012, the study of their higher order cores revealed, in 2019, the selection of appropriate areas in which the precursory minima βmin of the fluctuations β of the seismicity order parameter κ1 could be observed up to six months before all strong earthquakes above a certain threshold. The eastern Mediterranean region was studied in 2019, where all earthquakes of magnitude M≥7.1 were found to be preceded by βmin without any false alarm. Combining these results with the method of nowcasting earthquakes, introduced in 2016, for seismic risk estimation, here, we show that the epicenter of an impending strong earthquake can be estimated. This is achieved by employing—at the time of observing the βmin—nowcasting earthquakes in a square lattice grid in the study area and by averaging, self-consistently, the results obtained for the earthquake potential score. This is understood in the following context: The minimum βmin is ascertained to almost coincide with the onset of Seismic Electric Signals activity, which is accompanied by the development of long range correlations between earthquake magnitudes in the area that is a candidate for a mainshock.

Seismicity Patterns Prior to the Thessaly (Mw6.3) Strong Earthquake on 3 March 2021 in Terms of Multiresolution Wavelets and Natural Time Analysis

On 3 March 2021, a strong, shallow earthquake of moment magnitude, Mw6.3, occurred in northern Thessaly (Central Greece). To investigate possible complex correlations in the evolution of seismicity in the broader area of Central Greece before the Mw6.3 event, we apply the methods of multiresolution wavelet analysis (MRWA) and natural time (NT) analysis. The description of seismicity evolution by critical parameters defined by NT analysis, integrated with the results of MRWA as the initiation point for the NT analysis, forms a new framework that may possibly lead to new universal principles that describe the generation processes of strong earthquakes. In the present work, we investigate this new framework in the seismicity prior to the Mw6.3 Thessaly earthquake. Initially, we apply MRWA to the interevent time series of the successive regional earthquakes in order to investigate the approach of the regional seismicity at critical stages and to define the starting point of the natural time domain. Then, we apply the NT analysis, showing that the regional seismicity approached criticality a few days before the occurrence of the Mw6.3 earthquake, when the κ1 natural time parameter reached the critical value of κ1 = 0.070.

Magnitude-dependent epidemic-type aftershock sequences model for earthquakes

We propose a version of the pure temporal epidemic type aftershock sequences (ETAS) model: the ETAS model with correlated magnitudes. As for the standard case, we assume the Gutenberg-Richter law to be the probability density for the magnitudes of the background events. Instead, the magnitude of the triggered shocks is assumed to be probabilistically dependent on that of the relative mother events. This probabilistic dependence is motivated by some recent works in the literature and by the results of a statistical analysis made on some seismic catalogs [Spassiani and Sebastiani, J. Geophys. Res. 121, 903 (2016)10.1002/2015JB012398]. On the basis of the experimental evidences obtained in the latter paper for the real catalogs, we theoretically derive the probability density function for the magnitudes of the triggered shocks proposed in Spassiani and Sebastiani and there used for the analysis of two simulated catalogs. To this aim, we impose a fundamental condition: averaging over all the magnitudes of the mother events, we must obtain again the Gutenberg-Richter law. This ensures the validity of this law at any event's generation when ignoring past seismicity. The ETAS model with correlated magnitudes is then theoretically analyzed here. In particular, we use the tool of the probability generating function and the Palm theory, in order to derive an approximation of the probability of zero events in a small time interval and to interpret the results in terms of the interevent time between consecutive shocks, the latter being a very useful random variable in the assessment of seismic hazard.

Multifractal analysis of time series generated by discrete Ito equations

In this study, we show that discrete Ito equations with short-tail Gaussian marginal distribution function generate multifractal time series. The multifractality is due to the nonlinear correlations, which are hidden in Markov processes and are generated by the interrelation between the drift and the multiplicative stochastic forces in the Ito equation. A link between the range of the generalized Hurst exponents and the mean of the squares of all averaged net forces is suggested.

Minima of the fluctuations of the order parameter of global seismicity

It has been recently shown [N. V. Sarlis, Phys. Rev. E 84, 022101 (2011) and N. V. Sarlis and S.-R. G. Christopoulos, Chaos 22, 023123 (2012)] that earthquakes of magnitude M greater or equal to 7 are globally correlated. Such correlations were identified by studying the variance κ1 of natural time which has been proposed as an order parameter for seismicity. Here, we study the fluctuations of this order parameter using the Global Centroid Moment Tensor catalog for a magnitude threshold Mthres = 5.0 and focus on its behavior before major earthquakes. Natural time analysis reveals that distinct minima of the fluctuations of the order parameter of seismicity appear within almost five and a half months on average before all major earthquakes of magnitude larger than 8.4. This phenomenon corroborates the recent finding [N. V. Sarlis et al., Proc. Natl. Acad. Sci. U.S.A. 110, 13734 (2013)] that similar minima of the seismicity order parameter fluctuations had preceded all major shallow earthquakes in Japan. Moreover, on the basis of these minima a statistically significant binary prediction method for earthquakes of magnitude larger than 8.4 with hit rate 100% and false alarm rate 6.67% is suggested.

Natural time analysis of the Centennial Earthquake Catalog

By using the most recent version (1900-2007) of the Centennial Earthquake Catalog, we examine the properties of the global seismicity. Natural time analysis reveals that the fluctuations of the order parameter κ(1) of seismicity exhibit for at least three orders of magnitude a characteristic feature similar to that of the order parameter for other equilibrium or non-equilibrium critical systems-including self-organized critical systems. Moreover, we find non-trivial magnitude correlations for earthquakes of magnitude greater than or equal to 7.

Magnitude correlations in global seismicity

By employing natural time analysis, we analyze the worldwide seismicity and study the existence of correlations between earthquake magnitudes. We find that global seismicity exhibits nontrivial magnitude correlations for earthquake magnitudes greater than M(w) 6.5.

Nonextensivity and natural time: The case of seismicity

Nonextensive statistical mechanics, pioneered by Tsallis, has recently achieved a generalization of the Gutenberg-Richter law for earthquakes. This remarkable generalization is combined here with natural time analysis, which enables the distinction of two origins of self-similarity, i.e., the process' memory and the process' increments infinite variance. By using also detrended fluctuation analysis for the detection of long-range temporal correlations, we demonstrate the existence of both temporal and magnitude correlations in real seismic data of California and Japan. Natural time analysis reveals that the nonextensivity parameter q , in contrast to some published claims, cannot be considered as a measure of temporal organization, but the Tsallis formulation does achieve a satisfactory description of real seismic data for Japan for q=1.66 when supplemented by long-range temporal correlations.