Premium
Numerical evaluation of solute dispersion and dilution in unsaturated heterogeneous media
Author(s) -
Cirpka Olaf A.,
Kitanidis Peter K.
Publication year - 2002
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2001wr001262
Subject(s) - infiltration (hvac) , monte carlo method , dilution , vadose zone , soil science , mechanics , porous medium , mathematics , geotechnical engineering , environmental science , hydrology (agriculture) , statistics , soil water , geology , physics , meteorology , thermodynamics , porosity
We investigate flow and solute transport in a heterogeneous unsaturated soil under gravity‐driven flow conditions using a direct simulation Monte Carlo method. We account for the cross correlation between the heterogeneous fields of unsaturated flow parameters representing a silty loam. We determine through particle tracking the statistical moments of the travel time distribution for nonpoint injection of a conservative tracer. The one‐particle mean and variance characterize the mean arrival time and spreading of solute integrated over an entire horizontal cross section of the domain, whereas the two‐particle semivariogram characterizes the spreading observed at individual points. From these quantities we derive macrodispersivities and effective dispersivities. The latter reflect solute dilution. Owing to the nonlinearity of unsaturated flow the dispersivity values at a given depth vary with the mean infiltration rate. The two‐particle semivariogram hardly differs from the one‐particle variance of travel time at low infiltration rates, whereas dilution and macrodispersion rates differ significantly at high infiltration rates. We compare the dispersivities determined through particle tracking with approximations from linear stochastic theories based on the covariance of the Eulerian velocity field [ Dagan , 1988; Fiori and Dagan , 2000]. These linear stochastic theories slightly underestimate both dispersivity types. We present parametric models for the two types of dispersivities as function of the infiltration rate and the depth.