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.
Despite a multitude of small catchment studies, we lack a deep understanding of how variations in critical zone architecture lead to variations in hydrologic states and fluxes. This study characterizes hydrologic dynamics of 15 catchments of the U.S. Critical Zone Observatory (CZO) network where we hypothesized that our understanding of subsurface structure would illuminate patterns of hydrologic partitioning. The CZOs collect data sets that characterize the physical, chemical, and biological architecture of the subsurface, while also monitoring hydrologic fluxes such as streamflow, precipitation, and evapotranspiration. For the first time, we collate time series of hydrologic variables across the CZO network and begin the process of examining hydrologic signatures across sites. We find that catchments with low baseflow indices and high runoff sensitivity to storage receive most of their precipitation as rain and contain clay‐rich regolith profiles, prominent argillic horizons, and/or anthropogenic modifications. In contrast, sites with high baseflow indices and low runoff sensitivity to storage receive the majority of precipitation as snow and have more permeable regolith profiles. The seasonal variability of water balance components is a key control on the dynamic range of hydraulically connected water in the critical zone. These findings lead us to posit that water balance partitioning and streamflow hydraulics are linked through the coevolution of critical zone architecture but that much work remains to parse these controls out quantitatively.
more » « less- NSF-PAR ID:
- 10388227
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 57
- Issue:
- 3
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Precipitation is the primary driver of hydrological models, and its spatial and temporal variability have a great impact on water partitioning. However, in data‐sparse regions, uncertainty in precipitation estimates is high and the sensitivity of water partitioning to this uncertainty is unknown. This is a particular challenge in drylands (semi‐arid and arid regions) where the water balance is highly sensitive to rainfall, yet there is commonly a lack of in situ rain gauge data. To understand the impact of precipitation uncertainty on the water balance in drylands, here we have performed simulations with a process‐based hydrological model developed to characterize the water balance in arid and semi‐arid regions (DRYP: DRYland water Partitioning model). We performed a series of numerical analyses in the Upper Ewaso Ng'iro basin, Kenya driven by three gridded precipitation datasets with different spatio‐temporal resolutions (IMERG, MSWEP, and ERA5), evaluating simulations against streamflow observations and remotely sensed data products of soil moisture, actual evapotranspiration, and total water storage. We found that despite the great differences in the spatial distribution of rainfall across a climatic gradient within the basin, DRYP shows good performance for representing streamflow (KGE >0.6), soil moisture, actual evapotranspiration, and total water storage (
r >0.5). However, the choice of precipitation datasets greatly influences surface (infiltration, runoff, and transmission losses) and subsurface fluxes (groundwater recharge and discharge) across different climatic zones of the Ewaso Ng'iro basin. Within humid areas, evapotranspiration does not show sensitivity to the choice of precipitation dataset, however, in dry lowland areas it becomes more sensitive to precipitation rates as water‐limited conditions develop. The analysis shows that the highest rates of precipitation produce high rates of diffuse recharge in Ewaso uplands and also propagate into runoff, transmission losses and, ultimately focused recharge, with the latter acting as the main mechanism of groundwater recharge in low dry areas. The results from this modelling exercise suggest that care must be taken in selecting forcing precipitation data to drive hydrological modelling efforts, especially in basins that span a climatic gradient. These results also suggest that more effort is required to reduce uncertainty between different precipitation datasets, which will in turn result in more consistent quantification of the water balance. -
more » « less
Abstract Silicon stable isotope ratios (
30Si) of over 150 stream water samples were measured during seven storm events in six small critical zone observatory (CZO) catchments spanning a wide range in climate (sub‐humid to wet, tropical) and lithology (granite, volcanic, and mixed sedimentary). Here we report a cross‐site analysis of this dataset to gain insight into stream 30Si variability across low‐order catchments and to identify potential climate (i.e., runoff), hydrologic, lithologic, and biogeochemical controls on observed stream Si chemical and isotopic signatures. Event‐based 30Si exhibit variability both within and across sites (−0.22‰ to +2.27‰) on the scale of what is observed globally in both small catchments and large rivers. Notably, each site shows distinct 30Si signatures that are preserved even after normalization for bedrock composition. Successful characterization of observed cross‐site behavior requires the merging of two distinct frameworks in a novel combined model describing both non‐uniform fluid transit time distributions and multiple fractionating pathways in application to low‐order catchments. The combined model reveals that site‐specific architecture (i.e., biogeochemical reaction pathways and hydrologic routing) regulates stream silicon export signatures even when subject to extreme precipitation events. -
more » « less
Abstract Trees shape the critical zone and modulate terrestrial water storage yet observed streamflow responses to forest cover change vary. Differences in catchment area, soil water storage, management practices, tree species, and climate are among the many explanations proposed for heterogeneous hydrologic responses. We addressed evidence for the hypothesis that mean annual temperature (MAT) and the phase shift between precipitation and enhanced vegetation index (EVI) peaks,
θ , explain a significant amount of the variation in hydrologic response to forest cover loss. We selected 50 catchments with daily streamflow records spanning eight nations and seven climate regions. Categorical clustering of catchments was performed with MAT, θ, minimum EVI, catchment area, and percentage forest loss. Similar storm event runoff ratio responses to deforestation were best clustered by MAT andθ . High MAT tropical monsoonal catchments (Brazil, Myanmar, and Liberia) exhibited minimal evidence of increasing runoff ratios (increases observed in 9% of catchments). Low MAT subarctic, cold semi-arid, and humid continental catchments (US, Canada, and Estonia) showed consistent runoff increases around the time of snowmelt (94%). The deforestation runoff responses of temperate and subtropical catchments with Mediterranean, humid, and oceanic climates depended strongly onθ . We observe increased runoff following forest loss in a majority of catchments (90%) where precipitation peaks followed peak growing season (max EVI) (US). In contrast, where precipitation peaks preceded the growing season (South Africa and Australia) there was less evidence of increased runoff (25% of catchments). This research supports the strategic implementation of native forest conservation or restoration for simultaneously mitigating the effects of global climate change and regional or local surface runoff. -
more » « less
Abstract High‐elevation mountain catchments are often subject to large climatic and topographic gradients. Therefore, high‐density hydrogeochemical observations are needed to understand water sources to streamflow and the temporal and spatial behaviour of flow paths. These sources and flow paths vary seasonally, which dictates short‐term storage and the flux of water in the critical zone (CZ) and affect long‐term CZ evolution. This study utilizes multiyear observations of chemical compositions and water residence times from the Santa Catalina Mountains Critical Zone Observatory, Tucson, Arizona to develop and evaluate competing conceptual models of seasonal streamflow generation. These models were tested using endmember mixing analysis, baseflow recession analysis, and tritium model “ages” of various catchment water sources. A conceptual model involving four endmembers (precipitation, soil water, shallow, and deep groundwater) provided the best match to observations. On average, precipitation contributes 39–69% (55 ± 16%), soil water contributes 25–56% (41 ± 16%), shallow groundwater contributes 1–5% (3 ± 2%), and deep groundwater contributes ~0–3% (1 ± 1%) towards annual streamflow. The mixing space comprised two principal planes formed by (a) precipitation‐soil water‐deep groundwater (dry and summer monsoon season samples) and (b) precipitation‐soil water‐shallow groundwater (winter season samples). Groundwater contribution was most important during the wet winter season. During periods of high dynamic groundwater storage and increased hydrologic connectivity (i.e., spring snowmelt), stream water was more geochemically heterogeneous, that is, geochemical heterogeneity of stream water is storage‐dependent. Endmember mixing analysis and3H model age results indicate that only 1.4 ± 0.3% of the long‐term annual precipitation becomes deep CZ groundwater flux that influences long‐term deep CZ development through both intercatchment and intracatchment deep groundwater flows.
