Premium
Convective‐Dispersive Stream Tube Model for Field‐Scale Solute Transport: I. Moment Analysis
Author(s) -
Toride Nobuo,
Leij Feike J.
Publication year - 1996
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/sssaj1996.03615995006000020004x
Subject(s) - skewness , log normal distribution , non equilibrium thermodynamics , dispersion (optics) , field (mathematics) , moment (physics) , porous medium , mechanics , thermodynamics , physics , hydrology (agriculture) , chemistry , statistics , mathematics , geology , porosity , classical mechanics , optics , geotechnical engineering , organic chemistry , pure mathematics
Field‐scale solute transport is typically difficult to model due to the complexity and heterogeneity of flow and transport in natural soils. The stream tube model attempts to stochastically describe transport across the field for relatively short travel distances by viewing the field as a series of independent vertical soil columns. This study investigates the stream tube model with the chemical equilibrium and nonequilibrium convection‐dispersion equation (CDE) for local‐scale transport. A bivarlate (joint) lognormal probability density function was used for three pairs of random transport parameters: (i) the dispersion coefficient, D , and the pore‐water velocity, v ; (ii) the distribution coefficient for linear adsorption, K d , and v ; and (iii) the first‐order rate coefficient for nonequilibrium adsorption, α, and v. Expressions for travel time moments as a rsult of a Dirac input were derived to characterize field‐scale transport according to the stream tube model. The mean breakthrough time for the field‐scale flux‐averaged concentration, ĉ f , was found to be identical to that for the deterministic CDE. Variability in D has generally a minor effect on solute spreading compared with variability in v . Spreading of reactive solutes increased for negatively correlated v and K d , even if the variability in K d was relatively small, while nonequilibrium adsorption further increased spreading. If α was variable, a negative correlation between v and α enhanced the skewness of the breakthrough curve for ĉ f while spreading was independent of the correlation between α and v .