Effects of the Extended Water Retention Curve on Coupled Heat and Water Transport in the Vadose Zone

Understanding and simulating coupled heat and water transfer appropriately in the shallow subsurface is of vital significance for accurate prediction of soil evaporation that would improve the coupling between land surface and atmosphere. The theory of Philip and de Vries (1957) and its extensions (de Vries, 1958; Milly, 1982), although physically incomplete, are still adopted successfully to describe the coupled heat and water movement in field soils. However, the adsorptive water retention, which was ignored in Philip and de Vries theory and its extensions for characterizing soil hydraulic parameters, was shown to be non-negligible for soil moisture and evaporation flux calculation in dry field soils based on a recent synthetic analysis (Mohanty and Yang, 2013). In this study, we attempt to comprehensively investigate the effects of full range water retention curve on coupled heat and water transport simulation with a focus on soil moisture content, temperature and soil evaporative flux, based on two synthetic (sand and loam) and two field sites (Riverside, California and Audubon, Arizona) analysis. The results of synthetic sand and loam numerical modeling showed that when neglecting the adsorptive water retention, the resulting simulated soil water content would be larger, and the evaporative flux would be lower, respectively, compared to that obtained by the full range water retention curve mode. The simulated temperature did not show significant difference with or without accounting for adsorptive water retention. The evaporation underestimation when neglecting the adsorptive water retention is mainly caused by isothermal hydraulic conductivity underprediction. These synthetic findings were further corroborated by the Audubon, Arizona field site experimental results. The results from Riverside, California field experimental site showed that the soil surface can reach very dry status, although the soil profile below the drying front is not dry, which also to some extent justifies the necessity of employing full range water retention function in such generally not quite dry scenarios.