skip to main content

Title: The Climatic Water Balance and Topography Control Spatial Patterns of Atmospheric Demand, Soil Moisture, and Shallow Subsurface Flow

Catchment hydrometeorology and the organization of shallow subsurface flow are key drivers of active contributing areas and streamflow generation. However, understanding how the climatic water balance and complex topography contribute to these processes from hillslope to catchment scales remains difficult. We compared time series of vapor pressure deficits and soil moisture to the climatic water balance and topographic variables across six zero‐order catchments in the Lubrecht Experimental Forest (Montana, USA). We then evaluated how local hydrometeorology (volumetric water content and atmospheric vapor pressure deficit) affected the spatial occurrence of shallow subsurface flow. Generalized linear mixed model analysis revealed significant, temporally stable (monthly and seasonal average) patterns of hydrometeorology that can be predicted by the topographic wetness index and the dynamic climatic water deficit (CWD = potential evapotranspiration − actual evapotranspiration). Intracatchment patterns were significantly correlated to the topographic wetness index, while intercatchment patterns were correlated to spatiotemporal variance in the CWD during each time period. Spatial patterns of shallow subsurface flow were related to the hydrometeorological conditions of the site. We observed persistent shallow subsurface flow in convergent hillslope positions, except when a catchment was positioned in locations with high CWDs (low elevations and southerly aspects). Alternatively, we observed persistent subsurface flow across all hillslope positions (even 70‐m upslope from the hollow) when catchments were positioned in locations with especially low CWDs (northerly aspects and high elevations). These results highlight the importance of considering the superposition of the catchment‐scale climatic water balance and hillslope‐scale topography when characterizing hydrometeorology and shallow subsurface flow dynamics.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Page Range / eLocation ID:
p. 2370-2389
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Primary production is the entry point of energy and carbon into ecosystems, but modeling responses of primary production to “environmental stress” (i.e., reductions of primary production from nonoptimal environmental conditions) remains a key challenge and source of uncertainty in our understanding of Earth's carbon cycle. Here we develop an approach for estimating annual “environmental stress” from tree rings based on the proportion of the optimal diameter growth rate (from species‐specific allometric equations) that is realized in a given year. We assessed climatic, topographic, and soil drivers of environmental stress, as well as their interactions, using both empirical model experiments and linear mixed effect models. Climate gradients and interannual climate variability dominated spatial and temporal variability of environmental stress in much of the western United States, where the tree‐ring environmental stress index was positively correlated with antecedent climatic water balance (precipitation minus potential evapotranspiration) and negatively correlated with temperature and vapor pressure deficit. Excluding topographic and soil information from empirical models reduced their ability to capture spatial gradients in environmental stress, particularly in the eastern United States, where growth was not as strongly limited by climate. Mean climate conditions and topographic characteristics had significant interaction effects with the climatic water balance, indicating an increasing importance of winter moisture for warmer and drier sites and as elevation and topographic wetness index increased. These results suggest that including effects of antecedent climate (particularly in dry regions) and site topographic and soil characteristics could improve parameterization of environmental stress effects in primary production models.

    more » « less
  2. Concave hillslopes, also known as swales, play a pivotal role in facilitating the transfer of water and nutrients from hillslopes to riparian zones in headwater catchments by serving as pathways for surface and subsurface flow. Swales have long been recognized as hydrologic wet spots that contribute significantly to catchment runoff. However, long-term in situ observations have rarely confirmed this conventional understanding. In this study, we present a comprehensive analysis of three years of daily soil moisture measurements taken at multiple depths across 33 sites within a forested catchment to delineate key determinants of soil moisture variation at the catchment scale. Contrary to conventional understanding, our findings indicate that swales may represent areas of lower soil moisture, with those located on sunny hillslopes experiencing a higher frequency of dry conditions, amounting to 13.31% of the time during the study period, when compared to other topographical features such as slopes, valleys, and ridges. Topography alone accounted for less than half of the variation of soil moisture across the catchment. Whereas non-topographic factors (such as vegetation, soil porosity, soil temperature, and soil depth) strongly influenced soil moisture distribution in time and space. To better capture the nuances of soil moisture dynamics, we propose an enhancement to the Topographic Wetness Index (TWI) through the incorporation of non-topographic variables, which improves the model's representation of soil moisture variations. Our findings elucidate the complex interactions between topographic and non-topographic factors in shaping soil moisture dynamics at the catchment scale, offering valuable insights for refining process-based hydrological models and contributing to a deeper understanding of the spatial and temporal variability of soil moisture in headwater catchments. 
    more » « less
  3. 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.

    more » « less
  4. Abstract

    Groundwater flow direction within the critical zone of headwater catchments is often assumed to mimic land surface topographic gradients. However, groundwater hydraulic gradients are also influenced by subsurface permeability contrasts, which can result in variability in flow direction and magnitude. In this study, we investigated the relationship between shallow groundwater flow direction, surface topography, and the subsurface topography of low permeability units in a headwater catchment at the Hubbard Brook Experimental Forest (HBEF), NH. We continuously monitored shallow groundwater levels in the solum throughout several seasons in a well network (20 wells of 0.18–1.1 m depth) within the upper hillslopes of Watershed 3 of the HBEF. Water levels were also monitored in four deeper wells, screened from 2.4 to 6.9 m depth within glacial drift of the C horizon. We conducted slug tests across the well network to determine the saturated hydraulic conductivity (Ksat) of the materials surrounding each well. Results showed that under higher water table regimes, groundwater flow direction mimics surface topography, but under lower water table regimes, flow direction can deviate as much as 56 degrees from surface topography. Under these lower water table conditions, groundwater flow direction instead followed the topography of the top of the C horizon. The interquartile range ofKsatwithin the C horizon was two orders of magnitude lower than within the solum. Overall, our results suggest that the land surface topography and the top of the C horizon acted as end members defining the upper and lower bounds of flow direction variability. This suggests that temporal dynamics of groundwater flow direction should be considered when calculating hydrologic fluxes in critical zone and runoff generation studies of headwater catchments that are underlain by glacial drift.

    more » « less
  5. Abstract

    The seasonal evolution of the ground thermal regime in cold regions influences hydrologic flow paths, soil biogeochemistry, and hillslope geomorphology. In mountain environments, steep topography produces strong gradients in solar insolation, vegetation, and snowpack dynamics that lead to large differences in soil temperature over short distances, suggesting a need for high‐resolution, process‐based models that quantify the influence of topography. We present soil temperature and snow depth results from a coupled thermo‐hydrologic model compared to field observations from Gordon Gulch, a seasonally snow‐covered montane catchment in the Colorado Front Range in the Boulder Creek Critical Zone Observatory. The field site features two instrumented hillslopes with opposing aspects: Despite the persistent snowpack on the north‐facing slope, seasonally frozen ground is more prevalent there than the south‐facing slope, which experiences significantly higher incoming radiation that prevents the persistence of frozen ground. A novel modeling framework is developed by coupling a surface energy balance model incorporating solar radiation and snowpack processes to an existing subsurface model (PFLOTRAN‐ICE). The coupled model is used to reproduce strong aspect‐controlled differences in soil temperature and snow depth evident from observations during water years 2013–2016, including a higher incidence of frozen ground under the north‐facing slope. Representation of the snowpack and its insulating effects significantly improves soil temperature estimates on the north‐facing slope, particularly the duration of soil freezing in the spring, which is underestimated by 1–2 months without including the snowpack.

    more » « less