Abstract

Abstract. During the 2012 seismic sequence of the Emilia region (northern Italy), the earthquake ground motion in the epicentral area featured longer duration and higher velocity than those estimated by empirical-based prediction equations typically adopted in Italy. In order to explain these anomalies, we (1) build up a structural and geophysical 3-D digital model of the crustal sector involved in the sequence, (2) reproduce the earthquake ground motion at some seismological stations through physics-based numerical simulations and (3) compare the observed recordings with the simulated ones. In this way, we investigate how the earthquake ground motion in the epicentral area is influenced by local stratigraphy and geological structure buried under the Po Plain alluvium. Our study area covers approximately 5000 km2 and extends from the right Po River bank to the Northern Apennine morphological margin in the N–S direction, and between the two chief towns of Reggio Emilia and Ferrara in the W–E direction, involving a crustal volume of 20 km thickness. We set up the 3-D model by using already-published geological and geophysical data, with details corresponding to a map at scale of 1:250 000. The model depicts the stratigraphic and tectonic relationships of the main geological formations, the known faults and the spatial pattern of the seismic properties. Being a digital vector structure, the 3-D model can be easily modified or refined locally for future improvements or applications. We exploit high-performance computing to perform numerical simulations of the seismic wave propagation in the frequency range up to 2 Hz. In order to get rid of the finite source effects and validate the model response, we choose to reproduce the ground motion related to two moderate-size aftershocks of the 2012 Emilia sequence that were recorded by a large number of stations. The obtained solutions compare very well to the recordings available at about 30 stations in terms of peak ground velocity and signal duration. Snapshots of the simulated wavefield allow us to attribute the exceptional length of the observed ground motion to surface wave overtones that are excited in the alluvial basin by the buried ridge of the Mirandola anticline. Physics-based simulations using realistic 3-D geomodels show eventually to be effective for assessing the local seismic response and the seismic hazard in geologically complex areas.

Highlights

  • Computer-aided three-dimensional (3-D) geological modeling (e.g., Mallet, 2002) is becoming an increasingly important tool in geoscience studies for both the management of natural resources and the prevention of natural disasters. 3-D geological modeling allows the combination of multidisciplinary data in the shaping and visualization of the current knowledge of the geological structures and allows integration with new data or interpretations, as they become available (Calcagno, 2015)

  • The adequateness of the FPSM3D code in this kind of applications is demonstrated in the works by Chaljub et al (2015) and Maufroy et al (2015), aimed to estimate the accuracy of a number of numerical methods currently used for physicsbased predictions of earthquake ground motion in 3-D models of sedimentary basins

  • In order to investigate whether the ER3D model is able to reproduce the peculiar features of the observed earthquake ground motion, we performed a comparison between the ground motion recorded by 29 seismological stations deployed in the study area during the 2012 seismic sequence – see Fig. 1 and Table 4 – and the numerically predicted ground motion at the same locations

Read more

Summary

Introduction

Computer-aided three-dimensional (3-D) geological modeling (e.g., Mallet, 2002) is becoming an increasingly important tool in geoscience studies for both the management of natural resources and the prevention of natural disasters. 3-D geological modeling allows the combination of multidisciplinary data in the shaping and visualization of the current knowledge of the geological structures and allows integration with new data or interpretations, as they become available (Calcagno, 2015). An emblematic case of such deviations occurred during the 2012 Emilia seismic crisis, when unexpectedly long duration and large peak ground velocity (PGV) characterized the earthquake ground motion at some sites in the epicentral area (Priolo et al, 2012; Castro et al, 2013; Luzi et al, 2013; Barnaba et al, 2014; De Nardis et al, 2014) Those deviations have been attributed to the complexity of the geological structure beneath the Po Plain, which features a very large and deep alluvial basin bounded by two largely buried thrust-and-fold systems, the Northern Apennine chain in the south and the southern Alpine ridge in the north, respectively (Boccaletti et al, 1985). The computations were run using the HPC resources of the CINECA consortium in Bologna

The structural and geophysical 3-D model of central Emilia
Geological and seismotectonic setting of the study area
Integration of geological data in the 3-D digital model
Physical properties of geological formations
The FPSM3D code
Setup for the computations
Comparisons between numerical predictions and data
The simulated earthquake ground motion
Comparison with the empirical earthquake ground motion
Wave propagation across a vertical profile
Findings
Conclusions

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.