| Abstract: | Solute leaching in unsaturated soil is influenced by the variability in hydraulic functions (water retention and conductivity) that govern the flow process. Variability in measured soil hydraulic functions of a coarse-, medium- and fine-textured soil group was quantified with the scaling theory of similar media. Solute leaching in these soils was calculated with Monte Carlo simulation assuming, successively, hydraulic conductivity, K, volumetric water content, 0, and pressure head, h, to be constant. In addition to variability in hydraulic functions, variability in the solute retardation factor was also taken into account. To examine this effect five solutes were considered: a conservative solute (chloride), a non-retarded solute subject to decay (nitrate), a retarded solute that does not decay (cadmium) and two organic solutes which are retarded but have different sorption and decay parameters (the pesticide atrazine and a chlorinated hydrocarbon). The numerical results obtained with Monte Carlo simulation were in a number of instances verified with analytical solutions. The three soil groups distinguished showed considerable differences in vulnerability for leaching of the five solutes, emphasizing the importance of the effect of variability in soil hydraulic functions when studying solute leaching. Numerical and analytical results showed good agreement. Therefore, in relatively simple situations analytical solutions are attractive. However, in complicated situations, analytical solutions are cumbersome and numerical solutions are the only realistic alternative. |
