USGS Science⌕ Search

SEARCH · USGS Science

Results for “Current Climate Change Reports”

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

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2024

Participatory engagement to reduce communication gaps

Underserved communities, especially those in coastal areas in Puerto Rico, face significant threats from natural hazards such as hurricanes and rising sea levels. Limited funding hinders the investment in costly mitigation measures, increasing exposure to natural disasters. Providing coastal resources and data products through effective communication mechanisms is fundamental to improving the well-being of these underserved coastal communities. The overall objectives of the pilot effort to engage and connect with underserved coastal communities in Puerto Rico were the following: (1) compile a comprehensive database of the projects and resources relevant to natural hazards in Puerto Rico; (2) foster connections with Puerto Rican interested parties to better understand their priorities regarding coastal hazards and provide them with pertinent U.S. Geological Survey (USGS) resources; and (3) identify knowledge gaps to guide future USGS projects in Puerto Rico. Here we outline our participatory engagement framework and process, along with two specific resources developed with the information collected from this effort. These resources are available in English and Spanish and consist of user-friendly, non-technical information products. Among them are: (1) a website where users can learn about USGS research on landslides, hurricanes, earthquakes, water resources, coastal hazards, tsunamis, and ecosystem hazards and environmental contaminants, and (2) a geonarrative highlighting shoreline changes in Puerto Rico with sections on historical shoreline trends, hurricane impacts, and potential solutions that could help protect communities and mitigate coastal hazards. Continuing participatory engagement in future projects could enhance the accessibility and usability of natural hazards resources within the community.

San Juan↗
USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2025

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗
USGS2025

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗
USGS2024

Anaerobic biodegradation of perfluorooctane sulfonate (PFOS) and microbial community composition in soil amended with a dechlorinating culture and chlorinated solvents

Perfluorooctane sulfonate (PFOS), one of the most frequently detected per - and polyfluoroalkyl substances (PFAS) occurring in soil, surface water, and groundwater near sites contaminated with aqueous film-forming foam (AFFF), has proven to be recalcitrant to many destructive remedies, including chemical oxidation. We investigated the potential to utilize microbially mediated reduction (bioreduction) to degrade PFOS and other PFAS through addition of a known dehalogenating culture, WBC-2, to soil obtained from an AFFF-contaminated site. A substantial decrease in total mass of PFOS (soil and water) was observed in microcosms amended with WBC-2 and chlorinated volatile organic compound (cVOC) co-contaminants — 46.4 ± 11.0 % removal of PFOS over the 45-day experiment. In contrast, perfluorooctanoate (PFOA) and 6:2 fluorotelomer sulfonate (6:2 FTS) concentrations did not decrease in the same microcosms. The low or non-detectable concentrations of potential metabolites in full PFAS analyses, including after application of the total oxidizable precursor assay, indicated that defluorination occurred to non-fluorinated compounds or ultrashort-chain PFAS. Nevertheless, additional research on the metabolites and degradation pathways is needed. Population abundances of known dehalorespirers did not change with PFOS removal during the experiment, making their association with PFOS removal unclear. An increased abundance of sulfate reducers in the genus Desulfosporosinus (Firmicutes) and Sulfurospirillum (Campilobacterota) was observed with PFOS removal, most likely linked to initiation of biodegradation by desulfonation. These results have important implications for development of in situ bioremediation methods for PFAS and advancing knowledge of natural attenuation processes.

Science of the Total Environment↗
USGS2025

Dendroseismological investigation of redwood trees along the North Coast section of the San Andreas Fault

Sequoia sempervirens (coast redwood) tree rings have the potential to annually resolve late-Holocene earthquakes on the northern San Andreas Fault based on direct (e.g., physical damage) and indirect (e.g., co-seismic environmental change) impacts, but scarcity of suitable samples and challenges crossdating this long-lived species have limited progress. More precise dating of the pre-1906 (penultimate) earthquake can improve hazard assessment and understanding of rupture segmentation. We target old trees (maximum >815 yr) along the North Coast section of the fault (increment cores via rope-climbing, 11 living trees; plunge cuts, 23 stumps) and employ complementary disturbance detection methods including radial-growth averaging (tree- and series-level), cataloging anatomical indicators (e.g., traumatic resin ducts, TRD), and dating structural components (e.g., reiterated trunks, leans). Multi-centennial ring-width chronologies at Fort Ross (1569−2023) and Gualala (1397−2023) promote continued study with incomplete crossdating limiting utilization of some series. Growth pulses (reductions, releases) and TRD dispersed across the record reflect dynamic environments that obfuscate detection of earthquake signals. The 1906 earthquake did not leave strong signatures on most trees, and when it did, within-tree response varied from normal presentation to discoloration, TRD, and missing rings. Synchrony of indicators at both locations identified 1678−1680 (6 of 15 trees) and 1698−1700 (8 of 16 trees) as the strongest disturbances among dated rings in the time range of the penultimate earthquake, peaking at 1698 (15.7 % of possible growth and anatomical indicators), but the triggering mechanisms for these events are unknown.

California↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS2025

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗
USGS2025

Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration

Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model-data integration workflow to address uncertainty in groundwater response to beaver-induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning-based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down-valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down-valley underflow, from 2% during dry pond periods to 20% during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down-valley underflow, minimizes the effect of beaver-induced inundation on water quality. We emphasize the need for field-scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response.

Colorado↗
USGS2026

Surface variable‐based machine learning for scalable arsenic prediction in undersampled areas

In the United States, private wells are not federally regulated, and many households do not test for Arsenic (As). Chronic exposure is linked with multiple health outcomes, and risk can change sharply over short distances and with well depth. Coarse maps or sparse sampling often miss exceedances. Most existing models operate at ∼1 km resolution and use groundwater chemistry or detailed geologic logs, which limits their use in undersampled areas where improved guidance is most needed. We overcome these limitations by developing a machine learning model for Minnesota, USA, that predicts As exposure risk using only surficial variables from remote sensing and global data sets. Variables related to surface water hydrology and geomorphology are selected based on mechanistic links that control redox conditions and As mobilization. Local training was essential, and surficial geology variables that are more sensitive to local conditions were needed to maximize model accuracy. The resulting complete model was sufficiently sensitive to generate accurate and detailed risk maps and depth profiles of As concentrations above the 10 μg/L maximum contaminant level. Accuracy depended on local training data density. We identified a training data density of 0.07 wells/km 2 as a practical target for stable county-level performance. Maps of exceedance probabilities highlight priority areas for testing that are particularly important in rural communities that have received less sampling. These results support public health action by guiding where to install wells and where to test them, how much new sampling is needed, and where treatment outreach is most urgent.

Minnesota↗
USGS2024

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware↗
USGS2024

Spatial differences in soil nutrients along a hydrographic gradient on floodplains in Dongting Lake

The spatial heterogeneity of soil nutrients is crucial for the water bird and whole floodplain wetland ecosystem in large lakes, and it is influenced by the dramatic water level changes and sedimentation progress in West Dongting Lake (WDL). Soil samples were collected at various soil depths along the Yuan River and Li River that feed into WDL. The concentrations of soil total organic carbon (TOC), total nitrogen (TN), total phosphorus (TP), and soil grain size were tested. The stoichiometric ratios of C, N, P, and the mean value of soil grain size ( Mz ) were calculated. The differences of soil TOC, TN, TP and the stoichiometric ratio at different sites and soil depths were compared. Linear regression was used to explore the relationships of Mz and nutrient concentrations, and relationships between TOC, TN, and TP. Redundancy analysis was used to explore the relationship between soil nutrients, heavy metal concentrations, and plant community diversity. The results showed that the distributions of soil TOC, TN, and TP concentrations differed across regions in west Dongting Lake along the Yuan and Li Rivers. Total organic carbon concentration differed at different sedimentation depths. Soil grain size showed negative effect with soil TOC, TN, and TP concentrations in this region. Plant community diversity correlated positively with soil TOC and negatively with Hg. West Dongting Lake was N limited despite the high wet deposition of N. It could potentially be attributed to the insufficient presence of aerobic environments for microbes during intermittent flooding of the floodplain, coupled with feeble mineralization. This study can provide valuable insights for the conservation of water bird habitats and wetland ecosystems.

West Dongting Lake↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗