Crustal structure of northern Italy from the ellipticity of Rayleigh waves

Abstract

Northern Italy is a diverse geological region, including the wide and thick Po Plain sedimentary basin, which is bounded by the Alps and the Apennines. The seismically slow shallow structure of the Po Plain is difficult to retrieve with classical seismic measurements such as surface wave dispersion, yet the detailed structure of the region greatly affects seismic wave propagation and hence seismic ground shaking.

Here we invert Rayleigh wave ellipticity measurements in the period range 10–60 s for 95 stations in northern Italy using a fully non linear approach to constrain vertical and density profiles of the crust beneath each station. The ellipticity of Rayleigh wave ground motion is primarily sensitive to shear-wave velocity beneath the recording station, which reduces along-path contamination effects. We use the 3D layering structure in MAMBo, a previous model based on a compilation of geological and geophysical information for the Po Plain and surrounding regions of northern Italy, and employ ellipticity data to constrain and density within its layers. We show that ellipticity data from ballistic teleseismic wave trains alone constrain the crustal structure well. This leads to MAMBo-E, an updated seismic model of the region’s crust that inherits information available from previous seismic prospection and geological studies, while fitting new seismic data well. MAMBo-E brings new insights into lateral heterogeneity in the region’s subsurface. Compared to MAMBo, it shows overall faster seismic anomalies in the region’s Quaternary, Pliocene and Oligo-Miocene layers and better delineates the seismic structures of the Po Plain at depth. Two low velocity regions are mapped in the Mesozoic layer in the western and eastern parts of the Plain, which seem to correspond to the Monferrato sedimentary basin and to the Ferrara-Romagna thrust system, respectively.