Search

The preparation, initiation and occurrence dynamics of earthquakes in Italy are governed by several frequently unknown physical mechanisms and parameters. Understanding these mechanisms is crucial for developing new techniques and approaches for earthquake monitoring and hazard assessments. Here, we develop a first-order numerical model simulating quasi-static crustal interseismic loading, coseismic brittle episodic dislocation and post-seismic relaxation for extensional and compressional earthquakes in Italy based on a common framework of lithostatic and tectonic forces. Our model includes an upper crust, where the fault is locked, and a deep crust, where the fault experiences steady shear. The results indicate that during the interseismic phase, the contrasting behaviour between the upper locked fault segment and lower creeping fault segment generates a stretched volume at depth in the hangingwall via extensional tectonics while a contracted volume forms via compressional tectonics. The interseismic stress and strain gradients invert at the coseismic stage, with the interseismic dilated volume contracting during the coseismic stage, and vice versa. Moreover, interseismic stress gradients promote coseismic gravitational subsidence of the hangingwall for normal fault earthquakes and elastic uplift for reverse fault earthquakes. Finally, the post-seismic relaxation is characterized by further ground subsidence and uplift for normal and reverse faulting earthquakes, respectively, which is consistent with the faulting style. The fault is the passive feature, with slipping generating the seismic waves, whereas the energy activating the movement is stored mostly in the hangingwall volume. The main source of energy for normal faulting and thrust is provided by the lithostatic load and elastic load, respectively. [ABSTRACT FROM AUTHOR]

Stimulated by the recent advances in computational tools for the simulation of seismic wave propagation problems in realistic geological configurations, this paper presents a 3D physicsbased numerical analysis of near-source ground motion during the MW 6.0 2012 May 29 earthquake in the Po Plain, Northern Italy. To reproduce with sufficient accuracy some of the most peculiar features of the near-source strong-motion records and of the spatial variability of damage distribution, this study required a sequence of investigations, starting from the analysis of a wide set of near-source records, to the calibration of an improved kinematic seismic source model, up to the development of a 3D numerical model of the portion of the Po Plain interested by the earthquake. The latter includes the basin geometry, characterized by sediment thickness sharply varying from few tens of metres to some kilometres. The spatial resolution of the numerical model is suitable to propagate frequencies up to about 1.5 Hz. Numerical simulations were performed using the open-source high-performance code SPEED, based on the Discontinuous Galerkin Spectral Elements method. The 3D numerical model, coupled with the updated slip distribution along the rupturing fault, proved successful to reproduce with good agreement, measured through quantitative goodness-of-fit criteria, the most relevant features of the observed ground motion. These include: (i) the large fault normal velocity peaks at the near-source stations driven by updip directivity effects; (ii) the small-scale variability at short distance from the source, resulting in the out-of-phase motion at stations separated by only 3 km distance; (iii) the propagation of prominent trains of surface waves, especially in the Northern direction; (iv) the map of earthquake-induced ground uplift with maximum values of about 10 cm, in substantial agreement with geodetic measurements and (v) the two-lobed pattern of the peak ground velocity map, well correlated with the distribution of macroseismic intensity. [ABSTRACT FROM AUTHOR]