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.

On 2 December 2020, an earthquake with a magnitude of M_w 4.5 occurred near the city of Thiva (Greece). The aftershock sequence, triggered by ruptures on or near the Kallithea fault, continued until January 2021. Seven months later, new seismic activity began a few kilometers west of the initial events, with the swarm displaying a general trend of spatiotemporal migration toward the east–southeast until the middle of 2022. In order to understand the physical and statistical pattern of the swarm, the seismicity was relocated using HypoDD, and the magnitude of completeness was determined using the frequency–magnitude distribution. In order to define the existence of spatiotemporal seismicity clusters in an objective way, the DBSCAN clustering algorithm was applied to the 2020–2022 Thiva earthquake sequence. The extracted clusters permit the analysis of the spatiotemporal scaling properties of the main clusters using the Non-Extensive Statistical Physics (NESP) approach, providing detailed insights into the nature of the long-term correlation of the seismic swarm. The statistical pattern observed aligns with a Q-exponential distribution, with q_D values ranging from 0.7 to 0.8 and q_T values from 1.44 to 1.50. Furthermore, the frequency–magnitude distributions were analyzed using the fragment–asperity model proposed within the NESP framework, providing the non-additive entropic parameter (q_M). The results suggest that the statistical characteristics of earthquake clusters can be effectively interpreted using NESP, highlighting the complexity and non-additive nature of the spatiotemporal evolution of seismicity. In addition, the analysis of the properties of the seismicity clusters extracted using the DBSCAN algorithm permits the suggestion of possible physical mechanisms that drive the evolution of the two main and larger clusters. For the cluster that activated first and is located in the west–northwest part, an afterslip mechanism activated after the 2 September 2021, M 4.0 events seems to predominately control its evolution, while for the second activated cluster located in the east–southeast part, a normal diffusion mechanism is proposed to describe its migration pattern. Concluding, we can state that in the present work the application of the DBSCAN algorithm to recognize the existence of any possible spatiotemporal clustering of seismicity could be helping to provide detailed insight into the statistical and physical patterns in earthquake swarms.

Significant seismic activity has been witnessed in the area of Ridgecrest (Southern California) over the past 40 years, with the largest being the M_w 5.8 event on 20 September 1995. In July 2019, a strong earthquake of M_w 7.1, preceded by a M_w 6.4 foreshock, impacted Ridgecrest. The mainshock triggered thousands of aftershocks that were thoroughly documented along the activated faults. In this study, we analyzed the spatiotemporal variations of the frequency–magnitude distribution in the area of Ridgecrest using the fragment–asperity model derived within the framework of non-extensive statistical physics (NESP), which is well-suited for investigating complex dynamic systems with scale-invariant properties, multi-fractality, and long-range interactions. Analysis was performed for the entire duration, as well as within various time windows during 1981–2022, in order to estimate the q_M parameter and to investigate how these variations are related to the dynamic evolution of seismic activity. In addition, we analyzed the spatiotemporal q_M value distributions along the activated fault zone during 1981–2019 and during each month after the occurrence of the M_w 7.1 Ridgecrest earthquake. The results indicate a significant increase in the q_M parameter when large-magnitude earthquakes occur, suggesting the system’s transition in an out-of-equilibrium phase and its preparation for seismic energy release.