The numerical analysis of water and solute movement in scale heterogeneous profiles

A computer based numerical method is presented for the analysis of water and solute movement in unsaturated heterogeneous porous materials. Such a method is necessary since, for those field studies where solute movement is of concern, the soil profiles under consideration are invariably heterogeneous. The numerical analysis is based on a general one-dimensional finite difference soil water flow model which includes a numerical technique combining the concepts of scale heterogeneity with an interpolative soil water hysteresis model. An explicit finite difference solute movement subroutine is incorporated into the unsaturated flow model to describe the transport of nonreactive solutes. A velocity dependent longitudinal dispersion coefficient is used in the solution of the hydrodynamic dispersion equation. The resulting hysteretic scale heterogeneous solute movement model permits the study of solute dynamics during infiltrating and redistribution in realistically complex spatially varying soil profiles. Results are presented for the leaching of both coarse grading to fine and fine grading to coarse sand profiles. Both vertical and horizontal profiles are studied using either a constant flux or a constant concentration input boundary condition. The four cases studied demonstrate the versatility of the numerical method and emphasise the substantial differences in transport behavior that can arise between heterogeneous and homogeneous profiles.Now with BHP Petroleum Pty. Ltd., GPO Box 1911R, Melbourne, Vic. 3001, Australia.