Anion and cation leaching through large undisturbed soil cores under different flow regimes. 2 Simulation results

In the first part of this study, the movement of anions and cations during sequential leaching of CaSO 4 and KBr solutions under different flow regimes in 2 large undisturbed cores was simulated using a mobile–immobile convection–dispersion model (MIM). By using parameters from the literature or independently measured, and with 75% of the adsorption/exchange sites allocated to the immobile water volume, the model was able to simulate the leaching of bromide and sulfate in both cores, even when flow rates differed by more than 2 orders of magnitude. The leaching of cations (Ca 2+ and K + ) in the slow flow core, where cation exchange occurred, was also simulated reasonably well. But the simulation of cation leaching in the fast flow core, where the surface-applied solutes appeared very rapidly in the effluent, was poor. Sensitivity analysis of the model output showed the most critical parameters were the mobile water volume m, the fraction of the adsorption sites associated with the mobile water f s, and the mobile–immobile transfer coefficient α. Caution is needed when choosing these parameters for modelling purposes. The values used here gave good simulations of cation and anion leaching at flow rates of 3–20 mm/h, but were not satisfactory at very high flow rates (350 mm/h). Keywords: flow rates, ion exchange, mobile– immobile water, preferential flow, modelling. Australian Journal of Soil Research 37(4) 727 - 742 Full text doi:10.1071/SR98003 © CSIRO 1999