Open AccessUse of the Richards equation in land surface parameterizations
Author(s) -
Lee Deborah H.,
Abriola Linda M.
Publication year - 1999
Publication title -
journal of geophysical research: atmospheres
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.67
H-Index - 298
eISSN - 2156-2202
pISSN - 0148-0227
DOI - 10.1029/1999jd900951
Subject(s) - richards equation , infiltration (hvac) , surface runoff , soil science , groundwater , environmental science , loam , water content , moisture , surface water , hydrology (agriculture) , drainage , geology , soil water , geotechnical engineering , meteorology , environmental engineering , physics , ecology , biology
Accurately modeling infiltration and soil moisture within land surface parameterization schemes (LSPs) of coupled land surface‐atmosphere models is essential for producing realistic simulations of energy and moisture fluxes and for partitioning precipitation into infiltration, surface runoff, and drainage to groundwater. This report compares simulations of soil moisture, runoff, infiltration, and drainage to groundwater for a bare clay loam using three approaches: a finite difference solution of the vertically integrated Richards equation (an approach commonly used in LSPs), a highly resolved (spatially and temporally) finite element solution of Richards equation, and an analytical kinematic wave solution of Richards equation. Comparisons show that depth‐averaged soil moisture simulated using the vertically integrated Richards equation is only similar to those of the finite element solution for vertical spatial discretizations finer than those employed by most state‐of‐the‐art land surface‐atmosphere transfer schemes. The vertically integrated Richards equation overpredicts soil moisture in the near‐surface soil column and underpredicts drainage to groundwater. The infiltration formulation is found to be critical in partitioning precipitation into runoff, soil moisture, and drainage. Different infiltration formulations and vertical spatial discretizations may partly explain the very different land surface moisture and energy fluxes reported by the LSPs evaluated as part of the Project for Intercomparison of Land Surface Parameterization Schemes (PILPS) Phase 2(b) experiment.