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.
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- PAR ID:
- 10462804
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 33
- Issue:
- 5
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 759-774
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Flood irrigation is globally one of the most used irrigation methods. Typically, not all water that is applied during flood irrigation is consumed by plants or lost to evaporation. Return flow, the portion of applied water from flood irrigation that returns back to streams either via surface or subsurface flow, can constitute a large part of the water balance. Few studies have addressed the connection between vertical and lateral subsurface flows and its potential role in determining return flow pathways due to the difficulty in observing and quantifying these processes at plot or field scale. We employed a novel approach, combining induced polarization, time‐lapse electrical resistivity tomography, and time‐lapse borehole nuclear magnetic resonance, to identify flow paths and quantify changes in soil hydrological conditions under nonuniform application of flood irrigation water. We developed and tested a new method to track the wetting front in the subsurface using the full range of inverted resistivity values. Antecedent soil moisture conditions did not play an important role in preferential flow path activation. More importantly, boundaries between lithological zones in the soil profile were observed to control preferential flow pathways with subsurface run‐off occurring at these boundaries when saturation occurred. Using the new method to analyse time‐lapse resistivity measurements, we were able to track the wetting front and identify subsurface flow paths. Both uniform infiltration and preferential lateral flows were observed. Combining three geophysical methods, we documented the influence of lithology on subsurface flow processes. This study highlights the importance of characterizing the subsurface when the objective is to identify and quantify subsurface return flow pathways under flood irrigation.
-
more » « less
Abstract 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
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 Understanding land use/land cover (LULC) effects on tropical soil infiltration is crucial for maximizing watershed scale hydro‐ecosystem services and informing land managers. This paper reports results from a multiyear investigation of LULC effects on soil bulk infiltration in steep, humid tropical, and lowland catchments. A rainfall simulator applied water at measured rates on 2 × 6 m plots producing infiltration through structured, granulated, and macroporous Ferralsols in Panama's central lowlands. Time‐lapse electrical resistivity tomography (ERT) helped to visualize infiltration depth and bulk velocity. A space‐for‐time substitution methodology allowed a land‐use history investigation by considering the following: (a) a continuously heavy‐grazed cattle pasture, (b) a rotationally grazed traditional cattle pasture, (c) a 4‐year‐old (y.o.) silvopastoral system with nonnative improved pasture grasses and managed intensive rotational grazing, (d) a 7 y.o. teak (
Tectona grandis ( ), (g) an approximately 30 y.o. secondary succession forest, and (h) a >100 y.o. secondary succession forest. Within a land cover, unique plot sites totalled two at (a), (c), (d), (e), and (g); three at (b); and one at (f) and (h). Our observations confirmed measured infiltration scale dependency by comparing our 12 m2plot‐scale measurements against 8.9 cm diameter core‐scale measurements collected by others from nearby sites. Preferential flow pathways (PFPs) significantly increased soil infiltration capacity, particularly in forests greater than or equal to 10 y.o. Time‐lapse ERT observations revealed shallower rapid bulk infiltration and increased rapid lateral subsurface flow in pasture land covers when compared with forest land covers and highlighted how much subsurface flow pathways can vary within the Ferralsol soil class. Results suggest that LULC effects on PFPs are the dominant mechanism by which LULC affects throughfall partitioning, runoff generation, and flow pathways.Coffea canephora
