Ion transport through unsaturated soils: field experiments and regional simulations

Abstract

This paper describes the movement of anions and cations through soils at the regional scale using block-scale and regional simulations of one-dimensional ion transport through cultivated soils. The simulations were based on field experiments in a region of about 10 km2 in Lower Saxony, Germany. Transport was modelled with the convection–dispersion equation, and the cation exchange was described using the Gapon equation. We evaluated the spatial variation of cation exchange parameters, obtained estimates valid at the block scale, and simulated the one-dimensional transport of anions and cations. The movement of anions and cations was simulated over blocks using effective transport parameters calculated from local transport parameters. The approach led to a good agreement between measured and predicted concentrations of Br–, Na+, K+, Ca2+ and Mg2+ on four different 1 ha blocks. However, the mean concentrations of K+ in the soil solution in the uppermost horizons could not be described satisfactorily by the model. For the regional simulations, transport and exchange parameters were estimated by block kriging. All variograms of the exchange parameters were spatially structured with correlation lengths varying from 100 m to 300 m. Results of the regional simulations imply that Cl– and K+ were transported substantially deeper in the southern part than in the northern part of the area. The transport depth of the ions strongly depended on the pore water velocities. The simulation of solute transport to the water table showed the influence of the depth of water table on the estimated travel times, superimposing the influence of the transport parameters in the region. The results of the regional simulations also emphasize the importance of careful fertilization, especially in regions with shallow water tables such as in the north of the area.