USGS Science⌕ Search

FIND YOUR NEXT DISCOVERY

Results for “Geophysical Research Letters”

Original records, connected by a shared subject.

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

230 recordsLinked to original sources
USGS2026

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado↗
USGS2025

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California↗
USGS2025

S/P amplitude ratios with Distributed Acoustic Sensing and application to earthquake focal mechanisms

Distributed acoustic sensing (DAS), which transforms a fiber optic cable into an array of high frequency strainmeters, has the potential to help us characterize earthquakes with a dense sampling of measurements. While earthquake focal mechanisms are frequently determined using P -wave polarities and S / P amplitude ratios with inertial seismometers, the dense sampling of DAS over potentially large portions of the focal sphere can aid our solutions. Here, we consider ∼200 regional earthquakes and thousands of S / P measurements on a DAS cable colocated with a network of inertial seismometers near Arcata, California. We demonstrate the S / P ratio measurements made on DAS are similar to those made on the inertial seismometers and can be used to constrain focal mechanism solutions.

California↗
USGS2026

Decadal shifts in groundwater age detected by environmental tracers across California, USA

Groundwater age offers important insight into recharge, storage, and contamination risk. Although models predict age changes can be driven by pumping and climate variability, direct observational evidence remains limited. Here, we analyzed paired environmental tracer suites (tritium, carbon-14, and tritiogenic helium-3) collected a decade apart from 268 wells across California to assess the prevalence of groundwater age transience. Travel-time distribution models and statistical tests indicated age transience at 29% of sites, occurring most often in agricultural regions, such as the San Joaquin Valley and Southern Coast Ranges, where large carbon-14 changes coincided with substantial nitrate and chloride shifts. Sites with tritiogenic helium-3 data showed more frequent age transience, underscoring the value of multi-tracer data sets. These results provide the first regional evidence of widespread groundwater age change and a method for detecting changing water balances with implications for groundwater sustainability and water quality.

California↗
USGS2026

Asynchronous landslide seasonality across the United States

Mid-range landslide outlooks can facilitate weather-related landslide preparedness and disaster response planning, but seasonal landslide activity remains poorly quantified at continental scales. Leveraging >55,000 reported landslides from across the United States (U.S.), we used circular statistics to quantify landslide seasonality in 67 National Weather Service County Warning Areas (CWAs). We found regional differences in landslide season timing and duration, with transitions between domains variably corresponding to climate class or river basin. We assessed differences in seasonality by movement type for slides, flows, and falls, detecting apparent, but uncertain, differences between slide and fall seasonalities in 27 of 35 (77%) of CWAs with both types reported. In the Pacific Northwest, where long records exist, we found a credible shift toward a later mean landslide season in western Washington from 1990 to 2020, but no trend in western Oregon. Our results can provide emergency planners a resource to assess seasonal landslide probability nationwide.

Geophysical Research Letters↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS2025

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California↗
USGS2025

A northeast-dipping zone of low frequency earthquakes at the southern edge of Cascadia subduction

Tectonic tremor monitoring occasionally detects events in an anomalous zone in southern Cascadia, 50–100 km west of the main tremor band, near the expected southern edge of the subducting Gorda slab at the Mendocino triple junction. To investigate the geometry and temporal behavior of this tremor, we examine its constituent low-frequency earthquakes (LFEs) by developing 27 stacked LFE waveform templates that we use to detect events from 2018 to 2024. We then relocate LFE sources together with regional seismicity. We find that LFE hypocenters form a northeast-dipping alignment at 22–29 km depth, extending eastward from a zone of micro-earthquakes, ∼15 km south of the southern edge of Gorda slab seismicity. These LFE families exhibit small bursts of activity every few days. Considering the strong world-wide association of tremor and LFEs with high slip-rate, plate-bounding faults, we hypothesize these LFEs may demark the southern edge of Cascadia subduction.

California↗
USGS2026

Dynamic coupling between faulting, rifting and magmatism during 2021-2025 unrest on Reykjanes Peninsula, Iceland

Interactions among faulting, earthquakes, and eruptions are fundamental to plate tectonics and hazard forecasting yet rarely observed along mid-ocean ridges. On Iceland's Reykjanes Peninsula, seismotectonic–volcanic unrest resumed after nearly 800-year hiatus, providing an opportunity to observe these interactions during 2021–2025 activity. By integrating high-resolution seismicity, focal mechanisms, satellite geodesy, surface deformation, and eruption data, we document ∼4 m of total extension accommodated through 14 rifting episodes. The largest, in 2023, involved graben reactivation and diking, with seismic swarms and earthquake faulting that matched the surface ruptures, where strike-slip faulting preceded normal-faulting earthquakes and extension. The accrued extension was released by extension fractures triggered by magma accumulation. Long-term observations show no correlation between erupted magma volume, seismicity, and crustal extension. This highlights dynamic relation between rifting, faulting, and magmatism in transtensional settings and their implications for hazard assessment.

Reykjanes Peninsula↗
USGS2025

Improved prediction of postfire debris flows through rainfall anomaly maps

Predicting where runoff-generated debris flows might occur during rainfall on steep, recently burned terrain is challenging. Studies of mass-movement processes in unburned areas indicate that event locations are well-predicted by rainfall anomaly, R* , in which peak observed rainfall is normalized by local rainfall climatology. Here, we use remote and field methods to map debris flows triggered within the 2020 Dolan Fire burn area in coastal California, demonstrate that a short-duration R* metric predicts debris-flow occurrence more effectively than absolute peak intensity or longer-duration rainfall metrics, and show that incorporating an R* criterion into an existing debris-flow likelihood model can reduce false positive predictions and improve accuracy. We test R * at three other climatically distinct fires in California, demonstrating its utility for mapping likely debris-flow locations in different climates. We also consider how R* can benefit postfire debris-flow prediction given recent increases in climatological variability within individual burn perimeters.

Callifornia↗
USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗
USGS2025

Dating a medieval tsunami with uranium-series techniques on Caribbean corals

Uranium-series dates from coral boulders constrain the timing of a medieval tsunami from the Puerto Rico Trench. Previously reported evidence for this tsunami includes hundreds of coral boulders that came to rest hundreds of meters inland on Anegada, British Virgin Islands. New U-series dates on these coral boulders provide limiting dates for the tsunami. The narrowest limits were by dating interior bands of a coral that retains the hemispherical form of a living coral colony, and which include adjustments for the number of annual density band couplets between the dated samples and the boulder exteriors. By those limits, the tsunami dates between 1381 and 1391 CE, and likely occurred during summer or fall. The tsunami is important as the only reported sign that the eastern Puerto Rico Trench has produced a great earthquake. The dating may aid in defining the earthquake source and in communicating tsunami hazards.

British Virgin Islands, Puerto Rico↗
USGS2025

Unveiling coseismic deformation from differenced legacy aerial photography and modern lidar topography: The 1983 M6.9 Borah Peak earthquake, Idaho, USA

The 1983 M6.9 Borah Peak, Idaho, earthquake is one of the largest historical normal fault earthquakes in the western United States. We quantified meter-scale vertical change along the 35 km-long rupture using topographic differencing of 1966 aerial imagery and 2019 lidar-derived data. The initial differencing results are largely obscured by horizontal and vertical georeferencing errors and flight-line stripes. Our error corrections are designed to be insensitive to the coseismic deformation and reduced error by 50%. We calculated vertical separation and resolved a maximum of 2.02 ± 0.46 m at Doublespring Pass. Our vertical separation measurements are generally consistent with those from prior studies using field data and post-earthquake topographic data. However, the differencing measurements are a few decimeters lower than these prior measurements, indicating that differencing can isolate historical from prehistoric earthquake deformation. Our study demonstrates that revisiting historical earthquakes can provide new insights into the magnitude and patterns of coseismic deformation.

Idaho↗
USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗
USGS2026

The role of groundwater in contributing to surface water salinization in the Upper Colorado River Basin

Freshwater salinization impacts the availability of water for human use and ecosystem needs worldwide. It has been estimated that total dissolved solids (TDS) in the Colorado River Basin cause $350 million/year in damages and substantial resources are devoted to reducing TDS loading to streams. This study describes the development and application of coupled watershed models that enable TDS source tracking through the subsurface and across the landscape at a seasonal timestep for 35 years in the Upper Colorado River Basin. Results indicate that, on average, 75% of TDS loading to streams originates as baseflow, and 50% of loading is lagged in delivery by longer than one season. Snowmelt was identified as a dominant process controlling the transport of lagged TDS to streams. This approach informs when and where TDS mitigation efforts may be effective in a watershed that serves as a critical water supply for the southwestern United States.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2025

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗