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: Shallow and local or deep and regional? Inferring source groundwater characteristics across mainstem riverbank discharge faces
Abstract Riverbank groundwater discharge faces are spatially extensive areas of preferential seepage that are exposed to air at low river flow. Some conceptual hydrologic models indicate discharge faces represent the spatial convergence of highly variable age and length groundwater flowpaths, while others indicate greater consistency in source groundwater characteristics. Our detailed field investigation of preferential discharge points nested across mainstem riverbank discharge faces was accomplished by: (1) leveraging new temperature‐based recursive estimation (extended Kalman Filter) modelling methodology to evaluate seasonal, diurnal, and event‐driven groundwater flux patterns, (2) developing a multi‐parameter toolkit based on readily measured attributes to classify the general source groundwater flowpath depth and flowpath length scale, and, (3) assessing whether preferential flow points across discharge faces tend to represent common or convergent groundwater sources. Five major groundwater discharge faces were mapped along the Farmington River, CT, United States using thermal infrared imagery. We then installed vertical temperature profilers directly into 39 preferential discharge points for 4.5 months to track vertical discharge flux patterns. Monthly water chemistry was also collected at the discharge points along with one spatial synoptic of stable isotopes of water and dissolved radon gas. We found pervasive evidence of shallow groundwater sources at the upstream discharge faces along a wide valley section with deep bedrock, as primarily evidenced by pronounced diurnal discharge flux patterns. Discharge flux seasonal trends and bank storage transitions during large river flow events provided further indication of shallow, local sources. In contrast, downstream discharge faces associated with near surface cross cutting bedrock exhibited deep and regional source flowpath characteristics such as more stable discharge patterns and temperatures. However, many neighbouring points across discharge faces had similar discharge flux patterns that differed in chloride and radon concentrations, indicating the additional effects of localized flowpath heterogeneity overprinting on larger scale flowpath characteristics.  more » « less
Award ID(s):
1824820
PAR ID:
10478244
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Hydrological Processes
Volume:
37
Issue:
7
ISSN:
0885-6087
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. As the climate warms and dry periods become more extreme, shallow groundwater discharge is generally becoming a less reliable source of streamflow while deep groundwater discharge remains a more resilient source. The implications of shifts in the relative balance of shallow and deep groundwater discharge sources are profound in gaining streams. These different sources exert critical controls on stream temperature and water quality as influenced by legacy groundwater contaminant transport. Groundwater discharge flux rates over time were used for the inference of source groundwater characteristics to prominent riverbank groundwater discharge faces along the mainstem Farmington River, CT USA. To estimate groundwater discharge rates, we deployed sediment temperature loggers (iButton #DS1922L, Maxim Integrated, Inc., San Jose, CA, USA) in vertical profilers installed directly into mapped preferential groundwater discharge points across extensive riverbank discharge face features.Temperature data contained in this release were collected from June 24 to November 5, 2020, at 40 distinct discharge point riverbank locations, similar to those described by Barclay et al. (2022) and Briggs et al. (2022). Saturated sediment thermal conductivity and heat capacity were measured in-situ with a TEMPOS Thermal Property Analyzer (TEMPOS, Meter Group, Inc., Pullman, WA, USA) at multiple points across each riverbank discharge face to aid in estimating groundwater discharge flux rates. Barclay, J. R., Briggs, M. A., Moore, E. M., Starn, J. J., Hanson, A. E. H., & Helton, A. M. (2022). Where groundwater seeps: Evaluating modeled groundwater discharge patterns with thermal infrared surveys at the river-network scale. Advances in Water Resources, 160. https://doi.org/10.1016/j.advwatres.2021.104108 Briggs, M. A., Jackson, K. E., Liu, F., Moore, E. M., Bisson, A., & Helton, A. M. (2022). Exploring Local Riverbank Sediment Controls on the Occurrence of Preferential Groundwater Discharge Points. Water, 14(1). https://doi.org/10.3390/w14010011 
    more » « less
  2. We used spatial data from previously mapped preferential groundwater discharges throughout the Farmington River watershed in Connecticut and Massachusetts (https://doi.org/10.5066/P915E8JY) to guide water sample collection at known locations of groundwater discharging to surface water. In 2017 and 2019 - 2021, samples were collected during general river baseflow conditions (July ? November, less than 30.9 cms mean daily discharge (USGS gage 01189995, statistics 2010-2022) when the riverbank discharge points were exposed. We collected a suite of dissolved constituents and stable isotopes of water directly in the shallow saturated sediments of active points of discharge, and coincident stream chemical samples were also collected adjacent to locations of groundwater discharge. Data collected includes nutrients (NO3, NH4, Cl, SO4, PO4, dissolved organic carbon (DOC), and total nitrogen (TN)), greenhouse gases (CO2, CH4, and N2O), dissolved gases (N2, dissolved oxygen (DO)), conductivity, sediment characteristics, temperature, and spatial information. This dataset includes 2 main files: 1) Farmington_Chemistry_2017_2021.csv contains attribute information for each biogeochemical constituent collected at preferential groundwater discharges along the Farmington River network. 2)Farmington_Temporal_Cl_Rn_Iso_2020.csv contain attribute information for source characteristic data (Chloride, Radon, Isotope) collected at locations of repeat sampling at 5 groundwater seep faces along the Farmington River (Alsop and Rainbow Island). 
    more » « less
  3. This is the data used to create the published study of the same name published in the Journal of Hydrology in 2022 (10.1016/j.jconhyd.2022.104068). Shallow (<30 m) reducing groundwater commonly contains abundant dissolved arsenic (As) in Bangladesh. We hypothesize that dissolved As in iron (Fe)-rich groundwater discharging to rivers is trapped onto Fe(III)-oxyhydroxides which precipitate in shallow riverbank sediments under the influence of tidal fluctuations. Therefore, the goal of this study is to compare the calculated mass of sediment-bound As that would be sequestered from dissolved groundwater As that discharges through riverbanks of the Meghna River to the observed mass of As trapped within riverbank sediments. To calculate groundwater discharge, a Boussinesq aquifer analytical groundwater flow model was developed and constrained by cyclical seasonal fluctuations in hydraulic heads and river stages observed at three sites along a 13 km reach in central Bangladesh. At all sites, groundwater discharges to the river year-round but most of it passes through an intertidal zone created by ocean tides propagating upstream from the Bay of Bengal in the dry season. The annualized groundwater discharge per unit width at the three sites ranges from 173 to 891 m2/yr (average 540 m2/yr). Assuming that riverbanks have been stable since the Brahmaputra River avulsed far away from this area 200 years ago and dissolved As is completely trapped within riverbank sediments, the mass of accumulated sediment As can be calculated by multiplying groundwater discharge by ambient aquifer As concentrations measured in 1969 wells. Across all sites, the range of calculated sediment As concentrations in the riverbank is 78–849 mg/kg, which is higher than the observed concentrations (17–599 mg/kg). This discovery supports the hypothesis that the dissolved As in groundwater discharge to the river is sufficient to account for the observed buried deposits of As along riverbanks. 
    more » « less
  4. Groundwater discharge to rivers takes many forms, including preferential groundwater discharge points (PDPs) along riverbanks that are exposed at low flows, with multi-scale impacts on aquatic habitat and water quality. The physical controls on the spatial distribution of PDPs along riverbanks are not well-defined, rendering their prediction and representation in models challenging. To investigate the local riverbank sediment controls on PDP occurrence, we tested drone-based and handheld thermal infrared to efficiently map PDP locations along two mainstem rivers. Early in the study, we found drone imaging was better suited to locating tributary and stormwater inflows, which created relatively large water surface thermal anomalies in winter, compared to PDPs that often occurred at the sub-meter scale and beneath riparian tree canopy. Therefore, we primarily used handheld thermal infrared imaging from watercraft to map PDPs and larger seepage faces along 12-km of the fifth-order Housatonic River in Massachusetts, USA and 26-km of the Farmington River in Connecticut, USA. Overall, we mapped 31 riverbank PDPs along the Housatonic reach that meanders through lower permeability soils, and 104 PDPs along the Farmington reach that cuts through sandier sediments. Riverbank soil parameters extracted at PDP locations from the Soil Survey Geographic (SSURGO) database did not differ substantially from average bank soils along either reach, although the Farmington riverbank soils were on average 5× more permeable than Housatonic riverbank soils, likely contributing to the higher observed prevalence of PDPs. Dissolved oxygen measured in discharge water at these same PDPs varied widely, but showed no relation to measured sand, clay, or organic matter content in surficial soils indicating a lack of substantial near-surface aerobic reaction. The PDP locations were investigated for the presence of secondary bank structures, and commonly co-occurred with riparian tree root masses indicating the importance of localized physical controls on the spatial distribution of riverbank PDPs. 
    more » « less
  5. Estuaries are dynamic environments at the intersection of rivers and oceans. They provide important ecosystem services to coastal communities but have also been impacted by human activities. Understanding their physical functioning improves our ability to manage and restore estuaries. This work focuses on two aspects of estuarine dynamics that have not been a major focus of estuarine physical oceanography: cross-channel currents and water temperature. Cross-channel currents, although generally smaller in magnitude than the main along-channel flow, can influence the distribution of properties in estuaries. We investigate the patterns and magnitudes of cross-channel flow in the Hudson River estuary, finding distinct flow characteristics during spring and neap tides due to changes in stratification. We determine the contributions of channel curvature and Coriolis to these flows and quantify their impact on the along-channel momentum balance. We next explore water temperature dynamics in the Yaquina Bay estuary, using both observations and modeling to understand the relative importance of conservative mixing of river and ocean water and surface heat flux to temperature variability at seasonal, subseasonal, diurnal, and semidiurnal timescales. We collected two and a half years of temperature measurements throughout the estuary and found a strong seasonal cycle as well as subseasonal and high-frequency variability. We found that residence time was important in determining the importance of surface heat flux. High-frequency variability depended strongly on the along-channel temperature gradient and diurnal variance in particular was affected by the annual cycle in the phasing of the tides and the sun as well as heat exchange with intertidal sediments. To further investigate water temperature dynamics, we developed a three-dimensional hydrodynamic model of the Yaquina Bay estuary using the Finite Volume Community Ocean Model (FVCOM). We validated the model with observations and found that accurate boundary conditions, surface heat flux, and residence times were essential to reproducing observed temperature variability. We turned off surface heat flux in the model to show that surface heat flux contributes to temperature variability across timescales, with major impacts on the seasonal cycle as well as the along-channel temperature gradient. We also changed the estuary's boundary conditions, including sea level, ocean temperature, river temperature and discharge, and atmospheric conditions, to represent conditions at the end of the century under a high-emissions climate scenario. This simulation produced year-round estuary warming relative to the present, with consequences for estuarine habitat suitability for species such as oysters, eelgrass, and salmon. This work contributes to our improved understanding of the complex estuarine environment, which is central to maintaining the benefits of estuaries to humanity and the broader coastal ecosystem in a changing world. 
    more » « less