scholarly journals Estimating Ground-Motion Models via Bayesian Inference using Stan

2021 ◽  
Author(s):  
Nicolas Kuehn ◽  
Peter Stafford

We provide a simple introduction to the estimation of ground-motion models via Bayesian inference and the probabilistic programming language Stan.We show one ca implement a simple ground-motion model in Stan, and how can run the program from the computer environment R.We also show how one can access the results, and plot summaries of estimated parameters.A large number of different Stan models for the development https://github.com/pstafford/StanGMMTutorial.

2021 ◽  
pp. 875529302110552
Author(s):  
Silvia Mazzoni ◽  
Tadahiro Kishida ◽  
Jonathan P Stewart ◽  
Victor Contreras ◽  
Robert B Darragh ◽  
...  

The Next-Generation Attenuation for subduction zone regions project (NGA-Sub) has developed data resources and ground motion models for global subduction zone regions. Here we describe the NGA-Sub database. To optimize the efficiency of data storage, access, and updating, data resources for the NGA-Sub project are organized into a relational database consisting of 20 tables containing data, metadata, and computed quantities (e.g. intensity measures, distances). A database schema relates fields in tables to each other through a series of primary and foreign keys. Model developers and other users mostly interact with the data through a flatfile generated as a time-stamped output of the database. We describe the structure of the relational database, the ground motions compiled for the project, and the means by which the data can be accessed. The database contains 71,340 three-component records from 1880 earthquakes from seven global subduction zone regions: Alaska, Central America and Mexico, Cascadia, Japan, New Zealand, South America, and Taiwan. These data were processed on a component-specific basis to minimize noise effects in the data and remove baseline drifts. Provided ground motion intensity measures include peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations for a range of oscillator periods.


2008 ◽  
Vol 24 (1) ◽  
pp. 299-317 ◽  
Author(s):  
Jack W. Baker ◽  
Nirmal Jayaram

Ground motion models (or “attenuation relationships”) describe the probability distribution of spectral acceleration at an individual period, given a set of predictor variables such as magnitude and distance, but they do not address the correlations between spectral acceleration values at multiple periods or orientations. Those correlations are needed for several calculations related to seismic hazard analysis and ground motion selection. Four NGA models and the NGA ground motion database are used here to measure these correlations, and predictive equations are fit to the results. The equations are valid for periods from 0.01 seconds to 10 seconds, versus similar previous equations that were valid only between 0.05 and 5 seconds and produced unreasonable results if extrapolated. Use of the new NGA ground motion database also facilitates a first study of correlations from intra- and inter-event residuals. Observed correlations are not sensitive to the choice of accompanying ground motion model, and intra-event, inter-event, and total residuals all exhibit similar correlation structure. A single equation is thus applicable for a variety of correlation predictions. A simple example illustrates the use of the proposed equations for one hazard analysis application.


2021 ◽  
Author(s):  
Nicolas Kuehn

A nonrgodic ground-motion model explicitly takes systematic local source, path and site effects on the predicted ground-motion into account. With an increasing number of ground-motion records, it s possible to estimate these effects. Landwehr et al. (2016) proposed a varying coefficient model as a tool to estimate nonergodic ground-motion models, based on Gaussian processes. Gaussian processes are computationally expen- sive, so for large data sets, some approximations have to be made. Here, we compare different Bayesian implementations of varying coefficient models, using the probabilistic programming language Stan (Carpenter et al., 2017), and the integrated nested Laplace approximation (INLA) (Rue et al., 2009). The models are used to fit nonergodic models on the California subset of the NGA-West2 data set (Ancheta et al., 2014). We find that both implementations lead to very similar results, both in the estimated parameters and the predicted ground motions.


Author(s):  
Fabio Sabetta ◽  
Antonio Pugliese ◽  
Gabriele Fiorentino ◽  
Giovanni Lanzano ◽  
Lucia Luzi

AbstractThis 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.


Sign in / Sign up

Export Citation Format

Share Document