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.
Keywords: Bayesian Monte Carlo; HYDRUS; Maximum Likelihood estimation; Numerical Solution; Richards equation (RE); Soil moisture.