Array Analyses of Low-Frequency (0.1-0.5 Hz) Ambient Noise in Central Italy

In the framework of the seismological studies related to the activity of the Alto Tiberina Fault (ATF), a seismic array composed by 9 stations was deployed in the vicinity of Gubbio, central Italy

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

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

Surgical treatment of gastrointestinal stromal tumors of the duodenum. A literature review

Background: Gastrointestinal stromal tumors (GIST) are the most frequent mesenchymal tumours in the digestive tract. The duodenal GIST (dGIST) is the rarest subtype, representing only 4–5% of all GIST, but up to 21% of the resected ones. The diagnostic and therapeutic management of dGIST may be difficult due to the rarity of this tumor, its anatomical location, and the clinical behavior that often mimic a variety of conditions; moreover, there is lack of consent for their treatment. This study has evaluated the scientific literature to provide consensus on the diagnosis of dGIST and to outline possible options for surgical treatment. Methods: An extensive research has been carried out on the electronic databases MEDLINE, Scopus, EMBASE and Cochrane to identify all clinical trials that report an event or case series of dGIST. Results: Eighty-six studies that met the inclusion criteria were identified with five hundred forty-nine patients with dGIST: twenty-seven patients were treated with pancreatoduodenectomy and ninety-six with only local resection (segmental/wedge resections); in four hundred twenty-six patients it is not possible identify the type of treatment performed (pancreatoduodenectomy or segmental/wedge resections). Conclusions: dGISTs are a very rare subset of GISTs. They may be asymptomatic or may involve symptoms of upper GI bleeding and abdominal pain at presentation. Because of the misleading clinical presentation the differential diagnosis may be difficult. Tumours smaller than 2 cm have a low biological aggressiveness and can be followed annually by endoscopic ultrasound. The biggest ones should undergo radical surgical resection (R0). In dGIST there is no uniformly adopted surgical strategy because of the low incidence, lack of experience, and the complex anatomy of the duodenum. Therefore, individually tailored surgical approach is recommended. R0 resection with 1–2 cm clear margin is required. Lymph node dissection is not recommended due to the low incidence of lymphatic metastases. Tumor rupture should be avoided

The complex architecture of the 2009 MW 6.1 L’Aquila normal fault system (Central Italy) as imaged by 64,000 high-resolution aftershock locations.

On April 6th 2009, a MW 6.1 normal faulting earthquake struck the axial area of the Abruzzo region in Central Italy. We present high-precision hypocenter locations of an extraordinary dataset composed by 64,000 earthquakes recorded at a very dense seismic network of 60 stations operating for 9 months after the main event. Events span in magnitude (ML) between -0.9 to 5.9, reaching a completeness magnitude of 0.7. The dataset has been processed by integrating an accurate automatic picking procedure together with cross-correlation and double-difference relative location methods. The combined use of these procedures results in earthquake relative location uncertainties in the range of a few meters to tens of meters, comparable/lower than the spatial dimension of the earthquakes themselves). This data set allows us to image the complex inner geometry of individual faults from the kilometre to meter scale. The aftershock distribution illuminates the anatomy of the en-echelon fault system composed of two major faults. The mainshock breaks the entire upper crust from 10 km depth to the surface along a 14-km long normal fault. A second segment, located north of the normal fault and activated by two Mw>5 events, shows a striking listric geometry completely blind. We focus on the analysis of about 300 clusters of co-located events to characterize the mechanical behavior of the different portions of the fault system. The number of events in each cluster ranges from 4 to 24 events and they exhibit strongly correlated seismograms at common stations. They mostly occur where secondary structures join the main fault planes and along unfavorably oriented segments. Moreover, larger clusters nucleate on secondary faults located in the overlapping area between the two main segments, where the rate of earthquake production is very high with a long-lasting seismic decay

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

The deep structure of the Larderello-Travale geothermal field (Italy) from integrated, passive seismic investigations

AbstractWe report the preliminary results from a project (GAPSS-Geothermal Area Passive Seismic Sources), aimed at testing the resolving capabilities of passive exploration methods on a well-known geothermal area, namely the Larderello-Travale Geothermal Field (LTGF). Located in the western part of Tuscany (Italy), LTGF is the most ancient geothermal power field of the world. GAPSS consisted of up to 20 seismic stations deployed over an area of about 50 x 50 Km. During the first 12 months of measurements, we located more than 2000 earthquakes, with a peak rate of up to 40 shocks/day. Preliminary results from analysis of these signals include: (i) analysis of Shear-Wave-Splitting from local earthquake data, from which we determined the areal distribution of the most anisotropic regions; (ii) local-earthquake travel-time tomography for both P- and S-wave velocities; (iii) telesismic receiver function aimed at determining the high-resolution (<0.5km) S-velocity structure over the 0-20km depth range, and seismic anisotropy using the decomposition of the angular harmonics of the RF data-set; (iv) S-wave velocity profiling through inversion of the dispersive characteristics of Rayleigh waves from earthquakes recorded at regional distances. After presenting results from these different analyses, we eventually discuss their potential application to the characterisation and exploration of the investigated area

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