Improving the Numerical Solution of Soil Moisture-Based Richards Equation for Land Models with a Deep or Shallow Water Table

Abstract The soil moisture–based Richards equation is widely used in land models for weather and climate studies, but its numerical solution using the mass-conservative scheme in the Community Land Model is found to be deficient when the water table is within the model domain. Furthermore, these deficiencies cannot be reduced by using a smaller grid spacing. The numerical errors are much smaller when the water table is below the model domain. These deficiencies were overlooked in the past, most likely because of the more dominant influence of the free drainage bottom boundary condition used by many land models. They are fixed here by explicitly subtracting the hydrostatic equilibrium soil moisture distribution from the Richards equation. This equilibrium distribution can be derived at each time step from a constant hydraulic (i.e., capillary plus gravitational) potential above the water table, representing a steady-state solution of the Richards equation. Furthermore, because the free drainage condition h...

[1]  Yaping Shao,et al.  On the Choice of Soil Hydraulic Models in Land-Surface Schemes , 1999 .

[2]  J. Holton.  An introduction to dynamic meteorology , 2004 .

[3]  Ying Fan,et al.  Incorporating water table dynamics in climate modeling: 1. Water table observations and equilibrium water table simulations , 2007 .

[4]  Elfatih A. B. Eltahir,et al.  Representation of Water Table Dynamics in a Land Surface Scheme. Part I: Model Development , 2005 .

[5]  Roger A. Pielke,et al.  Coupled Atmosphere–Biophysics–Hydrology Models for Environmental Modeling , 2000 .

[6]  Reed M. Maxwell,et al.  Development of a Coupled Land Surface and Groundwater Model , 2005 .

[7]  R. Koster,et al.  A catchment-based approach to modeling land surface processes in a general circulation model , 2000 .

[8]  Zong-Liang Yang,et al.  The Project for Intercomparison of Land Surface Parameterization Schemes (PILPS): Phases 2 and 3 , 1993 .

[9]  G. Hornberger,et al.  Empirical equations for some soil hydraulic properties , 1978 .

[10]  Zhenghui Xie,et al.  A new parameterization for surface and groundwater interactions and its impact on water budgets with the variable infiltration capacity (VIC) land surface model , 2003 .

[11]  J. D. Tarpley,et al.  The multi‐institution North American Land Data Assimilation System (NLDAS): Utilizing multiple GCIP products and partners in a continental distributed hydrological modeling system , 2004 .

[12]  S. P. Neuman,et al.  SATURATED-UNSATURATED SEEPAGE BY FINITE ELEMENTS , 1973 .

[13]  G. Hornberger,et al.  A Statistical Exploration of the Relationships of Soil Moisture Characteristics to the Physical Properties of Soils , 1984 .

[14]  R. Betts,et al.  Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model , 2000, Nature.

[15]  Zong-Liang Yang,et al.  A simple TOPMODEL-based runoff parameterization (SIMTOP) for use in global climate models , 2005 .

[16]  Dara Entekhabi,et al.  Equivalent steady soil moisture profile and the time compression approximation in water balance modeling , 1994 .

[17]  P. Milly,et al.  A mass-conservative procedure for time-stepping in models of unsaturated flow , 1985 .

[18]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[19]  P. J. Ross,et al.  Modeling Soil Water and Solute Transport—Fast, Simplified Numerical Solutions , 2003 .

[20]  Praveen Kumar,et al.  Topographic Influence on the Seasonal and Interannual Variation of Water and Energy Balance of Basins in North America , 2001 .

[21]  J. Famiglietti,et al.  Multiscale modeling of spatially variable water and energy balance processes , 1994 .

[22]  Keith W. Oleson,et al.  Simulation of Global Land Surface Conditions from 1948 to 2004. Part I: Forcing Data and Evaluations , 2006 .

[23]  Peter E. Thornton,et al.  Technical Description of the Community Land Model (CLM) , 2004 .

[24]  Zhongbo Yu,et al.  On continental-scale hydrologic simulations with a coupled hydrologic model , 2006 .

[25]  Lehua Pan,et al.  A Transformed Pressure Head‐Based Approach to Solve Richards' Equation for Variably Saturated Soils , 1995 .

[26]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .

[27]  D. Randall,et al.  A Revised Land Surface Parameterization (SiB2) for Atmospheric GCMS. Part I: Model Formulation , 1996 .

[28]  Ann Henderson-Sellers,et al.  Biosphere-atmosphere transfer scheme(BATS) version 1e as coupled to the NCAR community climate model , 1993 .

[29]  R. Dickinson,et al.  The Community Land Model and Its Climate Statistics as a Component of the Community Climate System Model , 2006 .

[30]  John D. Albertson,et al.  Temporal dynamics of soil moisture variability: 1. Theoretical basis , 2003 .

[31]  Peter E. Thornton,et al.  An Improved Canopy Integration Scheme for a Land Surface Model with Prognostic Canopy Structure , 2007 .

[32]  Scott A. Isard,et al.  A Soil Moisture Climatology of Illinois , 1994 .

[33]  Randel Haverkamp,et al.  Parameter constraints on closed-form soilwater relationships , 1992 .

[34]  Gordon B. Bonan,et al.  Evaluating Aspects of the Community Land and Atmosphere Models (CLM3 and CAM3) Using a Dynamic Global Vegetation Model , 2006 .

[35]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[36]  E. Davidson,et al.  The role of deep roots in the hydrological and carbon cycles of Amazonian forests and pastures , 1994, Nature.