Procedure to Obtain Analytical Solutions of One-Dimensional Richards’ Equation for Infiltration in Two-Layered Soils
Publication: Journal of Hydrologic Engineering
Volume 21, Issue 7
Abstract
Richards’ equation (RE) constitutes the most used physically based formula for transient infiltration modeling in unsaturated soils, and it can be solved by using analytical or numerical procedures. Analytical methods require that certain assumptions should be made regarding the closed-form equations that are derived; if the assumptions can be considered reasonable, then the analytical procedures can represent a simple and practical tool (except some cases in which complicated boundary conditions could make the application of the methodology more difficult). Numerical methods do not have the strong limitations of available analytic procedures, and for this reason many numerical schemes are proposed in technical literature to solve RE. Nevertheless, numerical techniques could require high computational costs for application in large areas with respect to analytical procedures. For the frequent case of layered unsaturated soils, analytical solutions were derived by only imposing equal values of capillarity rise for different layers, and this is a clear questionable assumption. A procedure for one dimension is proposed in this work in order to obtain analytical solutions for any soil property: different values of capillarity rise of considered layers, any variation of saturated hydraulic conductivity with depth. The obtained results highlight a good capability of the proposed methodology to reproduce all the considered benchmark cases, constituted by analytical solutions reported in technical literature and application of a software.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
The contents of this paper include activities of a research collaboration funded by the Calabria Region, through the University of Calabria (ROP ESP 2007/2013—IV Axis Human Capital—Operative Objective M2—Action d.5) and the International Centre for Geohazards (ICG) at the Norwegian Geotechnical Institute (NGI). Their support is acknowledged with thanks.
References
Aristizabal, E., Velez, J. I., Martinez, H. E., and Jaboyedoff, M. (2015). “SHIA_Landslide: A distributed conceptual and physically based model to forecast the temporal and spatial occurrence of shallow landslides triggered by rainfall in tropical and mountainous basins.” Landslides, 1–21.
Babu, D. K. (1976). “Infiltration analysis and perturbation methods—3. Vertical infiltration.” Water Resour. Res., 12(5), 1019–1024.
Baca, R. G., Chung, J. N., and Mulla, D. J. (1997). “Mixed transform finite element method for solving the nonlinear equation for flow in variably saturated porous media.” Int. J. Numer. Methods Fluids, 24(5), 441–455.
Barontini, S., Ranzi, R., and Bacchi, B. (2007). “Water dynamics in a gradually nonhomogeneous soil described by the linearized Richards equation.” Water Resour. Res., 43(8), W08411.
Basha, H. A. (2011a). “Infiltration models for semi-infinite soil profiles.” Water Resour. Res., 47(8), W08516.
Basha, H. A. (2011b). “Infiltration models for soil profiles bounded by a water table.” Water Resour. Res., 47(10), W10527.
Baum, R. L., Godt, J. W., and Savage, W. Z. (2010). “Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration.” J. Geophys. Res., 115(F3), 1–26.
Baum, R. L., Savage, W. Z., and Godt, J. W. (2008). “TRIGRS—A FORTRAN program for transient rainfall infiltration and grid-based regional slope stability analysis, version 2.0.”, U.S. Geological Survey, Denver Publishing Service Center, Denver.
Bergamaschi, L., and Putti, M. (1999). “Mixed finite element and Newton-type linearizations for the solution of Richards’ equation.” Int. J. Numer. Methods Eng., 45, 1025–1046.
Broadbridge, P., Knight, J. H., and Rogers, C. (1988). “Constant rate infiltration into a bounded profile: Solutions of a nonlinear model.” Soil Sci. Soc. Am. J., 52(6), 1526–1533.
Broadbridge, P., and Rogers, C. (1990). “Exact solutions for vertical drainage and redistribution in soils.” J. Eng. Math., 24(1), 25–43.
Broadbridge, P., and White, I. (1988). “Constant rate rainfall infiltration: A versatile nonlinear model. 1. Analytic solution.” Water Resour. Res., 24(1), 145–154.
Brooks, R. H., and Corey, A. T. (1964). “Hydraulic properties of porous media.”, Colorado State Univ., Fort Collins, CO.
Brunone, B., Ferrante, M., Romano, N., and Santini, A. (2003). “Numerical simulations of one-dimensional infiltration into layered soils with the Richards equation using different estimates of the interlayer conductivity.” Vadose Zone J., 2(2), 193–200.
Bunsri, T., Sivakumar, M., and Hagare, D. (2008). “Numerical modeling of tracer transport in unsaturated porous media.” J. Appl. Fluid Mech., 1(1), 62–70.
Camporese, M., Ferraris, S., Putti, M., Salandin, P., and Teatini, P. (2006). “Hydrological modeling in swelling/shrinking peat soils.” Water Resour. Res., 42(6), W06420.
Casulli, V., and Zanolli, P. (2010). “A nested Newton type algorithm for finite volume methods solving Richards equation in mixed form.” Siam J. Sci. Comput., 32(4), 2255–2273.
Celia, M. A., Bouloutas, E. T., and Zarba, R. L. (1990). “A general mass-conservative numerical solution of the unsaturated flow equation.” Water Resour. Res., 26(7), 1483–1496.
Chen, H. X., and Zhang, L. M. (2014). “A physically-based distributed cell model for predicting regional rainfall-induced shallow slope failures.” Eng. Geol., 176, 79–92.
Cho, S. E. (2009). “Infiltration analysis to evaluate the surficial stability of two-layered slopes considering rainfall characteristics.” Eng. Geol., 105(1–2), 32–43.
Chounet, L. M., Hilhorst, D., Jouron, C., Kelanemer, Y., and Nicolas, P. (1999). “Simulation of water flow and heat transfer in soils by means of a mixed finite element method.” Adv. Water Resour., 22(5), 445–460.
Clement, T. P., William, R. W., and Molz, F. J. (1994). “A physically based two-dimensional, finite difference algorithm for modeling variably saturated flow.” J. Hydrol., 161(1–4), 71–90.
Cordano, E., and Rigon, R. (2013). “A mass-conservative method for the integration of the two-dimensional groundwater (Boussinesq) equation.” Water Resour. Res., 49(2), 1058–1078.
Dehotin, J., Vazquez, R. F., Braud, I., Debionne, S., and Viallet, P. (2011). “Modeling of hydrological processes using unstructured and irregular grids: 2D groundwater application.” J. Hydrol. Eng., 108–125.
Endrizzi, S., Gruber, S., Dall’Amico, M., and, Rigon, R. (2014). “GEOtop 2.0: Simulating the combined energy and water balance at and below the land surface accounting for soil freezing, snow cover and terrain effects.” Geosci. Model Dev., 7, 2831–2857.
Farthing, M. W., Kees, C. E., Coffey, T. S., Kelley, C. T., and Miller, C. T. (2003). “Efficient steady state solution techniques for variably saturated groundwater flow.” Adv. Water Resour., 26(8), 833–849.
Gardner, W. (1958). “Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table.” Soil Sci., 85(4), 228–232.
Geo-Slope International. (2012). Seep/W for finite element seepage analysis, users manual, Calgary, AB, Canada.
Ghezzehei, T. A., Kneafsey, T. J., and Su, G. W. (2007). “Correspondence of the Gardner and van Genuchten-Mualem relative permeability function parameters.” Water Resour. Res., 43(10), W10417.
Giraldez, J. V., and Sposito, G. (1985). “Infiltration in swelling soil.” Water Resour. Res., 21(1), 33–44.
Kavetski, D., Binning, P., and Sloan, S. W. (2001). “Adaptive backward Euler time stepping with truncation error control for numerical modeling of unsaturated fluid flow.” Int. J. Numer. Methods Eng., 53(6), 1301–1322.
Kees, C. E., and Miller, C. T. (2002). “Higher order time integration methods for two-phase flow.” Adv. Water Resour., 25(2), 159–177.
Lima-Vivancos, V., and Voller, V. R. (2004). “Two numerical methods for modelling variably saturated flow in layered media.” Vadose Zone J., 3(3), 1031–1037.
Lu, Z., Zhang, D., and Robinson, B. A. (2007). “Explicit analytical solutions for one-dimensional steady state flow in layered, heterogeneous unsaturated soils under random boundary conditions.” Water Resour. Res., 43(9), W09413.
Miller, C. T., Abhishek, C., and Farthing, M. W. (2005). “A spatially and temporally adaptive solution of Richards’ equation.” Adv. Water Resour., 29(4), 525–545.
Milly, P. C. D. (1985). “A mass conservation procedure for time-stepping in models of unsaturated flow.” Adv. Water Resour., 8(1), 32–36.
Paniconi, C., and Putti, M. (1994). “A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems.” Water Resour. Res., 30(12), 3357–3374.
Parlange, J.-Y. (1971). “Theory of water movement in soils: 2. One dimensional infiltration.” Soil Sci., 111(3), 170–174.
Parlange, M. B., Prasad, S. N., Parlange, J. -Y., and Romkens, M. J. M. (1992). “Extension of the Heaslet-Alksne technique to arbitrary soil water diffusivities.” Water Resour. Res., 28(10), 2793–2797.
Philip, J. R. (1969). “Theory of infiltration.” Adv. Hydrosci., 5, 215–296.
Philip, J. R. (1991). “Hillslope infiltration: Planar slopes.” Water Resour. Res., 27(1), 109–117.
Prasad, S. N., and Romkens, M. J. M. (1982). “An approximate integral solution of vertical infiltration under changing boundary conditions.” Water Resour. Res., 18(4), 1022–1028.
Press, W. H., Flannery, B. P., Teukolsky, S. A., and Vetterling, W. T. (1986). Numerical recipes, the art of scientific computing, Cambridge University Press, New York.
Richards, L. A. (1931). “Capillary conduction of liquids through porous mediums.” Physics, 1(5), 318–333.
Romano, N., Brunone, B., and Santini, A. (1998). “Numerical analysis of one-dimensional unsaturated flow in layered soils.” Adv. Water Resour., 21(4), 315–324.
Ross, P. J. (2003). “Modeling soil water and solute transport—Fast, simplified numerical solutions.” Agron. J., 95(6), 1352–1361.
Rucker, D. F., Warrick, A. W., and Ferré, T. P. A. (2005). “Parameter equivalence for the Gardner and van Genuchten soil hydraulic conductivity functions for steady vertical flow with inclusions.” Adv. Water Resour., 28(7), 689–699.
Sander, G. C., Cunning, I. F., Hogarth, W. L., and Parlange, J.-Y. (1991). “Exact solution for nonlinear, nonhysteretic redistribution in vertical soil of finite depth.” Water Resour. Res., 27(7), 1529–1536.
Sander, G. C., Parlange, J.-Y., Kuhnel, V., Hogarth, W. L., Lockington, D., and O’Kane, J. P. J. (1988). “Exact nonlinear solution for constant flux infiltration.” J. Hydrol., 97(3–4), 341–346.
Simoni, S., Zanoti, F., Bertoldi, G., and Rigon, R. (2007). “Modelling the probability of occurrence of shallow landslides and channelized debris flow using GEOtop_FS.” Hydrol. Processes, 22(4), 532–545.
Simunek, J., Sejna, M., Saito, H., Sakai, M., and van Genuchten, M. T. (2008a). “The HYDRUS-1D software package for simulating the movement of water, heat, and multiple solutes in variably saturated media, version 4.0, HYDRUS software series 3.” Dept. of Environmental Sciences, Univ. of California Riverside, Riverside, CA.
Simunek, J., van Genuchten, M. T., and Sejna, M. (2008b). “Development and applications of the HYDRUS and STANMOD software packages and related codes.” Vadose Zone J., 7(2), 587–600.
Srivastava, R., and Yeh, T. C. J. (1991). “Analytical solutions for one-dimensional, transient infiltration toward the water table in homogeneous and layered soils.” Water Resour. Res., 27(5), 753–762.
Su, N. (2010). “Theory of infiltration: Infiltration into swelling soils in a material coordinate.” J. Hydrol., 395(1–2), 103–108.
Su, N. (2014). “Mass-time and space-time fractional partial differential equations of water movement in soils: Theoretical framework and application to infiltration.” J. Hydrol., 519, 1792–1803.
Tabuada, M. A., Rego, Z. J. C., Vachaud, G., and Pereira, L. S. (1995). “Two-dimensional irrigation under furrow irrigation: Modeling, its validation and applications.” Agric. Water Manage., 27(2), 105–123.
van Genuchten, M. T. (1980). “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Sci. Soc. Am J., 44(5), 892–898.
Varado, N., Braud, I., Ross, P. J., and Haverkamp, R. (2006). “Assessment of an efficient numerical solution of the 1D Richards equation on bare soil.” J. Hydrol., 323(1–4), 244–257.
Witelski, T. P. (1997). “Perturbation analysis for wetting fronts in Richards’ equation.” Transp. Porous Media, 27(2), 121–134.
Witelski, T. P. (2005). “Motion of wetting fronts moving into partially pre-wet soil.” Adv. Water Resour., 28(10), 1133–1141.
Wu, L. Z., Zhang, L. M., and Huang, R. Q. (2012). “Analytical solution to 1D coupled water infiltration and deformation in two-layer unsaturated soils.” Int. J. Numer. Anal. Methods Geomech., 36(6), 798–816.
Yeh, T.-C. J. (1989). “One-dimensional steady-state infiltration in heterogeneous soils.” Water Resour. Res., 25(10), 2149–2158.
Yuan, F., and Lu, Z. (2005). “Analytical solutions for vertical flow in unsaturated, rooted soils with variable surface fluxes.” Vadose Zone J., 4(4), 1210–1218.
Zhan, T. L. T., and Ng, C. W. W. (2004). “Analytical analysis of rainfall infiltration mechanism in unsaturated soils.” Int. J. Geomech., 273–284.
Information & Authors
Information
Published In
Journal of Hydrologic Engineering
Volume 21 • Issue 7 • July 2016
Copyright
© 2016 American Society of Civil Engineers.
History
Received: Feb 24, 2015
Accepted: Dec 7, 2015
Published online: Mar 21, 2016
Published in print: Jul 1, 2016
Discussion open until: Aug 21, 2016
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.