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]

The dynamic behaviour and the seismic vulnerability of the ancient civic clock tower of Amatrice, dramatically damaged by the last shocks sequence of 2016 that occurred in Central Italy, is studied in this paper by means of advanced 3D numerical analyses with the Non-Smooth Contact Dynamics Method (NSCD). Thus, a discontinuous approach was used to assess the dynamic properties and the vulnerability of the masonry structure, through large displacements regulated by the Signorini’s law, concerning the impenetrability between the rigid bodies, and by the Coulomb’s law, regarding the dry-friction model. Afterward, different values were assigned to the friction coefficient of the models and real seismic shocks were then applied in nonlinear analyses. The major purpose of this study is to highlight that relevant data on the real structural behaviour of historical masonry may provided through advanced numerical analyses. The comparison between the results of numerical simulation and survey of the existing crack pattern of the bell tower permitted to validate the selected approach. Finally, from results and conclusions provided by this case study, it is possible to affirm that the used methodology may be applied to a wide variety of historical masonry structure in Europe. [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]