An official website of the United States government
Abstract Spatially integrated transport models have been applied widely to model hydrologic transport. However, we lack simple and process‐based theoretical tools to predict the transport closures—transit time distributions (TTDs) and StorAge Selection (SAS) functions. This limits our ability to infer characteristics of hydrologic systems from tracer observations and to make first‐order estimates of SAS functions in catchments where no tracer data is available. Here we present a theoretical framework linking TTDs and SAS functions to hydraulic groundwater theory at the hillslope scale. For hillslopes where the saturated hydraulic conductivity declines exponentially with depth, analytical solutions for the closures are derived that can be used as hypotheses to test against data. In the simplest form, the hillslope SAS function resembles a uniform or exponential distribution (corresponding to flow pathways in the saturated zone) offset from zero by the storage in the unsaturated zone that does not contribute to discharge. The framework is validated against nine idealized virtual hillslopes constructed using a 2‐D Richards equation‐based model, and against data from tracer experiments in two artificial hillslopes. Modeled internal age, life expectancy, and transit time structures reproduce theoretical predictions. The experimental data also support the theory, though further work is needed to account for the effects of time‐variability. The shape and tailing of TTDs and their power spectra are discussed. The theoretical framework yields several dimensionless numbers that can be used to classify hillslope scale flow and transport dynamics and suggests distinct water age structures for high or low Hillslope number.
more »
« less
Direct Observation of Hillslope Scale StorAge Selection Functions in Experimental Hydrologic Systems: Geomorphologic Structure and Preferential Discharge of Old Water
Abstract Spatially integrated water transport dynamics at the hillslope scale have rarely been observed directly, and underlying physical mechanisms of those dynamics are poorly understood. We present time‐variable transit time distributions and StorAge Selection (SAS) functions for a 28 days tracer experiment conducted at the Landscape Evolution Observatory, Biosphere 2, the University of Arizona, AZ, USA. The observed form of the SAS functions is concave, meaning that older water in the hillslope was preferentially discharged than younger water. The concavity is, in part, explained by the relative importance of advective and diffusive water dynamics and by the geomorphologic structure of the hillslopes. A simple numerical examination illustrates that, for straight plan‐shaped hillslopes, the saturated zone SAS function is concave when the hillslope Péclet (Pe) number is large (and thus when the advective water dynamics are more pronounced). We also investigated the effect of hillslope planform geometry on the saturated zone SAS function using a model and found that the more convergent the plan shape is, the more concave the SAS function is. A numerical examination indicates that the unsaturated zone SAS function is concave for straight and convergent hillslopes when the soil thickness is uniform. The concavity of those subcomponent SAS functions signifies that the hillslope scale SAS function is concave for straight or convergent plan shape hillslopes when the hillslope Pe number is high.
more »
« less
- PAR ID:
- 10443517
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 3
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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.more » « less
-
Amplified warming of high latitudes and rapid thaw of frozen ground threaten permafrost carbon stocks. The presence of permafrost modulates water infiltration and flow, as well as sediment transport, on soil-mantled slopes, influencing the balance of advective fluvial processes to diffusive processes on hillslopes in ways that are different from temperate settings. These processes that shape permafrost landscapes also impact the carbon stored on soil-mantled hillslopes via temperature, saturation, and slope stability such that carbon stocks and landscape morphometry should be closely linked. We studied 69,000 headwater basins between 25° and 90 °N to determine whether the thermal state of the soil sets the balance between hillslope (diffusive) and fluvial (advective) erosion processes, as evidenced by the density of the channel networks (i.e., drainage density) and the proportion of convex to concave topography (hillslopes and river valleys, respectively). Watersheds within permafrost regions have lower drainage densities than regions without permafrost, regardless of watershed glacial history, mean annual precipitation, and relief. We find evidence that advective fluvial processes are inhibited in permafrost landscapes compared to their temperate counterparts. Frozen soils likely inhibit channel development, and we predict that climate warming will lower incision thresholds to promote growth of the channel network in permafrost landscapes. By demonstrating how the balance of advective versus diffusive processes might shift with future warming, we gain insight into the mechanisms that shift these landscapes from sequestering to exporting carbon.more » « less
-
Abstract Both hydrological and geophysical data can be used to calibrate hillslope hydrologic models. However, these data often reflect hydrological dynamics occurring at disparate spatial scales. Their use as sole objectives in model calibrations may thus result in different optimum hydraulic parameters and hydrologic model behavior. This is especially true for mountain hillslopes where the subsurface is often heterogeneous and the representative elementary volume can be on the scale of several m3. This study explores differences in hydraulic parameters and hillslope‐scale storage and flux dynamics of models calibrated with different hydrological and geophysical data. Soil water content, groundwater level, and two time‐lapse electrical resistivity tomography (ERT) data sets (transfer resistance and inverted resistivity) from two mountain hillslopes in Wyoming, USA, are used to calibrate physics‐based surface–subsurface hydrologic models of the hillslopes. Calibrations are performed using each data set independently and all data together resulting in five calibrated parameter sets at each site. Model predicted hillslope runoff and internal hydrological dynamics vary significantly depending on the calibration data set. Results indicate that water content calibration data yield models that overestimate near‐surface water storage in mountain hillslopes. Groundwater level calibration data yield models that more reasonably represent hillslope‐scale storage and flux dynamics. Additionally, ERT calibration data yield models with reasonable hillslope runoff predictions but relatively poor predictions of internal hillslope dynamics. These observations highlight the importance of carefully selecting data for hydrologic model calibration in mountain environments. Poor selection of calibration data may yield models with limited predictive capability depending on modeling goals and model complexity.more » « less
-
Abstract How does hillslope structure (e.g., hillslope shape and permeability variation) regulate its hydro‐geochemical functioning (flow paths, solute export, chemical weathering)? Numerical reactive transport experiments and particle tracking were used to answer this question. Results underscore the first‐order control of permeability variations (with depth) on vertical connectivity (VC), defined as the fraction of water flowing into streams from below the soil zone. Where permeability decreases sharply and VC is low, >95% of water flows through the top 6 m of the subsurface, barely interacting with reactive rock at depth. High VC also elongates mean transit times (MTTs) and weathering rates. VC however is less of an influence under arid climates where long transit times drive weathering to equilibrium. The results lead to three working hypotheses that can be further tested.H1:The permeability variations with depth influence MTTs of stream water more strongly than hillslope shapes; hillslope shapes instead influence the younger fraction of stream water more.H2:High VC arising from high permeability at depths enhances weathering by promoting deeper water penetration and water‐rock interactions; the influence of VC weakens under arid climates and larger hillslopes with longer MTTs.H3:VC regulates chemical contrasts between shallow and deep waters (Cratio) and solute export patterns encapsulated in the power law slope b of concentration‐discharge (CQ) relationships.Higher VC leads to similar shallow versus deep water chemistry (Cratio∼1) and more chemostatic CQ patterns. Although supporting data already exist, these hypotheses can be further tested with carefully designed, co‐located modeling and measurements of soil, rock, and waters. Broadly, the importance of hillslope subsurface structure (e.g., permeability variation) indicate it is essential in regulating earth surface hydrogeochemical response to changing climate and human activities.more » « less
