Predictive-descriptive models for gas and solute diffusion coefficients in variably saturated porous media coupled to pore-size distribution: IV. Solute diffusivity and the liquid phase impedance factor

The solute diffusion coefficient in soil (Dp) and its dependency on volumetric soil-water content (theta) and matric potential (-psi) governs release of dissolved chemicals from low-permeable zones of contaminated soils. Further inspired by the intriguing early history of research findings on mineral nutrition of plants and the key role of ion diffusion played in this, we present a model directly linking the ion (solute) diffusion coefficient in variably saturated soil with both theta and psi. In literature, a linear relationship between the liquid phase impedance factor f(1), (ratio of relative solute diffusion coefficient, Dp/D-0 to theta) and theta is typically observed. Assuming a linear impedance factor equation (LIFE), we couple the LIFE with the Campbell soil-water retention model to relate f(1) and Dp/D-0 to soil-water matric potential. For loamy and clayey soils, the LIFE-Campbell model in agreement with measured data shows an almost linear relationship between f(1) and pF [= log(-psi)] between pF 1.7 and 4.2 (from near field capacity to plant wilting point moisture contents), while the relation between Dp/D-0 and pF is more nonlinear. We also evaluate links between LIFE threshold water content theta(th) (where solute diffusion ceases because of disconnected or thin water films), LIFE slope H (related to pore network tortuosity), and basic soil physical parameters. Linking theta(th) to a newly defined soil physical characteristic (FCvol; the volumetric content of clay plus organic matter) showed to be highly promising. For 21 soils, the value Of theta(th) (m(3) H2O m(-3) soil) is close to FCvol (m(3) fines m(-3) soil), with theta(th) typically > FCvol for compacted or aggregated soils, and theta(th) = 0.8 FCvol for seven nonaggregated and noncompacted soils with between 5 and 40% clay. Furthermore, theta(th) exhibits a nonlinear relationship with the Philips soil-water sorptivity (S). The LIFE slope H is typically close to unity but is also affected by soil aggregation and compaction. The LIFE-Campbell type model is useful to illustrate effects of pore-size distribution and soil type on solute diffusivity and the impedance factor of the soil liquid phase.