Linking quickflow response to subsurface state can improve our understanding of runoff processes that drive emergent catchment behaviour. We investigated the formation of non‐linear quickflows in three forested headwater catchments and also explored unsaturated and saturated storage dynamics, and likely runoff generation mechanisms that contributed to threshold formation. Our analyses focused on two reference watersheds at the Coweeta Hydrologic Laboratory (CHL) in western North Carolina, USA, and one reference watershed at the Susquehanna Shale Hills Critical Zone Observatory (SHW) in Central Pennsylvania, USA, with available hourly soil moisture, groundwater, streamflow, and precipitation time series over several years. Our study objectives were to characterise (a) non‐linear runoff response as a function of storm characteristics and antecedent conditions, (b) the critical levels of shallow unsaturated and saturated storage that lead to hourly flow response, and (c) runoff mechanisms contributing to rapidly increasing quickflow using measurements of soil moisture and groundwater. We found that maximum hourly rainfall did not significantly contribute to quickflow production in our sites, in contrast to prior studies, due to highly conductive forest soils. Soil moisture and groundwater dynamics measured in hydrologically representative areas of the hillslope showed that variable subsurface states could contribute to non‐linear runoff behaviour. Quickflow generation in watersheds at CHL were dominated by both saturated and unsaturated pathways, but the relative contributions of each pathway varied between catchments. In contrast, quickflow was almost entirely related to groundwater fluctuations at SHW. We showed that co‐located measurements of soil moisture and groundwater supplement threshold analyses providing stronger prediction and understanding of quickflow generation and indicate dominant runoff processes.
Soils derived from different lithologies and their controls on preferential flow remain underexplored in forested landscapes. In the same lithology, the propensity for preferential flow occurrence at different hillslope positions also remains largely elusive. By utilizing a soil moisture response time method, we compared preferential flow occurrence between a shale site (Shale Hills, silt loam soils) and a sandstone site (Garner Run, sandy loam soils) at four hillslope positions: ridge‐top, North‐ and South‐facing mid‐slopes and toe slope, for over 2 years. The catchments are neighbouring and covered by temperate forest. For the four hillslope positions, Shale Hills had higher preferential flow frequencies compared to Garner Run. Between these two catchments, the South‐facing mid‐slope sites showed the highest contrasts in preferential flow frequency (33.5% of events at Shale Hills vs. 8.8% at Garner Run) while the ridge‐top sites showed the lowest contrasts (18.7 vs. 13.2%). Additionally, over the unfrozen period, for seven out of eight monitoring sites, drier antecedent conditions tended to be more favourable for preferential flows to occur, with significant (
- NSF-PAR ID:
- 10452265
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 34
- Issue:
- 22
- ISSN:
- 0885-6087
- Format(s):
- Medium: X Size: p. 4295-4306
- Size(s):
- ["p. 4295-4306"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract More above‐ground biomass (kg m−2) grows in the northern Appalachian Mountains (USA) in forests on shale than on sandstone at all landscape positions other than ridgetops. This has been tentatively attributed to physical (rather than chemical) attributes of the substrates, such as elevation, particle size, and water capacity. However, shales have generally similar phosphorus (P) concentrations to sandstones and, in the Valley and Ridge province, they erode more quickly. This led us to hypothesize that faster replenishment of the lithogenic nutrient P in shale soils through erosion + soil production could instead control the differences in biomass. To test this, soils and foliage from 10 sites on shales and sandstones in the northern Appalachians from roughly the same elevation and aspect were analysed. We discovered that, when controlling for location, concentrations of bioavailable P in soils and P in foliage were higher and P resorbed from senescing red oak leaves was lower on slower‐eroding sandstone than on faster‐eroding shale. Lower resorption generally can be attributed to lower P limitation for trees. Further investigation of weathering and erosion on one of the sandstone–shale pairs within a larger, paired watershed study revealed that the differences in P concentrations in biomass and foliage between lithologies likely developed because sandstones act as ‘collectors’ that trap nutrients from residual and exogenous sources, while shales erode quickly and thus promote production of soil from bedrock that releases P to ecosystems. We concluded that the combined effects of differential rates of dust collection and erosion results in roughly equal biomass growing on sandstone and shale ridgetops. This work emphasizes the balance between a landscape's capacity to collect dust versus produce soil in controlling bioavailability of nutrients.
-
null (Ed.)Understanding and predicting catchment responses to a regional disturbance is difficult because catchments are spatially heterogeneous systems that exhibit unique moderating characteristics. Changes in precipitation composition in the Northeastern U.S. is one prominent example, where reduction in wet and dry deposition is hypothesized to have caused increased dissolved organic carbon (DOC) export from many northern hemisphere forested catchments; however, findings from different locations contradict each other. Using shifts in acid deposition as a test case, we illustrate an iterative “process and pattern” approach to investigate the role of catchment characteristics in modulating the steam DOC response. We use a novel dataset that integrates regional and catchment-scale atmospheric deposition data, catchment characteristics and co-located stream Q and stream chemistry data. We use these data to investigate opportunities and limitations of a pattern-to-process approach where we explore regional patterns of reduced acid deposition, catchment characteristics and stream DOC response and specific soil processes at select locations. For pattern investigation, we quantify long-term trends of flow-adjusted DOC concentrations in stream water, along with wet deposition trends in sulfate, for USGS headwater catchments using Seasonal Kendall tests and then compare trend results to catchment attributes. Our investigation of climatic, topographic, and hydrologic catchment attributes vs. directionality of DOC trends suggests soil depth and catchment connectivity as possible modulating factors for DOC concentrations. This informed our process-to-pattern investigation, in which we experimentally simulated increased and decreased acid deposition on soil cores from catchments of contrasting long-term DOC response [Sleepers River Research Watershed (SRRW) for long-term increases in DOC and the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO) for long-term decreases in DOC]. SRRW soils generally released more DOC than SSHCZO soils and losses into recovery solutions were higher. Scanning electron microscope imaging indicates a significant DOC contribution from destabilizing soil aggregates mostly from hydrologically disconnected landscape positions. Results from this work illustrate the value of an iterative process and pattern approach to understand catchment-scale response to regional disturbance and suggest opportunities for further investigations.more » « less
-
more » « less
Abstract Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment production within a landscape. Here, we investigate saprolite and weathered bedrock properties across a ridge‐valley system in the Northern California Coast Ranges, USA, where topography varies with slope aspect such that north‐facing slopes have thicker soils and are more densely vegetated than south‐facing slopes. We use active source seismic refraction surveys to extend observations made in boreholes to the hillslope scale. Seismic velocity models across several ridges capture a high velocity gradient zone (from 1,000 to 2,500 m/s) located ∼4–13 m below ridgetops that coincides with transitions in material strength and chemical depletion observed in boreholes. Comparing this transition depth across multiple north‐ and south‐facing slopes, we find that the thickness of saprolite does not vary with slope aspects. Additionally, seismic survey lines perpendicular and parallel to bedding planes reveal weathering profiles that thicken upslope and taper downslope to channels. Using a rock physics model incorporating seismic velocity, we estimate the total porosity of the saprolite and find that inherited fractures contribute a substantial amount of pore space in the upper 6 m, and the lateral porosity structure varies strongly with hillslope position. The aspect‐independent weathering structure suggests that the contemporary critical zone structure at Rancho Venada is a legacy of past climate and vegetation conditions.
-
more » « less
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 (
D 50 = 0.1–0.2 cm) both before and after the wildfire. In steeper catchments with cliffs, channel sediment is consistently coarse prior to fire (D 50 = 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.
