Abstract In drylands, runoff during storms redistributes water and nutrients from bare soil areas to vegetated patches, subsidizing vegetation with additional resources. The extent of this redistribution depends on the interplay between surface roughness and permeability; greater permeability in vegetated patches promotes run‐on to vegetation, but greater surface roughness diverts runoff, producing tortuous flow paths that bypass vegetation. Here, this interplay is examined in virtual experiments using the 2D Saint Venant Equations to measure runoff connectivity. Flowpaths are delineated using tracers advected by the flow. Distances between tracer sources and sinks along flowpaths measure hydrologic connectivity at two lengthscales: connectivity to the hillslope outlet and within‐slope source‐sink connectivity. Differences between these connectivity lengthscales indicate how flow may “by‐pass” vegetated patches within hillslopes. At the hillslope scale, a derived power‐law relation between the runoff coefficient and outlet connectivity describes hillslope water losses, providing a foundation for identifying landscapes likely to shed water.
more »
« less
Hillslope to channel hydrologic connectivity in a dryland ecosystem
Abstract Hydrologic connectivity refers to the processes and thresholds leading to water transport across a landscape. In dryland ecosystems, runoff production is mediated by the arrangement of vegetation and bare soil patches on hillslopes and the properties of ephemeral channels. In this study, we used runoff measurements at multiple scales in a small (4.67 ha) mixed shrubland catchment of the Chihuahuan Desert to identify controls on and thresholds of hillslope‐channel connectivity. By relating short‐ and long‐term hydrologic records, we also addressed whether observed changes in outlet discharge since 1977 were linked to modifications in hydrologic connectivity. Hillslope runoff production was controlled by the maximum rainfall intensity occurring in a 30‐min interval (I30), with small‐to‐negligible effects of antecedent surface soil moisture, vegetation cover, or slope aspect. AnI30threshold of nearly 10 mm/h activated runoff propagation from the shrubland hillslopes and through the main ephemeral channel, whereas anI30threshold of about 16 mm/h was required for discharge from the catchment outlet. Since storms rarely exceedI30, full hillslope‐channel connectivity occurs infrequently in the mixed shrubland, leading to <2% of the annual precipitation being converted into outlet discharge. Progressive decreases in outlet discharge since 1977 could not be explained by variations in precipitation metrics, includingI30, or the process of woody plant encroachment. Instead, channel modifications from the buildup of sediment behind measurement flumes may have increased transmission losses and reduced outlet discharge. Thus, alterations in channel properties can play an important role in the long‐term (45‐year) variations of rainfall–runoff dynamics of small desert catchments.
more »
« less
- Award ID(s):
- 2025166
- PAR ID:
- 10477062
- Publisher / Repository:
- Wiley Online Library
- Date Published:
- Journal Name:
- Ecosphere
- Volume:
- 14
- Issue:
- 11
- ISSN:
- 2150-8925
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Debris flows are powered by sediment supplied from steep hillslopes where soils are often patchy and interrupted by bare‐bedrock cliffs. The role of patchy soils and cliffs in supplying sediment to channels remains unclear, particularly surrounding wildfire disturbances that heighten debris‐flow hazards by increasing sediment supply to channels. Here, we examine how variation in soil cover on hillslopes affects sediment sizes in channels surrounding the 2020 El Dorado wildfire, which burned debris‐flow prone slopes in the San Bernardino Mountains, California. We focus on six headwater catchments (<0.1 km2) where hillslope sources ranged from a continuous soil mantle to 95% bare‐bedrock cliffs. At each site, we measured sediment grain size distributions at the same channel locations before and immediately following the wildfire. We compared results to a mixing model that accounts for three distinct hillslope sediment sources distinguished by local slope thresholds. We find that channel sediment in fully soil‐mantled catchments reflects hillslope soils (D50 = 0.1–0.2 cm) both before and after the wildfire. In steeper catchments with cliffs, channel sediment is consistently coarse prior to fire (D50 = 6–32 cm) and reflects bedrock fracture spacing, despite cliffs representing anywhere from 5% to 95% of the sediment source area. Following the fire, channel sediment size reduces most (5‐ to 20‐fold) in catchments where hillslope sources are predominantly soil covered but with patches of cliffs. The abrupt fining of channel sediment is thought to facilitate postfire debris‐flow initiation, and our results imply that this effect is greatest where bare‐bedrock cliffs are present but not dominant. A patchwork of bare‐bedrock cliffs is common in steeplands where hillslopes respond to channel incision by landsliding. We show how local slope thresholds applied to such terrain aid in estimating sediment supply conditions before two destructive debris flows that eventually nucleated in these study catchments in 2022.more » « less
-
Abstract. Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of three shale-underlain headwater catchments located in Pennsylvania, USA (the forested Shale Hills Critical Zone Observatory), and Wales, UK (the peatland-dominated Upper Hafren and forest-dominated Upper Hore catchments in the Plynlimon forest), dissimilar concentration–discharge (C–Q) behaviors are best explained by contrasting landscape distributions of soil solution chemistry – especially dissolved organic carbon (DOC) – that have been established by patterns of vegetation and soil organic matter (SOM). Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Furthermore, concentration–discharge relationships of non-chemostatic solutes changed following tree harvest in the Upper Hore catchment in Plynlimon, while no changes were observed for chemostatic solutes, underscoring the role of vegetation in regulating the concentrations of certain elements in the stream. These results indicate that differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where SOM is dominantly in lowlands (e.g., Shale Hills), we infer that non-chemostatic elements associated with organic matter are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), these non-chemostatic elements are released later during rainfall events. The distribution of SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.more » « less
-
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
-
Abstract. Catchment-scale hydrological studies on drylands are lacking because of thescarcity of consistent data: observations are often available at the plotscale, but their relevance for the catchment scale remains unclear. Adatabase of 24 years of stream gauge discharge and homogeneoushigh-resolution radar data over the eastern Mediterranean allows us to describe the properties of floods over catchments spanning from desert toMediterranean climates, and we note that the data set is mostly of moderateintensity floods. Comparing two climatic regions, desert and Mediterranean,we can better identify specific rainfall-runoff properties. Despite the large differences in rainfall forcing between the two regions, the resulting unitpeak discharges and runoff coefficients are comparable. Rain depth andantecedent conditions are the most important properties to shape floodresponse in Mediterranean areas. In deserts, instead, storm core propertiesdisplay a strong correlation with unit peak discharge and, to a lesser extent,with runoff coefficient. In this region, an inverse correlation with meancatchment annual precipitation suggests also a strong influence of localsurface properties. Preliminary analyses suggest that floods in catchmentswith wet headwater and dry lower section are more similar to desertcatchments, with a strong influence of storm core properties on runoffgeneration.more » « less