skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Methods for Quantifying Interactions Between Groundwater and Surface Water
Driven by the need for integrated management of groundwater (GW) and surface water (SW), quantification of GW–SW interactions and associated contaminant transport has become increasingly important. This is due to their substantial impact on water quantity and quality. In this review, we provide an overview of the methods developed over the past several decades to investigate GW–SW interactions. These methods include geophysical, hydrometric, and tracer techniques, as well as various modeling approaches. Different methods reveal valuable information on GW–SW interactions at different scales with their respective advantages and limitations. Interpreting data from these techniques can be challenging due to factors like scale effects, heterogeneous hydrogeological conditions, sediment variability, and complex spatiotemporal connections between GW and SW. To facilitate the selection of appropriate methods for specific sites, we discuss the strengths, weaknesses, and challenges of each technique, and we offer perspectives on knowledge gaps in the current science.  more » « less
Award ID(s):
1752995
PAR ID:
10612080
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Annual Reviews Environmental Resources
Date Published:
Journal Name:
Annual Review of Environment and Resources
Volume:
49
Issue:
1
ISSN:
1543-5938
Page Range / eLocation ID:
623 to 653
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This study attempts to integrate a Surface Water (SW) model Soil and Water Assessment Tool (SWAT) with an existing steady-state, single layer, unconfined heterogeneous aquifer Analytic Element Method (AEM) based Ground Water (GW) model, named Bluebird AEM engine, for a comprehensive assessment of SW and GW resources and its management. The main reason for integrating SWAT with the GW model is that the SWAT model does not simulate the distribution and dynamics of GW levels and recharge rates. To overcome this issue, often the SWAT model is coupled with the numerical GW model (either using MODFLOW or FEFLOW), wherein the spatial and temporal patterns of the interactions are better captured and assessed. However, the major drawback in integrating the two models (SWAT with—MODFLOW/FEM) is its conversion from Hydrological Response Unit’s (HRU)/sub-basins to grid/elements. To couple them, a spatial translation system is necessary to move the inputs and outputs back and forth between the two models due to the difference in discretization. Hence, for effective coupling of SW and GW models, it may be desirable to have both models with a similar spatial discretization and reduce the need for rigorous numerical techniques for solving the PDEs. The objective of this paper is to test the proof of concept of integrating a distributed hydrologic model with an AEM model at the same spatial units, primarily focused on surface water and groundwater interaction with a shallow unconfined aquifer. Analytic Element Method (AEM) based GW models seem to be ideal for coupling with SWAT due to their innate character to consider the HRU, sub-basin, River, and lake boundaries as individual analytic elements directly without the need for any further discretization or modeling units. This study explores the spatio-temporal patterns of groundwater (GW) discharge rates to a river system in a moist-sub humid region with SWAT-AEM applied to the San Jacinto River basin (SJRB) in Texas. The SW-GW interactions are explored throughout the watershed from 2000–2017 using the integrated SWAT-AEM model, which is tested against stream flow and GW levels. The integrated SWAT-AEM model results show good improvement in predicting the stream flow (R2 = 0.65–0.80) and GW levels as compared to the standalone SWAT model. Further, the integrated model predicted the low flows better compared to the standalone SWAT model, thus accounting for the SW-GW interactions. Almost 80% of the stream network experiences an increase in groundwater discharge rate between 2000 and 2017 with an annual average GW discharge rate of 1853 Mm3/year. The result from the study seems promising for potential applications of SWAT-AEM coupling in regions with considerable SW-GW interactions. 
    more » « less
  2. Abstract Groundwater/surface‐water (GW/SW) exchange and hyporheic processes are topics receiving increasing attention from the hydrologic community. Hydraulic, chemical, temperature, geophysical, and remote sensing methods are used to achieve various goals (e.g., inference of GW/SW exchange, mapping of bed materials, etc.), but the application of these methods is constrained by site conditions such as water depth, specific conductance, bed material, and other factors. Researchers and environmental professionals working on GW/SW problems come from diverse fields and rarely have expertise in all available field methods; hence there is a need for guidance to design field campaigns and select methods that both contribute to study goals and are likely to work under site‐specific conditions. Here, we present the spreadsheet‐based GW/SW‐Method Selection Tool (GW/SW‐MST) to help practitioners identify methods for use in GW/SW and hyporheic studies. The GW/SW‐MST is a Microsoft Excel‐based decision support tool in which the user selects answers to questions about GW/SW‐related study goals and site parameters and characteristics. Based on user input, the tool indicates which methods from a toolbox of 32 methods could potentially contribute to achieving the specified goals at the site described. 
    more » « less
  3. Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and downstream of the fractured hillslope (non- fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone were high, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. We interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow. 
    more » « less
  4. Abstract Several studies have focused on the importance of river bathymetry (channel geometry) in hydrodynamic routing along individual reaches. However, its effect on other watershed processes such as infiltration and surface water (SW)‐groundwater (GW) interactions has not been explored across large river networks. Surface and sbsurface processes are interdependent, therefore, errors due to inaccurate representation of one watershed process can cascade across other hydraulic or hydrologic processes. This study hypothesizes that accurate bathymetric representation is not only essential for simulating channel hydrodynamics but also affects subsurface processes by impacting SW‐GW interactions. Moreover, quantifying the effect of bathymetry on surface and subsurface hydrological processes across a river network can facilitate an improved understanding of how bathymetric characteristics affect these processes across large spatial domains. The study tests this hypothesis by developing physically based distributed models capable of bidirectional coupling (SW‐GW) with four configurations with progressively reduced levels of bathymetric representation. A comparison of hydrologic and hydrodynamic outputs shows that changes in channel geometry across the four configurations has a considerable effect on infiltration, lateral seepage, and location of water table across the entire river network. For example, when using bathymetry with inaccurate channel conveyance capacity but accurate channel depth, peak lateral seepage rate exhibited 58% error. The results from this study provide insights into the level of bathymetric detail required for accurately simulating flooding‐related physical processes while also highlighting potential issues with ignoring bathymetry across lower order streams such as spurious backwater flow, inaccurate water table elevations, and incorrect inundation extents. 
    more » « less
  5. Abstract It is important to understand how point measurements across spatially heterogeneous ecosystems are scaled to represent these systems. Stream biogeochemistry presents an illustrative example because water quality concerns within stream networks and recipient water bodies motivate heterogeneous watershed studies. Measurements of the stream water‐groundwater (SW‐GW) interface (i.e., the shallow stream subsurface) are well‐documented for point‐scale sampling density measurements (i.e., cm2–m2features), but poorly characterized for network‐scale sampling density measurements (i.e., km2; stream reaches and networks). Sampling the SW‐GW interface is more time and labor intensive than surface water sampling, meaning sample point selection must be made with care for network‐scale analyses. In this study, we endeavor to determine which of two common spatial sampling schemes is appropriate for characterizing SW‐GW interface biogeochemistry across a third‐order stream network, focusing on dissolved organic carbon. The first scheme, called Local Sampling, focuses on characterizing small‐scale (< 10 m2) variability produced by the local physical and biogeochemical heterogeneity, with fewer points across the stream network. The second scheme, called Longitudinal Sampling, has approximately the same number of measurements distributed over many more points across the stream network with less local variability characterization. This comparison reveals that selection of a Local Sampling versus a Longitudinal Sampling scheme influences the biogeochemical pattern interpretation at the stream network scale. Additionally, this study found that increasing observation efforts at the local scale added limited information for reach‐ to network‐scale biogeochemical patterns, suggesting that emphasis should be placed on characterizing variability across broader spatial scales with the Longitudinal Sampling approach. 
    more » « less