b map evaluation and on-fault stress state for the Antakya 2023 earthquakes.

The analysis of on-fault seismicity can enlighten the current stress state on the fault itself. Its definition is relevant to individuate fault patches that have not released all the accumulated stress even after the occurrence of a high magnitude earthquake. We use the b value to characterize the stress state on the fault of the Antakya 2023 main events, being b inversely proportional to the stress. The small magnitude seismicity occurring on the maximum slip fault-patches does not allow the b value estimation. This represents a strong indication that the maximum slip zone released most of the stress previously accumulated. Conversely, the lowest b values are located at the bends of the faults and close to the nucleation zone suggesting that, there, still exists not released stress implying that it could be reactivated in the future.

Statistics of seismicity to investigate the Campi Flegrei caldera unrest.

The knowledge of the dynamic of the Campi Flegrei calderic system is a primary goal to mitigate the volcanic risk in one of the most densely populated volcanic areas in the world. From 1950 to 1990 Campi Flegrei suffered three bradyseismic crises with a total uplift of 4.3 m. After 20 years of subsidence, the uplift started again in 2005 accompained by a low increment of the seismicity rate. In 2012 an increment in the seismic energy release and a variation in the gas composition of the fumaroles of Solfatara (in the central area of the caldera) were recorded. Since then, a slow and progressive increase in phenomena continued until today. We analyze the INGV - Osservatorio Vesuviano seismic catalogue of Campi Flegrei from 2000 to 2020 in order to look for any variation in the seismic parameters and compare them with geochemical monitored ones. A remarkable correlation between independent variables of earthquake cumulative number, CO/CO2 values and vertical ground deformation reveals a likely common origin. Moreover the correlation between all the variables here analysed enlightens that the same origin can cause the temporal behavior of all these variables. We interpret the seismological, geochemical and geodetic observable in terms of the injection of magmatic fluids into the hydrothermal system or its pressurization.

Multiple-time scaling and universal behavior of the earthquake interevent time distribution.

The interevent time distribution characterizes the temporal occurrence in seismic catalogs. Universal scaling properties of this distribution have been evidenced for entire catalogs and seismic sequences. Recently, these universal features have been questioned and some criticisms have been raised. We investigate the existence of universal scaling properties by analyzing a Californian catalog and by means of numerical simulations of an epidemic-type model. We show that the interevent time distribution exhibits a universal behavior over the entire temporal range if four characteristic times are taken into account. The above analysis allows us to identify the scaling form leading to universal behavior and explains the observed deviations. Furthermore, it provides a tool to identify the dependence on the mainshock magnitude of the c parameter that fixes the onset of the power law decay in the Omori law.

Forecasting of the first hour aftershocks by means of the perceived magnitude.

The majority of strong earthquakes takes place a few hours after a mainshock, promoting the interest for a real time post-seismic forecasting, which is, however, very inefficient because of the incompleteness of available catalogs. Here we present a novel method that uses, as only information, the ground velocity recorded during the first 30 min after the mainshock and does not require that signals are transferred and elaborated by operational units. The method considers the logarithm of the mainshock ground velocity, its peak value defined as the perceived magnitude and the subsequent temporal decay. We conduct a forecast test on the nine M  ≥ 6 mainshocks that have occurred since 2013 in the Aegean area. We are able to forecast the number of aftershocks recorded during the first 3 days after each mainshock with an accuracy smaller than 18% in all cases but one with an accuracy of 36%.

On-off intermittency in earthquake occurrence.

The clustered occurrence of earthquakes is viewed as an intermittent phenomenon, interpreting the clusters of events as chaotic bursts combined to the Poissonian occurrence of background seismicity. In particular, we suggest that it can be interpreted as an example of on-off intermittency. This kind of intermittency is parameter driven and exhibits certain universal statistical properties. The study of a Californian catalogue allows to interpret earthquake occurrence as an on-off intermittent phenomenon. Our results suggest the existence of a branching mechanism in earthquake occurrence well explained by epidemic type models.

Mechanical origin of aftershocks.

Aftershocks are the most striking evidence of earthquake interactions and the physical mechanisms at the origin of their occurrence are still intensively debated. Novel insights stem from recent results on the influence of the faulting style on the aftershock organisation in magnitude and time. Our study shows that the size of the aftershock zone depends on the fault geometry. We find that positive correlations among parameters controlling aftershock occurrence in time, energy and space are a stable feature of seismicity independently of magnitude range and geographic areas. We explain the ensemble of experimental findings by means of a description of the Earth Crust as an heterogeneous elastic medium coupled with a Maxwell viscoelastic asthenosphere. Our results show that heterogeneous stress distribution in an elastic layer combined with a coupling to a viscous flow are sufficient ingredients to describe the physics of aftershock triggering.

Spatial organization of foreshocks as a tool to forecast large earthquakes.

An increase in the number of smaller magnitude events, retrospectively named foreshocks, is often observed before large earthquakes. We show that the linear density probability of earthquakes occurring before and after small or intermediate mainshocks displays a symmetrical behavior, indicating that the size of the area fractured during the mainshock is encoded in the foreshock spatial organization. This observation can be used to discriminate spatial clustering due to foreshocks from the one induced by aftershocks and is implemented in an alarm-based model to forecast m > 6 earthquakes. A retrospective study of the last 19 years Southern California catalog shows that the daily occurrence probability presents isolated peaks closely located in time and space to the epicenters of five of the six m > 6 earthquakes. We find daily probabilities as high as 25% (in cells of size 0.04 × 0.04deg(2)), with significant probability gains with respect to standard models.

Role of static stress diffusion in the spatiotemporal organization of aftershocks.

We investigate the spatial distribution of aftershocks, and we find that aftershock linear density exhibits a maximum that depends on the main shock magnitude, followed by a power law decay. The exponent controlling the asymptotic decay and the fractal dimensionality of epicenters clearly indicate triggering by static stress. The nonmonotonic behavior of the linear density and its dependence on the main shock magnitude can be interpreted in terms of diffusion of static stress. This is supported by the power law growth with exponent H approximately 0.5 of the average main-aftershock distance. Implementing static stress diffusion within a stochastic model for aftershock occurrence, we are able to reproduce aftershock linear density spatial decay, its dependence on the main shock magnitude, and its evolution in time.

Influence of time and space correlations on earthquake magnitude.

A crucial point in the debate on the feasibility of earthquake predictions is the dependence of an earthquake magnitude from past seismicity. Indeed, while clustering in time and space is widely accepted, much more questionable is the existence of magnitude correlations. The standard approach generally assumes that magnitudes are independent and therefore in principle unpredictable. Here we show the existence of clustering in magnitude: earthquakes occur with higher probability close in time, space, and magnitude to previous events. More precisely, the next earthquake tends to have a magnitude similar but smaller than the previous one. A dynamical scaling relation between magnitude, time, and space distances reproduces the complex pattern of magnitude, spatial, and temporal correlations observed in experimental seismic catalogs.

Universality in solar flare and earthquake occurrence.

Earthquakes and solar flares are phenomena involving huge and rapid releases of energy characterized by complex temporal occurrence. By analyzing available experimental catalogs, we show that the stochastic processes underlying these apparently different phenomena have universal properties. Namely, both problems exhibit the same distributions of sizes, interoccurrence times, and the same temporal clustering: We find after flare sequences with power law temporal correlations as the Omori law for seismic sequences. The observed universality suggests a common approach to the interpretation of both phenomena in terms of the same driving physical mechanism.