This paper presents a pseudo-deterministic catchment runoff model based on the Richards equation model - the governing equation for subsurface flow. The subsurface flow in a catchment is described here by two-dimensional variably saturated flow (unsaturated and saturated). The governing equation is the Richards equation with a slight modification of the time derivative term, as considered e.g. by Neuman. The nonlinear nature of this problem appears in the unsaturated zone only, so it was possible to make use of adaptive domain decomposition algorithm. However delineating of the saturated zone boundary is a nonlinear computationally expensive issue. The simple one-dimensional Boussinesq equation was used here as a rough estimator of the saturated zone boundary. With this estimate the adaptive domain decomposition could always start with an optimal subdomain split, and thus it is now possible to avoid solving huge systems of linear equations in the initial iteration level.
With this measure it is possible to construct an efficient two-dimensional pseudodeterministic catchment runoff model. Finally, the model is tested against real data originating from the Modrava 2 experimental catchment, Czech Republic.
Whipkey R. Z. Subsurface storm-flow from forested slopes, International Association of Scientific Hydrology, Vol. 10, No. 2, 1965, pp. 74-85.
Fletcher P. W. The hydrologic function of forest soil in watershed management, Journal of Forestry, Vol. 50, 1952, pp. 359-362.
Freeze R. A. Three-dimensional, transient, Saturated-unsaturated flow in a groundwater basin, Water Resources Research, Vol. 7, No. 2, 1971, pp. 347-366.
Kao C. , Bouarfa S., Zimmer D. Steady state analysis of unsaturated flow above a shallow water-table aquifer drained by ditches, Journal of Hydrology, Vol. 250, No. 4, 2001, pp. 122-133.
Kosugi K. , Uchida T., Mizuyama T. Numerical calculation of soil pipe flow and its effect on water dynamics in a slope, Hydrological Processes, Vol. 18, No. 4, 2004, pp. 777-789.
Tani M. Analysis of runoff storage relationships to evaluate the runoff-buffering potential of a sloping permeable domain, Journal of Hydrology, Vol. 360, No. 4, 2008, pp. 132-146.
Loritz R. , Hassler S. K., Jackisch C., Allrogen N., Schaik L., Wienhöfer J., Zehe E. Picturing and modelling catchments by representative hillslopes, Hydrol. Earth Syst. Sci. Discuss, (in press)
Bishop J. M. , Callaghan M. V., Cey E. E., Bentley L. R. Measurement and simulation of subsurface tracer migration to tile drains in low permeability, macroporous soil, Water Resour. Res, Vol. 51, No. 6, 2015, pp. 3956-3981.
Wienhöfer J. , Zehe, E. Predicting subsurface storm flow response of a forested hillslope - the role of connected flow paths, Hydrol. Earth Syst. Sci, Vol. 18, No. 1, 2014, pp. 121-138.
Klaus J. , Zehe, E. Modeling rapid flow response of a tile-drained field site using a 2D physically based model, assessment of ‘equifinal’ model setups, Hydrol. Process, Vol. 24, No. 12, 2010, pp. 1595-1609.
Troch P. A. , Paniconi C., van Loon E. E. Hillslope-storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes, 1, Formulation and characteristic response, Water Resources Research, Vol. 39, No. 11, 2003, pp. 1-10.
Paniconi C. , Troch P. A., van Loon E. E., Hilberts A. G. L. Hillslope-storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes: 2, Intercomparison with a three-dimensional Richards equation model, Water Resources Research, Vol. 39, No. 11, 2003, pp. 11-17.
Kuraz M. , Mayer P., Leps M, Trpkosova D. An adaptive time discretization of the classical and the dual porosity model of Richards’ equation, Journal of Computational and Applied Mathematics, Vol. 233, No. 12, 2010, pp. 3167-3177.
Kuraz M. , Mayer P. Algorithms for solving Darcian flow in structured porous media, Acta Polytechnica, Vol. 53, No. 4, 2013, pp. 347-358.
Šolín P. , Kuráž M. Solving the nonstationary Richards equation with adaptive hp-FEM, Advances in Water Resources, Vol. 34, No. 9, 2011, pp. 1062-1081.
Kuraz M. , Mayer P., Havlicek V., Pech P. Domain decomposition adaptivity for the Richards equation, Computing, Vol. 95, Suppl. 1, 1913, pp. 501–519.
Kuraz M. , Mayer P., Pech P. Solving the nonlinear Richards equation model with adaptive domain decomposition, Journal of Computational and Applied Mathematics, Vol. 270, 2014, pp. 2-14.
Kuraz M. , Mayer P., Pech P. Solving the nonlinear and nonstationary Richards equation with multilevel adaptive domain decomposition (dd adaptivity), Applied Mathematics and Computation, Vol. 267, No. C, 2015, pp. 207–222.
Richards L. A. Capillary conduction of liquids through porous mediums, Journal of Applied Physics, Vol. 1, No. 5, 1931, pp. 318-334.
Huyakorn P. S. , Thomas S. D., Thompson B. M. Techniques for making finite elements competitive in modeling flow in variably saturated porous media, Water Resources Research, Vol. 20, No. 8, 1984, pp. 1099-1115.
Neuman S. P. , Witherspoon P. A. Finite element method analyzing steady seepage with a free surface, Water Resources Research, Vol. 6, No. 3, 1970, pp. 889-897.
Tocci M. D. , Kelley C. T., Miller C. T. Accurate and economical solution of the pressure head form of the Richards equation by the method of lines, Advances in Water Resources, Vol. 20, No. 1, 1997, pp. 1-14.
van Genuchten M. T. H. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Journal of Soil Science, Vol. 44, No. 5, 1980, pp. 892-898.
Mualem Y. A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resources Research, No. 12, 3, 1976-513.
Wooding R. A. , Chapman T. G. Groundwater flow over a sloping impermeable layer, Journal of Geophysical Research, Vol. 71, No. 12, 1966, pp. 2895-2902.
Hartani T. , Zimmer D., Lesaffre B. Drainage of sloping lands with variable recharge: analytical formulas and model development, Journal of Irrigation and Drainage Engineering, Vol. 127, No. 1, 2001, pp. 8-15.
Kuraz M. , Mayer P. DRUtES - an open source library for coupled nonlinear convectiondiffusion- reaction problems, http://www.drutes.org, (last visited 1 March 2016).
Kuraz M. , Mayer P., Havlicek V., Pech P., Pavlasek J. Dual permeability variably saturated flow and contaminant transport modeling of a nuclear waste repository with capillary barrier protection, Applied Mathematics and Computation, Vol. 219, No. 13, 2013, pp. 7127-7138.
Kruis J. , Koudelka T., Krejůí T. Salt transport in porous materials and its efficient computer implementation, Pollack Periodica, Vol. 9, No. 2, 2014, pp. 111-122.
Brož J. , Kruis J. Selection strategy for fixing nodes in FETI-DP method, Pollack Periodica, Vol. 5, No. 2, 2010, pp. 57-68.
Oudin L. , Hervieu F., Michel C., Perrin C., Andréassian V., Anctil F., Loumagne C. Which potential evapotranspiration input for a lumped rainfall-runoff model? Part 2, Towards a simple and efficient potential evapotranspiration model for rainfall-runoff modeling, Journal of Hydrology, Vol. 303, No. 1–4, 2005, pp. 290–306.
Kuraz M. , Havlicek V., Jacka L. Identifying the hydraulic parameters of Modrava catchment, Journal of Hydrology and Hydromechanics, (under review).