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

• Published in: Hydrological processes 2016-07, Vol.30 (14), p.2438-2450
• Main Authors: Ameli, A. A., McDonnell, J. J., Bishop, K.
• Format: Article
• Language: English
• Subjects: Catchment area; Catchments; Computational fluid dynamics; Depth; Distribution; Environmental Management; Exact solutions; exponential decline in saturated hydraulic conductivity with depth; Flow paths; Fluid flow; Gradients; Hydraulic conductivity; Hydraulics; Hydrology; Hypotheses; integrated flow and transport model; Interpolation; Mass balance; Mathematical models; Miljöledning; Model testing; Numerical schemes; Oceanografi, hydrologi, vattenresurser; Oceanography, Hydrology, Water Resources; Residence time; Robustness (mathematics); Saturated-Unsaturated flow; semi analytical model; Storm seepage; Subsoils; Subsurface flow; subsurface flow pathline; subsurface flow pathline; Saturated‐Unsaturated flow; Transit time; transit time distribution; Unsaturated flow; Vertical distribution; Water depth; Water flow
• ISSN: 0885-6087; 1099-1085; 1099-1085
• DOI: 10.1002/hyp.10777

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.