Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The Adaptive Behavior of a Soccer Team: An Entropy-Based Analysis
Next Article in Special Issue
Non-Extensive Statistical Analysis of Energetic Particle Flux Enhancements Caused by the Interplanetary Coronal Mass Ejection-Heliospheric Current Sheet Interaction
Previous Article in Journal
Creation of Numerical Constants in Robust Gene Expression Programming
Previous Article in Special Issue
High Density Nodes in the Chaotic Region of 1D Discrete Maps
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tsallis Entropy Index q and the Complexity Measure of Seismicity in Natural Time under Time Reversal before the M9 Tohoku Earthquake in 2011

by
Panayiotis A. Varotsos
1,2,
Nicholas V. Sarlis
1,2,* and
Efthimios S. Skordas
1,2
1
Section of Solid State Physics, Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos, 157 84 Athens, Greece
2
Solid Earth Physics Institute, Department of Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos, 157 84 Athens, Greece
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(10), 757; https://doi.org/10.3390/e20100757
Submission received: 24 July 2018 / Revised: 23 September 2018 / Accepted: 30 September 2018 / Published: 2 October 2018
(This article belongs to the Special Issue Theoretical Aspects of Kappa Distributions)

Abstract

:
The observed earthquake scaling laws indicate the existence of phenomena closely associated with the proximity of the system to a critical point. Taking this view that earthquakes are critical phenomena (dynamic phase transitions), here we investigate whether in this case the Lifshitz–Slyozov–Wagner (LSW) theory for phase transitions showing that the characteristic size of the minority phase droplets grows with time as t 1 / 3 is applicable. To achieve this goal, we analyzed the Japanese seismic data in a new time domain termed natural time and find that an LSW behavior is actually obeyed by a precursory change of seismicity and in particular by the fluctuations of the entropy change of seismicity under time reversal before the Tohoku earthquake of magnitude 9.0 that occurred on 11 March 2011 in Japan. Furthermore, the Tsallis entropic index q is found to exhibit a precursory increase.

1. Introduction

Earthquakes exhibit complex correlations in time, space and magnitude and have been the object of a multitude of studies [1,2,3,4,5,6,7]. The observed earthquake scaling laws (e.g., [8]) are considered to indicate the existence of phenomena closely associated with the proximity of the system to a critical point (e.g., [9]). Accepting this view that earthquakes are critical phenomena (dynamic phase transitions, where a mainshock is the new phase), the quantity by which one can identify the approach of a dynamical system to the state of criticality is termed order parameter. This parameter in the frame of a new time domain termed natural time χ [10] is just the quantity κ 1 = χ 2 χ 2 , as explained in Section 2 (cf. pp. 249–254 of Ref. [11]). The following two key properties have been shown [10,11] for the probability density function (PDF) P ( κ 1 ) of the κ 1 values in an earthquake catalog: First, different seismic areas σ P ( κ 1 ) versus ( μ κ 1 ) / σ —where μ stands for the mean of the κ 1 values and σ for their standard deviation—fall on a non-Gaussian universal curve which has a left exponential tail, showing that an extreme fluctuation may be orders of magnitude more probable than it would be if Gaussianity was valid, pointing to the existence of extreme events. Second, the PDF P ( κ 1 ) versus κ 1 before strong earthquakes exhibits a bimodal feature. This, for example, has been observed before the M7.3 Landers and before the M7.1 Hector Mine earthquakes that occurred in Southern California in 1992 and 1999, respectively (see pp. 274 and 278 of Ref. [11]). Another very recent example studied in Ref. [12] is the case of the M8.2 earthquake that occurred in the Chiapas region in Mexico on 7 September 2017 to which we now turn.
In the case of Mexico, the seismicity has been studied in natural time in Ref. [13] in the six tectonic regions Baja California (BC), Jalisco-Colima (J), Michoacán(M), Guerrero (G), Oaxaca (O) and Chiapas (CH) of the Mexican Pacific Coast shown in Figure 2 of Ref. [12]. This study showed that only for earthquakes in the regions CH, G and O a bimodal feature appears in the PDF P ( κ 1 ) versus κ 1 (see their Figure 3a). Among these three regions, the first one, i.e., CH, the PDF P ( κ 1 ) vs κ 1 of which is shown in Figure 3 of Ref. [12], had the highest probability for an extreme fluctuation (large earthquake) as can be seen by comparing their left exponential tail of the σ P ( κ 1 ) versus ( μ κ 1 ) / σ plot depicted in Figure 4 of Ref. [12], where the results of all six regions are shown. These results reveal that in principle extreme events in the Chiapas region were expected from the natural time analysis, as actually happened upon the occurrence of the M8.2 earthquake in this region.
In addition, by employing natural time analysis in Ref. [12], it was found that the entropy change Δ S under time reversal of the seismicity in the Chiapas area exhibited a clear minimum on 14 June 2017, which signaled that a major event was impending there as actually happened almost three months later with the occurrence of the M8.2 earthquake on 7 September 2017. This is so, because Δ S constitutes a key measure that may identify when the system approaches the critical point (dynamic phase transition) [11]. For example, Δ S has been applied [14] for the identification of the impending sudden cardiac death risk. Furthermore, it has been shown that Δ S provides a useful tool [15] to investigate the predictability of the Olami–Feder–Christensen (OFC) model for earthquakes [16], which is probably [17] the most studied non-conservative self-organized criticality model. In particular, it was found that Δ S exhibits a clear minimum [11] (or maximum if we define [15] Δ S S S instead of Δ S S S used here, see Section 2) before a large avalanche in the OFC model, which corresponds to a large earthquake. Furthermore, in a more recent study [18], by analyzing the seismicity during the six-year period 2012–2017 in natural time in the Chiapas region where the M8.2 earthquake occurred, we found that, on the same date as above, i.e., 14 June 2017, the complexity measure Λ i (see Section 2) associated with the fluctuations of the entropy change under time reversal exhibits an abrupt increase along with a simultaneous increase of the Tsallis entropic index q [19,20,21,22,23,24].The temporal variations of the latter entropic index before strong earthquakes arouse an intense international interest and has been studied by several workers [25,26,27,28,29,30,31,32].
Recapitulating the aforementioned recent studies [12] related with the M8.2 earthquake that occurred in Mexico on 7 September 2017, which is Mexico’s largest earthquake in more than a century, we can say that upon employing natural time analysis we found that almost three months before its occurrence the following precursory behavior was identified: The entropy change under time reversal exhibits a minimum [12] along with increased fluctuations of the entropy change under time reversal as well as by a simultaneous increase of the Tsallis entropic index q [18]. It is the main scope of this paper to investigate whether precursory behavior existed also before the super-giant M9 Tohoku earthquake that occurred in Japan on 11 March 2011. In addition, here we investigate whether the seminal work in 1961 by Lifshitz and Slyozov [33] and independently by Wagner [34], on phase transitions is applicable to this super-giant earthquake which can be considered, as mentioned above in the first paragraph, as a dynamic phase transition. We stress, however, that the present work differs essentially from our previous recent studies related with the M8.2 earthquake in Mexico, because—beyond the significant difference as far as their magnitudes is concerned—the following holds: These recent studies on Mexico’s earthquake have been carried out by analyzing in natural time the seismicity in Chiapas region during 2012–2017 (where the M8.2 earthquake occurred) since the aforementioned work in 2013 by Ramírez-Rojas and Flores-Márquez, which analyzed in natural time the seismicity in the six tectonic regions of Mexico, forwarded arguments that in the Chiapas region extreme events (i.e., large earthquakes) have been expected [13]. Unfortunately, such a work providing an estimate of the region where the M9 Tohoku earthquake was going to occur had not been published before its occurrence in 2011. (Such an estimate of the epicentral area was published recently [35] as explained below in Section 5.) In view of the lack of such an information at that time, natural time analysis is made here for the seismicity that occurred in the whole Japanese area N 25 46 E 125 148 for the period 1 January 1984–11 March 2011, the day of the M9 Tohoku earthquake. In other words, here we are going to answer the following question: If a super-giant earthquake is going to occur in a large area such as Japan, is it possible to identify a precursory behavior and estimate when the system approaches the critical point (earthquake occurrence)? This is answered in light of the fact that it is nowadays widely accepted, as for example stated by Holliday et al. [9], that the observed scaling laws indicate the existence of phenomena closely associated with the proximity of the system to a critical point as already mentioned. In particular, here we mainly focus our investigation on the complexity measure Λ i associated with the fluctuations of the entropy change under time reversal and on the Tsallis entropic index q which still continues to arouse international interest [36] by analyzing in natural time the seismicity occurring allover that area.
This paper is structured as follows: A summary of natural time analysis is given in Section 2, while the data along with the procedure for their analysis are presented in Section 3. Our results are described in Section 4 and a discussion follows in Section 5. Finally, our conclusions are summarized in Section 6.

2. Natural Time Analysis Background

Natural time analysis, introduced in the beginning of the 2000s [37,38,39,40,41], uncovers unique dynamic features hidden behind the time series of complex systems. For seismicity, in a time series comprising N earthquakes, the natural time χ k = k / N serves as an index for the occurrence of the kth earthquake. This index together with the energy Q k released during the kth earthquake of magnitude M k , i.e., the pair ( χ k , Q k ) , is studied in natural time analysis. Equivalently, one studies the pair ( χ k , p k ) , where
p k = Q k n = 1 N Q n
stands for the normalized energy released during the kth earthquake. The variance of χ weighted for p k , labeled by κ 1 , is given by [10,11,37,40,41,42]
κ 1 = k = 1 N p k ( χ k ) 2 k = 1 N p k χ k 2
where Q k , and hence p k , for earthquakes is estimated through the usual relation [43]
Q k 10 1.5 M k
In natural time analysis, an order parameter for seismicity has been proposed. In particular, it has been explained in Ref. [10] (see also pp. 249–253 of Ref. [11]) that the quantity κ 1 given by Equation (2) can be considered as an order parameter for seismicity since it changes abruptly when a mainshock (the new phase) occurs, and in addition the feature of its fluctuations resemble those in other non-equilibrium and equilibrium critical systems. Note that at least six earthquakes are needed for obtaining reliable κ 1 [10]. It has been found that the quantity κ 1 is a key parameter that enables recognition of the complex dynamical system under study entering the critical stage [11,37,38,39]. In short, Varotsos et al. [44] (see also p. 343 of Ref. [11]) showed that κ 1 becomes equal to 0.070 at the critical state for a variety of dynamical models. In Table 8.1 of Ref. [11], one can find a compilation of 14 cases in which the condition κ 1 = 0.070 has been ascertained. (This was also later confirmed in the analyses of very recent experimental results in Japan by Hayakawa and coworkers [45,46].) In addition, natural time has been recently employed by Turcotte and coworkers [47,48,49,50] as a basis for a new method to estimate the current level of seismic risk.
The entropy S in natural time is defined as the derivative with respect to d of the fluctuation function χ d χ d at d = 1 , which results in:
S χ ln χ χ ln χ
where f ( χ ) = k = 1 N p k f ( χ k ) . It is dynamic entropy and its value upon considering the time reversal T, i.e., T p k = p N k + 1 , is labeled by S . The value of S is, in general, different from S (e.g., [14]) (see also Ref. [51] and references therein), and thus S satisfies the conditions to be “causal”. The physical meaning of the change of entropy Δ S S S in natural time under time reversal has been discussed in Refs. [11,14,51]. The entropy S is a dynamic entropy, as mentioned, that exhibits [52] concavity, positivity, and Lesche stability [53,54].
Using a moving window of length i (number of consecutive events) sliding through the time series of L consecutive earthquakes 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 S i is obtained. By considering the standard deviation σ ( Δ S i ) of the time series of Δ S i S i ( S ) i , we define [11,18,55] 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 has been selected [18] to correspond to the standard deviation σ ( Δ S 100 ) of the time series of Δ S i of i = 100 events. This complexity measure quantifies how the statistics of Δ S i time series changes upon increasing the scale from 100 events to a longer scale, e.g., i = 10 3 events. The calculations are carried out by means of a window of length i (= number of successive earthquakes) sliding, each time by one earthquake, through the whole time series. The entropies S and S , and therefrom their difference Δ S i , are calculated each time, thus we also form a new time series consisting of successive Δ S i values and the complexity measure Λ i is determined according to its definition given in Equation (5).

3. Data and Analysis

The Japan Meteorological Agency (JMA) seismic catalogue was used (e.g., see Refs. [35,56]). We considered all earthquakes of magnitude M ≥ 3.5 from 1984 until the Tohoku earthquake occurrence on 11 March 2011 within the area 25 –46 N, 125 –148 E. The energy of earthquakes was obtained from the JMA magnitude M after converting [57] to the moment magnitude M w [43]. Setting a threshold M = 3.5 to assure data completeness, there exist 47,204 earthquakes in the area under discussion. Thus, we have on the average ∼150 earthquakes per month for the area considered.
The time evolution of the complexity measure Λ i is studied for a number of scales i of the seismicity with M ≥ 3.5 occurring in the whole area of Japan during the aforementioned 27-year period by choosing proper scales i as follows: We consider that recent investigations by means of natural time analysis showed that there exists the following interconnection between precursory low frequency (≤1 Hz) electric signals, termed Seismic Electric Signals (SES) (e.g., see Refs. [58,59]) and seismicity as follows [60]: The fluctuations β of the order parameter κ 1 of seismicity exhibit a minimum labeled β m i n when we observe the initiation of series of consecutive SES termed SES activities [51,61,62] exhibiting critical behavior [38,40,41] and have lead times ranging from a few weeks to around 5 1 2 months [11]. In addition, beyond this simultaneous appearance of SES activity and seismicity, Varotsos et al. [60] showed that these two phenomena are also linked closely in space, that gave the possibility of a reliable estimation of the epicentral area of an impending major earthquake. This has been subsequently confirmed in [35] for all major mainshocks of magnitude 7.6 or larger that occurred in Japan during 1984–2011 including the case of the M9 Tohoku earthquake. (Notably, before the latter earthquake, the minimum β m i n observed was the deepest during the whole period 1984–2011 investigated [56].) We return to this important point in Section 5. Before the initiation of the SES activity, and hence before β m i n , a stage has been detected in which the temporal correlations between earthquake magnitudes exhibit an anti-correlated behavior [63] while after its initiation long range correlations prevail between earthquake magnitudes. Thus, a significant change in the temporal correlations between earthquake magnitudes occurs when comparing the two stages that correspond to the periods before and just after the initiation of an SES activity. Since this change may be captured by the time evolution of Δ S i , we start our investigation of Δ S i from the scale of i 10 3 events, which is of the order of the number of seismic events M ≥ 3.5 that occur during a period around the maximum lead time of SES activities.

4. Results

We first present the results of our investigation concerning the complexity measure Λ i associated with the fluctuations of the entropy change under time reversal.
In Figure 1a–f, we plot the Λ i values for example for the time scales i = 2000, 3000 and 4000 events versus the conventional time from 1 January 1984 until the occurrence of the M9 Tohoku earthquake on 11 March 2011. In addition, all earthquakes of magnitude 7.0 or larger are also plotted in the same figure with vertical lines ending at circles read in the right scale. An inspection of this figure reveals that on 22 December 2010, i.e., almost two and a half months before the M9 earthquake an abrupt increase of the Λ i values for all the three scales is observed. This happens upon the occurrence of a M7.8 earthquake with an epicenter at 27.05 N 143.94 E [35]. This can be better visualized in Figure 2, which is a three-month excerpt of Figure 1f in expanded horizontal scale from 10 December 2010 to 11 May 2011 and depicts for the reader’s better inspection the increase Δ Λ i of the Λ i values after the occurrence of the aforementioned M7.8 earthquake on 22 December 2010.
A close inspection of Figure 2 shows that all three complexity measures Λ 2000 , Λ 3000 and Λ 4000 (cf. these symbols stand for the Λ i values at the scales i = 2000, 3000 and 4000, respectively) show a strong and abrupt increase on 22 December 2010 and after the occurrence of the aforementioned M7.8 earthquake exhibit a scaling behavior of the form
Δ Λ i = A ( t t 0 ) c
where the exponent c is independent of i with a value very close to 1 / 3 , while the pre-factors A are proportional to i (see Figure 3) and t 0 is approximately 0.2 days after the M7.8 earthquake occurrence. Equation (6) conforms to the seminal work by Lifshitz and Slyozov [33] and independently by Wagner [34] on phase transitions which shows that the time growth of the characteristic size of the minority phase droplets grows with time as t 1 / 3 .
To further elucidate the claim that an increase of Λ i is possibly associated with strong earthquakes, let us now further study the obvious increases (which are the most significant ones in Figure 1) of Λ i on 2 November 1989 and on 15 January 1993, as can be seen in Figure 1b, and investigate whether they exhibit a scaling behavior similar to that found above for the abrupt Λ i increase on 22 December 2010. Concerning the increase of Λ i on 2 November 1989, we give in Figure 4 the log-log plot of the changes Δ Λ i of the complexity measures Λ 2000 , Λ 3000 and Λ 4000 versus the elapsed time ( t t 0 ) in days from the occurrence time of a M7.1 earthquake on the same date increased by 0.022 days. Furthermore, for readers’ convenience, we also plot in Figure 4 the black straight line of slope 0.5 in order to easily visualize that the value of c for the three complexity measures Λ 2000 , Λ 3000 and Λ 4000 is close to 0.5, thus being distinctly different from the value c = 1 / 3 predicted by the LSW theory. In addition, an inspection of Figure 4 reveals that the pre-factors are not proportional to i (e.g., see that for i = 3000 events the green line is higher than the blue line that corresponds to i = 4000 events). In other words, the increase of Λ i associated with the earthquake on 2 November 1989 does not obey Equation (6), thus strongly deviating from LSW theory of phase transitions. By the same token as in Figure 4, we now give in Figure 5 the log-log plot that corresponds to the changes Δ Λ i of the complexity measures Λ 2000 , Λ 3000 and Λ 4000 versus the elapsed time ( t t 0 ) in days from the M7.5 earthquake occurrence on 15 January 1993 increased by 0.014 days. An inspection of this figure also reveals that the exponent in Figure 5 is close to 0.5—thus differing from the value c = 1 / 3 of LSW theory—and that the prefactors A are not proportional to the scale i, for example see the red crosses corresponding to the scale i = 2000 which do not practically differ from the green symbols of the larger scale i = 3000.
We now turn to the results obtained by means of non-extensive statistical mechanics [22], pioneered by Tsallis [19,24], which provides a framework for the study of complex systems in their non-equilibrium stationary states, systems with (multi)fractal and self-similar structures, long-range interacting systems, etc. This has found application [20,21] in the physics of earthquakes and especially in the description of the asperities in the faults on which earthquakes occur through the Tsallis entropic index q. Based on the earthquake magnitude distribution, one can deduce the [7,23] entropic index q and study how it varies with time as we approach a strong earthquake (for a recent review on this interesting aspect, see Ref. [32]). Figure 6 depicts the q-value versus conventional time during 1984–2011 as it is estimated [7] for several sliding windows of i = 1000, 2000, 3000, 4000 and 5000 consecutive earthquakes for M ≥ 3.5 in the Japanese area N 25 46 E 125 148 . We observe that, before the occurrence of the M9 Tohoku earthquake, the q-value exhibits an abrupt increase upon the occurrence of the M7.8 earthquake on 22 December 2010. A three-month excerpt of Figure 6 is given in expanded time scale in Figure 7 from 10 December 2010 until the Tohoku earthquake occurrence on 11 March 2011. By the same token as in Figure 3, Figure 4 and Figure 5, we now give in Figure 8 a log-log plot of the changes Δ q of the Tsallis entropic index versus the time ( t t 0 ) in days elapsed from the M7.8 earthquake occurrence on 22 December 2010 increased by approximately 0.2 days. An inspection of Figure 8 shows that the exponent c in Equation (6) is around 1 / 3 , as in LSW theory, but the prefactors A are not proportional to i.

5. Discussion

A close inspection of Figure 2 reveals the following interesting feature of the complexity measure Λ i : On 9 March 2011, a M7.3 foreshock occurred, and next day on 10 March 2011 all three complexity measures Λ 2000 , Λ 3000 and Λ 4000 exhibited a simultaneous decrease. What is the origin of this decrease, since next day, i.e., on 11 March 2011, the M9.0 mainshock occurred?
The following comments are now in order, based on a recent review [64] on the procedure on identifying the occurrence time of an impending major earthquake by means of natural time analysis. This can be achieved by analyzing in natural time the earthquakes in the candidate area. To apply this procedure, we need two important pieces of information. First, we need to know when we should set the natural time equal to zero and start the analysis. This is the time at which the system enters the critical stage. Second, we need an estimation of the candidate epicentral area. If geoelectrical measurements are available, both pieces of information become available upon the recording of an SES activity, because its initiation marks the time when the system enters the critical stage, and in addition the SES data provide an estimation of the epicentral area of the impending mainshock. On the other hand, if geoelectrical data are lacking, we make use of the following two recent findings by means of natural time analysis: first, the fluctuations β of the order parameter of seismicity in a large area exhibit a minimum a few months before a major earthquake almost simultaneously with the initiation of an SES activity [60]; and, second, a spatiotemporal study of this minimum unveils an estimate of the epicentral area of the impending major earthquake [35]. The application of this procedure to the determination of the occurrence time of the M9 Tohoku earthquake by Varotsos et al. [64] was made as follows:
As for the starting time of the natural time analysis of seismicity, they chose the date of 5 January 2011 which is the date of the appearance of the minimum of the fluctuations of the order parameter of seismicity before this major earthquake reported by Sarlis et al. [56]. This, which remarkably is the deepest minimum ever observed during the period investigated (1984–2011), almost coincides with the initiation of an SES activity since anomalous magnetic field variations appeared in the vertical component during the period 4–14 January 2011 at two measuring sites (Esashi and Mizusawa) lying at epicentral distances of around 130 km [65,66,67]. As for the estimation of the epicentral location of the impending mainshock without making use of SES data, Sarlis et al. [35] worked as follows: By dividing the entire Japanese region N 25 46 E 125 148 into small areas, a calculation of the fluctuations of κ 1 of seismicity is carried out on them. Some small areas show a minimum of the fluctuations almost simultaneously with the minimum in the entire Japanese region (on 5 January 2011) and such small areas cluster within a few hundred km from the actual epicenter, thus leading to an estimate of the candidate epicentral area. (Such an estimate of the epicentral area cannot be made of course in absence of the knowledge of the minimum of the fluctuations of the order parameter of seismicity on 5 January 2011. This absence happens as explained in the next paragraph if we do not consider the M7.8 earthquake on 22 December 2010.) A computation of the κ 1 values of seismicity in that area was made by starting from 5 January 2011. The results deduced by Varotsos et al. [64] are reproduced here in Figure 9 for M t h r e s = 4.2 to 5.0. Recalling that at least six earthquakes are needed [10] for obtaining reliable κ 1 value, which happens in the candidate epicentral area on 16 February 2011 for M t h r e s = 4.2, they depicted in their Figure 6a the computed κ 1 values during the last four weeks before the M9 Tohoku earthquake occurrence. This figure clearly showed that the condition κ 1 = 0.070 was not satisfied for all magnitude thresholds at least until the M7.3 earthquake on 9 March 2011. Here, we plot in expanded time scale in Figure 9 the κ 1 values of seismicity from 00:00 LT on 9 March 2011 until the Tohoku earthquake occurrence. This figure reveals that the condition κ 1 = 0.070 (which signals that the mainshock is going to occur within the next few days or so) is fulfilled for M t h r e s = 4.2–5.0 in the morning of 10 March 2011 upon the occurrence of the earthquakes from 08:36 to 13:14 LT, i.e., almost one day before the Tohoku earthquake, see the gray shaded area in Figure 9. (Such a determination of the occurrence time of Tohoku earthquake cannot be achieved of course in absence of the knowledge of the date of the minimum of the fluctuations of the order parameter of seismicity on 5 January 2011, because from this date we started the computation as mentioned. This absence happens if we do not consider the M7.8 earthquake on 22 December 2010, as is explained in the next paragraph.) In this Figure, we have also inserted the values of the change Δ Λ 2000 , Δ Λ 3000 and Δ Λ 4000 (depicted in Figure 2) of the complexity measures but multiplied by a factor of four to better visualize their variations. Quite interestingly, we now clearly see that they exhibit a decrease just after the gray shaded area in Figure 9 where the condition κ 1 = 0.070 was fulfilled, thus signaling that the mainshock was approaching.
We now investigate what happens with the observation of the minimum of the fluctuations of the order parameter of seismicity if we do not consider the entire Japanese region N 25 46 E 125 148 and select for example the area N 28 46 E 125 148 which does not include the epicenter (27.05°N,143.94°E) of the M7.8 earthquake on 22 December 2010. Figure 10 depicts the fluctuations β 200 of the order parameter of seismicity when a natural time window comprising 200 events (M 3.5 ) is sliding through the JMA catalog of the area N 28 46 E 125 148 from 1 January 1984 to the M9 Tohoku earthquake occurrence. An inspection of Figure 10 shows that the deepest minimum of β 200 during this period is not observed on 5 January 2010, as found in Refs. [42,56,63], when the investigation is extended to the entire Japanese region N 25 46 E 125 148 . The absence of this minimum reflects that neither the epicentral area nor the occurrence time of the M9 Tohoku earthquake can be estimated if we do not consider the M7.8 earthquake on 22 December 2010.

6. Main Conclusions

Almost two and a half months before the M9 Tohoku earthquake occurrence, i.e., on 22 December 2010, the following two facts were observed:
First, the complexity measure Λ i associated with the fluctuations of the entropy change of seismicity in natural time under time reversal exhibited an abrupt increase, which conforms to the Lifshitz–Slyozov–Wagner theory for phase transitions showing that the characteristic size of the minority phase droplets exhibits a scaling behavior in which time growth has the form A ( t t 0 ) 1 / 3 where the prefactors A are proportional to the scale i, while the exponent (1/3) is independent of i.
Second, the Tsallis entropic index q shows a simultaneous increase which interestingly exhibits the same exponent (1/3) (but the prefactors A are not proportional to the scale i).
Beyond the aforementioned two facts, we also found that the complexity measures Λ 2000 , Λ 3000 and Λ 4000 exhibit a decrease just after the fulfillment of the condition κ 1 = 0.070 almost one day before the Tohoku earthquake occurrence.

Author Contributions

All authors contributed equally to this research.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Telesca, L.; Lapenna, V.; Macchiato, M. Spatial variability of the time-correlated behaviour in Italian seismicity. Earth Planetary Sci. Let. 2003, 212, 279–290. [Google Scholar] [CrossRef]
  2. Huang, Q. Seismicity changes prior to the Ms8.0 Wenchuan earthquake in Sichuan, China. Geophys. Res. Lett. 2008, 35, L23308. [Google Scholar] [CrossRef]
  3. Lennartz, S.; Livina, V.N.; Bunde, A.; Havlin, S. Long-term memory in earthquakes and the distribution of interoccurrence times. EPL 2008, 81, 69001. [Google Scholar] [CrossRef]
  4. Telesca, L.; Lovallo, M. Non-uniform scaling features in central Italy seismicity: A non-linear approach in investigating seismic patterns and detection of possible earthquake precursors. Geophys. Res. Lett. 2009, 36, L01308. [Google Scholar] [CrossRef]
  5. Lennartz, S.; Bunde, A.; Turcotte, D.L. Modelling seismic catalogues by cascade models: Do we need long-term magnitude correlations? Geophys. J. Int. 2011, 184, 1214–1222. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, Q. Retrospective investigation of geophysical data possibly associated with the Ms8.0 Wenchuan earthquake in Sichuan, China. J. Asian Earth Sci. 2011, 41, 421–427. [Google Scholar] [CrossRef]
  7. Telesca, L. Maximum Likelihood Estimation of the Nonextensive Parameters of the Earthquake Cumulative Magnitude Distribution. Bull. Seismol. Soc. Am. 2012, 102, 886–891. [Google Scholar] [CrossRef]
  8. Turcotte, D.L. Fractals and Chaos in Geology and Geophysics, 2nd ed.; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  9. Holliday, J.R.; Rundle, J.B.; Turcotte, D.L.; Klein, W.; Tiampo, K.F.; Donnellan, A. Space-Time Clustering and Correlations of Major Earthquakes. Phys. Rev. Lett. 2006, 97, 238501. [Google Scholar] [CrossRef] [PubMed]
  10. Varotsos, P.A.; Sarlis, N.V.; Tanaka, H.K.; Skordas, E.S. Similarity of fluctuations in correlated systems: The case of seismicity. Phys. Rev. E 2005, 72, 041103. [Google Scholar] [CrossRef] [PubMed]
  11. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Natural Time Analysis: The New View of Time. Precursory Seismic Electric Signals, Earthquakes and other Complex Time-Series; Springer-Verlag: Berlin/Heidelberg, Germany, 2011. [Google Scholar] [CrossRef]
  12. 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]
  13. Ramírez-Rojas, A.; Flores-Márquez, E. Order parameter analysis of seismicity of the Mexican Pacific coast. Phys. A 2013, 392, 2507–2512. [Google Scholar] [CrossRef]
  14. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S.; Lazaridou, M.S. Identifying sudden cardiac death risk and specifying its occurrence time by analyzing electrocardiograms in natural time. Appl. Phys. Lett. 2007, 91, 064106. [Google Scholar] [CrossRef]
  15. 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]
  16. 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]
  17. 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]
  18. 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 Mexico M8.2 Earthquake on 7 September 2017. Entropy 2018, 20, 477. [Google Scholar] [CrossRef]
  19. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  20. Sotolongo-Costa, O.; Posadas, A. Fragment-Asperity Interaction Model for Earthquakes. Phys. Rev. Lett. 2004, 92, 048501. [Google Scholar] [CrossRef] [PubMed]
  21. 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]
  22. Tsallis, C. Introduction to Nonextensive Statistical Mechanics; Springer: Berlin, Germany, 2009. [Google Scholar] [CrossRef]
  23. 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]
  24. 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]
  25. Telesca, L. A non-extensive approach in investigating the seismicity of L’Aquila area (central Italy), struck by the 6 April 2009 earthquake (ML = 5.8). Terra Nova 2010, 22, 87–93. [Google Scholar] [CrossRef]
  26. Telesca, L. Tsallis-Based Nonextensive Analysis of the Southern California Seismicity. Entropy 2011, 13, 1267–1280. [Google Scholar] [CrossRef] [Green Version]
  27. Vallianatos, F.; Sammonds, P. Evidence of non-extensive statistical physics of the lithospheric instability approaching the 2004 Sumatran-Andaman and 2011 Honshu mega-earthquakes. Tectonophysics 2013, 590, 52–58. [Google Scholar] [CrossRef]
  28. Vallianatos, F.; Michas, G.; Papadakis, G. Non-extensive and natural time analysis of seismicity before the Mw6.4, October 12, 2013 earthquake in the South West segment of the Hellenic Arc. Phys. A Stat. Mech. Appl. 2014, 414, 163–173. [Google Scholar] [CrossRef]
  29. Papadakis, G.; Vallianatos, F.; Sammonds, P. A Nonextensive Statistical Physics Analysis of the 1995 Kobe, Japan Earthquake. Pure Appl. Geophys. 2015, 172, 1923–1931. [Google Scholar] [CrossRef]
  30. Papadakis, G.; Vallianatos, F.; Sammonds, P. Non-extensive statistical physics applied to heat flow and the earthquake frequency-magnitude distribution in Greece. Physica A: Statistical Mechanics and its Applications 2016, 456, 135–144. [Google Scholar] [CrossRef]
  31. Papadakis, G.; Vallianatos, F. Non-extensive statistical physics analysis of earthquake magnitude sequences in North Aegean Trough, Greece. Acta Geophys. 2017, 65, 555–563. [Google Scholar] [CrossRef]
  32. Vallianatos, F.; Michas, G.; Papadakis, G. Nonextensive Statistical Seismology: An Overview. In Complexity of Seismic Time Series: Measurement and Application; Elsevier Science: Amsterdam, The Netherlands, 2018; pp. 25–59. [Google Scholar] [CrossRef]
  33. Lifshitz, I.; Slyozov, V. The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. Solids 1961, 19, 35–50. [Google Scholar] [CrossRef]
  34. Wagner, C. Theorie der Alterung von Niederschlägen durch Umlösen (Ostwald-Reifung). Zeitschrift Elektrochemie Berichte Bunsengesellschaft Physikalische Chemie 1961, 65, 581–591. [Google Scholar] [CrossRef]
  35. Sarlis, N.V.; Skordas, E.S.; Varotsos, P.A.; Nagao, T.; Kamogawa, M.; Uyeda, S. Spatiotemporal variations of seismicity before major earthquakes in the Japanese area and their relation with the epicentral locations. Proc. Natl. Acad. Sci. USA 2015, 112, 986–989. [Google Scholar] [CrossRef] [PubMed]
  36. Livadiotis, G. Thermodynamic origin of kappa distributions. EPL (Europhys. Lett.) 2018, 122, 50001. [Google Scholar] [CrossRef]
  37. 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]
  38. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Long-range correlations in the electric signals that precede rupture. Phys. Rev. E 2002, 66, 011902. [Google Scholar] [CrossRef] [PubMed]
  39. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Seismic Electric Signals and Seismicity: On a tentative interrelation between their spectral content. Acta Geophys. Pol. 2002, 50, 337–354. [Google Scholar]
  40. 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]
  41. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Long-range correlations in the electric signals the precede rupture: Further investigations. Phys. Rev. E 2003, 67, 021109. [Google Scholar] [CrossRef] [PubMed]
  42. Varotsos, P.; Sarlis, N.; Skordas, E. Scale-specific order parameter fluctuations of seismicity in natural time before mainshocks. EPL (Europhys. Lett.) 2011, 96, 59002. [Google Scholar] [CrossRef]
  43. Kanamori, H. Quantification of Earthquakes. Nature 1978, 271, 411–414. [Google Scholar] [CrossRef]
  44. 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]
  45. Potirakis, S.M.; Asano, T.; Hayakawa, M. 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. Entropy 2018, 20, 199. [Google Scholar] [CrossRef]
  46. Potirakis, S.M.; Schekotov, A.; Asano, T.; Hayakawa, M. Natural time analysis on the ultra-low frequency magnetic field variations prior to the 2016 Kumamoto (Japan) earthquakes. J. Asian Earth Sci. 2018, 154, 419–427. [Google Scholar] [CrossRef]
  47. Rundle, J.B.; Turcotte, D.L.; Donnellan, A.; Grant Ludwig, L.; Luginbuhl, M.; Gong, G. Nowcasting earthquakes. Earth Space Sci. 2016, 3, 480–486. [Google Scholar] [CrossRef] [Green Version]
  48. Rundle, J.B.; Luginbuhl, M.; Giguere, A.; Turcotte, D.L. Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities. Pure Appl. Geophys. 2018, 175, 647–660. [Google Scholar] [CrossRef]
  49. Luginbuhl, M.; Rundle, J.B.; Hawkins, A.; Turcotte, D.L. Nowcasting Earthquakes: A Comparison of Induced Earthquakes in Oklahoma and at the Geysers, California. Pure Appl. Geophys. 2018, 175, 49–65. [Google Scholar] [CrossRef]
  50. Luginbuhl, M.; Rundle, J.B.; Turcotte, D.L. Natural Time and Nowcasting Earthquakes: Are Large Global Earthquakes Temporally Clustered? Pure Appl. Geophys. 2018, 175, 661–670. [Google Scholar] [CrossRef]
  51. 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]
  52. 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]
  53. Lesche, B. Instabilities of Renyi entropies. J. Stat. Phys. 1982, 27, 419. [Google Scholar] [CrossRef]
  54. Lesche, B. Renyi entropies and observables. Phys. Rev. E 2004, 70, 017102. [Google Scholar] [CrossRef] [PubMed]
  55. Sarlis, N.V.; Christopoulos, S.R.G.; Bemplidaki, M.M. Change ΔS of the entropy in natural time under time reversal: Complexity measures upon change of scale. EPL (Europhy. Lett.) 2015, 109, 18002. [Google Scholar] [CrossRef]
  56. Sarlis, N.V.; Skordas, E.S.; Varotsos, P.A.; Nagao, T.; Kamogawa, M.; Tanaka, H.; Uyeda, S. Minimum of the order parameter fluctuations of seismicity before major earthquakes in Japan. Proc. Natl. Acad. Sci. USA 2013, 110, 13734–13738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Tanaka, H.K.; Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. A plausible universal behaviour of earthquakes in the natural time-domain. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2004, 80, 283–289. [Google Scholar] [CrossRef] [Green Version]
  58. Varotsos, P.; Alexopoulos, K. Physical Properties of the variations of the electric field of the earth preceding earthquakes, I. Tectonophysics 1984, 110, 73–98. [Google Scholar] [CrossRef]
  59. Varotsos, P.; Alexopoulos, K. Physical Properties of the variations of the electric field of the earth preceding earthquakes, II. Tectonophysics 1984, 110, 99–125. [Google Scholar] [CrossRef]
  60. 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]
  61. Varotsos, P.; Lazaridou, M. Latest aspects of earthquake prediction in Greece based on Seismic Electric Signals. Tectonophysics 1991, 188, 321–347. [Google Scholar] [CrossRef]
  62. Varotsos, P.; Alexopoulos, K.; Lazaridou, M. Latest aspects of earthquake prediction in Greece based on Seismic Electric Signals, II. Tectonophysics 1993, 224, 1–37. [Google Scholar] [CrossRef]
  63. 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. Space Phys. 2014, 119, 9192–9206. [Google Scholar] [CrossRef] [Green Version]
  64. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Identifying the occurrence time of an impending major earthquake: A review. Earthq. Sci. 2017, 30, 209–218. [Google Scholar] [CrossRef]
  65. Xu, G.; Han, P.; Huang, Q.; Hattori, K.; Febriani, F.; Yamaguchi, H. Anomalous behaviors of geomagnetic diurnal variations prior to the 2011 off the Pacific coast of Tohoku earthquake (Mw9.0). J. Asian Earth Sci. 2013, 77, 59–65. [Google Scholar] [CrossRef]
  66. Han, P.; Hattori, K.; Xu, G.; Ashida, R.; Chen, C.H.; Febriani, F.; Yamaguchi, H. Further investigations of geomagnetic diurnal variations associated with the 2011 off the Pacific coast of Tohoku earthquake (Mw 9.0). J. Asian Earth Sci. 2015, 114, 321–326. [Google Scholar] [CrossRef]
  67. Han, P.; Hattori, K.; Huang, Q.; Hirooka, S.; Yoshino, C. Spatiotemporal characteristics of the geomagnetic diurnal variation anomalies prior to the 2011 Tohoku earthquake (Mw 9.0) and the possible coupling of multiple pre-earthquake phenomena. J. Asian Earth Sci. 2016, 129, 13–21. [Google Scholar] [CrossRef]
Figure 1. (af) Plot of the complexity measure Λ i versus the conventional time for the scales i = 2000 (red), 3000 (blue) and 4000 events (green) from 1 January 1984 until the M9 Tohoku earthquake on 11 March 2011. The vertical lines ending at circles depict the magnitudes (M ≥ 7) of earthquakes read in the right scale.
Figure 1. (af) Plot of the complexity measure Λ i versus the conventional time for the scales i = 2000 (red), 3000 (blue) and 4000 events (green) from 1 January 1984 until the M9 Tohoku earthquake on 11 March 2011. The vertical lines ending at circles depict the magnitudes (M ≥ 7) of earthquakes read in the right scale.
Entropy 20 00757 g001
Figure 2. An almost three-month excerpt of Figure 1 in expanded time scale which shows the change Δ Λ i of Λ i values versus the conventional time after the occurrence of the M7.8 earthquake on 22 December 2010 with an epicenter at 27.05 N 143.94 E. Note that, after the M7.3 foreshock that occurred on 9 March 2011, a decrease of the Δ Λ i appears.
Figure 2. An almost three-month excerpt of Figure 1 in expanded time scale which shows the change Δ Λ i of Λ i values versus the conventional time after the occurrence of the M7.8 earthquake on 22 December 2010 with an epicenter at 27.05 N 143.94 E. Note that, after the M7.3 foreshock that occurred on 9 March 2011, a decrease of the Δ Λ i appears.
Entropy 20 00757 g002
Figure 3. Log-log plot of the change Δ Λ i of the complexity measures Λ 2000 (red), Λ 3000 (green) and Λ 4000 (blue) versus the elapsed time ( t t 0 ) in days since the establishment of scaling behavior after the occurrence of the M7.8 earthquake on 22 December 2010. The value of t 0 is approximately 0.2 days measured from the M7.8 earthquake occurrence and c = 1 / 3 .
Figure 3. Log-log plot of the change Δ Λ i of the complexity measures Λ 2000 (red), Λ 3000 (green) and Λ 4000 (blue) versus the elapsed time ( t t 0 ) in days since the establishment of scaling behavior after the occurrence of the M7.8 earthquake on 22 December 2010. The value of t 0 is approximately 0.2 days measured from the M7.8 earthquake occurrence and c = 1 / 3 .
Entropy 20 00757 g003
Figure 4. The same as Figure 3, but for the M7.1 earthquake on 2 November 1989 with an epicenter at 39.86 N 143.05 E. The value of t 0 is 0.022 days measured from the M7.1 earthquake occurrence and c = 0.5 .
Figure 4. The same as Figure 3, but for the M7.1 earthquake on 2 November 1989 with an epicenter at 39.86 N 143.05 E. The value of t 0 is 0.022 days measured from the M7.1 earthquake occurrence and c = 0.5 .
Entropy 20 00757 g004
Figure 5. The same as Figure 3, but for the M7.5 earthquake on 15 January 1993 with an epicenter at 42.92 N 144.35 E. The value of t 0 is 0.014 days measured from the M7.5 earthquake occurrence and c = 0.5 .
Figure 5. The same as Figure 3, but for the M7.5 earthquake on 15 January 1993 with an epicenter at 42.92 N 144.35 E. The value of t 0 is 0.014 days measured from the M7.5 earthquake occurrence and c = 0.5 .
Entropy 20 00757 g005
Figure 6. (af) Plots of the values of the Tsallis entropic index q at several scales i = 1000, 2000, 3000, 4000 and 5000; events as shown by the colors in the inset.
Figure 6. (af) Plots of the values of the Tsallis entropic index q at several scales i = 1000, 2000, 3000, 4000 and 5000; events as shown by the colors in the inset.
Entropy 20 00757 g006
Figure 7. An almost three-month excerpt of Figure 6 after the occurrence of the M7.8 earthquake on 22 December 2010.
Figure 7. An almost three-month excerpt of Figure 6 after the occurrence of the M7.8 earthquake on 22 December 2010.
Entropy 20 00757 g007
Figure 8. Log-log plot of the change Δ q of the values of the Tsallis entropic index q versus the elapsed time ( t t 0 ) in days since the establishment of scaling behavior after the occurrence of the M7.8 earthquake on 22 December 2010. The value of t 0 is 0.2 days measured from the M7.8 earthquake occurrence and c = 1 / 3 .
Figure 8. Log-log plot of the change Δ q of the values of the Tsallis entropic index q versus the elapsed time ( t t 0 ) in days since the establishment of scaling behavior after the occurrence of the M7.8 earthquake on 22 December 2010. The value of t 0 is 0.2 days measured from the M7.8 earthquake occurrence and c = 1 / 3 .
Entropy 20 00757 g008
Figure 9. The κ 1 values as well as the values of the change Δ Λ 2000 , Δ Λ 3000 and Δ Λ 4000 of the complexity measures versus the conventional time since 00:00 LT on 9 March 2011 until the M9 Tohoku earthquake occurrence. The shaded area marks the period in the morning of 10 March 2011 during which the condition κ 1 = 0.070 is fulfilled.
Figure 9. The κ 1 values as well as the values of the change Δ Λ 2000 , Δ Λ 3000 and Δ Λ 4000 of the complexity measures versus the conventional time since 00:00 LT on 9 March 2011 until the M9 Tohoku earthquake occurrence. The shaded area marks the period in the morning of 10 March 2011 during which the condition κ 1 = 0.070 is fulfilled.
Entropy 20 00757 g009
Figure 10. The fluctuations β 200 of the order parameter of seismicity when a window comprising 200 events is sliding through the JMA catalog (M 3.5 ) in the area N 28 46 E 125 148 which does not contain the epicenter (27.05 N,143.94 E) of the M7.8 earthquake on 22 December 2010 for the periods: (a) 1 January 1984–1 January 2000; and (b) 1 January 2000 until the M9 Tohoku earthquake occurrence.
Figure 10. The fluctuations β 200 of the order parameter of seismicity when a window comprising 200 events is sliding through the JMA catalog (M 3.5 ) in the area N 28 46 E 125 148 which does not contain the epicenter (27.05 N,143.94 E) of the M7.8 earthquake on 22 December 2010 for the periods: (a) 1 January 1984–1 January 2000; and (b) 1 January 2000 until the M9 Tohoku earthquake occurrence.
Entropy 20 00757 g010

Share and Cite

MDPI and ACS Style

Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Tsallis Entropy Index q and the Complexity Measure of Seismicity in Natural Time under Time Reversal before the M9 Tohoku Earthquake in 2011. Entropy 2018, 20, 757. https://doi.org/10.3390/e20100757

AMA Style

Varotsos PA, Sarlis NV, Skordas ES. Tsallis Entropy Index q and the Complexity Measure of Seismicity in Natural Time under Time Reversal before the M9 Tohoku Earthquake in 2011. Entropy. 2018; 20(10):757. https://doi.org/10.3390/e20100757

Chicago/Turabian Style

Varotsos, Panayiotis A., Nicholas V. Sarlis, and Efthimios S. Skordas. 2018. "Tsallis Entropy Index q and the Complexity Measure of Seismicity in Natural Time under Time Reversal before the M9 Tohoku Earthquake in 2011" Entropy 20, no. 10: 757. https://doi.org/10.3390/e20100757

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