Geology topics

Jonathan P. Stewart

Publications and source records attributed to Jonathan P. Stewart.

At least 19 recordsLinked to original sources

USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗

USGS2024

Data-driven adjustments for combined use of NGA-East hard-rock ground motion and site amplification models

Model development in the Next Generation Attenuation-East (NGA-East) project included two components developed concurrently and independently: (1) earthquake ground-motion models (GMMs) that predict the median and aleatory variability of various intensity measures conditioned on magnitude and distance, derived for a reference hard-rock site condition with an average shear-wave velocity in the upper 30 m ( V S30 ) = 3000 m/s; and (2) a site amplification model that modifies intensity measures for softer site conditions. We investigate whether these models, when used in tandem, are compatible with ground-motion recordings in central and eastern North America (CENA) using an expanded version of the NGA-East database that includes new events from November 2011 (end date of NGA-East data curation) to April 2022. Following this expansion, the data set has 187 events, 2096 sites, and 16,272 three-component recordings, although the magnitude range remains limited (∼4 to 5.8). We compute residuals using 17 NGA-East GMMs and three data selection criteria that reflect within-CENA regional variations in ground-motion attributes. Mixed-effects regression of the residuals reveals a persistent pattern in which ground motions are overpredicted at short periods (0.01–0.6 s, including peak ground acceleration (PGA)) and underpredicted at longer periods. These misfits are regionally variable, with the Texas–Oklahoma–Kansas region having larger absolute misfits than other parts of CENA. Two factors potentially influencing these misfits are (1) differences in the site amplification models used to adjust the data to the reference condition during NGA-East GMM development relative to CENA amplification models applied since the 2018 National Seismic Hazard Model (NSHM), and (2) potential bias in simulation-based factors used to adjust ground motions from the hard-rock reference condition to a V S30 = 760 m/s condition. We provide adjustment factors and their epistemic uncertainties and discuss implications for applications.

Kansas, Oklahoma, Texas↗

USGS2023

Findings from a decade of ground motion simulation validation research and a path forward

Simulated ground motions have the potential to advance seismic hazard assessments and structural response analyses, particularly for conditions with limited recorded ground motions such as large magnitude earthquakes at short source-to-site distances. However, rigorous validation of simulated ground motions is needed for hazard analysts, practicing engineers, or regulatory bodies to be confident in their use. A decade ago, validation exercises were mainly limited to comparisons of simulated-to-observed waveforms and median values of spectral accelerations for selected earthquakes. The Southern California Earthquake Center (SCEC) Ground Motion Simulation Validation (GMSV) group was formed to increase coordination between simulation modelers and research engineers with the aim of devising and applying more effective methods for simulation validation. Here, we summarize what has been learned in over a decade of GMSV activities, principally reflecting the views of the SCEC research community but also extending our findings and suggestions for a path forward to broader United States and worldwide simulation validation efforts. We categorize different validation methods according to their approach and the metrics considered. Two general approaches are to compare validation metrics from simulations to those from historical records or to those from semi-empirical models. Validation metrics are categorized into ground motion characteristics and structural responses. We discuss example validation studies that have been impactful in the past decade and suggest future research directions. Key lessons learned are that validation is application-specific, our outreach and dissemination need improvement, and much validation-related research remains unexplored.

Earthquake Spectra↗

USGS2023

High-pass corner frequency selection for implementation in the USGS automated ground motion processing tool

Earthquake ground motion processing for next-generation attenuation (NGA) projects required human inspection to select high-pass corner frequencies ( f cHP ), which is time-intensive and subjective. With growth in the number of recordings per event and interest in enhancing repeatability, we sought to develop automated procedures for f cHP selection. These procedures consider signal-to-noise ratio (SNR) and non-physical features in the displacement time series that indicate high- and/or low-frequency noise effects. The procedures are implemented in a US Geological Survey software package (gmprocess). We extend previous procedures for SNR-based corner frequency selection to also check for low-frequency artifacts in the displacement record using a polynomial fit to improve f cHP selection. We evaluate the performance of the SNR and polynomial fit criteria using recordings from the 2020 M5.1 North Carolina and the 2013 M4.7 Southern Ontario earthquakes. Data processed with the SNR-only criteria can have low f cHP and displacement drift; the displacement check increases f cHP and reduces displacement drift.

Conference Paper↗

USGS2023

Panel review of Ground Motion Characterization Model in 2023 NSHM

The 2023 National Seismic Hazard Model (NSHM; Petersen et al., 2023) has two major components – a seismic source characterization (SSC) model and a ground motion characterization (GMC) model. The US Geological Survey (USGS) established separate panels to review and provide input on these two models. Both panels are advisory, meaning that they provide input on technical issues for consideration by the USGS NSHM team, but they do not have decision making authority. Here, we report on the activities and recommendations of the Ground Motion Characterization Panel, made up of the authors of this review. Final modeling decisions are presented in separate USGS documents, including Petersen et al., (2023) and Moschetti et al., (2023). Where modeling decisions depart from our recommendations, the rationale is explained in those publications.

Report↗

USGS2022

Comparisons of the NGA-Subduction ground motion models

In this article, ground-motion models (GMMs) for subduction earthquakes recently developed as part of the Next Generation Attenuation-Subduction (NGA-Sub) project are compared. The four models presented in this comparison study are documented in their respective articles submitted along with this article. Each of these four models is based on the analysis of the large NGA-Sub database. Three of the four current models are developed for a global version as well as separate regionalized models. The fourth model was developed based on earthquakes only from Japan, and as such is applicable only for Japan. As part of this comparison study, a general discussion on the parameterization of the four models and the regionalization of the three models is provided. The specific strengths and or weaknesses or the technical decisions and justifications of any one model are not part of this comparison. A selected suite of deterministic attenuation curves and spectra are presented for the models along with a selected suite of currently used subduction models. A limited number of comparisons are presented in this article with a larger number of comparisons and the digital values provided in the electronic attachment. In addition to these scenario calculation comparisons, the results from a standard probabilistic seismic hazard analysis (PSHA) for two sites located in the Pacific Northwest Region in the state of Washington are presented. These calculations highlight the potential impact of using the new GMMs. Based on the comparisons presented here, a general understanding of these new GMMs can be obtained with the expectation that the implementation of a specific seismic hazard study should incorporate similar and additional comparisons and sensitivity studies pertinent to the site of interest.

Earthquake Spectra↗

USGS2022

Validating predicted site response in sedimentary basins from 3D ground motion simulations

We introduce procedures to validate site response in sedimentary basins as predicted using ground motion simulations. These procedures aim to isolate contributions of site response to computed intensity measures relative to those from seismic source and path effects. In one of the validation procedures, simulated motions are analyzed in the same manner as earthquake recordings to derive non-ergodic site terms. This procedure compares the scaling with sediment isosurface depth of simulated versus empirical site terms (the latter having been derived in a separate study). A second validation procedure utilizes two sets of simulations, one that considers three-dimensional (3D) basin structure and a second that utilizes a one-dimensional (1D) representation of the crustal structure. Identical sources are used in both procedures, and after correcting for variable path effects, differences in ground motions are used to estimate site amplification in 3D basins. Such site responses are compared to those derived empirically to validate both the absolute levels and the depth scaling of site response from 3D simulations. We apply both procedures to southern California in a manner that is consistent between the simulated and empirical data (i.e. by using similar event locations and magnitudes). The results show that the 3D simulations overpredict the depth-scaling and absolute levels of site amplification in basins. However, overall patterns of site amplification with depth are similar, suggesting that future calibration may be able to remove observed biases.

Earthquake Spectra↗

USGS2022

Performance of NGA-East GMMs and site amplification models relative to CENA ground motions

We investigate bias in ground motions predicted for Central and Eastern North America (CENA) using ground motion models (GMMs) combined with site amplification models developed in the NGA-East project. Bias is anticipated because of de-coupled procedures used in the development of the GMMs and site amplification models. The NGA-East GMMs were mainly calibrated by adjusting CENA data to a reference site condition using a site amplification model appropriate for active tectonic regions. Hence, these GMMs are likely biased relative to the CENA reference site condition (3000 m/sec shear wave velocity). Moreover, the NGA-East site amplification model recommended for hazard applications contains a simulation-based term for amplification between the reference condition and time-averaged shear wave velocity V S30 =760 m/sec, which is uncertain and has not been calibrated against data from sites with that reference condition. Using the NGA-East dataset, we apply mixed-effects residual analysis and identify that period-dependent bias in 5% damped response spectral acceleration is present across a wide frequency range, but is strongest (i.e., overestimating by a factor of 2) at short oscillator periods <0.2 sec. Ongoing work to remedy this bias consists of expanding the NGA-East dataset with more recent recordings and enhanced metadata, particularly regarding site conditions.

Conference Paper↗

USGS2022

Data resources for NGA-subduction project

A relational database was developed over a five-year period to support ground motion model (GMM) development for the Next Generation Attenuation-Subduction (NGA-Sub) project. The relational database has components that interact according to a database schema, including a source and path component used to describe attributes of seismic sources in global subduction regions and to compute source-to-site distances, a site component that describes attributes of sites where recordings have been made, and a ground motion component. The source component of the database has information for 1880 earthquakes, mainly from the following regions: the Pacific Northwest region of North America, Alaska and the Aleutian Islands, Japan, Taiwan, New Zealand, South America, Central America, and Mexico. Of the 1880 earthquakes, 88 have finite fault models (FFMs) from the literature that were systematically reviewed, distilled to one more rectangular shapes, and trimmed according to procedures based on percentage of total slip. For earthquakes without FFMs, a simulation routine is used to represent finite fault effects required for distance calculations. This simulation routine was adjusted and made more uniform in its application than in prior NGA projects. All earthquakes are classified as interface, intraslab, shallow crustal, or outer rise, using uniform protocols developed for this project. All earthquakes are also assigned class designations adapted from a prior NGA project for active regions, that allows foreshock, mainshock, and aftershock events to be distinguished. The site component of the database is described in a companion paper (Ahdi et al. 2020 [1]). The ground motion component of the database consists of median – and maximum – horizontal component peak parameters (peak ground acceleration, PGA and peak ground velocity, PGV) and pseudo-spectral accelerations (PSa) at 111 oscillator periods and 11 damping ratios. Response spectra were also computed for the vertical component. Fourier amplitude spectra (FAS) and duration metrics were also computed. The ground motion recordings were obtained from collaborating organizations world-wide as uncorrected (Vol 1) digital recordings, that were corrected (componentspecific low – and high – pass filters and baseline correction, as needed) following Pacific Earthquake Engineering Research Center (PEER)/NGA protocols. The relational database operates on each of these (and other) database components to dynamically draw relevant parameters into a single file, known as a flatfile, that is used by researchers engaged in GMM development. The flatfiles used in model development are being published with the NGA-Sub GMMs as products of the NGA-Sub project.

Extramural-Authored Publication Paper↗

USGS2021

Relational database for horizontal‐to‐vertical spectral ratios

Frequency‐dependent horizontal‐to‐vertical spectral ratios (HVSRs) of Fourier amplitudes from three‐component recordings can provide useful information for site response modeling. However, such information is not incorporated into most ground‐motion models, including those from Next‐Generation Attenuation projects, which instead use the time‐averaged shear‐wave velocity ( ⁠ V S "> V S VS ⁠ ) in the upper 30 m of the site and sediment depth terms. To facilitate utilization of HVSR, we developed a publicly accessible relational database. This database is adapted from a similar repository for V S "> V S VS data and provides microtremor‐based HVSR data (mHVSR) and supporting metadata, but not parameters derived from the data. Users can interact with the data directly within a web portal that contains a graphical user interface (GUI) or through external tools that perform cloud‐based computations. Within the database GUI, the median horizontal‐component mHVSR can be plotted against frequency, with the mean and mean ± one standard deviation (representing variability across time windows) provided. Using external interactive tools (provided as a Jupyter Notebook and an R script), users can replot mHVSR (as in the database) or create polar plots. These tools can also derive parameters of potential interest for modeling purposes, including a binary variable indicating whether an mHVSR plot contains peaks, as well as the fitted properties of those peaks (frequencies, amplitudes, and widths). Metadata are also accessible, which includes site location, details about the instruments used to make the measurements, and data processing information related to windowing, antitrigger routines, and filtering.

Seismological Research Letters↗

USGS2021

NGA-Subduction research program

This article summarizes the Next Generation Attenuation (NGA) Subduction (NGA-Sub) project, a major research program to develop a database and ground motion models (GMMs) for subduction regions. A comprehensive database of subduction earthquakes recorded worldwide was developed. The database includes a total of 214,020 individual records from 1,880 subduction events, which is by far the largest database of all the NGA programs. As part of the NGA-Sub program, four GMMs were developed. Three of them are global subduction GMMs with adjustment factors for up to seven worldwide regions: Alaska, Cascadia, Central America and Mexico, Japan, New Zealand, South America, and Taiwan. The fourth GMM is a new Japan-specific model. The GMMs provide median predictions, and the associated aleatory variability, of RotD50 horizontal components of peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral acceleration (PSA) at oscillator periods ranging from 0.01 to 10 s. Three GMMs also quantified “within-model” epistemic uncertainty of the median prediction, which is important in regions with sparse ground motion data, such as Cascadia. In addition, a damping scaling model was developed to scale the predicted 5%-damped PSA of horizontal components to other damping ratios ranging from 0.5% to 30%. The NGA-Sub flatfile, which was used for the development of the NGA-Sub GMMs, and the NGA-Sub GMMs coded on various software platforms, have been posted for public use.

Earthquake Spectra↗

USGS2021

Ergodic site response model for subduction zone regions

We present an ergodic site response model with regional adjustments for use with subduction zone ground-motion models. The model predicts site amplification of peak ground acceleration, peak ground velocity, and 5% damped pseudo-spectral accelerations of the orientation-independent horizonal component for oscillator periods from 0.01 to 10 s. The model depends on the time-averaged shear-wave velocity in the upper 30 m ( V S 30 ), basin depth, and region and is independent of subduction earthquake type. It has three components: a linear site-amplification term in the form of V S 30 -scaling, a nonlinear term that depends on V S 30 and shaking intensity parameterized by peak ground acceleration at the reference-rock velocity condition of 760 m/s, and a basin sediment-depth term for Japan and Cascadia conditioned on the depth to the 2.5 km/s shear-wave velocity isosurface ( Z 2.5 ). A global V S 30 -scaling model is provided along with regional adjustments for Japan, Taiwan, South America, Alaska, and Cascadia. The nonlinear model is global, with a functional form that has often been used to fit nonlinear responses inferred from simulations, but here we calibrate it empirically. Relative to a prior model for shallow earthquakes in active tectonic regions, our subduction zone global V S 30 -scaling is comparable at short periods (<1.0 s) but weaker at long periods, while the nonlinear site response is generally less pronounced but extends to lower levels of shaking. Basin depth models are conditioned on the difference of the actual Z 2.5 and a V S 30 -conditioned mean Z 2.5 . Sites with positive differential depths have increased long-period site responses and decreased short-period responses, with the opposite occurring for negative differential depths.

Earthquake Spectra↗

USGS2021

NGA-subduction global ground motion models with regional adjustment factors

We develop semi-empirical ground motion models (GMMs) for peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral accelerations for periods from 0.01 to 10 s, for the median orientation-independent horizontal component of subduction earthquake ground motion. The GMMs are applicable to interface and intraslab subduction earthquakes in Japan, Taiwan, Mexico, Central America, South America, Alaska, the Aleutian Islands, and Cascadia. The GMMs are developed using a combination of data inspection, data regression with respect to physics-informed functions, ground-motion simulations, and geometrical constraints for certain model components. The GMMs capture observed differences in source and path effects for interface and intraslab events, conditioned on moment magnitude, rupture distance, and hypocentral depth. Site effect and aleatory variability models are shared between event types. Regionalized GMM components include the model constant (that controls ground motion amplitude), anelastic attenuation, magnitude-scaling break point, linear site response, and sediment depth terms. We develop models for the aleatory between-event variability ( τ ) "> ( τ ) , within-event variability ( ϕ ) "> ( ϕ ) , single-station within-event variability ( ϕ SS ) "> ( ϕ SS ) , and site-to-site variability ( ϕ S 2 S ) "> ( ϕ S 2 S ) . Ergodic analyses should use the median GMM and aleatory variability computed using the between-event and within-event variability models. An analysis incorporating non-ergodic site response should use the median GMM at the reference shear-wave velocity condition, a site-specific site response model, and aleatory variability computed using the between-event and single-station within-event variability models. Epistemic uncertainty in the median model is represented by standard deviations on the regional model constants, which facilitates scaled-backbone representations of model uncertainty in hazard analyses.

Earthquake Spectra↗

USGS2021

Horizontal-to-vertical spectral ratios from California sites: Open-source database and data interpretation to establish site parameters

Frequency-dependent horizontal-to-vertical spectral ratios (HVSR) of Fourier amplitudes from three-component recordings can provide information on one or more site resonant frequencies and relative levels of amplification at those frequencies. Such information is potentially useful for predicting site amplification but is not present in site databases that have been developed over the last 15–20 years for the Next-Generation Attenuation (NGA) projects, which instead use the time-averaged shear-wave velocity (VS) in the upper 30 m of the site (VS30) as the primary site parameter and are supplemented with basin depth terms where available. As a consequence, HVSR parameters are also not used in NGA ground motion models. In order for HVSR-based parameters to be used in future versions of site databases, a publicly accessible repository of this information is needed. We adapt a relational database developed to archive and disseminate VS data to also include HVSR. The database provides relevant microtremor-based HVSR data (mHVSR) and supporting metadata. We consider the most relevant data to be the frequency-dependent mHVSR, where the horizontal is taken as the median component and also as a function of horizontal azimuth (referred to as polar plots). Relevant metadata includes site location information, details about the equipment used to make the measurements, and processing details related to windowing, anti-trigger routines, and filtering. We describe the database schema developed to organize and present this information.

California↗

USGS2020

A ground‐motion prediction model for shallow crustal earthquakes in Greece

Using a recently completed database of uniformly processed strong‐motion data recorded in Greece, we derive a ground‐motion prediction model (GMPM) for horizontal‐component peak ground velocity, peak ground acceleration, and 5% damped pseudoacceleration response spectra, at 105 periods ranging from 0.01 to 10 s. The equations were developed by modifying a global GMPM, to account for more rapid attenuation and weaker magnitude scaling in the Greek ground motions than in the global GMPM. Our GMPM is calibrated using the Greek data for distances up to 300 km, magnitudes from 4.0 to 7.0, and time‐averaged 30 m shear‐wave velocities from 150 to 1200    m / s "> 1200 m / s 1200 m/s ⁠ . The GMPM has important attributes for hazard applications including magnitude scaling that extends the range of applicability to M "> M M 8.0 and nonlinear site response. These features are possible because they are well constrained by data in the global GMPM from which our model is derived. An interesting feature of the Greek data, also observed previously in studies of mid‐magnitude events (6.1–6.5) in Italy, is that they are substantially overpredicted by the global GMPM, which may be a repeatable regional feature, but may also be influenced by soil–structure interaction. This bias is an important source of epistemic uncertainty that should be considered in hazard analysis.

Bulletin of the Seismological Society of America↗

USGS2020

Liquefaction and related ground failure from July 2019 Ridgecrest earthquake sequence

The 2019 Ridgecrest earthquake sequence produced a 4 July M 6.5 foreshock and a 5 July M 7.1 mainshock, along with 23 events with magnitudes greater than 4.5 in the 24 hr period following the mainshock. The epicenters of the two principal events were located in the Indian Wells Valley, northwest of Searles Valley near the towns of Ridgecrest, Trona, and Argus. We describe observed liquefaction manifestations including sand boils, fissures, and lateral spreading features, as well as proximate non‐ground failure zones that resulted from the sequence. Expanding upon results initially presented in a report of the Geotechnical Extreme Events Reconnaissance Association, we synthesize results of field mapping, aerial imagery, and inferences of ground deformations from Synthetic Aperture Radar‐based damage proxy maps (DPMs). We document incidents of liquefaction, settlement, and lateral spreading in the Naval Air Weapons Station China Lake US military base and compare locations of these observations to pre‐ and postevent mapping of liquefaction hazards. We describe liquefaction and ground‐failure features in Trona and Argus, which produced lateral deformations and impacts on several single‐story masonry and wood frame buildings. Detailed maps showing zones with and without ground failure are provided for these towns, along with mapped ground deformations along transects. Finally, we describe incidents of massive liquefaction with related ground failures and proximate areas of similar geologic origin without ground failure in the Searles Lakebed. Observations in this region are consistent with surface change predicted by the DPM. In the same region, geospatial liquefaction hazard maps are effective at identifying broad percentages of land with liquefaction‐related damage. We anticipate that data presented in this article will be useful for future liquefaction susceptibility, triggering, and consequence studies being undertaken as part of the Next Generation Liquefaction project.

California↗

USGS2020

Risk-targeted alternatives to deterministic ground motion caps in U.S. seismic provisions

Since their inception over 20 years ago, the maximum considered earthquake ground motion maps in U.S. building codes have capped probabilistic values with deterministic ground motions from characteristic earthquakes on known active faults. This practice has increasingly been called into question both because of spatially non-uniform risk levels that are produced (risk being higher mainly in coastal California) and practical difficulties in defining characteristic earthquakes from recent earthquake rupture forecast models. We describe two proposals developed to enable phase-out of deterministic caps. One approach modestly increases collapse risk targets nationwide based on recent information on return periods of characteristic earthquakes on major central and eastern U.S. seismic sources; adoption of this approach would remove the perceived need for caps in California. The second approach uses geographically varying collapse risk targets, being higher near the highly active faults in California and unchanged elsewhere. Neither approach was adopted for the 2020 National Earthquake Hazards Reduction Program recommended seismic Provisions for new building structures, but they are described in a Part 3 document to accompany the Provisions and Commentary .

Earthquake Spectra↗

USGS2018

Development of a United States community shear wave velocity profile database

We present the details of a multi-institutional effort to develop an open-access shear-wave velocity ( V S ) profile database (PDB), which will include a public repository for V S profile data in the United States. V S profiles are an essential resource for ground motion modeling and other applications. The minimum requirements for a site to be included in the database are in situ geophysical V S measurements and location metadata (geodetic coordinates and elevation). Other information is included as available, including geotechnical logs, penetration resistance, laboratory test data, ground water elevation, and P-wave velocity profiles. The project is currently at the stage of data collection (over 4500 V S profiles) and prototype data model development. The database will be presented as an online map-based interface with downloadable V S profile and metadata information. This paper serves as a progress report to the geotechnical and earthquake engineering communities, as we seek community engagement and support.

Conference Paper↗