Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first‐order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated‐soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.
Mountain block systems are critical to water resources and have been heavily studied and modeled in recent decades. However, due to lack of field data, there is little consistency in how models represent the mountain block subsurface. While there is a large body of research on subsurface heterogeneity, few studies have evaluated the effect that common conceptual choices modelers make in mountainous systems have on simulated hydrology. Here we simulate the hydrology of a semi‐idealized headwater catchment using six common conceptual models of the mountain block subsurface. These scenarios include multiple representations of hydraulic conductivity decaying with depth, changes in soil depth with topography, and anisotropy. We evaluate flow paths, discharge, and water tables to quantify the impact of subsurface conceptualization on hydrologic behavior in three dimensions. Our results show that adding higher conductivity layers in the shallow subsurface concentrates flow paths near the surface and increases average saturated flow path velocities. Increasing heterogeneity by adding additional layers or introducing anisotropy increases the variance in the relationship between the age and length of saturated flow paths. Discharge behavior is most sensitive to heterogeneity in the shallow subsurface layers. Water tables are less sensitive to layering than they are to the overall conductivity in the domain. Anisotropy restricts flow path depths and controls discharge from storage but has little effect on governing runoff. Differences in the response of discharge, water table depth, and residence time distribution to subsurface representation highlight the need to consider model applications when determining the level of complexity that is needed.
more » « less- Award ID(s):
- 1806383
- PAR ID:
- 10448445
- 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
-
more » « less
Abstract -
How does the physical and chemical structure of the Critical Zone (CZ), defined as the zone from treetops to the bottom of groundwater, govern its hydro-biogeochemical functioning? Multiple lines of evidence from past and newly emerging research have prompted the shallow and deep partitioning concentration-discharge (C-Q) hypothesis. The hypothesis states that in-stream C-Q relationships are shaped by distinct source waters from flow paths at different depths. Base flows are often dominated by deep groundwater and mostly reflect groundwater chemistry, whereas high flows are often dominated by shallow soil water and thus mostly reflect soil water chemistry. The contrasts between shallow soil water versus deeper groundwater chemistry shape stream solute export patterns. In this context, the vertical connectivity that regulates the shallow and deep flow partitioning is essential in determining chemical contrasts, biogeochemical reaction rates in soils and parent rocks, and ultimately solute export patterns. This talk will highlight insights gleaned from multiple lines of recent studies that include collation of water chemistry data from soils, rocks, and streams in intensively monitored watersheds, meta-analysis of stream chemistry data at the continental scale, and integrated reactive transport modeling at the hillslope and watershed scales. The hypothesis underscores the importance of subsurface vertical structure and connectivity relative to the extensively studied horizontal connectivity. It also alludes to the potential of using streams as mirrors for subsurface water chemistry, and the potential of using C-Q relationships to infer flow paths and biogeochemical reaction rates and the response of earth’s subsurface to climate and human perturbations. Broadly, this simple conceptual framework links CZ subsurface structure to its functioning under diverse climate, geology, and land cover conditions.more » « less
-
more » « less
Abstract Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds.
-
more » « less
Abstract Dissolved organic and inorganic carbon (DOC and DIC) influence water quality, ecosystem health, and carbon cycling. Dissolved carbon species are produced by biogeochemical reactions and laterally exported to streams via distinct shallow and deep subsurface flow paths. These processes are arduous to measure and challenge the quantification of global carbon cycles. Here we ask: when, where, and how much is dissolved carbon produced in and laterally exported from the subsurface to streams? We used a catchment‐scale reactive transport model, BioRT‐HBV, with hydrometeorology and stream carbon data to illuminate the “invisible” subsurface processes at Sleepers River, a carbonate‐based catchment in Vermont, United States. Results depict a conceptual model where DOC is produced mostly in shallow soils (3.7 ± 0.6 g/m2/yr) and in summer at peak root and microbial respiration. DOC is flushed from soils to the stream (1.0 ± 0.2 g/m2/yr) especially during snowmelt and storms. A large fraction of DOC (2.5 ± 0.2 g/m2/yr) percolates to the deeper subsurface, fueling deep respiration to generate DIC. DIC is exported predominantly from the deeper subsurface (7.1 ± 0.4 g/m2/yr, compared to 1.3 ± 0.3 g/m2/yr from shallow soils). Deep respiration reduces DOC and increases DIC concentrations at depth, leading to commonly observed DOC flushing (increasing concentrations with discharge) and DIC dilution patterns (decreasing concentrations with discharge). Surprisingly, respiration processes generate more DIC than weathering in this carbonate‐based catchment. These findings underscore the importance of vertical connectivity between the shallow and deep subsurface, highlighting the overlooked role of deep carbon processing and export.
-
more » « less
Abstract A portion of water not consumed by crops during flood irrigation can flow back across the surface or through the subsurface to adjacent surface water bodies and streams as return flow. Few studies have directly addressed subsurface processes governing return flow and the importance of structural complexity on hydrologic process representation. It is challenging to measure and model these subsurface flow paths using traditional hydrologic observations. In this study, we assess the impact of subsurface structural complexity on vadose zone flow representation in a two‐dimensional transport model by varying structural complexity derived from background geophysical data. We assessed four model structures each with three soil types of homogeneous hydrologic properties, two of which were evaluated with and without an anisotropy factor. Wetting front arrival times, derived from time‐lapse electrical resistivity measurements during flood irrigation field experiments, were used to evaluate the different representations of soil profile structures. These data indicated both vertical and lateral preferential flow in the subsurface during flood irrigation. Inclusion of anisotropy in the saturated hydraulic conductivity field improved the ability to model subsurface hydrologic behavior when flow processes shifted from uniform to heterogeneous flow, as occurs with lateral subsurface return flow under flood irrigation driven by a large pressure gradient. This reduced the need for detailed spatial discretization to represent these observed subsurface flow processes. The resulting simple three‐layer model structure was better able to model both the vertical and lateral flow processes than a more complex geospatial structure, suggesting that overinterpretation of smoothed inverted profiles could lead to misrepresentation of the subsurface structure.
