Seismic anisotropy reveals focused mantle flow around the Calabrian slab (Southern Italy)

SKS splitting at the Calabrian subduction zone, with delay times (δt) up to 3s, reveals the presence of a strong anisotropic fabric. Fast directions (ϕ) are oriented NNE-SSW in the Calabrian Arc (C.A.) and rotate NNW-SSE to the north following the arcuate shape of the subducting plate. We interpret the trench-parallel ϕ as local-scale mantle flow driven by the retrograde motion of the slab; the absence of trench perpendicular ϕ below the Southern Apennines (S.A.) excludes the presence of a seismically detectable return flow at its NE edge. This may be due to the relative youth and limited width of the S.A. slab tear. A possible return flow is identified farther north at the boundary of the S.A. and Central Apennines. Different and weaker anisotropy is present below the Apulian Platform (A.P.). This implies that the influence of the slab rollback in the sub-slab mantle is limited to less then 100 km from the slab

SKS splitting in Southern Italy: anisotropy variations in a fragmented subduction zone.

In this paper we present a collection of good quality shear wave splitting measurements in Southern Italy. In addition to a large amount of previous splitting measurements, we present new data from 15 teleseisms recorded from 2003 to 2006 at the 40 stations of the CAT/SCAN temporary network. These new measurements provide additional constraints on the anisotropic behaviour of the study region and better define the fast directions in the southern part of the Apulian Platform. For our analysis we have selected wellrecorded SKS phases and we have used the method of Silver and Chan to obtain the splitting parameters: the azimuth of the fast polarized shear wave (φ) and delay time (δt). Shear wave splitting results reveal the presence of a strong seismic anisotropy in the subduction system below the region. Three different geological and geodynamic regions are characterized by different anisotropic parameters. The Calabrian Arc domain has fast directions oriented NNE–SSW and the Southern Apennines domain has fast directions oriented NNW–SSE. This rotation of fast axes, following the arcuate shape of the slab, is marked by a lack of resolved measurements which occurs at the transition zone between those two domains. The third domain is identified in the Apulian Platform: here fast directions are oriented almost N–S in the northern part and NNE–SSW to ENE–WSW in the southern one. The large number of splitting parameters evaluated for events coming from different back-azimuth allows us to hypothesize the presence of a depth-dependent anisotropic structure which should be more complicated than a simple 2 layer model below the Southern Apennines and the Calabrian Arc domains and to constrain at 50 km depth the upper limit of the anisotropic layer, at least at the edge of Southern Apennines and Apulian Platform. We interpret the variability in fast directions as related to the fragmented subduction system in the mantle of this region. The trench-parallel φ observed in Calabrian Arc and in Southern Apennines has its main source in the asthenospheric flow below the slab likely due to the pressure induced by the retrograde motion of the slab itself. The pattern of φ in the Apulian Platform does not appear to be the direct result of the rollback motion of the slab, whose influence is limited to about 100 km from the slab. The anisotropy in the Apulian Platform may be related to an asthenospheric flow deflected by the complicated structure of the Adriatic microplate or may also be explained as frozen-in lithospheric anisotropy

Seismic Anisotropy beneath Northern Victoria Land from SKS Splitting Analysis

Abstract. Teleseismic data recorded by temporary and permanent stations located in the Northern Victoria Land region are analysed in order to identify the presence and location of seismic anisotropy. We work on data recorded by 24 temporary seismographic stations deployed between 1993 and 2000 in different zones of the Northern Victoria Land, and by the permanent very broad-band stations TNV located near the Italian Base M. Zucchelli. The temporary networks monitored an area extending from Terra Nova Bay towards the South beyond the David Glacier and up to the Indian Ocean northward. To better constrain our study, we also provide an analysis of data recorded by TNV in the same period of time and we take into account also SKS shear wave splitting measurements performed by Barruol and Hoffman (1999) on data recorded by DRV. This study, to be considered as preliminary, reveals the presence of seismic anisotropy below the study region, with a mainly NW-SE fast velocity direction below the Terra Nova Bay area and rather large delay times, that mean a deep rooted anisotropic layer

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