Numerical simulation of unsaturated flow along preferential pathways: implications for the use of mass balance calculations for isotope storm hydrograph separation

An objective common to many watershed studies is to separate storm hydrographs into two components: water that was present in the watershed prior to a storm event (soil moisture and groundwater) and water which fell on the watershed during the storm. To use this approach, a number of assumptions must be made including that the composition of water in the soil moisture and groundwater reservoirs are constant and known. The objective of this paper is to show that in settings where flow and transport are dominated by preferential pathways for flow, steady state mass balance calculations for quantitative hydrograph separation may be in error. We present field data from a site where flow and transport are dominated by preferential pathways (relict fractures in saprolite of sedimentary rocks) which indicate that the δ18O content of the water in the unsaturated and shallow saturated zones is not constant over the course of a storm event. We use a numerical model to further explore the interactions between the fractures and surrounding matrix. Both the field data and modeling results indicate that the δ18O of the previous storm event(s) has a strong influence on water in the fractures. On the time scale of a storm event, only the water in the matrix immediately surrounding the fracture mixes with water in the fracture, while the bulk of the matrix is isolated from fracture flow. The spatial and temporal heterogeneity of the δ18O in the subsurface and the isolation of the most of the matrix water from flow in fractures make the measurement of a singular δ18O value for subsurface reservoirs problematic and the assumption of a constant value doubtful. Since most near-surface geologic materials have preferential flow paths, we suggest that quantitative hydrograph separation using mass balance techniques is not possible in most situations. Future field and modeling investigations using the approach outlined here could be designed to explore the important temporal and spatial scales of variability in watersheds, and lead to a more quantitative approach to storm hydrograph separation.

[1]  H. Elsenbeer,et al.  Hydrologic pathways and stormflow hydrochemistry at South Creek, northeast Queensland , 1994 .

[2]  Piotr Maloszewski,et al.  On the theory of tracer experiments in fissured rocks with a porous matrix , 1985 .

[3]  Hydrology of a forested hillslope during storm events , 1990 .

[4]  M. Sklash,et al.  The Role Of Groundwater In Storm Runoff , 1979 .

[5]  Jeffrey J. McDonnell,et al.  Modeling Base Flow Soil Water Residence Times From Deuterium Concentrations , 1991 .

[6]  G. Wilson,et al.  Field-scale transport from a buried line source in variably saturated soil , 1993 .

[7]  David E. Radcliffe,et al.  Soil structure development and preferential solute flow , 1999 .

[8]  D. A. Lietzke,et al.  Hydrogeochemical processes controlling subsurface transport from an upper subcatchment of Walker Branch watershed during storm events. 1. Hydrologic transport processes , 1991 .

[9]  P. Mulholland.  Hydrometric and stream chemistry evidence of three storm flowpaths in Walker Branch Watershed , 1993 .

[10]  Dale D. Huff,et al.  Use of specific conductance and contact time relations for separating flow components in storm runoff , 1979 .

[11]  D. Dewalle,et al.  Three-component tracer model for stormflow on a small Appalachian forested catchment , 1988 .

[12]  E. Sudicky,et al.  Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media , 1996 .

[13]  P. Mulholland,et al.  Use of radon-222 and calcium as tracers in a three-end-member mixing model for streamflow generation on the West Fork of Walker Branch Watershed , 1993 .

[14]  R. J. Luxmoore,et al.  Estimating Macroporosity in a Forest Watershed by use of a Tension Infiltrometer1 , 1986 .

[15]  P. Mulholland,et al.  Hydrogeochemical Response of a Forested Watershed to Storms: Effects of Preferential Flow Along Shallow and Deep Pathways , 1990 .

[16]  J. McDonnell,et al.  Hydrograph Separation Using Continuous Open System Isotope Mixing , 1995 .

[17]  W. Hendershot,et al.  Separating streamflow into groundwater, solum and upwelling flow and its implications for hydrochemical modelling , 1993 .

[18]  D. Genereux,et al.  Quantifying uncertainty in tracer‐based hydrograph separations , 1998 .

[19]  R. J. Luxmoore,et al.  Micro-, Meso-, and Macroporosity of Soil , 1981 .

[20]  Glenn V. Wilson,et al.  Physical and chemical controls of preferred path flow through a forested hillslope , 1990 .

[21]  George F. Pinder,et al.  Determination of the ground‐water component of peak discharge from the chemistry of total runoff , 1969 .