Abstract

Abstract. This paper presents WHETGEO and its 1D deployment: a new physically based model simulating the water and energy budgets in a soil column. The purpose of this contribution is twofold. First, we discuss the mathematical and numerical issues involved in solving the Richardson–Richards equation, conventionally known as the Richards equation, and the heat equation in heterogeneous soils. In particular, for the Richardson–Richards equation (R2) we take advantage of the nested Newton–Casulli–Zanolli (NCZ) algorithm that ensures the convergence of the numerical solution in any condition. Second, starting from numerical and modelling needs, we present the design of software that is intended to be the first building block of a new customizable land-surface model that is integrated with process-based hydrology. WHETGEO is developed as an open-source code, adopting the object-oriented paradigm and a generic programming approach in order to improve its usability and expandability. WHETGEO is fully integrated into the GEOframe/OMS3 system, allowing the use of the many ancillary tools it provides. Finally, the paper presents the 1D deployment of WHETGEO, WHETGEO-1D, which has been tested against the available analytical solutions presented in the Appendix.

Highlights

  • The Earth’s critical zone (CZ) is defined as the heterogeneous near-surface environment in which complex interactions involving rock, soil, water, air, and living organisms regulate the natural habitat and determine the availability of life-sustaining resources (National Research Council, 2001)

  • The paper presents the 1D deployment of WHETGEO, WHETGEO-1D, which has been tested against the available analytical solutions presented in the Appendix

  • – WHETGEO makes available the NCZ algorithm, which has a priori convergence ensured for any choice of time step and for a great variety of boundary and initial conditions

Read more

Summary

Introduction

The Earth’s critical zone (CZ) is defined as the heterogeneous near-surface environment in which complex interactions involving rock, soil, water, air, and living organisms regulate the natural habitat and determine the availability of life-sustaining resources (National Research Council, 2001). Central to simulating the processes in the CZ is the study of soil moisture dynamics (Clark et al, 2015a; Tubini et al, 2021b). In the following we suggest that studying the CZ requires tools that are not yet readily available to researchers; we propose one of our own. These tools should be flexible enough to allow the quick embedding of advancements in science

Setting up the water budget
The three or four worlds
The necessary coupling with the energy budget
Heat transport
Informatics
Organization and scope
General issues of the R2 equation
The discretization of the R2 equation
Surface boundary condition
Heat transport numerics issues
The discretization of the heat equation
Driving the heat equation with the surface energy budget
Algorithms
Design and deployment of WHETGEO-1D
Design requirements
The software organization
Generic programming at work
Information for users and developers
Workflow for users
Inputs and outputs
Workflow for developers
Conclusions
Homogeneous soil
Layered soil
Infiltration excess
Saturation excess
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call