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 soilbedrock 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.Subsurface stormflow response to precipitation is also a threshold-like process [e.g., Peters et al., 2003;Tromp-van Meerveld and McDonnell, 2006b] that depends upon the microtopography, mesotopography, and macrotopography-not of the surface topography but the topography of the restricting layer at depth Key Points: Richards' and overland flow approaches similarly simulate subsurface stormflow Infiltration excess subsurface stormflow (SSF) can occur along the bedrock surface SSF and overland flow show similar behaviors in terms of the role of their lower boundaries Correspondence to: A. A. Ameli, aaameli@uwaterloo.ca Citation: Ameli, A. A., J. R. Craig, and J. J. McDonnell (2015), Are all runoff processes the same? Numerical experiments comparing a Darcy-Richards solver to an overland flowbased approach for subsurface storm runoff simulation, Water Resour. Water Resources Research PUBLICATIONS (e.g., soil-bedrock interface, soi...

Groundwater Subsidy From Headwaters to Their Parent Water Watershed: A Combined Field‐Modeling Approach

“…This reflects a system type similar to the theory proposed by Berghuijs and Kirchner () where groundwater can be much older than the streamflow when most of the streamflow originates from the highly conductive soil layer above the bedrock groundwater which produces subsurface stormflow (cf. Ameli et al, ). Such an interpretation is supported by the recent work of Gabrielli ().…”

Does Wetland Location Matter When Managing Wetlands for Watershed‐Scale Flood and Drought Resilience?

“…Ameli et al. ). In the groundwater–surface water interaction model, the water level at each intact wetland was considered to be equal to the median of the long‐term water level in the wetland simulated using the surface flow routing model during the eight‐month study period; this water level was then used as the constant head boundary condition at the intact wetland.…”

“…In addition, given the small rainfall to evapotranspiration ratio in the Nose Creek watershed, typical of a dry continental climate, there is a small likelihood of the occurrence of lateral water movement in the shallow unsaturated zone through subsurface stormflow, which is typically observed in humid forested landscapes with permeable soils and a relatively large rainfall to evapotranspiration ratio (c.f. Ameli et al 2015). In the groundwater-surface water interaction model, the water level at each intact wetland was considered to be equal to the JAWRA JOURNAL OF THE AMERICAN WATER RESOURCES ASSOCIATION median of the long-term water level in the wetland simulated using the surface flow routing model during the eight-month study period; this water level was then used as the constant head boundary condition at the intact wetland.…”

The soil moisture velocity equation

“…Downer and Ogden [] found that “effective” parameters could be used to overcome this effect, but these changes in the parameter values resulted in too much infiltration as the soil surface approached saturation. Two and three‐dimensional Richards' equation simulations of even small catchments at appropriate spatial resolutions require very long run times [ Ameli et al ., ], when the solution domain is required to meet the REV assumption.…”