Upper Mantle Domains beneath Central-Southern Italy: Petrological, Geochemical and Geophysical Constraints

—The Italian peninsula shows high complexity of the mantle-crust system and of the Plio-Quaternary magmatism. The lithospheric thickness has remarkable lateral variations from about 110 km to about 30 km. Intermediate and deep-focus earthquakes indicate the presence of a lithospheric slab under the Aeolian-Calabrian area and at the southern end of Campania. Much less extensive intermediate-depth seismicity characterizes the Roman-Tuscany region, where the existence of a relic slab has been hypothesized. The deep seismicity in the southern Tyrrhenian Sea is associated with active calcalkaline to shoshonitic volcanism in the Aeolian arc. Alkaline potassic volcanism occurs in central Italy, and potassic lamproitic magmatism coexists with crustal anatectic and various types of hybrid rocks in the Tuscany area.¶The parallelism between changing magmatism and variation of the structure of the crust-mantle system makes central-southern Italy a key place where petrological and geophysical data can be used to work out an integrated model of the structure and composition of the upper mantle. Beneath Tuscany the upper mantle has been affected by intensive subduction-related metasomatism. This caused the formation of phlogopite-rich veins that cut through residual spinel-harzburgite and dunite. These veins, possibly partially molten, may explain the unusually soft mechanical properties that are detected just below the Moho. In the Roman Province, the upper mantle is formed by a relatively thin lid (the mantle part of the lithosphere) and by metasomatic fertile peridotite, probably connected with the upraise of an asthenospheric mantle wedge above the Apennines subduction zone. Geochemical data indicate that metasomatism, though still related to subduction, had different characteristics and age than in Tuscany. In the eastern sector of the Aeolian arc and in the Neapolitan area, the upper mantle appears to be distinct from the Roman and Tuscany areas and is probably formed by fertile peridotite contaminated by the presently active subduction of the Ionian Sea floor.¶The overall picture is that of a mosaic of various mantle domains that have undergone different evolutionary history in terms of both metasomatism and pre-metasomatic events. The coexistence side by side of these sectors is a key factor that has to be considered by models of the geodynamic evolution of the Central Mediterranean area.     

Neo-deterministic seismic hazard assessment in North Africa

North Africa is one of the most earthquake-prone areas of the Mediterranean. Many devastating earthquakes, some of them tsunami-triggering, inflicted heavy loss of life and considerable economic damage to the region. In order to mitigate the destructive impact of the earthquakes, the regional seismic hazard in North Africa is assessed using the neo-deterministic, multi-scenario methodology (NDSHA) based on the computation of synthetic seismograms, using the modal summation technique, at a regular grid of 0.2?×?0.2°. This is the first study aimed at producing NDSHA maps of North Africa including five countries: Morocco, Algeria, Tunisia, Libya, and Egypt. The key input data for the NDSHA algorithm are earthquake sources, seismotectonic zonation, and structural models. In the preparation of the input data, it has been really important to go beyond the national borders and to adopt a coherent strategy all over the area. Thanks to the collaborative efforts of the teams involved, it has been possible to properly merge the earthquake catalogues available for each country to define with homogeneous criteria the seismogenic zones, the characteristic focal mechanism associated with each of them, and the structural models used to model wave propagation from the sources to the sites. As a result, reliable seismic hazard maps are produced in terms of maximum displacement (D _max), maximum velocity (V _max), and design ground acceleration.     

1D velocity structure of the Po River plain (Northern Italy) assessed by combining strong motion and ambient noise data

Strong ground motions recorded on the sedimentary deposits of the Po River alluvial plain during the Emilia (Northern Italy) Mw 5.7 earthquake of May 29, 2012 are used to assess the vertical profile of shear-wave velocity above the limestone basement. Data were collected by a linear array installed for site effect studies after the Mw 5.9 mainshock of May 20, 2012. The array stations, equipped with both strong and weak motion sensors, are aligned in the South–North direction, at distances ranging from 1 to 26 km from the epicenter. The vertical components of ground motion show very distinctive, large-amplitude, low-frequency dispersive wave trains. Wavelet analysis yields group-velocity dispersion curve in the 0.2–0.7 Hz frequency band. The availability of a long ambient noise record allows estimates of the site resonance frequency along with its stability among stations. The joint inversion of dispersion of surface waves and ellipticity curves derived from ambient noise H/V allows extending investigations down to the sediment-limestone interface, at a depth of about 5,000 m. Our results add new information about the velocity structure at a scale that is intermediate between the local scale already investigated by other authors with small-aperture arrays using ambient noise and the regional scale inferred from modeling of seismogram waveforms recorded at hundreds of kilometers from the source.     

Forest structure of a maple old-growth stand: a case study on the Apennines mountains (Southern Italy)

In Europe, very small forest areas can be considered to be old-growth, and they are mainly located in Eastern Europe. The typical structures of old growth forests infrequently occur in Mediterranean mountainous environments, since they have been affected by human activities for centuries. This study focused on a remote and almost pure Italian maple stand located in southern Italy, which has not been managed for long time due to its inaccessibility. The effects of natural evolution on the forest stand were evaluated through the analysis of the spatial and chronological structure and the regeneration patterns, then estimating the amounts and quality of deadwood occurrence. Across the whole stand, all the trees with DBH (diameter at breast height) larger than 50 cm (LLT, large living trees) were measured (DBH and height) and age was also determined through a dendrochronological approach. The diameters observed ranged between 50 and 145 cm with ages of 120 to 250 years. The Latham index calculated for trees within the sample plot highlighted a multilayered canopy with a dominant layer of large living trees (age > 120 years). The size-class distribution of stems had a reverse-J shape, and basal area was 52 m² ha^-1. Deadwood was exclusively constituted by standing dead trees and CWD and its volume was on average 31 m³ ha^-1.Pure Italian maple forests are generally rare in Europe, and it was unexpected to find a forest stand characterized by a so complex structure with old growth attributes. The study of complex forest stand, even if small, could give precious information on the forest evolution, clarifying also diverse auto-ecological traits of tree species that usually are not common in our forests.     

Neo-Deterministic and Probabilistic Seismic Hazard Assessments: a Comparison over the Italian Territory

Estimates of seismic hazard obtained using the neo-deterministic approach (NDSHA) and the probabilistic approach (PSHA) are compared for the Italian territory. The NDSHA provides values larger than those given by the PSHA in areas where large earthquakes are observed and in areas identified as prone to large earthquakes, but lower values in low-seismicity areas. These differences suggest the adoption of the flexible, robust and physically sound NDSHA approach to overcome the proven shortcomings of PSHA, thus allowing for a reliable seismic hazard estimation, especially for those areas characterized by a prolonged quiescence, i.e. in tectonically active sites where events of only moderate size have occurred in historical times.     

Seismic response of L��Aquila downtown from comparison between 2D synthetics spectral ratios of SH, P-SV and Rayleigh waves and observations of the 2009 earthquake sequence

This is the first part of a study on the seismic response of the L’Aquila city using 2D simulation and experimental data. We have studied two velocity-depth models with the aim of outlining the behavior of a velocity reversal in the top layer, which is associated with the stiff Brecce de L’Aquila unit (BrA). In this setting, the SMTH model is topped by a layer with about 2:1 impedance contrast with the underlying layer while the NORV model has no velocity reversal. We have simulated the propagation of SH and P-SV wavefields in the range 0–10 Hz for incidence 0°–90°. Earthquake spectral ratios of the horizontal and vertical components at six sites in L’Aquila downtown are compared to corresponding synthetics spectral ratios. The vertical component of P-SV synthetics enables us to investigate a remarkable amplification effect seen in the vertical component of the recorded strong motion. Sites AQ04 and AQ05 are best matched by synthetics from the NORV model while FAQ5 and AQ06 have a better match with synthetics spectral ratios from the SMTH model. All simulations show this behavior systematically, with horizontal and near-horizontal incident waves predicting the overall pattern of matches more clearly than vertical and near-vertical incidence. The model inferences are in agreement with new geological data reporting lateral passages in the top layer from the stiff BrA to softer sediments. Matches are good in terms of frequency of the first amplification peak and of spectral amplitude: the horizontal components have spectral ratio peaks predominantly at 0.5 Hz in the simulations and at 0.7 Hz in the data, both with amplitudes of 4, while the vertical component spectral ratios reach values of 6 at frequencies of about 1 Hz in both data and simulations. The vertical component spectral ratios are very well matched using Rayleigh waves with incidence at 90°. The NORV model without the velocity reversal predicts spectral ratio peaks for the horizontal components at frequencies up to 6 Hz. The reversal of velocity acts as a low-pass frequency filter on the horizontal components reducing the amplification effect of the sediment filled valley.     

Synthetic tsunami mareograms for realistic oceanic models

We study how the tsunami mode is generated by a scaled double-couple seismic source, and how it propagates in realistic oceanic models. The method developed and used is the direct extension to tsunami waves propagating in multilayered oceanic media of the well-known Haskell method. The most intensive tsunamis may be expected from sources located within the crust in the deep-water parts of the ocean. The extension to laterally heterogeneous structures shows that, if the thickness of the ocean liquid layer diminishes, the maximum amplitude of the tsunami wave train increases.     

Extension of global creepex definition (Ms-mb) to local studies (Md-ML): the case of the Italian region

Creepex, defined as the systematic deviation of the magnitude of a single earthquake from the linear orthogonal regression between local magnitude M_L and coda duration magnitude M_d, calculated for the whole region, is used as a measure of the frequency content of the seismic sources in the Italian region. Predominantly high-frequency events are found in the two areas of Quaternary tectonic shortening in North-Central Italy and in the Calabrian Arc. This result, confirmed by two independent statistical tests, is in agreement with the global pattern obtained from the study of the regression between body-wave magnitude, m_b, and surface-wave magnitude M_s: systematic shift to high frequencies in the energy release of seismic sources located in subduction zones and to low frequencies in zones of spreading. The analysis of the correlation between the patterns of heat flow and of seismic source spectral properties indicates that these source properties, in general, do not reflect thermal conditions in the lithosphere, but rather represent the result of tectonic processes.     

Neo-deterministic definition of seismic input for residential seismically isolated buildings

This paper deals with the neo-deterministic definition of the seismic input in the municipality of Nimis (Italy), aimed at the design of residential seismically isolated buildings. The seismic input is defined by the computation of realistic synthetic seismograms considering different levels of detail for the earthquake source, both for a bedrock model and taking into account the specific site conditions. The horizontal response spectrum, calculated in the centre of the municipality by modelling the most dangerous source, advises against the construction of a building with a fixed base, but it is compatible with the seismic isolation, and it has been, therefore, used for the design of a residential seismically isolated building. The maximum displacement for the isolation system has been estimated about 17 cm, a value much lower than that provided by the code design response spectrum (28 cm). The importance of the realistic modelling, which permits the generalization of empirical observations by means of, physically sound, theoretical considerations, is evident, as it allows for the optimisation of the structural design with respect to the site of interest.