Journal cover Journal topic
Annales Geophysicae An interactive open-access journal of the European Geosciences Union
Journal topic
Ann. Geophys., 37, 315-324, 2019
https://doi.org/10.5194/angeo-37-315-2019

Special issue: Satellite observations for space weather and geo-hazard

Ann. Geophys., 37, 315-324, 2019
https://doi.org/10.5194/angeo-37-315-2019

Regular paper 10 May 2019

Regular paper | 10 May 2019

# Phenomena preceding major earthquakes interconnected through a physical model

Phenomena preceding major earthquakes
Panayiotis A. Varotsos1,2, Nicholas V. Sarlis1,2, and Efthimios S. Skordas1,2 Panayiotis A. Varotsos et al.
• 1Department of Physics, Section of Solid State Physics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos 157 84, Athens, Greece
• 2Department of Physics, Solid Earth Physics Institute, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos 157 84, Athens, Greece
Abstract

The analysis of earthquake time series in a new time domain termed natural time enables the uncovering of hidden properties in time series of complex systems and has been recently employed as the basis of a method to estimate seismic risk. Natural time also enables the determination of the order parameter of seismicity, which is a quantity by means of which one can identify when the system approaches the critical point (the mainshock occurrence is considered the new phase). Applying this analysis, as an example, to the Japanese seismic data from 1 January 1984 until the super-giant M 9 Tōhoku earthquake on 11 March 2011, we find that almost 3 months before its occurrence the entropy change of seismicity under time reversal is minimized on 22 December 2010, which signals an impending major earthquake. On this date the order parameter fluctuations of seismicity exhibit an abrupt increase. This increase is accompanied by various phenomena; e.g., from this date the horizontal GPS azimuths start to become gradually oriented toward the southern direction, while they had random orientation during the preceding period. Two weeks later, a minimum of the order parameter fluctuations of seismicity appears accompanied by anomalous Earth magnetic field variations and by full alignment of the orientations of GPS azimuths southwards leading to the most intense crust uplift. These phenomena are discussed and found to be in accordance with a physical model which seems to explain on a unified basis anomalous precursory changes observed either in ground-based measurements or in satellite data.

1 Introduction

Almost 8 years have passed since the Tōhoku earthquake (EQ) that occurred on 11 March 2011 in Japan with magnitude (M) 9.0, the largest magnitude ever recorded in Japan. In the meantime, independent research groups reported anomalous precursory variations of quantities like the geomagnetic field , seismicity upon analyzing it in natural time , and Earth's surface displacements measured by Global Positioning System (GPS) . In the 1980s, a short-term earthquake prediction method was introduced based on the observation of seismic electric signals (SESs), which are low-frequency transient changes in the electric field of the Earth preceding EQs . Several SESs recorded within a short time are termed SES activity . Major EQs are preceded by intense SES activities accompanied by evident Earth magnetic field variations mainly recorded on the z component (; ). This method was motivated by a physical model for SES generation, which also foresees that some additional transient multidisciplinary phenomena should be simultaneously generated and ended well before the EQ rupture as schematically shown in Fig. 1a. This was in direct contrast to other proposed precursory mechanisms (Rikitake1981) that usually exhibit anomalous behavior becoming more intense upon approaching the EQ failure as seen in Fig. 1b. The scope of the present study is twofold: first, investigate whether the transient phenomena foreseen by this model actually appeared before the M 9 Tōhoku EQ, having been observed by ground-based measurements or by GPS data; second, report other precursory phenomena that appeared before the Tōhoku EQ almost simultaneously with the transient phenomena that were expected on the basis of this model.

Figure 1Schematic diagram showing the two distinct approaches proposed for a precursory behavior. (a) The case of the physical model for SES generation, which differs greatly from other suggested mechanisms (Rikitake1981) in which the anomalous precursory behavior becomes more intense upon approaching the EQ occurrence (b).

The aforementioned physical model for SES generation, termed the pressure stimulated polarization current (PSPC) model (; summarized also in , and ), suggests the following (see Fig. 2): in the Earth, electric dipoles always exist due to lattice imperfections (point and linear defects; e.g., see ) in the ionic constituents of rocks. In the future focal region of an EQ, where the electric dipoles have initially random orientations (Fig. 2c), the stress, σ, starts to gradually increase due to an excess stress disturbance (Fig. 2a). Let us call this stage A hereafter. When this gradually increasing stress reaches a critical value (σcr), the electric dipoles exhibit a cooperative orientation (Fig. 2e) resulting in the emission of a transient SES (Fig. 2b) with current density j. We call this stage B. pointed out that the PSPC model is unique among other models in that SESs would be generated spontaneously during the gradual increase in stress without requiring any sudden change in stress such as microfracturing (or faulting).

Figure 2Schematic diagram of the physical model proposed for the SES generation. (a) Before an EQ, the stress σ gradually increases in the focal area versus the time t towards reaching a critical value σcr. (b) When σ reaches σcr a transient electric signal is emitted that constitutes an SES. (c) Random orientation of the electric dipoles at small stress. (d) Partial orientation at an intermediate stress σ(<σcr). (e) Cooperative orientation of the electric dipoles when σ = σcr.

Observations of SES activities in Japan , in China (see , and references therein; for example see the geoelectric field changes depicted in Fig. 2b of , that started almost 50 d before the Ms 8.0 Wenchuan EQ in 2008), and in Mexico (see p. 220 of Ramírez-Rojas et al.2011) and in California (where magnetic field variations similar to those associated with the SES activities in Greece have been reported; e.g., see ; see also ) have shown that their lead time lies in the range from a few weeks to a few months or so, in agreement with earlier observations in Greece (where the maximum lead time observed is 5.5 months or so; e.g., see ). Hence, the SES observations in various EQ-prone areas reveal that before the occurrence of major EQs there is a crucial timescale (from a few weeks to around a few months or so), in which the critical stress σcr is attained and changes in other associated physical quantities should become also detectable .

This paper is structured as follows: in the next section, i.e., Sect. 2, we present the anomalous variations of multidisciplinary nature observed before the M 9 Tōhoku EQ by independent research groups, while our own findings obtained by natural time analysis of the seismicity of Japan are given in Sect. 3. In the subsequent Sect. 4, we investigate whether the observed precursory variations are in accordance with the PSPC model, and in the final section, i.e., Sect. 5, we summarize our conclusions.

2 Anomalous variations of multidisciplinary nature observed by independent groups before the 2011 M 9 Tōhoku EQ

## 2.1 Earth's magnetic field variations

found anomalous behavior of geomagnetic diurnal variations mainly in the vertical component at the Esashi station (ESA) located at about 135 km from the M 9 Tōhoku EQ focal zone for about 10 d, i.e., 4–14 January 2011. They analyzed geomagnetic data of a 3-year period, i.e., from 1 January 2010 to 31 December 2012, by computing ratios of diurnal variation range between the target station ESA and the remote reference station Kakioka about 300 km from the EQ epicenter. To validate this finding, further investigations were reported by after analyzing geomagnetic data of 16-year long-term observations in Japanese stations. They again found that the mean values of the ratios of the diurnal variations in the vertical component showed a clear anomaly exceeding the statistical threshold during the aforementioned period, i.e., 4–14 January 2011, and in addition they emphasized that this anomaly is unique in over 16 years. This has been further validated by the most recent study of , who analyzed geomagnetic data of long-term observations at 17 stations in Japan. They found that the above unique anomaly in the vertical component has also been observed at a second station at Mizusawa (MIZ) in the Tōhoku region, which is about 20 km southwest of the ESA station. This fact that both ESA and MIZ show clear anomalies at the same time suggests that the anomaly cannot be the result of observation system error or artificial noises . Furthermore, this anomaly cannot be attributed to magnetic storms since it has been observed during a period in which no moderate–strong magnetic storms were recorded.

## 2.2 Earth's surface displacements

Daily resolution data retrieved from the 1243 GPS stations in Japan were utilized by to expose surface displacements before the M 9 Tōhoku EQ. They applied the method proposed by on filtering long-term plate movements, short-term noise, and frequency-dependent (i.e., semiannual and annual) variations from the three-component GPS data for all stations. The N–S and E–W components were utilized to compute the orientations of the horizontal azimuths, termed GPS azimuths.

Figure 3Schematic diagram that compiles the multidisciplinary changes before the Tōhoku EQ. Period of observations: from 12 December 2010 until 23 January 2011. The black text describes the observations by means of residual GPS data, the brown text describes those obtained by means of Earth's magnetic field variations, the green text describes those obtained by detrended fluctuation analysis (DFA), the blue text describes those obtained by natural time analysis, and the purple text describes those obtained by anomalous groundwater and radon changes.

In general, the residual surface displacements are random . However, as depicted by black text in Fig. 3, southward movements became evident on 5 January 2011, i.e., 65 d before the Tōhoku EQ. Other changes before and after 5 January 2011 have also been observed as follows: while during the period 12–22 December 2010 random orientations of GPS azimuths prevailed, a gradual alignment toward the southern direction started on 22 December 2010 and continued until around 5 January 2011, accompanied with a gradual uplift of the crust. The most intense crust uplift was observed approximately on 5 January 2011 together with the full alignment of GPS azimuths southwards (more details are given in Sect. 4).

## 2.3 Changes in the level and temperature of confined groundwater

reported that anomalous groundwater changes started 3 months before the Tōhoku EQ. In particular, groundwater level and temperature decreased almost simultaneously in a 2000 m well at a spa, Goyo-onsen, in Iwate Prefecture, 155 km northwest of the epicenter. This simultaneous decrease occurred only once in the 3.5-year records when considering that the recordings started in this source since October 2007. emphasized that Tōhoku EQ is the only EQ that was preceded by anomalous changes in both water level and temperature. The exact date of the initiation of this phenomenon is not mentioned explicitly by , because the measurements were not made continuously but were taken intermittently and irregularly (the average interval between them being 8 d). They plotted, however, the consecutive measurements versus conventional time in their Fig. 1 (period 2007–2012) in which one can read that the phenomenon initiated around 20 December 2010, which agrees with what they state, i.e., around 3 months before Tōhoku EQ.

In addition, reported that, according to radon concentration measurements in the groundwater in the Izu Peninsula (at a distance about 500 km from the epicenter), an increase started almost 3 months before Tōhoku EQ. This increase occurred only this time during a 35-year observation.

3 Precursory changes observed by means of natural time analysis of Japanese seismicity

Natural time analysis uncovers important hidden properties in time series of complex systems and has been recently employed by Turcotte and coworkers as the basis of a new methodology (nowcasting) to estimate the current seismic risk level .

## 3.1 Natural time analysis background

In a time series comprising N EQs, the natural time for the occurrence of the kth EQ of energy Qk is defined as ${\mathit{\chi }}_{k}=k/N$. In natural time analysis, we study the evolution of the pair (χk,pk), where

$\begin{array}{}\text{(1)}& {p}_{k}={Q}_{k}/\sum _{n=\mathrm{1}}^{N}{Q}_{n}\end{array}$

denotes the normalized energy released during the kth EQ. Qk and hence pk for earthquakes are estimated through the relation (Kanamori1978)

$\begin{array}{}\text{(2)}& {Q}_{k}\propto {\mathrm{10}}^{\mathrm{1.5}{M}_{k}}.\end{array}$

It is widely accepted that the observed earthquake scaling laws indicate the existence of phenomena closely associated with the proximity of the system to a critical point (the mainshock is the new phase). In particular, it has been indicated by that it seems possible that systems that operate persistently near a threshold of instability are in some way like thermodynamic systems near critical points (EQ can be regarded as a stick-slip frictional instability of a preexisting fault). The order parameter of seismicity is the quantity by which one can identify the approach of the dynamical system to a critical point (the mainshock is the new phase). It was argued by that the variance

$\begin{array}{}\text{(3)}& {\mathit{\kappa }}_{\mathrm{1}}=〈{\mathit{\chi }}^{\mathrm{2}}〉-〈\mathit{\chi }{〉}^{\mathrm{2}}\end{array}$

of natural time χ weighted for pk, given by

$\begin{array}{}\text{(4)}& {\mathit{\kappa }}_{\mathrm{1}}=\sum _{k=\mathrm{1}}^{N}{p}_{k}\left({\mathit{\chi }}_{k}{\right)}^{\mathrm{2}}-{\left(\sum _{k=\mathrm{1}}^{N}{p}_{k}{\mathit{\chi }}_{k}\right)}^{\mathrm{2}},\end{array}$

may serve as an order parameter of seismicity.

The entropy S in natural time is defined by

$\begin{array}{}\text{(5)}& S\equiv 〈\mathit{\chi }\mathrm{ln}\mathit{\chi }〉-〈\mathit{\chi }〉\mathrm{ln}〈\mathit{\chi }〉,\end{array}$

where the brackets $〈\mathrm{\dots }〉\equiv \sum \left(\mathrm{\dots }\right){p}_{k}$ denote averages with respect to the distribution pk, i.e., $〈f\left(\mathit{\chi }\right)〉\equiv \sum f\left({\mathit{\chi }}_{k}\right){p}_{k}$. Notably, the functional given by Eq. (5) has been shown to exhibit positivity, concavity, and experimental stability, which are the three requirements in order to be characterized as entropic functional. Furthermore, note that the entropy S is a dynamic entropy depending on the sequential order of the events and not simply a statistical entropy (e.g., Shannon entropy) (see Varotsos et al.2005a). Upon considering time reversal $\stackrel{\mathrm{^}}{T}$, i.e., $\stackrel{\mathrm{^}}{T}{p}_{k}={p}_{N-k+\mathrm{1}}$, the value S changes to a value S:

$\begin{array}{ll}& {S}_{-}=\sum _{k=\mathrm{1}}^{N}{p}_{N-k+\mathrm{1}}\left(\frac{k}{N}\right)\mathrm{ln}\left(\frac{k}{N}\right)\\ \text{(6)}& & \phantom{\rule{1em}{0ex}}-\left(\sum _{k=\mathrm{1}}^{N}\frac{k}{N}{p}_{N-k+\mathrm{1}}\right)\mathrm{ln}\left[\sum _{l=\mathrm{1}}^{N}\frac{l}{N}{p}_{N-l+\mathrm{1}}\right].\end{array}$

The physical meaning of the entropy change $\mathrm{\Delta }S\equiv S-{S}_{-}$ in natural time under time reversal is discussed in and .

Using a moving window of length i (number of events) sliding through the time series of L consecutive events, the entropy in natural time is determined for each position $j=\mathrm{1},\mathrm{2},\mathrm{\dots },L-i$ of the sliding window. Thus, a time series of Si is obtained. By considering the standard deviation σSi) of the time series of $\mathrm{\Delta }{S}_{i}\equiv {S}_{i}-\left({S}_{-}{\right)}_{i}$, we define the complexity measure Λi:

$\begin{array}{}\text{(7)}& {\mathrm{\Lambda }}_{i}=\frac{\mathit{\sigma }\left(\mathrm{\Delta }{S}_{i}\right)}{\mathit{\sigma }\left(\mathrm{\Delta }{S}_{\mathrm{100}}\right)},\end{array}$

where the denominator has been selected to correspond to the standard deviation σS100) of the time series of ΔSi of i=100 events (e.g., Ramírez-Rojas et al.2018).

ΔS constitutes a key measure that may identify when the system approaches the critical point (dynamic phase transition). For example, ΔS has been applied for the identification of the time of an impending sudden cardiac death risk . Furthermore, it has been used for the study of the predictability of the Olami–Feder–Christensen (OFC) model for earthquakes , which is probably the most studied non-conservative self-organized criticality (SOC) model. The OFC model originated by a simplification of the Burridge and Knopoff spring–block model by mapping it into a non-conservative cellular automaton simulating the earthquake's behavior and introducing dissipation in the family of SOC systems. In particular, it was found that ΔS exhibits a clear minimum (or maximum if we define $\mathrm{\Delta }S\equiv {S}_{-}-S$ instead of $\mathrm{\Delta }S\equiv S-{S}_{-}$; e.g., see ) before a large avalanche in the OFC model, which corresponds to a large earthquake. For example, by analyzing the seismicity during the period 2012–2017 in natural time in the Chiapas region of Mexico where the M 8.2 earthquake occurred on 7 September 2017, we observed that the entropy change ΔS of seismicity under time reversal was minimized almost 3 months before and in particular on 14 June 2017.

## 3.2 Results from natural time analysis of seismicity

Interesting results have been recently obtained upon analyzing the Japan seismic catalog in natural time and computing the fluctuations of κ1. To compute the κ1 fluctuations, we use a sliding natural time window comprising the number i of EQs that would occur on average in a few months or so . We then calculate the average value μ(κ1) and the standard deviation σ(κ1) of the ensemble of κ1 obtained. The quantity

$\begin{array}{}\text{(8)}& {\mathit{\beta }}_{i}\equiv \mathit{\sigma }\left({\mathit{\kappa }}_{\mathrm{1}}\right)/\mathit{\mu }\left({\mathit{\kappa }}_{\mathrm{1}}\right)\end{array}$

is defined as the variability of κ1 . The time evolution of the β value can then be pursued by sliding the excerpt i through the EQ catalog, and the corresponding minimum value is labeled βmin. The following key results have been obtained.

found that the fluctuations β of κ1 of seismicity exhibited a clearly detectable minimum approximately at the time of the initiation of a pronounced SES activity recorded by around 2 months before the volcanic–seismic swarm activity in 2000 in the Izu Islands region, Japan.

analyzed the Japan seismic catalog in natural time from 1 January 1984 to 11 March 2011, and their results showed that the fluctuations β of κ1 of seismicity exhibited distinct minima βmin a few months before all the shallow earthquakes of magnitude 7.6 or larger that occurred during this 27-year period in the Japanese area ${N}_{\mathrm{25}}^{\mathrm{46}}{E}_{\mathrm{125}}^{\mathrm{148}}$. Among these minima, the minimum before the M 9 Tōhoku EQ observed at around 5 January 2011 was the deepest. Subsequently, found that the spatiotemporal variations of βmin enable the estimation of the epicentral area of the impending mainshock for all these EQs of magnitude 7.6 or larger.

In addition, focused on the minima βmin preceding all magnitude 8 (and 9) class EQs in the Japanese area from 1 January 1984 to 11 March 2011 and applied detrended fluctuation analysis (DFA) to the earthquake magnitude time series. DFA has been established as a standard method to investigate long-range correlations in non-stationary time series in diverse fields (e.g., see Peng et al.1993, 1994; Telesca and Lovallo2009) including the study of geomagnetic data associated with the M 9.0 Tōhoku EQ . The results of DFA are described in terms of the so-called DFA exponent , hereafter labeled α (α=0.5 means random, α greater than 0.5 means long-range correlations, and α less than 0.5 means anti-correlations). The following three main features have been identified : the minima βmin are observed during periods when long-range temporal correlations between EQ magnitudes have been developed since the corresponding DFA exponent is α>0.5. Before (bef) the minima βmin there exists a stage in which an evident anti-correlated behavior appears showing a minimum αmin,bef in the DFA exponent markedly smaller than 0.5. Finally, after (aft) the minima βmin, the long-range correlations break down to an almost random behavior (α=0.5), possibly turning to anti-correlation exhibiting a minimum αmin,aft with ${\mathit{\alpha }}_{\mathrm{min},\mathrm{aft}}<\mathrm{0.5}$. These three main features of temporal correlations between EQ magnitudes can be visualized in Fig. 3 for the case of the M 9 Tōhoku EQ. The deepest minimum βmin was observed around 5 January 2011 during a period in which long-range correlations between EQ magnitudes prevail with a DFA exponent α>0.5. Before this minimum an anti-correlated behavior was identified on 22 December 2010 with ${\mathit{\alpha }}_{\mathrm{min},\mathrm{bef}}=\mathrm{0.35}$. After βmin the long-range correlations break down on 13 January 2011 to an almost random behavior with α=0.5, and subsequently the behavior turned to anti-correlation on 23 January 2011 with ${\mathit{\alpha }}_{\mathrm{min},\mathrm{aft}}=\mathrm{0.42}$. As explained in the next section, the main features of the temporal correlations between EQ magnitudes obtained by DFA appeared simultaneously with distinct phases of crustal deformation identified by GPS measurements described in the previous section.

Recently we have shown that almost 3 months before the M 9 Tōhoku earthquake, i.e., on 22 December 2010, the following additional facts have been observed: first, the complexity measure Λi associated with the fluctuations of the entropy change of seismicity under time reversal exhibited an abrupt increase which conforms to the seminal work by Lifshitz and Slyozov and independently by 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\left(t-{t}_{\mathrm{0}}{\right)}^{\mathrm{1}/\mathrm{3}}$ . It was also found that the increase ΔΛi of Λi follows the latter form and that the prefactors A are proportional to the scale i, while the exponent (1∕3) is independent of i . Second, the Tsallis entropic index q (Tsallis1988) shows a simultaneous increase which interestingly exhibits the same exponent (1∕3) . Third, a minimum ΔSmin of the change ΔS of the entropy of seismicity in the entire Japanese region under time reversal was found by , who also demonstrated that the probability of obtaining such a minimum by chance is approximately 3 %, thus showing that it is statistically significant. In addition, the robustness of the appearance of this minimum on 22 December 2010 upon changing the EQ depth, the EQ magnitude threshold, and the size of the area investigated has been documented . Such a minimum is of precursory nature, signaling that a large EQ is impending according to the natural time analysis of the OFC model as mentioned in Sect. 3.1. Fourth, studying the fluctuations β of κ1 of seismicity in the entire Japanese region ${N}_{\mathrm{25}}^{\mathrm{46}}{E}_{\mathrm{125}}^{\mathrm{148}}$ versus the conventional time from 1 January 1984 until the Tōhoku EQ occurrence on 11 March 2011, we find a large fluctuation of β upon the occurrence of the M 7.8 earthquake near Chichijima on 22 December 2010. This finding has also been checked for several lengths from i=150 to 500 events, which also revealed the following : upon increasing i it is observed (see Figs. 2b and 4e of Sarlis et al.2013) that the increase Δβi of the βi fluctuation on 22 December 2010 becomes distinctly larger – obeying the interrelation $\mathrm{\Delta }{\mathit{\beta }}_{i}=\mathrm{0.5}\mathrm{ln}\left(i/\mathrm{114.3}\right)$ – which does not happen (see Fig. 4a–d of Sarlis et al.2013) for the increases in the β fluctuations upon the occurrences of all other shallow EQs in Japan of magnitude 7.6 or larger during the period from 1 January 1984 to the time of the M 9 Tōhoku EQ. Hence, the β fluctuation on 22 December 2010 accompanying the minimum ΔSmin is unique. Its presence is of paramount importance for the validity of the physical model that will be discussed in the next section.

4 Compatibility of the observed precursory phenomena with the PSPC physical model

Let us now discuss the multidisciplinary observations described in the previous two sections (and compiled in Fig. 3) that preceded the M 9 Tōhoku EQ. As we shall see all these observations are directly evidenced from the PSPC model, except probably of the anomalous changes in the level and temperature of confined groundwater.

A striking fact is that around 5 January 2011 the phenomenon of aligned orientations of the GPS azimuths occurred almost simultaneously with two other phenomena, i.e., the initiation of the anomalous Earth magnetic field variations and the minimum βmin of the fluctuations of the order parameter of seismicity. This is strikingly reminiscent of the mechanism of the emission of SES activity (stage B of PSPC model) in which, upon reaching σ=σcr, the electric dipoles exhibit cooperative orientation that also reflects alignment of the horizontal GPS azimuths (black in Fig. 3), as expected by , leading to the most intense crust uplift. Such an SES emission is directly evidenced by the observed anomalous variations of the Earth's magnetic field mainly in the vertical component (brown in Fig. 3). This emission is consistent with the observation of the deepest minimum βmin around 5 January 2011 of the fluctuations of the order parameter of seismicity (blue in Fig. 3) in view of the up-to-date experimental results mentioned above that an SES activity initiates almost simultaneously with both the observation of βmin and the establishment of long-range temporal correlation between earthquake magnitudes.

Another striking fact is that the aforementioned simultaneous appearance of the three phenomena around 5 January 2011 has been preceded by a stage of an evident anti-correlated behavior between earthquake magnitudes since it was found that ${\mathit{\alpha }}_{\mathrm{min},\mathrm{bef}}=\mathrm{0.35}$ upon the occurrence on 22 December 2010 of the M 7.8 EQ in southern Japan at 27.05 N, 143.94 E. On the same date, the horizontal GPS azimuths, which were initially random, started to become gradually oriented toward the southern direction probably due to an excess stress disturbance. This may also originate the simultaneous appearance of various phenomena, including the large abrupt increase in the order parameter fluctuations along with an abrupt increase in the complexity measure Λi of the change in the entropy of seismicity under time reversal (recall that this change is then minimized); an increase in the Tsallis entropic index q (Tsallis1988); and start of groundwater anomalous changes, i.e., groundwater level drop, temperature decrease, and increase in radon concentration. This corresponds to stage A of the PSPC model, according to which an excess stress disturbance starts gradually increasing until reaching σcr.

After the occurrence of βmin at around 5 January 2011, the intense crust uplift was gradually mitigated and the orientations of GPS azimuths returned to random around 13 January 2011, thus agreeing with the DFA exponent α=0.5 . The behavior turned to anti-correlation around 23 January 2011 with DFA exponent ${\mathit{\alpha }}_{\mathrm{min},\mathrm{aft}}=\mathrm{0.42}$, and a shift of earthquake-related stress disturbance was observed , where westward movements replaced the southward ones; i.e., the orientations of the residual displacements were realigned along the western direction and the crust depressed. After this change on 23 January 2011 the stress disturbance gradually approached the threshold of the fault rupture, and the orientations of the residual displacements became random again , in agreement with the DFA exponent of the earthquake magnitude time series being close to 0.5 until around 10 February 2011 (see Fig. 5 of Varotsos et al.2014), which indicates random behavior. This fact that the Tōhoku EQ occurred after the emergence of an almost random behavior did not come as a surprise since it is strikingly reminiscent of similar findings in other complex time series as follows: in the case of electrocardiograms, for example, the long-range temporal correlations that characterize the healthy heart rate variability break down for individuals at high risk of sudden cardiac death (SCD), and this is often accompanied by the emergence of uncorrelated randomness (SCD could be viewed as a critical phenomenon; e.g., see ).

5 Summary and conclusions

Several phenomena of multidisciplinary nature preceded the M 9 Tōhoku EQ that occurred on 11 March 2011. Leaving aside the details, these phenomena were mainly accumulated around two dates, i.e., 22 December 2010 and 5 January 2011, which concur with the two stages A and B of the PSPC physical model, respectively. These phenomena include the following.

• a.

Around the date 22 December 2010:

• 1.

The entropy change of seismicity under time reversal is minimized along with increased fluctuations (since Λi increases).

• 2.

There is an increase in the fluctuations of the order parameter of seismicity.

• 3.

The DFA exponent decreased to the value ${\mathit{\alpha }}_{\mathrm{min},\mathrm{bef}}=\mathrm{0.35}$, which is the lowest observed during the period 1984–2011 of our study, pointing to an evident anti-correlated behavior in the earthquake magnitude time series.

• 4.

The horizontal GPS azimuths started to become gradually oriented toward the southern direction (while they had random orientations during the preceding period 12–22 December 2010).

• 5.

Anomalous changes in the groundwater started (level drop, temperature decrease, and probably increase in radon concentration).

• 6.

Increase in the Tsallis entropic index q.

• b.

Around the date 5 January 2011:

• 1.

Unprecedented minimum βmin of the fluctuations of the order parameter of seismicity.

• 2.

Anomalous magnetic field variations started (which, according to Maxwell equations, should be accompanied by a strong SES activity).

• 3.

Full alignment of the orientations of the GPS azimuths southwards accompanied by the most intense crust uplift.

• 4.

Long-range temporal correlations in the earthquake magnitude time series.

All the above phenomena were observed to begin and end well before the M 9 Tōhoku EQ occurrence as schematically shown in Fig. 1a in accordance with the PSPC model (Fig. 2).

Data availability
Data availability.

Author contributions
Author contributions.

PAV, NVS, and ESS designed the research, performed the research, analyzed the results, and reviewed the manuscript.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Review statement
Review statement.

This paper was edited by Georgios Balasis and reviewed by two anonymous referees.

References

Bernardi, A., Fraser-Smith, A. C., McGill, P. R., and Villard, O. G.: ULF magnetic field measurements near the epicenter of the Ms 7.1 Loma Prieta earthquake, Phys. Earth Planet. In., 68, 45–63, https://doi.org/10.1016/0031-9201(91)90006-4, 1991. a

Burridge, R. and Knopoff, L.: Model and theoretical seismicity, B. Seismol. Soc. Am., 57, 341–371, 1967. a

Carlson, J. M., Langer, J. S., and Shaw, B. E.: Dynamics of earthquake faults, Rev. Mod. Phys., 66, 657–670, https://doi.org/10.1103/RevModPhys.66.657, 1994. a, b

Chen, C.-H., Yeh, T.-K., Liu, J.-Y., Wang, C.-H., Wen, S., Yen, H.-Y., and Chang, S.-H.: Surface Deformation and Seismic Rebound: Implications and Applications, Surv. Geophys., 32, 291–313, https://doi.org/10.1016/j.jseaes.2013.11.009, 2011. a, b

Chen, C.-H., Wen, S., Liu, J.-Y., Hattori, K., Han, P., Hobara, Y., Wang, C.-H., Yeh, T.-K., and Yen, H.-Y.: Surface displacements in Japan before the 11 March 2011 M 9.0 Tohoku-Oki earthquake, J. Asian Earth Sci., 80, 165–171, https://doi.org/10.1016/j.jseaes.2013.11.009, 2014. a, b, c, d, e, f, g

Fan, Y.-Y., Du, X.-B., Zlotnicki, J., Tan, D.-C., An, Z.-H., Chen, J.-Y., Zheng, G.-L., Liu, J., and Xie, T.: The Electromagnetic Phenomena Before the Ms8.0 Wenchuan Earthquake, Chinese J. Geophys., 53, 997–1010, https://doi.org/10.1002/cjg2.1570, 2010. a

Fraser-Smith, A. C., Bernardi, A., McGill, P. R., Ladd, M. E., Helliwell, R. A., and Villard, O. G.: Low-frequency magnetic-field measurements near the epicenter of the Ms-7.1 Loma-Prieta earthquake, Geophys. Res. Lett., 17, 1465, https://doi.org/10.1029/GL017i009p01465, 1990. a

Goldberger, A. L., Amaral, L. A. N., Hausdorff, J. M., Ivanov, P. C., Peng, C.-K., and Stanley, H. E.: Fractal dynamics in physiology: Alterations with disease and aging, P. Natl. Acad. Sci. USA, 99, 2466–2472, https://doi.org/10.1073/pnas.012579499, 2002. a

Han, P., Hattori, K., Xu, G., Ashida, R., Chen, C.-H., Febriani, F., and 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., 114, 321–326, https://doi.org/10.1016/j.jseaes.2015.02.022, 2015. a, b

Han, P., Hattori, K., Huang, Q., Hirooka, S., and 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., 129, 13–21, https://doi.org/10.1016/j.jseaes.2016.07.011, 2016. a, b

Holliday, J. R., Rundle, J. B., Turcotte, D. L., Klein, W., Tiampo, K. F., and Donnellan, A.: Space-Time Clustering and Correlations of Major Earthquakes, Phys. Rev. Lett., 97, 238501, https://doi.org/10.1103/PhysRevLett.97.238501, 2006. a

Huang, Q.: Retrospective investigation of geophysical data possibly associated with the Ms8.0 Wenchuan earthquake in Sichuan, China, J. Asian Earth Sci., 41, 421–427, https://doi.org/10.1016/j.jseaes.2010.05.014, 2011. a

Kanamori, H.: Quantification of Earthquakes, Nature, 271, 411–414, https://doi.org/10.1038/271411a0, 1978. a

Lifshitz, I. and Slyozov, V.: The kinetics of precipitation from supersaturated solid solutions, J. Phys. Chem. Solids, 19, 35–50, https://doi.org/10.1016/0022-3697(61)90054-3, 1961. a

Luginbuhl, M., Rundle, J. B., Hawkins, A., and Turcotte, D. L.: Nowcasting Earthquakes: A Comparison of Induced Earthquakes in Oklahoma and at the Geysers, California, Pure Appl. Geophys., 175, 49–65, https://doi.org/10.1007/s00024-017-1678-8, 2018a. a

Luginbuhl, M., Rundle, J. B., and Turcotte, D. L.: Natural Time and Nowcasting Earthquakes: Are Large Global Earthquakes Temporally Clustered?, Pure Appl. Geophys., 175, 661–670, https://doi.org/10.1007/s00024-018-1778-0, 2018b. a

Molchanov, O. A. and Hayakawa, M.: Generation of ULF electromagnetic emissions by microfracturing, Geophys. Res. Lett., 22, 3091–3094, https://doi.org/10.1029/95GL00781, 1995. a

Olami, Z., Feder, H. J. S., and Christensen, K.: Self-organized criticality in a continuous, nonconservative cellular automaton modeling earthquakes, Phys. Rev. Lett., 68, 1244–1247, https://doi.org/10.1103/physrevlett.68.1244, 1992. a

Orihara, Y., Kamogawa, M., and Nagao, T.: Preseismic Changes of the Level and Temperature of Confined Groundwater related to the 2011 Tohoku Earthquake, Sci. Rep., 4, 6907, https://doi.org/10.1038/srep06907, 2014. a, b, c, d

Peng, C.-K., Buldyrev, S. V., Goldberger, A. L., Havlin, S., Simons, M., and Stanley, H. E.: Finite-size effects on long-range correlations: Implications for analyzing DNA sequences, Phys. Rev. E, 47, 3730–3733, https://doi.org/10.1103/physreve.47.3730, 1993. a, b, c

Peng, C.-K., Buldyrev, S. V., Havlin, S., Simons, M., Stanley, H. E., and Goldberger, A. L.: Mosaic organization of DNA nucleotides, Phys. Rev. E, 49, 1685–1689, https://doi.org/10.1103/physreve.49.1685, 1994. a, b, c

Ramírez-Rojas, A., Telesca, L., and Angulo-Brown, F.: Entropy of geoelectrical time series in the natural time domain, Nat. Hazards Earth Syst. Sci., 11, 219–225, https://doi.org/10.5194/nhess-11-219-2011, 2011. a

Ramírez-Rojas, A., Flores-Márquez, E. L., Sarlis, N. V., and 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, 20, 477, https://doi.org/10.3390/e20060477, 2018. a, b

Ramos, O., Altshuler, E., and Måløy, K. J.: Quasiperiodic Events in an Earthquake Model, Phys. Rev. Lett., 96, 098501, https://doi.org/10.1103/physrevlett.96.098501, 2006. a

Rikitake, T.: Current research in earthquake prediction, D. Reidler Publishing Company Dordrecht Holland, 1981. a, b

Rong, Y., Wang, Q., Ding, X., and Huang, Q.: Non-uniform scaling behavior in Ultra-Low-Frequency (ULF) geomagnetic signals possibly associated with the 2011 M9.0 Tohoku earthquake, Chinese J. Geophys., 55, 3709, https://doi.org/10.6038/j.issn.0001-5733.2012.11.018, 2012. a

Rundle, J. B., Turcotte, D. L., Donnellan, A., Grant Ludwig, L., Luginbuhl, M., and Gong, G.: Nowcasting earthquakes, Earth Space Sci., 3, 480–486, https://doi.org/10.1002/2016EA000185, 2016. a

Rundle, J. B., Luginbuhl, M., Giguere, A., and Turcotte, D. L.: Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities, Pure Appl. Geophys., 175, 647–660, https://doi.org/10.1007/s00024-017-1720-x, 2018. a

Sarlis, N. and Varotsos, P.: Magnetic field near the outcrop of an almost horizontal conductive sheet, J. Geodyn., 33, 463–476, https://doi.org/10.1016/S0264-3707(02)00008-X, 2002. a

Sarlis, N., Skordas, E., and Varotsos, P.: The change of the entropy in natural time under time-reversal in the Olami-Feder-Christensen earthquake model, Tectonophysics, 513, 49–53, https://doi.org/10.1016/j.tecto.2011.09.025, 2011. a, b

Sarlis, N. V., Skordas, E. S., Varotsos, P. A., Nagao, T., Kamogawa, M., Tanaka, H., and Uyeda, S.: Minimum of the order parameter fluctuations of seismicity before major earthquakes in Japan, P. Natl. Acad. Sci. USA, 110, 13734–13738, https://doi.org/10.1073/pnas.1312740110, 2013. a, b, c, d, e

Sarlis, N. V., Skordas, E. S., Varotsos, P. A., Nagao, T., Kamogawa, M., and Uyeda, S.: Spatiotemporal variations of seismicity before major earthquakes in the Japanese area and their relation with the epicentral locations, P. Natl. Acad. Sci. USA, 112, 986–989, https://doi.org/10.1073/pnas.1422893112, 2015. a, b, c

Sarlis, N. V., Skordas, E. S., and Varotsos, P. A.: A remarkable change of the entropy of seismicity in natural time under time reversal before the super-giant M9 Tohoku earthquake on 11 March 2011, Europhys. Lett., 124, 29001, https://doi.org/10.1209/0295-5075/124/29001, 2018a. a, b

Sarlis, N. V., Skordas, E. S., Varotsos, P. A., Ramírez-Rojas, A., and Flores-Márquez, E. L.: Natural time analysis: On the deadly Mexico M8.2 earthquake on 7 September 2017, Physica A, 506, 625–634, https://doi.org/10.1016/j.physa.2018.04.098, 2018b. a, b

Telesca, L. and 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., 36, L01308, https://doi.org/10.1029/2008GL036247, 2009. a

Tsallis, C.: Possible generalization of Boltzmann-Gibbs statistics, J. Stat. Phys., 52, 479–487, https://doi.org/10.1007/BF01016429, 1988. a, b

Tsunomori, F. and Tanaka, H.: Anomalous change of groundwater radon concentration monitored at Nakaizu well in 2011, Radiat. Meas., 60, 35–41, https://doi.org/10.1016/j.radmeas.2013.11.006, 2014. a

Uyeda, S., Nagao, T., Orihara, Y., Yamaguchi, T., and Takahashi, I.: Geoelectric potential changes: Possible precursors to earthquakes in Japan, P. Natl. Acad. Sci. USA, 97, 4561–4566, https://doi.org/10.1073/pnas.97.9.4561, 2000. a

Uyeda, S., Hayakawa, M., Nagao, T., Molchanov, O., Hattori, K., Orihara, Y., Gotoh, K., Akinaga, Y., and Tanaka, H.: Electric and magnetic phenomena observed before the volcano-seismic activity in 2000 in the Izu Island Region, Japan, P. Natl. Acad. Sci. USA, 99, 7352–7355, https://doi.org/10.1073/pnas.072208499, 2002. a, b

Uyeda, S., Kamogawa, M., and Tanaka, H.: Analysis of electrical activity and seismicity in the natural time domain for the volcanic-seismic swarm activity in 2000 in the Izu Island region, Japan, J. Geophys. Res., 114, B02310, https://doi.org/10.1029/2007JB005332, 2009a. a, b

Uyeda, S., Nagao, T., and Kamogawa, M.: Short-term earthquake prediction: Current status of seismo-electromagnetics, Tectonophysics, 470, 205–213, https://doi.org/10.1016/j.tecto.2008.07.019, 2009b. a

Varotsos, P.: The Physics of Seismic Electric Signals, TERRAPUB, Tokyo, 2005. a

Varotsos, P.: Point defect parameters in β-PbF2 revisited, Solid State Ionics, 179, 438–441, https://doi.org/10.1016/j.ssi.2008.02.055, 2008. a

Varotsos, P. and Alexopoulos, K.: Physical Properties of the variations of the electric field of the earth preceding earthquakes, I, Tectonophysics, 110, 73–98, https://doi.org/10.1016/0040-1951(84)90059-3, 1984a. a, b

Varotsos, P. and Alexopoulos, K.: Physical Properties of the variations of the electric field of the earth preceding earthquakes, II, Tectonophysics, 110, 99–125, https://doi.org/10.1016/0040-1951(84)90060-X, 1984b. a, b

Varotsos, P. and Alexopoulos, K.: Thermodynamics of Point Defects and their Relation with Bulk Properties, North Holland, Amsterdam, 1986. a, b

Varotsos, P. and Lazaridou, M.: Latest aspects of earthquake prediction in Greece based on Seismic Electric Signals, Tectonophysics, 188, 321–347, https://doi.org/10.1016/0040-1951(91)90462-2, 1991. a, b

Varotsos, P., Alexopoulos, K., Nomicos, K., and Lazaridou, M.: Earthquake prediction and electric signals, Nature, 322, 120, https://doi.org/10.1038/322120a0, 1986. a

Varotsos, P., Alexopoulos, K., and Lazaridou, M.: Latest aspects of earthquake prediction in Greece based on Seismic Electric Signals,II, Tectonophysics, 224, 1–37, https://doi.org/10.1016/0040-1951(93)90055-O, 1993. a, b

Varotsos, P., Sarlis, N., and Skordas, E.: Scale-specific order parameter fluctuations of seismicity in natural time before mainshocks, Europhys. Lett., 96, 59002, https://doi.org/10.1209/0295-5075/96/59002, 2011. a, b

Varotsos, P. A., Sarlis, N. V., and Skordas, E. S.: Attempt to distinguish electric signals of a dichotomous nature, Phys. Rev. E, 68, 031106, https://doi.org/10.1103/PhysRevE.68.031106, 2003a. a

Varotsos, P. A., Sarlis, N. V., and Skordas, E. S.: Electric Fields that “arrive” before the time derivative of the magnetic field prior to major earthquakes, Phys. Rev. Lett., 91, 148501, https://doi.org/10.1103/PhysRevLett.91.148501, 2003b. a

Varotsos, P. A., Sarlis, N. V., Skordas, E. S., and Lazaridou, M. S.: Entropy in Natural Time Domain, Phys. Rev. E, 70, 011106, https://doi.org/10.1103/physreve.70.011106, 2004. a, b, c

Varotsos, P. A., Sarlis, N. V., Skordas, E. S., and Lazaridou, M. S.: Natural entropy fluctuations discriminate similar-looking electric signals emitted from systems of different dynamics, Phys. Rev. E, 71, 011110, https://doi.org/10.1103/physreve.71.011110, 2005a. a, b

Varotsos, P. A., Sarlis, N. V., Tanaka, H. K., and Skordas, E. S.: Some properties of the entropy in the natural time, Phys. Rev. E, 71, 032102, https://doi.org/10.1103/physreve.71.032102, 2005b. a

Varotsos, P. A., Sarlis, N. V., Tanaka, H. K., and Skordas, E. S.: Similarity of fluctuations in correlated systems: The case of seismicity, Phys. Rev. E, 72, 041103, https://doi.org/10.1103/physreve.72.041103, 2005c. a

Varotsos, P. A., Sarlis, N. V., Skordas, E. S., and Lazaridou, M. S.: Identifying sudden cardiac death risk and specifying its occurrence time by analyzing electrocardiograms in natural time, Appl. Phys. Lett., 91, 064106, https://doi.org/10.1063/1.2768928, 2007. a, b, c

Varotsos, P. A., Sarlis, N. V., and 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, https://doi.org/10.1007/978-3-642-16449-1, 2011. a, b, c, d, e, f, g, h, i, j

Varotsos, P. A., Sarlis, N. V., Skordas, E. S., and Lazaridou, M. S.: Seismic Electric Signals: An additional fact showing their physical interconnection with seismicity, Tectonophysics, 589, 116–125, https://doi.org/10.1016/j.tecto.2012.12.020, 2013.  a

Varotsos, P. A., Sarlis, N. V., and Skordas, E. S.: Study of the temporal correlations in the magnitude time series before major earthquakes in Japan, J. Geophys. Res.-Space, 119, 9192–9206, https://doi.org/10.1002/2014JA020580, 2014. a, b, c, d, e, f

Varotsos, P. A., Sarlis, N. V., and 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, 20, 757, https://doi.org/10.3390/e20100757, 2018. a, b, c

Varotsos, P. A., Sarlis, N. V., and Skordas, E. S.: Natural time analysis: Important changes of the order parameter of seismicity preceding the 2011 M9 Tohoku earthquake in Japan, Europhys. Lett., 125, 69001, https://doi.org/10.1209/0295-5075/125/69001, 2019. a, b

Wagner, C.: Theorie der Alterung von Niederschlägen durch Umlösen (Ostwald-Reifung), Zeitschrift für Elektrochemie, Berich. Bunsen. Gesell., 65, 581–591, 1961. a

Xu, G., Han, P., Huang, Q., Hattori, K., Febriani, F., and 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., 77, 59–65, https://doi.org/10.1016/j.jseaes.2013.08.011, 2013. a, b