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.
Time‐lapse electrical resistivity tomography (ERT) data are increasingly used to inform the hydrologic dynamics of mountainous environments at the hillslope scale. Despite their popularity and recent advancements in hydrogeophysical inversion methods, few studies have shown how time‐lapse ERT data can be used to determine hydraulic parameters of subsurface water flow models. This study uses synthetic and field‐collected, hillslope‐scale, time‐lapse ERT data to determine subsurface hydraulic properties of a two‐layer, physics‐based, 2‐D vertical flow model with predefined layer and boundary locations. Uncoupled and coupled hydrogeophysical inversion methods are combined with a fine‐earth fraction optimization scheme to reduce the number of parameters needing calibration and interpret the influence of the hydraulic parameters on the hydrologic model predictions. Inversions of synthetic ERT data recover the prescribed fine‐earth fraction bulk density to within 0.1 g cm−3. Field‐collected ERT data from a mountain hillslope result in hydrologic model dynamics that are consistent with previous studies and measured water content data but struggle to capture measured groundwater levels. The uncoupled hydrogeophysical inversion method is more sensitive to changes in hydraulic parameter values of the lower hydrologic model layer than the coupled hydrogeophysical inversion method. Time series of minimum objective function value simulations indicate that periodically collected ERT data may recover hydraulic parameters to a similar level of uncertainty as daily ERT data. Using simple hydrologic model domains within hydrogeophysical inversions shows promise for providing reasonable hydrologic predictions while maintaining relatively simple calibration schemes and should be explored further in future studies.
more » « less- Award ID(s):
- 1818550
- NSF-PAR ID:
- 10419210
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 4
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approachmore » « less
Abstract The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min−1/2for sites with a subsurface dominated by saprolitic material to 0.156 m min−1/2for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.
-
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.
-
more » « less
Abstract The hydrologic dynamics and geomorphic evolution of watersheds are intimately coupled—runoff generation and water storage are controlled by topography and properties of the surface and subsurface, while also affecting the evolution of those properties over geologic time. However, the large disparity between their timescales has made it difficult to examine interdependent controls on emergent hydrogeomorphic properties, such as hillslope length, drainage density, and extent of surface saturation. In this study, we develop a new model coupling hydrology and landscape evolution to explore how runoff generation affects long‐term catchment evolution, and analyze numerical results using a nondimensional scaling framework. We focus on hydrologic processes dominating in humid climates where storm runoff primarily arises from shallow subsurface flow and from precipitation on saturated areas. The model solves hydraulic groundwater equations to predict the water‐table elevation given prescribed, constant groundwater recharge. Water in excess of the subsurface capacity for transport becomes overland flow, which generates shear stress on the surface and may detach and transport sediment. This affects the landscape form that in turn affects runoff generation. We show that (a) four dimensionless parameters describe the possible steady state landscapes that coevolve under steady recharge; (b) hillslope length increases with increasing transmissivity relative to the recharge rate; (c) three topographic metrics—steepness index, Laplacian curvature, and topographic index—together provide a basis for interpreting landscapes that have coevolved with runoff generated via shallow subsurface flow. Finally we discuss the possibilities and limitations for quantitative comparisons between the model results and real landscapes.
-
The belowground architecture of the critical zone (CZ) consists of soil and rock in various stages of weathering and wetness that acts as a medium for biological growth, mediates chemical reactions, and controls partitioning of hydrologic fluxes. Hydrogeophysical imaging provides unique insights into the geometries and properties of earth materials that are present in the CZ and beyond the reach of direct observation beside sparse wellbores. An improved understanding of CZ architecture can be achieved by leveraging the geophysical measurements of the subsurface. Creating categorical models of the CZ is valuable for driving hydrologic models and comparing belowground architectures between different sites to interpret weathering processes. The CZ architecture is revealed through a novel comparison of hillslopes by applying facies classification in the elastic-electric domain driven by surface-based hydrogeophysical measurements. Three pairs of hillslopes grouped according to common geologic substrates — granite, volcanic extrusive, and glacially altered — are classified by five different hydrofacies classes to reveal the relative wetness and weathering states. The hydrofacies classifications are robust to the choice of initial mean values used in the classification and noncontemporaneous timing of geophysical data acquisition. These results will lead to improved interdisciplinary models of CZ processes at various scales and to an increased ability to predict the hydrologic timing and partitioning. Beyond the hillslope scale, this enhanced capability to compare CZ architecture can also be exploited at the catchment scale with implications for improved understanding of the link between rock weathering, hydrochemical fluxes, and landscape morphology.more » « less
