## Identifying the occurrence time of an impending mainshock: A very recent case.

P. A. Varotsos, 1, 2, ∗ N. V. Sarlis, 1, 2 E. S. Skordas, 1, 2 and M. S. Lazaridou-Varotsos 2

1 Section of Solid State Physics, Physics Department, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos 157 84, Athens, Greece 2 Solid Earth Physics Institute, Physics Department, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos 157 84, Athens, Greece

The procedure by means of which the occurrence time of an impending mainshock can be identified by analyzing in natural time the seismicity in the candidate area subsequent to the recording of a precursory Seismic Electric Signals (SES) activity is reviewed. Here, we report the application of this procedure to an Mw5.4 mainshock that occurred in Greece on 17 November 2014 and was strongly felt in Athens. This mainshock (which is pretty rare since it is the strongest in that area for more than half a century) was preceded by an SES activity recorded on 27 July 2014 and the results of the natural time analysis reveal that the system approached the critical point (mainshock occurrence) early in the morning on 15 November 2014. SES activities that have been recently recorded are also presented. Furthermore, in a Note we discuss the case of the Mw5.3 earthquake that was also strongly felt in Athens on 19 July 2019 (Parnitha fault).

PACS numbers: 05.40.-a, 05.45.Tp, 91.30.Dk, 89.75.-k

## I. INTRODUCTION

Earthquakes (EQs) in general exhibit complex correlations in time, space, and magnitude M which have been investigated by several authors [1-20]. The earthquake scaling laws [21](Turcotte 1997) indicate the existence of phenomena closely associated with the proximity of the system to a critical point, e.g., see Holliday et al. [4]. Here, we take this view that mainshocks are (non-equilibrium) critical phenomena. Major EQs are preceded by transient changes of the electric field of the Earth termed Seismic Electric Signals (SES) [22, 23](Varotsos and Alexopoulos 1984a,b). A series of such signals recorded within a short time are called SES activities [24-27], the average lead time of which is of the order of a few months [28]. It has been suggested that SES are emitted when the stress in the focal area of the impending mainshock reaches a critical value [22, 23, 29, 30]. This suggestion is strengthened by the finding [31] that the fluctuations of the order parameter of seismicity defined in the frame of natural time analysis (see the next section) minimize upon the initiation of an SES activity exhibiting long range temporal correlations[32]. Such minima of the fluctuations of the order parameter of seismicity have been identified before all major (M ≥ 7.6) EQs in Japan[33, 34]. The identification of the occurrence time of an impending mainshock within a short time window is a challenge. This becomes possible when employing a procedure that combines SES data and natural time analysis of the seismicity [35-40]. In short, the initiation of the SES activity marks the time when the system enters the critical stage and then the

∗ Electronic address: pvaro@otenet.gr

natural time analysis of the subsequent seismicity in the candidate area (which is determined on the basis of SES data, e.g., see Varotsos [26]) identifies when the system approaches the critical point, i.e., the mainshock occurrence, e.g., see Fig.1 of Huang et al. [41]. It is one of the aims of this paper to report a characteristic application of this procedure, which refers to an SES activity that was followed by a pronounced Mw=5.4 mainshock in Greece on 17 November 2014, which is pretty rare as explained later.

## II. SUMMARY OF THE PROCEDURE TO IDENTIFY THE OCCURRENCE TIME OF AN IMPENDING MAINSHOCK

Let us first summarize the natural time analysis[36] in the case of seismicity: In a time series comprising N EQs, the natural time χ k = k/N serves as an index for the occurrence of the k-th EQ. The combination of this index with the energy Q k released during the k-th EQ of magnitude M , i.e., k the pair ( χ , Q k k ), is studied in natural time analysis. Alternatively, one studies the pair ( χ , p k k ), where p k = Q / k ∑ N n =1 Q n stands for the normalized energy released during the k-th EQ. It has been found that the variance of χ weighted for p k , designated by κ 1 , which is given by[28, 35-37, 42, 43]

$$\stackrel { \infty } { \text{ES} } \\ \underset { \substack { \leftarrow \\ \text{the} } } { \leftarrow } \rightarrow \\ \text{the} \rightarrow \rightarrow \left ( \chi \right ) - \left < \chi \right > ^ { 2 } = \sum _ { k = 1 } ^ { N } p _ { k } ( \chi _ { k } ) ^ { 2 } - \left ( \sum _ { k = 1 } ^ { N } p _ { k } \chi _ { k } \right ) ^ { 2 }. \quad ( 1 )$$

plays a prominent role in natural time analysis. In particular, κ 1 may serve as an order parameter for seismicity[38] and it has been empirically observed[28, 35-40, 44] that κ 1 of the seismicity in the candidate area

above a magnitude threshold M thres subsequent to an SES activity becomes equal to 0.070 when approaching the critical point ( mainshock occurrence ). Note that Q k , and hence p k , for earthquakes is estimated through the usual relation[45]: Q k ∝ 10 1 5M . k .

Upon the recording of an SES activity, one can estimate an area A within which the impending mainshock is expected to occur. The magnitude M of the expected EQ is estimated through the relation log 10 ( ∆ V L ) ≈ 0 3 . M + const , e.g., see [24], where for a given measuring dipole . of length L and a given seismic area the SES amplitude ∆ V/L is found from the anomalous variation ∆ V of the potential difference between the corresponding two electrodes. When area A reaches criticality, one expects in general that all its subareas have also reached criticality simultaneously. At that time, therefore, the evolution of seismicity in each of its subareas is expected to result in κ 1 values close to 0.070. Assuming equi-partition of probability among the subareas[39], the distribution Prob( κ 1 ) of the κ 1 values of all subareas should be peaked at around 0.070 exhibiting also magnitude threshold invariance. This usually occurs a few days to around one week before the mainshock, thus it enables the prediction of the occurrence time of major EQs with time window of the order of a week or less.

## III. APPLICATION TO A PRONOUNCED SEISMIC ACTIVITY IN GREECE

The SES activity shown in Fig.1(a) was recorded on 27 July 2014 at Keratea (KER) geoelectrical station, the location of which is depicted with the red bullet in Fig.1(b). On the basis of the selectivity map of this station (i.e, the map showing all seismic areas in the past that gave rise to SES recorded at this station, e.g., see Varotsos and Lazaridou [24] and the ratio of the SES components the candidate area was determined [46]. This is depicted here by the rectangle in Fig.1(b) as was designated in the uppermost right part of Fig.2 of the paper uploaded by Sarlis et al. [46] on 7 August 2014.

Wenow proceed to the natural time analysis of the seismicity subsequent to the aforementioned SES activity at KERwithin the candidate area N(37.7-39.0)E(22.6-24.2). The EQ catalogue of the Institute of Geodynamics of the National Observatory of Athens available on 2 February 2015 at http://www.gein.noa.gr/services/current\_ catalogue.php was used, e.g. see [47, 48]. Figure 2(a) depicts Prob( κ 1 ) versus κ 1 of seismicity for M thres =2.8 (the data used are compiled in Table 1 of [49]) for the period after 27 October 2014, i.e., almost three weeks before the mainshock occurrence on 17 November 2014. During this period six smaller EQs occurred and we observe that Prob( κ 1 ) maximizes at κ 1 =0.070 upon the occurrence of the last EQ, i.e., the ML=2.8 EQ at 01:01 UTC on 15 November 2014. It is remarkable that the same behavior is observed in Fig.2(b) where in the computation of the κ 1 values we discarded from the seismicity of the

FIG. 1: (color online) . (a) The SES activity of dichotomous nature recorded at the Keratea (KER) geoelectrical station of the SES telemetric network. (b)The predicted epicental area designated by the rectangle on a map in which the location of the KER station (red bullet) is shown along with that of other geolectrical stations Lamia (LAM), Loutraki (LOU) and Pirgos (PIR) (black bullets). The epicenters of the strongest EQs in Greece (Mw ≥ 6.5) during the last decade are also shown with stars. The central station of the SES telemetric network is located at Athens (ATH, black square).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-c577b1e4-8307-4d5f-83c3-eff738600d42.jpg)

candidate area N(37.7-39.0)E(22.6-24.2) the EQs that occurred within the subarea N(37.7-38.3)E(22.6-23.3). This is consistent with the fact that the latter subarea constitutes the preliminary selectivity map of the LOU station, see Fig.1(b), which however did not show any SES activity simultaneously with the one initiated on 27 July 2014 at KER station (alternatively, the area resulting from the subtraction of the above two areas could have been announced as a candidate area for the impending mainshock). To assure that this behavior exhibits also magnitude threshold invariance, we repeated the calculation that resulted in Fig.2(b), but for low magnitude thresholds (so that to have a large number of EQs). In particular, Figs.2(c), (d), (e), (f) depict the corresponding results for M thres =1.8, 1.9, 2.0, and 2.1, respectively, which do show that Prob( κ 1 ) versus κ 1 exhibit local maximum at κ 1 =0.070 upon the occurrence of the aforementioned EQ on 15 November 2014 (the seismic data used in order to obtain Figs.2(b) to 2(f) are given in Table 2 of [49]). Actually, almost three days later, i.e, at 23:05 UTC on 17 November 2014, the Mw(USGS)=5.4 EQ occurred with an epicenter at 38.67oN,23.39oE (followed by a smaller Mw(USGS)=5.1 EQ at 23:09 UTC with epicenter at 38.68oN,23.24oE). It should be mentioned that EQs of such magnitude occur there very rarely. In particular, no EQ with Mw(USGS) ≥ 5.4 took place within the coordinates N(38.3-39.0)E(23.0-23.8) since 1965. In view of this very rare occurrence, it is interesting to study this case in the future by employing an approach[50] which uses SES and a neural network (trained by relevant data of earlier cases) to predict the magnitude and the occurrence time of the forthcoming EQ.

FIG. 2: How the histograms of Prob( κ 1 ) versus κ 1 evolve event by event in the natural time analysis of the seismicity subsequent to the initiation of the SES activity depicted in Fig.1(a). In each panel, the magnitude threshold (M thres ) used in the calculation is also depicted. For details on the exact (sub)areas within the rectangle of Fig.1(b) considered in each panel see Section 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-76b44492-89a2-4bd2-97a5-7987a8678965.jpg)

## IV. CONCLUSIONS

A pronounced Mw(USGS)=5.4 EQ was strongly felt at Athens, Greece, on 17 November 2014. This is pretty rare since it is the strongest EQ that occurred in that area since 1965. The procedure based on natural time analysis of the seismicity subsequent to an SES activity recorded on 27 July 2014 at the KER station close to Athens revealed that the system approached the critical point (mainshock occurrence) just a few days before, i.e., on 15 November 2014.

More recent SES activities: Despite severe experimental difficulties during the current period, it seems that an SES activity of more or less similar polarity has recently been recorded at KER on 15 March 2020 (Fig.4). Natural time analysis of the subsequent seismicity in the area designated by the rectangle in Fig.1(b) was carried out. (Remarkably, such an analysis has just been reported ( www. nature.com/articles/s41598-020-59333-4 ) as being a powerful tool to detect the onset of acceleration as an early warning of an impending failure.) The results of this analysis have been described in the two

previous versions of this preprint and followed by the EQs that occurred on 3 September, 11 September, and 2 December 2020 with epicenters at 38.17 N23.99 E, o o 38.12 N23.18 E, o o and 38.33 N23.46 E o o of magnitude Ms(ATH)=4.8, Ms(ATH)=4.7, and Ms(ATH)=4.9, respectively.

Around 18:00 UTC on 8 December 2020 an SES activity appeared on KER station (Fig. 5) obeying the properties described in the main text and the Section 7.2 of Ref.[28]. Beyond the results obtained before the ML(ATH)=6.0 earthquake on 3 March 2021 we report the following: By applying the same analysis in natural time of the subsequent seismicity in the area described above we find that upon the occurrence of the ML(ATH)=2.4 earthquake at 04:15 UTC on 18 March 2021 with an epicenter at 38.30 N23.69 E we find that o o Prob( κ 1 ) versus κ 1 exhibit maximum at κ 1 =0.070 upon considering M thres =2.0 and 2.2 as shown in Figs. 6A and B, respectively. The same holds upon the occurrence of the ML(ATH)=2.1 earthquake at 01:24 UTC on 18 March 2021 with an epicenter at 38.16 N22.91 E as o o shown in Figs 7A and B upon considering M thres =2.1 and 2.2.

In continuation of this analysis in the area described above -see Fig. 1(b)we find that Prob( κ 1 ) versus κ 1 exhibit maximum at κ 1 =0.070, for M thres =2.3 and 2.4, upon the occurrence at 08:02 UTC on 24 March 2021 of the ML(ATH)=2.8 EQ with an epicenter at 38.30 N23.68 E (Fig. 8). In addition, a principal peak at o o κ 1 =0.070 is observed upon the occurrence at 06:15 UTC on 25 March 2021 (see Fig. 9 for M thres =3.0) of the ML(ATH)=3.5 EQ with an epicenter at 38.76 N23.41 E o o as well as a maximum at κ 1 =0.070, for M thres =2.1 and 2.2, upon the occurrence at 04:15 UTC on 28 March 2021 of a ML(ATH)=2.1 EQ with an epicenter at 38.76 o N23.38 E. o

On 2 June 2021 an SES activity was recorded at KER geoelectrical station (Fig. 10) as in Fig. 3 . In the natural time analysis of the subsequent seismicity in the area described above -see Fig. 1(b)the critical condition κ 1 = 0 070 was found to exhibit magnitude threshold . invariance in the range M thres =2.9 to M thres =3.5, upon the occurrence at 10:08 UTC on 12 September 2021 of the ML(ATH)=3.6 EQ with an epicenter at 38.07 N23.76 E. o o

On 15 November 2021 an additional SES activity was recorded at KER station (Fig. 11) pointing to the continuation of the process of approaching criticality in the aforementioned region under investigation.

Actually, a ML(ATH)=4.0, i.e., Ms(ATH)=4.5, EQ occurred at 16:13 UTC on 26 December 2021 with an epicenter at 38.10 N23.14 E o o lying inside the expected rectangular area of Fig. 1(b) felt also in Athens. This EQ was preceded at 17:13 UTC on 25 December 2021 by an ML(ATH)=2.4 event with an epicenter at 38.31 N23.41 E upon the occurrence of which the condio o tion κ 1 = 0 070 was fulfilled for . M thres =1.7 and 1.8 (Fig. 12). Strikingly, the Ms(ATH)=4.5 EQ was also followed by an ML(ATH)=2.1 event at 20:39 UTC on 26 Decem-

FIG. 3: (color online) The SES activity recorded at KER geolectrical station on 19 June 2019.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-fc76a29f-6458-44fe-9b5b-41ad36a00883.jpg)

ber 2021 with an epicenter at 38.52 N23.95 E upon the o o occurrence of which the condition κ 1 = 0 070 was also . fulfilled for M thres =1.7 and 1.8 (Fig. 13) pointing to the conclusion that the approach to the critical point is still in progress.

On 25 February 2022, upon the occurrence of the ML(ATH)=3.3 EQ at 17:19 UTC with an epicenter at 38.49 N23.61 E, we found that the criticality condition o o κ 1 = 0 070 has been obeyed exhibiting magnitude thresh-. old invariance in the range M thres =2.0 to M thres =2.5 (Fig.14), which signals the approach of the system to the critical point.

From 5 April 2022 until 8 April 2022, the most wide magnitude threshold invariance of the criticality condition κ 1 = 0 070 (from . M thres =2.3 to M thres =3.2, see Fig. 15) concerning the region under investigation has been observed.

KERATEA STATION 15/03/2020

15

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-6f9e845b-ceec-4b8b-8d4c-c7ece0d1f5e4.jpg)

✁

Mar

✁

2020

FIG. 4: The SES activity recorded at KER geolectrical station on 15 March 2020.

FIG. 5: The SES activity at KER station (raw data) on 8 December 2020

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-96ba087c-db44-4a8f-b279-6e0b9cd4da34.jpg)

)

κ

1

Prob(

A

0.5

0.4

0.3

0.2

0.1

0.00

0.02

0.04

0.06

κ

1

0.08

0.10

0.12

120

130

No of EQs after SES

0.14

110

140

150

FIG. 6: The same as Figure 2 but after the SES activity at KER station on 8 December 2020, upon the occurrence of the ML(ATH)=2.4 EQ on 18 March 2021 with an epicenter at 38.30 N23.69 E by considering M o o thres =2.0 (A) and 2.2 (B), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-74b30c13-7add-4b54-988d-474dfff0db35.jpg)

100

90

80

70

60

50

40

30

20

10

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-1e90d43c-007c-4807-aaa9-7c9468ec441f.jpg)

FIG. 7: The same as Figure 2 but after the SES at KER station on 8 December 2020 upon the occurrence of the ML(ATH)=2.1 EQ on 18 March 2021 with an epicenter at 38.16 N22.91 E by considering M o o thres =2.1 (A) and 2.2 (B), respectively.

FIG. 8: The same as Figure 2 but after the SES at KER station on 8 December 2020, for M thres =2.3 (A) and 2.4 (B), upon the occurrence of the ML(ATH)=2.8 EQ at 08:02 UTC on 24 March 2021 with an epicenter at 38.30 N23.68 E. o o

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-06c834fd-bdb8-4395-a0c1-d2c6de0a1d2f.jpg)

FIG. 9: The same as Figure 2 but after the SES at KER station on 8 December 2020 upon the occurrence of the ML(ATH)=3.5 EQ at 06:15 UTC on 25 March 2021 with an epicenter at 38.76 o N23.41 E for M o thres =3.0. A principal peak at κ 1 =0.070 is also observed for M thres =3.1 and 3.2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-7f0a8e53-9c4b-423a-bec0-590a16f7f5c8.jpg)

FIG. 10: The SES activity recorded at KER geolectrical station on 2 June 2021.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-6885b19d-9b72-4f8c-ab94-20cbcd2add18.jpg)

FIG. 11: The SES activity recorded at KER geolectrical station on 15 November 2021. Only the intense colors correspond to the current data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-e8b0f080-3e2e-4a22-b2d5-8679f87f0272.jpg)

)

κ

1

Prob(

A

0.5

0.4

0.3

0.2

0.1

0.00

0.02

0.04

0.06

κ

1

0.08

0.10

0.12

100

110

40

70

60

No of EQs after SES

0.14

90

120

130

FIG. 12: The same as Figure 2 but after the SES at KER station on 15 November 2021 (Fig.11) upon the occurrence of the ML(ATH)=2.4 EQ at 17:13 UTC on 25 December 2021 with an epicenter at 38.31 o N23.41 E for M o thres =1.7(A) and 1.8(B).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-8c5ad6b5-3c12-45fd-94a8-7703a2e59924.jpg)

80

50

30

20

10

FIG. 13: The same as Figure 2 but after the SES at KER station on 15 November 2021 (Fig.11) upon the occurrence of the ML(ATH)=2.1 EQ at 20:39 UTC on 26 December 2021 with an epicenter at 38.52 o N23.95 E for M o thres =1.7(A) and 1.8(B).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-6a32a235-1934-4bcf-a836-3aaca1d1c45f.jpg)

FIG. 14: The same as Figure 2 but after the SES at KER station on 15 November 2021 (Fig.11) upon the occurrence of the ML(ATH)=3.3 EQ at 17:19 UTC on 25 February 2022 with an epicenter at 38.49 N23.61 E for M o o thres =2.0(a), 2.1(b), 2.2(c), 2.3(d), 2.4(e), and 2.5(f).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-cd9dc518-8827-4c0c-8358-4aa4e8540fb9.jpg)

FIG. 15: The same as Figure 2 but after the SES at KER station on 15 November 2021 (Fig.11) for M thres =2.3(a), 2.4(b), 2.5(c), 2.6(d), 2.7(e), 2.8(f), 2.9(g), 3.0(h), 3.1(i) and 3.2(j) upon the occurrence of the EQs: ML(ATH)=2.5 at 21:27 UTC on 7 April 2022 for (a) and (b), ML(ATH)=3.0 at 11:08 UTC on 8 April 2022 for (c), ML(ATH)=3.5 at 9:43 UTC on 8 April 2022 for (d) and (g), ML(ATH)=2.9 at 17:06 UTC on 5 April 2022 for (e), ML(ATH)=2.9 at 9:04 UTC on 8 April 2022 for (f), ML(ATH)=3.5 at 15:46 UTC on 7 April 2022 for (h), ML(ATH)=3.1 at 0:55 UTC on 7 April 2022 for (i), and ML(ATH)=3.8 at 11:36 UTC on 7 April 2022 for (j).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-9289380d-b3c1-466c-9904-16b1ecec6c2d.jpg)

Note added on 1 August 2019. In the main text of the previous version of this paper[49], it has been reported that a pronounced M (USGS)=5.4 w earthquake (EQ) -or ML(ATH)=5.2 EQ, thus Ms(ATH)=ML(ATH)+0.5=5.7 EQ -which was strongly felt at Athens, Greece at 23:05 UTC on 17 November 2014 with an epicenter at 38.64 o N 23.40 E has been preo ceded by an SES activity at Keratea (KER) geoelectrical station uploaded in the arXiv [46] almost three months before, i.e., on 7 August 2014 (only if the expected EQ magnitude Ms(ATH) estimated from the amplitude of the SES activity is comparable with or larger than 6.0, quick report is uploaded before the EQ occurrence as explained in the subsection 7.2 of Ref. [28]). It has been followed by an EQ of equal magnitude (ML=5.2) almost four minutes later, i.e., at 23:09 UTC on 17 November 2014, practically at the same epicenter, i.e., at 38.64 o N 23.41 E. o

Here, we report that at 11:13 UTC on 19 July 2019, a M (USGS)=5.3 EQ -or ML(ATH)=5.1- was also stongly w felt at Athens with an epicenter at 38.12 N 23.53 E. It o o has been preceded by an SES activity at KER which can be seen in Fig.3. An inspection of this figure shows that it had a different ratio of the SES components compared to the SES activity depicted in Fig.1(a) that preceded the previous EQ in 2014 mentioned above. This explains

- [1] D. Sornette, Critical Phenomena in the Natural Sciences: Chaos, Fractals, Selforganization, and Disorder: Concepts and Tools (Springer-Verlag, Berlin Heidelberg, 2000).
- [2] A. Corral, Phys. Rev. Lett. 92 , 108501 (2004).
- [3] J. Davidsen and M. Paczuski, Phys. Rev. Lett. 94 , 048501 (2005).
- [4] J. R. Holliday, J. B. Rundle, D. L. Turcotte, W. Klein, K. F. Tiampo, and A. Donnellan, Phys. Rev. Lett. 97 , 238501 (2006).
- [5] A. Saichev and D. Sornette, Phys. Rev. Lett. 97 , 078501 (2006).
- [6] J. F. Eichner, J. W. Kantelhardt, A. Bunde, and S. Havlin, Phys. Rev. E 75 , 011128 (2007).
- [7] S. Lennartz, V. N. Livina, A. Bunde, and S. Havlin, EPL (Europhysics Letters) 81 , 69001 (2008).
- [8] S. Lennartz, A. Bunde, and D. L. Turcotte, Geophys. J. Int. 184 , 1214 (2011).
- [9] E. Lippiello, L. de Arcangelis, and C. Godano, Phys. Rev. Lett. 103 , 038501 (2009).
- [10] E. Lippiello, C. Godano, and L. de Arcangelis, Geophys. Res. Lett. 39 , L05309 (2012).
- [11] N. V. Sarlis, E. S. Skordas, and P. A. Varotsos, Phys. Rev. E 80 , 022102 (2009).
- [12] N. V. Sarlis, E. S. Skordas, and P. A. Varotsos, Phys. Rev. E 82 , 021110 (2010).
- [13] L. Telesca and M. Lovallo, Geophys. Res. Lett. 36 , L01308 (2009).
- [14] M. Bottiglieri, L. de Arcangelis, C. Godano, and E. Lippiello, Phys. Rev. Lett. 104 , 158501 (2010).
- [15] L. Telesca, Physica A 389 , 1911 (2010).

why (e.g., see Varotsos and Lazaridou [24]) the recent 19 July 2019 EQ occurred at a different region of the SES selectivity map of the measuring station depicted in Fig. 1(b). Another important difference between these two cases is that the recent event was followed almost 1 hour later, i.e., at 12:11 UTC on 19 July 2019, approximately at the same epicenter, i.e., at 38.10 o N 23.58 E, o by a smaller EQ with ML=4.3. The study of the evolution of the seismic activity is made by analysing in natural time the events occurring in the candidate area (Fig. 1(b)) subsequent to the SES activity recorded at KER on 19 June 2019 (Fig. 3) by means of the procedure developed in Section II. In view of the complexity of this procedure, which is the most accurate, one may alternatively rely -but only approximately- on the upper time chart depicted in Fig.28 of Ref.[24], which is explained in simple words in p.35 of Varotsos et al. [51].

## Acknowledgments

We gratefully acknowledge the continuous supervision and technical support of the geoelectrical stations of the SES telemetric network by Vasilis Dimitropoulos, Spyros Tzigkos and George Lampithianakis.

- [16] N. V. Sarlis, Phys. Rev. E 84 , 022101 (2011).
- [17] Q. Huang and X. Ding, Bull. Seismol. Soc. Am. 102 , 1878 (2012).
- [18] N. V. Sarlis and S.-R. G. Christopoulos, Chaos 22 , 023123 (2012).
- [19] P. Varotsos, N. Sarlis, and E. Skordas, EPL 96 , 59002 (2011).
- [20] P. Varotsos, N. Sarlis, and E. Skordas, EPL 99 , 59001 (2012).
- [21] D. L. Turcotte, Fractals and Chaos in Geology and Geophysics (Cambridge University Press, Cambridge, 1997), 2nd ed.
- [22] P. Varotsos and K. Alexopoulos, Tectonophysics 110 , 73 (1984).
- [23] P. Varotsos and K. Alexopoulos, Tectonophysics 110 , 99 (1984).
- [24] P. Varotsos and M. Lazaridou, Tectonophysics 188 , 321 (1991).
- [25] P. Varotsos, K. Alexopoulos, and M. Lazaridou, Tectonophysics 224 , 1 (1993).
- [26] P. Varotsos, The Physics of Seismic Electric Signals (TERRAPUB, Tokyo, 2005).
- [27] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Chaos 19 , 023114 (2009).
- [28] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Natural Time Analysis: The new view of time. Precursory Seismic Electric Signals, Earthquakes and other Complex Time-Series (Springer-Verlag, Berlin Heidelberg, 2011).
- [29] P. Varotsos and K. Alexopoulos, Thermodynamics of Point Defects and their Relation with Bulk Properties (North Holland, Amsterdam, 1986).

- [30] P. A. Varotsos, N. V. Sarlis, E. S. Skordas, S. Uyeda, and M. Kamogawa, EPL 92 , 29002 (2010).
- [31] P. A. Varotsos, N. V. Sarlis, E. S. Skordas, and M. S. Lazaridou, Tectonophysics 589 , 116 (2013).

[32]

P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, J. Geo- phys. Res.:

Space Physics

119

,

9192 (2014).

- [33] N. V. Sarlis, E. S. Skordas, P. A. Varotsos, T. Nagao, M. Kamogawa, H. Tanaka, and S. Uyeda, Proc. Natl. Acad. Sci. USA 110 , 13734 (2013).
- [34] N. V. Sarlis, E. S. Skordas, P. A. Varotsos, T. Nagao, M. Kamogawa, and S. Uyeda, Proc. Natl. Acad. Sci. USA 112 , 986 (2015).
- [35] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Practica of Athens Academy 76 , 294 (2001).
- [36] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Phys. Rev. E 66 , 011902 (2002).
- [37] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Acta Geophysica Polonica 50 , 337 (2002).
- [38] P. A. Varotsos, N. V. Sarlis, H. K. Tanaka, and E. S. Skordas, Phys. Rev. E 72 , 041103 (2005).
- [39] N. V. Sarlis, E. S. Skordas, M. S. Lazaridou, and P. A. Varotsos, Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 84 , 331 (2008).
- [40] N. V. Sarlis, Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 89 , 438 (2013).
- [41] Q. Huang, H. Ren, D. Zhang, and Y. J. Chen, Proc. Jpn Acad. Ser. B Phys. Biol. Sci. 91 , 17 (2015).
- [42] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Phys. Rev. E 67 , 021109 (2003).
- [43] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, Phys. Rev. E 68 , 031106 (2003).
- [44] P. A. Varotsos, N. V. Sarlis, E. S. Skordas, and M. S. Lazaridou, J. Appl. Phys. 103 , 014906 (2008).
- [45] H. Kanamori, Nature 271 , 411 (1978).
- [46] N. V. Sarlis, P. V. Varotsos, and E. S. Skordas, arXiv e-prints arXiv:1301.7634v7 (2014), 1301.7634v7.
- [47] G. Chouliaras, N. S. Melis, G. Drakatos, and K. Makropoulos, Advances in Geosciences 36 , 7 (2013).
- [48] A. Mignan and G. Chouliaras, Seismological Research Letters 85 , 657 (2014).
- [49] P. A. Varotsos, N. V. Sarlis, E. S. Skordas, S.-R. G. Christopoulos, and M. S. Lazaridou-Varotsos, Earthquake Science 28 , 215 (pages 7) (2015).
- [50] M. Moustra, M. Avraamides, and C. Christodoulou, Expert Systems with Applications 38 , 15032 (2011).
- [51] P. Varotsos, M. Lazaridou, K. Eftaxias, G. Antonopoulos, J. Makris, and J. Kopanas, in The Critical Review of VAN: Earthquake Prediction from Seismic Electric Signals , edited by S. J.Lighthill (World Scientific, Singapore, 1996), pp. 29-76.
- [52] Note added on 8 June 2022: The evolution in the lower time chart of p. 65 of Ref.[51] describing the case of the M6.5 EQ that occurred almost 27 years ago (commented on by the Times of London on the next day of EQ occurrence, see p.142 of Ref. [53]) seems to be similar to that followed after a probable SES activity recorded at VOL on 3 May 2022.
- [53] Note added on 25 November 2022: In a recent work, Christopoulos et al. [ Appl. Sci. 12 , 7496 (2022)] by employing natural time analysis of the order parameter fluctuation minima β min showed that with statistical significant of the order of 10 -5 , the occurrence time of M ≥ 8 5 . EQs can be determined with a maximum lead time of 9 months with outstanding discrimination. This holds also

for the M 9 Tohoku EQ that occurred in Japan on 11 March 2011, which has been preceded by the deepest β min value in Japan on around 5 January 2011 that appeared almost simultaneously with an SES activity, i.e., by a lead time of around two months which is smaller than the maximum lead time of ≈ 5 1 2 months observed to date for SES activities in Greece, where -during the last 40 year operation of the VAN telemetric networkno EQ of magnitude significantly larger than M 7 has yet occurred. The latter does not preclude of course that, for appreciably larger EQs, we may observe appreciably larger values of the lead time of SES activities (i.e., up to 9 months, or so, as mentioned above). This reflects that we should extend the computation of the κ 1 value (for the identification of the occurrence time of an impending EQ with an accuracy of a few days to 1 week or so, as for example for the current case of our update on 8 June 2022, see also p. 64 of Ref. [51]) until around 9 months after the SES activities' recording.

- [54] Note added on 12 January 2023: An SES activity has been recorded at Assiros (ASS) station (close to Thessaloniki) on 8 January 2023 as shown in Fig. 16. The corresponding SES selectivity map of ASS station is depicted in Fig. 17.
- [55] Note added on 7 February 2023: On 2 February 2023, upon the occurrence of the ML(ATH)=2.1 EQ at 20:51 UTC with an epicenter at 38.99 o N24.18 E, we found that o the criticality condition κ 1 = 0 070 . has been obeyed inside the selectivity map of ASS station (Fig. 17) exhibiting magnitude threshold invariance in the range M thres =1.7 to M thres =1.9, which signals the approach of the system to the critical point.
- [56] Note added on 16 March 2023: In continuation, at 02:46 and 02:49 UTC on 16 March 2023, upon the occurrence of the two EQs of magnitude ML(ATH)=2.4 and ML(ATH)=2.3 with epicenters at 39.09 N23.18 E and o o 40.36 N24.07 E, respectively, o o and using the thresholds M thres =2.1, 2.3, and 2.4, it seems that the critical point is approached.
- [57] Note added on 8 May 2023: Remarkable consecutive facts (i.e., fulfillment of the criticality condition κ 1 = 0 070 . in the vicinity of the critical point) emerged during the last few days upon the occurrence of the following EQs: (A) on 5 May 2023 11:35 UTC at 39.98 N23.10 E, o o a ML(ATH)=2.6 for the thresholds M thres =2.1 and 2.4, (B) on 8 May 2023 03:48 UTC at 40.07 N24.46 E, o o a ML(ATH)=2.0 for the threshold M thres =2.0, and (C) on 8 May 2023 06:50 UTC at 40.96 N22.36 E, a o o ML(ATH)=2.3 for the thresholds M thres = 2 2 and 2.3. . These EQs occurred within the shaded regions c and a in Fig. 17. According to pages 1 and 13 of Ref.[26], the regions a, b, c refer to seismic areas from which SES activities have been recorded at ASS and seems probable that they are interconnected to a larger area as depicted in Fig.4.3.5b in p. 104 of Ref.[26]. Attention is drawn to the point that the latter event (C) just occurred in the area c.
- [58] Note added on 22 July 2023: Following Ref. [53], upon the occurrence of a ML(ATH)=3.8 EQ with an epicenter at 39.53 o N23.67 E, we find that the criticality condition o κ 1 = 0 070 has been obeyed for . a multitude of magnitude thresholds including M thres =2.4, 2.9, 3.1, 3.2, 3.3, and 3.4. It is challenging that at three of them, i.e., M thres =2.4, 3.1, and 3.2, the plot Prob( κ 1 ) versus κ 1

is maximized at κ 1 = 0 070 pointing to an imminent ap-. proach of the critical point, a fact that has been observed for the first time.

- [59] Note added on 11 August 2023: The aforementioned important fact observed for the first time on 22 July 2023, i.e., the plot Prob( κ 1 ) versus κ 1 is maximized at the ciriticality condition κ 1 = 0 070, . was again observed upon the occurrence of the following two events: the ML(ATH)=3.6 EQ on 10 August 2023 at 12:39:56 UTC with an epicenter at 39.51 N23.67 E o o and the ML(ATH)=3.7 EQ on 11 August 2023 at 13:20:56 UTC with an epicenter at 39.51 N23.69 E. o o We note that for the second case the criticality condition holds for M thres =3.1, 3.2, 3.3, 3.4, and 3.5 pointing to an imminent approach of the critical point.
- [60] M. S. Lazaridou-Varotsos, Earthquake Prediction by Seismic Electric Signals. The success of the VAN method over thirty years (Springer-Verlag, Berlin Heidelberg, 2013).
- [61] Note added on 5 June 2024: Almost one hour after the occurrence of an EQ felt at Athens of local magnitude ML(ATH)=2.7 with epicenter at 38.26 N23.86 E (2km o o SSW of Kalamos in Attica), bay like SES were recorded from approximately 22:30 UTC on 30 May 2024 until

02:40 UTC on 31 May 2024 at KER, see Fig. 18. Their feature is strikingly reminiscent of that observed again at Keratea station almost four decades ago in particular from approximately 19:52 UTC on 8 September 1986 until around 05:22 UTC on 9 September 1986 followed by the two destructive EQs (Ms(ATH)=6.2 and 5.6) on 13 September 1986 and 15 September 1986 whose epicenters were located almost 15km from Kalamata city. It should be noticed that other stations like the one located at Nafplion, see Fig.1 of Varotsos and Alexopoulos [22], placed at smaller epicentral distance, did not show any SES, thus indicating the extremely important role of SES selectivity effect as discussed by Varotsos and Alexopoulos [62].

- [62] P. Varotsos and K. Alexopoulos, Tectonophysics 136 , 335 (1987).
- [63] Note added on 30 August 2024: We clarify that a bay like SES depicted in Fig. 18 may precede any EQ located within the KER SES selectivity map (Fig.1b) of magnitude comparable with Ms(ATH) ≈ 6.0 that may occur within the time window of an SES activity from a few weeks to around 5 1 2 months or so.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-fe104f12-428a-4fe9-a5ea-c5bec3f446c1.jpg)

8-Jan-2023

FIG. 16: The SES activity recorded at ASS geolectrical station on 8 January 2023.

✁

✁

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-0ce48109-1944-43ce-af5f-6c1c3c05b29b.jpg)

✁

✁

✁

FIG. 17: Map showing the operating geoelectric VAN stations (blue triangles). The ASS selectivity map is bounded by the red dashed-dotted line according to Ref. [40].

FIG. 18: The SES recorded at channels A αντ -E 2 , A αντ -1 N , A αντ -∆ , 4 A 3 -∆ , 4 A αντ -A 3 of KER geolectrical station from 22:30 UTC on 30 May 2024 until 02:40 UTC on 31 May 2024.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00279v23-figure-88ee6786-4859-48a3-a75d-431d942044d5.jpg)
P. A. Varotsos N. V. Sarlis E. S. Skordas M. S. Lazaridou-Varotsos
Identifying the occurrence time of an impending mainshock: A very recent case
The procedure by means of which the occurrence time of an impending mainshock
can be identified by analyzing in natural time the seismicity in the candidate
area subsequent to the recording of a precursory Seismic Electric Signals (SES)
activity is reviewed. Here, we report the application of this procedure to an
Mw5.4 mainshock that occurred in Greece on 17 November 2014 and was strongly
felt in Athens. This mainshock (which is pretty rare since it is the strongest
in that area for more than half a century) was preceded by an SES activity
recorded on 27 July 2014 and the results of the natural time analysis reveal
that the system approached the critical point (mainshock occurrence) early in
the morning on 15 November 2014. SES activities that have been recently
recorded are also presented. Furthermore, in a Note we discuss the case of the
Mw5.3 earthquake that was also strongly felt in Athens on 19 July 2019
(Parnitha fault).