More Like this

1. Ground‐Motion Amplification in Cook Inlet Region, Alaska, from Intermediate‐Depth Earthquakes, Including the 2018 Mw 7.1 Anchorage Earthquake

   https://doi.org/10.1785/0220190179  

   Moschetti, Morgan P.; Thompson, Eric M.; Rekoske, John; Hearne, Michael G.; Powers, Peter M.; McNamara, Daniel E.; Tape, Carl (November 2019, Seismological Research Letters)

   null (Ed.)

   Abstract We measure pseudospectral and peak ground motions from 44 intermediate‐depth Mw≥4.9 earthquakes in the Cook Inlet region of southern Alaska, including those from the 2018 Mw 7.1 earthquake near Anchorage, to identify regional amplification features (0.1–5  s period). Ground‐motion residuals are computed with respect to an empirical ground‐motion model for intraslab subduction earthquakes, and we compute bias, between‐, and within‐event terms through a linear mixed‐effects regression. Between‐event residuals are analyzed to assess the relative source characteristics of the Cook Inlet earthquakes and suggest a difference in the scaling of the source with depth, relative to global observations. The within‐event residuals are analyzed to investigate regional amplification, and various spatial patterns manifest, including correlations of amplification with depth of the Cook Inlet basin and varying amplifications east and west of the center of the basin. Three earthquake clusters are analyzed separately and indicate spatial amplification patterns that depend on source location and exhibit variations in the depth scaling of long‐period basin amplification. The observations inform future seismic hazard modeling efforts in the Cook Inlet region. More broadly, they suggest a greater complexity of basin and regional amplification than is currently used in seismic hazard analyses. 

   more » « less
2. The 2023 US 50-State National Seismic Hazard Model: Overview and implications

   https://doi.org/10.1177/87552930231215428  

   Petersen, Mark D; Shumway, Allison M; Powers, Peter M; Field, Edward H; Moschetti, Morgan P; Jaiswal, Kishor S; Milner, Kevin R; Rezaeian, Sanaz; Frankel, Arthur D; Llenos, Andrea L; et al (February 2024, Earthquake Spectra)

   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. 

   more » « less
3. Estimating site amplification variability in Yangon, Myanmar, from a dense nodal seismic array

   https://doi.org/10.1177/87552930251352843  

   Islam, Md Mohimanul; Persaud, Patricia; Thant, Myo; Win, Khin; Sandvol, Eric (November 2025, Earthquake Spectra)

   We explore a novel acquisition geometry that can be used to estimate the linear component of site amplification using a dense nodal seismic network installed in Yangon, Myanmar’s largest city. The city is surrounded by several seismically active faults, including the Sagaing Fault, which is capable of generating M_w > 7.0 earthquakes. As part of the Irrawaddy delta system, this densely populated city sits on young water-saturated alluvium that is likely to amplify earthquake ground motions. Assessing site response is crucial for understanding the seismic hazard potential to minimize the loss of property and lives. Using a dense seismic array comprised of 110 three-component nodes, we estimated the frequency-dependent site amplification pattern of Yangon from regional (Lg) and local (Sg) seismic phases. Since this acquisition geometry is not sensitive to Q or geometric spreading, this approach provides a fast and cost-effective way to estimate the linear component of site response as a function of frequency. Our Lg and Sg site response results identify regions with high site amplification that have significantly greater seismic hazard risks for regional and local distance earthquakes. We observed consistent site response characteristics between both Lg and Sg phases. Site amplification patterns correlate well with the surficial geology and subsurface structure beneath the city. De-amplification is observed across all frequencies at stations located above an anticlinal structure composed of older Pliocene rocks (i.e. the Irrawaddy Formation). Conversely, highly amplified areas correspond to younger Pleistocene to recent alluvial plains consisting of loose, unconsolidated alluvium. We found a dominant horizontal-to-vertical spectral ratio (HVSR) peak at ∼1.0 Hz from ambient noise, likely corresponding to the thickness of unconsolidated sediments. We suggest that the growing number of nodal networks worldwide can be used to estimate frequency-dependent site amplification, addressing key data gaps in seismic hazard assessment. 

   more » « less
4. Ground-Motion Simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, Earthquake

   https://doi.org/10.1785/0220250333  

   Boyd, Oliver S; Bozdağ, Ebru; Kehoe, Haiyang L; Moschetti, Morgan P (January 2026, Seismological Research Letters)

   Abstract Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume 1D Earth structure. In this study, we use a 3D seismic velocity model to better constrain the depth and focal mechanism of the 5 April 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280 km wide by 260 km long by 77 km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey (USGS) National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the Next Generation Attenuation-East Project (NGA-East) ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic GMM, with a 76% reduction of residual variance in observed ground motions averaged over 3, 5, 7, and 10 s periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the USGS regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard. 

   more » « less
5. A Discontinuous Galerkin Method for Simulating 3D Seismic Wave Propagation in Nonlinear Rock Models: Verification and Application to the M _w 7.8 2015 Gorkha, Nepal Earthquake

   https://doi.org/10.1029/2025JB031378  

   Niu, Zihua; Gabriel, Alice‐Agnes; Wolf, Sebastian; Ulrich, Thomas; Lyakhovsky, Vladimir; Igel, Heiner (July 2025, Journal of Geophysical Research: Solid Earth)

   Abstract The nonlinear mechanical responses of rocks and soils to seismic waves play an important role in earthquake physics, influencing ground motion from source to site. Continuous geophysical monitoring, such as ambient noise interferometry, has revealed co‐seismic wave speed reductions extending tens of kilometers from earthquake sources. However, the mechanisms governing these changes remain challenging to model, especially at regional scales. Using a nonlinear damage model constrained by laboratory experiments, we develop and apply an open‐source 3D discontinuous Galerkin method to simulate regional co‐seismic wave speed changes during the 2015 M_w7.8 Gorkha earthquake. We find pronounced spatial variations of co‐seismic wave speed reduction, ranging from <0.01% to >50%, particularly close to the source and within the Kathmandu Basin, while disagreement with observations remains. The most significant reduction occurs within the sedimentary basin and varies with basin depths, whereas wave speed reductions correlate with the fault slip distribution near the source. By comparing ground motions from simulations with elastic, viscoelastic, elastoplastic, and nonlinear damage rheologies, we demonstrate that the nonlinear damage model effectively captures low‐frequency ground motion amplification due to strain‐dependent wave speed reductions in soft sediments. We verify the accuracy of our approach through comparisons with analytical solutions and assess its scalability on high‐performance computing systems. The model shows near‐linear strong and weak scaling up to 2,048 nodes, enabling efficient large‐scale simulations. Our findings provide a physics‐based framework to quantify nonlinear earthquake effects and emphasize the importance of damage‐induced wave speed variations for seismic hazard assessment and ground motion predictions. 

   more » « less