Exckss salts may be removed from soil by leaching, but ponding water on the soil's surface and allowing infiltration requires large quantities of water. During such leaching water flows preferentially through macropores between aggregates, while the flow within aggregates is much less. Consequently, solute within aggregates is removed much more slowly, thus decreasing overall leaching efficiency. For this reason intermittent ponding can be more efficient because it allows time for solute to diffuse to the surfaces of aggregates during the rest period and subsequently be removed in macropore flow. We explored solute transport in aggregated soils under intermittent leaching in three ways: theoretically, by laboratory experiments on columns of porous ceramic spheres as analogues of aggregates, and by simulation. Solute movement during displacement is described by the mobile-immobile convectiondispersion equation. During the rest period flow ceases, and solute redistributes within the aggregates by diffusion, the key variable being the effective diffusion coefficient, D,, of the solute in the aggregates, and longitudinally by diffusion within macropores (though this was ignored in the simulation). We estimated D, for our porous spheres from observations of solute outflow into finite volumes of stirred distilled water. The theory was validated against experiments on saturated columns for different aggregatesize distributions, flow velocities, and displacement and rest periods, with most parameters estimated independently. Experiments and simulations showed that water savings of 25% were possible under our laboratory conditions, increasing as aggregate size, flow velocity and duration of rest period increased. The potential of intermittent leaching in the field is considered. Nomenclature a = sphere radius d ds q t =time v z = depth of leaching water applied = depth of soil to be leached = volumetric liquid flux density (= Darcy velocity) = average pore-water velocity (= q/@,) = space co-ordinate (positive downwards) 167 @ b B 1 =0.144721n (7). where bl has various values depending on @, Values of bl for various @ values are given in Table 1 of Rao et al. (1980a) C = solute concentration in solution C, = solute concentration in solution at equilibrium C, = solute concentration in solution between the spheres, the mobile-water region Ci, = solute concentration in solution within the spheres, Ci, = average solute concentration in solution (averaged C*im = overall average of Ci, Cinp = concentration of solution entering the column Co D, Do K, L Li Pi S So T V = outflow volume the immobile-water region over the sphere volume) within the spheres -= initial concentration of solution within the spheres = effective diffusion coefficient in porous spheres = ionic diffusion coefficient in free water = hydrodynamic dispersion coefficient for mobile = length of the porous column = imaginary length of a (semi-infinite) column = mass proportion of spheres of radius ai = amount of solute remaining in the column =initial amount of sol...