Anisotropic seismic structure of the lithosphere beneath the Adriatic coast of Italy constrained with mode‐converted body waves

PS converted waves observed near Ancona on the Adriatic coast of central Italy, as revealed by teleseismic receiver functions (RFs), vary with earthquake back‐azimuth and epicentral distance in a manner consistent with a 1‐D anisotropic seismic structure. Using reflectivity calculations, we develop a profile of anisotropic seismic velocity through the Adriatic lithosphere at this locality. We infer crustal thickness of ∼45 km. Anisotropy within the crust appears at ∼15‐km depth, suggesting a decollement between the subducting Adriatic lithosphere and the overriding crustal wedge. Lineation of inferred rock fabric is compatible with simple shear in ENE‐WSW direction. In the upper mantle, we infer an anisotropic layer at 80–90 km depth. If caused by olivine crystals alignment, the nearly north‐south lineation of the inferred rock fabric would be consistent with some nearby shear‐wave splitting observations. This anisotropic layer may be related to mantle deformation induced by the rollback of Adriatic lithosphere

S wave Splitting in Central Apennines (Italy): anisotropic parameters in the crust during seismic sequences

In this work, we reviewed the main anisotropic results obtained in the last two decades along the Central Apennines. Moreover, we improved this database, with new results coming from the seismicity that occurred in the Montereale area, between 2009 and 2017, which corresponds to a spatio-temporal gap in the previously analyzed datasets. The examined papers concerned both seismic sequences (as Colfiorito in 1997, Pietralunga in 2010, L’Aquila in 2009, Amatrice in 2016) and background seismicity (as the 2000-2001 Città di Castello experiment). The whole of the collected results shows a general NW-SE fast shear wave direction consistent with both the orientation of the extensional active Quaternary and inherited compressive fault systems, focal mechanisms and local stress field. Also, we observed a more intense anisotropy strength (normalized delay time > 0.006 s/km) nearby the strongest events (M > 5), all concentrated in the hanging-wall of the activated fault systems. In fact, this area is deeply affected by the surrounding rock volume perturbations that, in turn, have altered both the local stress field and crustal fracturing network. The most common anisotropic interpretative models that could explain our results are 1) the stress-induced anisotropy according to the Extensive-Dilatancy Anisotropy (EDA) model where the anisotropic pattern is related to the local stress variation and most of the variability is visible in time; 2) the tectonic-controlled anisotropy according to the Structural-Induced Anisotropy (SIA) model where the anisotropic pattern is related to the major structural features and most of the variability is visible only in space. As reported by the examined studies in Central Apennines the possibility to discriminate between stress and structural anisotropy is quite complex in a region where the directions of the extensional regime, the in situ horizontal maximum stress, the strike of major faults, both active and inherited coincide. Generally, in this review, we noted an overlap and mixture of the two aforementioned mechanisms and, just through a temporal analysis, made in the Montereale area, we supposed a predominant stressinduced anisotropy only in rock volumes where anisotropic parameter variations have been detected

The Pollino seismic sequence: Can shear wave anisotropy monitoring help earthquakes forecast?

Since the late the late ’60s-early ’70s era seismologists started developed theories that included variations of the elastic property of the Earth crust and the state of stress and its evolution crust prior to the oc- currence of a large earthquake. Among the others the theory of the dilatancy (Scholz et al., 1973): when a rock is subject to stress, the rock grains are shifted generating micro-cracks, thus the rock itself in- creases its volume. Inside the fractured rock, fluid saturation and pore pressure play an important role in earthquake nucleation, by modulating the effective stress. Thus measuring the variations of wave speed and of anisotropic parameter in time can be highly informative on how the stress leading to a major fault failure builds up. In 80s and 90s such kind of research on earthquake precursor slowed down and the priority was given to seismic hazard and ground motions studies, which are very important since these are the basis for the building codes in many countries. Today we have dense and sophisticated seismic networks to measure wave-fields characteristics: we archive continuous waveform data recorded at three components broad-band seismometers, we almost routinely obtain high resolution ear- thquake locations. Therefore we are ready to start to systematically look at seismic-wave propagation properties to possibly reveal short-term variations in the elastic properties of the Earth crust. In active fault areas and volcanoes, tectonic stress variation influences fracture field orientation and fluid migration processes, whose evolution with time can be monitored through the measurement of the anisotropic pa- rameters ( Piccinini et al., 2006). Through the study of S waves anisotropy it is therefore potentially possible to measure the presence, migration and state of the fluid in the rock traveled by seismic waves, thus providing a valuable route to understanding the seismogenic phenomena and their precursors (Crampin & Gao, 2010). Variations of anisotropic parameter and of the ratio between the compressional (P-wave) and the shear (S-wave) seismic velocities, the Vp/Vs (Nur, 1972) have been recently observed and measured during the preparatory phase of a major earthquake (Lucente et al. 2010). Here we show the anisotropic parameters at station MMN during the Pollino seismic sequence 2010-2013

Peeking inside the mantle structure beneath the Italian region through SKS shear wave splitting anisotropy: a review

Over the years, seismic anisotropy characterization has become one of the most popular methods to study and understand the Earth’s deep structures. Starting from more than 20 years ago, considerable progress has been made to map the anisotropic structure beneath Italy and the Central Mediterranean area. In particular, several past and current international projects (such as RETREAT, CAT/SCAN, CIFALPS, CIFALPS-2, AlpArray) focused on retrieving the anisotropic structure beneath Italy and surrounding regions, promoting advances in the knowledge of geological and geodynamical setting of this intriguing area. All of these studies aimed at a better understanding the complex and active geodynamic evolution of both the active and remnant subduction systems characterising this region and the associated Apennines, Alps and Dinaric belts, together with the Adriatic and Tyrrhenian basins. The presence of dense high-quality seismic networks, permanently run by INGV and other institutions, and temporary seismic stations deployed in the framework of international projects, the improvements in data processing and the use of several and even more sophisticated methods proposed to quantify the anisotropy, allowed to collect a huge amount of anisotropic parameters. Here a collection of all measurements done on core refracted phases are shown and used as a measure of mantle deformation and interpreted into geodynamic models. Images of anisotropy identify well-developed mantle flows around the sinking European and Adriatic slabs, recognised by tomographic studies. Slab retreat and related mantle flow are interpreted as the main driving mechanism of the Central Mediterranean geodynamics

CRUSTAL FRACTURING FIELD AS REVEALED BY SEISMIC ANISOTROPY IN THREE SEISMOGENIC AREAS OF THE APENNINIC CHAIN

In the last three years, we developed, tested and improved an automatic analysis code to calculate the shear wave splitting parameters, fast polarization direction (φ) and delay time (∂t). The code is a set of MatLab scripts able to retrieve crustal anisotropy parameters from three-component seismic recording of local earthquakes using horizontal component cross-correlation method. The analysis procedure consists in choosing an appropriate frequency range, that better highlights the signal containing the shear waves, and a length of time window on the seismogram centred on the S arrival (the temporal window contains at least one cycle of S wave). The code was compared to other two automatic analysis code (SPY and SHEBA) and tested on three Italian areas (Val d’Agri, Tiber Valley and L’Aquila surrounding) along the Apennine mountains. For each region we used the anisotropic parameters resulting from the automatic computation as a tool to determine the fracture field geometries connected with the active stress field. The anisotropic fast directions are used to define the active stress field (EDA model), finding a general consistence between fast direction and main stress indicators (focal mechanism and borehole break-out). The magnitude of delay time is used to define the fracture field intensity finding higher value in the volume where micro-seismicity occurs. Furthermore we studied temporal variations of anisotropic parameters in order to explain if fluids play an important role in the earthquake generation process. The close association of anisotropic parameters variations and seismicity rate changes supports the hypothesis that the background seismicity is influenced by the fluctuation of pore fluid pressure in the rocks

Seismic measurements to reveal short-term variations in the elastic properties of the Earth crust

Since the late the late ’60s-early ’70s era seismologists started developed theories that included variations of the elastic property of the Earth crust and the state of stress and its evolution crust prior to the occurrence of a large earthquake. Among the others the theory of the dilatancy (Scholz et al., 1973): when a rock is subject to stress, the rock grains are shifted generating micro-cracks, thus the rock itself increases its volume. Inside the fractured rock, fluid saturation and pore pressure play an important role in earthquake nucleation, by modulating the effective stress. Thus measuring the variations of wave speed and of anisotropic parameter in time can be highly informative on how the stress leading to a major fault failure builds up. In 80s and 90s such kind of research on earthquake precursor slowed down and the priority was given to seismic hazard and ground motions studies, which are very important since these are the basis for the building codes in many countries. Today we have dense and sophisticated seismic networks to measure wave-fields characteristics: we archive continuous waveform data recorded at three components broad-band seismometers, we almost routinely obtain high resolution earthquake locations. Therefore we are ready to start to systematically look at seismic-wave propagation properties to possibly reveal short-term variations in the elastic properties of the Earth crust. One seismological quantity which, since the ‘70s, is recognized to be diagnostic of the level of fracturation and/or of the pore pressure in the rock, hence of its state of stress, is the ratio between the compressional (P-wave) and the shear (S-wave) seismic velocities, the Vp/Vs (Nur, 1972; Kisslinger and Engdahl, 1973). Variations of this ratio have been recently observed and measured during the preparatory phase of a major earthquake (Lucente et al. 2010). In active fault areas and volcanoes, tectonic stress variation influences fracture field orientation and fluid migration processes, whose evolution with time can be monitored through the measurement of the anisotropic pa- rameters (Miller and Savage, 2001; Piccinini et al., 2006). Through the study of S waves anisotropy it is therefore potentially possible to measure the presence, migration and state of the fluid in the rock traveled by seismic waves, thus providing a valuable route to understanding the seismogenic phenomena and their precursors (Crampin & Gao, 2010). In terms of determination of Earth crust elastic properties, recent studies (Brenguier et al., 2008; Chen et al., 2010; Zaccarelli et al., 2011) have shown how it is possible to estimate the relative variations in the wave speed through the analysis of the crosscorrelation of ambient seismic noise. In this paper we analyze in detail two seismological methods dealing with shear wave splitting and seismic noise cross correlation: a short historical review, their theoretical bases, the problems, learnings, limitations and perspec- tives. Moreover we discuss the results of these methods already applied on the data recorded in the L’Aquila region, before and after the destructive earthquake of April 6th 2009, represent their self an interesting case study

SEISMIC ANISOTROPY AND ITS RELATION WITH FAULTS AND STRESS FIELD IN THEVAL D'AGRI (SOUTHERN ITALY).

Shear-wave splitting is measured at 17 seismic stations deployed in the Val DAgri by INGV, which recorded local back-ground seismicity from May 2005 to June 2006 . The splitting results suggest the presence of an anisotropic upper crust (max hypocentral depth 15.5 km). The dominant fast polarisation direction strikes NW-SE parallel to the Apennines orogen and is approximately parallel to the maximum horizontal stress in the region and also parallel to the strike of the main normal faults in the Val DAgri. The size of the delay times, average is 0.1 second suggests 4.5% shear-wave velocity anisotropy. At stations located at the North West portion of the deployment average delay times are larger on the order of 0.2s. These parameters agree with an interpretation of seismic anisotropy in terms of the Extensive-Dilatancy Anisotropy model which considers the rock volume to be pervaded by fluid-saturated microcracks aligned by the active stress field. We cannot completely rule out the contribution of aligned macroscopic fractures as the cause of the shear wave anisotropy even if the parallel shear-wave polarisations we found are diagnostic of transverse isotropy with a horizontal axis of symmetry. This symmetry is commonly explained by parallel stress-aligned microcracks

IMAGING THE ACTIVE STRESS FIELD OF THREE SEISMOGENIC AREAS ALONG THE APENNINES AS REVEALED BY CRUSTAL ANISOTROPY

During the last decades, the study of seismic anisotropy has provided useful information for the interpretation and evaluation of the stress field and active crustal deformation. Seismic anisotropy can yield valuable information on upper crustal structure, fracture field, and presence of fluid-saturated rocks crossed by shear waves. Several studies worldwide demonstrate that seismic anisotropy is related to stress-aligned, filled-fluid micro-cracks (EDA model). An automatic analysis code, “Anisomat+”, was developed, tested and improved to calculate the anisotropic parameters: fast polarization direction (φ) and delay time (∂t). Anisomat+ has been compared to other two automatic analysis codes (SPY and SHEBA) and tested on three zones of the Apennines (Val d’Agri, Tiber Valley and L’Aquila surroundings). The anisotropic parameters, resulting from the automatic computation, have been interpreted to determine the fracture field geometries; for each area, we defined the dominant fast direction and the intensity of the anisotropy, interpreting these results in the light of the geological and structural setting and of two anisotropic interpretative models, proposed in the literature. In the first one, proposed by Zinke and Zoback, the local stress field and cracks are aligned by tectonics phases and are not necessarily related to the presently active stress field. Therefore the anisotropic parameters variations are only space-dependent. In the second, EDA model, and its development in the APE model fluid-filled micro-cracks are aligned or ‘opened’ by the active stress field and the variation of the stress field might be related to the evolution of the pore pressure in time; therefore in this case the variation of the anisotropic parameters are both space- and time- dependent. We recognized that the average of fast directions, in the three selected areas, are oriented NW-SE, in agreement with the orientation of the active stress field, as suggested by the EDA model, but also, by the proposed by Zinke and Zoback model; in fact, NW-SE direction corresponds also to the strike of the main fault structures in the three study regions. The mean values of the magnitude of the normalized delay time range from 0.005 s/km to 0.007 s/km and to 0.009 s/km, respectively for the L'Aquila (AQU) area, the High Tiber Valley (ATF) and the Val d'Agri (VA), suggesting a 3-4% of crustal anisotropy. In each area are also examined the spatial and temporal distribution of anisotropic parameters, which lead to some innovative observations, listed below. 1) The higher values of normalized delay times have been observed in those zones where most of the seismic events occur. This aspect was further investigated, by evaluating the average seismic rate, in a time period, between years 2005 and 2010, longer than the lapse of time, analyzed in the anisotropic studies. This comparison has highlighted that the value of the normalised delay time is larger where the seismicity rate is higher. 2) In the Alto Tiberina Fault area the higher values of normalised delay time are not only related to the presence of a high seismicity rate but also to the presence of a tectonically doubled carbonate succession. Therefore, also the lithology, plays a important role in hosting and preserving the micro-fracture network responsible for the anisotropic field. 3) The observed temporal variations of anisotropic parameters, have been observed and related to the fluctuation of pore fluid pressure at depth possibly induced by different mechanisms in the different regions, for instance, changes in the water table level in Val D’Agri, occurrence of the April 6th Mw=6.1 earthquake in L’Aquila.Since these variations have been recognized, it is possible to affirm that the models that better fit the results, both in term of fast directions and of delay times, seems to be EDA and APE models

Reply to Comment by A. Argnani on "Geometry of the Deep Calabrian Subduction from Wide‐Angle Seismic Data and 3‐D Gravity Modeling"

Keypoints This contribution is a reply on a comment submitted by A. Argnani. The alternate interpretation of the wide-angle seismic model is discussed. The Alfeo Fault system is proposed to be the current location of STEP fault. Abstract Andrea Argnani in his comment on Dellong et al., 2020 (Geometry of the deep Calabrian subduction (Central Mediterranean Sea) from wide‐angle seismic data and 3D gravity modeling), proposes an alternate interpretation of the wide-angle seismic velocity models presented by Dellong et al., 2018 and Dellong et al., 2020 and proposes a correction of the literature citations in these paper. In this reply, we discuss in detail all points raised by Andrea Argnani

Passive seismology and deep structure in central Italy

n the last decade temporary teleseismic transects have become a powerful tool for investigating the crustal and upper mantle structure. In order to gain a clearer picture of the lithosphere-asthenosphere structure in peninsular Italy, between 1994 and 1996, we have deployed three teleseismic transects in northern, central, and southern Apennines, in the framework of the project GeoModAp (European Community contract EV5V-CT94–0464). Some hundreds of teleseisms were recorded at each deployment which lasted between 3 and 4 months. Although many analyses are still in progress, the availability of this high quality data allowed us to refine tomographic images of the lithosphere-asthenosphere structure with an improved resolution in the northern and central Apennines, and to study the deformation of the upper mantle looking at seismic anisotropy through shear-wave splitting analysis. Also, a study of the depth and geometry of the Moho through the receiver function technique is in progress. Tomographic results from the northernmost 1994 and the central 1995 teleseismic experiments confirm that a high-velocity anomaly (HVA) does exist in the upper 200–250 km and is confined to the northern Apenninic arc. This HVA, already interpreted as a fragment of subducted lithosphere is better defined by the new temporary data, compared to previous works, based only on data from permanent stations. No clear high-velocity anomalies are detected in the upper 250 km below the central Apennines, suggesting either a slab window due to a detachment below southern peninsular Italy, or a thinner, perhaps continental slab of Adriatic lithosphere not detectable by standard tomography. We found clear evidence of seismic anisotropy in the uppermost mantle, related to the main tectonic processes which affected the studied regions, either NE–SW compressional deformation of the lithosphere beneath the mountain belt, or arc-parallel asthenospheric flow (both giving NW–SE fast polarization direction), and successive extensional deformation ( E–W trending) in the back-arc basin of northern Tyrrhenian and Tuscany. Preliminary results of receiver function studies in the northern Apennines show that the Moho depth is well defined in the Tyrrhenian and Adriatic regions while its geometry underneath the mountain belt is not yet well constrained, due to the observed high complexity.Published479-4934T. Sismicità dell'ItaliaJCR Journa