Are all runoff processes the same? Numerical experiments comparing a Darcy-Richards solver to an overland flow-based approach for subsurface storm runoff simulation

A. A. Ameli*, J. R. Craig, J. J. McDonnell

*Corresponding author for this work

Research output: Contribution to journalArticle

19 Citations (Scopus)

Abstract

Hillslope runoff theory is based largely on the differentiation between infiltration excess overland flow, saturation excess overland flow, and subsurface stormflow. Here we explore to what extent a 2-D friction-based overland flow model is useful for predicting hillslope-scale subsurface stormflow, posited here as phenomenologically the same as infiltration excess at depth. We compare our results to a 3-D variably saturated Darcy-Richards subsurface solver for individual rainfall runoff events. We use field data from the well-studied Panola Mountain Experimental hillslope in Georgia USA. Our results show that the two models are largely indistinguishable in terms of their ability to simulate the hillslope hydrograph magnitude and timing for a range of slopes and rainfall depths. Furthermore, we find that the descriptive ability of the overland flow model is comparable to the variably saturated subsurface flow model in terms of its ability to represent the spatial distribution of subsurface stormflow and infiltration across the soil-bedrock interface. More importantly, these results imply that the physics of infiltration excess subsurface stormflow at the soil-bedrock interface is similar to infiltration excess overland flow at the soil surface, in terms of detention storage, loss along the lower boundary, and threshold-like activation at the larger hillslope scale. Given the phenomenological similarity of overland flow and subsurface stormflow and the fact that overland flow model predictions are considerably faster to run (particularly as slope and rainfall depth increase), these findings imply that new forms of hillslope-scale subsurface storm runoff predictions may be possible with the knowledge of bedrock permeability and limited soil information. Finally, this work suggests that the role of soil mantle vis-a-vis subsurface stormflow is mainly as a filter that delays the development of patches of saturation along the bedrock surface. Our model results show that simple realizations of soil based on a few soil depth measurements can possibly be enough to characterize this filtering effect.

Original languageEnglish
Pages (from-to)10008-10028
Number of pages21
JournalWater Resources Research
Volume51
Issue number12
Early online date30 Dec 2015
DOIs
Publication statusPublished - Dec 2015

Keywords

  • subsurface stormflow
  • soil-bedrock interface
  • overland flow
  • Darcy-Richards' equation
  • diffusion wave equation
  • bedrock permeability and topography

Cite this

Are all runoff processes the same? Numerical experiments comparing a Darcy-Richards solver to an overland flow-based approach for subsurface storm runoff simulation. / Ameli, A. A.; Craig, J. R.; McDonnell, J. J.

In: Water Resources Research, Vol. 51, No. 12, 12.2015, p. 10008-10028.

Research output: Contribution to journalArticle

@article{939153e9e0f940a7984f9a99ba1fcf1f,
title = "Are all runoff processes the same? Numerical experiments comparing a Darcy-Richards solver to an overland flow-based approach for subsurface storm runoff simulation",
abstract = "Hillslope runoff theory is based largely on the differentiation between infiltration excess overland flow, saturation excess overland flow, and subsurface stormflow. Here we explore to what extent a 2-D friction-based overland flow model is useful for predicting hillslope-scale subsurface stormflow, posited here as phenomenologically the same as infiltration excess at depth. We compare our results to a 3-D variably saturated Darcy-Richards subsurface solver for individual rainfall runoff events. We use field data from the well-studied Panola Mountain Experimental hillslope in Georgia USA. Our results show that the two models are largely indistinguishable in terms of their ability to simulate the hillslope hydrograph magnitude and timing for a range of slopes and rainfall depths. Furthermore, we find that the descriptive ability of the overland flow model is comparable to the variably saturated subsurface flow model in terms of its ability to represent the spatial distribution of subsurface stormflow and infiltration across the soil-bedrock interface. More importantly, these results imply that the physics of infiltration excess subsurface stormflow at the soil-bedrock interface is similar to infiltration excess overland flow at the soil surface, in terms of detention storage, loss along the lower boundary, and threshold-like activation at the larger hillslope scale. Given the phenomenological similarity of overland flow and subsurface stormflow and the fact that overland flow model predictions are considerably faster to run (particularly as slope and rainfall depth increase), these findings imply that new forms of hillslope-scale subsurface storm runoff predictions may be possible with the knowledge of bedrock permeability and limited soil information. Finally, this work suggests that the role of soil mantle vis-a-vis subsurface stormflow is mainly as a filter that delays the development of patches of saturation along the bedrock surface. Our model results show that simple realizations of soil based on a few soil depth measurements can possibly be enough to characterize this filtering effect.",
keywords = "subsurface stormflow , soil-bedrock interface , overland flow , Darcy-Richards' equation , diffusion wave equation , bedrock permeability and topography",
author = "Ameli, {A. A.} and Craig, {J. R.} and McDonnell, {J. J.}",
note = "Acknowledgments We thank Ali Sarhai, Mahsa Jessri, Willemijn Appels, Ilja van Meerveld, Chris Graham, Jessica Liggett, Tommaso Moramarco, Sven Frei, and Luisa Hopp for useful conversations. We appreciate the feedback of anonymous reviewers and associate editor of Water Resources Research journal. We thank Aquanty, Inc. for the provision of the HydroGeoSphere software used in this experiment. This research was funded by a 2014 NSERC Discovery grant and NSERC Accelerator grant to J. J. M. The data used for all simulations are available upon request from the corresponding author.",
year = "2015",
month = "12",
doi = "10.1002/2015WR017199",
language = "English",
volume = "51",
pages = "10008--10028",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "12",

}

TY - JOUR

T1 - Are all runoff processes the same? Numerical experiments comparing a Darcy-Richards solver to an overland flow-based approach for subsurface storm runoff simulation

AU - Ameli, A. A.

AU - Craig, J. R.

AU - McDonnell, J. J.

N1 - Acknowledgments We thank Ali Sarhai, Mahsa Jessri, Willemijn Appels, Ilja van Meerveld, Chris Graham, Jessica Liggett, Tommaso Moramarco, Sven Frei, and Luisa Hopp for useful conversations. We appreciate the feedback of anonymous reviewers and associate editor of Water Resources Research journal. We thank Aquanty, Inc. for the provision of the HydroGeoSphere software used in this experiment. This research was funded by a 2014 NSERC Discovery grant and NSERC Accelerator grant to J. J. M. The data used for all simulations are available upon request from the corresponding author.

PY - 2015/12

Y1 - 2015/12

N2 - Hillslope runoff theory is based largely on the differentiation between infiltration excess overland flow, saturation excess overland flow, and subsurface stormflow. Here we explore to what extent a 2-D friction-based overland flow model is useful for predicting hillslope-scale subsurface stormflow, posited here as phenomenologically the same as infiltration excess at depth. We compare our results to a 3-D variably saturated Darcy-Richards subsurface solver for individual rainfall runoff events. We use field data from the well-studied Panola Mountain Experimental hillslope in Georgia USA. Our results show that the two models are largely indistinguishable in terms of their ability to simulate the hillslope hydrograph magnitude and timing for a range of slopes and rainfall depths. Furthermore, we find that the descriptive ability of the overland flow model is comparable to the variably saturated subsurface flow model in terms of its ability to represent the spatial distribution of subsurface stormflow and infiltration across the soil-bedrock interface. More importantly, these results imply that the physics of infiltration excess subsurface stormflow at the soil-bedrock interface is similar to infiltration excess overland flow at the soil surface, in terms of detention storage, loss along the lower boundary, and threshold-like activation at the larger hillslope scale. Given the phenomenological similarity of overland flow and subsurface stormflow and the fact that overland flow model predictions are considerably faster to run (particularly as slope and rainfall depth increase), these findings imply that new forms of hillslope-scale subsurface storm runoff predictions may be possible with the knowledge of bedrock permeability and limited soil information. Finally, this work suggests that the role of soil mantle vis-a-vis subsurface stormflow is mainly as a filter that delays the development of patches of saturation along the bedrock surface. Our model results show that simple realizations of soil based on a few soil depth measurements can possibly be enough to characterize this filtering effect.

AB - Hillslope runoff theory is based largely on the differentiation between infiltration excess overland flow, saturation excess overland flow, and subsurface stormflow. Here we explore to what extent a 2-D friction-based overland flow model is useful for predicting hillslope-scale subsurface stormflow, posited here as phenomenologically the same as infiltration excess at depth. We compare our results to a 3-D variably saturated Darcy-Richards subsurface solver for individual rainfall runoff events. We use field data from the well-studied Panola Mountain Experimental hillslope in Georgia USA. Our results show that the two models are largely indistinguishable in terms of their ability to simulate the hillslope hydrograph magnitude and timing for a range of slopes and rainfall depths. Furthermore, we find that the descriptive ability of the overland flow model is comparable to the variably saturated subsurface flow model in terms of its ability to represent the spatial distribution of subsurface stormflow and infiltration across the soil-bedrock interface. More importantly, these results imply that the physics of infiltration excess subsurface stormflow at the soil-bedrock interface is similar to infiltration excess overland flow at the soil surface, in terms of detention storage, loss along the lower boundary, and threshold-like activation at the larger hillslope scale. Given the phenomenological similarity of overland flow and subsurface stormflow and the fact that overland flow model predictions are considerably faster to run (particularly as slope and rainfall depth increase), these findings imply that new forms of hillslope-scale subsurface storm runoff predictions may be possible with the knowledge of bedrock permeability and limited soil information. Finally, this work suggests that the role of soil mantle vis-a-vis subsurface stormflow is mainly as a filter that delays the development of patches of saturation along the bedrock surface. Our model results show that simple realizations of soil based on a few soil depth measurements can possibly be enough to characterize this filtering effect.

KW - subsurface stormflow

KW - soil-bedrock interface

KW - overland flow

KW - Darcy-Richards' equation

KW - diffusion wave equation

KW - bedrock permeability and topography

U2 - 10.1002/2015WR017199

DO - 10.1002/2015WR017199

M3 - Article

VL - 51

SP - 10008

EP - 10028

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 12

ER -