Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Seismological Research Letters”

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 271 records · Page 15Linked to original sources
USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗
USGS2024

Summary of the discussions during 2023 SSA topical meeting on “Future Directions for Physics-Based Ground Motion Modeling”

The Seismological Society of America (SSA) topical conference, Future Directions for Physics‐Based Ground Motion Modeling, was held in Vancouver, Canada, on 10–13 October 2023, co‐sponsored by the Seismological Society of Japan and co‐chaired by Annemarie Baltay of the U.S. Geological Survey and Hiroshi Kawase of Kyoto University. This meeting brought together many researchers and practitioners interested in modeling, observing, and utilizing ground‐motion models (GMMs). Scientists gathered to discuss complex kinematic and dynamic rupture simulation approaches, empirical representations of the earthquake source, site and path effects, physical modeling of the recording site, challenges for model extrapolation, and overall prediction accuracy and simulation validation. The four‐day meeting included many posters as well as oral presentations, with each session followed by lively discussion sections, upon which we report here.

Seismological Research Letters↗
USGS2024

Preface to focus section on new frontiers and advances in global seismology

Over the last century, many of the fundamental advances in our understanding of the solid Earth have been underpinned by seismic observations recorded on long‐running networks of globally distributed seismic instruments (e.g., Agnew et al. , 1976; Romanowicz et al. , 1984; Hanka and Kind, 1994; Peterson and Hutt, 2014; Ringler et al. , 2022a). During this time, seismic data quality and the speed of dissemination have improved substantially from early analog paper records to digital, very broadband data transmitted in near‐real time (Steim, 2015) and rapidly archived in online data repositories with...

Seismological Research Letters↗
USGS2024

How does the onset of offset influence geologic slip rates?

Geologic slip rates are typically based on the displacement accrued by a geomorphic or stratigraphic feature and the age of the offset feature. Because slip rates are commonly calculated by dividing the displacement of a faulted marker by its age, they contain two open time intervals: the elapsed time between the age of an offset feature and the age of the earthquake that displaced the feature, and the time between the present‐day and the most recent earthquake. Here, we explore the influence of including unconstrained open intervals in geologic slip rate calculations. We test the degree to which these open intervals affect geologic slip rates and their uncertainties, and we find that their influence depends primarily on mean earthquake recurrence intervals (RIs). Slip rates on faults with longer RIs, such as the Wasatch fault, can be greatly influenced by an increase of up to 20% when accounting for open intervals. In contrast, slip rates on faults with shorter RIs, such as the San Andreas fault, are only slightly influenced by the assumption that slip rates calculated over open intervals approximate those calculated over closed intervals. Our analyses indicate that faults with moderate slip rates (∼0.2–5 mm/yr) are sensitive to both open interval effects themselves, as well as methods to quantify and account for these effects. We re‐evaluate how slip rates are calculated and defined in displacement–time space using published deformation records. We explore the utility of assigning a probability distribution to the initiation of offset of the oldest faulted feature and the timing of the most recent earthquake (MRE). We find that calculating geologic slip rates without using probability distributions that capture the timing of the MRE and the onset of offset of the oldest faulted feature, especially on slow‐to‐moderate slip rate faults, can lead to systematic underestimation of average geologic slip rates.

Seismological Research Letters↗
USGS2024

Recent expansion of the Cascades Volcano Observatory geophysical network at Mount Rainier for improved volcano and lahar monitoring

The U.S. Geological Survey Cascades Volcano Observatory (CVO) recently expanded its continuous geophysical monitoring at Mount Rainier, an active stratovolcano in Washington state. CVO monitors volcanoes in Oregon, Washington, and Idaho to characterize volcanic systems and detect unrest. Mount Rainier has a history of large lahar occurrences in the Holocene, including at least one that may not have been associated with volcanic activity. Pierce County, Washington, is one of the areas most at risk from large lahars. In the 1990s, CVO collaborated with Pierce County to install the Rainier lahar detection system (RLDS), an automated system designed to detect large lahars in high‐risk drainages and mitigate hazards to heavily populated areas. The system was designed to detect lahars within 5–10 min of their occurrence and alert authorities of the need to evacuate populated low‐lying areas before lahar arrival. In addition, CVO and the Pacific Northwest Seismic Network (PNSN) maintained and expanded a network of seismic and geodetic monitoring stations on and near the edifice to provide adequate volcano monitoring capabilities. Since 2016, CVO has worked to upgrade the existing RLDS and to expand its capabilities into other drainages around Mount Rainier. This expansion includes installation of 25 new broadband seismic stations with many including infrasound along high‐risk drainages, as well as support for equipment upgrades at existing PNSN and CVO volcano monitoring sites. All stations transmit continuous, near‐real‐time data with dramatically improved spatial coverage for volcano monitoring and lahar hazard mitigation compared to the previous system.

Washington↗
USGS2024

Fault-dike-magma interactions inferred from transcrustal conical structures under Akutan volcano

Volcano monitoring and eruption forecasting requires accurate characterization of transcrustal magmatic structures in order to place volcanic unrest in context within the system where it occurs. Structural imaging using local seismicity is limited to seismogenic depths. Here, we exploit arrivals in teleseismic receiver functions that change polarity with backazimuth to image two surfaces beneath Akutan volcano in the Aleutian arc. The two surfaces delineate an upper to midcrustal inverted conical volume that deepens and thickens away from the volcanic center, with thicknesses of 3 to 13 km. The top of the volume is at depths of 2-3 km below sea level at distances of ~5-15 km from the caldera center. The bottom is at depths of 7-15 km at the same distances, and the cone’s thickness increases outward from ~5 to ~10 km depth. The signal is best fit by a volume with anisotropy with fast symmetry planes that dip outwards from the center and downward increases in shear velocity at both interfaces. The upper boundary coincides with the top of Akutan’s volcano-tectonic (VT) seismogenic zone, with the VT seismicity exhibiting outward dipping planar features matching the anisotropic fast plane orientation within the volume. The bottom of the anisotropic volume is below the termination depth of the majority of the VT seismicity and is therefore likely associated with the brittle-ductile transition. Long-period (LP) events associated previously with magma movement are concentrated below the anisotropic VT volume. Because of the strong spatial association with VT seismicity, we interpret the volume as consisting of concentric outward dipping faults/dikes that align the seismogenic response to stress changes from magmatic processes. Our observations map this volume independent of the present-day seismicity distribution and thus provide a spatially more complete image of the magmatic system.

Alaska↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
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

On algorithmically determined versus traditional macroseismic intensity assignments

The utility of macroseismic data, defined as the effects of earthquakes on humans and the built environment, has been increasingly recognized following the advent of online systems that now produce unprecedented volumes of macroseismic intensity information. Contributed reports from the U.S. Geological Survey “Did You Feel It?” (DYFI) system ( Wald et al ., 1999 ) are used to generate intensity values with an algorithm based on seminal work by Dengler and Dewey (1998) . The algorithm was developed initially to reproduce intensity values assigned by expert opinion using questionnaire results collected by telephone survey. In this article, I discuss reasons why intensity values from (self‐selected) DYFI responses can differ from values that would be assigned by expert opinion given more complete data from randomly selected participants. For example, with the data used by Dengler and Dewey (1998) , intensities near 4 could be determined from the percentage of people who felt shaking in each town. With less spatially rich data from self‐selected participants, this percentage often cannot be determined reliably. Audible noises are key additional diagnostic criteria for modified Mercalli intensity (MMI) 4, but, although the DYFI system includes a question about noise, following Dengler and Dewey (1998) , the DYFI algorithm does not include a noise indicator. At the upper end of the scale, as defined the DYFI algorithm yields a maximum intensity value of 9.05, nominally corresponding to peak ground acceleration of 75% g . These and other factors can result in DYFI values that are low compared to traditional MMI values assigned using expert opinion, even absent factors that can bias traditional MMI assignments. Modern ground‐motion intensity conversion equations determined using DYFI intensities are expected to be appropriate for DYFI intensities, but the results of this study suggest that biases may be introduced if DYFI and traditional intensities are assumed to be interchangeable.

Seismological Research Letters↗
USGS2024

Preface to the focus section on volcano monitoring in the Americas

From the Andes to the Aleutian Islands, the Americas are rich with volcanism that spans a diverse range of tectonic settings, eruptive styles, levels of activity, and hazards. Over the past 120 yr, the Americas have witnessed catastrophic volcanic eruptions that have significantly impacted nearby populations. Notable events include the 8 May 1902, pyroclastic density current from Mount Pelée in Martinique, which resulted in the loss of over 28,000 lives in Saint‐Pierre, with only one or two survivors (Lacroix, 1904), and the 1985 eruption of Nevado del Ruiz in Colombia, which triggered a lahar and left an estimated...

Seismological Research Letters↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2024

Distinguishing natural sources from anthropogenic events in seismic data

As seismic data are increasingly used to investigate a diverse range of subsurface phenomena beyond regular fast-rupturing earthquakes (Peng and Gomberg, 2010; Beroza and Ide, 2011), it is important to acknowledge that human-generated ground vibrations may be mistaken for naturally generated subsurface processes (Larose et al., 2015; Li et al., 2018). Correct discrimination of natural processes from anthropogenic noise is especially pressing given the trend in seismic detection research toward automated algorithms and machine learning methods (Yoon et al., 2015; Kong et al., 2019;Mousavi and Beroza, 2022) and the growth in seismic data collection in new environments such as urban and industry settings (e.g., Díaz et al.,2017).

Seismological Research Letters↗
USGS2025

Seismic intensities, ground motions and basin amplification from the 2022 Mw 6.4 Ferndale, California intraslab earthquake

On December 20, 2022, a Mw 6.4 earthquake occurred at a depth of 18 km within the subducting Gorda plate in the Mendocino Triple Junction (MTJ), one of the most seismically active regions of the contiguous United States, causing widespread damage to local communities. Here we document the seismic intensities, ground motions, and basin amplification effects recorded by this earthquake across Humboldt County as part of an ongoing scientific effort to understand subduction zone earthquake hazards. Modified Mercalli Intensity (MMI) values from our post-earthquake field survey report shaking intensities as high as VIII (Severe). Strong ground motion data from 54 seismic stations were processed to calculate amplitude and frequency content parameters across Northern California. The maximum calculated geometric mean PGA and PGV are 1053 cm/s 2 and 52 cm/s (respectively), both recorded within the 3 km deep Eel River sedimentary basin. Comparisons with four published Ground Motion Prediction Equations indicate that PGA and PGV measurements align with expected attenuation-distance patterns for intraslab earthquakes of this nature. Within the Eel River Basin, ground motions for frequencies above 1 Hz are amplified, with respect to reference stations, by factors greater than 2. Our findings suggest that peak ground motions were mainly caused by sedimentary basin site-effects within the Eel River basin, although rupture directivity may have also increased ground motion amplitudes locally. While the Ferndale area is primarily impacted by shallow intraslab earthquakes, our results also raise questions about site-response and basin amplification hazards from a potential megathrust earthquake. The high seismicity rates of the southernmost Cascades call for stronger regional preparedness and improved strategies to mitigate the effects of such a large-scale disaster.

California↗
USGS2025

Surface rupture from an aftershock: Remote observations from the January 2024 Wushi (Aykol), China, earthquakes

On 22 January 2024, a moment magnitude ( M w ⁠⁠ ) 7.0 earthquake occurred in the southern Tian Shan region of western China, followed by an M w 5.7 aftershock on 29 January. The M w 7.0 mainshock occurred on a northwest‐dipping fault plane, at a hypocentral depth of 13 km. Both Interferometric Synthetic Aperture Radar (InSAR) and optical satellite imagery observations indicate a lack of surface rupture for the M w 7.0 earthquake. However, mapping and subpixel correlation of high‐resolution optical satellite imagery and InSAR observations spanning the 29 January M w 5.7 thrust‐faulting aftershock support the interpretation of fault surface rupture. We observed ∼4.8 km of discontinuous surface rupture associated with the M w 5.7 aftershock. The surface rupture trace trends subparallel to the northeast–southwest bedding orientation within the limb of a syncline. Subpixel image correlation and digital elevation model (DEM) differencing indicate mean on‐fault offsets of 1 m, with a peak of ∼1.4 m, along the fault. Given these observations, we interpret the M w 5.7 aftershock produced surface rupture on an unmapped flexural‐slip fault within a syncline in the hanging wall of the Maidan fault or on a back thrust associated with the Maidan fault. Pre‐existing fault features along the surface rupture trace in the landscape from pre‐earthquake DEMs indicate that the fault that ruptured in the M w 5.7 aftershock may have ruptured in the past. Surface rupture only from an aftershock in a seismic sequence is a critical observation that highlights the need to understand secondary or distributed faulting in larger earthquakes for seismic and fault displacement hazard analysis.

Tian Shan region↗
USGS2025

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗