Improving the numerical solution of soil moisture-based Richards equation for land models with a deep or shallow water table

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 has serious deficiencies, a new bottom boundary condition based on the equilibrium soil moisture distribution at each time step is proposed that also provides an effective and direct coupling between groundwater and surface water.