New insights into crustal structure, Cenozoic magmatism, CO2 degassing and seismogenesis in the southern Apennines and Irpinia region from local earthquake tomography

We present high-resolution Vp and Vp/Vs models of the southern Apennines (Italy) computed using local earthquakes recorded from 2006 to 2011 with a graded inversion scheme that progressively resolves the crustal structure, from the large scale of the Apennines belt to the local scale of the normal-fault system. High-Vp bodies defined in the upper and mid crust under the external Apennines are interpreted as extensive mafic intrusions revealing anorogenic magmatism episodes that broadened on the Adriatic domain during Paleogene. Under the mountain belt, a low-Vp region, annular to the Neapolitan volcanic district, indicates the existence of a thermal/fluid anomaly in the mid crust, coinciding with a shallow Moho and diffuse degassing of deeply derived CO2. In the belt axial zone, low Vp/Vs gas-pressurized rock volumes under the Apulian carbonates correlate to high heat flow, strong CO2-dominated gas emissions of mantle origin and shallow carbonate reservoirs with pressurized CO2 gas caps. We hypothesize that the pressurized fluid volumes located at the base of the active fault system influence the rupture process of large normal-faulting earthquakes, like the 1980 Mw6.9 Irpinia event, and that major asperities are confined within the high-Vp Apulian carbonates. This study confirms once more that pre-existing structures of the Pliocene Apulian belt controlled the rupture propagation during the Irpinia earthquake. The main shock broke a 30 km long, NE-dipping seismogenic structure, whereas delayed ruptures (both the 20 s and the 40 s sub-events) developed on antithetic faults, reactivating thrust faults located at the eastern edge of the Apulian belt

A new view of Italian seismicity using 20 years of instrumental recordings

Abstract In this paper, we show the seismicity of the past 20 years that occurred in Italy and surrounding regions. Hypocentral locations have been obtained by using P- and S-wave arrival times from the INGV national and several regional permanent seismic networks. More than 48,000 events, selected from an original data set of about 99,780, are used to reconstruct the most complete seismic picture of the Italian region so far. The seismicity distribution allows inference on seismotectonics of this complex region of subduction versus continental collision. Our results clearly reveal the geometry of the Adria and the Ionian subduction and a continuous normal fault belt in the upper crust, following the Apennines mountain range. The depth of the seismogenic layer is computed from the cut-off of seismicity at depth and shows large variations along and across the seismic active regions. Earthquakes are generated by the different velocity of slab retreat and the subsequent asthenospheric upwelling. D 2004 Elsevier B.V. All rights reserved

Heterogeneities along the 2009 L’Aquila normal fault inferred by the b-value distribution

In this study we map the distribution of the b-value of the Gutenberg-Richter law—as well as complementary seismicity parameters—along the fault responsible for the 2009 MW 6.1 L'Aquila earthquake. We perform the calculations for two independent aftershock sub-catalogs, before and after a stable magnitude of completeness is reached. We find a substantial spatial variability of the b-values, which range from 0.6 to 1.3 over the fault plane. The comparison between the spatial distribution of the b-values and the main-shock slip pattern shows that the largest slip occurs in normal-to-high b-values portion of the fault plane, while low b-value is observed close to the main-shock nucleation. No substantial differences are found in the b-value computed before and after the main-shock struck in the small region of the fault plane populated by foreshocks

Insights into post-emplacement lava flow dynamics at Mt. Etna volcano from 2016 to 2021 by synthetic aperture radar and multispectral satellite data

Post-emplacement dynamics of lava flows is governed by several factors such as poroelastic deformation of the substrate; gravity-induced repacking and rearrangement of the vesicle-bearing fluid lava and other void spaces by superposed flows; lava densification processes; viscoelastic strain relaxation of the ground caused by the lava load; thermal cooling and contraction of the solid lava; and discrete motion of surface blocks. Here we investigate post-emplacement lava flow dynamics at the Mt. Etna volcano, and we infer on the possible causes by exploiting optical and radar satellite data. Synthetic aperture radar data from Sentinel-1 satellite mission provided high-resolution horizontal and vertical displacement rates and displacement time series of the lava flows emplaced on the Mt. Etna volcano summit from January 2016 to July 2021. Sentinel-2 multispectral data allowed to identify the lava flows boundaries emplaced during the December 2018 and May 2019 paroxysms. Finally, high resolution COSMO-SkyMed radar data allowed to account for the topographic changes generated by the lava emplacement by means of stereo radargrammetry technique. Such an unprecedented dataset provided a full picture of the lava flow dynamics, whose kinematics is governed lava cooling, which in turn produce thermal contraction of the lava body and viscous compaction of the underlying substrate. Both phenomena act at different periods, being the thermal contraction predominant for recent lava flows. Downslope sliding is also invoked, especially for recent lava flows emplaced on high slope areas

Rilocalizzazione di terremoti, modelli 3D di Vp, Vp/Vs e Qp nell’area geotermica di Larderello.

L'area geotermica di Larderello è caratterizzata da un elevato flusso di calore (tra 120 e 1000 mW/m2) e geologicamente caratterizzata da bacini postorogenici riempiti da depositi neogenici, risultato di un diffuso assottigliamento crostale accompagnato da risalite di materiale dal mantello superiore fino a profondità di pochi chilometri

Shallow subduction beneath Italy

This paper presents a velocity model of the Italian (central Mediterranean) lithosphere in unprecedented detail. The model is derived by inverting a set of 166,000 Pg and Pn seismic wave arrival times, restricted to the highest-quality data available. The tomographic images reveal the geometry of the subduction-collision system between the European, Adriatic, and Tyrrhenian plates, over a larger volume and with finer resolution than previous studies. We find two arcs of low-Vp anomalies running along the Alps and the Apennines, describing the collision zones of underthrusting continental lithospheres. Our results suggest that in the Apennines, a significant portion of the crust has been subducted below the mountain belt. From the velocity model we can also infer thermal softening of the crustal wedge above the subducting Adriatic plate. In the Tyrrhenian back-arc region, strong and extensive low-Vp anomalies depict upwelling asthenospheric material. The tomographic images also allow us to trace the boundary between the Adriatic and the Tyrrhenian plates at Moho depth, revealing some tears in the Adriatic-Ionian subducting lithosphere. The complex lithospheric structure described by this study is the result of a long evolution; the heterogeneities of continental margins, lithospheric underthrusting, and plate indentation have led to subduction variations, slab tears, and asthenospheric upwelling at the present day. The high-resolution model provided here greatly improves our understanding of the central Mediterranean’s structural puzzle. The results of this study can also shed light on the evolution of other regions experiencing both oceanic and continental subduction

Automatic seismic phase picking and consistent observation error assessment: application to the Italian seismicity

Accuracy of seismic phase observation and consistency of timing error assessment define the quality of seismic waves arrival times. High-quality and large data sets are prerequisites for seismic tomography to enhance the resolution of crustal and upper mantle structures. In this paper we present the application of an automated picking system to some 600 000 seismograms of local earthquakes routinely recorded and archived by the Italian national seismic network. The system defines an observation weighting scheme calibrated with a hand-picked data subset and mimics the picking by an expert seismologist. The strength of this automatic picking is that once it is tuned for observation quality assessment, consistency of arrival times is strongly improved and errors are independent of the amount of data to be picked. The application to the Italian local seismicity documents that it is possible to automatically compile a precise, homogeneous and large data set of local earthquake Pg and Pn arrivals with related polarities. We demonstrate that such a data set is suitable for high-precision earthquake location, focal mechanism determination and high-resolution seismic tomograph

Flank sliding: A valve and a sentinel for paroxysmal eruptions and magma ascent at Mount Etna, Italy

Rising magma, dike intrusions, and flank collapse are observed at many volcanoes worldwide, but how they interact is still poorly documented. Extensive synthetic aperture radar interferometry and continuous global positioning system observations captured a sharp dike intrusion at Mount Etna, Italy, during the 2018 paroxysm that triggered a vigorous seaward sliding of the eastern flank connected with brittle failure and deep magmatic resourcing. We propose a feedback process between flank acceleration and magma intrusion that derives from the interaction between the long- and short-term deformation of the volcano. The flank sliding acts as a valve that modulates the emplacement and eruption of magma within the shallow system. Rapid flank acceleration could potentially evolve into sudden collapses and seismic release at shallow depth. In turn, flank slip events could act as a sentinel for changes in magma depth and paroxysmal eruptions at Mount Etna

From 3D to 4D passive seismic tomography: The sub-surface structure imaging of the Val d’Agri region, southern Italy

Local earthquakes (passive seismic) tomography (LET) is a well established tool for the imaging of the sub-surface structure. Alternative to active seismics, the main advantages of using natural sources are the better sounding in deeper portions of the upper crust, the relatively low cost, and the direct availability of S-waves. The main drawback is the achievable model resolution, which is limited by the density of the seismic network and the distribution of elastic sources, rather than the elastic wave frequency. Recently, 4D variations (in space and time) of velocity anomalies have been recognized in active volcanoes (Patanè et al., 2006) and normal faulting systems and ascribed to the medium response to transient geological processes, like dyke intrusions or fluid pressure increase on fault planes. In this paper we show how LET contributes to the imaging of the upper crust in a very attractive region like the Val d’Agri in southern Italy, which hosts both significant oil fields and seismogenic structures. We show that LET allows to improve the definition of the crust structure, at depths larger than those sampled by conventional seismic profiles, and detect the space-time dependency of elastic properties in response to local variations of fluid pressur

Discovering geothermal supercritical fluids: a new frontier for seismic exploration

Exploiting supercritical geothermal resources represents a frontier for the next generation of geothermal electrical power plant, as the heat capacity of supercritical fluids (SCF),which directly impacts on energy production, is much higher than that of fluids at subcritical conditions. Reconnaissance and location of intensively permeable and productive horizons at depth is the present limit for the development of SCF geothermal plants. We use, for the first time, teleseismic converted waves (i.e. receiver function) for discovering those horizons in the crust. Thanks to the capability of receiver function to map buried anisotropic materials, the SCF-bearing horizon is seen as the 4km-depth abrupt termination of a shallow, thick, ultra-high (>30%) anisotropic rock volume, in the center of the Larderello geothermal field. The SCF-bearing horizon develops within the granites of the geothermal field, bounding at depth the vapor-filled heavily-fractured rock matrix that hosts the shallow steam-dominated geothermal reservoirs. The sharp termination at depth of the anisotropic behavior of granites, coinciding with a 2 km-thick stripe of seismicity and diffuse fracturing, points out the sudden change in compressibility of the fluid filling the fractures and is a key-evidence of deep fluids that locally traversed the supercritical conditions. The presence of SCF and fracture permeability in nominally ductile granitic rocks open new scenarios for the understanding of magmatic systems and for geothermal exploitation