Abstract:A two-layer model based on the integrated form of Richards' equation (RE) was recently developed to simulate the soil water movement in the roots layer and the vadose zone with a relatively shallow and dynamic water table. The model simulates thickness-averaged volumetric water content and matric suction as opposed to point values and was numerically verified for three soil textures using HYDRUS as a benchmark. However, the strengths and limitations of the two-layer model and its performance in stratified soils and under actual field conditions have not been tested. This study further examined the two-layer model using two numerical verification experiments and, most importantly, tested its performance at site-level under actual, highly variable hydroclimate conditions. Moreover, model parameters were estimated and uncertainty and sources of errors were quantified using a Bayesian framework. First, the two-layer model was evaluated for 231 soil textures under varying soil layer thicknesses with a uniform soil profile. Second, the two-layer model was assessed for stratified conditions where the top and bottom soil layers have contrasting hydraulic conductivities. The model was evaluated by comparing soil moisture and flux estimates to those from the HYDRUS model. Last, a case study of model application using data from a Soil Climate Analysis Network (SCAN) site was presented. Bayesian Monte Carlo (BMC) method was implemented for model calibration and quantifying sources of uncertainty under real hydroclimate and soil conditions. For a homogeneous soil profile, the two-layer model generally had excellent performance in estimating volumetric water content and fluxes, while the model performance slightly declined with increasing layer thickness and coarser textured soils. The model configurations regarding layer thicknesses and soil textures that generate accurate soil moisture and flux estimations were further suggested. With the two layers of contrasting permeability, model-simulated soil moisture contents and fluxes agreed well with those computed by HYDRUS, indicating that the two-layer model accurately handles the water flow dynamics around the layer interface. In the field application, given the highly variable hydroclimate conditions, the two-layer model combined with the BMC method showed good agreement with the observed average soil moisture of the root zone and the vadose zone below (RMSE <0.021 during calibration and <0.023 during validation periods). The contribution of parametric uncertainty to the total model uncertainty was too small compared to other sources. The numerical tests and the site level application showed that the two-layer model can reliably simulate thickness-averaged soil moisture and estimate fluxes in the vadose zone under various soil and hydroclimate conditions. Results also indicated that the BMC method could be a robust framework for vadose zone hydraulic parameters identification and model uncertainty estimation.

Method of Minimum Flux in Saturation Layer for Calculating Stable Water Infiltration Through Layered Soil

Extending the Integrated Richards Equation Method to New Boundary Conditions: Model Description and Validation

Evaluating The Water Holding Capacity Of Multilayer Soil Profiles Using Hydrus-1d And Multi-Criteria Decision Analysis

An Easily Implemented Agro-Hydrological Procedure with Dynamic Root Simulation for Water Transfer in the Crop–soil System: Validation and Application

Semi-analytical solution of pore-water pressure in unsaturated ground and infinite slope considering highly nonlinear soil hydraulic properties

A novel approximate solution to slope rainfall infiltration

Functional Extremum Solution and Parameter Estimation for One-Dimensional Vertical Infiltration using the Brooks-Corey Model

Methodology to simulate unsaturated zone hydrology in Storm Water Management Model (SWMM) for green infrastructure design and evaluation

Efficient Bayesian Inverse Modeling of Water Infiltration in Layered Soils

Soil moisture migration model based on flux-concentration relation under steady rainfall

Numerical Simulation of High-pressure Infiltration Resistance of One-dimensional Soil Column

A two-layer numerical model of soil moisture dynamics: Model development

NUMERICAL SIMULATIONS OF WATER INTERCEPTING AND DRAINAGE ASSOCIATED WITH AIR FILLING FOR LANDSLIDE TREATMENT

Two-Layer numerical model of soil moisture dynamics: Model assessment and Bayesian uncertainty estimation

Infiltration-runoff Model for Layered Soils Considering Air Resistance and Unsteady Rainfall

Water Field Distribution Characteristics under Slope Runoff and Seepage Coupled Effect Based on the Finite Element Method

Rainfall infiltration model considering gas resistance effect in layered soil

Approximate water vertical infiltration extremum path in unsaturated soils

Tracing the soil water flux in variably saturated zones using temperature data

Feasibility Investigation of Improving the Modified Green–Ampt Model for Treatment of Horizontal Infiltration in Soil

A modified model predicting the hydraulic properties of variably saturated soil accounting for film and capillary flow