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 289 records · Page 16Linked to original sources
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗
USGS2025

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina↗
USGS2025

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California↗
USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗
USGS2025

Revisiting seismological discoveries of the inner core

Seismology has been used as a tool for understanding the current physical properties of the interior of the Earth and its dynamic evolution with remarkable success over the last century. Much of this progress is due to the ever‐expanding set of high‐quality quantitative observations of teleseismic waveforms recorded at seismographic stations worldwide. In this work, we revisit historical seismological studies that helped first identify a core distinct in physical properties from the overlying mantle, followed by the detection of an inner core that was eventually verified to be solid based on normal‐mode eigenperiods. Along with a brief overview of past studies of the Earth’s inner core, we examine the reproducibility of these results and discuss how historical data compare against modern observations. After accounting for past normal‐mode misidentifications, we confirm that introducing a solid inner core is required to afford significant improvements in fits to both radial modes and core‐sensitive spheroidal overtones. Strong shear dissipation in the inner core of the radial reference Earth model, REM1D ( ⁠⁠ Q u = 89.54 ), fits the reference datasets of both normal‐mode eigenperiods and quality factors accounting for physical dispersion. Because a liquid region would only have bulk dissipation, a narrow range of low Q u values that are preferred by the reference datasets affords additional evidence of a solid inner core. In addition, we find that there is little systematic bias in the timing accuracy of historical data, although large variances exist. Investigations into the temperature, composition, and evolution of the inner core, as well as the reproducibility of past studies, can benefit from the reconciliation of historical and modern seismological datasets.

Seismological Research Letters↗
USGS2026

Quantifying seismic properties of a river channel at Mount Rainier for use in debris flow monitoring and analysis

Theoretical models that relate debris flow properties to their seismic signature suggest seismic methods may be used to remotely characterize properties of these events. However, the complexity of debris flow sources and poorly constrained material properties in near‐surface environments limit our ability to determine important attributes of debris flows, such as volume and particle size distribution, from seismic records alone. In this study, we explore the sensitivity of debris flow seismic signals to subsurface seismic characteristics using an established debris flow seismicity model and subsurface properties derived from active‐source measurements in the Tahoma Creek stream channel at Mount Rainier, United States. Using refraction and multichannel analysis of surface waves, we estimate 1D primary and secondary wave velocity profiles to a depth of 11 m and calculate the frequency‐varying phase velocity ( ⁠⁠ v c ) and Rayleigh‐wave quality factor (Q R ⁠⁠ ) for frequencies between 9 and 50 Hz. We find that the Tahoma Creek stream channel has low v c ⁠ , varying with frequency between 226 and 434 m/s, and is highly attenuating, with Q R estimated below 13.3 at all analyzed frequencies. We model the seismic signal of a hypothetical debris flow in Tahoma Creek using our measured values and compare the results against the same model but varying v c and Q R over a range of frequency‐independent, single values commonly used in the literature when measurements are not available ( ⁠⁠ v c = 250-750 m/s, Q R = 3-33 ⁠ ). We find that the power spectral densities of the modeled debris flows vary by orders of magnitude within the subset of our test values, highlighting the benefits of measuring material properties when using modeled debris flow seismic signals for quantitative monitoring.

Washington↗
USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Making many out of one: Synthetic geologic deformation model distributions for use in USGS NSHM25‐PRVI Puerto Rico-U.S. Virgin Island update

A key use‐case of geologic slip rates is within deformation models used in probabilistic seismic hazard analyses. Field‐derived geologic slip rates have formed the cornerstone of deformation models in such applications for decades. Recent advancements in seismic hazard analyses have expanded the use of faults for which geologic slip rates are not well constrained using categorical slip rate estimates. Because of these advancements, application of a geologic deformation model for use in 2025 U.S. Geological Survey National Seismic Hazard Model Puerto Rico‐U.S. Virgin Islands (NSHM25‐PRVI) proved challenging due to: (1) a lack of field‐based geologic slip rates, and (2) a lack of epistemic uncertainty distributions within a broad range of estimated slip rates. Preliminary versions of the NSHM25‐PRVI model sampled these slip rate bins in a coincident manner along preferred and extreme value branches, which yielded untenable correlations in mean hazard results. To minimize the influence of correlated uncertainties amid these challenges, we develop a synthetic epistemic uncertainty distribution for deformation rate on each crustal fault. Each fault has a weighting schema across four possible distribution shapes: uniform, normal, triangular favoring local minima, and triangular favoring local maxima. The synthetic distributions are then sampled several times for each logic tree branch. The results provide a more realistic distribution of rates across the study region as compared with using correlated extrema sampling. This exploration of our method in a small region like PRVI can pave the way for larger‐scale, more complicated applications (e.g., western United States).

Puerto Rico, US Virgin Islands↗
USGS2025

Late Pleistocene kinematics of the Great Southern Puerto Rico Fault Zone, Puerto Rico

Several onshore faults in southern Puerto Rico have recently been recognized as Quaternary active. However, the kinematics of these faults, particularly any lateral component, remain largely unconstrained. It is difficult to characterize low strain‐rate faults, partially due to extensive erosional and anthropogenic landscape modification, steep relief, and frequent landsliding, limiting the preservation of geomorphic features that could serve as recorders of fault motion. Here, we constrain the kinematics along sections of the Great Southern Puerto Rico Fault Zone (GSPRFZ) on the southern coastal plain of Puerto Rico. We integrate ∼1‐m‐resolution light detection and ranging (lidar)‐derived topography, historical air photos, and field mapping to identify a series of ∼50–1200‐m‐long fault scarps and lineaments that trend northwest–southeast and extend for ≥25 km across the southern coastal plain. Fault scarps are primarily south facing, cut across topography, and displace Quaternary deposits and landforms. We document multiple offset geomorphic markers, including channel thalwegs and interfluves formed in deposits previously mapped as Quaternary piedmont alluvial plain. We observe both vertical (south‐side‐down) and right‐lateral meter‐scale displacements, which indicate that the GSPRFZ accommodates right‐lateral oblique motion in the late Pleistocene, consistent with northeast motion of the Puerto Rico and the Virgin Islands microplate away from the Hispaniola block.

Puerto Rico↗
USGS2026

Opportunities and challenges for incorporating lithium-ion batteries into geophysical power systems

Batteries play a critical role at seismic and other geophysical stations in ensuring that the equipment receives uninterrupted power to make continuous observations of the Earth. Often sited in remote environments under a wide range of environmental conditions, these stations have almost ubiquitously used conventional lead‐acid batteries. However, recent advances in rechargeable lithium‐ion battery technology have led to the development of lithium‐iron phosphate (LFP) batteries that can replace lead‐acid batteries in energy storage systems. In the past three years, LFP battery production has grown over an order of magnitude, reducing their cost and increasing their availability. In this article, we examine the benefits and challenges associated with incorporating several types of LFP batteries into the power systems of geophysical stations. We review how fundamental differences in battery chemistry enable LFP batteries to weigh less than half of comparable lead‐acid batteries and may have longer service lives but have the added complexity of being managed by an integrated battery monitoring system (BMS) microcontroller. Through a controlled heat chamber experiment, we demonstrate that LFP batteries can be ideally suited for use in high‐temperature environments by offering longer service lives than lead‐acid batteries. However, for safety reasons, the BMS prevents most LFP batteries from being charged at temperatures <0°C, which complicates their use in colder regions and can lead to equipment damage if the charging system is not appropriately designed for this battery type. Finally, we address some of the unique challenges associated with adopting lithium batteries, including unexpected substantial (≤14 W) power consumption from the BMS, safety concerns, and shipping restrictions. It is our goal to help inform station operators and geophysicists about potential benefits and issues that can arise from using LFP batteries at seismic and geophysical stations in a range of environments.

Seismological Research Letters↗
USGS2025

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2025

Reflections on a trio of North American earthquakes in 1925

In 1925, three moderately large damaging earthquakes occurred in North America over four months: the 28 February (local time; LT) M 6.2 Charlevoix, 27 June (LT) M 6.6 Montana, and 29 June M 6.5 Santa Barbara earthquakes. The centennial anniversaries of these events motivated this retrospective consideration focused on the ground motions generated by the three events, including a reconsideration of early intensity assignments for the Montana earthquake. At the time, these three earthquakes appeared to support the arguments of some geologists who downplayed the severity of seismic hazard in southern California relative to other parts of the country. Some of the arguments advanced at that time, for example that Los Angeles “has the least to fear from ‘Acts of God’ of any city under the American flag,” ( Hill, 1928 ) sound naïve if not laughable now, but a comparison of well‐constrained shaking distributions for the three earthquakes reveals the dramatic difference in wave propagation efficiency in western versus eastern North America (ENAM), which leads to moderate ENAM events being felt to much larger distances. At M 6.2, the 1925 Charlevoix earthquake was a notably large event in ENAM. This earthquake was the largest event in eastern Canada since 1870 and caused damage in the epicentral region in addition to towns as far away as 200 km, with felt shaking extending over 1000 km. In contrast, felt shaking from the Santa Barbara earthquake barely extended beyond ∼200 km. Compiling published intensity distributions for larger ENAM earthquakes, we show that perceptible earthquake shaking is not uncommon in ENAM over century time scales, but experience with weakly felt shaking may incline people to downplay potential earthquake risk.

California, Montana, Quebec↗
USGS2026

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗