USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geochronology, Thermochronology and Time Scale Research”

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,663 records · Page 28Linked to original sources
USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗
USGS2024

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Geological Society, London, Special Publications↗
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2026

Permafrost thaw controls iron flux from wetlands and sulfide-bearing rocks to Arctic rivers and streams

Recent warming has caused widespread iron mobilization into Arctic waterbodies that degrades ecosystems and threatens natural resources. Yet, understanding where and when iron flux occurs remains limited. Here, we investigate iron loading across regional to local scales in Arctic Alaska using climate, water chemistry, and borehole data together with mapped geology and permafrost presence. We show that both anoxic microbial iron reduction and acid rock drainage from iron-sulfide oxidation mobilize iron. Iron influx is strongly associated with lowland wetlands, sulfide-rich upland bedrock, and near-surface permafrost. Acid rock drainage chemistry correlates very strongly with the depth of seasonal thaw above permafrost from the previous year, indicating a one-year lag. These findings clarify the spatial and temporal dynamics of Arctic river rusting, provide a mechanistic understanding of the phenomenon, and may allow anticipation of its occurrence and assessment of its implications for aquatic ecosystem health and subsistence resources under ongoing climate change.

Alaska↗
USGS2025

Divide and conquer: Separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process – the probability of the existence of a seismic phase (detection probability) and the probability of correctly identifying the phase arrival time (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labeled phase arrival time, and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks) and the height of the peaks are taken as “probability” of the picks. Here, we show that this “probability” represents neither the detection nor the timing probabilty because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multi-class classifier to obtain a probability mass function to address the timing probability. This new approach makes the deep learning-based phase picking process more interpretable and gives us options to logically control seismic monitoring workflows.

Geophysical Journal International↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS2026

Tringa flavipes (Lesser Yellowlegs) from separate breeding sites subdivides the Prairie Pothole Region in space and time during southbound migration

Some staging regions support multiple groups of the same migratory species, each of which may use the region differently. Characterizing the ways, in which separate groups use such regions can therefore help to identify vulnerabilities during this sensitive period of the annual cycle. The Prairie Pothole Region (PPR) is a massive wetland complex in the northern Great Plains of North America used by ∼11 million shorebirds during migration. The PPR has been heavily modified by agriculture and is experiencing varied effects of global climate change, threatening the health of the shorebirds that rely on it. Here, we used 6 seasons of southbound tracking data of Tringa flavipes (Lesser Yellowlegs)—a long-distance migratory shorebird species with an estimated population decline of 63% over the last 4 decades—from 9 sites across their breeding range to explore differences in migratory behavior within this important staging region. We found that 75% of tracked individuals used the region during southbound migration, and T. flavipes from different breeding sites detoured 110–875 km from their most direct migratory route to access the PPR. Individuals that arrived later stayed longer and made more stops within the region than those that arrived early. Individuals originating from different breeding sites also displayed spatial and temporal segregation within the region: T. flavipes from southwest and central Alaska relied heavily on the northwestern PPR, while those from Canada used the central and southeastern portions of the PPR. Finally, timing of use varied among groups, but the southeastern PPR became increasingly important over the course of the southbound migratory window, as other wetlands likely dried out. Our study highlights the portions of the PPR of critical importance to migrating T. flavipes and the diversity of ways, in which different groups from within the same species can use a single staging region.

Alaska, Manitoba, Northwest Territories, Ontario, ↗
USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗
USGS2026

Satellite time series analysis to quantify changing climax ciénegas using a state and transition model approach

Ciénegas are rare wetlands in arid landscapes of the North American Southwest, historically providing critical ecological and hydrological functions but increasingly threatened by changing climate and land use pressures. This study quantifies changes in ciénega condition and floodplain dynamics using a state-and-transition model (STM) informed by expert knowledge and remote sensing. Key factors include woody plant encroachment, water availability, and soil aggradation. We mapped 31 ciénegas with high-resolution imagery and analyzed Landsat data (1985–2023) to assess vegetation health and moisture using the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Infrared Index (NDII). Results show substantial interannual variability in phenology, water stress, and soil moisture, with regional drying and elevation strongly influencing ciénega resilience. We classified ciénegas into three functional states—healthy, desiccated, and dormant—and mapped their 2023 condition. Trend analyses indicate most ciénegas exhibit greening despite drought, though localized variability underscores the need for site-specific management. None are in a stable climax (reference) state; rather, they transition among states in response to external drivers. Increasing woody plant cover and surface drying, likely linked to declining regional water tables, favor deep-rooted species over wetland grasses—a pattern mirrored in adjacent control plots. Spatially explicit analysis revealed intra-ciénega variability often masked by aggregated data, highlighting the importance of high-resolution monitoring. Seasonal and long-term trends provide context for understanding ciénega dynamics, including degradation and restoration pathways. This study emphasizes the importance of groundwater conservation and demonstrates how remote sensing supports long-term monitoring. The STM framework offers a practical tool for adaptive management to sustain freshwater resources in arid environments.

Arizona, New Mexico↗
USGS2025

Aftershocks in stress shadows are inconsistent with modeled static Coulomb stress changes

Aftershock triggering is commonly attributed to increases in static Coulomb stress. In some areas, termed "stress shadows", a decrease in Coulomb stress is predicted to suppress earthquake occurrence. However, aftershocks are often observed in the modeled stress shadows. We examine several hypotheses that attempt to reconcile these shadow aftershocks with the static Coulomb stress change model: (1) they appear to be in shadows because of inaccuracy in the stress change calculations, (2) they occur on faults of unusual orientation which actually experienced increased Coulomb stress, (3) they occur on faults with different frictional properties, not modeled well by Coulomb stress, and (4) they are secondary aftershocks triggered by prior aftershocks or afterslip. When tested on the 2016 Mw7.0 Kumamoto, Japan, and 2019 Mw7.1 Ridgecrest, California, aftershock sequences, none of these hypotheses can explain the majority of the shadow aftershocks, and taken together these hypotheses can explain only about half of these aftershocks. This implies that Coulomb stress modeling that lacks small-scale fault zone heterogeneity might be inadequate to fully capture the true static stress changes and/or that other physical triggering models are needed, for example transient processes such as delayed triggering by dynamic stress changes from the passing seismic waves.

California↗
USGS2026

A diatom-based quantitative sea-ice proxy for the Bering and Chukchi seas

Sea ice affects Earth's climate system on both regional and global scales. Its incorporation into climate can be used to achieve more accurate predictions of future climate. However, instrumental records of sea-ice concentration do not extend earlier than 1978. In an effort to extend this record, we constructed a proxy using the generalized additive model based on relative abundances of five easy-to-identify diatom species found in sediment samples across the Bering and Chukchi seas. Here we present the first quantitative diatom-based sea-ice proxy developed for Beringia. The developed proxy has been applied to two sediment cores in the Bering Sea ranging from 0 to 25.7 ka (HLY0204 51JPC) and 369 to 430 ka (IODP Exp 323 Site U1345) and one in the Chukchi Sea ranging from 2.7 to 10 ka (HLY0204 24JPC). The obtained reconstructions of sea-ice concentrations are similar, but not identical to previously published qualitative and nearby records based on other proxies. Because our results are quantitative, they can be incorporated into regional climate models. The proxy is publicly available as an R Shiny application (app) and can be applied to any diatom count from marine sediments in the region.

Bering Sea, Chukchi Sea↗
USGS2026

Keeping pace: Migrating sea lampreys adjust swim speed in response to shifting water velocity to maintain a constant ground speed

Upstream migrating fish must balance energetic costs with environmental and ecological constraints. Theoretical models predict that fish can conserve energy by adjusting swim speed in response to changing water velocity, but detailed field tests of these predictions remain rare. We used fine-scale acoustic telemetry and hydrodynamic modeling to track adult sea lamprey (Petromyzon marinus) migrating through a river and evaluated how movement speeds varied with environmental conditions and individual traits. Swim speed increased with faster water velocities, while ground speed remained constant, supporting the hypothesis that lampreys regulate swim effort to maintain a fixed migration pace. The average ground speed aligned with a predicted optimum of ∼1 body length per second. Both swim and ground speeds declined slightly in deeper water, suggesting a trade-off between movement efficiency and predator avoidance. These findings provide field-based support for a theoretical model of energy optimization and underscore how migratory fish adaptively adjust behavior to local hydrodynamic conditions. Our results highlight the value of using high-resolution tracking and environmental modeling to reveal behavioral tactics that underlie efficient migration strategies.

Michigan↗
USGS2025

Introduction to the special section on improving measurements of earthquake source parameters

Earthquake source parameters such as magnitude, seismic moment, source dimension, stress drop, and radiated energy are fundamental to understanding earthquake physics, and are also key ingredients in earthquake ground‐motion modeling, rupture simulation, and statistical seismology. However, the uncertainties in these parameters estimated from the radiated seismic wavefield are large due to variability in approaches, including site and attenuation characterizations, and so estimates for an individual earthquake made by different studies can vary greatly. Estimating spectral source parameters remains a popular topic, due to a combination of their intrinsic importance and their apparent ease of measurement, but also a controversial one, due to many sources of variability and large uncertainty. The available methods coupled with necessary parameter choices and assumptions in the analysis make it challenging for researchers to apply methods or understand the reliability in results or reported source parameters. This Special Section on Improving Measurements of Earthquake Source Parameters showcasing comparisons between methods and studies seeks to alleviate some of these difficulties to help the community identify the important components and trade‐offs of decomposing recorded seismograms into their source, path, and site components.

Bulletin of the Seismological Society of America↗
USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗
USGS2026

Variability and consistency in wildfire susceptibility: Insights from a national compilation

Background Wildfire risk in the United States is rising and remains a land management priority. The quantitative wildfire risk assessment (QWRA) framework integrates fuels, topography, weather and values at risk to estimate the potential change in value from wildfire. Within this, response functions (RFs) represent how values respond to fire intensity. These are often based on expert judgment, but variation across assessments is unclear. Aims This study uses data from the US Geological Survey (USGS) Wildfire Hazard and Risk Assessment Clearinghouse to characterize consistency and variation across categories and contexts. Methods We applied descriptive statistics to summarize RFs, using tables, box-and-whisker plots and heat maps stratified by highly valued resource or asset (HVRA) category and spatial scale. Key results RFs and value definitions vary, especially for ecosystem-related resources. Some functions, such as for buildings in the wildland–urban interface (WUI), translate well across contexts, while others require more input. Conclusions Some functions are broadly transferable, while others need customization. This analysis provides references and starting points for improvement to RFs in QWRAs. Interpretations Expanding the clearinghouse and dataset and building more transparency in expert elicitation can build trust among communities, agencies and end-users, and can support efficient use of limited resources to mitigate wildfire risk.

International Journal of Wildland Fire↗
USGS2026

Groundwater spatial variability within an atoll island: Assessing shallow aquifer heterogeneity with geophysical and physicochemical measurements

This study examines the spatial variability of shallow groundwater on Dhigelaabadhoo Island using electromagnetic induction surveys, groundwater monitoring, and sediment analyses. The research reveals how variations in island morphology—such as differences in elevation, reef flat width, and sediment composition—affect the spatial distribution of groundwater lenses and the overall aquifer dynamics. Saltwater intrusion is especially pronounced in low elevated areas, with narrow reef flat plate and areas where higher hydraulic conductivity—driven by the presence of coarser sediments—is observed, whereas regions characterized by finer sediments, higher elevation, and wider reef flat plates tend to support more symmetric and less saline groundwater lenses. The geophysical investigations reveal that tidal oscillations alter groundwater movement by markedly changing water levels and conductivity, thereby underscoring the critical need to account for temporal dynamics in atoll coastal aquifer systems and the importance of integrating tidal dynamics into the aquifer zone. The findings highlight the significant role of intrinsic morphological and external hydrodynamic factors in shaping groundwater distribution on atoll islands, offering critical insights for sustainable freshwater resource management.

Dhigelaabadhoo Island↗
USGS2024

Taking heat (downstream): Simulating groundwater and thermal equilibrium controls on annual paired air–water temperature signal transport in headwater streams

Headwater stream temperature often exhibits spatial variation at the kilometer-scale, but the relative importance of the underlying hydrogeological processes and riverine perturbations remains poorly understood. In this study, we investigated the relative importance of groundwater (GW) and other processes on downstream annual stream temperature signal characteristics using deterministic heat budget model (HFLUX) scenarios within an idealized stream reach representative of mountainous forested conditions. We summarized annual stream thermal regimes from the relationship of paired sinusoidal air and water temperature signals (amplitude ratio, phase lag, and mean ratio). Results showed that downstream changes in annual temperature depended on the thermal gradient between water and the hypothetical equilibrium temperature (where all heat fluxes sum to zero). GW inflow, riparian shading, and the boundary input signal were the most significant factors affecting downstream annual water temperature signals, while flow volume and channel dimensions impacted how quickly annual temperature signals changed. Effects of GW were dominated by advective rather than conductive heat exchange processes, but conduction played a larger role when GW input was more spatially diffuse. Our results indicated several mechanisms by which local processes may affect stream thermal resilience to disturbances and can help guide management of wildfire and climate change.

Journal of Hydrology↗