A convective-dispersive-adsorptive flow model for solute transport in soils. I. Model description and some simulations

Abstract

Solute transport predictions are usually based on mathematical (analytical and numerical) solutions of partial differential equations describing convective-dispersive-adsorptive flow and ion-soil interactions. In this paper computer simulations based on a numerical model of solute transport in soil columns have been used to illustrate the significance of flow dynamics and adsorption-desorption parameters in determining the shape and position of solute breakthrough curves. Approximating non-linear isotherms of the form S = KCM, where S is adsorption phase concentration, Cis solution phase concentration, and K and M are constants, by linear isotherms of the form S = KC + å where E is the intercept, to facilitate the use of analytical solutions. has been shown to result in significant errors when the retardation factor fails to account for the intercept å (when å is zero, the retardation factor equals 1 + âK/è, where â is the bulk density and è is water content). Using the dynamic form of Freundlich adsorption, the importance of non-linearity in the equilibrium isotherms is demonstrated. Computer simulations have also been used to demonstrate the extreme importance of intermittent flow, the simultaneous occurrence of adsorption dynamics, ionic fixation and hysteresis during instantaneous adsorption and desorption in determining the transport of solutes in soils.