USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science and Technology - Water”

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,230 records · Page 5Linked to original sources
USGS2026

Thermal detectability of subsurface water ice on Mars: A comparative analysis for the Subsurface Water Ice Mapping (SWIM) Project

We have developed a new global map of shallow ground ice distribution, SWIM23, based on Mars Global Surveyor Thermal Emission Spectrometer data and made systematic comparisons between this new map and two similar, previously developed data products. We have explored the origin of differences between the three ice maps by detailing technical and procedural differences in their development, by making global pixel-by-pixel comparisons, and by carrying out a series of one-dimensional thermal simulations to explore fundamental physical limitations of thermal ice-detection techniques. These efforts and the production of a composite thermal ice-consistency map supported integration of multiple geophysical data products relevant to ice detection in the upper meter of the Martian regolith into the larger Mars Subsurface Water Ice Mapping project. Our work also highlights fundamental physical limitations to thermal ice detection as a technique, particularly the rapid fall-off in ice detection sensitivity at depths >30 cm, which produces maximum uncertainty in the presence and depth of ice within regions preferred for potential human landing sites. A future Mars orbiter mission designed to detect ice and support crewed landing site selection in the midlatitude region should give payload priority to an instrument capable of probing the 1–5 m depth range (i.e., a high-frequency radar), over a next-generation thermal spectrometer, which is unlikely to offer clarity on ice table depths or lateral continuity of the ice table in the locations of highest interest.

Planetary Science Journal↗
USGS2025

Using probability difference to compare streamflow information of alternatives for efficient operation of monitoring networks

Efficient operation of streamflow monitoring networks requires investments in technology and labor that provide the greatest benefits from available resources. Economic analyses comparing the costs and benefits from different types of alternatives for monitoring have not been practical to implement. Streamflow information provides a generic measure of benefits that can be incorporated into operational decisions as an objective for monitoring networks. A methodology for comparing how accuracy, monitoring period, and monitoring instead of modeling affects streamflow information is developed from information-theoretic approaches for network design but contributes three novel features: (1) a probability-difference model for conditional probability of monotonically paired variables, (2) explicit discounting of unverified information that may exceed the accuracy of streamflow records, and (3) run analysis to account for non-stationarity in streamflow probabilities. Application of the methodology to the U.S. Geological Survey streamflow monitoring network indicates the value of monitoring period to reduce the uncertainty of streamflow probabilities and, thus, increase streamflow information. The methodology has important limitations, particularly for sites with non-perennial streamflow, but demonstrates that probability difference could be used to evaluate operational alternatives to increase the efficiency of monitoring networks.

PLOS Water↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
USGS2026

Baseflow and snowmelt sustained streamflow in the Upper Colorado River Basin, 1986-2020

The Upper Colorado River Basin (UCRB) faces substantial water availability limitations. Although most streamflow originates as snowmelt, the partitioning of snowmelt between surface runoff and groundwater recharge and subsequent groundwater discharge to streams is highly uncertain. On average, over half of the streamflow in the UCRB is estimated to originate from groundwater discharge to streams, highlighting the importance of baseflow in sustaining surface water. However, the historical patterns of baseflow and streamflow, along with their variability over space and time and their specific sources, remain unknown at the basin scale. This study addresses those gaps by characterizing the sources and transport pathways of both baseflow and streamflow in the UCRB at a seasonal timestep from 1986 to 2020, including the lagged delivery of subsurface water to streams beyond the current season, using coupled models of baseflow and streamflow. Between 1986 and 2020, on average 63% of UCRB streamflow originated from baseflow. About half of this baseflow took longer than one season to reach streams, and outside the snowmelt season, baseflow was the dominant source of streamflow. Snowmelt was a key source of both baseflow and streamflow. Current season snowmelt contributed 33% of streamflow via runoff, and 22% of the 29% of streamflow that originated as current season baseflow via subsurface flow to streams. Over the study period, baseflow index (BFI) declined in headwaters and increased at mid-elevations. Springtime increases in BFI demonstrate the increasingly important role baseflow plays in water supply. Identifying the sources, locations, and timing of water that contributed to the UCRB outlet can inform management of water resources in the basin.

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

Irrigated agriculture influences selenium levels in an endangered marsh bird

Selenium bioaccumulation in aquatic food webs poses risks to wildlife, particularly in wetlands receiving irrigation runoff. The Salton Sea, California’s largest lake, is primarily sustained by agricultural drainage. This drainage creates wetland habitat along the lakeshore that many bird species depend on, including the federally endangered Yuma Ridgway’s rail ( Rallus obsoletus yumanensis ). However, these marshes may pose an ecological trap – attracting rails despite high selenium exposure. We captured rails during the 2020–2023 breeding seasons and compared rail selenium levels within three types of marshes (fed with irrigation runoff, Colorado River water, or groundwater). We collected blood, breast feathers, and head feathers of rails in all three water sources for selenium comparisons. We tagged adult rails with GPS transmitters to locate nests and foraging locations where we collected eggshells, unhatched eggs, and prey. We assessed selenium exposure by collecting multiple prey species commonly eaten by rails in all three water sources. Selenium concentrations varied among sampling locations. Selenium concentrations in most sample types were predominately influenced by water source and marsh inflow velocity (sometimes in combination with marsh size). Distance to inflow, however, did not influence selenium concentrations in any sample type. Selenium concentrations were highest in agricultural-fed marshes compared to river-fed and spring-fed marshes. Increased marsh inflow velocities resulted in lower selenium concentrations. Given the risk of an ecological trap, our results suggest that supplementing wetlands with Colorado River water could mitigate selenium bioaccumulation in Yuma Ridgway’s rails.

California↗
USGS2025

Relations of groundwater quality to long-term surface disposal of produced water near the Midway-Sunset and Buena Vista Oil Fields, California, USA

Contamination of groundwater by oil-field fluids in proximity to oil and gas development has been an issue of concern to water users and regulators given long histories of development and legacy disposal practices. A robust set of geochemical tracers including petroleum hydrocarbon compounds, thermogenic gases, inorganic ion concentrations, stable isotopes, radioactive isotopes, and noble gases were used to assess if oil-field fluids mixed with groundwater near the Midway-Sunset and Buena Vista Oil Fields in California, USA. Results show evidence of mixing of oil-field fluids with groundwater within the study area from either anthropogenic or natural processes. Produced water plumes associated with modern surface disposal facilities, used since the late 1950s, extend up to 1.5 km and currently remain within the boundaries of the oil fields. Plumes associated with earlier routing of produced water down natural drainages and in large retention structures (Midway Basin and Sunset Basin) near the Buena Vista Lake Bed are present in groundwater east of the oil fields. Based on geochemical tracer evidence, aerial imagery, and aerial electromagnetic surveys, these legacy plumes reach the western portion of the Central Valley aquifer system, an important groundwater resource for agricultural and domestic supply. The legacy plume associated with Sunset Basin may further be detected downgradient in deeper groundwater beneath the southern extent of the Buena Vista Lake Bed based on the presence of thermogenic gases and petroleum hydrocarbon compounds.

California↗
USGS2025

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona↗
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

Multiple machine-learning estimation of groundwater levels and trends for the regional Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer provides irrigation, public, and domestic water supplies across the south-central United States. Declining groundwater levels require improved characterization of changing conditions. Traditional potentiometric-surface mapping does not use all available water-level data or quantify uncertainty. To address these limitations, we developed a data-driven multiple machine-learning (MML) framework delivered through two open-source R packages. The covMRVAgen1 software assembles covariates to 155,960 monthly groundwater levels from 57,695 wells; the mmlMRVAgen1 software trains Cubist and Random Forest models, blends them, and makes 1-kilometer gridded predictions of monthly potentiometric surfaces for the period January 1980–December 2022. The MML approach provides a methodological foundation for region-scale spatiotemporal groundwater prediction and uncertainty quantification, generating 90-percent prediction limits with appropriate empirical coverage. Model performance is acceptable, with a root-mean-square error of about 4.2 feet, standard deviation of 24.82 feet, and a normalized Nash–Sutcliffe efficiency of 0.973.

Arkansas, Illinois, Louisiana, Mississippi, Missou↗
USGS2025

A review of abrupt permafrost thaw: Definitions, usage, and a proposed conceptual framework

Purpose of Review We review how ‘abrupt thaw’ has been used in published studies, compare these definitions to abrupt processes in other Earth science disciplines, and provide a definitive framework for how abrupt thaw should be used in the context of permafrost science. Recent Findings We address several aspects of permafrost systems necessary for abrupt thaw to occur and propose a framework for classifying permafrost processes as abrupt thaw in the future. Based on a literature review and our collective expertise, we propose that abrupt thaw refers to thaw processes that lead to a substantial persistent environmental change within a few decades. Abrupt thaw typically occurs in ice-rich permafrost but may be initiated in ice-poor permafrost by external factors such as hydrologic change (i.e., increased streamflow, soil moisture fluctuations, altered groundwater recharge) or wildfire. Summary Permafrost thaw alters greenhouse gas emissions, soil and vegetation properties, and hydrologic flow, threatening infrastructure and the cultures and livelihoods of northern communities. The term ‘abrupt thaw’ has emerged in scientific discourse over the past two decades to differentiate processes that rapidly impact large depths of permafrost, such as thermokarst, from more gradual, top-down thaw processes that impact centimeters of near-surface permafrost over years to decades. However, there has been no formal definition for abrupt thaw and its use in the scientific literature has varied considerably. Our standardized definition of abrupt thaw offers a path forward to better understand drivers and patterns of abrupt thaw and its consequences for global greenhouse gas budgets, impacts to infrastructure and land-use, and Arctic policy- and decision-making.

Current Climate Change Reports↗
USGS2025

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗
USGS2025

Submarine groundwater discharge creates cold‐water refugia that can mitigate exposure of heat stress in nearshore corals

Coral reef mortality around the world is accelerating due to human activities and rising sea temperatures that cause bleaching, which is expected to become more frequent. Our ability to predict which corals will be most resilient, however, remains limited due to insufficient information characterizing nearshore temperature and habitat conditions. In this study, we examine how submarine groundwater discharge (SGD) reduces nearshore water temperatures and exposure of corals to heat stress, complementing the understanding that SGD can adversely affect coral when it contains elevated nutrient concentrations. Data from fixed nearshore sensors and vertical depth profiles along ~100 km of the western shoreline of the Island of Hawai’i from 2003 to 2014 demonstrate that submarine groundwater discharge (SGD) can reduce nearshore water temperatures by 1 °C–5°C and create estuarine-like conditions with salinities as low as 20 PSU, where the prevalent coral species, Pocillopora meandrina, Porites lobata , and Montipora capitata , thrive. Time-series temperature records reveal that exposure to high ambient ocean temperatures, which are known to initiate bleaching events, are reduced up to 5%–46% of the time. Coral health surveys indicated coral bleaching in response to moderately high annual temperatures in 2010 and 2011, with more colonies affected farther from cold, SGD-fed waters. Synthesis of these results, along with coral response data following the more extreme marine heat wave of 2014–2015, demonstrates lower coral loss and greater coral recovery near groundwater seeps, particularly those with higher flux and influence on reducing nearshore water temperatures. Our results demonstrate that SGD may therefore provide a beneficial ecosystem service and enhance coral reef resilience, particularly where human-related nutrient additions to groundwater can be mitigated. The implications of our findings are relevant across tropical coasts where groundwater inputs can be substantial, such as the Caribbean and Indo-Pacific, and contribute to improving our understanding of coral sensitivity to gradients in temperature and nutrient stress. Improved management of groundwater resources could thus be vital to local–regional strategies for mitigating future heat stress.

Hawaii↗
USGS2026

Estimating paleotemperature using stable isotopes of soil-formed phyllosilicates from paleosols: A review

Fossilized soils, or paleosols, contain soil-formed phyllosilicates whose stable isotopic compositions may be used to calculate paleotemperature and thus reconstruct ancient terrestrial environments. Though paleosols are common in the geologic record, the use of phyllosilicates as paleotemperature proxies is limited in the literature owing to difficulties with selecting optimal paleosols, isolation from non-clay minerals and organic materials, mixtures of phyllosilicates in natural samples, wide variations of chemical compositions for phyllosilicates, and limited to undefined equilibrium fractionation factors between phyllosilicates-water. Here, we address these challenges by examining and comparing methods used for sample selection, mineral isolation, pretreatments, mineral identification, conventional and developing methods for oxygen and hydrogen isotopic analyses, and determination of phyllosilicate-water equilibrium fractionation factors, concluding with recommendations for best approaches for paleotemperature estimation. Additionally, we discuss how to identify and avoid detrital phyllosilicates, the impacts of diagenesis, comparison of stable isotope and non-isotope paleosol paleotemperature proxies, and challenges and opportunities for broadly using paleosols as paleoclimate archives. With ongoing efforts to refine this multi-faceted paleotemperature approach, the stable isotope geochemistry of soil-formed phyllosilicates continues to be an invaluable proxy system, enhancing our understanding of terrestrial paleoenvironments and paleoclimate.

Kansas↗
USGS2025

A geospatial analysis of water-quality threats from orphan wells in principal and secondary aquifers of the United States

Throughout the history of oil and gas production in the United States, millions of wells have been drilled for exploration and energy production. Hundreds of thousands of unplugged wells are no longer actively producing and are currently under orphan status, with no responsible party obligated for plugging. Orphan wells can pose threats to water resources by providing pathways for contaminants such as hydrocarbons and brines to migrate into water-supply aquifers. In this study, we investigate the potential threats to groundwater resources posed by orphan wells at the national scale. Water-quality data is extremely sparse in relation to orphan wells nationally and may not be suitable for identifying contamination from oil and gas development. We used geospatial and statistical methods to evaluate which principal and secondary aquifer systems may be most susceptible to contamination from orphan wells. Analysis involved three sets of susceptibility factors including: 1) factors related to the number and density of orphan wells; 2) factors that can threaten well integrity and contribute to transport of contaminants; and 3) factors related to groundwater withdrawal rates and the affected populations/communities in the event of water quality disturbances. From a dataset of 117,672 documented orphan wells, 64,203 fall within a principal aquifer system, while the remainder fall within a secondary aquifer system. By assessing the combination of well integrity and hydrogeologic factors within these aquifer systems, five groupings of principal aquifers were identified, where groups ranged from aquifer systems with high numbers of orphan wells, multiple well integrity threats and high withdrawals, to aquifers with a relatively low number of orphan wells, limited well integrity threats and minimal water use. Three regions of the country emerge containing aquifers with higher susceptibility to contamination from orphan oil and gas wells. These regions include 1) The Appalachian Basin (including the Pennsylvanian Aquifer System), 2) The Gulf Coast Aquifers (including the Coastal Lowlands Aquifer system) and 3) The California Aquifers (including the California Coastal Basin Aquifer system). This work is the first multivariate geospatial investigation of orphan wells and groundwater resources on a national scale, and sheds light on which aquifers are most susceptible to groundwater contamination from orphan wells.

Science of the Total Environment↗
USGS2025

Evaluation of daily stream temperature predictions (1979-2021) across the contiguous United States using a spatiotemporal aware machine learning algorithm

Stream temperature controls a variety of physical and biological processes that affect ecosystems, human health, and economic activities. We used 42 years (1979–2021) of data to predict daily summary statistics of stream temperature across >50,000 stream reaches in the contiguous United States using a recurrent graph convolution network. We comprehensively documented the performance – both across all reaches and by stream type (e.g., reservoir or groundwater influence) – as a baseline for future improvement. The model showed reach-level RMSE of <2 °C with 90 % prediction intervals that contain 90.7 % of observations. We also assessed how the model captured variability in ecologically relevant metrics (e.g., R 2 for annual 7-day maximum = 0.76; R 2 for days exceeding 25 °C = 0.75). This model does not outperform state-of-the-art machine learning efforts (e.g., RMSE ≤1.5 °C) due to a limited input set but does provide the most spatially complete modeling to date to support water availability assessments.

contiguous United States↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗