Search

[No Abstract]

Few studies have quantified the differences between celerity and velocity of hillslope water flow and explained the processes that control these differences. Here, we asses these differences by combining a 24-day hillslope sprinkling experiment with a spatially explicit hydrologic model analysis. We focused our work on Watershed 10 at the H. J. Andrews Experimental Forest in western Oregon. Celerities estimated from wetting front arrival times were generally much faster than average vertical velocities of δ2H. In the model analysis, this was consistent with an identifiable effective porosity (fraction of total porosity available for mass transfer) parameter, indicating that subsurface mixing was controlled by an immobile soil fraction, resulting in the attenuation of the δ2H input signal in lateral subsurface flow. In addition to the immobile soil fraction, exfiltrating deep groundwater that mixed with lateral subsurface flow captured at the experimental hillslope trench caused further reduction in the δ2H input signal. Finally, our results suggest that soil depth variability played a significant role in the celerity–velocity responses. Deeper upslope soils damped the δ2H input signal, while a shallow soil near the trench controlled the δ2H peak in lateral subsurface flow response. Simulated exit time and residence time distributions with our hillslope hydrologic model showed that water captured at the trench did not represent the entire modeled hillslope domain; the exit time distribution for lateral subsurface flow captured at the trench showed more early time weighting.

The Hydropedograph Toolbox has been developed to provide a set of standardized tools for analyzing soil moisture time series in an efficient and consistent manner. This toolbox contains various modules that permit the exploration and visualization of key soil hydrological parameters and processes using multi-depth real-time soil moisture monitoring datasets. This includes statistical summary, soil water release curve, preferential flow occurrence, hydraulic redistribution, and the relationship between soil moisture and soil temperature. After describing this toolbox, this paper demonstrates the utility of this toolbox in a case study from the Shale Hills Critical Zone Observatory in USA. The case study illustrates the topographic impacts on soil moisture dynamics along a hillslope transect, and quantifies the frequency of the occurrence of preferential flow, diel fluxes of water, and seasonal storage dynamics. It is expected that such a toolbox, with continued enhancements in the future and wide applications across diverse landscapes, can facilitate the advancement of comparative hydrology and hydropedology.

The spatial patterns of groundwater recharge on hillslopes with a thin soil mantle overlying bedrock are poorly known. Complex interactions between vertical percolation of water through the soil, permeability contrasts between soil and bedrock and lateral redistribution of water result in large spatial variability of water moving into the bedrock. Here, we combine new measurements of saturated hydraulic conductivity of soil mantle and bedrock of the well-studied Panola Mountain experimental hillslope with previously collected (sub)surface topography and soil depth data to quantify the factors affecting the spatial pattern of bedrock groundwater recharge. We use geostatistical characteristics of the measured permeability to generate spatial fields of saturated hydraulic conductivity for the entire hillslope. We perform simulations with a new conceptual model with these random fields and evaluate the resulting spatial distribution of groundwater recharge during individual rainstorms and series of rainfall events. Our simulations show that unsaturated drainage from soil into bedrock is the prevailing recharge mechanism and accounts for 60% of annual groundwater recharge. Therefore, soil depth is a major control on the groundwater recharge pattern through available storage capacity and controlling the size of vertical flux. The other 40% of recharge occurs during storms that feature transient saturation at the soil-bedrock interface. Under these conditions, locations that can sustain increased subsurface saturation because of their topographical characteristics or those with high bedrock permeability will act as hotspots of groundwater recharge when they receive lateral flow. Copyright © 2015 John Wiley & Sons, Ltd.

The difference between celerity and velocity of hillslope water flow is poorly understood. We assessed these differences by combining a 24-day hillslope sprinkling experiment with a spatially explicit hydrologic model analysis. We focused our work at Watershed 10 at the H. J. Andrews Experimental Forest in western Oregon. δ2H label was applied at the start of the sprinkler experiment. Maximum event water (δ2H labeled water) contribution was 26 % of lateral subsurface flow at 20 h. Celerities estimated from wetting front arrival times were generally much faster (on the order of 10–377 mm h−1) than average vertical velocities of δ2H (on the order of 6–17 mm h−1). In the model analysis, this was consistent with an identifiable effective porosity (fraction of total porosity available for mass transfer) parameter, indicating that subsurface mixing was controlled by an immobile soil fraction, resulting in an attenuated δ2H in lateral subsurface flow. Furthermore, exfiltrating bedrock groundwater that mixed with lateral subsurface flow captured at the experimental hillslope trench caused further reduction in the δ2H input signal. Our results suggest that soil depth variability played a significant role in the velocity-celerity responses. Deeper upslope soils damped the δ2H input signal and played an important role in the generation of the δ2H breakthrough curve. A shallow soil (~ 0.30 m depth) near the trench controlled the δ2H peak in lateral subsurface flow response. Simulated exit time and residence time distributions with the hillslope hydrologic model were consistent with our empirical analysis and provided additional insights into hydraulic behavior of the hillslope. In particular, it showed that water captured at the trench was not representative for the hydrological and mass transport behavior of the entire hillslope domain that generated total lateral subsurface flow, because of different exit time distributions for lateral subsurface flow captured at the trench and total lateral subsurface flow.

Despite the widespread acceptance of hydrologic importance, controls on the initiation of preferential flow in natural soil profiles and the frequency of its occurrence at different times of year remain elusive. This study determined the controls and frequency of preferential flow occurrence in the Shale Hills Critical Zone Observatory. Soil moisture profiles and precipitation were monitored at 10 sites along a topographic gradient for >3 yr, encompassing 175 precipitation events. For each event and each site, the flow regime was classified as either preferential flow, sequential flow, or nondetectable flow based on the sequence of soil moisture response at various depths within the same site. Preferential flow here specifically refers to out-of-sequence soil moisture response, with a deeper horizon responding to precipitation earlier than a shallower horizon. Indices describing antecedent precipitation, precipitation characteristics, precipitation timing, and initial soil moisture were examined to determine the characteristics of events that resulted in preferential flow vs. those that resulted in sequential flow. Analyses showed that preferential flow was common throughout the catchment, occurring during 17 to 54% of the 175 events at each of the 10 monitored sites. Preferential flow occurred in at least one site during 90% of the 175 events. While the frequency of preferential flow appeared insensitive to topographic position, the controls on preferential flow initiation varied with landscape position. Analysis of subsets of the time series data showed that while the frequency of preferential flow can be determined from 1 yr of real-time monitoring, the controls on preferential flow require much longer (≥3 yr) monitoring to be reliably identified.

Climate change is expected to increase the intensity of precipitation, but our ability to model the consequences for soil respiration are limited by a lack of data from soils that are saturated and draining. In this study, we used large intact soil columns (28 × 30 cm) to 1) quantify changes in CO 2 flux as soils drain from saturated conditions, and 2) to determine which soil water metrics best predict instantaneous maximum CO 2 flux. The columns were from three agricultural landscape positions that vary in soil properties. We simulated water table fluctuations that were observed at the field site (and predicted to increase in future climate scenarios) by flooding the columns from bottom to surface and then allowing the columns to drain for 96 h while monitoring volumetric soil water content (VWC), water filled pore space (WFPS), water content normalized to field capacity, matric potential, and CO 2 flux. Mean cumulative CO 2 flux was 4649 mg CO 2 ―C m − 2 96 h − 1 . Regardless of landscape position, CO 2 flux rates exhibited a single maximum slightly below saturation, near field capacity. This result suggests that many field studies have not captured soil respiration rates when water availability is optimum for heterotrophic respiration. Across landscape positions, matric potential was the most consistent indicator of instantaneous maximum CO 2 flux, with maximum fluxes occurring within the narrow range of − 0.15 to − 4.89 kPa. In contrast, instantaneous maximum CO 2 flux rates occurred between 95 and 131% of water content normalized to field capacity, 72–97% WFPS, and 29–45% VWC. Thus, our data suggest that instantaneous maximum CO 2 flux rates, a key parameter in ecosystem models, can be predicted across an agricultural landscape with diverse soils if matric potential is used as a water scalar.

Recent studies have called for a new unifying hydrological theory at the hillslope and watershed scale, emphasizing the importance of coupled process understanding of the interactions between hydrology, ecology, pedology, geochemistry and geomorphology. The Biosphere 2 Hillslope Experiment is aimed at tackling this challenge and exploring how climate, soil and vegetation interact and drive the evolution of the hydrologic hillslope behavior. A set of three large-scale hillslopes (18 m by 33 m each) will be built in the climate-controlled experimental biome of the Biosphere 2 facility near Tucson, Arizona, USA. By minimizing the initial physical complexity of these hillslopes, the spontaneous formation of flow pathways, soil spatial heterogeneity, surface morphology and vegetation patterns can be observed over time. This paper documents the hydrologic design process for the Biosphere 2 Hillslope Experiment, which was based on design principles agreed upon among the Biosphere 2 science community. Main design principles were that the hillslopes should promote spatiotemporal variability of hydrological states and fluxes, facilitate transient lateral subsurface flow without inducing overland flow and be capable of supporting vegetation. Hydrologic modeling was used to identify a hillslope configuration (geometry, soil texture, soil depth) that meets the design objectives. The recommended design for the hillslopes consists of a zero-order basin shape with a 10 degree overall slope, a uniform soil depth of 1 m and a loamy sand soil texture. The sensitivity of the hydrologic response of this design to different semi-arid climate scenarios was subsequently tested. Our modeling showed that the timing of rainfall in relation to the timing of radiation input controls the spatiotemporal variability of moisture within the hillslope and the generation of lateral subsurface flow. The Biosphere 2 Hillslope Experiment will provide an excellent opportunity to test hypotheses, observe emergent patterns and advance the understanding of interactions.