Order Parameter and Entropy of Seismicity in Natural Time before Major Earthquakes: Recent Results

Comparing prediction efficiency in the BTW and Manna sandpiles

The state-of-the-art in the theory of self-organized criticality reveals that a certain inactivity precedes extreme events, which are located on the tail of the event probability distribution with respect to their sizes. The existence of the inactivity allows for the prediction of these events in advance. In this work, we explore the predictability of the Bak-Tang-Wiesenfeld (BTW) and Manna models on the square lattice as a function of the lattice length. For both models, we use an algorithm that forecasts the occurrence of large events after a fall in activity. The efficiency of the prediction can be universally described in terms of the event size divided by an appropriate power-law function of the lattice length. The power-law exponents are projected to be 2.75 and 3 for the Manna and BTW models respectively. The scaling with the exponent 2.75 is known for collapsing of the entire size-frequency relationship in the Manna model. However, the correspondence between events on different lattices in the BTW model requires a variety of exponents where 3 is the largest. This indicates that in thermodynamic limit, prediction does exist in the Manna but not in the BTW model, at least based on inactivity. The difference in the universality classes may underline the difference in the prediction.

A Comprehensive Numerical Model for Reservoir-Induced Earthquake Risk Assessment

The assessment of seismic risk and the prevention of earthquake occurrences during reservoir operation present significant challenges in terms of accurate determination. This study aims to address this issue by developing a numerical model. The primary objective is to estimate the vulnerability of different fault types to reservoir impoundment. This model integrates essential parameters such as fault dip and the relative orientation between the reservoir and potential earthquakes, and it is structured within a risk framework using polar coordinates. Through comprehensive computations, we evaluate the alterations in elastic stress and fluid pore pressure resulting from water impoundment. This is achieved by employing a fully coupled two-dimensional poroelastic approach. Furthermore, our model incorporates relevant seismic data to enhance its accuracy. The findings of our study underscore that the critical factor lies in the fault's precise positioning with respect to the reservoir. The risk associated with a fault is contingent upon both its location and its orientation, emphasizing the importance of these factors in determining hazardous zones.

Complexity in Geophysical Time Series of Strain/Fracture at Laboratory and Large Dam Scales: Review

One of the interesting directions of complexity theory is the investigation of the synchronization of mechanical behavior of large-scale systems by weak forcing, which is one of manifestations of nonlinearity/complexity of a system. The effect of periodic weak mechanical or electromagnetic forcing leading to synchronization was studied on the laboratory load-spring system as well as on a big dam's strain data. Due to synchronization, the phase space structure of the forced system strongly depends on the weak forcing intensity-determinism show itself in the recurrence of definite states of the forced system. The nonlinear dynamics of tilts/strains/seismicity near grand dams reflect both the complexity of the mentioned time series, connected with the natural agents (regional and local geodynamics), which were presented even before dam erection, as well as the effects of the water level (WL) variation in the reservoir, which is a quasi-periodic forcing superimposed on the natural geodynamic background. Both these effects are documented by the almost half-century of observations at the large Enguri Dam. The obtained data on the dynamics of strain/seismicity near a large dam can be used for the assessment of the possible risks, connected with the abrupt change of routine dynamics of construction.

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.

Natural Time Analysis of Global Seismicity

Natural time analysis enables the introduction of an order parameter for seismicity, which is just the variance of natural time χ, κ1=⟨χ2⟩−⟨χ⟩2. During the last years, there has been significant progress in the natural time analysis of seismicity. Milestones in this progress are the identification of clearly distiguishable minima of the fluctuations of the order parameter κ1 of seismicity both in the regional and global scale, the emergence of an interrelation between the time correlations of the earthquake (EQ) magnitude time series and these minima, and the introduction by Turcotte, Rundle and coworkers of EQ nowcasting. Here, we apply all these recent advances in the global seismicity by employing the Global Centroid Moment Tensor (GCMT) catalog. We show that the combination of the above three milestones may provide useful precursory information for the time of occurrence and epicenter location of strong EQs with M≥8.5 in GCMT. This can be achieved with high statistical significance (p-values of the order of 10−5), while the epicentral areas lie within a region covering only 4% of that investigated.