The prevailing opinion on the temporal origin of water in a hillslope stormflow hydrograph is that the pre-event water represents a dominant fraction. Such conclusion is usually based on hydrograph separation techniques using stable water isotopes (or other conservative tracers) in conjunction with a mass balance approach. In this study, a two-dimensional dual-continuum model was used to study preferential flow of water and transport of Oxygen-18 (O-18) in a vertical cross-section of a hillslope located in a temperate spruce forest. The effects of hydrodynamic mixing and the spatiotemporal variability of isotopic signatures on estimated pre-event/event water fractions in the hillslope discharge were studied by means of numerical simulation experiments. Pre-event and event water contributions to hillslope stormflow were evaluated using a two-component mass balance approach combined with the 2D flow and transport simulations involving real as well as synthetic O-18 signatures. Long-term simulations of O-18 transport in the hillslope segment were compared with the observed O-18 content in soil water and in the hillslope effluent. The results of the long-term simulations indicated significant mixing of pre-event and event water occurring near the subsurface trench and in the soil above the soil–bedrock interface where the transfer of O-18 from the soil matrix to the preferential pathways takes place. Despite the dominant role of preferential flow in the generation of hillslope stormflow, the pre-event water formed 52–84% of total subsurface stormflow. The mass balance approach failed in partitioning the hillslope discharge into the pre-event/event water components for two thirds of the selected rainfall–runoff episodes due to similar natural isotopic signatures of pre-event and event water. The analysis showed that spatially and temporally variable exchange of O-18 between the soil matrix and preferential pathways exerted a primary control on the estimates of the temporal origin of water in the hillslope runoff. It was demonstrated that the degree of hydrodynamic mixing in the flow domain played an important role in the interpretation of the isotope-based hydrograph separation.
Read full abstract