USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Seismological Record”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23Linked to original sources
USGS2025

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California↗
USGS2018

The 19 September 2017 M 7.1 Puebla‐Morelos earthquake: Spectral ratios confirm Mexico City zoning

One important element of understanding basin response to strong shaking is the analysis of spectral ratios, which may provide information about the dominant frequency of ground motion at specific locations. Spectral ratios computed from accelerations recorded by strong‐motion stations in Mexico City during the mainshock of the 19 September 2017 M "> M M 7.1 Puebla‐Morelos earthquake reveal predominate periods consistent with those mapped in the 2004 Mexican seismic design code. Furthermore, the predominant periods thus computed validate those studies using mainshock and aftershock recordings of the handful strong‐motion stations that recorded the 19 September 1985 M "> M M 8.1 Michoacán earthquake. Even though the number of stations in each of the zones (zones I, II, IIIa, b, c, and d) is not the same, they still allow confirmation of site frequencies (periods) attributable to the specific zones (particularly those in zones IIIa, b, c, and d). Spectral ratios are computed with two different methods: (1) horizontal to horizontal (H/H) ratio of smoothed amplitude spectrum of a horizontal channel in direction X of a station with respect to the smoothed amplitude spectrum of the horizontal channel in the same X direction of a reference stiff soil (or rock) station and (2) horizontal to vertical (H/V) ratio (or also known as the Nakamura method) of both horizontal (H) and vertical (V) channels of the same station. We show a comparison of the identified frequencies (periods) derived by both methods and find they are very similar and in good agreement with those indicated in the zoning maps of Mexico City in the 2004 seismic design code.

Mexico City↗
USGS2019

Characteristics and spatial variability of wind noise on near-surface broadband seismometers

By coupling with the ground, wind causes ground motion that appears on seismic records as noise across a wide bandwidth. This wind-generated noise can drown out important features such as small earthquakes and prevent observation of normal modes from large earthquakes. Because the wind field is heterogeneous at local scales due to structures, diurnal heating, and topography, wind-induced seismic noise may be different on seismometers installed just meters apart. We have investigated the spatial variability of wind-induced noise using two weather sensors separated by approximately ~100 m and co-located with one deep borehole and four near-surface broadband seismometers. We found that at longer periods (>5 s), increasing wind speed causes increases in noise on the horizontal components of seismometers. While this has been previously observed, we also measured a γ2-coherence of less than 0.2 between the wind speed, wind direction, and the pressure recorded by our weather stations. We further observed a loss of coherence between the vertical components of our seismometers from 8 s to 20 s period. The amplitude of the drop-in coherence appears to depend on the substrate surrounding the seismometer. Based on two previously-developed theoretical models, we found that the local material surrounding the sensor could be amplifying the wind-generated noise. We also investigated the frequency dependence of wind-induced noise and found that the dominant source of high-frequency seismic noise at some sites could be anthropogenic rather than induced by wind. Additionally, we estimated the linear relationship between the root mean squares (RMS) of wind speed and RMS seismic velocity for all sensors, finding substantial variability between different installments. A more detailed understanding of the complex processes by which wind-induced noise is generated can inform the installation of sensors and the development of methods for mitigation of these effects, thus improving the overall quality of seismic data.

New Mexico↗
USGS2019

Ground-motion attenuation in the Sacramento-San Joaquin delta, California, from 14 Bay Area earthquakes, including the 2014 M 6.0 South Napa earthquake

Peak ground motions (acceleration and velocity) radiated by earthquakes in the San Francisco Bay area and recorded within the Sacramento–San Joaquin Delta generally attenuate faster with distance than the Next Generation Attenuation-West2 ground-motion prediction equations (GMPEs). We evaluate the attenuation for a wide set of paths into the Delta by analyzing recorded ground motions from fourteen 4 ≤ M < 7 earthquakes located on major Bay area faults: the San Andreas, Calaveras, Hayward, West Napa, and Green Valley faults. We select stations within azimuthal ranges of 38°–114° into the Delta and calculate the residuals of the peak ground motions relative to the Boore et al. (2014) GMPEs. We then fit the natural log of these peak ground acceleration and peak ground velocity residuals for each earthquake to the function a−krγ, in which a is an event term and krγ is the differential attenuation. Although there is some variation in the differential attenuation obtained for each earthquake, the peak ground motions from most of the 14 events attenuate faster than predicted by the Boore et al. (2014) GMPEs. The differential attenuation does not appear to depend on azimuth or magnitude of the earthquake; however, earthquake depth may have an effect. Our results suggest that attenuation models for the Delta can be significantly improved through regionalization, although this regionalization will increase the model complexity and the epistemic uncertainty.

California↗
USGS2021

Effect of fixing earthquake depth in ShakeAlert algorithms on performance for intraslab earthquakes

We investigate whether assuming a fixed shallow depth in the ShakeAlert network‐based earthquake early warning system is sufficient to produce accurate ground‐motion based alerts for intraslab earthquakes. ShakeAlert currently uses a fixed focal depth of 8 km to estimate earthquake location and magnitude. This is an appropriate way to reduce computational costs without compromising alert accuracy in California, where earthquakes typically occur on shallow crustal faults. In the Pacific Northwest (PNW), however, the most common moderate‐magnitude events occur within the subducting Juan de Fuca slab at depths between ∼35 and 65 km. Using a dataset of seismic recordings from 37 M w 4.5+ intraslab earthquakes from the PNW and Chile, we replay events through the Earthquake Point‐Source Integrated Code and eqInfo2GM algorithms to estimate source parameters and compute modified Mercalli intensity (MMI) alert threshold contours. Each event is replayed twice—once using a fixed 8 km depth and a second time using the actual catalog earthquake depth. For each depth scenario, we analyze MMI III and IV contours using various performance metrics to determine the number of correctly alerted sites and measure warning times. We determine that shallow depth replays are more likely to produce errors in location estimates of greater than 50 km if the event is located outside of a seismic network. When located within a seismic network, shallow and catalog depth replays have similar epicenter estimates. Results show that applying catalog earthquake depth does not improve the accuracy of magnitude estimates or MMI alert threshold contours, or increase warning times. We conclude that using a fixed shallow earthquake depth for intraslab earthquakes will not significantly impact alert accuracy in the PNW.

California, Oregon, Washington↗
USGS2011

Testing long-period ground-motion simulations of scenario earthquakes using the Mw 7.2 El Mayor-Cucapah mainshock: Evaluation of finite-fault rupture characterization and 3D seismic velocity models

Using a suite of five hypothetical finite-fault rupture models, we test the ability of long-period ( T >2.0 s) ground-motion simulations of scenario earthquakes to produce waveforms throughout southern California consistent with those recorded during the 4 April 2010 M w 7.2 El Mayor-Cucapah earthquake. The hypothetical ruptures are generated using the methodology proposed by Graves and Pitarka (2010) and require, as inputs, only a general description of the fault location and geometry, event magnitude, and hypocenter, as would be done for a scenario event. For each rupture model, two Southern California Earthquake Center three-dimensional community seismic velocity models (CVM-4m and CVM-H62) are used, resulting in a total of 10 ground-motion simulations, which we compare with recorded ground motions. While the details of the motions vary across the simulations, the median levels match the observed peak ground velocities reasonably well, with the standard deviation of the residuals generally within 50% of the median. Simulations with the CVM-4m model yield somewhat lower variance than those with the CVM-H62 model. Both models tend to overpredict motions in the San Diego region and underpredict motions in the Mojave desert. Within the greater Los Angeles basin, the CVM-4m model generally matches the level of observed motions, whereas the CVM-H62 model tends to overpredict the motions, particularly in the southern portion of the basin. The variance in the peak velocity residuals is lowest for a rupture that has significant shallow slip (<5 km depth), whereas the variance in the residuals is greatest for ruptures with large asperities below 10 km depth. Overall, these results are encouraging and provide confidence in the predictive capabilities of the simulation methodology, while also suggesting some regions in which the seismic velocity models may need improvement.

California↗
USGS2017

Modeling of high‐frequency seismic‐wave scattering and propagation using radiative transfer theory

This is a study of the nonisotropic scattering process based on radiative transfer theory and its application to the observation of the M 4.3 aftershock recording of the 2008 Wells earthquake sequence in Nevada. Given a wide range of recording distances from 29 to 320 km, the data provide a unique opportunity to discriminate scattering models based on their distance‐dependent behaviors. First, we develop a stable numerical procedure to simulate nonisotropic scattering waves based on the 3D nonisotropic scattering theory proposed by Sato (1995). By applying the simulation method to the inversion of M 4.3 Wells aftershock recordings, we find that a nonisotropic scattering model, dominated by forward scattering, provides the best fit to the observed high‐frequency direct S waves and S ‐wave coda velocity envelopes. The scattering process is governed by a Gaussian autocorrelation function, suggesting a Gaussian random heterogeneous structure for the Nevada crust. The model successfully explains the common decay of seismic coda independent of source–station locations as a result of energy leaking from multiple strong forward scattering, instead of backscattering governed by the diffusion solution at large lapse times. The model also explains the pulse‐broadening effect in the high‐frequency direct and early arriving S waves, as other studies have found, and could be very important to applications of high‐frequency wave simulation in which scattering has a strong effect. We also find that regardless of its physical implications, the isotropic scattering model provides the same effective scattering coefficient and intrinsic attenuation estimates as the forward scattering model, suggesting that the isotropic scattering model is still a viable tool for the study of seismic scattering and intrinsic attenuation coefficients in the Earth.

Bulletin of the Seismological Society of America↗
USGS2024

Reference 1D seismic velocity models for volcano monitoring and imaging: Methods, models, and applications

Seismic velocity models of the crust are an integral part of earthquake monitoring systems at volcanoes. 1D models that vary only in depth are typically used for real‐time hypocenter determination and serve as critical reference models for detailed 3D imaging studies and geomechanical modeling. Such models are usually computed using seismic tomographic methods that rely on P ‐ and S ‐wave arrival‐time picks from numerous earthquakes recorded at receivers around the volcano. Traditional linearized tomographic methods that jointly invert for source locations, velocity structure, and station corrections depend critically on having reasonable starting values for the unknown parameters, are susceptible to local misfit minima and divergence, and often do not provide adequate uncertainty information. These issues are often exacerbated by sparse seismic networks, inadequate distributions of seismicity, and/or poor data quality common at volcanoes. In contrast, modern probabilistic global search methods avoid these issues only at the cost of increased computation time. In this article, we review both approaches and present example applications and comparisons at several volcanoes in the United States, including Mount Hood (Oregon), Mount St. Helens (Washington), the Island of Hawai’i, and Mount Cleveland (Alaska). We provide guidance on the proper usage of these methods as relevant to challenges specific to volcano monitoring and imaging. Finally, we survey‐published 1D P ‐wave velocity models from around the world and use them to derive a generic stratovolcano velocity model, which serves as a useful reference model for comparison and when local velocity information is sparse.

Seismological Research Letters↗
USGS2022

Improved resolution across the Global Seismographic Network: A new era in low-frequency seismology

The Global Seismographic Network (GSN)—a global network of ≈150 very broadband stations—is used by researchers to study the free oscillations of the Earth (≈0.3–10 mHz) following large earthquakes. Normal‐mode observations can provide information about the radial density and anisotropic velocity structure of the Earth (including near the core–mantle boundary), but only when signal‐to‐noise ratios at very low frequencies are sufficiently high. Most normal‐mode observations in the past three decades have been made using Streckeisen STS‐1 vault seismometers. However, these sensors are no longer being manufactured or serviced. Candidate replacement sensors, the Streckeisen STS‐6 and the Nanometrics T‐360GSN, have been recently installed in boreholes, postholes, and vaults at several GSN stations and GSN testbeds. In this study, we examine normal‐mode spectra following three M w "> M w 8 earthquakes in 2021 and from one M w "> M w 8.2 earthquake in 2014 to evaluate the change in GSN low‐frequency performance on the vertical component. From this analysis, we conclude that the number of GSN stations capable of resolving normal modes following M w "> M w 8 earthquakes has nearly doubled since 2014. The improved observational capabilities will help better understand the radial velocity and density estimates of the Earth.

The Seismic Record↗
USGS2013

Crustal structure and fault geometry of the 2010 Haiti earthquake from temporary seismometer deployments

Haiti has been the locus of a number of large and damaging historical earthquakes. The recent 12 January 2010 Mw 7.0 earthquake affected cities that were largely unprepared, which resulted in tremendous losses. It was initially assumed that the earthquake ruptured the Enriquillo Plantain Garden fault (EPGF), a major active structure in southern Haiti, known from geodetic measurements and its geomorphic expression to be capable of producing M 7 or larger earthquakes. Global Positioning Systems (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data, however, showed that the event ruptured a previously unmapped fault, the Léogâne fault, a north‐dipping oblique transpressional fault located immediately north of the EPGF. Following the earthquake, several groups installed temporary seismic stations to record aftershocks, including ocean‐bottom seismometers on either side of the EPGF. We use data from the complete set of stations deployed after the event, on land and offshore, to relocate all aftershocks from 10 February to 24 June 2010, determine a 1D regional crustal velocity model, and calculate focal mechanisms. The aftershock locations from the combined dataset clearly delineate the Léogâne fault, with a geometry close to that inferred from geodetic data. Its strike and dip closely agree with the global centroid moment tensor solution of the mainshock but with a steeper dip than inferred from previous finite fault inversions. The aftershocks also delineate a structure with shallower southward dip offshore and to the west of the rupture zone, which could indicate triggered seismicity on the offshore Trois Baies reverse fault. We use first‐motion focal mechanisms to clarify the relationship of the fault geometry to the triggered aftershocks.

Bulletin of the Seismological Society of America↗
USGS2013

The Mw 5.8 Mineral, Virginia, earthquake of August 2011 and aftershock sequence: constraints on earthquake source parameters and fault geometry

The M w 5.8 earthquake of 23 August 2011 (17:51:04 UTC) (moment, M0 5.7×10 17 N·m) occurred near Mineral, Virginia, within the central Virginia seismic zone and was felt by more people than any other earthquake in United States history. The U.S. Geological Survey (USGS) received 148,638 felt reports from 31 states and 4 Canadian provinces. The USGS PAGER system estimates as many as 120,000 people were exposed to shaking intensity levels of IV and greater, with approximately 10,000 exposed to shaking as high as intensity VIII. Both regional and teleseismic moment tensor solutions characterize the earthquake as a northeast‐striking reverse fault that nucleated at a depth of approximately 7±2 km. The distribution of reported macroseismic intensities is roughly ten times the area of a similarly sized earthquake in the western United States (Horton and Williams, 2012). Near‐source and far‐field damage reports, which extend as far away as Washington, D.C., (135 km away) and Baltimore, Maryland, (200 km away) are consistent with an earthquake of this size and depth in the eastern United States (EUS). Within the first few days following the earthquake, several government and academic institutions installed 36 portable seismograph stations in the epicentral region, making this among the best‐recorded aftershock sequences in the EUS. Based on modeling of these data, we provide a detailed description of the source parameters of the mainshock and analysis of the subsequent aftershock sequence for defining the fault geometry, area of rupture, and observations of the aftershock sequence magnitude–frequency and temporal distribution. The observed slope of the magnitude–frequency curve or b‐value for the aftershock sequence is consistent with previous EUS studies (b=0.75), suggesting that most of the accumulated strain was released by the mainshock. The aftershocks define a rupture that extends between approximately 2–8 km in depth and 8–10 km along the strike of the fault plane. Best‐fit modeling of the geometry of the aftershock sequence defines a rupture plane that strikes N36°E and dips to the east‐southeast at 49.5°. Moment tensor solutions of the mainshock and larger aftershocks are consistent with the distribution of aftershock locations, both indicating reverse slip along a northeast–southwest striking southeast‐dipping fault plane.

Virginia↗
USGS2022

High‐frequency rupture processes of the 2014 Mw 8.2 Iquique and 2015 Mw 8.3 Illapel, Chile, earthquakes determined from strong‐motion recordings

Strong‐motion recordings of the 2014 M w "> M w Mw 8.2 Iquique and 2015 M w "> M w Mw 8.3 Illapel, Chile, earthquakes were analyzed to determine rupture propagation and the location, timing, and strength of subevents that produce most of the high‐frequency (≥1 Hz) ground motions. A moving window,cross‐correlation analysis of recordings from a local dense array, band‐pass filtered at 1 Hz, directly shows that the Iquique earthquake ruptured to the southeast over a distance of about 60 km. Array analysis of lower frequency energy (0.03–0.1 Hz) indicates that it occurred updip of the high‐frequency rupture. A methodology was developed for inverting the envelopes of acceleration records (1–5 Hz) to map high‐frequency source factors on the rupture zone and was applied to the two earthquakes. Waveforms of M w "> M w Mw 6 earthquakes were used as empirical Green’s functions in the inversions. High‐frequency subevents within the two M w "> M w Mw 8 earthquakes were located at depths ranging from 25 to 55 km and mostly occurred downdip of the peak slip of these earthquakes. Fourier spectral ratios of the Iquique mainshock with respect to M w "> M w Mw 5–6 aftershocks were fit to determine their stress drops. The stress drops were roughly constant from M w "> M w Mw 5 to 8 at 10–20 MPa. A compound rupture model is described in which subevents occur in areas of spatially heterogeneous strength and stress on the rupture, and produce the high‐frequency radiated energy of the overall earthquake, but are not located in the areas of peak slip. The stress drop of the overall earthquake is shown to equal the root mean square stress drop of subevents averaged over the rupture area.

Iquique, Illapel↗
USGS2011

88 hours: The U.S. Geological Survey National Earthquake Information Center response to the March 11, 2011 Mw 9.0 Tohoku earthquake

The M 9.0 11 March 2011 Tohoku, Japan, earthquake and associated tsunami near the east coast of the island of Honshu caused tens of thousands of deaths and potentially over one trillion dollars in damage, resulting in one of the worst natural disasters ever recorded. The U.S. Geological Survey National Earthquake Information Center (USGS NEIC), through its responsibility to respond to all significant global earthquakes as part of the National Earthquake Hazards Reduction Program, quickly produced and distributed a suite of earthquake information products to inform emergency responders, the public, the media, and the academic community of the earthquake's potential impact and to provide scientific background for the interpretation of the event's tectonic context and potential for future hazard. Here we present a timeline of the NEIC response to this devastating earthquake in the context of rapidly evolving information emanating from the global earthquake-response community. The timeline includes both internal and publicly distributed products, the relative timing of which highlights the inherent tradeoffs between the requirement to provide timely alerts and the necessity for accurate, authoritative information. The timeline also documents the iterative and evolutionary nature of the standard products produced by the NEIC and includes a behind-the-scenes look at the decisions, data, and analysis tools that drive our rapid product distribution.

Tohoku↗
USGS2006

Long-period effects of the Denali earthquake on water bodies in the Puget Lowland: Observations and modeling

Analysis of strong-motion instrument recordings in Seattle, Washington, resulting from the 2002 Mw 7.9 Denali, Alaska, earthquake reveals that amplification in the 0.2-to 1.0-Hz frequency band is largely governed by the shallow sediments both inside and outside the sedimentary basins beneath the Puget Lowland. Sites above the deep sedimentary strata show additional seismic-wave amplification in the 0.04- to 0.2-Hz frequency range. Surface waves generated by the Mw 7.9 Denali, Alaska, earthquake of 3 November 2002 produced pronounced water waves across Washington state. The largest water waves coincided with the area of largest seismic-wave amplification underlain by the Seattle basin. In the current work, we present reports that show Lakes Union and Washington, both located on the Seattle basin, are susceptible to large water waves generated by large local earthquakes and teleseisms. A simple model of a water body is adopted to explain the generation of waves in water basins. This model provides reasonable estimates for the water-wave amplitudes in swimming pools during the Denali earthquake but appears to underestimate the waves observed in Lake Union.

Bulletin of the Seismological Society of America↗
USGS2016

Simulation of earthquake ground motions in the eastern United States using deterministic physics‐based and site‐based stochastic approaches

Earthquake ground‐motion recordings are scarce in the central and eastern United States (CEUS) for large‐magnitude events and at close distances. We use two different simulation approaches, a deterministic physics‐based method and a site‐based stochastic method, to simulate ground motions over a wide range of magnitudes. Drawing on previous results for the modeling of recordings from the 2011 M w 5.8 Mineral, Virginia, earthquake and using the 2001 M w 7.6 Bhuj, India, earthquake as a tectonic analog for a large magnitude CEUS event, we are able to calibrate the two simulation methods over this magnitude range. Both models show a good fit to the Mineral and Bhuj observations from 0.1 to 10 Hz. Model parameters are then adjusted to obtain simulations for M w 6.5, 7.0, and 7.6 events in the CEUS. Our simulations are compared with the 2014 U.S. Geological Survey weighted combination of existing ground‐motion prediction equations in the CEUS. The physics‐based simulations show comparable response spectral amplitudes and a fairly similar attenuation with distance. The site‐based stochastic simulations suggest a slightly faster attenuation of the response spectral amplitudes with distance for larger magnitude events and, as a result, slightly lower amplitudes at distances greater than 200 km. Both models are plausible alternatives and, given the few available data points in the CEUS, can be used to represent the epistemic uncertainty in modeling of postulated CEUS large‐magnitude events.

Bulletin of the Seismological Society of America↗
USGS2014

Refinements to the Graves and Pitarka (2010) Broadband Ground Motion Simulation Method

This brief article describes refinements to the Graves and Pitarka (2010) broadband ground motion simulation methodology (GP2010 hereafter) that have been implemented in version 14.3 of the SCEC Broadband Platform (BBP). The updated version of our method on the current SCEC BBP is referred to as GP14.3. Our simulation technique is a hybrid approach that combines low-­‐frequency and high-­‐frequency motions computed with different methods into a single broadband response. The separate low-­‐ and high-­‐frequency components have traditionally been called “deterministic” and “stochastic”, respectively; however, this nomenclature is an oversimplification. In reality, the low-­‐frequency approach includes many stochastic elements, and likewise, the high-­‐frequency approach includes many deterministic elements (e.g., Pulido and Kubo, 2004; Hartzell et al., 2005; Liu et al., 2006; Frankel, 2009; Graves and Pitarka, 2010; Mai et al., 2010). While the traditional terminology will likely remain in use by the broader modeling community, in this paper we will refer to these using the generic terminology “low-­‐frequency” and “high-­‐ frequency” approaches. Furthermore, one of the primary goals in refining our methodology is to provide a smoother and more consistent transition between the low-­‐ and high-­‐ frequency calculations, with the ultimate objective being the development of a single unified modeling approach that can be applied over a broad frequency band. GP2010 was validated by modeling recorded strong motions from four California earthquakes. While the method performed well overall, several issues were identified including the tendency to over-­‐predict the level of longer period (2-­‐5 sec) motions and the effects of rupture directivity. The refinements incorporated in GP14.3 are aimed at addressing these issues with application to the simulation of earthquakes in Western US (WUS). These refinements include the addition of a deep weak zone (details in following section) to the rupture characterization and allowing perturbations in the correlation of rise time and rupture speed with the specified slip distribution. Additionally, we have extended the parameterization of GP14.3 so that it is also applicable for simulating Eastern North America (ENA) earthquakes. This work has been guided by the comprehensive set of validation studies described in Goulet and Abrahamson (2014) and Dreger et al. (2014). The GP14.3 method shows improved performance relative to GP2010, and we direct the interested reader to Dreger et al. (2014) for a detailed assessment of the current methodology. In this paper, we concentrate on describing the modifications in more detail, and also discussing additional refinements that are currently being developed.

California↗
USGS2013

The attenuation of Fourier amplitudes for rock sites in eastern North America

We develop an empirical model of the decay of Fourier amplitudes for earthquakes of M 3–6 recorded on rock sites in eastern North America and discuss its implications for source parameters. Attenuation at distances from 10 to 500 km may be adequately described using a bilinear model with a geometric spreading of 1/R 1.3 to a transition distance of 50 km, with a geometric spreading of 1/R 0.5 at greater distances. For low frequencies and distances less than 50 km, the effective geometric spreading given by the model is perturbed using a frequency‐ and hypocentral depth‐dependent factor defined in such a way as to increase amplitudes at lower frequencies near the epicenter but leave the 1 km source amplitudes unchanged. The associated anelastic attenuation is determined for each event, with an average value being given by a regional quality factor of Q=525f 0.45 . This model provides a match, on average, between the known seismic moment of events and the inferred low‐frequency spectral amplitudes at R=1 km (obtained by correcting for the attenuation model). The inferred Brune stress parameters from the high‐frequency source terms are about 600 bars (60 MPa), on average, for events of M>4.5.

Bulletin of the Seismological Society of America↗
USGS2022

Earthquake early warning for estimating floor shaking levels of tall buildings

This article investigates methods to improve earthquake early warning (EEW) predictions of shaking levels for residents of tall buildings. In the current U.S. Geological Survey ShakeAlert EEW system, regions far from an epicenter will not receive alerts due to low predicted ground‐shaking intensities. However, residents of tall buildings in those areas may still experience significant shaking due to the acceleration amplification caused by tall buildings’ dynamic behavior, as recently experienced by residents of the 52‐story building in downtown Los Angeles (DTLA) during the 2019 M 7.1 Ridgecrest earthquake. Using more than 400 recorded response data acquired from 77 instrumented buildings in California, here we compare the Federal Emergency Management Agency (FEMA) P‐58 and American Society of Civil Engineers (ASCE) 7‐16 simplified equations for peak floor acceleration (PFA), finding that the ASCE estimation is close to the median of data recorded in large and long‐distance events, whereas the current FEMA estimation is not suitable. In the second part of this article, four instrumented tall buildings in DTLA are extensively studied, and the performance of the simplified and response spectrum (RS) methods giving both an estimation of the free‐field horizontal peak ground acceleration (PGA) and pseudospectral acceleration is evaluated. The results show that the RS method is as accurate as the response history analysis as long as the ground‐motion RS is accurate, whereas the ASCE 7‐16 prediction is conservative. However, when ground‐motion RS or PGA is estimated for DTLA using a ground‐motion model (GMM), the performance of the RS method significantly degrades due to underestimation by the GMM at long periods. The results of this study imply that a nonergodic GMM, which may give more accurate prediction in Los Angeles, could improve the results for PFA when the building’s behavior is dominated by a few long‐period fundamental modes, as is the case for the 52‐story building in DTLA.

California↗