USGS Science⌕ Search
Home/Archive/USGS

USGS · 70022302

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

P. M. Barlow·L.A. DeSimone·A.F. Moench
View original source ↗Export metadata ↓

Abstract

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to streamstage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream-aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river-aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.Analytical step-response functions, developed for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream-aquifer systems. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank seepage rates and bank storage.

Explore related subjects

Journal of HydrologyArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

P. M. Barlow, L.A. DeSimone, A.F. Moench. 2000. Aquifer response to stream-stage and recharge variations. II. Convolution method and applications. https://doi.org/10.1016/s0022-1694(00)00176-1

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Reservoir releases and land cover interact to drive event-scale nitrate export in a large agricultural basin

Understanding the drivers of nitrate export in rivers is critical for developing effective nutrient management strategies. However, few studies have explored event-scale drivers of export in large river basins with human modifications like reservoirs. Here, we analyzed nitrate concentration-discharge (C-Q) relationships from 215 events at the outlet of the Kansas River Basin, USA (155,690 km 2 ) from 2014 to 2022 to (1) characterize event-scale nitrate export behaviors in a large agricultural basin, and (2) determine how different event characteristics are linked to these behaviors. We found that C-Q behaviors varied greatly, with 60% of events producing nitrate enrichment (n = 130) and 40% of events producing nitrate dilution (n = 85). These behaviors were correlated with complex spatial interlinkages between climate and land cover: across most of the basin, nitrate enrichment was correlated with drier antecedent conditions, but in wetter areas with higher proportions of urban and forested land cover, enrichment was more strongly correlated with precipitation magnitude/intensity. This difference in hydroclimatic controls on nitrate export might be related to differential distributions of nitrate sources within these land covers. Upstream of major reservoirs, however, neither variable was strongly correlated with C-Q behavior, suggesting that nitrate attenuation within reservoirs decouples event-scale concentration signals in upstream waters from those downstream. Reservoir outflows had variable impacts on C-Q behavior, reflecting reservoir-specific variations in nitrate attenuation efficiency. Together, these results identify specific complex interactions between hydroclimate, land cover, reservoir positioning, and individual reservoir properties that control event-scale nitrate export from large basins.

Colorado, Kansas, Nebraska↗
USGS2026

Deep learning error post-processing improves stochastic watershed modeling

Hydrologic extremes, including floods and droughts, pose substantial societal risks that are expected to intensify with climate change. Deterministic watershed models (DWMs) remain a mainstay for modeling these extremes, but lack explicit representation of uncertainty, limiting their utility for risk-informed planning. Stochastic watershed models (SWMs) address this limitation by generating ensembles of streamflow via models of observed DWM residuals. However, most SWMs struggle with the complex dependence between DWM residuals and the underlying hydrologic state, which can complicate stochastic simulations under nonstationary climates. Deep learning (DL) models, whether used as standalone models or post-processors for process-based DWMs, offer a pathway to address this challenge by reducing conditional dependence. In this study, we evaluate SWMs applied to seven models: three process-based models (PRMS, Hymod, and HBV), their hybrid process-DL counterparts, and a pure DL DWM, focusing on daily simulations and extremes under both historical conditions and synthetic climate change scenarios. Results for a case study watershed in Massachusetts show that SWMs applied to hybrid or pure DL DWMs consistently outperform those applied to process-based DWMs. However, an SWM applied to the pure DL model exhibits weaknesses at low flows for this study basin, underscoring the value of hybrid approaches. Extending the analysis across 73 additional basins demonstrates that these improvements are robust and generalizable statewide. This work highlights the potential of a DL-enhanced stochastic watershed modeling framework to advance hydrologic risk prediction under changing climate conditions, offering a scalable methodology for integrating uncertainty into watershed modeling for long-term planning.

Journal of Hydrology↗
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↗