Continental degassing of helium in an active tectonic setting (northern Italy): the role of seismicity

In order to investigate the variability of helium degassing in continental regions, its release from rocks and emission into the atmosphere, here we studied the degassing of volatiles in a seismically active region of northern Italy (MwMAX = 6) at the Nirano-Regnano mud volcanic system. The emitted gases in the study area are CH4–dominated and it is the carrier for helium (He) transfer through the crust. Carbon and He isotopes unequivocally indicate that crustal-derived fluids dominate these systems. An high-resolution 3-dimensional reconstruction of the gas reservoirs feeding the observed gas emissions at the surface permits to estimate the amount of He stored in the natural reservoirs. Our study demonstrated that the in-situ production of 4He in the crust and a long-lasting diffusion through the crust are not the main processes that rule the He degassing in the region. Furthermore, we demonstrated that micro-fracturation due to the field of stress that generates the local seismicity increases the release of He from the rocks and can sustain the excess of He in the natural reservoirs respect to the steady-state diffusive degassing. These results prove that (1) the transport of volatiles through the crust can be episodic as function of rock deformation and seismicity and (2) He can be used to highlight changes in the stress field and related earthquakes

The shallow boreholes at The AltotiBerina near fault Observatory (TABOO; northern Apennines of Italy)

Abstract. As part of an interdisciplinary research project, funded by the European Research Council and addressing the mechanics of weak faults, we drilled three 200–250 m-deep boreholes and installed an array of seismometers. The array augments TABOO (The AltotiBerina near fault ObservatOry), a scientific infrastructure managed by the Italian National Institute of Geophysics and Volcanology. The observatory, which consists of a geophysical network equipped with multi-sensor stations, is located in the northern Apennines (Italy) and monitors a large and active low-angle normal fault. The drilling operations started at the end of 2011 and were completed by July 2012. We instrumented the boreholes with three-component short-period (2 Hz) passive instruments at different depths. The seismometers are now fully operational and collecting waveforms characterised by a very high signal to noise ratio that is ideal for studying microearthquakes. The resulting increase in the detection capability of the seismic network will allow for a broader range of transients to be identified

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

Radiography of a normal fault system by 64,000 high-precision earthquake locations: The 2009 L'Aquila (central Italy) case study

We studied the anatomy of the fault system where the 2009 L'Aquila earthquake (M_W 6.1) nucleated by means of ~64 k high-precision earthquake locations spanning 1 year. Data were analyzed by combining an automatic picking procedure for P and S waves, together with cross-correlation and double-difference location methods reaching a completeness magnitude for the catalogue equal to 0.7 including 425 clusters of similar earthquakes. The fault system is composed by two major faults: the high-angle L'Aquila fault and the listric Campotosto fault, both located in the first 10 km of the upper crust. We detect an extraordinary degree of detail in the anatomy of the single fault segments resembling the degree of complexity observed by field geologists on fault outcrops. We observe multiple antithetic and synthetic fault segments tens of meters long in both the hanging wall and footwall along with bends and cross fault intersections along the main fault and fault splays. The width of the L'Aquila fault zone varies along strike from 0.3 km where the fault exhibits the simplest geometry and experienced peaks in the slip distribution, up to 1.5 km at the fault tips with an increase in the geometrical complexity. These characteristics, similar to damage zone properties of natural faults, underline the key role of aftershocks in fault growth and co-seismic rupture propagation processes. Additionally, we interpret the persistent nucleation of similar events at the seismicity cutoff depth as the presence of a rheological (i.e., creeping) discontinuity explaining how normal faults detach at depth

Radiography of a normal fault system by 64,000 high-precision earthquake locations: The 2009 L’Aquila (central Italy) case study

We studied the anatomy of the fault system where the 2009 L’Aquila earthquake (MW 6.1) nucleated by means of ~64 k high-precision earthquake locations spanning 1 year. Data were analyzed by combining an automatic picking procedure for P and S waves, together with cross-correlation and double-difference location methods reaching a completeness magnitude for the catalogue equal to 0.7 including 425 clusters of similar earthquakes. The fault system is composed by two major faults: the high-angle L’Aquila fault and the listric Campotosto fault, both located in the first 10 km of the upper crust. We detect an extraordinary degree of detail in the anatomy of the single fault segments resembling the degree of complexity observed by field geologists on fault outcrops. We observe multiple antithetic and synthetic fault segments tens of meters long in both the hanging wall and footwall along with bends and cross fault intersections along the main fault and fault splays. The width of the L’Aquila fault zone varies along strike from 0.3 km where the fault exhibits the simplest geometry and experienced peaks in the slip distribution, up to 1.5 km at the fault tips with an increase in the geometrical complexity. These characteristics, similar to damage zone properties of natural faults, underline the key role of aftershocks in fault growth and co-seismic rupture propagation processes. Additionally, we interpret the persistent nucleation of similar events at the seismicity cutoff depth as the presence of a rheological (i.e., creeping) discontinuity explaining how normal faults detach at depth

Analysis and interpretation of the impact of missense variants in cancer

Large scale genome sequencing allowed the identification of a massive number of genetic variations, whose impact on human health is still unknown. In this review we analyze, by an in silico-based strategy, the impact of missense variants on cancer-related genes, whose effect on protein stability and function was experimentally determined. We collected a set of 164 variants from 11 proteins to analyze the impact of missense mutations at structural and functional levels, and to assess the performance of state-of-the-art methods (FoldX and Meta-SNP) for predicting protein stability change and pathogenicity. The result of our analysis shows that a combination of experimental data on protein stability and in silico pathogenicity predictions allowed the identification of a subset of variants with a high probability of having a deleterious phenotypic effect, as confirmed by the significant enrichment of the subset in variants annotated in the COSMIC database as putative cancer-driving variants. Our analysis suggests that the integration of experimental and computational approaches may contribute to evaluate the risk for complex disorders and develop more effective treatment strategie

Aseismic deformation associated with an earthquake swarm in the northern Apennines (Italy)

Analyzing the displacement time series from continuous GPS (cGPS) with an Independent Component Analysis, we detect a transient deformation signal that correlates both in space and time with a seismic swarm activity (maximum M_w=3.69 ± 0.09) occurred in the hanging wall of the Altotiberina normal fault (Northern Apennines, Italy) in 2013–2014. The geodetic transient lasted ∼6 months and produced a NW-SE trending extension of ∼5.3 mm, consistent with the regional tectonic regime. The seismicity and the geodetic signal are consistent with slip on two splay faults in the Altotiberina fault (ATF) hanging wall. Comparing the seismic moment associated with the geodetic transient and the seismic events, we observe that seismicity accounts for only a fraction of the measured geodetic deformation. The combined seismic and aseismic slip decreased the Coulomb stress on the locked shallow portion of the ATF, while the transition region to the creeping section has been loaded

An automatically generated high-resolution earthquake catalogue for the 2016–2017 Central Italy seismic sequence, including P and S phase arrival times

The 2016–2017 central Italy earthquake sequence began with the first main shock near the town of Amatrice on August 24 (Mw 6.0), and was followed by two subsequent large events near Visso on October 26 (Mw 5.9) and Norcia on October 30 (Mw 6.5), plus a cluster of four events with Mw > 5.0 within few hours on 18 January 2017. The affected area had been monitored before the sequence started by the permanent Italian National Seismic Network (RSNC), and was enhanced during the sequence by temporary stations deployed by the National Institute of Geophysics and Volcanology and the British Geological Survey. By the middle of September, there was a dense network of 155 stations, with a mean separation in the epicentral area of 6–10 km, comparable to the most likely earthquake depth range in the region. This network configuration was kept stable for an entire year, producing 2.5 TB of continuous waveform recordings. Here we describe how this data was used to develop a large and comprehensive earthquake catalogue using the Complete Automatic Seismic Processor (CASP) procedure. This procedure detected more than 450 000 events in the year following the first main shock, and determined their phase arrival times through an advanced picker engine (RSNI-Picker2), producing a set of about 7 million P- and 10 million S-wave arrival times. These were then used to locate the events using a non-linear location (NLL) algorithm, a 1-D velocity model calibrated for the area, and station corrections and then to compute their local magnitudes (ML). The procedure was validated by comparison of the derived data for phase picks and earthquake parameters with a handpicked reference catalogue (hereinafter referred to as ‘RefCat’). The automated procedure takes less than 12 hr on an Intel Core-i7 workstation to analyse the primary waveform data and to detect and locate 3000 events on the most seismically active day of the sequence. This proves the concept that the CASP algorithm can provide effectively real-time data for input into daily operational earthquake forecasts, The results show that there have been significant improvements compared to RefCat obtained in the same period using manual phase picks. The number of detected and located events is higher (from 84 401 to 450 000), the magnitude of completeness is lower (from ML 1.4 to 0.6), and also the number of phase picks is greater with an average number of 72 picked arrival for a ML = 1.4 compared with 30 phases for RefCat using manual phase picking. These propagate into formal uncertainties of ±0.9 km in epicentral location and ±1.5 km in depth for the enhanced catalogue for the vast majority of the events. Together, these provide a significant improvement in the resolution of fine structures such as local planar structures and clusters, in particular the identification of shallow events occurring in parts of the crust previously thought to be inactive. The lower completeness magnitude provides a rich data set for development and testing of analysis techniques of seismic sequences evolution, including real-time, operational monitoring of b-value, time-dependent hazard evaluation and aftershock forecasting

Fault structure and slip localization in carbonate-bearing normal faults: An example from the Northern Apennines of Italy

Carbonate-bearing normal faults are important structures for controlling fluid flow and seismogenesis within the brittle upper crust. Numerous studies have tried to characterize fault zone structure and earthquake slip processes along carbonate-bearing faults. However, due to the different scales of investigation, these studies are not often integrated to provide a comprehensive fault image. Here we present a multi-scale investigation of a normal fault exhumed from seismogenic depths. The fault extends for a length of 10 km with a maximum width of about 1.5 km and consists of 5 sub-parallel and interacting segments. The maximum displacement (370e650 m) of each fault segment is partitioned along sub-parallel slipping zones extending for a total width of about 50 m. Each slipping zone is characterized by slipping surfaces exhibiting different slip plane phenomena. Fault rock development is controlled by the protolith lithology. In massive limestone, moving away from the slip surface, we observe a thin layer (<2 cm) of ultracataclasite, cataclasite (2e10 cm) and fault breccia. In marly limestone, the fault rock consists of a cataclasite with hydrofractures and smectite-rich pressure solution seams. At the micro-nanoscale, the slip surface consists of a continuous and thin (<300 mm) layer composed of coarse calcite grains (~5e20 mm in size) associated with sub-micrometer grains showing fading grain boundaries, voids and/or vesicles, and suggesting thermal decomposition processes. Micrometer-sized calcite crystals show nanoscale polysynthetic twinning affected by the occurrence of subgrain boundaries and polygonalized nanostructures. Investigations at the kilometres-tens of meter scale provide fault images that can be directly compared with high-resolution seismological data and when combined can be used to develop a comprehensive characterization of seismically active fault structures in carbonate lithologies. Micro and nanoscale investigations along the principal slipping zone suggest that different deformation processes, including plastic deformation and thermal decomposition, were active during seismic slip

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