Numerical investigation of saturated source area behavior at the small catchment scale

The objective of this research is to explore the relationship between small catchment properties and the temporal growth and decay of saturated source areas (SSA). A simple physics-based hydrologic model, which we call the Sandbox model, is developed for this purpose. A thorough sensitivity analysis is undertaken to evaluate model response to variations in model parameters. Sandbox model output is compared to that from the semi-distributed conceptual model, TOPMODEL, a model with a wide spread acceptance. Plotting the temporal evolution of the extent of saturated source area versus catchment average soil water content during a number of wetting and drying cycles shows a wide variety of trajectories or hysteretic loops. A parametric analysis was performed to quantify the impact of hypothetical catchment properties on the relationship between saturated area extent and basin-average soil water content, revealing hysteretic behavior. It is shown that this hysteresis is the result of changes in groundwater table slope, which is usually less than, and seldom equal to, the land-surface slope in non-saturated areas.

[1]  R. D. Black,et al.  Partial Area Contributions to Storm Runoff in a Small New England Watershed , 1970 .

[2]  R. Allan Freeze,et al.  Role of subsurface flow in generating surface runoff: 2. Upstream source areas , 1972 .

[3]  K. Beven,et al.  Shenandoah Watershed Study: Calibration of a Topography‐Based, Variable Contributing Area Hydrological Model to a Small Forested Catchment , 1985 .

[4]  William T. Sloan,et al.  A physics‐based function for modeling transient groundwater discharge at the watershed scale , 2000 .

[5]  M. Anderson,et al.  A laboratory model to investigate the soil moisture conditions on a draining slope , 1977 .

[6]  R. Healy,et al.  Documentation of computer program VS2D to solve the equations of fluid flow in variably saturated porous media , 1987 .

[7]  R. Healy,et al.  Simulation of solute transport in variably saturated porous media with supplemental information on modifications to the US Geological Survey's computer program VS2D , 1990 .

[8]  E. O'Loughlin.  Prediction of Surface Saturation Zones in Natural Catchments by Topographic Analysis , 1986 .

[9]  V. Singh,et al.  Computer Models of Watershed Hydrology , 1995 .

[10]  E. O'Loughlin,et al.  Saturation regions in catchments and their relations to soil and topographic properties , 1981 .

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

[12]  C. Paniconi,et al.  Hillslope‐storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes: 1. Formulation and characteristic response , 2003 .

[13]  Fred L. Ogden,et al.  GSSHA: Model To Simulate Diverse Stream Flow Producing Processes , 2004 .

[14]  M. Wigmosta,et al.  A distributed hydrology-vegetation model for complex terrain , 1994 .

[15]  J. Hewlett,et al.  Moisture and energy conditions within a sloping soil mass during drainage , 1963 .

[16]  Tammo S. Steenhuis,et al.  SMoRMod – A GIS-integrated Rainfall-runoff Model , 1996 .

[17]  Walter C. Boughton,et al.  Systematic Procedure for Evaluating Partial Areas of Watershed Runoff , 1990 .

[18]  Keith Beven,et al.  TOPMODEL : a critique. , 1997 .

[19]  Daniel Hillel,et al.  Soil and water , 1971 .

[20]  P. Troch,et al.  Analytical solutions to a hillslope-storage kinematic wave equation for subsurface flow , 2002 .

[21]  Dara Entekhabi,et al.  Basin hydrologic response relations to distributed physiographic descriptors and climate , 2001 .

[22]  R. Wooding,et al.  A hydraulic model for the catchment-stream problem: I. Kinematic-wave theory , 1965 .

[23]  F. Ogden,et al.  Saturated area formation on nonconvergent hillslope topography with shallow soils: A numerical investigation , 2000 .

[24]  C. Paniconi,et al.  Hillslope‐storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes: 2. Intercomparison with a three‐dimensional Richards equation model , 2003 .

[25]  I. Iorgulescu,et al.  Validation of TOPMODEL on a small Swiss catchment , 1994 .

[26]  V. T. Chow,et al.  Hydrodynamic Modeling of Two-Dimensional Watershed Flow , 1973 .

[27]  Paul S. P. Cowpertwait,et al.  Further developments of the neyman‐scott clustered point process for modeling rainfall , 1991 .

[28]  Walter J. Rawls,et al.  Green‐ampt Infiltration Parameters from Soils Data , 1983 .

[29]  P. E. O'connell,et al.  An introduction to the European Hydrological System — Systeme Hydrologique Europeen, “SHE”, 1: History and philosophy of a physically-based, distributed modelling system , 1986 .

[30]  Christopher J. Duffy,et al.  A TWO-STATE INTEGRAL-BALANCE MODEL FOR SOIL MOISTURE AND GROUNDWATER DYNAMICS IN COMPLEX TERRAIN , 1996 .

[31]  R. Betson,et al.  What is watershed runoff , 1964 .

[32]  D. Entekhabi,et al.  Groundwater-surface water interaction and the climatic spatial patterns of hillslope hydrological response , 1999 .

[33]  Keith Beven,et al.  The sensitivity of hydrological models to spatial rainfall patterns: an evaluation using observed data , 1994 .

[34]  K. Beven,et al.  Testing a physically-based flood forecasting model (TOPMODEL) for three U.K. catchments , 1984 .

[35]  K. Beven,et al.  The Relationship of Catchment Topography and Soil Hydraulic Characteristics to Lake Alkalinity in the , 1989 .

[36]  Dara Entekhabi,et al.  Hillslope and Climatic Controls on Hydrologic Fluxes , 1995 .

[37]  Keith Beven,et al.  NEW METHOD DEVELOPED FOR STUDYING FLOW ON HILLSLOPES , 1996 .

[38]  Ben Chie Yen,et al.  A laboratory study of surface runoff due to moving rainstorms , 1969 .