skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Conditioned Simulation of Ground-Motion Time Series at Uninstrumented Sites Using Gaussian Process Regression
ABSTRACT Ground-motion time series are essential input data in seismic analysis and performance assessment of the built environment. Because instruments to record free-field ground motions are generally sparse, methods are needed to estimate motions at locations with no available ground-motion recording instrumentation. In this study, given a set of observed motions, ground-motion time series at target sites are constructed using a Gaussian process regression (GPR) approach, which treats the real and imaginary parts of the Fourier spectrum as random Gaussian variables. Model training, verification, and applicability studies are carried out using the physics-based simulated ground motions of the 1906 Mw 7.9 San Francisco earthquake and Mw 7.0 Hayward fault scenario earthquake in northern California. The method’s performance is further evaluated using the 2019 Mw 7.1 Ridgecrest earthquake ground motions recorded by the Community Seismic Network stations located in southern California. These evaluations indicate that the trained GPR model is able to adequately estimate the ground-motion time series for frequency ranges that are pertinent for most earthquake engineering applications. The trained GPR model exhibits proper performance in predicting the long-period content of the ground motions as well as directivity pulses.  more » « less
Award ID(s):
2025310
PAR ID:
10321638
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
BSSA
Date Published:
Journal Name:
Bulletin of the Seismological Society of America
Volume:
112
Issue:
1
ISSN:
0037-1106
Format(s):
Medium: X Other: BSSA
Sponsoring Org:
National Science Foundation
More Like this
  1. Following an earthquake, ground motion time series are needed to carry out site-specific nonlinear response history analysis. However, the number of currently available recording instruments is sparse; thus, the ground motion time series at uninstrumented sites must be estimated. Tamhidi et al. developed a Gaussian process regression (GPR) model to generate ground motion time series given a set of recorded ground motions surrounding the target site. This GPR model interpolates the observed ground motions’ Fourier Transform coefficients to generate the target site’s Fourier spectrum and the corresponding time series. The robustness of the optimized hyperparameter of the model depends on the surrounding observation density. In this study, we carried out sensitivity analysis and tuned the hyperparameter of the GPR model for various observation densities. The 2019 M7.1 Ridgecrest and 2020 M4.5 South El Monte earthquake data sets recorded by the Community Seismic Network and California Integrated Seismic Network in Southern California are used to demonstrate the process. To provide a tool to quantify the uncertainty of the generated motions, a methodology to develop realizations of ground motion time series is also incorporated. The results illustrate that the uncertainty of the generated motions is lower at longer periods. It is shown that the observation density in the proximity of the target site plays a vital role in both error and uncertainty reduction of the generated time series. To demonstrate the concept, the effect of additional observations from combined recording networks is investigated. 
    more » « less
  2. Long‐period ground motions from large (Mw≥7.0) subduction‐zone earthquakes are a real threat for large‐scale human‐made structures. The Nankai subduction zone, Japan, is expected to host a major megathrust earthquake in the near future and has therefore been instrumented with offshore and onshore permanent seismic networks. We use the ambient seismic field continuously recorded at these stations to simulate the long‐period (4–10 s) ground motions from past and future potential offshore earthquakes. First, we compute impulse response functions (IRFs) between an ocean‐bottom seismometer of the Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) network, which is located offshore on the accretionary wedge, and 60 onshore Hi‐net stations using seismic interferometry by deconvolution. As this technique only preserves the relative amplitude information of the IRFs, we use a moderate Mw 5.5 event to calibrate the amplitudes to absolute levels. After calibration, the IRFs are used together with a uniform stress‐drop source model to simulate the long‐period ground motions of the 2004 Mw 7.2 intraplate earthquake. For both events, the residuals of the 5\\% damped spectral acceleration (SA) computed from the horizontal and vertical components of the observed and simulated waveforms exhibit almost no bias and acceptable uncertainties. We also compare the observed SA values of the Mw 7.2 event to those from the subduction‐zone BC Hydro ground‐motion model (GMM) and find that our simulations perform better than the model. Finally, we simulate the long‐period ground motions of a hypothetical Mw 8.0 subduction earthquake that could occur along the Nankai trough. For this event, our simulations generally exhibit stronger long‐period ground motions than those predicted by the BC Hydro GMM. This study suggests that the ambient seismic field recorded by the ever‐increasing number of ocean‐bottom seismometers can be used to simulate the long‐period ground motions from large megathrust earthquakes. 
    more » « less
  3. ABSTRACT CyberShake is a high-performance computing workflow for kinematic fault-rupture and earthquake ground-motion simulation developed by the Statewide California Earthquake Center to facilitate physics-based probabilistic seismic hazard assessment (PSHA). CyberShake exploits seismic reciprocity for wave propagation by computing strain green tensors along fault planes, which in turn are convolved with rupture models to generate surface seismograms. Combined with a faultwide hypocentral variation of each simulated rupture, this procedure allows for generating ground-motion synthetics that account for realistic source variability. This study validates the platform’s kinematic modeling of physics-based seismic wave propagation simulations in Southwest Iceland as the first step toward migrating CyberShake from its original study region in California. Specifically, we have implemented CyberShake workflows to model 2103 fault ruptures and simulate the corresponding two horizontal components of ground-motion velocity on a 5 km grid of 625 stations in Southwest Iceland. A 500-yr-long earthquake rupture forecast consisting of 223 hypothetical finite-fault sources of Mw 5–7 was generated using a physics-based model of the bookshelf fault system of the Southwest Iceland transform zone. For each station, every reciprocal simulation uses 0–1 Hz Gaussian point sources polarized along two horizontal grid directions. Comparison of the results in the form of rotation-invariant synthetic pseudoacceleration spectral response values at 3, 4, and 5 s periods are in good agreement with the Icelandic strong motion data set and a suite of empirical Bayesian ground-motion prediction equations (GMPEs). The vast majority of the physics-based simulations fall within one standard deviation of the mean GMPE predictions, previously estimated for the area. At large magnitudes for which no data exist in Iceland, the synthetic data set may play an important role in constraining GMPEs for future applications. Our results comprise the first step toward comprehensive and physics-based PSHA for Southwest Iceland. 
    more » « less
  4. ABSTRACT Earthquake ground motions in the vicinity of receivers couple with the atmosphere to generate pressure perturbations that are detectable by infrasound sensors. These so-called local infrasound signals traverse very short source-to-receiver paths, so that they often exhibit a remarkable correlation with seismic velocity waveforms at collocated seismic stations, and there exists a simple relationship between vertical seismic velocity and pressure time series. This study leverages the large regional network of infrasound sensors in Alaska to examine local infrasound from several light to great Alaska earthquakes. We estimate seismic velocity time series from infrasound pressure records and use these converted infrasound recordings to compute earthquake magnitudes. This technique has potential utility beyond the novelty of recording seismic velocities on pressure sensors. Because local infrasound amplitudes from ground motions are small, it is possible to recover seismic velocities at collocated sites where the broadband seismometers have clipped. Infrasound-derived earthquake magnitudes exhibit good agreement with seismically derived values. This proof-of-concept demonstration of computing seismic magnitudes from infrasound sensors illustrates that infrasound sensors may be utilized as proxy vertical-component seismometers, making a new data set available for existing seismic techniques. Because single-sensor infrasound stations are relatively inexpensive and are becoming ubiquitous, this technique could be used to augment existing regional seismic networks using a readily available sensor platform. 
    more » « less
  5. Ground motion selection has become increasingly central to the assessment of earthquake resilience. The selection of ground motion records for use in nonlinear dynamic analysis significantly affects structural response. This, in turn, will impact the outcomes of earthquake resilience analysis. This paper presents a new ground motion clustering algorithm, which can be embedded in current ground motion selection methods to properly select representative ground motion records that a structure of interest will probabilistically experience. The proposed clustering-based ground motion selection method includes four main steps: 1) leveraging domain-specific knowledge to pre-select candidate ground motions; 2) using a convolutional autoencoder to learn low-dimensional underlying characteristics of candidate ground motions’ response spectra – i.e., latent features; 3) performing k-means clustering to classify the learned latent features, equivalent to cluster the response spectra of candidate ground motions; and 4) embedding the clusters in the conditional spectra-based ground motion selection. The selected ground motions can represent a given hazard level well (by matching conditional spectra) and fully describe the complete set of candidate ground motions. Three case studies for modified, pulse-type, and non-pulse-type ground motions are designed to evaluate the performance of the proposed ground motion clustering algorithm (convolutional autoencoder + k-means). Considering the limited number of pre-selected candidate ground motions in the last two case studies, the response spectra simulation and transfer learning are used to improve the stability and reproducibility of the proposed ground motion clustering algorithm. The results of the three case studies demonstrate that the convolutional autoencoder + k-means can 1) achieve 100% accuracy in classifying ground motion response spectra, 2) correctly determine the optimal number of clusters, and 3) outperform established clustering algorithms (i.e., autoencoder + k-means, time series k-means, spectral clustering, and k-means on ground motion influence factors). Using the proposed clustering-based ground motion selection method, an application is performed to select ground motions for a structure in San Francisco, California. The developed user-friendly codes are published for practical use. 
    more » « less