skip to main content

Title: Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approach

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
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Hydrological Processes
Page Range / eLocation ID:
p. 759-774
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. 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) plantation, (e) an approximately 10 y.o. secondary succession forest, (f) a 12 y.o. coffee plantation(Coffea canephora), (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.

    more » « less
  5. Claypan soils cover approximately 10 million acres across several states in the central United States. The soils are characterized by a highly impermeable clay layer within the profile that impedes water flow and root growth. While some claypan soils can be productive, they must be carefully managed to avoid reductions to crop productivity due to root restrictions, water, and nutrient limitations. Clay soils are usually resistant to erosion but may exacerbate erosion of the silt-loam topsoil. Soil production potential is the capacity of soil to produce at a given level (yield per acre). The productive capacity is tied to soil characteristics, which can be highly variable within a field. In this project, we have used imagery analysis to study the aerial images and terrain of fields during different productive times of the year to identify where soil samples should be collected for more discrete analysis. Soil samples provide valuable information; however, the amount of data obtained from a relatively small area within a field does not provide sufficient information to delineate the subsurface characteristics. To address the limitations of sampling, we have also employed the use of yield maps collected from commercial yield monitors on production-scale combines and surface electrical conductivity measurements (Sassenrath and Kulesza, 2017). Soil conductivity is a measurement of how well a representative volume of soil conducts electricity. Soil conductivity is a function of the soil clay content, moisture content, and other measurable soil properties (Kitchen et al., 2003); as such, it has become a valuable tool for mapping in-field variability. The main advantage of a soil conductivity measurement is that the entire surface of a field can be imaged. The disadvantage of a soil conductivity measurement is that data are only collected near the surface (10 – 30 inches) and the measurements are relative measurements. This means that the conductivity mappers can identify changes in soil properties, but they cannot directly tell researchers what caused these changes. Electrical resistivity tomography (ERT) is a popular near-surface geophysical measurement for geophysical and engineering applications. The term “near-surface” generally means down to around 30 feet in the subsurface. Electrical resistivity is the reciprocal measurement of electrical conductivity; therefore, both systems measure differences in the same soil properties. ERT measurements are different than surface electrical conductivity measurements because ERT collects a “slice” of data into the subsurface, as opposed to only changes at the surface area. Relative measurements, similar to those collected in an electrical conductivity survey, are collected; however, in ERT studies the data are mathematically inverted to yield the true electrical resistivity of the soil with depth. This allows an interpretation of the changing soil properties with depth to reduce the required amount of sampling. A disadvantage of an ERT survey is that the data acquisition is stationary so mapping an entire field is not feasible. We have used a coupled process of imagery and terrain analysis, yield maps, and electrical conductivity measurements to guide the locations of ERT surveys in this project (Tucker-Kulesza et al. 2017). 
    more » « less