Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Transactions of the American Fisheries Society”

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.

1,248 records · Page 4Linked to original sources
USGS2026

Comparative crop yield forecasting using satellite-derived biophysical and agro-climatic predictors in Sub-Saharan Africa

Timely and accurate crop yield forecasting is central to food security early warning systems, particularly in climate-vulnerable regions. While operational forecasting frameworks commonly rely on precipitation and vegetation indices such as NDVI, their ability to provide actionable lead time remains limited. Here, we evaluate the added value of satellite-derived biophysical Essential Climate Variables (ECVs): Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FAPAR), for forecasting millet yield in Burkina Faso (BF) and maize yield in South Africa (ZA) and Malawi (MW). Using Random Forest models, we quantify forecast skill across the growing season at both national and sub-national scales. Results show that LAI and FAPAR provide effective forecast lead times of approximately 4 months in BF, 2 months in ZA, and up to 6 months in MW relative to harvest. At peak performance, Mean Absolute Percentage Error (MAPE) reaches 19.8% (LAI) and 23.8% (FAPAR) in BF, 12.0% and 9.8% in ZA, and 21.8% and 20.8% in MW, respectively. Across countries, biophysical parameters often outperform NDVI and precipitation, particularly in arid and semi-arid regions. At the sub-national level, LAI and FAPAR enable classification of administrative units into high and moderate-skill forecast units, revealing strong spatial heterogeneity linked to crop dominance. However, forecast skill declines where the target crop is not the dominant type, highlighting an important limitation for operational deployment. Overall, the findings suggest that satellite-derived biophysical parameters can provide earlier and more spatially resolved yield signals than commonly used predictors, with potential to improve the timeliness and effectiveness of food security early warning systems.

Remote Sensing Applications: Society and Environme↗
USGS2025

Estimated average annualized tsunami losses for the United States

Tsunami hazards are substantial threats to coastal communities across the United States (U.S.) and its territories. U.S. states and territories collaborate through the National Tsunami Hazard Mitigation Program (NTHMP) to develop their own tsunami-hazard information for outreach and evacuation planning. An effort to curate this tsunami-hazard information to support comprehensive risk analysis at the national level has not yet been completed. In support of this effort, the Federal Emergency Management Agency (FEMA) collaborated with the NTHMP, the National Oceanic and Atmospheric Administration (NOAA) and the U.S. Geological Survey (USGS) starting in 2023. This collaboration included the collection and analysis of existing tsunami hazard data and methods in the U.S. Tsunami subject matter experts identified and selected scientifically defensible methods for estimating the risks to buildings and populations in coastal communities. These efforts may support decision making regarding resilience policies, priorities, strategies and funding levels. Tsunamis can be triggered by earthquakes, subaerial or submarine landslides, volcanic eruptions, glacial calving, near-earth objects, weather or other events. These events can cause severe destruction, injuries, and loss of life due to powerful currents and flooding. Tsunamis pose a substantial threat to the western United States and all U.S. territories, as described below. ■ Hawaii is threatened by distant tsunamis due to its central location in the Pacific Ocean basin and has a history of local events. ■ Alaska, particularly the Aleutian Islands, faces local tsunami threats due to proximity to the Alaska-Aleutian Subduction Zone, as well as distant tsunamis from around the Pacific Ocean basin. ■ The western coast of the U.S. is threatened by distant tsunamis from around the Pacific Ocean basin and local source tsunamis from earthquakes generated within the Cascadia Subduction Zone in the Pacific Northwest. ■ American Samoa faces local tsunami threats from earthquakes generated in the nearby Tonga Trench, as well as distant tsunami threats. ■ Guam and the Commonwealth of the Northern Mariana Islands are threatened by local tsunamis from the nearby Mariana Subduction Zone, as well as distant sources from around the Pacific Ocean Basin. ■ Puerto Rico and the United States Virgin Islands are threatened by multiple local and distant tsunami sources, such as the Puerto Rico Trench (PRT), given their location in the complex seismic region of the Caribbean Sea. Several historical events stand out because of their catastrophic impacts. ■ In the Pacific Northwest, the 1700 Cascadia earthquake caused a tsunami that affected coastal Native American communities, though the extent of the damage is not fully documented (Ludwin, et al., 2005). ■ In Puerto Rico, the 1918 earthquake triggered a tsunami that caused $77 million in damage in 2022 dollars and 116 fatalities, primarily along the western coast (Coffman et al., 1982). ■ The 1946 Aleutian Islands earthquake triggered a massive tsunami that devastated Hilo, Hawaii, killing 158 people and resulting in approximately $375 million in damage (adjusted to 2022 dollars) (Fisher et al., 2023). ■ The 1964 Alaska earthquake (M 9.2) generated tsunamis that caused severe destruction in some communities across Alaska, Oregon, and California. This disaster led to a total of 124 fatalities and approximately $2.9 billion in property damage (adjusted to 2022 dollars) (Brocher et al., 2014) (Alaska Science Center, 2024). ■ In American Samoa, a tsunami generated by the 2009 Samoa earthquake (Mw 8.1) caused widespread devastation, resulting in 34 confirmed fatalities (Apatu et al., 2013) and economic losses exceeding $160 million (adjusted to 2022 dollars) (DHS, 2011). More recent events, including the 2010 Chile earthquake, the 2011 Japan earthquake, and the 2022 Tonga volcanic eruption, resulted in millions of dollars in damage to numerous ports and harbors in the U.S. South Pacific territories, Hawaii, and along the west coast of the U.S. (Lynett, et al., 2022) (Wilson, et al., 2013). Since these events, the expansion of the built environment in lowlying areas along the coast has increased the exposure of buildings and people, thereby further escalating community risk from tsunamis. This report provides a comprehensive national assessment of earthquake-generated tsunami risk. It does not include impacts from tsunamis generated by landslides, volcanic eruptions, glacial calving, near-earth objects, weather, or other events. This study is based on the best available hazard data from the U.S. Pacific Coast (California, Oregon and Washington), Alaska, Hawaii, U.S. Pacific Territories (American Samoa, Guam and Commonwealth of the Northern Mariana Islands) and Caribbean Territories (Puerto Rico and United States Virgin Islands). Tsunami risks associated with states along the East Coast, Gulf Coast, and Great Lakes are not included in this study because Hazus 6.1 software (FEMA 2024a) does not currently include the ability to analyze tsunami risk in those states. Once modeling capabilities and tsunami hazard data become available for additional states, FEMA may incorporate these data into future editions of this study.

Alaska, California, Hawaii Oregon, Washington↗
USGS2025

Do Graviquakes exist?

The “Graviquake” model, proposed in 2015 as an alternative to the elastic dislocation model, posits that normal faults are passive features dominated by coseismic gravitational collapse into a dilated crustal wedge, and that normal faulting is fundamentally distinct from strike‐slip and reverse faulting. Developed using finite‐element modeling before the 2016 central Apennines earthquake sequence, the model was revamped based on interpreted Differential Interferometric Synthetic Aperture Radar data from these events and used as evidence for a gravitational collapse episode. However, this interpretation relies on miscalculated elevation changes and is not corroborated by independent geophysical and seismological observations. Our analysis exposes fundamental flaws in the Graviquake model. By assuming that faults are passive players, it underrepresents the dynamic role of strain accumulation and release in rocks adjacent to faults. The hypothesized rapid expulsion of overpressurized fluids appears inconsistent with observed diffusion rates and lacks supporting seismological evidence. Part of the uplifted–subsided volume imbalance is likely an artifact arising from data processing, and in part is a transient effect due to the delayed response of the lower crust. Moment tensor analyses detect no isotropic components indicative of gravitational collapse, and observed ground motion and stress‐drop levels remain fully consistent with elastic dislocation theory. In addition, finite‐element modeling of normal faulting replicates observed surface deformation without invoking a collapsing wedge. The Graviquake model proposes a representation of normal‐faulting mechanics that differs significantly from established models and observations. Gravity does play a role in normal faulting, but the elastic dislocation theory remains the definitive framework of fault mechanics. Reinterpreting the 2016 earthquakes as a cascade of gravitational episodes, based on incorrect data processing and modeling, fails to substantiate the Graviquake hypothesis. Persistence in advocating this model could mislead seismic hazard assessment and undermine our understanding of normal faulting.

Bulletin of the Seismological Society of America↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗
USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

New Mexico↗
USGS2025

Site response and wave propagation effects in the eastern United States

Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.

eastern United States↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS2026

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2025

Rare earth element-mineralized carbonatite in the Bear Lodge Alkaline Complex, USA—Ore genesis implications from fluid inclusion characterization

Rare earth element (REE) resources of the Bear Lodge Alkaline Complex, Wyoming, are hosted in variably leached carbonatite spatially related to diatreme breccia pipes. We investigated the genesis of REE and lesser-known gold resources through fluid inclusion analysis of carbonatite, fluorite breccia, and smoky quartz vein samples. Physicochemical characteristics of inclusion-trapped fluids were evaluated using petrography, microthermometry, Raman spectroscopy, decrepitate mound analysis, energy-dispersive spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and noble gas isotope analysis. Microthermometry results reveal three fluid types that affected carbonatite dikes within deeper zones that escaped near-surface, ore-grade REE enrichment: (1) high-temperature (330–432°C) magmatic fluid captured in fine-grained calcite; (2) REE-enriched alkali bicarbonate-sulfate brine; and (3) low-temperature (117–182°C) diluted magmatic or meteoric water. Multiphase brine-like inclusions contain burbankite, nahcolite, strontianite, celestine and alkali sulfate daughter crystals, linking them to early burbankite mineralization. Peripheral smoky quartz and fluorite occurrences at Smith Ridge, 1.5 km from the central carbonatite dike swarm, contain primary inclusions that are Cl-poor and rich in Na-HCO 3 -SO 4 , similar to secondary and pseudosecondary inclusions in carbonatites. Helium isotopes reveal a MORB-like source for carbonatite samples and an older crust signature at Smith Ridge, consistent with the proximal ridge-top exposures of Archean granite. Results from this fluid inclusion study coupled with previous studies of carbonatite mineral paragenesis, show that light REEs (LREEs) were not mobilized great distances. Instead, burbankite crystallized within carbonatite from alkali bicarbonate fluids. With sodium retained in early burbankite, outward-emanating fluids enriched in potassium relative to sodium (higher K:Na) resulted in potassium–ferric iron metasomatism of silicate host rocks. This alkali fractionation was accompanied by fractionation of LREEs and heavy REEs (HREEs), with LREEs dominating the central carbonatite resources. In contrast, areas of peripheral REE mineralization at Bear Lodge are commonly characterized by higher HREE:LREE ratios. The K:Na ratio of associated fenites or alteration assemblages could be indicative of early crystallized burbankite in carbonatites and REE fractionation processes potentially leading to areas of concentrated HREEs with greater supply vulnerabilities.

South Dakota, Wyoming↗
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↗
USGS2026

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America↗
USGS2026

Detecting earthquakes in noisy real-time GNSS data with deep learning for improved PGD magnitude estimation

To disseminate accurate and useful warnings, earthquake early warning (EEW) systems must quickly determine the size and location of an earthquake to estimate expected shaking. Traditional seismic‐based algorithms tend to underestimate the true magnitudes of large earthquakes, a phenomenon known as magnitude saturation. This limitation motivated the recent inclusion of Global Navigation Satellite Systems (GNSS) data into the U.S. Geological Survey’s ShakeAlert EEW system with the Geodetic First Approximation of Size and Time (GFAST) algorithm because GNSS data do not saturate with large ground motions. However, the noise levels of GNSS data are very high compared with traditional seismic data, which obscures P ‐wave arrivals and can result in less accurate magnitude estimations if displacement amplitudes are low, such as for lower magnitude earthquakes or large source–station distances. In this study, we develop a deep‐learning model that detects earthquakes in GNSS data and use the Ridgecrest, California, earthquake sequence as a case study to demonstrate how the model could act as a filter to reduce the amount of low‐quality data that enters an algorithm like GFAST. To preserve our limited real earthquake data for model inference, we generated a training dataset composed of >700,000 synthetic displacement waveforms. We combined the synthetic waveforms with real‐time GNSS noise to produce realistically noisy training waveforms and then tested our model on additional synthetic data and performed inference using the real data that were held back. We discuss the performance of our trained model on both the unseen synthetic data and real inference data. Our model can be used to selectively filter only high‐quality data where an earthquake signal is observed for input into an algorithm like GFAST (outperforming a simple signal‐to‐noise ratio–based filter) to reduce the error in GFAST’s real‐time earthquake magnitude estimations.

California↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2025

An empirical Green’s function approach for isolating directivity effects in earthquake ground-motion amplitudes

In this study, we apply an empirical Green’s function (eGf) method within a ground‐motion modeling framework to mitigate trade‐offs between source, path, and site effects. Many physical processes contribute to spatial variations in observed ground motions, including earthquake radiation pattern, directivity, variable path attenuation, and site effects. Current nonergodic ground‐motion models use spatially varying coefficients for path and site effects, but they do not address trade‐offs with complex earthquake source effects. To quantify the influence of directivity on ground‐motion amplitudes, we use records from multiple smaller earthquakes with epicenters near that of a larger event. We use these small magnitude events as eGfs and estimate repeatable path and site effects at individual stations, assuming that the average adjustments are not controlled by directivity. We adjust residuals from the larger earthquake using the eGf terms, isolating effects related to the rupture. This method clearly enhances the observed broadband directivity observed in the 2022 M 5.1 and 2007 M 5.4 Alum Rock earthquake ground motions, reinforcing the conclusion that their ruptures were unilateral. For the 2004 M 6.0 Parkfield earthquake, we find a bilateral rupture model better fits the data because variations in rupture velocity, slip rate, and slip distribution seem to have a stronger effect on the ground motions than rupture direction alone. Applying eGf adjustments reduces the standard deviation of the rupture models over the three earthquakes by 32% on average and by up to 57% for the 2022 Alum Rock earthquake, confirming we have effectively removed repeatable effects related to the wave propagation path and site response. We propose a novel measure of the frequency‐dependent directivity amplification strength as the reduction in ground‐motion residual variability gained by fitting a directivity model; for the three earthquakes considered, this parameter varies between 25% and 75%, indicating that directivity can strongly influence ground motions and should be considered in ground‐motion modeling.

California↗