Simulation of water flow and solute transport in free-drainage lysimeters and field soils with heterogeneous structures

Lysimeters are valuable for studying the fate and transport of chemicals in soil. Large-scale field lysimeters are used to assess pesticide behaviour and radionuclide transport, and are assumed to represent natural field conditions better than laboratory columns. Field lysimeters are usually characterized by a free-draining lower boundary. As a result, the hydraulic gradient is disrupted, and leachate cannot be collected until the bottom of the lysimeter becomes saturated. We compared heterogeneously structured, free-drainage lysimeters and field soils with respect to water flow and solute transport. Numerical simulations were carried out in a two-dimensional heterogeneous sandy soil under unsaturated water flow conditions with the CHAIN_2D code. Three different soil structures (isotropic, horizontal, and vertical) were generated, and Miller-Miller similitude was used to scale the hydraulic properties of the soil. The results showed that ponding occurs at the bottom of the lysimeter for the three soil structures and that it occurred faster and was more pronounced with the vertical structure (preferential flow effect). Breakthrough curves of a conservative solute (bromide) showed that solutes are moving faster in the field than in the lysimeters. Fewer differences between lysimeters and field soils were found with the horizontal soil structure than with the isotropic and vertical structures.