Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Article”

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.

979 records · Page 21Linked to original sources
USGS2024

A multidisciplinary approach that considers occurrence, geochemistry, bioavailability, and toxicity to prioritize critical minerals for environmental research

Critical minerals (or critical elements) are minerals or elements that are essential to global security and development and have supply chains vulnerable to disruption. In general, knowledge of the environmental behavior and health effects of critical elements is needed to support the development of safe and environmentally responsible supplies. This knowledge includes identifying potential consequences of increased critical element production and use, alternative critical element sources such as mine wastes, and adverse effects of critical elements on ecosystem condition and organismal health. Here we identify significant data gaps in the understanding of critical elements in surficial and aquatic environments, and the need, given the large number of commodities (50) identified on the 2022 critical minerals list for the United States, for an approach to prioritize them for study of their environmental fate and effects. We propose a multidisciplinary approach for this prioritization, considering measures of occurrence, geochemistry, bioavailability, and toxicity. We describe relatively easy-to-obtain metrics for each of these topic areas and demonstrate the utility of this integrated prioritization approach using indium and zinc as examples. This approach facilitates prioritizing research with a focus on those critical elements that are most mobile in the environment, bioavailable, toxic, or simply lacking data in these categories.

Environmental Science & Technology↗
USGS2025

Grand Canyon landslide-dam and paleolake triggered by the Meteor Crater impact at 56 ka

This paper hypothesizes that the Meteor Crater impact in Arizona, USA, 56,000 years ago triggered landslides in Grand Canyon that dammed the Colorado River and formed Nankoweap paleolake. This is compatible with shock and earthquake physics for the impact that infer a M5.4 seismic event, attenuated to an effective magnitude of M3.5 at Grand Canyon. Results that support the hypothesis include radiocarbon dating of driftwood and luminescence dating of associated slack-water lake sediments that are preserved in caves up to 60 m above the modern Colorado River. Radiocarbon ages from two locations, including Stanton’s Cave, date the driftwood as 55.25 ± 2.44 ka (n = 4). Sediments associated with the driftwood gave a luminescence age of 56.00 ± 6.39 ka (n = 2). These six Grand Canyon dates, and three published ages for the Meteor Crater impact, show statistically indistinguishable results that support the hypothesis for a geologically instantaneous series of events with a mean age of 55.60 ± 1.30 ka. This work highlights the value of radiocarbon dating near the limits of the technique, integration of multiple dating methods, and seismic and landslide hazards associated with meteorite impacts in regions of extreme topography like Grand Canyon.

Arizona↗
USGS2026

Per- and polyfluoroalkyl substances in waters associated with oil and gas development in the Denver Basin

Use of per- and polyfluoroalkyl substances (PFAS) in the petroleum industry could be a cause for concern due to the large volumes of produced water (PW) generated during oil and gas extraction, the reuse of these wastes in water-stressed regions, and adverse health outcomes related to PFAS exposures. However, PW PFAS characterization is nearly absent in the literature, and hydraulic fracturing (HF) chemical disclosures often omit the identities of additives as proprietary. Here we evaluate PFAS in PW samples from three petroleum wells in the Denver Basin during their first year of production. Total concentrations of targeted PFAS (Σ 40 PFAS) were < 35 ng/L in PW samples, with short-chain PFAS like perfluorobutanoic acid persisting throughout the sampled duration. Analysis of freshwater inputs for hydraulic fracturing (Σ 40 PFAS ~ 113 ng/L) and mixed fracture fluid (Σ 40 PFAS ~ 69 ng/L) indicated much of the targeted PFAS content was derived from the input water, and not from HF additives, however samples subjected to oxidation indicated the presence of PFAS precursors that would not be detected by targeted analysis. This study highlights that while PFAS content is low in the studied PWs, the potential for redistribution of PFAS in the environment may be a consideration for reuse applications.

Colorado, Nebraska, Wyoming↗
USGS2025

Cave records reveal recent origin of North America’s deepest canyon

We explore how and when Hells Canyon, North America’s deepest river gorge (~2,400 m deep), formed, addressing these fundamental questions first posed by W. Lindgren [ The Gold Belt of the Blue Mountains of Oregon (1901)]. Existing hypotheses about the canyon’s formation and timing of incision remain speculative due to a lack of direct constraints and geomorphic analysis in the canyon. Herein, we combine cosmogenic nuclide dating of cave-bound river deposits, river profile analysis, and numerical modeling to provide the first direct age constraints and systematic analysis of incision processes at work in Hells Canyon. Our study reveals a significant drainage capture triggered rapid incision at ~2.1 ± 1.0 Ma, establishing the Snake River’s modern route into the Columbia River system. The increased drainage area and subsequent increase in stream power resulted in the rapid incision of Hells Canyon and the formation of tributary knickpoints (KPs) that decrease in elevation away from the capture location. Cosmogenic dating of cave deposits indicates incision rates increased from ~0.01 to ~0.16 mm y −1 . Numerical modeling of the stream capture supports these observations, demonstrating how abrupt drainage area increase drives rapid river incision. Our findings from Hells Canyon provide a well-constrained example of how drainage capture can dramatically shape the evolution of a major river gorge.

Idaho, Oregon, Washington↗
USGS2026

Landslide volume estimation from seismic waveform features using a global catalog of seismogenic events

Seismic waves generated by large, rapid landslides encode information about the source and can be analyzed rapidly following an event. Even remote landslides can pose hazards to downstream communities, so rapid detection and characterization using existing seismic monitoring networks could be beneficial. In this study, we expand on past regionally limited work by presenting a globally applicable method for estimating landslide volume from seismic features that could be integrated into future landslide seismic monitoring frameworks. We train the model using multivariable linear regression and five seismic features derived from recordings of 129 landslide events of a range of styles and locations with independently estimated volumes. We present two preferred models, one that combines long‐period (LP) and high‐frequency (HF) features and one for use on smaller landslides without observed long‐period signals. We find that our best‐performing model, applicable to landslides larger than 100,000 m 3 with signals containing observable long‐period energy, requires only two features: LP (20–100 s) absolute maximum amplitude and HF (1–5 Hz) rise time (time between the signal start and the maximum envelope amplitude) and has an R 2 score of 0.79. This model predicted volumes within one order of magnitude for 55 out of 58 events. We find that this combination mitigates the trade‐off between mass and acceleration and other variations of landslide style that limit methods based on amplitude alone.

Seismological Research Letters↗
USGS2025

Kiloyear cycles of carbonate and Mg-silicate replacement at Von Damm hydrothermal vent field

The Von Damm vent field (VDVF) on the Mid-Cayman Rise in the Caribbean Sea is unique among modern hydrothermal systems in that the chimneys and mounds are almost entirely composed of talc. We analyzed samples collected in 2020 and report that in addition to disordered talc of variable crystallinity, carbonates are a major class of mineral at VDVF. The carbonate minerals include aragonite, calcite, magnesium-rich calcite, and dolomite. Talc and carbonate mineral textures indicate that, rather than replacing volcanic host rock, they precipitate from the mixing of hydrothermal fluids and seawater at the seafloor, occurring in chimneys and surrounding rubble. Alternating precipitation of this mineral assemblage is pervasive, with carbonate minerals typically being succeeded by talc, and with indications that in some cases talc and carbonate minerals replace one another. Stable carbon isotopic data indicate the carbonate minerals originate from the mixing of seawater and hydrothermal fluid, which is supported by U-Th data. Radiocarbon calcite ages and talc 234 U- 230 Th isochron ages indicate mineral ages spanning over thousands to tens of thousands of years. Analyses of these samples illustrate a dynamic system that transitions from carbonate-dominated to Mg-silicate−dominated precipitation over time scales of thousands of years. Our observations raise questions regarding the eventual fate of seafloor precipitates and whether carbonate and silicate minerals in such settings are sequestered and represented in the rock record.

Geology↗
USGS2026

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States↗
USGS2026

Rewilding the upper Klamath River Basin: Rapid recolonization of Chinook salmon following the world’s largest dam removal

Habitat fragmentation is a global driver of decline in biodiversity, yet evaluations of ecological response following re-connection remain rare. In 2024, four hydroelectric dams on the Klamath River were removed, reconnecting over 640 km of anadromous habitat and marking the largest dam removal project completed to date. We used a stationary high-resolution multi-beam imaging SONAR to quantify adult Chinook salmon ( Oncorhynchus tshawytscha ) passage into newly reopened reaches during the first two migratory seasons following dam removal. An estimated 7,742 (95% CI: 7,702–7,778) Chinook salmon migrated upstream in 2024 and 13,310 (95% CI: 12,876–13,733) in 2025, representing 18–19% of fish returning to the Klamath River Basin. Within this period, Chinook salmon recolonized 88% of their documented historical distribution, expanding into 345 km of reconnected habitats and reaching elevations of up to 1,250 m. Recolonization required fish to swim through dam removal sites, complex hydraulic corridors and engineered fishways, demonstrating the propensity for rapid recolonization at basin scale despite these challenges. These findings provide the first quantitative assessment of salmon response to the world’s largest dam removal and demonstrate how reconnection can catalyze large-scale recolonization in historically fragmented river networks.

California, Oregon↗
USGS2024

Potential climate and human water-use effects on water-quality trends in a semiarid, western U.S. watershed: Fountain Creek, Colorado, USA

Nutrients, total dissolved solids (TDS), and trace elements affect the suitability of water for human and natural needs. Here, trends in such water-quality constituents are analyzed for 1999–2022 for eight nested monitoring sites in the 24,000 km 2 Fountain Creek watershed in Colorado, USA, by using the weighted regressions on time, discharge, and season (WRTDS) methodology. Fountain Creek shares characteristics with other western U.S. watersheds: (1) an expanding but more water-efficient population, (2) a heavy reliance on imported water, (3) a semiarid climate trending towards warmer and drier conditions, and (4) shifts of water from agricultural to municipal uses. The WRTDS analysis found both upward and downward trends in the concentrations of nutrients that reflected possible shifts in effluent management, instream uptake, and water conservation by a watershed population that grew by about 40%. Selenium, other trace elements, and TDS can pose water-quality challenges downstream and their concentrations were found to have a downwards trend. Those trends could be driven by either a warming and drying of the local climate or decreased agricultural irrigation, as both would reduce recharge and subsequent mobilization from natural geologic sources via groundwater discharge. The patterns illustrate how changes in climate and water use may have affected water quality in Fountain Creek and demonstrate the patterns to look for in other western watersheds.

Colorado↗
USGS2025

Accurate simulation of flow through dipping aquifers with MODFLOW 6 using enhanced cell connectivity

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.

Groundwater↗
USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗
USGS2026

From landslide susceptibility to risk assessment in the conterminous U.S.

Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.

Engineering Geology↗
USGS2008

Late Jurassic igneous rocks in south-central Arizona and north-central Sonora: Magmatic accompaniment of crustal extension

In Middle Jurassic time, the region of south-central Arizona and north-central Sonora was part of a continental magmatic arc. In early Late Jurassic time, this magmatic arc gradually gave way to a extensional regime. Late Jurassic rifting was accompanied by magmatism, typically small in volume relative to the preceding arc phase, and with mildly alkaline compositions. These igneous rocks, the 158- to 146-Ma Ko Vaya Suite, make up much or most of the Comobabi, Artesa, Quijotoa, and Brownell Mountains and Sierra del Cobre. The bimodal Ko Vaya Suite constitutes volcanic to shallow plutonic complexes comprising chiefly trachyandesite and volcanic wacke; monzodiorite; rhyolite porphyry; distinctive, compositionally and texturally heterogeneous, quartzpoor granite (the Ko Vaya Granite); A-type perthite granite and quartz syenite; aphyric intrusive rhyolite; and locally abundant hematite veins. Alteration, particularly potassic alteration, is widespread. Minimally altered samples of granite and quartz syenite have the composition of trachyte or trachydacite and alkali rhyolite. The Ko Vaya Suite is nonconformably overlain by the latest Jurassic Sand Wells Formation, approximately correlative with the Glance Conglomerate, the basal unit of the Bisbee Group. The Ko Vaya Suite is in some ways analogous to the coeval Independence dike swarm of eastern California, but less diffuse. We infer that the Ko Vaya Suite represents especially intense or focused Late Jurassic extensional magmatism.

Arizona, Sonora↗
USGS2026

From model to action: Identifying the gaps in coastal marsh models for decision making

Coastal salt marsh modeling projects are designed to answer questions about where salt marshes may migrate, how they can migrate, and how they will be impacted by changing water levels. While this information is vital for climate-resilient marsh conservation, translating these findings into actionable results remains a challenge. The authors examined the actionability of marsh modeling efforts in the United States by analyzing data collected from semi-structured interviews of marsh model users ( n = 24) across two ongoing projects along the U.S. East Coast and the Gulf of America (Gulf of Mexico). By qualitatively analyzing the interview data, the authors found that the tasks of practitioners who use coastal marsh model outputs fall into three major themes: (1) marsh restoration, (2) planning with uncertainty, and (3) conserving habitat for marsh-reliant species. For each of these themes, the authors identify unmet needs, including high spatial resolution information for local planning, accessible descriptions of uncertainty to increase user confidence, and the incorporation of human dimensions data (e.g., human alterations to the landscape such as impoundments and culverts) for a comprehensive understanding of the coastal salt marsh. Marsh modeling projects should strategize how to fulfil these unmet user needs so that marsh modeling outputs can better support decisions related to marsh conservation and restoration.

Gulf of America, United States East Coast↗
USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2024

A new water temperature modeling approach to predict thermal habitat suitability for nonnative cichlids in Florida rivers

As global temperatures increase, the spatiotemporal arrangement of thermal habitats in Florida rivers may shift, creating the potential for greater dispersal and establishment of nonnative tropical freshwater fishes. To understand how water temperature changes may affect the spatial distribution of these nonnative species, more effective water temperature prediction models are necessary. Currently, most models employ either a generalized air–water temperature relationship or require expensive and complicated tools to measure hydrometeorological factors (e.g. groundwater input). Thus, we developed a novel modeling approach that is accurate, accessible, and cost-effective in allowing fisheries managers to project water temperatures in rivers across Central and North Florida. To characterize the potential for nonnative fishes to spread northward, we evaluated two hardy and abundant species currently found primarily in South Florida: Mayan Cichlid ( Mayaheros urophthalmus ) and Oscar ( Astronotus ocellatus ). Our results show an increase in thermally suitable winter days for both species in 10 of 11 rivers studied, consistent with predicted water temperature warming under 16 climate-change scenarios spanning different levels of air temperature warming (+1 °C, +2 °C, +3 °C, +4 °C) and precipitation/groundwater thermal sensitivity (0, 0.33, 0.66, 1). Considering resource limitations, fisheries managers can use our water temperature modeling approach to predict effects of climate change on Mayan Cichlid and Oscar survival, growth, and dispersal and take actions to manage potential northward movement of these species.

Florida↗
USGS2025

Gaps in water quality modeling of hydrologic systems

This review assesses gaps in water quality modeling, emphasizing opportunities to improve next-generation models that are essential for managing water quality and are integral to meeting goals of scientific and management agencies. In particular, this paper identifies gaps in water quality modeling capabilities that, if addressed, could support assessments, projections, and evaluations of management alternatives to support ecosystem health and human beneficial use of water resources. It covers surface water and groundwater quality modeling, dealing with a broad suite of physical, biogeochemical, and anthropogenic drivers. Modeling capabilities for six constituents (or constituent categories) are explored: water temperature, salinity, nutrients, sediment, geogenic constituents, and contaminants of emerging concern. Each constituent was followed through the coupled atmospheric-hydrologic-human system, with prominent modeling gaps described for a diverse array of relevant inputs, processes, and human activities. Commonly identified modeling gaps primarily fall under three types: (1) model gaps, (2) data gaps, and (3) process understanding gaps. In addition to potential solutions for addressing specific individual modeling limitations, some broad approaches (e.g., enhanced data collection and compilation, machine learning, reduced-complexity modeling) are discussed as ways forward for tackling multiple gaps. This gap analysis establishes a framework of diverse approaches that may support improved process representation, scale, and accuracy of models for a wide range of water quality issues.

Water↗
USGS2026

Predictions of anthropogenic background PFAS concentrations in soil and relation to bedrock lithology and groundwater quality

Detectable concentrations of per- and polyfluoroalkyl substances (PFAS) have been observed in soils in remote areas and presumably originate from atmospheric deposition. These anthropogenic background concentrations may enable some PFAS to leach to groundwater at levels that exceed regulatory criteria for drinking water. However, anthropogenic background soil concentrations and their connection to groundwater are not well characterized. We developed a boosted regression tree model to predict perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) concentrations in shallow soils across northern New England. Low soil pH was the most important predictor of elevated anthropogenic PFOS and PFOA concentrations in background soils, rather than potential PFAS sources, land use, or population density. Total organic carbon (TOC) was also an important predictor for PFOS soil concentrations. Model predictions indicate that 73% of the shallow soils within Maine, Vermont, and New Hampshire exceed New Hampshire’s Soil Remediation Standard for PFOS (0.5 ng/g) and 41% exceed the PFOA standard (0.4 ng/g). Analysis of soil model results and groundwater data suggests that areas with high soil pH are associated with higher groundwater detection frequencies, illustrating how areas with less retention in soil are, conversely, also areas with potentially greater groundwater vulnerability. Further analysis indicates that groundwater may be more vulnerable in calcareous lithologies.

Maine, New Hampshire, Vermont↗