Earthquake Nowcasting: Retrospective Testing in Greece.

Earthquake nowcasting (EN) is a modern method of estimating seismic risk by evaluating the progress of the earthquake (EQ) cycle in fault systems. EN evaluation is based on a new concept of time, termed 'natural time'. EN employs natural time, and uniquely estimates seismic risk by means of the earthquake potential score (EPS), which has been found to have useful applications both regionally and globally. Amongst these applications, here we focused on Greece since 2019, for the estimation of the EPS for the largest-magnitude events, MW(USGS) ≥ 6, that occurred during our study period: for example, the MW= 6.0 WNW-of-Kissamos EQ on 27 November 2019, the MW= 6.5 off-shore Southern Crete EQ on 2 May 2020, the MW= 7.0 Samos EQ on 30 October 2020, the MW= 6.3 Tyrnavos EQ on 3 March 2021, the MW= 6.0 Arkalohorion Crete EQ on 27 September 2021, and the MW= 6.4 Sitia Crete EQ on 12 October 2021. The results are promising, and reveal that the EPS provides useful information on impending seismicity.

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

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

Natural Time Analysis: The Area under the Receiver Operating Characteristic Curve of the Order Parameter Fluctuations Minima Preceding Major Earthquakes.

It has been reported that major earthquakes are preceded by Seismic Electric Signals (SES). Observations show that in the natural time analysis of an earthquake (EQ) catalog, an SES activity starts when the fluctuations of the order parameter of seismicity exhibit a minimum. Fifteen distinct minima-observed simultaneously at two different natural time scales and deeper than a certain threshold-are found on analyzing the seismicity of Japan from 1 January 1984 to 11 March 2011 (the time of the M9 Tohoku EQ occurrence) 1 to 3 months before large EQs. Six (out of 15) of these minima preceded all shallow EQs of magnitude 7.6 or larger, while nine are followed by smaller EQs. The latter false positives can be excluded by a proper procedure (J. Geophys. Res. Space Physics 2014, 119, 9192-9206) that considers aspects of EQ networks based on similar activity patterns. These results are studied here by means of the receiver operating characteristics (ROC) technique by focusing on the area under the ROC curve (AUC). If this area, which is currently considered an effective way to summarize the overall diagnostic accuracy of a test, has the value 1, it corresponds to a perfectly accurate test. Here, we find that the AUC is around 0.95 which is evaluated as outstanding.

Identifying the Occurrence Time of the Deadly Mexico M8.2 Earthquake on 7 September 2017.

It has been shown that some dynamic features hidden in the time series of complex systems can be unveiled if we analyze them in a time domain termed natural time. In this analysis, we can identify when a system approaches a critical point (dynamic phase transition). Here, based on natural time analysis, which enables the introduction of an order parameter for seismicity, we discuss a procedure through which we could achieve the identification of the occurrence time of the M8.2 earthquake that occurred on 7 September 2017 in Mexico in Chiapas region, which is the largest magnitude event recorded in Mexico in more than a century. In particular, we first investigated the order parameter fluctuations of seismicity in the entire Mexico and found that, during an almost 30-year period, i.e., from 1 January 1988 until the M8.2 earthquake occurrence, they were minimized around 27 July 2017. From this date, we started computing the variance of seismicity in Chiapas region and found that it approached the critical value 0.070 on 6 September 2017, almost one day before this M8.2 earthquake occurrence.

Spatiotemporal variations of seismicity before major earthquakes in the Japanese area and their relation with the epicentral locations.

Using the Japan Meteorological Agency earthquake catalog, we investigate the seismicity variations before major earthquakes in the Japanese region. We apply natural time, the new time frame, for calculating the fluctuations, termed ß, of a certain parameter of seismicity, termed κ1. In an earlier study, we found that ß calculated for the entire Japanese region showed a minimum a few months before the shallow major earthquakes (magnitude larger than 7.6) that occurred in the region during the period from 1 January 1984 to 11 March 2011. In this study, by dividing the Japanese region into small areas, we carry out the ß calculation on them. It was found that some small areas show ß minimum almost simultaneously with the large area and such small areas clustered within a few hundred kilometers from the actual epicenter of the related main shocks. These results suggest that the present approach may help estimation of the epicentral location of forthcoming major earthquakes.

Study in Natural Time of Geoelectric Field and Seismicity Changes Preceding the Mw6.8 Earthquake on 25 October 2018 in Greece.

A strong earthquake of magnitude M w 6.8 struck Western Greece on 25 October 2018 with an epicenter at 37.515 ∘ N 20.564 ∘ E. It was preceded by an anomalous geolectric signal that was recorded on 2 October 2018 at a measuring station 70 km away from the epicenter. Upon analyzing this signal in natural time, we find that it conforms to the conditions suggested for its identification as precursory Seismic Electric Signal (SES) activity. Notably, the observed lead time of 23 days lies within the range of values that has been very recently identified as being statistically significant for the precursory variations of the electric field of the Earth. Moreover, the analysis in natural time of the seismicity subsequent to the SES activity in the area candidate to suffer this strong earthquake reveals that the criticality conditions were obeyed early in the morning of 18 October 2018, i.e., almost a week before the strong earthquake occurrence, in agreement with earlier findings. Finally, when employing the recent method of nowcasting earthquakes, which is based on natural time, we find an earthquake potential score around 80%.

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

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.

Estimating the Compressibility of Osmium from Recent Measurements of Ir-Os Alloys under High Pressure.

Several fcc- and hcp-structured Ir-Os alloys have been recently studied up to 30 GPa at room temperature by means of synchrotron-based X-ray powder diffraction in diamond anvil cells. Using their bulk moduli, which increase with increasing osmium content, showing a deviation from linearity, and after employing a thermodynamical model, it was concluded that the bulk modulus for osmium is slightly smaller than that for diamond. Here, a similar conclusion is obtained upon employing an alternative model, thus strengthening the conclusion that osmium is the densest but not the most incompressible element. This is particularly interesting for Earth Sciences because it may be of key importance toward clarifying the anomalous elastic properties of the Earth's core.

Minimum of the order parameter fluctuations of seismicity before major earthquakes in Japan.

It has been shown that some dynamic features hidden in the time series of complex systems can be uncovered if we analyze them in a time domain called natural time χ. The order parameter of seismicity introduced in this time domain is the variance of χ weighted for normalized energy of each earthquake. Here, we analyze the Japan seismic catalog in natural time from January 1, 1984 to March 11, 2011, the day of the M9 Tohoku earthquake, by considering a sliding natural time window of fixed length comprised of the number of events that would occur in a few months. We find that the fluctuations of the order parameter of seismicity exhibit distinct minima a few months before all of the shallow earthquakes of magnitude 7.6 or larger that occurred during this 27-y period in the Japanese area. Among the minima, the minimum before the M9 Tohoku earthquake was the deepest. It appears that there are two kinds of minima, namely precursory and nonprecursory, to large earthquakes.

Natural time analysis of critical phenomena.

A quantity exists by which one can identify the approach of a dynamical system to the state of criticality, which is hard to identify otherwise. This quantity is the variance κ(1)(≡<χ(2)> - <χ>(2)) of natural time χ, where = Σp(k)f(χ(k)) and p(k) is the normalized energy released during the kth event of which the natural time is defined as χ(k) = k/N and N stands for the total number of events. Then we show that κ(1) becomes equal to 0.070 at the critical state for a variety of dynamical systems. This holds for criticality models such as 2D Ising and the Bak-Tang-Wiesenfeld sandpile, which is the standard example of self-organized criticality. This condition of κ(1) = 0.070 holds for experimental results of critical phenomena such as growth of rice piles, seismic electric signals, and the subsequent seismicity before the associated main shock.

Micro-scale, mid-scale, and macro-scale in global seismicity identified by empirical mode decomposition and their multifractal characteristics.

The magnitude time-series of the global seismicity is analyzed by the empirical mode decomposition giving rise to 14 intrinsic mode functions (IMF) and a trend. Using Hurst analysis one can identify three different sums of these IMFs and the trend which exhibit distinct multifractal behaviour and correspond to micro-, mid- and macro-scales. Their multifractal detrended fluctuation analysis reveals that the micro-scale time-series exhibits anticorrelated behaviour in contrast to the mid-scale one which is long-range correlated. Concerning the mid-scale one, in the range of 30 to 300 consecutive events the maximum entropy method power spectra indicates that it exhibits an 1/f α behaviour with α close to 1/3 which is compatible with the long-range correlations identified by detrended fluctuation analysis during periods of stationary seismicity. The results have been also verified to hold regionally for the earthquakes in Japan and shed light on the significance of the mid-scale of 30 to 300 events in the natural time analysis of global (and regional) seismicity. It is shown that when using the mid-scale time-series only, we can obtain results similar to those obtained by the natural time analysis of global seismicity when focusing on the prediction of earthquakes with M ≥ 8.4.