Universal Non-Extensive Statistical Physics Temporal Pattern of Major Subduction Zone Aftershock Sequences

Forecasting Strong Subsequent Earthquakes in Greece with the Machine Learning Algorithm NESTORE

Aftershocks of earthquakes can destroy many urban infrastructures and exacerbate the damage already inflicted upon weak structures. Therefore, it is important to have a method to forecast the probability of occurrence of stronger earthquakes in order to mitigate their effects. In this work, we applied the NESTORE machine learning approach to Greek seismicity from 1995 to 2022 to forecast the probability of a strong aftershock. Depending on the magnitude difference between the mainshock and the strongest aftershock, NESTORE classifies clusters into two types, Type A and Type B. Type A clusters are the most dangerous clusters, characterized by a smaller difference. The algorithm requires region-dependent training as input and evaluates performance on an independent test set. In our tests, we obtained the best results 6 h after the mainshock, as we correctly forecasted 92% of clusters corresponding to 100% of Type A clusters and more than 90% of Type B clusters. These results were also obtained thanks to an accurate analysis of cluster detection in a large part of Greece. The successful overall results show that the algorithm can be applied in this area. The approach is particularly attractive for seismic risk mitigation due to the short time required for forecasting.

Complexity of Recent Earthquake Swarms in Greece in Terms of Non-Extensive Statistical Physics

Greece exhibits the highest seismic activity in Europe, manifested in intense seismicity with large magnitude events and frequent earthquake swarms. In the present work, we analyzed the spatiotemporal properties of recent earthquake swarms that occurred in the broader area of Greece using the Non-Extensive Statistical Physics (NESP) framework, which appears suitable for studying complex systems. The behavior of complex systems, where multifractality and strong correlations among the elements of the system exist, as in tectonic and volcanic environments, can adequately be described by Tsallis entropy (S_q), introducing the Q-exponential function and the entropic parameter q that expresses the degree of non-additivity of the system. Herein, we focus the analysis on the 2007 Trichonis Lake, the 2016 Western Crete, the 2021-2022 Nisyros, the 2021-2022 Thiva and the 2022 Pagasetic Gulf earthquake swarms. Using the seismicity catalogs for each swarm, we investigate the inter-event time (T) and distance (D) distributions with the Q-exponential function, providing the q_T and q_D entropic parameters. The results show that q_T varies from 1.44 to 1.58, whereas q_D ranges from 0.46 to 0.75 for the inter-event time and distance distributions, respectively. Furthermore, we describe the frequency-magnitude distributions with the Gutenberg-Richter scaling relation and the fragment-asperity model of earthquake interactions derived within the NESP framework. The results of the analysis indicate that the statistical properties of earthquake swarms can be successfully reproduced by means of NESP and confirm the complexity and non-additivity of the spatiotemporal evolution of seismicity. Finally, the superstatistics approach, which is closely connected to NESP and is based on a superposition of ordinary local equilibrium statistical mechanics, is further used to discuss the temporal patterns of the earthquake evolution during the swarms.