Abstract
Issues related to the numerical simulation of moisture migration patterns in the unsaturated zone and in the vicinity of mature trees are explored in this paper. The research is based on the use of Richard's equation for unsaturated moisture flow incorporating a sink term. A numerical solution has been achieved via the finite-element method for spatial discretization along with a finite-difference time-marching scheme. An axisymmetric solution is developed to represent water uptake near an established tree. The approach adopted utilizes radial symmetry and assumes a linear distribution of water extraction rates with both depth and radius. The model has been validated by direct comparison with field measurements recorded (by others) for a mature lime tree located on a boulder clay subsoil. Non-linear hydraulic properties have been obtained from independent published data. A good correlation between field data and simulated results has been achieved. The simulation covers a full annual cycle starting from field capacity in winter, extending through a full spring–summer drying period and subsequent autumn recharge. It is believed that this is the first attempt to simulate the behaviour of an established tree over such a time-scale. This relatively straightforward approach is thought to be suitable for development and application to a range of geo-engineering problems (e.g. slope stability, shrinkage/heave prediction, etc).