Search USGSSearch

USGS · 70276582

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

Abstract

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 17.00324098022206° to 19.302954943559428° latitude; -68.32307657947791° to -64.67497609311363° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

BibTeXRIS

Allison M. Shumway, Kevin Ross Milner, Peter M. Powers, Morgan P. Moschetti, Jason M. Altekruse, Julie A. Herrick, Andrea L. Llenos, Kyle B. Withers, Edward H. Field, Brad T. Aagaard, Richard W. Briggs, Alexandra Elise Hatem, Kirstie Lafon Haynie, Andrew J. Michael, Jessica A. Thompson Jobe, Kishor S. Jaiswal, Brandon S. Clayton, Nico Luco, Mark D. Petersen, Sanaz Rezaeian, Thomas L. Pratt, Yuehua Zeng. 2026-04-28. 2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results. https://doi.org/10.1002/esp4.70055

Cite the original work for its findings. Save a collection to share your selection of sources.

Discover connections

Connections use source metadata and explicit phrase matches, not verified experimental comparisons.

KEEP EXPLORING

Related USGS reports

Earthquake stress drop, source spectra from laboratory experiments, and the apparent unsteadiness of frictional slip

Over a range of co-seismic total energy release, high sample rate slip, stress, and slip velocity measurements were made during frictional failure to determine seismic source properties. The near-field measurements were used to estimate stress drops using spectral approaches that are analogous to the standard seismological moment-corner determined Δσ and acceleration amplitude Δσ arms stress drops. Rapid frictional sliding produces slip velocity amplitude spectra that vary with frequency -1 so long as the energy dissipated by friction approaches the total energy release. Under these conditions Δσ and Δσ arms are similar in size. Since these are nearly direct measurements of fault slip velocity and because they produce flat acceleration spectra between the corner and an apparent fmax, friction in these experiments seems to produce dynamic acceleration and deceleration at all resolved frequencies. The physical origin of this unsteady motion is not known exactly but likely reflects rapid dynamic variations in fault strength or applied stress. On the other hand, when the total energy stored prior to failure is artificially made somewhat larger relative to that which can be dissipated by friction the velocity spectra vary with frequency -2 immediately below the corner. At these conditions Δσ and Δσ arms diverge. The excess stored energy is preferentially partitioned into radiated energy upon failure. Collectively the experiments imply that to produce unsteady, white noise accelerations that are observed for natural earthquakes, requires that at least 95% of the energy released does not escape the source region to be radiated to the far-field.

Bulletin of the Seismological Society of America

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş

A partially nonergodic ground-motion model for Fourier amplitude spectra for the San Francisco Bay area, California, USA

We develop a partially nonergodic ground-motion model (GMM) for Fourier amplitude spectra for the San Francisco Bay Area, California, USA, using the Bayless and Abrahamson (2019) GMM as a reference ergodic GMM and developing location-dependent adjustments to the predicted median and variance. We compile regional ground-motion data from moment magnitude (𝑀 w ) >3 earthquakes occurring during 2000–2022 for which magnitude information is available in the U.S. Geological Survey Comprehensive Catalog (Guy et al., 2015). The data set predominantly consists of records from 𝑀 w 3.5–4.5 earthquakes but includes three well-recorded 𝑀 w > 5 events. Ground-motion residuals are evaluated using the time-averaged shear-wave velocity in the top 30 m (𝑉 S30 ) from the California-specific map of Thompson et al. (2018) and basin-depth site parameters from the seismic velocity model of Aagaard and Hirakawa (2021). The 𝑉 S30 dependence and basin-depth scaling of the reference ergodic GMM of Bayless and Abrahamson (2019) are evaluated and modified with the updated data set. We compute maps of site adjustments using a varying-coefficient model that considers the spatial correlation structure and uncertainties at each observation location. The spatial covariance model is developed using ground-motion residuals that are standardized by the uncertainty model, which allows for consideration of the aleatory variability in developing the site adjustments. The covariance model is fit considering the means and standard deviations of the site terms at all locations. The use of partially nonergodic median adjustments results in modified variance components of the within-event variability. Due to the low number of large-magnitude earthquakes that control seismic hazard in the data set, we do not modify between-event variance; however, we present adjustments to site-to-site variability for use in partially nonergodic hazard assessments.

California