Geology ReportsSearch

SEARCH · Geology Reports

Results for “Journal of Hydrology”

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,659 records · Page 13Linked to original sources

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States

Wet meadow regeneration through restoration of biophysical feedbacks

Wet meadows are globally significant ecosystems that provide critical hydrological, ecological, and biogeochemical functions, yet their extent has declined dramatically due to land use changes and hydrologic alteration. These sedge-dominated wetlands exist at the drier end of the wetland gradient, maintained by shallow groundwater and periodic inundation. This paper is a global synthesis of the ecological, geomorphic, and hydrological dynamics of wet meadows, with an emphasis on alluvial systems, to inform effective restoration strategies. We compare wet meadows to other wetlands, classify them into palustrine, lacustrine, and alluvial types, then focus on alluvial wet meadows and discuss how their formation and persistence depend on ground and surface water interactions, sediment deposition and flow obstructions, all mediated by biological processes. In particular, we highlight the role of hydric graminoids in resisting erosion and maintaining soil cohesion, how beaver promote meadow persistence, and the significance of wet meadows as carbon sinks. We also present stratigraphic evidence demonstrating that incision, often triggered by anthropogenic activity or changing climate, is the primary mechanism of alluvial wet meadow degradation, resulting in water table decline and shifts in vegetation composition. Restoration requires reversing these incisional processes through techniques that elevate water tables, disperse flow and retain sediment—methods traditionally associated with either soil conservation or stream restoration. These include nature-based solutions that create obstructions such as beaver dams and their analogues, rock and wood-based obstructions and incision trench or gully filling and grading. Given their multifunctional value—including but not limited to flood attenuation, biodiversity support, and carbon sequestration—wet meadows warrant a focused restoration framework. This review advocates for a valley-floor scale restoration paradigm that integrates hydrological reconnection, sediment retention, and biological reinforcement to ensure long-term resilience of these systems in the face of changing climate and land use pressures.

Frontiers in Environmental Science

Dynamic baseflow storage estimates and the role of topography, geology and evapotranspiration on streamflow recession characteristics in the Neversink Reservoir Watershed, New York

Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.

New York

Peak-discharge frequency and potential extreme peak discharge for natural streams in the Brazos River basin, Texas

The 2-, 5-, 10-, 25-, 50-, and 100-year peak discharges were estimated for 186 streamflow-gaging stations with at least 8 years of data for natural streams in and near the Brazos River Basin, Texas. Multiple regression equations were developed to estimate peak-discharge frequency for the 2-, 5-, 10-, 25-, 50-, and 100-year recurrence intervals for each of three hydrologic regions that compose the Brazos River Basin. The equations for each region are a function of significant basin characteristics (explanatory variables). The significant explanatory variables among six that were tested are the contributing drainage area and stream slope for regions 1 and 2 and the contributing drainage area for region 3. For the three sets of equations, the coefficient of determination ranges from 0.59 to 0.93, and the standard error ranges from 0.184 to 0.391 log units. A larger coefficient of determination and a lower standard error generally are associated with the equations for hydrologic regions 2 and 3. Statistics from the regression analysis allow computation of the prediction interval associated with a given significance level for a peak-discharge frequency estimate. The regression equations can be used to estimate peak discharges for sites at, near, or away from sites with streamflow-gaging stations. The potential extreme peak-discharge curves as related to contributing drainage area were estimated for each of the three hydrologic regions from measured extreme peaks of record at 186 sites with streamflow-gaging stations and from measured extreme peaks at 37 sites without streamflow-gaging stations in and near the Brazos River Basin. The potential extreme peak-discharge curves generally are similar for hydrologic regions 1 and 2, and the curve for region 3 consistently is below the curves for regions 1 and 2, which indicates smaller peak discharges.

Texas

Hydrologic investigations and a preliminary conceptual model of the groundwater system at North Penn Area 1 Superfund Site, Souderton, Montgomery County, Pennsylvania

The U.S. Geological Survey (USGS) conducted hydrogeologic investigations, reviewed existing data, and developed a preliminary conceptual model of the groundwater system as part of technical support of the U.S. Environmental Protection Agency (EPA) at the North Penn Area 1 Superfund Site (hereafter, the NP1 Site) located within the Borough of Souderton in Montgomery County, Pennsylvania. Field work and monitoring took place during 2012–18. The area is underlain by sedimentary formations that form a fractured-rock aquifer used for drinking water and industrial supply. The EPA placed the Site on the National Priorities List in 1989, identifying tetrachloroethylene (PCE) and trichloroethylene (TCE) as contaminants of concern. During 2012–18, the USGS conducted field activities that included drilling an 82-foot (ft)-deep monitoring well (MG 2220) in 2016, reconstructing a 208-ft-deep former industrial production well (MG 668 [Granite Knitting Mill]), and collecting borehole geophysical and video logs and water levels from those and five additional wells, which ranged in depth from about 50 to 200 ft below land surface. Continuous water levels were collected during 2014–17, and a synoptic set of water levels were measured in April 2018 in the seven wells. The borehole geophysical logs (caliper, acoustic televiewer, natural gamma, single-point resistance, vertical flow, and fluid temperature and resistivity) and borehole video logs in the seven wells were evaluated to assess potential for lithologic correlation and to identify and describe water-bearing features, which included both low- and high-angle fractures and other openings oriented along dipping bedding planes, joints, or possible faults. Borehole geophysical logs collected by USGS in 1992 in a 300-ft-deep former production well near the Site were also evaluated. Few to no distinctive features were identified on geophysical logs (natural gamma and single-point resistance) that could be used for correlation, thus limiting this approach to determining local geologic structure. Extensive fracturing in the upper 62 ft of monitoring well MG 2220 indicates that the well was likely drilled through a zone of faulting, and other evidence of faulting is present in the area near the Site. Assessment of continuous water levels showed hydraulic connections among some wells as indicated by rising or falling water levels in response to changes in pumping rates at nearby wells. A map of water levels measured in April 2018 indicates potential for groundwater flow generally toward the stream to the south and southwest of the Site, but the limited water-level data are insufficient to describe vertical groundwater gradients or lateral gradients in any detail. Review of 1999–2022 volatile organic compound (VOC) monitoring data collected by the Pennsylvania Department of Environmental Protection for five monitoring wells indicates that the highest groundwater concentrations of PCE and TCE were found in samples from extraction well MG 2201 (S-1) downgradient from, and nearest to, the previously identified Site contaminant source area, and these concentrations fluctuated through time. PCE concentrations were higher than TCE concentrations in samples from all five monitoring wells and were much higher than TCE concentrations in samples from extraction well MG 2201 (S-1). Temporally variable recharge is a possible factor affecting observed fluctuations in PCE concentrations in groundwater samples from well extraction MG 2201 (S-1), as indicated by a general inverse relation between PCE concentrations and water levels in a nearby long-term observation well. The PCE concentration of 1,830 micrograms per liter (μg/L) in a May 2018 water sample from monitoring well MG 2220 was more than four times the PCE concentration of 444 μg/L in a December 2017 sample from the nearby extraction well MG 2201 (S-1), which is open to fewer fractures. Low concentrations of VOCs were measured in surface water at two stream sites downgradient from wells with the highest groundwater VOC concentrations at the Site, indicating that discharge of contaminated groundwater to the stream is likely. Development of a conceptual model of the groundwater system was constrained by limited data. In areas with no pumping, groundwater-flow directions generally are thought to be controlled by topography and geologic structure (bedding orientation) and likely to the south and southwest of the Site, with local flow directions affected by orientations of fractures, joints, and local faults. Additional investigations that could help improve the conceptual model of the groundwater system and help delineate the extent of groundwater contamination and its transport are discussed.

Pennsylvania

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York

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

Exploring landscape and geologic controls on spatial patterning of streambank groundwater discharge in a mixed land use watershed

Preferential groundwater discharge features along stream corridors are ecologically important at local and stream network scales, yet we lack quantification of the multiscale controls on the spatial patterning of groundwater discharge. Here we identify physical attributes that best explain variation in the presence and lateral extent of preferential groundwater discharges along two 5th order streams, the Housatonic and Farmington Rivers, and 32 1st to 4th order reaches across the Farmington River network. We mapped locations of preferential groundwater discharge exposed along streambanks using handheld thermal infrared cameras paired with high-resolution topographic and land use land cover datasets, surficial soil characteristic maps, and depth-to-bedrock geophysical measurements. The unconfined Housatonic River, MA, USA (12 km) had fewer discharge locations and less lateral extent (41 discharge locations with 38 m of active discharge/km of river) compared to the partially confined Farmington River, CT, USA (26 km; 169 discharge locations with 129 m of active discharge/km of river). Using a moving window analysis, we found along both rivers that discharge was more likely to occur where bank slopes were steeper, floodplain extent was narrower, and degree of confinement was higher. Along the Farmington River, groundwater discharge was more likely to occur where saturated hydraulic conductivity was higher and depth-to-bedrock was shallower. Among the 32 stream reaches surveyed (33.2 km of total stream length) within the Farmington River watershed, preferential discharge was observed in all but two stream reaches, varied from 0 to 25% of lateral extent along stream banks (mean = 6%), and was more likely to occur where stream reach slopes were steep, saturated hydraulic conductivity was high, and watershed urbanization was low. Our results show that, though both surface (e.g., topographic, land use land cover) and subsurface (e.g., soil characteristics, bedrock depth) factors control the prevalence of streambank preferential groundwater discharge, the dominant controls vary across valley settings and stream sizes.

Connecticut, Massachusetts

A novel drive-point multilevel system to investigate PFAS and other contaminants of global concern in the hyporheic zone of a wastewater effluent dominated stream

Contaminants found in treated wastewater discharged to streams, including pharmaceuticals and per- and polyfluoroalkyl substances (PFAS), are of global concern due to their deleterious effects on aquatic ecosystems and potential impacts to human health. Hyporheic zones have strong potential for contaminant attenuation. Assessing this potential requires collection of physical and biogeochemical data within the hyporheic zone. This study tested the applicability of a novel drive-point multilevel system (DP-MLS) for quantifying head profiles and characterizing contaminant concentrations in the hyporheic zone of a temperate region effluent dominated stream (EDS). DP-MLS, each with 4 ports, were installed in the stream bed at two sites, DS-1 and DS-2, 0.2 and 4.7 km downstream of the effluent outfall, respectively. Head profiles were measured and groundwater collected for analysis of pharmaceuticals and PFAS temporally over two years. The DP-MLS withstood rapid changes in stage, ice formation, and floating debris. Vertical hydraulic gradients (VHG) were generally upward but varied in magnitude indicating heterogeneity in hydraulic conductivity and variability in flow conditions. Upward VHG were also about 2X larger at DS-1 than at DS-2. Contaminant concentration profiles consistently showed penetration of pharmaceuticals and PFAS to 1 m below the bed at DS-2 while there was less penetration, lower groundwater concentrations, and more temporal variability in concentrations at DS-1. Integration of the physical and chemical data suggests weaker upwelling conditions at DS-2 are more easily reversed during periods of high stream stage, which could facilitate migration of wastewater contaminants into the bed. However, further studies incorporating other transport processes and reach scale dynamics are required to fully characterize these exchanges. Overall, this study demonstrates the efficacy of these novel DP-MLSs for characterization of the hyporheic zone and provides new insights into the occurrence, composition, and persistence of wastewater derived contaminants in the hyporheic zone of a well-studied EDS.

Iowa

Streamflow depletion caused by groundwater pumping: Fundamental research priorities for management-relevant science

Reductions in streamflow caused by groundwater pumping, known as “streamflow depletion,” link the hydrologic process of stream-aquifer interactions to human modifications of the water cycle. Isolating the impacts of groundwater pumping on streamflow is challenging because other climate and human activities concurrently impact streamflow, making it difficult to separate individual drivers of hydrologic change. In addition, there can be lags between when pumping occurs and when streamflow is affected. However, accurate quantification of streamflow depletion is critical to integrated groundwater and surface water management decision making. Here, we highlight research priorities to help advance fundamental hydrologic science and better serve the decision-making process. Key priorities include (a) linking streamflow depletion to decision-relevant outcomes such as ecosystem function and water users to align with partner needs; (b) enhancing partner trust and applicability of streamflow depletion methods through benchmarking and coupled model development; and (c) improving links between streamflow depletion quantification and decision-making processes. Catalyzing research efforts around the common goal of enhancing our streamflow depletion decision-support capabilities will require disciplinary advances within the water science community and a commitment to transdisciplinary collaboration with diverse water-connected disciplines, professions, governments, organizations, and communities.

Water Resource Research

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee

Transit time modeling framework for predicting freshwater salinization in urban catchments

The salinity of inland freshwaters is rising globally, particularly in urban watersheds where winter road deicers are widely applied. Attributing stream salinity dynamics to specific sources and transport pathways remains challenging due to episodic salt inputs, engineered drainage, and strong coupling between hydrology and subsurface storage. We present a modeling framework that couples climate-driven deicer build-up and wash-off with transient transit time distribution theory to simulate salt transport through drainage, interflow, and groundwater pathways. Applied to an urban watershed in Northern Virginia (USA), the model reproduces ten years of high-frequency stream salinity measurements across daily-to-decadal timescales. The calibrated model implies an average deicer application of 206 tonnes Cl −1 yr −1 , or roughly one 20 kg bag of rock salt person −1 yr −1 when normalized by the 20,000 people living in the watershed. In winter months, higher infiltration routes a large fraction of snowmelt and deicers into shallow subsurface pathways, enhancing vadose-zone and interflow contributions to stream salinity. Limited subsurface storage capacity and seasonal hydrologic turnover flush excess chloride from the vadose zone and groundwater during subsequent summer storms. By linking climate-driven deicer inputs, hydrologic connectivity, and stream water age, the framework provides a transferable basis for diagnosing and managing freshwater salinization in urban watersheds.

Virginia

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences

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

Hydrology and water and sediment quality at James Campbell National Wildlife Refuge near Kahuku, Island of Oahu, Hawaii

The James Campbell National Wildlife Refuge occupies two lowland marsh and pond complexes on the northern coastal plain of Oahu: the mostly natural ponds and wetlands of the Punamano Unit and the constructed ponds of the Kii Unit. The U.S. Fish and Wildlife Service manages the Refuge primarily to protect and enhance habitat for four endangered species of Hawaiian waterbirds. Kii Unit is fed by artesian wells and rainfall, whereas Punamano Unit is fed naturally by rainfall, runoff, and ground-water seepage. Streams drain from the uplands into lowland ditches that pass through Kii Unit on their way to the ocean. A high-capacity pump transfers water from the inner ditch terminus at Kii to the ocean outlet channel. Stormwaters also exit the inner ditch system over flood-relief swales near the outlet pump and through a culvert with a one-way valve. A hydrologic investigation was done from November 1996 through February 1998 to identify and quantify principal inflows and outflows of water to and from the Refuge, identify hydraulic factors affecting flooding, document ground-water/surface-water interactions, determine the adequacy of the current freshwater supply, and determine water and sediment quality. These goals were accomplished by installing and operating a network of stream-gaging stations, meteorology stations, and shallow ground-water piezometers, by computing water budgets for the two Refuge units, and by sampling and analyzing water and pond-bottom sediments for major ions, trace metals, and organic compounds. Streamflow during the study was dominated by winter stormflows, followed by a gradual recession of flow into summer 1997, as water that had been stored in alluvial fans drained to lowland ditches. Outflow at the ditch terminus in 1997 was 125 million gallons greater than measured inflow to the coastal plain, mainly reflecting gains from ground water along the ditches between outlying gages and the ditch terminus. Of the measured 1997 outflow, 98 percent was through the Kii outlet pump, with the outlet culvert valve only opening for brief periods during storms. Large volumes of stormflow overflowed the flood-relief swales unmeasured. The largest storm of the study, in November 1996, was estimated to have a flood frequency of about 3 to 4 years. Streamflow exceeded culvert capacity and overtopped Kamehameha Highway at Kalaeokahipa Stream and Hospital ditch. Slight overbank flooding in Kii ditch resulted strictly from high discharge. Minor overbank flooding farther out on the coastal plain probably was caused mainly by the small hydraulic gradients available to convey stormflows along the lowland ditches. Stormwaters flooded Kii ponds and flowed back upstream along Punamano ditch into Punamano marsh, introducing suspended sediment and possibly other contaminants to the Refuge. Two smaller storms in January 1997 resulted in smaller flows and no overbank flooding. The Kii outlet pump ran continuously for 7 days during the November 1996 storm and for 1 to 2 days during the January 1997 storms. During all three storms, the outlet culvert valve opened and the inner ditches overtopped the flood-relief swales, allowing free outflow of water from the inner ditch. Backwater effects hindered drainage during the January 1997 storms at Hospital ditch at Kamehameha Highway, and at Punamano ditch at Nudist Camp Road (where the backflow into Punamano marsh in November 1996 constituted an extreme backwater effect). A probable marine backwater effect was imposed at the ocean outlet ditch during the November 1996 storm through a combination of high spring tides and wave setup from large surf. Whether this backwater effect propagated upstream in the ditches to affect inland sites could not be determined conclusively. A sand plug may have built up in the ocean outlet channel before the November 1996 storm, but if so, it probably washed out prior to, or early in the storm, and was not present at the time of peak stage at inland sites. A season-long buildup of the sand plug in late 1997 was inferred from rising water levels in the outlet ditch. Seawater flows up the outlet channel or over the sand berm and into the outer ditch system on most high tides, and particularly during spring high tides. Ponds and ditches of the Refuge and surrounding lowlands have mud- and clay-lined bottoms that form an effective confining unit and inhibit interaction with an underlying shallow limestone aquifer. At Kii Unit, pond levels are higher than adjacent ditch levels and underlying ground-water levels, establishing lateral and downward head gradients that could foster seepage losses from the ponds. Regional ground-water discharge from the Koolau aquifer to the coastal-plain sediments is mostly diffuse, but is concentrated where ridges of Koolau Basalt plunge beneath coastal-plain sediments near Punamano Unit and at the head of Hospital ditch. Kii ditch gains brackish ground water downstream of Kamehameha Highway. Wastewater disposal from the sewage treatment plant adjoining Kii Unit poses little or no threat to Refuge habitat. Disposal is at six injection wells located 0.45 mi away at Kahuku, and the wells inject into confined limestone aquifers that do not extend to Kii Unit. The natural freshwater supply to Punamano Unit is adequate for maintaining the wildlife habitat, judging from stable pond levels and low salinities there. A monthly water budget for Punamano showed an apparent annual deficit in measurable flows in 1997, requiring unmeasured ground-water gains equalling 51 inches of water. The freshwater supply to Kii Unit is inadequate according to Refuge managers, because there is not enough water to manipulate levels adequately in the ponds during most of the year, and particularly during the driest months. This is confirmed by monthly deficits in the water budget for the Kii ponds during summer months. However, the Kii budget showed an annual surplus in measurable flows for 1997 equalling 24 inches of water. Unmeasured losses are required to explain the apparent annual surplus, such as discharge to the ditches through pond water-control structures and downward and/or lateral ground-water seepage. The apparent surplus at Kii is strictly hydrologic and is not a surplus in a management sense; it cannot be stored or used to supply the Refuge, but instead reflects losses from the system that render this amount of water unavailable for use. The budget year, 1997, was drier than normal (24 percent below long-term mean rainfall) and so the measured potential evaporation for 1997 was probably higher than the long-term mean. Few metals or organic compounds of potential concern were detected in pond and ditch waters and in pond-bottom sediments. Detected pesticides were at trace levels or just above minimum reporting limits. Exceptions that exceeded quality guidelines for freshwater sediment were copper and zinc in sediment from Kii ponds C and D, and copper in sediment from Punamano north pond. Therefore, urban and agricultural runoff probably have contributed little in the way of harmful metals or organic compounds to the Refuge, although the potential for such contribution remains from periodic flooding of the ponds by ditch stormflows. Salinity was low throughout most Refuge waters, qualifying as fresh to slightly brackish and suitable for the environmental needs of Refuge fauna. Higher salinities have been observed in ditches during past periods of sugarcane cultivation and saltwater aquaculture, however. Resumption of saltwater aquaculture could raise ditch salinities if saltwater effluents are disposed directly into the ditches, as they were in the past.

Hawaii