de Franco, Roberto

This paper presents an overview of the geophysical activities for the seismic microzonation of 138 municipalities belonging to four Italian regions (Abruzzo, Lazio, Marche and Umbria) that were severely damaged by the seismic sequence of Central Italy (August 2016–January 2017). This study is the result of a collaborative effort between research Institutions and professional geologists with the support of local Administrations and the Italian Civil Protection Department and sets an unprecedented large-scale example of geophysical investigations supporting detailed seismic microzonation studies. This manuscript presents the methodological approach adopted for the geophysical activities, including the technical protocols and procedures, the best practices, the final products and the results supporting a detailed microzonation study of III level. The first step of the study was the collection and critical review of all available geophysical and geological information for planning the new geophysical surveys (specifically their type and location), in order to assess the subsoil geometry and the seismic characterization of the areas under investigation. Integration with the newly acquired geophysical data allowed the identification of zones with homogeneous local seismic hazard as well as the reference seismo-stratigraphy for each area, defining for each geological unit the ranges of the relevant properties in seismic amplification studies: layering and thicknesses, density, P-wave and S-wave seismic velocity. We also present a few representative case studies illustrating the geophysical investigation for different geomorphological situations. These examples, together with the findings of the entire project, are discussed to point out the strength points and the criticalities, as well as the necessary requirements in the application of geophysical methods to detailed microzonation studies.

The formation of Cenozoic mountain belts in the Mediterranean realm was preceded by tens of millions of years of subduction, forming volcanic arcs, and frontal contractional systems. In addition, subduction usually involves slab rollback and formation of oceanic backarcs. Although such structure must have influenced the orogeny of Mediterranean mountain belts, no active analog has been mapped with modern crustal-scale seismic methods. Here, we study the entire Calabrian subduction system to map the structure resulting from Tethys lithosphere subduction and slab rollback, in a process that must be akin to that operating during a phase of the formation of the Mediterranean orogenic belts. We present a crustal-scale cross section of the entire Calabrian subduction system obtained from on- and off-shore wide-angle seismic data. The 2D P-wave velocity section shows spatially abrupt (<5 km of profile distance) structural and petrological transitions from the Ionian sedimentary wedge and Calabrian arc, to the rifted NW Calabrian margin, where the Quaternary Aeolian arc is emplaced. The margin, then, transitions northwards into the Marsili backarc region, where exhumed mantle and localized volcanism occurred during its formation. This complex structure implies rapid temporal and spatial changes between magmatic and amagmatic processes, and between compressional and extensional regimes during the evolution of this subduction system. We find that some terranes involved in the Alpine orogeny share petrological and tectonic similarities with some domains of the Calabrian subduction system. Based on the results of this study we propose the Calabrian Arc system as an analog for the subduction structuration that preceded the formation of Alpine orogenic systems.

The first mainshock (Mw 6.0) of the 2016 Central Italy seismic sequence, severely struck the Amatrice village and the surrounding localities. After a few days, some Italian Institutions, coordinated by the “Center for Seismic Microzonation and its applications”, carried out several preparatory activities for seismic microzonation of the area. A temporary seismic network was installed that monitored about 50 sites in epicentral area. The network produced a huge amount of records in a wide range of magnitude up to Mw 6.5. For about half of the recording stations, detailed site characterization was undertaken, encompassing single station noise measurements and S-wave velocity profiles. The geological and geophysical data together with the collected dataset of seismic signals were exploited to investigate the site response of selected stations. Significant amplifications are found in the correspondence of several sites that experienced a high level of damage (Imcs >IX), mainly at short and intermediate periods

In August 2016, a magnitude 6.0 earthquake struck Central Italy, starting a devastating seismic sequence, aggravated by other two events of magnitude 5.9 and 6.5, respectively. After the first mainshock, four Italian institutions installed a dense temporary network of 50 seismic stations in an area of 260 km2. The network was registered in the International Federation of Digital Seismograph Networks with the code 3A and quoted with a Digital Object Identifier ( https://doi.org/10.13127/SD/ku7Xm12Yy9 ). Raw data were converted into the standard binary miniSEED format, and organized in a structured archive. Then, data quality and completeness were checked, and all the relevant information was used for creating the metadata volumes. Finally, the 99 Gb of continuous seismic data and metadata were uploaded into the INGV node of the European Integrated Data Archive repository. Their use was regulated by a Memorandum of Understanding between the institutions. After an embargo period, the data are now available for many different seismological studies.

The three-dimensional pattern of elastic moduli (bulk modulus, Young modulus, shear modulus) of the upper crust (0-10 km depth) has been determined in the Friuli area (north-eastern Italy) from the 3D Vp, Vp/Vs and density structures. Firstly, 3D Pwave velocity and P to S velocity ratio were modeled by joint inversion for hypocentres and velocity structure. Then, we apply the tomographic inversion method of Sequential Integrated Inversion (SII) to recover the three dimensional density structure. The pattern of the elastic moduli is characterized by marked lateral and depth variations that reflect the geologic-structural heterogeneity of the area, produced by the superposition of several tectonic phases with different orientations of the principal axes of stress. The bulk (K), Young (E) and shear (G) moduli image a high rigidity body with an irregular shape, at 4-8 km depth. The body is characterized by G ≥ 3.2·1010 N·m-2, K ≥ 6.8·1010 N·m-2 and E ≥ 8.4·1010 N·m-2 and is associated to platform limestones and dolomitic rocks. The seismicity is mainly located along the sharp variations of the moduli pattern, in or adjacent to high rigidity zones. The most severe earthquakes (ML between 4.5 and 6.4), occurred in the study area from 1976 to the present day, are located in a transition zone from high to low rigidity patterns. Our interpretation is that the elastic moduli variations, closely related to variability in rock mechanical properties, influence the occurrence of earthquakes by processes of stress concentrations. The values of the elastic moduli recently obtained from laboratory measurements on the main lithologic units fall in the middlehigh range of the values obtained with the present investigation. Keywords: Seismic tomography, gravity anomalies, seismicity, elastic moduli, Friuli, NE Italy

The objective of this work is to perform a purely empirical assessment of the actual capabilities of the horizontal-to-vertical (H/V) spectral ratio technique to provide reliable and relevant information concerning site conditions and/or site amplification. This objective has been tackled through the homogeneous (re)processing of a large volume of earthquakes and ambient noise data recorded by different research teams in more than 200 sites located mainly in Europe, but also in the Caribbean and in Tehran. The original recordings were first gathered in a specific database with information on both the sites and recorded events. Then, for all sites close to an instrumented reference, average site-to-reference spectral ratios (“spectral ratio method” (SSR)) were derived in a homogeneous way (window selection, smoothing, signal-to-noise ratio threshold, averaging), as well as H/V ratios (“HVSRE–RF”) on earthquake recordings. H/V ratios were also obtained from noise recordings at each site (either specific measurements, or extracted from pre- or post-event noise windows). The spectral curves resulting from these three techniques were estimated reliable for a subset of 104 sites, and were thus compared in terms of fundamental frequency, amplitude and amplification bandwidth, exhibiting agreements and disagreements, for which interpretations are looked for in relation with characteristics of site conditions. The first important result consists in the very good agreement between fundamental frequencies obtained with either technique, observed for 81% of the analyzed sites. A significant part of the disagreements correspond to thick, low frequency, continental sites where natural noise level is often very low and H/V noise ratios do not exhibit any clear peak. The second important result is the absence of correlation between H/V peak amplitude and the actual site amplification measured on site-to-reference spectral ratios. There are, however, two statistically significant results about the amplitude of the H/V curve: the peak amplitude may be considered as a lower bound estimate of the actual amplification indicated by SSR (it is smaller for 79% of the 104 investigated sites), and, from another point of view, the difference in amplitude exhibits a questioning correlation with the geometrical characteristics of the sediment/basement interface: large SSR/HV differences might thus help to detect the existence of significant 2D or 3D effects.

A new experiment called SERAPIS (SEismic Reflection/Refraction Acquisition Project for Imaging complex volcanic Structures) has been planned and carried out, based on off-shore seismic energization and data acquisition on land and on sea-bottom. The experiment was performed in September, 2001 during which the vessel NADIR of IFREMER (equipped with 12, 16-liters airgun) produced more than 5000 air gun shots recorded at a sea-bottom seismograph array of 72 OBS and 62 stations installed on-land. Active seismic refraction DSS (Deep Seismic Soundings) acquired during the surveys conducted in 1980 and 1985 were recovered jointly with seismic data acquired in the Campi Flegrei area in the framework of the MareVes97 (an experiment devoted to the definition of the structure of the Somma-Vesuvio complex) offshore survey. The data set acquired during the SERAPIS experiment has been successfully used to infer 3D images of the volcanic structures of Campi Flegrei and Neapolitan bay. Active seismic waveforms and related P-picks (more than 90000 data) from the SERAPIS experiment are also available in the project data server.

We performed an analysis of refraction data recorded in Italy since 1968 in the frame of the numerous deep seismic sounding and wide-angle reflection/refraction projects. The aims of this study are to construct a parametric database including the recording geometric information relative to each profile, the phase pickings and the results of some kinematic analyses performed on the data, and to define a reference 1D velocity model for the Italian territory from all the available refraction data. As concerns the first goal, for each seismic section we picked the P-wave first-arrival-times, evaluated the uncertainties of the arrival-times pickings and determined from each travel time-offset curve the 1D velocity model. The study was performed on 419 seismic sections. Picking was carried out manually by an algorithm which includes the computation of three picking functions and the picking- error estimation. For each of the travel time-offset curves a 1D velocity model has been calculated. Actually, the 1D velocity-depth functions were estimated in three different ways which assume: a constant velocitygradient model, a varying velocity-gradient model and a layered model. As regards the second objective of this work, a mean 1D velocity model for the Italian crust was defined and compared with those used for earthquake hypocentre locations and seismic tomographic studies by different institutions operating in the Italian area, to assess the significance of the model obtained. This model can be used in future works as input for a next joint tomographic inversion of active and passive seismic data.

The effect of local site amplification has been recognized as an important factor in ground motion assessment and is nowadays frequently studied. A common approach to evaluate ground shaking is to first estimate ground motion parameters at rock sites and then to correct them introducing site transfer functions derived from experimental data and from numerical modelling. The site transfer functions to be used to modify ground motion evaluated at rock sites can be evaluated starting from strong motion, weak motion and microtremor data. According to the amount and quality of the available data the transfer functions are evaluated for specific sites, in order to be used as a punctual information, or as representative of an average local condition in selected areas. All the available geological and geotechnical data must be collected to put some constraints on the obtained results and to permit numerical modeling to be compared with experimental results. The obtained transfer functions can be introduced in scenario studies convolving rock seismograms by the pulse response of the upper layers for different situations considered as representative of the geology of the studied areas. The capability of describing local site effects is strongly affected by the amount of seismological, geophysical and geotechnical data available. This is particularly true if numerical modelling needs to be performed and if the contribution of non linear soil behaviour has to be taken into account. For the three areas investigated in the framework of the project, the different amount of available data and information guided the performed studies and the obtained results. For Colfiorito test site, the availability of strong motion data recorded during the largest events of the Umbria Marche sequence (1997-98) yields well constrained information for specific sites. For Città di Castello the collection of weak motion and microtremor data allowed to reconstruct the geometry of the sedimentary basin underlying the city and to define zones with homogenous site response where to evaluate site transfer functions in a 1D approximation, including non linear behaviour. For Val D’Agri area, the lack of seismic and geotechnical data did not allow to describe in detail the site response in the sedimentary basin. In this case some sample sites with a known uppermost geological structure were selected as representative of the seismic response of the basin. For them, microtremor data were collected to put some constraint on the transfer functions computed in a 1D approximation. In this case it was not possible to consider the non linear soil behaviour.

A site response experiment was performed in the basin of Città di Castello (a small town in Central Italy) in May 2001. This study is part of a project on the evaluation of seismic hazard in seismogenic areas funded by the Gruppo Nazionale Difesa dai Terremoti (GNDT). The experiment consisted of a dense fixed transect configuration with most of the stations recording in continuous mode, and several ambient noise measurements both in single station and in array configuration spread over the investigated area. The dense transect was composed of 26 seismic stations in a crosswise configuration with a maximum inter-station distance of 250 m. The stations were deployed in the southern part of the basin, from the eastern bedrock outcrop to the western edge, across the town. About 70 earthquakes were recorded during 10 days of deployment, generally low magnitude or regional events. We located 23 earthquakes and 17 of them were located using the waveform similarity approach at 4 stations outside the target area. These 4 stations were part of a dense temporary seismic network involved in a previous experiment of the same project, aimed at performing a high-resolution picture of the local seismicity. Delay analysis on the recorded waveforms allowed us to infer the basin geometry at depth and estimate the S-wave velocity of sediments. Moreover, we evaluated relative site response along the E-W transect by performing a standard spectral ratio. Amplification factors up to 9 are found inside the basin; at frequencies above 5 Hz stations closer to the edges show higher amplification, whereas stations located in the middle of the basin, where the alluvial sediments are thicker (CD11-CD14), show higher amplification below 5 Hz. We considered the average amplification in two frequency bands (1-5 Hz and 5-10 Hz), representative of the resonance frequency for 2-3 storey buildings and 1 storey houses,respectively. Our results suggest that the potential hazard for 2-3 storey buildings is higher in the center of the basin (amplification factor up to 6), and for 1 storey houses is higher at the edges (amplification factor up to 5).