Relatively small changes in the shape of the soil water retention curve near saturation can signi®cantly aect the results of numerical simulations of variably saturated ̄ow, including the performance of the numerical scheme itself in terms of stability and rate of convergence. In this paper, we use a modi®ed form of the van Genuchten±Mualem (VGM) soil hydraulic functions to account for a very ...