This work presents an up-to-date model for the simulation of non-stationary ground motions, including several novelties compared to the original study of Sabetta and Pugliese (Bull Seism Soc Am ...86:337–352, 1996). The selection of the input motion in the framework of earthquake engineering has become progressively more important with the growing use of nonlinear dynamic analyses. Regardless of the increasing availability of large strong motion databases, ground motion records are not always available for a given earthquake scenario and site condition, requiring the adoption of simulated time series. Among the different techniques for the generation of ground motion records, we focused on the methods based on stochastic simulations, considering the time- frequency decomposition of the seismic ground motion. We updated the non-stationary stochastic model initially developed in Sabetta and Pugliese (Bull Seism Soc Am 86:337–352, 1996) and later modified by Pousse et al. (Bull Seism Soc Am 96:2103–2117, 2006) and Laurendeau et al. (Nonstationary stochastic simulation of strong ground-motion time histories: application to the Japanese database. 15 WCEE Lisbon, 2012). The model is based on the S-transform that implicitly considers both the amplitude and frequency modulation. The four model parameters required for the simulation are: Arias intensity, significant duration, central frequency, and frequency bandwidth. They were obtained from an empirical ground motion model calibrated using the accelerometric records included in the updated Italian strong-motion database ITACA. The simulated accelerograms show a good match with the ground motion model prediction of several amplitude and frequency measures, such as Arias intensity, peak acceleration, peak velocity, Fourier spectra, and response spectra.
In this paper the site categorization criteria and the corresponding site amplification factors proposed in the 2021 draft of Part 1 of Eurocode 8 (2021-draft, CEN/TC250/SC8 Working Draft N1017) are ...first introduced and compared with the current version of Eurocode 8, as well as with site amplification factors from recent empirical ground motion prediction equations. Afterwards, these values are checked by two approaches. First, a wide dataset of strong motion records is built, where recording stations are classified according to 2021-draft, and the spectral amplifications are empirically estimated computing the site-to-site residuals from regional and global ground motion models for reference rock conditions. Second, a comprehensive parametric numerical study of one-dimensional (1D) site amplification is carried out, based on randomly generated shear-wave velocity profiles, classified according to the new criteria. A reasonably good agreement is found by both approaches. The most relevant discrepancies occur for the shallow soft soil conditions (soil category E) that, owing to the complex interaction of shear wave velocity, soil deposit thickness and frequency range of the excitation, show the largest scatter both in terms of records and of 1D numerical simulations. Furthermore, 1D numerical simulations for soft soil conditions tend to provide lower site amplification factors than 2021-draft, as well as lower than the corresponding site-to-site residuals from records, because of higher impact of non-linear (NL) site effects in the simulations. A site-specific study on NL effects at three KiK-net stations with a significantly large amount of high-intensity recorded ground motions gives support to the 2021-draft NL reduction factors, although the very limited number of recording stations allowing such analysis prevents deriving more general implications. In the presence of such controversial arguments, it is reasonable that a standard should adopt a prudent solution, with a limited reduction of the site amplification factors to account for NL soil response, while leaving the possibility to carry out site-specific estimations of such factors when sufficient information is available to model the ground strain dependency of local soil properties.
The Engineering Strong-Motion (ESM) flatfile is a parametric table which contains verified and reliable metadata and intensity measures of manually processed waveforms included in the ESM database. ...The flatfile has been developed within the Seismology Thematic Core Service of EPOS-IP (European Plate Observing System Implementation Phase) and it is disseminated throughout a web portal (
http://esm.mi.ingv.it/flatfile-2018/flatfile.php
) for research and technical purposes. The adopted criteria for flatfile compilation aim to collect strong motion data and related metadata in a uniform, updated, traceable and quality-checked way to develop Ground Motion Models (GMMs) for Probabilistic Seismic Hazard Assessment (PSHA) and engineering applications. In this paper, we present the characteristics of ESM flatfile in terms of recording, event and station distributions, and we discuss the most relevant features of the Intensity Measures (IMs) of engineering interest included in the table. The dataset for flatfile compilation includes 23,014 recordings from 2179 earthquakes and 2080 stations from Europe and Middle-East. The events are characterized by magnitudes in the range 3.5–8.0 and refer to different tectonics regimes, such as shallow active crustal and subduction zones. Intensity measures include peak and integral parameters and duration of each waveform. The spectral amplitudes of the (5% damping) acceleration and displacement response are provided for 36 periods, in the interval 0.01–10 s, as well as the 103 amplitudes of the Fourier spectrum for the frequency range 0.04–50 Hz. Several statistics are shown with reference to the most significant metadata for GMMs calibrations, such as moment magnitude, focal depth, several distance metrics, style of faulting and parameters for site characterization. Furthermore, we also compare and explain the most relevant differences between the metadata of ESM flatfile with those provided by the previous flatfile derived in RESORCE (Reference Database for Seismic Ground Motion in Europe) project.
This work offers a novel methodological framework to address the problem of generating data-driven earthquake shaking fields at different vibration periods, which are key to support decision making ...and civil protection planning. We propose to analyse the entire profiles of spectral accelerations and project their information content to unsampled locations in the system, based on the theory of Object Oriented Spatial Statistics. The proposed methodology combines a non-ergodic ground motion model with a fully functional model for the residual term, the latter consisting of (i) the spatially-varying systematic effects due to source, site and path, and (ii) the remaining aleatory error. The proposed methodology allows to generate multiple shaking scenarios conditioned on the data, jointly and consistently for all the vibration periods, overcoming the intrinsic limitations of existing multivariate approaches to the problem. The approach is tested on a vast dataset of ground motion records collected in the study-area of the Po Plain (Northern Italy), for which a region-specific fully non-ergodic GMM was previously calibrated. Our validation tests demonstrate the potentiality of the approach, which is capable to effectively simulate spectral acceleration profiles, while keeping the ability to capture the main physical features of ground motion patterns in the region.
ShakeMap is the tool to evaluate the ground motion effect of earthquakes in vast areas. It is useful to delimit the zones where the shaking is expected to have been most significant, for civil ...defense rapid response. From the earthquake engineering point of view, it can be used to infer the seismic actions on the built environment to calibrate vulnerability models or to define the reconstruction policies based on observed damage vs shaking. In the case of long-lasting seismic sequences, it can be useful to develop ShakeMap envelopes, that is, maps of the largest ground intensity among those from the ShakeMap of (selected) events of a seismic sequence, to delimit areas where the effects of the whole sequence have been of structural engineering relevance. This study introduces ShakeMap envelopes and discusses them for the central Italy 2016–2017 seismic sequence. The specific goals of the study are: (i) to compare the envelopes and the ShakeMap of the main events of the sequence to make the case for sequence-based maps; (ii) to quantify the exceedance of design seismic actions based on the envelopes; (iii) to make envelopes available for further studies and the reconstruction planning; (iv) to gather insights on the (repeated) exceedance of design seismic actions at some sites. Results, which include considerations of uncertainty in ShakeMap, show that the sequence caused exceedance of design hazard in thousands of square kilometers. The most relevant effects of the sequence are, as expected, due to the mainshock, yet seismic actions larger than those enforced by the code for structural design are found also around the epicenters of the smaller magnitude events. At some locations, the succession of ground-shaking that has excited structures, provides insights on structural damage accumulation that has likely taken place; something that is not accounted for explicitly in modern seismic design. The envelopes developed are available as supplemental material.
We present the results of a consistency check performed over a flatfile of accelerometric data extracted from the ITalian ACcelerometric Archive (ITACA), enriched with velocimetric records of events ...with magnitude M < 4.0. The flatfile, called ITACAext, includes 31,967 waveforms from 1709 shallow crustal earthquakes, in the magnitude range from 3.0 to 6.9, and occurred in the period of 1972–2019 in Italy. The consistency check is carried out by decomposing the residuals obtained from a reference ground motion model, for the ordinates of the 5% damped acceleration response spectra. The residual components are subsequently analyzed to identify a list of events, stations, and records that significantly deviate from the median trends predicted by the model. The results indicate that about 10% of events and stations are outliers, while only 1% of the waveforms present anomalous amplitudes. The asymmetrical azimuthal coverage of seismic stations around the epicenter is the most common issue that can affect the estimates of the repeatable event residual term. On the other hand, peculiarities in the site-response or wrong estimates of the soil parameters (i.e., the average shear-wave velocity in the first 30 m of the subsoil) are the main issues related to the repeatable station residuals. Finally, single records can show large residuals because of issues related to signal acquisition (e.g., multiple events, noisy records) or possible near-source effects (e.g., rupture directivity).
The aim of this paper is to identify the ground motion models (GMMs), applicable in active shallow crustal regions (ASCRs) and subduction zones (SZs), to be used for the new release of the seismic ...hazard model of Italy (MPS19) for peak ground acceleration and 12 ordinates of the acceleration response spectra (5% damping) in the period range 0.05–4 s. The steps to achieve such goal are: (1) a pre-selection of the GMMs that takes into account the suitability for the application to the Italian territory and the fulfillment of the new hazard model requirements; (2) the assessment of a proper scoring of the pre-selected GMMs using strong-motion data recorded in Italy; (3) the selection of the GMMs to be used in the hazard calculation. The final set of GMMs describes satisfactorily the epistemic uncertainty of the ground shaking process, privileging the simplicity and flexibility of the functional form. Finally, the weights of the selected GMMs are assigned combining the results of the scoring and the weights obtained through an experts’ elicitation.
This work describes the analysis of the strong-motion data from the Engineering Strong Motion database (ESM,
http://esm.mi.ingv.it
), aimed at: (1) extract a dataset of accelerometric waveforms ...recorded during the 2016–2017 Central Italy seismic sequence; (2) identify the recording stations to be used as reference sites for further seismological analysis; (3) select the records to be used as input for seismic microzonation of higher level at 137 municipalities. Firstly, a residual analysis is carried out on the extracted dataset to perform: (1) the quality check of the waveforms recorded by temporary networks installed soon after the occurrence of the first main shock (M 6.0, 24 August 2016); (2) the estimation of the site-to-site residual term for each recording station with the aim of recognising potential reference rock sites. Finally, the software REXELite, integrated within the ESM website, is adopted to select suites of spectrum-compatible accelerograms, that will be used as input for calculating site amplifications through 1D and 2D simulations at sites which suffered the greatest damage. The results of this work demonstrate the success of the synergy among Italian institutions. The setup of key infrastructures, such as emergency networks and data repositories, together with the knowledge developed during national projects, turned out to be successful in terms of timely intervention during the emergency phase and the planning of the post-emergency.
Near-source effects can amplify seismic ground motion, causing large demand to structures and thus their identification and characterization is fundamental for engineering applications. Among the ...most relevant features, forward-directivity effects may generate near-fault records characterized by a large velocity pulse and unusual response spectral shape amplified in a narrow frequency-band. In this paper, we explore the main statistical features of acceleration and displacement response spectra of a suite of 230 pulse-like signals (impulsive waveforms) contained in the NESS1 (NEar Source Strong-motion) flat-file. These collected pulse-like signals are analyzed in terms of pulse period and pulse azimuthal orientation. We highlight the most relevant differences of the pulse-like spectra compared to the ordinary (i.e., no-pulse) ones, and quantify the contribution of the pulse through a corrective factor of the spectral ordinates. Results show that the proposed empirical factors are able to capture the amplification effect induced by near-fault directivity, and thus they could be usefully included in the framework of probabilistic seismic hazard analysis to adjust ground-motion model (GMM) predictions.