Premium
Model for Nonreactive Solute Transport in Structured Soils with Continuous Preferential Flow Paths
Author(s) -
Wallach R.,
Steenhuis T. S.
Publication year - 1998
Publication title -
soil science society of america journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.836
H-Index - 168
eISSN - 1435-0661
pISSN - 0361-5995
DOI - 10.2136/sssaj1998.03615995006200040005x
Subject(s) - dimensionless quantity , soil water , outflow , matrix (chemical analysis) , mass transfer , chemistry , flow (mathematics) , reynolds number , mechanics , thermodynamics , soil science , geology , physics , meteorology , chromatography , turbulence
A mathematical model for solute movement in a structured soil with well‐defined and continuous preferential paths was developed. The model divides the soil profile into one mobile and one “stagnant” pore group. For well‐structured soils, the mobile pore group consists of a few well‐connected pores that conduct the nonreactive solute downward very rapidly. Only a narrow matrix layer of stagnant solute along the interface between the two pore groups takes part in solute exchange with the preferential paths. Due to differences in time scales between convection and chemical transfer, the rate of chemical exchange between the preferential path and the active matrix layer for short and moderate times after chemical application is controlled mainly by the preferential flow concentration and, to a lesser extent, by the concentration in the matrix active layer. This is expressed mathematically by multiplying the concentration in the active matrix layer in the mass balance equation by a small dimensionless parameter, ɛ ≪ 1. A regular small perturbations problem is then obtained whose solution is expressed as the sum of the zeroth‐ and first‐order approximations. The simulated breakthrough curve (BTC) is made up of piecewise linear lines and shows a good fit to the measured solute concentration in the outflow of undisturbed columns that have significant preferential flow. Due to the model's simplicity, the transport parameters are estimated directly by fitting the model output to the measured BTCs.