Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment

A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model.     

The ITalian ACcelerometric Archive (ITACA): processing of strong-motion data

The Italian strong-motion database was created during a joint project between Istituto Nazionale di Geofisica e Vulcanologia (INGV, Italian Institute for Geophysics and Vulcanology) and Dipartimento della Protezione Civile (DPC, Italian Civil Protection). The aim of the project was the collection, homogenization and distribution of strong motion data acquired in Italy in the period 1972–2004 by different institutions, namely Ente Nazionale per l’Energia Elettrica (ENEL, Italian electricity company), Ente per le Nuove tecnologie, l’Energia e l’Ambiente (ENEA, Italian energy and environment organization) and DPC. Recently the strong-motion data relative to the 23th December 2009, Parma (Mw = 5.4 and Mw = 4.9) and to the April 2009 L’Aquila sequences (13 earthquakes with 4.1 ≤ Mw ≤ 6.3) were included in the Italian Accelerometric Archive (ITACA) database (beta release). The database contains 7,038 waveforms from analog and digital instruments, generated by 1.019 earthquakes with magnitude up to 6.9 and can be accessed on-line at the web site . The strong motion data are provided in the unprocessed and processed versions. This article describes the steps followed to process the acceleration time series recorded by analogue and digital instruments. The procedures implemented involve: baseline removal, instrumental correction, band pass filtering with acausal filters, integration of the corrected acceleration in order to obtain velocity and displacement waveforms, computation of acceleration response spectra and strong motion parameters. This procedure is applied to each accelerogram and it is realised to preserve the low frequency content of the records.     

Sequence stratigraphy based on high-resolution seismic profiles in the late Pleistocene and Holocene deposits of the Venice area

The sequence-stratigraphic investigation by Very High-Resolution (VHR) seismic profiles allowed recognition of the detailed architecture of the late Pleistocene and Holocene succession of the Venice area. In this way deposits previously known by the analyses of scattered cores, mainly taken along the lagoon margin and the littoral strips, have been correlated at regional scale including the near offshore sector and the result has pointed out the lateral variability of the stratal architecture. Late Pleistocene deposits consist of an aggrading floodplain and fluvial channel fills accumulated during decreasing eustatic sea level, and they are coeval with offlapping forced regressive marine wedges in the Central Adriatic basin. The Holocene sequence is composed of three main seismic units separated by major stratal surfaces. Unit 1 (up to 9 m thick) is formed by channelized deposits separated by areas showing sub-horizontal and hummocky reflectors, and is bounded at the base by a surface that records prolonged conditions of subaerial exposure and at the top by a flatter surface resulting from erosion by marine processes. Deposits of Unit 1 are interpreted as estuarine and distributary channel fills, and back-barrier strata. Unit 2 is well distinguishable from Unit 1 only in the offshore area and at the barrier island bounding the Venice Lagoon, and is composed of a prograding marine wedge (up to 10 m thick) that interacts laterally with ebb tidal deltas. Unit 3 consists of a tidal channel complex and inlet deposits, which testify the evolution of the lagoon area. Tidal channels are entrenched in the lagoon mud flat (coeval with Units 1–2) and cut the Pleistocene–Holocene boundary in several places.Following current sequence-stratigraphic concepts, the Holocene sequence is composed of a paralic transgressive systems tract (TST) (Unit 1) overlying a sequence boundary (the Pleistocene–Holocene boundary) and overlain by a marine highstand systems tract (HST) (Unit 2) in seaward locations and by highstand lagoonal deposits landwards. TST and HST are separated by a downlap surface that is amalgamated with a wave ravinement surface in several places. Unit 3 is coeval with the upper part of Unit 2, and its development has been favoured by human interventions, which led to a transgression limited to the lagoon area.Local factors during the deposition, i.e. subsidence, sediment supply, physiography, and current/wave regimes, led to a significant lateral variability in the architecture of the Holocene sequence, as evidenced by the extreme thickness variation of the TST along both depositional strike and dip. The HST, instead, shows less pronounced strike variations in the stratal architecture. Also, present data clearly evidence that the human impact has a great relevance in influencing the late Holocene sedimentation.     

Is there a correlation between lithosphere structure and the statistical properties of seismicity?

The elastic parameters and the thickness of the lithospheric part of the mantle, the lid, are correlated with the seismicity of Italy and surrounding areas through the estimation of several statistical parameters. In particular, the intensity of the flow of seismic events, deprived of aftershocks, reveals concentrations associated with large gradients in the lithospheric thickness. The average depth of crustal shocks tends to be very shallow in areas of thin and soft (low shear-wave velocities) lid. The positive influence, probably reflecting the thermal state in the lithosphere, is very large in connection with the thin and soft lid and where strong gradients in lid thickness are observed, and is small in connection with a thick and hard lid (high shear-wave velocities).     

Source Parameters of Weak Crustal Earthquakes of the Vrancea Region from Short-period Waveform Inversion

High-frequency records of nine low magnitude shallow earthquakes of the Vrancea (Romania) seismic region are inverted for the seismic moment tensor (MT). An approach is suggested regarding how to obtain at least a rough estimate of the MT when the information on the structure of the crust is poor. Here simple 1-D layered models are used in the Greens function synthesis despite the fact that the structure of the region is undoubtedly very complex. Different 1-D models were used for different source-station paths to approximate lateral variations. Record of a station located on a ray path which crosses a structure differing substantially from a 1-D model may however bias the retrieved MT essentially. Therefore, we did not collectively process all records, but subsets of stations separately. We check the consistency of the MTs resulting from these individual bootstrap solutions, and reject those which differ substantially, assuming that this is due to the oversimplification of the forward modeling. Thereafter, the averaged moment tensor yielded by the consistent subsets is accepted. Moreover, the distribution of the T, N and P axes from the moment tensors used for averaging provides a rough estimate of their reliability.Following this simplistic procedure, we found seven acceptably constrained solutions among nine events processed. Their P axes are compared with the general trend of the stress in the area: three comply with it, the others should be attributed to the complex stress field active in the region.Acknowledgment This research has been made possible by MURST (40% and 60%), by UNESCO-IGCP project 414 Realistic Modelling of Seismic Input for Megacities and Large Urban Areas and NATO SfP 972266. Partial support issued from the Grant Agency of Czech Rep. (Grant 205/02/0383) and from the National Agency for Science, Technology and Innovation of Romania (Grant 6185 GR/2000).     

Three Decades of Seismic Activity at Mt. Vesuvius: 1972–2000

We analyse the seismic catalogue of the local earthquakes which occurred at Somma-Vesuvius volcano in the past three decades (1972–2000). The seismicity in this period can be described as composed of a background level, characterised by a low and rather uniform rate of energy release and by sporadic periods of increased seismic activity. Such relatively intense seismicity periods are characterised by energy rates and magnitudes progressively increasing in the critical periods. The analysis of the b value in the whole period evidences a well-defined pattern, with values of b progressively decreasing, from about 1.8 at the beginning of the considered period, to about 1.0 at present. This steady variation indicates an increasing dynamics in the volcanic system. Within this general trend it is possible to identify a substructure in the time sequence of the seismic events, formed by the alternating episodes of quiescence and activity. The analysis of the source moment tensor of the largest earthquakes shows that the processes at the seismic source are generally not consistent with simple double-couples, but that they are compatible with isotropic components, mostly indicating volumetric expansion. These components are shown to be statistically significant for most of the analysed events. Such focal mechanisms can be interpreted as the effect of explosion phenomena, possibly related to volatile exsolution from the crystallising magma. The availability of a reduced amount of high quality data necessary for the inversion of the source moment tensor, the still limited period of systematic observation of Vesuvius micro-earthquakes and, above all, the absence of eruptive events during such interval of time, cannot obviously permit the outlining of any formal premonitory signal. Nevertheless, the analysis reported in this paper indicates a progressively evolving dynamics, characterised by a generally increasing trend in the seismic activity in the volcanic system and by a significant volumetric component of recent major events, thus posing serious concern for a future evolution towards eruptive activity.