Anon-linear equilibrium adsorption model wasapplied to predict the movementot solutes through unsatruated soils. All parameters such as water contenet (θ), hydraulic conductivity (K), soilwater diffusivity (D), dispersion coefficient (D.), pore velocity (v) were treated as variables and the latter four also as functions of θ. Water content and solute concentration profiles were predicted with this model assumes the experimental conditions of Elrick et al. (1979) in which KCl was the solute and Brookstonclay loam the soil. Good agreement was found with experimental and predicted water content profileseven though curve fitting was not used in this case. For the solute concentratlon profiles only one parameter,I was varied to fit the experimental data. The best-fit value for I was found to be well wlthin acceptable limits. Although only first order approximations were used the effect of numerical dispersion was not accounte dispersion was not accounted for in the solution since tne value of D., was found not to be significant due to the low pore water velocity in the nusaturated case. Fairly sharp fronts were obtained for the solute concentration profile indicating low dispersion in the unsaturated case.