Search USGS⌕ Search

Geology topics

Kyle B. Withers

Publications and source records attributed to Kyle B. Withers.

17 recordsLinked to original sources

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

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.

Puerto Rico, U.S. Virgin Islands↗

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands↗

The 2023 U.S. National Seismic Hazard Model: Subduction ground motion models

The US Geological Survey National Seismic Hazard Models (NSHMs) are used to calculate earthquake ground-shaking intensities for design and rehabilitation of structures in the United States. The most recent 2014 and 2018 versions of the NSHM for the conterminous United States included major updates to ground-motion models (GMMs) for active and stable crustal tectonic settings; however, the subduction zone GMMs were largely unchanged. With the recent development of the next generation attenuation-subduction (NGA-Sub) GMMs, and recent progress in the utilization of “M9” Cascadia earthquake simulations, we now have access to improved models of ground shaking in the US subduction zones and the Seattle basin. The new NGA-Sub GMMs support multi-period response spectra calculations. They provide global models and regional terms specific to Cascadia and terms that account for deep-basin effects. This article focuses on the updates to subduction GMMs for implementation in the 2023 NSHM and compares them to the GMMs of previous NSHMs. Individual subduction GMMs, their weighted averages, and their impact on the estimated mean hazard relative to the 2018 NSHM are discussed. The updated logic trees include three of the new NGA-Sub GMMs and retain two older models to represent epistemic uncertainty in both the median and standard deviation of ground-shaking intensities at all periods of interest. Epistemic uncertainty is further represented by a three-point logic tree for the NGA-Sub median models. Finally, in the Seattle region, basin amplification factors are adjusted at long periods based on the state-of-the-art M9 Cascadia earthquake simulations. The new models increase the estimated mean hazard values at short periods and short source-to-site distances for interface earthquakes, but decrease them otherwise, relative to the 2018 NSHM. On softer soils, the new models cause decreases to the estimated mean hazard for long periods in the Puget Lowlands basin but increases within the deep Seattle portion of this basin for short periods relative to the 2018 NSHM.

Earthquake Spectra↗

Basin effects from 3D simulated ground motions in the Greater Los Angeles region for use in seismic-hazard analyses

We develop basin-depth-scaling models (i.e. “basin terms”) from the long-period ( ⁠ T ≥2s⁠ ) simulated ground motions of the Southern California Earthquake Center (SCEC) CyberShake project for use in seismic hazard analyses at sites within the sedimentary basins of southern California. Basin terms use the Next Generation Attenuation (NGA)-West-2 ground-motion models (GMMs) as reference models and use their functional forms with slight modifications. We investigate the use of two approaches to incorporate the time-averaged shear-wave velocity in the upper 30 m ( ⁠ V S30 ⁠ ) in these calculations and find that the use of site-specific and uniform V S30 has minor effects on the resulting basin terms for this data set. By centering the simulated ground motions on the basin terms, we separate the information from the simulations about absolute ground-motion level from information relating to the relative amplifications, such as the differences between shallow- and deep-basin sites. Recent observations from sedimentary basins of southern California indicate that additional amplification effect may persist at relatively shallow basin depths (i.e. the GMM basin terms should have positive values when differential depths, δ Z1⁠ , are near zero), and we present models for “centered” and “adjusted” basin-depth scaling models that reflect this potential. The simulation-modified GMMs are appropriate for crustal sources and for deep-basin sites ( ⁠δ Z1>0⁠ ) within basins of the Greater Los Angeles region, for the magnitudes and distances defined by each of the reference NGA-West-2 GMMs.

California↗

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra↗

The 2023 US National Seismic Hazard Model: Ground-motion characterization for the conterminous United States

We update the ground-motion characterization for the 2023 National Seismic Hazard Model (NSHM) for the conterminous United States. The update includes the use of new ground-motion models (GMMs) in the Cascadia subduction zone; an adjustment to the central and eastern United States (CEUS) GMMs to reduce misfits with observed data; an updated boundary for the application of GMMs for shallow, crustal earthquakes in active tectonic regions (i.e. western United States (WUS)) and stable continental regions (i.e. CEUS); and the use of improved models for the site response of deep sedimentary basins in the WUS and CEUS. Site response updates include basin models for the California Great Valley and for the Portland and Tualatin basins, Oregon, as well as long-period basin effects from three-dimensional simulations in the Greater Los Angeles region and in the Seattle basin; in the CEUS, we introduce a broadband (0.01- to 10-s period) amplification model for the effects of the passive-margin basins of the Atlantic and Gulf Coastal Plains. In addition, we summarize progress on implementing rupture directivity models into seismic hazard models, although they are not incorporated in the 2023 NSHM. We implement the ground-motion characterization for the 2023 NSHM in the US Geological Survey’s code for probabilistic seismic hazard analysis, nshmp-haz-v2 , and present the sensitivity of hazard to these changes. Hazard calculations indicate widespread effects from adjustments to the CEUS GMMs, from the incorporation of Coastal Plain amplification effects, and from the treatment of shallow-basin and out-of-basin sites in the San Francisco Bay Area and Los Angeles region, as well as locally important changes from subduction-zone GMMs, and from updated and new WUS basins.

conterminous United States↗

Framework for implementing damping scaling factors in U.S. Geological Survey National Seismic Hazard Models

Traditionally, probabilistic seismic hazard analysis (PSHA) has focused on calculating ground motion hazard curves for elastic, 5%-damped pseudo spectral accelerations, Sa(T,5%), which are used as the basis for engineering design parameters and targets for ground motion selection and modification. However, structures and geotechnical systems can exhibit a wide range of damping ratios both above and below the 5% level, depending on the construction material, structural system, nonstructural elements, or subsurface soil properties. When spectral parameters at such damping levels are required for certain applications, 5%-damped accelerations have traditionally been extracted from PSHA-based hazard curves and adjusted outside of the hazard integral using damping scaling factors (DSF) such as those from Newmark & Hall (1982). Recent advances in the development of more rigorous and comprehensive damping scaling models (e.g., Rezaeian et al., 2014; Rezaeian et al., 2021) have allowed for the modeling of means and standard deviations of DSFs as functions of earthquake source and path properties for crustal, intraslab, and subduction interface tectonic environments. These DSF models can be applied to ground motion model (GMM) estimates of Sa(T,5%) for a given earthquake rupture scenario to produce a corresponding mean and standard deviation Sa at a specified damping ratio β, Sa(T,β). In this study, the DSF models of Rezaeian et al. (2014) and Rezaeian et al. (2021) are implemented within the U.S. Geological Survey National Seismic Hazard Model (NSHM) PSHA framework to calculate probabilistic hazard curves for spectral accelerations at damping ratios from 0.5% to 30%. The DSF models are applied directly to the mean and standard deviation of Sa(T,5%) predictions from each GMM in the NSHM logic tree. Resulting hazard curves and uniform hazard and risk spectra for Sa(T,β) are presented for several geographic locations and compared with corresponding spectra estimated using current design practices by applying the same DSFs outside of the PSHA calculation. Key differences between the two methods for estimating Sa(T,β) are discussed, and potential strategies are presented for the implementation and usage of the hazard-consistent Sa(T,β) in building codes. Comparing the results to those from DSFs used in current design practices that are mainly based on Newmark & Hall (1982) is not explored in this study.

Conference Paper↗

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗

Automated detection of clipping in broadband earthquake records

Because the amount of available ground‐motion data has increased over the last decades, the need for automated processing algorithms has also increased. One difficulty with automated processing is to screen clipped records. Clipping occurs when the ground‐motion amplitude exceeds the dynamic range of the linear response of the instrument. Clipped records in which the amplitude exceeds the dynamic range are relatively easy to identify visually yet challenging for automated algorithms. In this article, we seek to identify a reliable and fully automated clipping detection algorithm tailored to near‐real‐time earthquake response needs. We consider multiple alternative algorithms, including (1) an algorithm based on the percentage difference in adjacent data points, (2) the standard deviation of the data within a moving window, (3) the shape of the histogram of the recorded amplitudes, (4) the second derivative of the data, and (5) the amplitude of the data. To quantitatively compare these algorithms, we construct development and holdout datasets from earthquakes across a range of geographic regions, tectonic environments, and instrument types. We manually classify each record for the presence of clipping and use the classified records. We then develop an artificial neural network model that combines all the individual algorithms. Testing on the holdout dataset, the standard deviation and histogram approaches are the most accurate individual algorithms, with an overall accuracy of about 93%. The combined artificial neural network method yields an overall accuracy of 95%, and the choice of classification threshold can balance precision and recall.

Seismological Research Letters↗

Spectral damping scaling factors for horizontal components of ground motions from subduction earthquakes using NGA-Subduction data

This article develops global models of damping scaling factors (DSFs) for subduction zone earthquakes that are functions of the damping ratio, spectral period, earthquake magnitude, and distance. The Next Generation Attenuation for subduction earthquakes (NGA-Sub) project has developed the largest uniformly processed database of recorded ground motions to date from seven subduction regions: Alaska, Cascadia, Central America and Mexico, South America, Japan, Taiwan, and New Zealand. NGA-Sub used this database to develop new ground motion models (GMMs) at a reference 5% damping ratio. We worked with the NGA-Sub project team to develop an extended database that includes pseudo-spectral accelerations (PSA) for 11 damping ratios between 0.5% and 30%. We use this database to develop parametric models of DSF for both interface and intraslab subduction earthquakes that can be used to adjust any subduction GMM from a reference 5% damping ratio to other damping ratios. The DSF is strongly influenced by the response spectral shape and the duration of motion; therefore, in addition to the damping ratio, the median DSF model uses spectral period, magnitude, and distance as surrogate predictor variables to capture the effects of the spectral shape and the duration of motion. We also develop parametric models for the standard deviation of DSF. The models presented in this article are for the RotD50 horizontal component of PSA and are compared with the models for shallow crustal earthquakes in active tectonic regions. Some noticeable differences arise from the considerably longer duration of interface records for very large magnitude events and the enriched high-frequency content of intraslab records, compared with shallow crustal earthquakes. Regional differences are discussed by comparing the proposed global models with the data from each subduction region along with recommendations on the applicability of the models.

Earthquake Spectra↗

A machine learning approach to developing ground motion models from simulated ground motions

We use a machine learning approach to build a ground motion model (GMM) from a synthetic database of ground motions extracted from the Southern California CyberShake study. An artificial neural network is used to find the optimal weights that best fit the target data (without overfitting), with input parameters chosen to match that of state-of-the-art GMMs. We validate our synthetic-based GMM with empirically based GMMs derived from the globally based Next Generation Attenuation West2 data set, finding near-zero median residuals and similar amplitude and trends (with period) of total variability. Additionally, we find that the artificial neural network GMM has similar bias and variability to empirical GMMs from records of the recent Ridgecrest event, which neither GMM has included in its formulation. As simulations continue to better model broadband ground motions, machine learning provides a way to utilize the vast amount of synthetically generated data and guide future parameterization of GMMs.

California↗

Combining dynamic rupture simulations with ground motion data to characterize seismic hazard from Mw 3-5.8 earthquakes in Oklahoma and Kansas

Many seismically active areas suffer from a lack of near‐source ground‐motion recordings, making ground‐motion prediction difficult at distances within &#x223C; 40 &#x2009;&#x2009; km "> ∼ 40 km ∼40 km from an earthquake. We aim to aid the development of near‐source ground‐motion prediction equations (GMPEs) by generating synthetic ground‐motion data via simulation. Building on previous work using point‐source moment tensor sources to simulate small ( ⁠ M w "> M w Mw 3–4) earthquakes for a target region encompassing north central Oklahoma and south central Kansas, we perform dynamic rupture simulations of earthquakes up to M w "> M w Mw 5.8. We introduce complexity into the rupture process by adding stochastically generated heterogeneity to initial stress conditions. Our simulated ground‐motion data are added to a catalog of recorded ground‐motion data to construct a composite recorded‐simulated ground‐motion catalog that we use to develop a GMPE for the target region. This procedure can be generalized and used for GMPE development in other regions with near‐source ground‐motion data scarcity, which could directly benefit critical applications such as the National Seismic Hazard Maps produced by the U.S. Geological Survey (USGS).

Kansas, Oklahoma↗

A suite of exercises for verifying dynamic earthquake rupture codes

We describe a set of benchmark exercises that are designed to test if computer codes that simulate dynamic earthquake rupture are working as intended. These types of computer codes are often used to understand how earthquakes operate, and they produce simulation results that include earthquake size, amounts of fault slip, and the patterns of ground shaking and crustal deformation. The benchmark exercises examine a range of features that scientists incorporate in their dynamic earthquake rupture simulations. These include implementations of simple or complex fault geometry, off‐fault rock response to an earthquake, stress conditions, and a variety of formulations for fault friction. Many of the benchmarks were designed to investigate scientific problems at the forefronts of earthquake physics and strong ground motions research. The exercises are freely available on our website for use by the scientific community.

Seismological Research Letters↗

Integrate urban‐scale seismic hazard analyses with the U.S. National Seismic Hazard Model

For more than 20 yrs, damage patterns and instrumental recordings have highlighted the influence of the local 3D geologic structure on earthquake ground motions (e.g., M "> M M 6.7 Northridge, California, Gao et al. , 1996 ; M "> M M 6.9 Kobe, Japan, Kawase, 1996 ; M "> M M 6.8 Nisqually, Washington, Frankel, Carver, and Williams, 2002 ). Although this and other local‐scale features are critical to improving seismic hazard forecasts, historically they have not been explicitly incorporated into the U.S. National Seismic Hazard Model (NSHM, national model and maps), primarily because the necessary basin maps and methodologies were not available at the national scale. Instead, the U.S. Geological Survey (USGS), its partners, and external groups developed urban seismic hazard maps (urban models and maps) that consider detailed site effects in local areas (e.g., Wong et al. , 2002 ; Cramer et al. , 2006 ; Frankel et al. , 2007 ; Graves et al. , 2011 ). The disconnect between the urban and national hazard models, however, means that the national models, which underlie U.S. building codes and other applications, do not make use of all of the scientific results informing earthquake ground‐shaking hazards. We recommend that future U.S. national seismic hazard assessment make use of all available regional information, including that in urban models. In this column, we describe the roles of and differences between the urban and national models, and discuss the obstacles to and benefits of integrating the urban models with the national model. Future progress on this issue will require further coordination and implementation efforts between the USGS and external partners.

Seismological Research Letters↗