An integral-balance nonlinear model to simulate changes in soil moisture,groundwater and surface runoff dynamics at the hillslope scale |
| |
| Institution: | 1. Department of Mathematics, The University of Iowa, Iowa City, IA 52242, USA;2. IIHR – Hydroscience & Engineering, The University of Iowa, Iowa City, IA 52242, USA;1. Consiglio Nazionale delle Ricerche, Istituto di Ricerca Sulle Acque, Via Salaria km 29.300, 00015 Monterotondo, RM, Italy;2. Dipartimento DICATAM, Università degli Studi di Brescia, Via Branze 43, 25123 Brescia, Italy;3. Consiglio Nazionale delle Ricerche, Istituto di Ricerca Sulle Acque, UOS Brugherio, Via del Mulino, 19, 20861 Brugherio, MB, Italy;1. Department of Earth Sciences, “Sapienza” University of Rome, P.le A. Moro 5, 00185 Roma, Italy;2. Department of Physics and Earth Sciences, University of Ferrara, Via Saragat 1, 44122 Ferrara, Italy;3. Department of Environmental, Biological and Pharmaceutical Sciences and Technologies, Second University of Naples, Via Vivaldi 43, 81100 Caserta, Italy |
| |
| Abstract: | We present a system of ordinary differential equations (ODEs) capable of reproducing simultaneously the aggregated behavior of changes in water storage in the hillslope surface, the unsaturated and the saturated soil layers and the channel that drains the hillslope. The system of equations can be viewed as a two-state integral-balance model for soil moisture and groundwater dynamics. Development of the model was motivated by the need for landscape representation through hillslopes and channels organized following stream drainage network topology. Such a representation, with the basic discretization unit of a hillslope, allows ODEs-based simulation of the water transport in a basin. This, in turn, admits the use of highly efficient numerical solvers that enable space–time scaling studies. The goal of this paper is to investigate whether a nonlinear ODE system can effectively replicate observations of water storage in the unsaturated and saturated layers of the soil. Our first finding is that a previously proposed ODE hillslope model, based on readily available data, is capable of reproducing streamflow fluctuations but fails to reproduce the interactions between the surface and subsurface components at the hillslope scale. However, the more complex ODE model that we present in this paper achieves this goal. In our model, fluxes in the soil are described using a Taylor expansion of the underlying storage flux relationship. We tested the model using data collected in the Shale Hills watershed, a 7.9-ha forested site in central Pennsylvania, during an artificial drainage experiment in August 1974 where soil moisture in the unsaturated zone, groundwater dynamics and surface runoff were monitored. The ODE model can be used as an alternative to spatially explicit hillslope models, based on systems of partial differential equations, which require more computational power to resolve fluxes at the hillslope scale. Therefore, it is appropriate to be coupled to runoff routing models to investigate the effect of runoff and its uncertainty propagation across scales. However, this improved performance comes at the expense of introducing two additional parameters that have no obvious physical interpretation. We discuss the implications of this for hydrologic studies across scales. |
| |
| Keywords: | Hillslope model Surface vs subsurface interactions Nonlinear ODE system Shale Hills watershed |
| 本文献已被 ScienceDirect 等数据库收录！ |
|
