skip to main content


Title: Integrating observations and models to determine the effect of seasonally frozen ground on hydrologic partitioning in alpine hillslopes in the Colorado Rocky Mountains, USA
Abstract

This study integrated spatially distributed field observations and soil thermal models to constrain the impact of frozen ground on snowmelt partitioning and streamflow generation in an alpine catchment within the Niwot Ridge Long‐Term Ecological Research site, Colorado, USA. The study area was comprised of two contrasting hillslopes with notable differences in topography, snow depth and plant community composition. Time‐lapse electrical resistivity surveys and soil thermal models enabled extension of discrete soil moisture and temperature measurements to incorporate landscape variability at scales and depths not possible with point measurements alone. Specifically, heterogenous snowpack thickness (~0–4 m) and soil volumetric water content between hillslopes (~0.1–0.45) strongly influenced the depths of seasonal frost, and the antecedent soil moisture available to form pore ice prior to freezing. Variable frost depths and antecedent soil moisture conditions were expected to create a patchwork of differing snowmelt infiltration rates and flowpaths. However, spikes in soil temperature and volumetric water content, as well as decreases in subsurface electrical resistivity revealed snowmelt infiltration across both hillslopes that coincided with initial decreases in snow water equivalent and early increases in streamflow. Soil temperature, soil moisture and electrical resistivity data from both wet and dry hillslopes showed that initial increases in streamflow occurred prior to deep soil water flux. Temporal lags between snowmelt infiltration and deeper percolation suggested that the lateral movement of water through the unsaturated zone was an important driver of early streamflow generation. These findings provide the type of process‐based information needed to bridge gaps in scale and populate physically based cryohydrologic models to investigate subsurface hydrology and biogeochemical transport in soils that freeze seasonally.

 
more » « less
Award ID(s):
1637686 2012730
NSF-PAR ID:
10367223
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Hydrological Processes
Volume:
35
Issue:
10
ISSN:
0885-6087
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Soil moisture is a key control on runoff generation and biogeochemical processes on hillslopes. Precipitation events can evoke different soil moisture responses with depth through the soil profile, and responses can differ among landscape positions along a hillslope. We sought to elucidate the nature of these responses by estimating changes in water content, response time between peak precipitation and peak soil moisture, and wetting front velocities for 43 storms at 45 locations on three adjacent hillslopes within a headwater catchment of the southern Appalachian Mountains (NC, USA). We used a multivariate modeling approach to quantify the relative influences and the predictability of soil moisture responses by a combination of landscape and storm characteristics. We quantified the lag correlations between hillslope mean soil moisture and catchment runoff to demonstrate how storm properties and hillslope‐scale characteristics may influence runoff at the catchment outlet. Soil moisture responses varied widely, and no consistent patterns were observed among response metrics laterally or vertically along hillslopes. In contrast to other studies, we found that the relative influence of hillslope properties and storm characteristics varied with soil moisture responses and during storms. Antecedent conditions and storm depths influenced the strength of lag correlations between soil moisture and runoff, whereas storm mean intensity was correlated with the lag times. These results highlight the utility of intensive observations for characterizing heterogeneity in soil moisture responses, suggesting, among other things, a need for better representation of the subsurface processes in rainfall‐runoff models. Identifying the relative importance of drivers can be beneficial in building parsimonious hydrological models.

     
    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

    Climate change threatens water resources in snowmelt‐dependent regions by altering the fraction of snow and rain and spurring an earlier snowmelt season. The bulk of hydrological research has focused on forecasting response in streamflow volumes and timing to a shrinking snowpack; however, the degree to which subsurface storage offsets the loss of snow storage in various alpine geologic settings, i.e. the hydrogeologic buffering capacity, is still largely unknown. We address this research need by assessing the affects of climate change on storage and runoff generation for two distinct hydrogeologic settings present in alpine systems: a low storage granitic and a greater storage volcanic hillslope. We use a physically based integrated hydrologic model fully coupled to a land surface model to run a base scenario and then three progressive warming scenarios, and account for the shifts in each component of the water budget. For hillslopes with greater water retention, the larger storage volcanic hillslope buffered streamflow volumes and timing, but at the cost of greater reductions in groundwater storage relative to the low storage granite hillslope. We found that the results were highly sensitive to the unsaturated zone retention parameters, which in the case of alpine systems can be a mix of matrix or fracture flow. The presence of fractures and thus less retention in the unsaturated zone significantly decreased the reduction in recharge and runoff for the volcanic hillslope in climate warming scenarios. This approach highlights the importance of incorporating physically based subsurface flow in to alpine hydrology models, and our findings provide ways forward to arrive at a conceptual model that is both consistent with geology and hydrologic principles. Copyright © 2016 John Wiley & Sons, Ltd.

     
    more » « less
  4. Abstract

    This paper presents a top–down approach for soil moisture and sap flux sampling design with the goal of understanding ecohydrologic response to interannual climate variation in the rain–snow transition watersheds. The design is based on a priori estimates of soil moisture and transpiration patterns using a physical distributed model, Regional Hydro‐Ecologic Simulation System (RHESSys). RHESSys was initially calibrated with existing snow depth and streamflow data. Calibrated model estimates of seasonal trajectories of snowmelt, root‐zone soil moisture storage, and transpiration were used to develop five hydrologic similarity indicators and map these at (30 m) patch scale across the study watershed. The partitioning around medoids‐clustering algorithm was then used to define six distinctive spatially explicit clusters based on the five hydrologic similarity indictors. A representative site within each cluster was identified for sampling. For each site, soil moisture sensors were installed at the 30‐ and 90‐cm depths and at the five soil pits and a sap flux sensor at the averaged‐size white fir tree for each site. The model‐based cluster analysis suggests that the elevation gradient and topographically driven flow drainage patterns are the dominant drivers of spatial patterns of soil moisture and transpiration. The comparison of model‐based calculated hydrological similarity indicators with measured‐data‐based values shows that spatial patterns of field‐sampled soil moisture data typically fell within uncertainty bounds of model‐based estimates for each cluster. There were however several notable exceptions. The model failed to capture the soil moisture and sap flux dynamics in a riparian zone site and in a site where lateral subsurface flow may not follow surface topography. Results highlight the utility of using a hypothesis driven sampling strategy, based on a physically based model, for efficiently providing new information that can drive both future measurements and strategic refinements to model inputs, parameters, or structure that might reduce these errors. Future research will focus on strategies for using of finer scale representations of microclimate, topography, vegetation, and soil properties to improve models.

     
    more » « less
  5. Abstract

    Frozen ground influences subsurface hydrology by reducing soil permeability, impeding infiltration, and inducing surface runoff. As the seasonal snowpack strongly controls the evolution of the ground thermal regime, it is necessary to represent the snowpack in models simulating frozen ground and its hydrologic consequences. Conventional understanding of the relationship between snowcover and frozen ground considers the snowpack as an insulator, shielding the ground from cold winter temperatures and thus inhibiting frozen ground. However, observations and modeling at Niwot Ridge, a seasonally snow‐covered alpine catchment in the headwaters of the Boulder Creek watershed, illustrate that snowpack cold content can promote and preserve frozen ground beneath the snowpack during the snowmelt season, unlike bare ground patches that thaw relatively rapidly due to exposure to solar radiation and warm spring temperatures. We present results from a coupled snowpack and subsurface thermo‐hydrologic model including freeze‐thaw processes that captures these contrasts reasonably well. Our results suggest that the cold snowpack at this site acts as a heat sink and promotes frozen ground, in contrast to the conventional understanding of the snowpack as an insulator. In the subalpine zone, infiltration simulations show that shallow freezing beneath snow‐covered ground has a much stronger effect on infiltration than shallow freezing beneath bare ground because the soil beneath the snow remained frozen during snowmelt, whereas bare patches thawed by the time they received excess snowmelt run‐on.

     
    more » « less