An official website of the United States government
Abstract Subterranean estuaries (STEs) form at the land‐sea boundary where groundwater and seawater mix. These biogeochemically reactive zones influence groundwater‐borne nutrient concentrations and speciation prior to export via submarine groundwater discharge (SGD). We examined a STE located along the York River Estuary (YRE) to determine if SGD delivers dissolved inorganic nitrogen (DIN) and phosphorus (DIP) to the overlying water. We assessed variations in STE geochemical profiles with depth across locations, times, and tidal stages, estimated N removal along the STE flow path, measured hydraulic gradients to estimate SGD, and calculated potential nutrient fluxes. Salinity, dissolved oxygen (DO), DIN, and DIP varied significantly with depth and season (p < 0.05), but not location or tidal stage. Ammonium dominated the DIN pool deep in the STE. Moving toward the sediment surface, ammonium concentrations decreased as nitrate and DO concentrations increased, suggesting nitrification. Potential sediment N removal rates mediated by denitrification were <8 mmoles N m−2 d−1. The total groundwater discharge rate was 38 ± 11 L m−2 d−1; discharge followed tidal and seasonal patterns. Net SGD nutrient fluxes were 0.065–3.2 and 0.019–0.093 mmoles m−2 d−1for DIN and DIP, respectively. However, microbial N removal in the STE may attenuate 0.58% to >100% of groundwater DIN. SGD fluxes were on the same order of magnitude as diffusive benthic fluxes but accounted for <10% of the nutrients delivered by fluvial advection in the YRE. Our results indicate the importance of STE biogeochemical transformations to SGD flux estimations and their role in coastal eutrophication and nutrient dynamics.
more »
« less
Improved Prediction of Management‐Relevant Groundwater Discharge Characteristics Throughout River Networks
Abstract Groundwater discharge zones connect aquifers to surface water, generating baseflow and serving as ecosystem control points across aquatic ecosystems. The influence of groundwater discharge on surface flow connectivity, fate and transport of contaminants and nutrients, and thermal habitat depends strongly on hydrologic characteristics such as the spatial distribution, age, and depth of source groundwater flow paths. Groundwater models have the potential to predict spatial discharge characteristics within river networks, but models are often not evaluated against these critical characteristics and model equifinality with respect to discharge processes is a known challenge. We quantify discharge characteristics across a suite of groundwater models with commonly used frameworks and calibration data. We developed a base model (MODFLOW‐NWT) for a 1,570‐km2watershed in the northeastern United States and varied the calibration data, control of river‐aquifer exchange directionality, and resolution. Most models (n = 11 of 12) fit similarly to calibration metrics, but patterns in discharge location, flow path depth, and subsurface travel time varied substantially. We found (1) a 15% difference in the percent of discharge going to first‐order streams, (2) threefold variations in flow path depth, and (3) sevenfold variations in the subsurface travel times among the models. We recalibrated three models using a synthetic discharge location data set. Calibration with discharge location data reduced differences in simulated discharge characteristics, suggesting an approach to improved equifinality based on widespread field‐based mapping of discharge zones. Our work quantifying variation across common modeling approaches is an important step toward characterizing and improving predictions of groundwater discharge characteristics.
more »
« less
- Award ID(s):
- 1824820
- PAR ID:
- 10448344
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 56
- Issue:
- 10
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Most terrestrial allochthonous organic matter enters river networks through headwater streams during high flow events. In headwaters, allochthonous inputs are substantial and variable, but become less important in streams and rivers with larger watersheds. As allochthonous dissolved organic matter (DOM) moves downstream, the proportion of less aromatic organic matter with autochthonous characteristics increases. How environmental factors converge to control this transformation of DOM at a continental scale is less certain. We hypothesized that the amount of time water has spent travelling through surface waters of inland systems (streams, rivers, lakes, and reservoirs) is correlated to DOM composition. To test this hypothesis, we used established river network scaling relationships to predict relative river network flow‐weighted travel time (FWTT) of water for 60 stream and river sites across the contiguous United States (3090 discrete samples over 10 water years). We estimated lentic contribution to travel times with upstream in‐network lake and reservoir volume. DOM composition was quantified using ultraviolet and visible absorption and fluorescence spectroscopy. A combination of FWTT and lake and reservoir volume was the best overall predictor of DOM composition among models that also incorporated discharge, specific discharge, watershed area, and upstream channel length. DOM spectral slope ratio (R2 = 0.77) and Freshness Index (R2 = 0.78) increased and specific ultraviolet absorbance at 254 nm (R2 = 0.68) and Humification Index (R2 = 0.44) decreased across sites as a function of FWTT and upstream lake volume. This indicates autochthonous‐like DOM becomes continually more dominant in waters with greater FWTT. We assert that river FWTT can be used as a metric of the continuum of DOM composition from headwaters to rivers. The nature of the changes to DOM composition detected suggest this continuum is driven by a combination of photo‐oxidation, biological processes, hydrologically varying terrestrial subsidies, and aged groundwater inputs.more » « less
-
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
-
Abstract Coupled groundwater flow and heat transport within hyporheic zones extensively affect water, energy, and solute exchange with surrounding sediments. The local and cumulative implications of this tightly coupled process strongly depend on characteristics of drivers (i.e., discharge and temperature of the water column) and modulators (i.e., hydraulic and thermal properties of the sediment). With this in mind, we perform a systematic numerical analysis of hyporheic responses to understand how the temporal variability of river discharge and temperature affect flow and heat transport within hyporheic zones. We identify typical time series of river discharge and temperature from gauging stations along the headwater region of Mississippi River Basin, which are characterized by different degrees of flow alteration, to drive a physics‐based model of the hyporheic exchange process. Our modeling results indicate that coupled groundwater flow and heat transport significantly affects the dynamic response of hyporheic zones, resulting in substantial differences in exchange rates and characteristic time scales of hyporheic exchange processes. We also find that the hyporheic zone dampens river temperature fluctuations increasingly with higher frequency of temperature fluctuations. This dampening effect depends on the system transport time scale and characteristics of river discharge and temperature variability. Furthermore, our results reveal that the flow alteration reduces the potential of hyporheic zones to act as a temperature buffer and hinders denitrification within hyporheic zones. These results have significant implications for understanding the drivers of local variability in hyporheic exchange and the implications for the development of thermal refugia and ecosystem functioning in hyporheic zones.more » « less
