Next Article in Journal
GHZ States as Tripartite PR Boxes: Classical Limit and Retrocausality
Next Article in Special Issue
Statistical Significance of Earth’s Electric and Magnetic Field Variations Preceding Earthquakes in Greece and Japan Revisited
Previous Article in Journal
Approach to Evaluating Accounting Informatization Based on Entropy in Intuitionistic Fuzzy Environment
Previous Article in Special Issue
Criticality Analysis of the Lower Ionosphere Perturbations Prior to the 2016 Kumamoto (Japan) Earthquakes as Based on VLF Electromagnetic Wave Propagation Data Observed at Multiple Stations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Complexity Measures Associated with the Fluctuations of the Entropy in Natural Time before the Deadly México M8.2 Earthquake on 7 September 2017

1
Departamento de Ciencias Básicas, Universidad Autónoma Metropolitana Azcapotzalco, Av. San Pablo 180, Mexico City 02200, Mexico
2
Instituto de Geofísica, Universidad Nacional Autonoma de Mexico, Mexico City 04510, Mexico
3
Section of Solid State Physics and Solid Earth Physics Institute, Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, 157 84 Zografos, Greece
*
Author to whom correspondence should be addressed.
Submission received: 22 May 2018 / Revised: 12 June 2018 / Accepted: 16 June 2018 / Published: 20 June 2018

Abstract

:
We analyse seismicity during the 6-year period 2012–2017 in the new time domain termed natural time in the Chiapas region where the M8.2 earthquake occurred, Mexico’s largest earthquake in more than a century, in order to study the complexity measures associated with fluctuations of entropy as well as with entropy change under time reversal. We find that almost three months before the M8.2 earthquake, i.e., on 14 June 2017, the complexity measure associated with the fluctuations of entropy change under time reversal shows an abrupt increase, which, however, does not hold for the complexity measure associated with the fluctuations of entropy in forward time. On the same date, the entropy change under time reversal has been previously found to exhibit a minimum [Physica A 506, 625–634 (2018)]; we thus find here that this minimum is also accompanied by increased fluctuations of the entropy change under time reversal. In addition, we find a simultaneous increase of the Tsallis entropic index q.

1. Introduction

A widespread opinion is that the two recent deadly Mexico earthquakes, i.e., M8.2 on 7 September 2017 and M7.1 on 19 September 2017, in Chiapas and Morelos regions, respectively, had an unusual cause. In particular, although most big Mexican earthquakes occur along the interface between the subducting Cocos plate and the North American plate, in the case of the M8.2 earthquake that struck the Chiapas state, the earthquake occurred within the Cocos plate itself. Thus, since some seismologists say that this type of faulting should not produce such large earthquakes [1], they characterize it as an “extremely strange” event [2]. As for the M7.1 earthquake, which struck central Mexico on 19 September, it occurred at 57 km depth [3] also within the Cocos plate near the northern limit of the Mexican flat-slab [4,5], where it begins to plunge beneath the North American plate. In short, the two quakes happened at two different spots within the Cocos tectonic plate and surprised seismologists [6].
In a previous paper [7], upon considering the analysis of seismicity in the new time domain termed natural time, we showed that the occurrence of the M8.2 earthquake, Mexico’s largest earthquake in more than a century, should not be considered unexpected. Specifically, this analysis led to the result that in the Chiapas region, where the M8.2 earthquake occurred, the probability for the occurrence of an extreme event was the highest compared to other regions in Mexico. Furthermore, in this region, the same analysis revealed that the entropy change ΔS under time reversal exhibited a pronounced minimum on 14 June 2017, i.e., almost 3 months before the occurrence of the M8.2 earthquake, which pointed to the conclusion that an extreme event was likely to take place there in view of the following: Upon considering the Olami-Feder-Christensen (OFC) model for earthquakes [8], which is probably [9] the most studied non-conservative (supposedly), self-organized criticality (SOC) model (see also [10]), we found that the value of the entropy change under time reversal shows a clear minimum [11,12] before a large avalanche, which corresponds to a large earthquake.
The main purpose of this paper is to investigate the complexity measures associated with the fluctuations of either the entropy defined in natural time and/or the entropy change under time reversal of the seismicity in the Chiapas region from 01 January 2012 until 20 October 2017. In addition, we study whether, beyond the above mentioned precursory ΔS minimum reported in [7], there exist also similar changes in the complexity measures investigated here. These measures have been found to be useful when the system approaches the critical point (dynamic phase transition) as, for example, is the case in an impending sudden cardiac death risk ([13,14], see also Section 9.4.1 of Ref. [11]).

2. Natural Time Analysis. The Entropy Defined in Natural Time and the Associated Complexity Measures

For a time series of N events, we define an index for the occurrence of the k-th event by χk = k/N, which we term natural time χ. In this analysis [11,15], we ignore the time intervals between consecutive events, but preserve their order and energy Qk. We then study the pairs (χk, Qk) or the pairs (χk, pk) where
p k = Q k k = 1 N Q k
is the normalized energy for the k-th event. The entropy S in natural time is defined [13,16,17] by
S = χ l n χ χ l n χ
where the brackets <…> ≡ ∑(…)pk denote averages with respect to the distribution pk, i.e., <f(χ)> ≡ ∑f(χk)pk. The entropy S is a dynamic entropy that exhibits [18] concavity, positivity, and Lesche stability [19,20].
Upon considering time reversal T ^ , i.e., T ^ pk = pN−k+1, the value of S changes to a value S:
S _ = k = 1 N k N ln ( k N ) p N k + 1 ( k = 1 N k N p N k + 1 ) l n ( l = 1 N l N p N l + 1 )
The physical meaning of the change of entropy ∆SS − S in natural time under time reversal has been discussed in References. [11,21,22].
Using a moving window of length i (number of consecutive events) sliding through the time series of L consecutive events, the entropy in natural time has been determined for each position j = 1, 2, …, L − i of the sliding window. Thus, a time series of Si has been constructed [13] whose standard deviation is designated by δSi. The study of the effect of the change of scale i on δSi is made [14] by means of the complexity measure
λ i =   δ S i δ S 100
where the denominator δS100 is arbitrarily selected to stand for the δS value of a short scale, i.e., 100 events, while the numerator corresponds to a longer scale, e.g., i = 103 events. If instead of δSi the standard deviation σSi) of the time series of ΔSiSi (S)i is used, we define [11,23] the complexity measure Λi,
Λ i =   σ ( Δ S i ) σ ( Δ S 100 ) ,
when a moving window of i consecutive events is sliding through the time series and the denominator σS100) is arbitrarily selected to correspond to the standard deviation σS100) of the time series of ΔSi of i = 100 events. In other words, this complexity measure quantifies how the statistics of ∆Si time series changes upon increasing the scale from 100 events to a longer scale e.g., i = 103 events.

3. Data and Analysis

The seismic data analysed in natural time come from the seismic catalog of the National Seismic Service (SSN) of the Universidad Nacional Autónoma de México (www.ssn.unam.mx) and cover the period 1 January 2012 until 20 October 2017. Since in this analysis Qk should be [11,15] proportional to the energy emitted during the k-th earthquake of magnitude Mk, we assumed Q k 10 1.5 M k [11]. To assure catalog completeness a magnitude threshold M ≥ 3.5 has been imposed. Furthermore, in order to investigate whether our results obey magnitude threshold invariance, we also repeat our calculations for M ≥ 4.0, as will be explained in the next Section.
Our calculations are made, as mentioned in the previous section, by means of a window of length i (= number of successive events) sliding each time by one event through the whole time series. The entropies S and S, and therefrom their difference ∆Si, are calculated each time; we thus also form a new time series consisting of successive ∆Si values. The complexity measures λi and Λi are determined according to their definitions given in Equations (4) and (5), respectively.

4. Results

In Figure 1a–c we plot the values calculated for the quantities Si, (S)i, and ∆Si, respectively, versus the conventional time for all M ≥ 3.5 earthquakes in the Chiapas region during the period 1 January 2012 to the date of occurrence of the M8.2 earthquake for the scales 102, 3 × 103, and 4 × 103 events. The study of the first scale (102 events) is needed for the calculation of the denominator of Equations (4) and (5), while the selection of the other scales (3 × 103 events and longer) was made for the following reasons, as also explained in Ref. [7]. Since ∼11,500 earthquakes (M ≥ 3.5) occurred in this area from 1 January 2012 until the occurrence of the M8.2 earthquake on 7 September 2017, we find an average of around 170 earthquakes per month. We take into account that recent investigations by means of natural time analysis revealed that the fluctuations of the order parameter of seismicity exhibit [24] a minimum when a series of precursory low frequency (≤0.1 Hz) electric signals, termed Seismic Electric Signals (SES) activity (e.g., see [22]) (whose lead time is up to around 5.5 months [11]), is initiated. While this minimum of the order parameter of seismicity is observed during a period in which long-range correlations prevail between earthquake magnitudes, another stage appears before this minimum in which the temporal correlations between earthquake magnitudes exhibit a distinctly different behaviour, i.e., an evident anticorrelated behaviour [25]. The significant change between these two stages in the temporal correlations between earthquake magnitudes is likely to be captured by the time evolution of ΔSi, hence we started the presentation of our study of ΔSi in Ref. [7] from the scale of i ∼ 103 events, i.e., around the maximum lead time of SES activities. This study led to the conclusion that at scales i = 3 × 103 or longer (e.g., 4 × 103 and 5 × 103 events), a pronounced minimum becomes evident at the date 14 June 2017 (when a M7 earthquake occurred, i.e., almost 3 months before the M8.2 earthquake that struck Mexico’s Chiapas state). Interestingly, this minimum of ΔSi was found to exhibit magnitude threshold invariance. Hence, in the following we will present our results for these scales and in addition focus our attention on what happens before the occurrence of the M7 earthquake on 14 June 2017 and before the M8.2 earthquake on 7 September 2017.
In particular, we will investigate the values of the complexity measures λi and Λi on the following 6 dates: On 1 June 2017 (i.e., almost two weeks before the occurrence of the M7 event on 14 June), on 14 June 2017 (upon the occurrence of the last event before the M7 earthquake on 14 June 2017), on 1 July 2017, on 1 August 2017, on 1 September 2017 (6 days before the M8.2 earthquake), and on 7 September 2017 (upon the occurrence of the last small event before the M8.2 earthquake).

4.1. Results for the Complexity Measure λi

We will now present our results for the complexity measure λi, which is solely associated with the fluctuations δSi of the entropy in forward time. In Figure 2, we plot the λi values versus the conventional time by considering all M ≥ 3.5 earthquakes from 1 January 2012 until the occurrence of the M8.2 earthquake. A close inspection of this figure in all scales investigated does not reveal any remarkable change before the occurrence of the M8.2 earthquake. The same conclusion is drawn when we plot in Figure 3 the λi values versus the scale i of the number of M ≥ 3.5 events that occurred in the Chiapas region from 1 January 2012 until the dates mentioned above before the two earthquakes, i.e., the M7 earthquake on 14 June 2017 and the M8.2 earthquake on 7 September 2017. Remarkably, the resulting values for all six dates coincide for each i value without showing any precursory variation. In other words, even when considering the λi value on 7 September 2017, upon the occurrence of a small event just before the M8.2 earthquake, the value at i = 5000 does not significantly differ from the other λi values that correspond, for example, to that of 1 June 2017, i.e., more than 3 months before the M8.2 earthquake occurrence.

4.2. Results for the Complexity Measure Λi

In Figure 4 we plot the Λi values versus the conventional time by considering all M ≥ 3.5 earthquakes in the Chiapas region from 1 January 2012 until the occurrence of the M8.2 earthquake. An inspection of this figure reveals that upon the occurrence of the M7 earthquake on 14 June 2017 an abrupt increase of Λi is observed in all three scales.
Furthermore, in Figure 5 we plot the Λi values versus the scale i of the number of M ≥ 3.5 events in the Chiapas region from 1 January 2012 until the six dates mentioned above before the two earthquakes, i.e., the M7 earthquake on 14 June 2017 and the M8.2 earthquake on 7 September 2017. We observe that for the scales i = 3000 events and larger an evident abrupt increase of the Λi value is observed upon the occurrence of the M7 earthquake on 14 June 2017. This date of the abrupt increase of Λi remains invariant upon changing the magnitude threshold, for example, see the result also plotted in Figure 5 by increasing the magnitude threshold to M ≥ 4.0 instead of M ≥ 3.5 and considering that the number of earthquakes then decreases by a factor of around 4.

5. Discussion

While the complexity measure Λi does exhibit an evident precursory change almost three months before the M8.2 earthquake, the complexity measure λi does not. This should not be considered as surprising in view of the following: Λi is associated with the fluctuations of the entropy change under time reversal, and we know that ∆Si is a key quantity to determine the time of an impending dynamic phase transition [20] as, for example, is the case in sudden cardiac death risk. This is also the case with earthquakes because the observed earthquake scaling laws (e.g., [26]) indicate the existence of phenomena closely associated with the proximity of the system to a critical point. Thus, taking the view that a strong earthquake is a critical phenomenon (dynamic phase transition), it may not come as a surprise that the fluctuations of ∆Si, and hence Λi, may serve for the estimation of the time of its occurrence. We note that in [12] the predictability of the OFC model was studied by using the entropy change under time reversal. It was found that the value of ΔSi exhibits a clear minimum (if we define ΔS = SS) [7] or maximum (if we define ΔS = SS) [12] before large avalanches in the OFC model, thus this minimum provides a decision variable for the prediction of a large avalanche. The prediction quality of this minimum was studied, following [27], by choosing the receiver operating characteristics (ROC) graph [28], see Figure 2 of [12]. The ROC graph was constructed for the OFC model, and it was found that the predictions made on the basis of ΔS are statistically significant, and, as concerns the origin of their predictive power, this should be attributed to the fact that ΔS is able to catch the ‘true’ time arrow.
To further study the increase of Λi possibly associated with strong earthquakes, we now investigate the following: At first glance, one may claim that the behavior of Λi observed on 14 June 2017 seems to be very similar to what happens around the end of 2015, and both changes seem to be related to earthquakes around these times. In order to examine whether such a similarity is true, we now depict in Figure 6a the complexity measures Λi versus the conventional time in an expanded time scale for M ≥ 3.5. A close inspection of this figure reveals that all the three complexity measures Λ3000, Λ4000, and Λ5000 (which stand for the Λi values at the scales i = 3000, 4000, and 5000, respectively) exhibit a strong and abrupt increase on 14 June 2017 (after the occurrence of a M7.0 earthquake). Such a strong and abrupt increase for all these three complexity measures is not observed at any other time. In particular, on 17 December 2015 after the occurrence of a M6.6 earthquake, only Λ3000 and Λ5000 exhibit a strong and abrupt increase, while the corresponding increase for Λ4000 is much smaller and obviously does not scale with i. The behavior of the increase ΔΛi (≡Λi(t) − Λi(tEQ)) of the complexity measure Λi upon the occurrence of a strong earthquake at tEQ is further studied in Figure 6b,c for the earthquakes of 14 June 2017 and 17 December 2015, respectively. An inspection of these two figures shows that while the abrupt increase on 14 June 2017 exhibits a scaling behavior of the form
Δ Λ i = A   ( t t 0 ) c ,
where the pre-factors A are proportional to i and the exponent c is independent of i with a value close to 0.33, this does not hold for the case of 17 December 2015. We also comment that quite interestingly there exists a striking similarity between the results shown in Figure 6b and the seminal work by Lifshitz and Slyozov [29], and independently by Wagner [30], who found a form like that of Equation (6) with an exponent close to 1/3 for the time growth of the characteristic size of the minority phase droplets in phase transitions (e.g., [31]).
Concerning the statistical significance of the abrupt increase of the complexity measure Λi three months before the M8.2 Chiapas earthquake as a precursor to strong earthquakes, we employ the most recent method of event coincidence analysis [32] as implemented by the CoinCalc package in R. More specifically, for every month we estimate the largest earthquake magnitude Mmax observed and compare it with a threshold M t h r e s t a r g e t , if M max > M t h r e s t a r g e t a strong earthquake event occurred in that month. The precursor event time series is composed by the abrupt increase in events of Λi (see Figure 4 and Figure 7a) and exhibits two events for i = 5000 and 3000 (one in December 2015 and one in June 2017), while for i = 4000 a single event (in June 2017). Figure 7 shows how the corresponding p-values vary upon increasing the threshold M t h r e s t a r g e t to obtain the real situation by chance. We observe that only in three out of the nine results shown, these p-values exceed 10%, while for the largest threshold M t h r e s t a r g e t the corresponding p-values are 4.8%, 2.8%, and 7.0% for i = 3000, 4000, and 5000, respectively, pointing to the statistical significance of the precursor under discussion.
Similar results to those obtained here have been found for other time series prior to large earthquakes in other places. As a characteristic example we mention the M9 Tohoku mega-earthquake that occurred in Japan on 11 March 2011. In particular, the complexity measure Λi, associated with the fluctuations of the entropy change of the seismicity in Japan (with M ≥ 3.5), was studied during the period from 2000 to 11 March 2011 and found to exhibit an abrupt increase upon the occurrence of a M7.8 earthquake on 22 December 2010. Remarkably, on the same date [25] the temporal correlations between earthquake magnitudes exhibited anticorrelated behavior with the lowest value (≈0.35) ever observed of the exponent in the detrended fluctuation analysis [33,34,35] we employed in that study. Details on the interconnection between the changes in the temporal correlations between earthquake magnitudes and the precursory ΔS minimum in Japan will be published shortly elsewhere.
Finally, we comment on the results obtained on the basis of the pioneering work of Tsallis [36], which introduced a non-additive entropy Sq (e.g., see [37,38]). This has found application in the physics of earthquakes and especially in the description of the asperities [39,40] in the faults on which earthquakes occur through the entropic index q. Based on the earthquake magnitude distribution one can obtain [41,42] entropic index q and study its variation with time as we approach a strong EQ (for a recent review see [43]). Figure 8 depicts the q-value versus conventional time as it is estimated [42] for sliding windows of W = 1000 and W = 250 consecutive earthquakes in the Chiapas region for M ≥ 3.5 and M ≥ 4.0, shown as red and the blue lines, respectively. We observe that before the occurrence of the M8.2 Chiapas earthquake, the q-value starts to grow gradually during May 2017 and exhibits an abrupt increase upon the occurrence of the earthquake on 14 June 2017 for both magnitude thresholds. This behavior is reminiscent to that observed before the 1995 Kobe earthquake in Japan (e.g., see Figure 3 of [44]).

6. Conclusions

It was found that the complexity measure λi studied during the period 2012–2017, which is solely associated with the fluctuations of entropy in forward time, does not exhibit any obvious precursory change before the M8.2 earthquake on 7 September 2017.
On the other hand, the complexity measure Λi associated with the fluctuations of entropy change under time reversal shows an evident increase on 14 June 2017, accompanied by an abrupt increase of the Tsallis entropic index q. Interestingly, on the same date reported by Sarlis et al. in [7], the entropy change under time reversal has been found to exhibit a minimum, as was found in the OFC model [12].

Author Contributions

Conceptualization, A.R.-R., E.L.F.-M., P.A.V., N.V.S.; Methodology, P.A.V., N.V.S.; Software, P.A.V., N.V.S., A.R.-R., E.L.F.-M.; Formal Analysis, N.V.S., A.R.-R., P.A.V., E.L.F.-M.; Investigation, A.R.-R., E.L.F.-M., P.A.V., N.V.S.; Data Curation, A.R.-R., E.L.F.-M.; Writing-Original Draft Preparation, A.R.-R., N.V.S.; Writing-Review & Editing, P.A.V., N.V.S.; Project Administration, A.R.-R., N.V.S.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mega, E.R. Deadly Mexico earthquake had unusual cause. Nature 2017, 549. [Google Scholar] [CrossRef]
  2. Wade, L. Unusual quake rattles Mexico. Science 2017, 357, 1084. [Google Scholar] [CrossRef] [PubMed]
  3. Witze, A. Pair of deadly Mexico quakes puzzles scientists. Nature 2017, 549. [Google Scholar] [CrossRef] [PubMed]
  4. Pérez-Campos, X.; Kim, Y.; Husker, A.; Davis, P.M.; Clayton, R.W.; Iglesias, A.; Pacheco, J.F.; Singh, S.K.; Manea, V.C.; Gurnis, M. Horizontal subduction and truncation of the Cocos Plate beneath central Mexico. Geophys. Res. Lett. 2008, 35, l18303. [Google Scholar] [CrossRef]
  5. Manea, V.; Manea, M.; Ferrari, L.; Orozco-Esquivel, T.; Valenzuela, R.; Husker, A.; Kostoglodov, V. A review of the geodynamic evolution of flat slab subduction in Mexico, Peru, and Chile. Tectonophysics 2016, 695 (Suppl. C), 27–52. [Google Scholar] [CrossRef]
  6. Witze, A. Deadly Mexico quakes not linked. Nature 2017, 549, 442. [Google Scholar] [CrossRef] [PubMed]
  7. Sarlis, N.V.; Skordas, E.S.; Varotsos, P.A.; Ramírez-Rojas, A.; Flores-Márquez, E.L. Natural time analysis: On the deadly Mexico M8.2 earthquake on 7 September 2017. Phys. A 2018, 506, 625–634. [Google Scholar] [CrossRef]
  8. Olami, Z.; Feder, H.J.S.; Christensen, K. Self-organized criticality in a continuous, nonconservative cellular automaton modeling earthquakes. Phys. Rev. Lett. 1992, 68, 1244–1247. [Google Scholar] [CrossRef] [PubMed]
  9. Ramos, O.; Altshuler, E.; Måløy, K.J. Quasiperiodic events in an earthquake model. Phys. Rev. Lett. 2006, 96, 098501. [Google Scholar] [CrossRef] [PubMed]
  10. Caruso, F.; Kantz, H. Prediction of extreme events in the OFC model on a small world network. Eur. Phys. J. 2011, 79, 7–11. [Google Scholar] [CrossRef]
  11. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Natural Time Analysis: The new view of time. In Precursory Seismic Electric Signals, Earthquakes and Other Complex Time-Series; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  12. Sarlis, N.; Skordas, E.; Varotsos, P. The change of the entropy in natural time under time-reversal in the Olami-Feder-Christensen earthquake model. Tectonophysics 2011, 513, 49–53. [Google Scholar] [CrossRef]
  13. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Entropy in natural time domain. Phys. Rev. E 2004, 70, 011106. [Google Scholar] [CrossRef] [PubMed]
  14. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Natural entropy fluctuations discriminate similar-looking electric signals emitted from systems of different dynamics. Phys. Rev. E 2005, 71, 011110. [Google Scholar] [CrossRef] [PubMed]
  15. Varotsos, P.; Sarlis, N.V.; Skordas, E.S.; Uyeda, S.; Kamogawa, M. Natural time analysis of critical phenomena. Proc. Natl. Acad. Sci. USA 2011, 108, 11361–11364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Spatio-temporal complexity aspects on the interrelation between seismic electric signals and seismicity. Pract. Athens Acad. 2001, 76, 294–321. [Google Scholar]
  17. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Attempt to distinguish electric signals of a dichotomous nature. Phys. Rev. E 2003, 68, 031106. [Google Scholar] [CrossRef] [PubMed]
  18. Varotsos, P.A.; Sarlis, N.V.; Tanaka, H.K.; Skordas, E.S. Some properties of the entropy in the natural time. Phys. Rev. E 2005, 71, 032102. [Google Scholar] [CrossRef] [PubMed]
  19. Lesche, B.J. Instabilities of Rényi entropies. Stat. Phys. 1982, 27, 419–422. [Google Scholar] [CrossRef]
  20. Lesche, B. Rényi entropies and observables. Phys. Rev. E 2004, 70, 017102. [Google Scholar] [CrossRef] [PubMed]
  21. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Identifying sudden cardiac death risk and specifying its occurrence time by analysing electrocardiograms in natural time. Appl. Phys. Lett. 2007, 91, 064106. [Google Scholar] [CrossRef]
  22. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Fluctuations, under time reversal, of the natural time and the entropy distinguish similar looking electric signals of different dynamics. J. Appl. Phys. 2008, 103, 014906. [Google Scholar] [CrossRef] [Green Version]
  23. Nicholas, V.; Sarlis, S.; Christopoulos, R.G.; Bemplidaki, M.M. Change ΔS of the entropy in natural time under time reversal: Complexity measure upon change of scale. Eur. Lett. 2015, 109, 18002. [Google Scholar] [CrossRef]
  24. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Seismic Electric Signals: An additional fact showing their physical interconnection with seismicity. Tectonophysics 2013, 589, 116–125. [Google Scholar] [CrossRef]
  25. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Study of the temporal correlations in the magnitude time series before major earthquakes in Japan. J. Geophys. Res. 2014, 119, 9192–9206. [Google Scholar] [CrossRef] [Green Version]
  26. Turcotte, D.L. Fractals and Chaos in Geology and Geophysics, 2nd ed.; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  27. Garber, A.; Hallerberg, S.; Kantz, H. Predicting extreme avalanches in self-organized critical Sandpiles. Phys. Rev. E 2009, 80, 026124. [Google Scholar] [CrossRef] [PubMed]
  28. Fawcett, T. An Introduction to ROC Analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  29. Lifshitz, I.M.; Slyozov, V.V. The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. Solids 1961, 19, 35–50. [Google Scholar] [CrossRef]
  30. Wagner, C. Theorie der Alterung von NiederschlagendurchUmlosen (Ostwald-Reifung). Z. Elektrochem. 1981, 65, 581–591. [Google Scholar]
  31. Bray, A.J. Theory of phase-ordering kinetics. Adv. Phys. 1994, 43, 357–459. [Google Scholar] [CrossRef] [Green Version]
  32. Siegmund, J.F.; Siegmund, N.; Donner, R.V. CoinCalc—A new R package for quantifying simultaneities of event series. Comput. Geosci. 2017, 98, 64–72. [Google Scholar] [CrossRef]
  33. Peng, C.K.; Buldyrev, S.V.; Havlin, S.; Simons, M.; Stanley, H.E.; Goldberger, A.L. Mosaic organization of DNA nucleotides. Phys. Rev. E 1944, 49, 1685–1689. [Google Scholar] [CrossRef]
  34. Peng, C.K.; Buldyrev, S.V.; Goldberger, A.L.; Havlin, S.; Mantegna, R.N.; Simons, M.; Stanley, H.E. Statisticalproperties of DNA sequences. Phys. A 1995, 221, 180–192. [Google Scholar] [CrossRef]
  35. Peng, C.K.; Havlin, S.; Stanley, H.E.; Goldberger, A.L. Quantification of scaling exponents and crossoverphenomena in nonstationary heartbeat time series. Chaos 1995, 5, 82–87. [Google Scholar] [CrossRef] [PubMed]
  36. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  37. Tsallis, C. Introduction to Nonextensive Statistical Mechanics; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  38. Tsallis, C. The Nonadditive Entropy Sq and Its Applications in Physics and Elsewhere: Some Remarks. Entropy 2011, 13, 1765–1804. [Google Scholar] [CrossRef] [Green Version]
  39. Sotolongo-Costa, O.; Posadas, A. Fragment-Asperity Interaction Model for Earthquakes. Phys. Rev. Lett. 2004, 92, 048501. [Google Scholar] [CrossRef] [PubMed]
  40. Silva, R.; França, G.S.; Vilar, C.S.; Alcaniz, J.S. Nonextensive models for earthquakes. Phys. Rev. E 2006, 73, 026102. [Google Scholar] [CrossRef] [PubMed]
  41. Sarlis, N.V.; Skordas, E.S.; Varotsos, P.A. Nonextensivity and natural time: The case of seismicity. Phys. Rev. E 2010, 82, 021110. [Google Scholar] [CrossRef] [PubMed]
  42. Telesca, L. Maximum Likelihood Estimation of the Nonextensive Parameters of the Earthquake Cumulative 375 Magnitude Distribution. Bull. Seismol. Soc. Am. 2012, 102, 886–891. [Google Scholar] [CrossRef]
  43. Vallianatos, F.; Michas, G.; Papadakis, G. Nonextensive Statistical Seismology: An Overview. In Complexity of Seismic Time Series; Chelidze, T., Vallianatos, F., Telesca, L., Eds.; Elsevier: New York, NY, USA, 2018; pp. 25–59. [Google Scholar]
  44. Papadakis, G.; Vallianatos, F.; Sammonds, P. A Nonextensive Statistical Physics Analysis of the 1995 Kobe, 379 Japan Earthquake. Pure Appl. Geophys. 2015, 172, 1923–1931. [Google Scholar] [CrossRef]
Figure 1. Plot of the entropies Si (a) and (S)i (b) as well as the entropy change ∆Si under time reversal (c) versus the conventional time for the three scales i = 102 (red), 3 × 103 (green), and 4 × 103 (blue) events when analysing all earthquakes with M ≥ 3.5. The vertical lines ending at circles depict the earthquake magnitudes, which are read in the right scale.
Figure 1. Plot of the entropies Si (a) and (S)i (b) as well as the entropy change ∆Si under time reversal (c) versus the conventional time for the three scales i = 102 (red), 3 × 103 (green), and 4 × 103 (blue) events when analysing all earthquakes with M ≥ 3.5. The vertical lines ending at circles depict the earthquake magnitudes, which are read in the right scale.
Entropy 20 00477 g001
Figure 2. Plot of the values of the complexity measure λi versus the conventional time that correspond to the scales i = 3 × 103 (green), 4 × 103 (blue), and 5 × 103 (cyan) events when considering all earthquakes in the Chiapas region with M ≥ 3.5 since 2012.
Figure 2. Plot of the values of the complexity measure λi versus the conventional time that correspond to the scales i = 3 × 103 (green), 4 × 103 (blue), and 5 × 103 (cyan) events when considering all earthquakes in the Chiapas region with M ≥ 3.5 since 2012.
Entropy 20 00477 g002
Figure 3. Plot of the complexity measure λi versus the scale i (number of events) for all M ≥ 3.5 earthquakes in the Chiapas region since 1 January 2012. The λi values are calculated for each i value at the following dates: 1 June 2017 (yellow solid circles), 14 June 2017 (cyan squares), 1 July 2017 (red plus), 1 August 2017 (blue star), 1 September 2017 (green cross), and 7 September 2017 (red circle, until the last event before the M8.2 earthquake on 7 September 2017).
Figure 3. Plot of the complexity measure λi versus the scale i (number of events) for all M ≥ 3.5 earthquakes in the Chiapas region since 1 January 2012. The λi values are calculated for each i value at the following dates: 1 June 2017 (yellow solid circles), 14 June 2017 (cyan squares), 1 July 2017 (red plus), 1 August 2017 (blue star), 1 September 2017 (green cross), and 7 September 2017 (red circle, until the last event before the M8.2 earthquake on 7 September 2017).
Entropy 20 00477 g003
Figure 4. Plot of the values of the complexity measure Λi versus the conventional time that correspond to the scales i = 3 × 103 (green), 4 × 103 (blue), and 5 × 103 (cyan) events when considering all earthquakes in the Chiapas region with M ≥ 3.5 since 2012.
Figure 4. Plot of the values of the complexity measure Λi versus the conventional time that correspond to the scales i = 3 × 103 (green), 4 × 103 (blue), and 5 × 103 (cyan) events when considering all earthquakes in the Chiapas region with M ≥ 3.5 since 2012.
Entropy 20 00477 g004
Figure 5. Plot of the complexity measure Λi versus the scale i (number of events) for all M ≥ 3.5 earthquakes (a) as well as for all M ≥ 4.0 earthquakes (b) in the Chiapas region since 1 January 2012. The Λi values are calculated for each i value at the following dates: 1 June 2017 (yellow solid circles), 14 June 2017 (cyan squares), 1 July 2017 (red plus), 1 August 2017 (blue star), 1 September 2017 (green cross) and 7 September 2017 (red circle, until the last event before the M8.2 earthquake on 7 September 2017).
Figure 5. Plot of the complexity measure Λi versus the scale i (number of events) for all M ≥ 3.5 earthquakes (a) as well as for all M ≥ 4.0 earthquakes (b) in the Chiapas region since 1 January 2012. The Λi values are calculated for each i value at the following dates: 1 June 2017 (yellow solid circles), 14 June 2017 (cyan squares), 1 July 2017 (red plus), 1 August 2017 (blue star), 1 September 2017 (green cross) and 7 September 2017 (red circle, until the last event before the M8.2 earthquake on 7 September 2017).
Entropy 20 00477 g005
Figure 6. (a) The same as in Figure 4 but in an expanded timescale; (b,c) plot of the change ΔΛi of the complexity measure Λi after the occurrence of the M7.0 earthquake on 14 June 2017 (b) and the M6.6 earthquake on 17 December 2015 (c) versus the time elapsed (tt0) in days. In both cases, the value of t0 has been selected approximately 30 min after earthquake occurrence.
Figure 6. (a) The same as in Figure 4 but in an expanded timescale; (b,c) plot of the change ΔΛi of the complexity measure Λi after the occurrence of the M7.0 earthquake on 14 June 2017 (b) and the M6.6 earthquake on 17 December 2015 (c) versus the time elapsed (tt0) in days. In both cases, the value of t0 has been selected approximately 30 min after earthquake occurrence.
Entropy 20 00477 g006aEntropy 20 00477 g006b
Figure 7. Results from the event coincidence analysis between the earthquake events of magnitude M > M t h r e s t a r g e t and the events of abrupt increase of Λi (see Figure 4 and Figure 6a). The probability to obtain by chance (p-value) the observed results when considering the events of abrupt increase of Λi as precursors to strong earthquakes versus the threshold M t h r e s t a r g e t for various scales i. The p-values have been obtained by using the compute code eca.es of the CoinCalc package [32] of R.
Figure 7. Results from the event coincidence analysis between the earthquake events of magnitude M > M t h r e s t a r g e t and the events of abrupt increase of Λi (see Figure 4 and Figure 6a). The probability to obtain by chance (p-value) the observed results when considering the events of abrupt increase of Λi as precursors to strong earthquakes versus the threshold M t h r e s t a r g e t for various scales i. The p-values have been obtained by using the compute code eca.es of the CoinCalc package [32] of R.
Entropy 20 00477 g007
Figure 8. The entropic index q versus conventional time for sliding windows W = 1000 (red) and W = 250 (blue) consecutive earthquakes in the Chiapas region since 2012 for M ≥ 3.5 and M ≥ 4.0, respectively.
Figure 8. The entropic index q versus conventional time for sliding windows W = 1000 (red) and W = 250 (blue) consecutive earthquakes in the Chiapas region since 2012 for M ≥ 3.5 and M ≥ 4.0, respectively.
Entropy 20 00477 g008

Share and Cite

MDPI and ACS Style

Ramírez-Rojas, A.; Flores-Márquez, E.L.; Sarlis, N.V.; Varotsos, P.A. The Complexity Measures Associated with the Fluctuations of the Entropy in Natural Time before the Deadly México M8.2 Earthquake on 7 September 2017. Entropy 2018, 20, 477. https://0-doi-org.brum.beds.ac.uk/10.3390/e20060477

AMA Style

Ramírez-Rojas A, Flores-Márquez EL, Sarlis NV, Varotsos PA. The Complexity Measures Associated with the Fluctuations of the Entropy in Natural Time before the Deadly México M8.2 Earthquake on 7 September 2017. Entropy. 2018; 20(6):477. https://0-doi-org.brum.beds.ac.uk/10.3390/e20060477

Chicago/Turabian Style

Ramírez-Rojas, Alejandro, Elsa Leticia Flores-Márquez, Nicholas V. Sarlis, and Panayiotis A. Varotsos. 2018. "The Complexity Measures Associated with the Fluctuations of the Entropy in Natural Time before the Deadly México M8.2 Earthquake on 7 September 2017" Entropy 20, no. 6: 477. https://0-doi-org.brum.beds.ac.uk/10.3390/e20060477

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop