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

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

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

Automatic seismic phase picking and consistent observation error assessment: application to 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 paperwe present the application of an automated picking system to some 600000 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 tomography

The 2009 L’Aquila (Central Italy) Seismic Sequence.

On April 6 (01:32 UTC) 2009 a MW 6.1 normal faulting earthquake struck the axial area of the Abruzzo region in Central Italy. The earthquake heavily damaged the city of L’Aquila and its surroundings, causing 308 casualties, 70,000 evacuees and incalculable losses to the cultural heritage. We present the geometry of the fault system composed by two main normal fault planes, reconstructed by means of seismicity distribution: almost 3000 events with ML≥1.9 occurred in the area during the 2009. The events have been located with a 1D velocity model we computed for the area by using data of the seismic sequence. The mainshock, located at around 9.3 km of depth beneath the town of L’Aquila, activated a 50° (+/- 3) SW-dipping and ~135° NW-trending normal fault with a length of about 16 km. The aftershocks activated the whole 10 km of the upper crust up to the surface. The geometry of the fault is coherent with the mapped San Demetrio-Paganica and Mt. Stabiata normal faults. The whole normal fault system that reached about 50 km of length by the end of December in the NW-trending direction, was activated within the first few days of the sequence when most of the energetic events occurred. The main shock fault plane was activated by a foreshock sequence culminated with a MW 4.0 on the 30th of March (13:38 UTC), showing extensional kinematic with a minor left lateral component. The second major structure, located to the north close to Campotosto village, is controlled by a MW 5.0 which occurred on the same day of the main shock (the 6th of April at 23:15 UTC) and by a MW 5.2 event (9th of April - 00:53 UTC). The fault plane shows a shallower dip angle with respect to the main fault plane, of about 35° with a tendency to flattening towards the deepest portion. Due to the lack of seismicity above 5 km depth, the connection between this structure and the mapped Monti della Laga fault is not straightforward. This northern segment is recognisable for about 12-14 km of length, always NW-trending and forming a right lateral step with the main fault plane. The result is a en-echelon system overlapping for about 6 km. Seismicity pattern also highlights the activation of numerous minor normal fault segments within the whole fault system. The deepest is located at around 13-15 km of depth, south of the L’Aquila mainshock, and it seems to be antithetic to the main fault plane