The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.

[1]  M. J. Abedini,et al.  Performance assessment of low-order versus high-order numerical schemes in the numerical simulation of aquifer flow , 2016 .

[2]  James R. Craig,et al.  Are all runoff processes the same? Numerical experiments comparing a Darcy‐Richards solver to an overland flow‐based approach for subsurface storm runoff simulation , 2015 .

[3]  H. Grip.  Sweden's first forest hydrology field study 1905–1926: contemporary relevance of inherited conclusions and data from the Rokliden Hillslope , 2015 .

[4]  J. McDonnell,et al.  Temporal dynamics of catchment transit times from stable isotope data , 2015 .

[5]  L. Hubert‐Moy,et al.  Upscaling Nitrogen Removal Capacity from Local Hotspots to Low Stream Orders’ Drainage Basins , 2015, Ecosystems.

[6]  J. McDonnell,et al.  Debates—The future of hydrological sciences: A (common) path forward? A call to action aimed at understanding velocities, celerities and residence time distributions of the headwater hydrograph , 2014 .

[7]  J. Craig,et al.  Semianalytical series solutions for three‐dimensional groundwater‐surface water interaction , 2014 .

[8]  Peter A. Troch,et al.  Separating physical and meteorological controls of variable transit times in zero‐order catchments , 2013 .

[9]  J. Craig,et al.  Series solutions for saturated–unsaturated flow in multi-layer unconfined aquifers , 2013 .

[10]  K. Beven,et al.  Integrated modeling of flow and residence times at the catchment scale with multiple interacting pathways , 2013 .

[11]  Keith Beven,et al.  Macropores and water flow in soils revisited , 2013 .

[12]  Mark Bakker,et al.  Semi-analytic modeling of transient multi-layer flow with TTim , 2013, Hydrogeology Journal.

[13]  Amvrossios C. Bagtzoglou,et al.  Methods for simulating solute breakthrough curves in pumping groundwater wells , 2012, Comput. Geosci..

[14]  G. Botter.  Catchment mixing processes and travel time distributions , 2012 .

[15]  J. Seibert,et al.  Water storage in a till catchment. I: Distributed modelling and relationship to runoff , 2011 .

[16]  J. Seibert,et al.  Water storage in a till catchment. II: Implications of transmissivity feedback for flow paths and turnover times , 2011 .

[17]  Xiao-Wei Jiang,et al.  A new analytical solution of topography‐driven flow in a drainage basin with depth‐dependent anisotropy of permeability , 2011 .

[18]  A. Wörman,et al.  The use of spectral analysis-based exact solutions to characterize topography-controlled groundwater flow , 2011 .

[19]  Keith Beven,et al.  Catchment travel time distributions and water flow in soils , 2011 .

[20]  V. Zlotnik,et al.  Effects of Multiscale Anisotropy on Basin and Hyporheic Groundwater Flow , 2011, Ground water.

[21]  Enrico Bertuzzo,et al.  Catchment residence and travel time distributions: The master equation , 2011 .

[22]  S. Dunn,et al.  Using lumped conceptual rainfall–runoff models to simulate daily isotope variability with fractionation in a nested mesoscale catchment , 2011 .

[23]  S. Ge,et al.  An analytical study on stagnation points in nested flow systems in basins with depth‐decaying hydraulic conductivity , 2011 .

[24]  M. Cardenas,et al.  Groundwater flow, transport, and residence times through topography‐driven basins with exponentially decreasing permeability and porosity , 2010 .

[25]  Doerthe Tetzlaff,et al.  Gamma distribution models for transit time estimation in catchments: Physical interpretation of parameters and implications for time‐variant transit time assessment , 2010 .

[26]  J. Davies,et al.  Scaling implications of a new multiple interacting pathways hillslope model , 2010 .

[27]  Doerthe Tetzlaff,et al.  Generality of fractal 1/f scaling in catchment tracer time series, and its implications for catchment travel time distributions , 2010 .

[28]  S. Ge,et al.  Simultaneous rejuvenation and aging of groundwater in basins due to depth‐decaying hydraulic conductivity and porosity , 2010 .

[29]  Enrico Bertuzzo,et al.  Transport in the hydrologic response: Travel time distributions, soil moisture dynamics, and the old water paradox , 2010 .

[30]  S. Ge,et al.  Effect of exponential decay in hydraulic conductivity with depth on regional groundwater flow , 2009 .

[31]  Jeffrey J. McDonnell,et al.  A new time‐space accounting scheme to predict stream water residence time and hydrograph source components at the watershed scale , 2009 .

[32]  S. Dunn,et al.  Using long-term data sets to understand transit times in contrasting headwater catchments , 2009 .

[33]  Aldo Fiori,et al.  Travel time distribution in a hillslope: Insight from numerical simulations , 2008 .

[34]  Ty P. A. Ferré,et al.  Steady state lateral water flow through unsaturated soil layers , 2008 .

[35]  J. Craig.  Analytical solutions for 2D topography-driven flow in stratified and syncline aquifers , 2008 .

[36]  R. Maxwell,et al.  Demonstrating fractal scaling of baseflow residence time distributions using a fully‐coupled groundwater and land surface model , 2008 .

[37]  Peter Salamon,et al.  A review and numerical assessment of the random walk particle tracking method. , 2006, Journal of contaminant hydrology.

[38]  Kellie B. Vaché,et al.  A process‐based rejectionist framework for evaluating catchment runoff model structure , 2006 .

[39]  H. Haitjema,et al.  Are Water Tables a Subdued Replica of the Topography? , 2005, Ground water.

[40]  M. Bakker,et al.  Two‐dimensional steady unsaturated flow through embedded elliptical layers , 2004 .

[41]  Kevin Bishop,et al.  Resolving the Double Paradox of rapidly mobilized old water with highly variable responses in runoff chemistry , 2004 .

[42]  A. Warrick,et al.  Steady infiltration from line sources into a layered profile , 2003 .

[43]  Jeffrey J. McDonnell,et al.  How does rainfall become runoff? A combined tracer and runoff transfer function approach , 2003 .

[44]  Kevin Bishop,et al.  Simulating interactions between saturated and unsaturated storage in a conceptual runoff model , 2003 .

[45]  Kevin Bishop,et al.  Groundwater dynamics along a hillslope: A test of the steady state hypothesis , 2003 .

[46]  R. Haggerty,et al.  Transport with multiple-rate exchange in disordered media. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  J. Kirchner,et al.  Catchment-scale advection and dispersion as a mechanism for fractal scaling in stream tracer concentrations , 2001 .

[48]  T. Gisiger.  Scale invariance in biology: coincidence or footprint of a universal mechanism? , 2001, Biological reviews of the Cambridge Philosophical Society.

[49]  J. Kirchner,et al.  Fractal stream chemistry and its implications for contaminant transport in catchments , 2000, Nature.

[50]  L. Nyberg.  WATER-FLOW PATH INTERACTIONS WITH SOIL HYDRAULIC-PROPERTIES IN TILL SOIL AT GARDSJON, SWEDEN , 1995 .

[51]  Rajesh Srivastava,et al.  Analytical solutions for one-dimensional, transient infiltration towards the water table in homogeneous and layered soils , 1991 .

[52]  A. Rodhe.  On the generation of stream runoff in till soils , 1989 .

[53]  K. Beven.  Infiltration into a class of vertically non-uniform soils , 1984 .

[54]  K. Beven,et al.  A physically based, variable contributing area model of basin hydrology , 1979 .

[55]  Malcolm G. Anderson,et al.  The role of topography in controlling throughflow generation , 1978 .

[56]  R. D. Harr,et al.  Water flux in soil and subsoil on a steep forested slope , 1977 .

[57]  D. Weyman,et al.  MEASUREMENTS OF THE DOWNSLOPE FLOW OF WATER IN A SOIL , 1973 .

[58]  J. Tóth.  A Theoretical Analysis of Groundwater Flow in Small Drainage Basins , 1963 .

[59]  J. Philip,et al.  The Theory of Infiltration , 1958 .

[60]  W. R. Gardner.  SOME STEADY‐STATE SOLUTIONS OF THE UNSATURATED MOISTURE FLOW EQUATION WITH APPLICATION TO EVAPORATION FROM A WATER TABLE , 1958 .

[61]  J. Philip,et al.  THE THEORY OF INFILTRATION: 4. SORPTIVITY AND ALGEBRAIC INFILTRATION EQUATIONS , 1957 .

[62]  A. Wörman,et al.  The Impact of Hydraulic Conductivity on Topography Driven Groundwater Flow , 2007 .

[63]  J. Seibert.  On TOPMODEL's ability to simulate groundwater dynamics , 1999 .

[64]  John D. Hewlett,et al.  Non-Point and Diffused Water Sources: A Variable Source Area Problem , 1975 .