Search USGS⌕ Search

SEARCH · Search USGS

Results for “Hydrologic Processes”

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.

536 records · Page 7Linked to original sources
USGS2026

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado↗
USGS2025

A spatiotemporal deep learning approach for predicting daily air-water temperature signal coupling and identification of key watershed physical parameters in a montane watershed

Seasonal shifts from runoff to groundwater dominance influence daily headwater stream temperatures, especially where local groundwater input is strong. This input buffers temperature during hot periods, supporting cold-water habitats. Recent studies use air–water temperature signal metrics to identify zones of strong stream–groundwater connectivity. While Previous studies used air–water signal ratios as proxies for groundwater influence but were limited to specific sites and periods, without dynamic forecasting. This study is the first to forecast daily A r as a spatiotemporal signal using a Graph Convolutional Network–Long Short-Term Memory (GCN-LSTM) model. The model was trained using hydroclimate data (air temperature, precipitation, shortwave radiation, streamflow) and watershed physical features (e.g., sand content, slope). Results showed high predictive skill, achieving R 2 (NSE, RMSE) of 0.86 (0.73, 0.0004) for one-day-ahead to 0.52 (0.50, 0.0009) for seven-days ahead forecasts. Prior studies often have not explicitly incorporated spatial hydrogeologic drivers, but this model explicitly incorporates them to assess their impact on A r forecasting and stream-groundwater connectivity. Feature analysis identified mean sand, elevation, slope, clay, and TWI as key predictors of A r . Stronger groundwater signals appeared in hillslopes, elevations, and tributaries, highlighting watershed influence on streamflow. However, limitations include reliance on historical air–water temperature patterns for training and limited representation of extreme climate conditions. Despite these limitations, unlike previous studies relying on measured in-situ stream and air temperature, this study forecasts A r directly from climate and physiographic features after training, avoiding in-situ data requirements. Findings aiding predictions of stream ecosystem resilience.

New York↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗
USGS1998

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

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2024

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

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2026

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

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report↗
USGS2026

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

Mapping karst groundwater flow paths and delineating recharge areas for springs in the Little Sequatchie and Pryor Cove watersheds, Tennessee

The Little Sequatchie River and Pryor Cove Branch, in southern Tennessee, drain the eastern escarpment of the Cumberland Plateau to the Sequatchie River near the southern end of the Sequatchie Valley. The Little Sequatchie River is the largest tributary to the Sequatchie River by drainage area, covering over 120 square miles. The hydrology of the two drainage areas has been largely altered by karst processes, which has caused the majority of the streams to sink at the contact between the Mississippian Pennington Formation and the underlying Mississippian Bangor Limestone. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service and Tennessee Department of Environment and Conservation, initiated a study in 2021 to map the karst groundwater pathways in both watersheds in order to delineate recharge areas for several springs. One of these springs, Sequatchie Cave, represents a significant habitat for two Species of Greatest Conservation Need, the Glyphopsyche sequatchie (Sequatchie caddisfly) and the federally endangered Marstonia ogmorhaphe (royal marstonia). Springs and springflow-dominated streams in the Little Sequatchie River valley and Pryor Cove also provide water for agricultural practices and serve as a drinking water source for nearby communities. During the study, a total of 25 dye injections were conducted over eight rounds from January 2022 through March 2023. Dye traces from these injections helped to delineate recharge areas for six major springs, ranging from 7.3 to 65.2 square miles in area. The majority of the dye traces remained subsurface (from sinkpoint to recovery site) for long distances, with karst groundwater travelling nearly 8 miles before resurfacing. The dye traces also had rapid traveltimes, often travelling hundreds to thousands of feet per hour. The goal of this project was to provide scientific data related to karst groundwater pathways and spring recharge areas to aid State and Federal agencies in making informed decisions to protect and preserve this unique and vulnerable karst system.

Tennessee↗
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↗
USGS2000

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

Hydrogeology of the Chickasaw National Recreation Area, Murray County, Oklahoma

The Travertine District (Park) of the Chickasaw National Recreation Area, operated and maintained by the National Park Service, is near the City of Sulphur in south-central Oklahoma. The Park was established in 1902 because of its unique hydrologic setting, which includes Rock Creek, Travertine Creek, numerous mineralized and freshwater springs, and a dense cover of riparian vegetation. Since the turn of the century several flowing artesian wells have been drilled within and adjacent to the Park. Discharge from many of these springs and the numbers of flowing wells have declined substantially during the past 86 years. To determine the cause of these declines, a better understanding of the hydrologic system must be obtained. The U.S. Geological Survey, in cooperation with the National Park Service, has appraised hydrologic information obtained for the Park from several studies conducted during 1902-87. The principal geologic units referred to in this report are the Arbuckle Group and the overlying Simpson Group. These rocks are of Upper Cambrian to Middle Ordovician age and are composed of dolomitic limestone, with some sandstones and shales in the Simpson Group. Surface geologic maps give a general understanding of the regional subsurface geology, but information about the subsurface geology within the Park is poor. The Simpson and Arbuckle aquifers are the principal aquifers in the study area. The two aquifers are not differentiated readily in some parts of the study area because of the similarity of the Simpson and Arbuckle rocks; thus, both water-bearing units are referred to frequently as the Arbuckle-Simpson aquifer. The aquifers are confined under the Park, but are unconfined east and south of the Park. Precipitation on the outcrop area of the Arbuckle aquifer northeast and east of the Park recharges the freshwater springs (Antelope and Buffalo Springs) near the east boundary of the Park. The source of water from mineralized springs located in the central part of the Park, and flowing wells within and north of the Park, is believed to be a mix of waters from rocks of the Arbuckle and Simpson Groups. The source of water from two highly mineralized springs, Bromide and Medicine, that ceased to flow in the early 1970’s is believed to be from the Simpson Group. Water-quality characteristics reflect the sources of ground water in the study area. The highly mineralized springs near the western end of the Park are a sodium chloride type with dissolved solids greater than 4,500 mg/L. The freshwater springs near the eastern end of the Park are a calcium bicarbonate type with total dissolved solids of less than 400 mg/L. Flow from the artesian wells has declined substantially during the past 86 years and the wells are estimated to currently discharge only about 10 percent of the total flow reported in 1939. The depletion is believed to be caused by a gradual lowering of the hydraulic head within the aquifer. The influence on the hydrologic system of local municipal and industrial pumping from the Arbuckle-Simpson aquifer is difficult to discern because the system is much more sensitive to precipitation than to pumpage. Ground-water levels and spring flows in this region respond rapidly to precipitation. The effects of withdrawals from the City of Sulphur and Oklahoma Gas and Electric Company power-plant water-well fields are not discernible at wells and springs. The hydrologic system may be influenced by pumping, particularly during extended dry periods of several years, but the impact of pumping on the system cannot be determined without further investigation.

Oklahoma↗
USGS2026

Effects of groundwater withdrawals for water bottling and municipal use, Wards Brook Valley, Maine and New Hampshire

Hydrologic models for the Wards Brook valley near Fryeburg, Maine were developed for historical (2016 – 2021) and hypothetical future conditions (2046 – 2065 and 2080 – 2099) to understand the effects of groundwater withdrawals for bottled water and municipal use on hydrologic conditions (stream base flows and groundwater levels). Analyses showed that the simulated base flows in Wards Brook were reduced because of pumping for both municipal water supplies and for water bottling, and about half of the total pumping impact on the base flows in Wards Brook was from the bottled water extraction. Simulated flows were greater than the minimum recommended streamflow of 2,180 cubic meters per day (400 gallons per minute) throughout the historical period. Simulated groundwater levels at two of three nearby ponds (Round Pond and Davis Pond) were minimally affected by pumping conditions, and effects were primarily from the municipal well closest to the ponds. Several estimates of future projected recharge were used to understand the potential effects of groundwater withdrawals on hydrologic conditions under multiple hypothetical climate conditions. Annual projected recharge rates in the mid- and late-21st century from two climate scenarios (stabilized greenhouse-gas emissions and high greenhouse-gas emissions) were similar to rates for 2016 – 2021. However, monthly recharge patterns for the future periods shifted toward more recharge in the winter months (December, January, and February) and less recharge in April, May, and October relative to 2016 – 2021. The lowest mean monthly base flows from the future emission scenarios all remain larger than the minimum recommended streamflow and indicate no long-term declines in flow relative to historical conditions. However, simulated base flows during hypothetical 3-year drought scenarios declined below minimum recommended streamflow during the summer months in the stabilized- and high-emission scenarios in the mid-21st century. Although water is generally plentiful in the Wards Brook valley, reduced pumping may be needed to maintain streamflows in Wards Brook under future climate conditions similar to modeled drought scenarios.

EarthArXiv↗
USGS2025

Quantitative subsurface characterization illuminates the origin of the Quaternary Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer (MRVA) is vital to U.S. food security and global agricultural supply. However, quantitative understanding of its Quaternary origin, architecture, and hydrologic function remains incomplete. Here we develop a three-dimensional hydrostratigraphic model to characterize the deposition of clay and silt, fine-medium sands, and graveliferous sands using lithologic data from 75,000 boreholes compiled across the Lower Mississippi Valley and a geostatistical method—interval kriging. We find that cyclic glacial entrenchments, evidenced by remnants of pre-Wisconsinan postglacial sediments, alongside geodynamic activities shaped the MRVA basal configuration. Stratal weakening from faulting and salt diapirism enhanced glacial incision and thereby produced abrupt aquifer thickening. We demarcate the top of graveliferous sands as the regional marker of the Pleistocene-Holocene transition. The MRVA hydrostratigraphy reveals hydrologic function and geologic controls on groundwater storage and quality, advancing the assessment of aquifer sustainability under a changing climate, with implications for alluvial aquifers globally.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗