While the Land Surface Models (LSM) community vigorously evaluates the accuracy of new conceptual models against observations (validation), the accuracy of model numerical implementation (verification) is seldom examined. In this work, we developed an open source numerical library for solving terrestrial multi-physics processes. The numerical library, which is coupled to the E3SM Land Model (ELM), presently supports the solution of a variably saturated flow model that includes coupled soil and plant (root, xylem) hydraulics and a thermal advection-diffusion model.
This work developed a standalone, extensible numerical library that can robustly solve terrestrial multi-physics processes at global scales. We provide several examples applied to coupled soil and plant hydraulics and soil thermal dynamics and demonstrate the value of the numerical library for verifying a range of multi-dimensional, multi-physics problems.
Current generation LSMs have several shortcomings, including the lack of robust numerical methods to solve discretized equations, monolithic software design that hinders testing of a process representation in isolation from other components, and absence of a flexible coupling framework to solve tightly couple multi-physics problems. While the LSM community vigorously engages in accessing the accuracy of the conceptual model to represent the real system (validation) through numerous Model Intercomparison Projects, the fidelity of the numerical implementation of the conceptual model (verification) is seldom evaluated. In this work, we develop a numerically robust standalone Multi-Physics Problem (MPP) library for solving terrestrial processes at global scales with support for flexible coupling strategies. The MPP library uses the Method of Manufactured Solutions (MMS) to verify the numerical implementation when the exact solution is unknown. The underlying numerical engine of the MPP library is PETSc, whose DMComposite subclass provides an interface for assembling individual physics components as the part of a global matrix for solving tightly coupled multi-physics problems. Physics formulations supported by the MPP library include the variably saturated Richards equation and a thermal advection-diffusion model. The MMS technique was applied to verify the MPP library for a range of single and multi-physics problems that were solved in one or three dimensions. The numerical library solves variably saturated flow that includes coupled soil and plant (root, xylem) hydraulics and a thermal advection-diffusion model. For all problems, the error between the simulated and manufactured solution reduced as mesh size decreased and expected order of accuracy metric was obtained. The MPP library has been coupled with the E3SM Land Model.