The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

Stormflow generation in a small rainforest catchment in the Luquillo Experimental Forest,Puerto Rico

Various complementary techniques were used to investigate the stormflow generating processes in a small headwater catchment in northeastern Puerto Rico. Over 100 samples were taken of soil matrix water, macropore flow, streamflow and precipitation, mainly during two storms of contrasting magnitude, for the analysis of calcium, magnesium, silicon, potassium, sodium and chloride. These were combined with hydrometric information on streamflow, return flow, precipitation, throughfall and soil moisture to distinguish water following different flow paths. Geo‐electric sounding was used to survey the subsurface structure of the catchment, revealing a weathering front that coincided with the elevation of the stream channel instead of running parallel to surface topography. The hydrometric data were used in combination with soil physical data, a one‐dimensional soil water model (VAMPS ) and a three‐component chemical mass‐balance mixing model to describe the stormflow response of the catchment. It is inferred that most stormflow travelled through macropores in the top 20 cm of the soil profile. During a large event, saturation overland flow also accounted for a considerable portion of the stormflow, although it was not possible to quantify the associated volume fully. Although the mass‐balance mixing model approach gave valuable information about the various flow paths within the catchment, it was not possible to distill the full picture from the model alone; additional hydrometric and soil physical evidence was needed to aid in the interpretation of the model results. Copyright © 2004 John Wiley & Sons, Ltd.     

Role of bedrock groundwater in the rainfall–runoff process in a small headwater catchment underlain by volcanic rock

The role of bedrock groundwater in rainfall–runoff processes is poorly understood. Hydrometric, tracer and subsurface water potential observations were conducted to study the role of bedrock groundwater and subsurface flow in the rainfall–runoff process in a small headwater catchment in Shiranui, Kumamoto prefecture, south‐west Japan. The catchment bedrock consists of a strongly weathered, fractured andesite layer and a relatively fresh continuous layer. Major chemical constituents and stable isotopic ratios of δ¹⁸O and δD were analysed for spring water, rainwater, soil water and bedrock groundwater. Temporal and spatial variation in SiO₂ showed that stream flow under the base flow condition was maintained by bedrock groundwater. Time series of three components of the rainstorm hydrograph (rainwater, soil water and bedrock groundwater) separated by end member mixing analysis showed that each component fluctuated during rainstorm, and their patterns and magnitudes differed between events. During a typical mid‐magnitude storm event, a delayed secondary runoff peak with 1·0 l s⁻¹ was caused by increase in the bedrock groundwater component, whereas during a large rainstorm event the bedrock groundwater component increased to ≈ 2·5 l s⁻¹. This research shows that the contribution of bedrock groundwater and soil water depends strongly on the location of the groundwater table, i.e. whether or not it rises above the soil–bedrock interface. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of bedrock permeability on subsurface stormflow and the water balance of a trenched hillslope at the Panola Mountain Research Watershed,Georgia, USA

The effect of bedrock permeability on subsurface stormflow initiation and the hillslope water balance is poorly understood. Previous hillslope hydrological studies at the Panola Mountain Research Watershed (PMRW), Georgia, USA, have assumed that the bedrock underlying the trenched hillslope is effectively impermeable. This paper presents a series of sprinkling experiments where we test the bedrock impermeability hypothesis at the PMRW. Specifically, we quantify the bedrock permeability effects on hillslope subsurface stormflow generation and the hillslope water balance at the PMRW. Five sprinkling experiments were performed by applying 882–1676 mm of rainfall over a ～5·5 m × 12 m area on the lower hillslope during ～8 days. In addition to water input and output captured at the trench, we measured transpiration in 14 trees on the slope to close the water balance. Of the 193 mm day^?1 applied during the later part of the sprinkling experiments when soil moisture changes were small, <14 mm day^?1 was collected at the trench and <4 mm day^?1 was transpired by the trees, with residual bedrock leakage of >175 mm day^?1 (91%). Bedrock moisture was measured at three locations downslope of the water collection system in the trench. Bedrock moisture responded quickly to precipitation in early spring. Peak tracer breakthrough in response to natural precipitation in the bedrock downslope from the trench was delayed only 2 days relative to peak tracer arrival in subsurface stormflow at the trench. Leakage to bedrock influences subsurface stormflow at the storm time‐scale and also the water balance of the hillslope. This has important implications for the age and geochemistry of the water and thus how one models this hillslope and watershed. Copyright © 2006 John Wiley & Sons, Ltd.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 John Wiley & Sons, Ltd.     

Hydrological responses to rainfall events including the extratropical cyclone Gloria in two contrasting Mediterranean headwaters in Spain; the perennial font del Regàs and the intermittent Fuirosos

Catchment hydrological responses to precipitation inputs, particularly during exceptionally large storms, are complex and variable, and our understanding of the associated runoff generation processes during those events is limited. Hydrological monitoring of climatically and hydrologically distinct catchments can help to improve this understanding by shedding light on the interplay between antecedent soil moisture conditions, hydrological connectivity, and rainfall event characteristics. This knowledge is urgently needed considering that both the frequency and magnitude of extreme precipitation events are increasing worldwide as a consequence of climate change. In autumn 2018, we installed water level sensors to monitor stream water and near-stream groundwater levels at two Mediterranean forest headwater catchments with contrasting hydrological regimes: Font del Regàs (sub-humid climate, perennial flow regime) and Fuirosos (semi-arid climate, intermittent flow regime). Both catchments are located in northeastern Spain, where the extratropical cyclone Gloria hit in January 2020 and left in ca. 65 h outstanding accumulated rainfalls of 424 mm in Font del Regàs and 230 mm in Fuirosos. During rainfall events of low mean intensity, hydrological responses to precipitation inputs at the semi-arid Fuirosos were more delayed and more variable than at the sub-humid Font del Regàs. We explain these divergences by differences in antecedent soil moisture conditions and associated differences in catchment hydrological connectivity between the two catchments, which in this case are likely driven by differences in local climate rather than by differences in local topography. In contrast, during events of moderate and high mean rainfall intensities, including the storm Gloria, precipitation inputs and hydrological responses correlated similarly in the two catchments. We explain this convergence by rapid development of hydrological connectivity independently of antecedent soil moisture conditions. The data set presented here is unique and contributes to our mechanistic understanding on how streams respond to rainfall events and exceptionally large storms in catchments with contrasting flow regimes.     

On the potential of MetOp ASCAT‐derived soil wetness indices as a new aperture for hydrological monitoring and prediction: a field evaluation over Luxembourg

In situ soil moisture data from the Bibeschbach experimental catchment in Luxembourg are used to evaluate relative surface soil moisture observed with the MetOp‐A Advanced Scatterometer (ASCAT). Filtered and bias‐corrected surface soil wetness indices (SWIs) derived from coarse‐resolution (25 km) C‐band scatterometer observations are shown to be highly correlated (r = 0.86) with catchment‐averaged soil moisture measured in the field. The combination of ASCAT and ENVISAT Advanced Synthetic Aperture Radar (ASAR) data sets yields high‐resolution (1 km) relative surface soil moisture that is equally well correlated with in situ measurements. It is concluded that for soil moisture monitoring applications at a catchment scale, the two soil moisture products are equivalent. The best correlation between the SWI derived from ASCAT and ASCAT‐ASAR with in situ soil moisture observations at ca. 5 cm depth is obtained with a characteristic time length parameter T equal to 288 h. These results suggest that satellite‐derived surface soil wetness may serve as proxy for soil storage that enables the monitoring of abrupt switches in river system dynamics to appear when an effective field capacity is exceeded and rapid subsurface stormflow is initiated. In catchments where soil moisture is the main controlling factor of rapid subsurface flow, MetOp ASCAT–derived SWI has the potential to monitor how a river system approaches a critical threshold. Copyright © 2011 John Wiley & Sons, Ltd.     

Runoff generation in zero-order precambrian shield catchments: The stormflow response of a heterogeneous landscape

Hydrometric and isotopic (oxygen-18) observations were used to delineate the runoff processes operating in several headwater catchments on the Precambrian Shield of Canada. The catchments comprise patches of conifer forest situated on thin soils among areas of lichen-covered granitic bedrock. Horton overland flow occurred from the lichen-bedrock areas in all precipitation events that exceeded 4–6 mm. Runoff from the forest stands occurred mainly as subsurface stormflow, but in some instances saturation overland flow was observed. The occurrence of saturation overland flow was controlled by the topography of the bedrock beneath the forest soils. The area contributing runoff and the pathway by which water was conveyed to the catchment outflow switched from the open lichen-bedrock areas producing overland flow on the rising limb of the storm hydrograph to the forest stands contributing subsurface stormflow on the recession limb of the hydrograph. The areal extent and position of the landscape units in the basin were important to the rate and magnitude of stormflow production. Runoff was generated from the catchments only during and immediately after snowmelt and/or rainfall events. The catchments were dry and/or frozen for about 70% of the year.     

Brixenbach research catchment: Quantification of runoff process proportions in a small Alpine catchment depending on soil moisture states and precipitation characteristics

The Brixenbach valley is a small Alpine torrent catchment (9.2 km², 820–1950 m a.s.l., 47.45°, 12.26°) in Tyrol, Austria. Intensive hydrological research in the catchment since more than 12 years, including a hydrogeological survey, pedological and land use mapping, measurements of precipitation, runoff, soil moisture and infiltration as well as the conduction of rainfall simulations, has contributed to understand the hydrological response of the catchment, its subcatchments and specific sites. The paper presents a synthesis of the research in form of runoff process maps for different soil moisture states and precipitation characteristics, derived with the aid of a newly developed Soil-hydrological model. These maps clearly visualize the differing runoff reaction of different subcatchments. The pasture dominated areas produce high surface flow rates during short precipitation events (1 h, 86 mm) with high rainfall intensity, whilst the forested areas often develop shallow subsurface flow. Dry preconditions lead to a slight reduction of surface flow, long rainfall events (24 h, 170 mm) to a dominance of deep subsurface flow and percolation.     

Run‐off formation in a humid,temperate headwater catchment using a combined hydrological,hydrochemical and isotopic approach (Jizera Mountains,Czech Republic)

The formation of baseflow and stormflow was examined in the 1.18 km² part of the headwater catchment Uhlí?ská, Jizera Mountains, Czech Republic, over the period 2007–2011, by means of run‐off data and environmental tracers ¹⁸O and SiO₂. The baseflow, computed using the digital filter approach BFLOW, contributes 67% to total streamflow and has a mean residence time of 12.3 months. It is formed by groundwater discharge from the valley deluviofluvial granitic sediments, in combination with soil water in weathered layers on hillslopes during rainfall and snowmelt periods. The prevailing source of the groundwater is the infiltration of snowmelt water. Analysis of 20 run‐off events and their hysteretic patterns demonstrated that the stormflow water has a residence time of about 4 months and is generated by preferential flow on hillslopes combined by soil matrix drainage. Because of slower flow in the soil matrix, the enrichment of pore water in SiO₂ is more pronounced. The stormflow and snowmelt water flowing via preferential pathways of upslope minerals soils pushes the pre‐event groundwater through the pathways in wetlands to the stream, and the wetland can be therefore considered as groundwater supplied. This mechanism has been found to be typical for the groundwater‐supplied headwater catchments of the Jizera Mountains and can be also assumed in other mountainous headwaters of the granitic massif in Central Europe. The main methodological contribution of this study are the residence time calculations stratified by baseflow and event flow, identifying run‐off components of different travel times to streams and linking them with geochemical run‐off sources. This achievement was possible because of a comprehensive dataset on hydrology, stable isotopes and silica hydrochemistry in all relevant run‐off generation components. This concept indicates that a possible long‐term change in snowmelt may affect the run‐off regime of headwater catchments to climate or land‐use changes. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatio‐temporal variability of piezometric response on two steep alpine hillslopes

Information on the main drivers of subsurface flow generation on hillslopes of alpine headwater catchments is still missing. Therefore, the dominant factors controlling the water table response to precipitation at the hillslope scale in the alpine Bridge Creek Catchment, Northern Italy, were investigated. Two steep hillslopes of similar size, soil properties and vegetation cover but contrasting topography were instrumented with 24 piezometric wells. Sixty‐three (63) rainfall‐runoff events were selected over three years in the snow‐free months to analyse the influence of rainfall depth, antecedent moisture conditions, hillslope topographic characteristics and soil depth on shallow water table dynamics. Piezometric response, expressed as percentage of well activation and water peak magnitude, was strongly correlated with soil moisture status, as described by an index combining antecedent soil moisture and rainfall depth. Hillslope topography was found to be a dominant control only for the convex‐divergent hillslope and during wet conditions. Timing of water table response depended primarily on soil depth and topographic position, with piezometric peak response occurring later and showing a greater temporal variability at the hillslope bottom, characterized by thicker soil. The relationship between mean hillslope water table level and standard deviation for all wells reflected the timing of the water table response at the different locations along the hillslopes. The outcomes of this research contribute to a better understanding of the controls on piezometric response at the hillslope scale in steep terrain and its role on the hydrological functioning of the study catchment and of other sites with similar physiographic characteristics. Copyright © 2014 John Wiley & Sons, Ltd.     

A modified subsurface stormflow model of hillsides in forest catchment

Subsurface flow plays an important role in forest catchment hydrological processes, for which a modified model is established in this paper. Firstly, by taking soil samples from Natural Preserve Forests of Changbai Mountain, two crucial parameters for subsurface flow, the saturated hydraulic conductivity and effective porosity, were measured in the laboratory. Secondly, submodels of the two parameters varying logarithmically with soil depth were established through regressive analysis. Then a modified subsurface stormflow model (Modified model) was founded by substituting the submodels into a storage–discharge model (Sloan's model), established by Sloan in 1983. Finally, to verify the Modified model, five rainfall events on a simulated hillside were carried out. The subsurface flow processes were simulated using the Modified model, Sloan's model and Robinson's model. The comparison of simulated subsurface stormflow processes using the three models respectively with measured ones showed that the Modified model obtained better accuracy for peak flow and total amount of subsurface stormflow than the other two models. Copyright © 2005 John Wiley & Sons, Ltd.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Suspended sediment transport in a small Mediterranean agricultural catchment

The aim of this study is to analyze suspended sediment transport in a Mediterranean agricultural catchment under traditional soil and water conservation practices. Field measurements were conducted in Can Revull, a small ephemeral catchment (1.03 km²) on the island of Mallorca. This study uses continuous turbidity records to analyse suspended sediment transport regimes, construct and interpret multiple regression models of total suspended sediment concentration (SSC) and of SSC related to stormflow discharge, and assess the sediment loads and yields of three hydrological years (2004–2005 to 2006–2007). An annual average SSC of 17.3 mg l^?1, with a maximum of 2270 mg l^?1, was recorded in the middle of the winter period when rainfall intensities are high and headwater slopes are ploughed and thus bare. Strong seasonal contrasts of baseflow dynamics associated with different degrees of dilution provide a large scatter in SSC and in the derived rating curves, reflecting that other factors control the supply of suspended sediment. Multiple regression models identify rainfall intensity as the most significant variable in sediment supply. However, under baseflow conditions, physical and biological processes generate sediment in the channel that is subsequently removed during high flow. In contrast, when baseflow is not present, rainfall intensity is the only process that supplies sediment to the channel, mostly from hillslopes. Considering the study period as average in terms of total annual rainfall and intensities, suspended sediment yields were an order of magnitude lower than those obtained in other Mediterranean catchments, a factor that can be related to the historical use of soil conservation practices. Copyright © 2009 John Wiley & Sons, Ltd.     

High-frequency monitoring of the occurrence of preferential flow on hillslopes and its relationship with rainfall features,soil moisture and landscape

ABSTRACT

High-frequency monitoring was conducted to quantify the frequency and controlling factors of preferential flow (PF) in a monsoon-influenced sub-humid mountainous catchment (6.48 km²) of Northern China. Rainfall was measured using nine bucket raingauges. Soil moisture probes were set up at 12 sites to observe the PF. Overall, 129 rainfall events were identified during the years 2014–2016. The average PF occurrence was 41%, which increased to 71% during heavy rainfall events (>20 mm) revealing a strong influence of the amount and intensity of rainfall. The study also revealed that the PF increased with antecedent soil moisture. Soil moisture was much higher on flat sites compared to sloping sites, providing evidence that the topography has a strong influence on rainfall infiltration and runoff which, subsequently, influence soil moisture variation and the occurrence of PF. Our findings provide valuable insights into the hydrological processes for studies in regions with similar environmental conditions.     

Soil moisture response to rainfall on the Chinese Loess Plateau after a long‐term vegetation rehabilitation

The Chinese Loess Plateau (CLP) is a unique Critical Zone with deep loess deposits, where soil moisture is primarily replenished by seasonal monsoon rainfall. However, the role of vegetation, coupled with complex topography, on rainwater infiltration on the CLP, especially after long‐term revegetation for controlling erosion, is inadequately quantified. Over the growing season of 2016, we monitored soil moisture at the 30‐min interval at 5 depths (10, 20, 40, 60, and 100 cm) in an afforested catchment and a nearby catchment with natural regrowth of grasses. Two monitoring sites were established in each catchment, one in the downhill gully and the other in the uphill slope. We found that vegetation, topography, and rainfall attributes together determined rainwater infiltration and soil moisture replenishment. An accumulated rainfall amount of 9 mm was required to trigger soil moisture response at 10‐cm depth at the 2 grassland sites and the forestland uphill‐slope site whereas 14 mm of rainfall was required for the forestland gully site covered by dense undergrowth and trees. Rainfall events with larger sums and higher peak intensities permitted rainwater infiltration to deeper soil depths. However, no rain recharged soil moisture to 100‐cm depth during the monitoring period. The forestland uphill‐slope site showed the deepest wetting depth (up to 60‐cm depth), fastest wetting‐front velocity (up to 4 cm/hr below 10‐cm depth), and the most significant soil moisture increase (up to 15% cm ³ cm^?3 increase at 10‐cm depth) after rainfall in the growing season. The grassland gully site had the highest soil water storage, whereas soil moisture was depleted the most at the forestland gully site. Findings of this study reveal the transient dynamics of soil moisture after rainfall on the CLP, which signifies the role of revegetation on rainwater infiltration in the loess Critical Zone.     

Characterizing the climatic water balance dynamics and different runoff components in a poorly gauged tropical forested catchment,Nicaragua

ABSTRACT

The water balance dynamics and runoff components of a tropical forested catchment (46?km²) on the southwestern Pacific coast of Nicaragua were studied combining hydrometry, geological characterization and hydrochemical and isotopic tracers (three-component hydrograph separation). The climatic water balance was estimated for 2010/11, 2011/12 and 2012/13 with net values of 811?mm year^-1, 782?mm year^-1 and –447?mm year^-1, respectively. Runoff components were studied at different spatial and temporal scales, demonstrating that different sources and temporal contributions are controlled by dominant landscape elements and antecedent rainfall. In forested sub-catchments, permeable soils, stratigraphy and steep slopes favour subsurface stormflow generation contributing 50% and 53% to total discharge. At catchment scale, landscape elements such as smooth slopes, wide valleys, deeper soils and water table allow groundwater recharge during rainfall events. Groundwater dominates the hydrograph (50% of total discharge) under dry prior conditions. However, low soil infiltration capacity generates a larger surface runoff component (42%) under wet prior conditions which dominates total discharge. Our results show that forested areas are important to reduce surface runoff and thus soil degradation, which is relevant for the design of water management plans.

Editor D. Koutsoyiannis Associate editor D. Gerten     

Testing the ‘two water worlds’ hypothesis under variable preferential flow conditions

Widespread observations of ecohydrological separation are interpreted by suggesting that water flowing through highly conductive soil pores resists mixing with matrix storage over periods of days to months (i.e., two ‘water worlds’ exist). These interpretations imply that heterogeneous flow can produce ecohydrological separation in soils, yet little mechanistic evidence exists to explain this phenomenon. We quantified the separation between mobile water moving through preferential flow paths versus less mobile water remaining in the soil matrix after free-drainage to identify the amount of preferential flow necessary to maintain a two water world's scenario. Soil columns of varying macropore structure were subjected to simulated rainfall of increasing rainfall intensity (26 mm h⁻¹, 60 mm h⁻¹, and 110 mm h⁻¹) whose stable isotope signatures oscillated around known baseline values. Prior to rainfall, soil matrix water δ²H nearly matched the known value used to initially wet the pore space whereas soil δ¹⁸O deviated from this value by up to 3.4‰, suggesting that soils may strongly fractionate ¹⁸O. All treatments had up to 100% mixing between rain and matrix water under the lowest (26 mm h⁻¹) and medium (60 mm h⁻¹) rainfall intensities. The highest rainfall intensity (110 mm h⁻¹), however, reduced mixing of rain and matrix water for all treatments and produced significantly different preferential flow estimates between columns with intact soil structure compared to columns with reduced soil structure. Further, artificially limiting exchange between preferential flow paths and matrix water reduced bypass flow under the most intense rainfall. We show that (1) precipitation offset metrics such as lc-excess and d-excess may yield questionable interpretations when used to identify ecohydrological separation, (2) distinct domain separation may require extreme rainfall intensities and (3) domain exchange is an important component of macropore flow.