VZJ
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 8 March 2006
Published in Vadose Zone J 5:184-203 (2006)
DOI: 10.2136/vzj2005.0024
© 2006 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Full Text Free
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Vanderborght, J.
Right arrow Articles by Vereecken, H.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Vanderborght, J.
Right arrow Articles by Vereecken, H.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Vanderborght, J.
Right arrow Articles by Vereecken, H.
Related Collections
Right arrow Solute Transport Models
Right arrow Stochastic Processes
Right arrow Vadose Zone Processes and Chemical Transport

Stochastic Continuum Transport Equations for Field-Scale Solute Transport

Overview of Theoretical and Experimental Results

Jan Vanderborght*, Roy Kasteel and Harry Vereecken

Agrosphere, ICG-IV, Forschungszentrum Jülich GmbH, D-52425 Jülich, Germany

Figure 1
View larger version (17K):

[in a new window]
 
Fig. 1. Normalized equivalent dispersivity due to spatial variability of local flow velocity, ({lambda}eq{lambda}dL)/({sigma}2f{gamma}1), as a function of normalized travel distance, x1/{gamma}1, in media with isotropic structures ({gamma}1 = {gamma}2 = {gamma}3), elongated structures parallel ({gamma}1/{gamma}2,3 = 5,) and perpendicular ({gamma}1/{gamma}2,3 = 0.2) to the direction of the mean flow. Full lines: {gamma}1/{lambda}dL = 500. Dashed line: {gamma}1/{lambda}dL = 50.

 

Figure 2
View larger version (19K):

[in a new window]
 
Fig. 2. Variance of stream tube velocities, {sigma}vs2 (red), and dispersivities that characterize: the spreading of an averaged breakthrough curve, equivalent dispersivity, {lambda}eq (black), the spreading of locally measured breakthrough curve, stream tube dispersivity, {lambda}s (blue), and the pore-scale dispersivity {lambda}d (green), in reference planes at different distances from the injection surface in an heterogeneous aquifer with {sigma}f2 = 1, {gamma}1 = {gamma}2 = 5 m, {gamma}3 = 1 m, {lambda}dL = 0.1 m, {lambda}dT = 0.01 m. Thick lines are first-order predictions; symbols refer to parameters derived from a numerical experiment in a realization of the conductivity field.

 

Figure 3
View larger version (52K):

[in a new window]
 
Fig. 3. Maps of simulated concentrations distributions in a heterogeneous two-dimensional saturated flow field ({sigma}f2 = 1, {gamma}1 = {gamma}2 = 10 cm, {lambda}dL = 0.2 cm, {lambda}dT = 0.02 cm, {theta}s = 0.5) of three different tracers: inert tracer (top panel), a nonlinearly sorbing reactive tracer (kf = 1, nf = 0.67), and a strongly nonlinearly sorbing tracer (kf = 0.6, nf = 0.3). Concentration maps are shown when the center of mass of the plume is at 40 cm below the injection surface (x1 = 60 cm in the plots).

 

Figure 4
View larger version (18K):

[in a new window]
 
Fig. 4. Breakthrough curves of flux averaged concentrations at two depths, 0.5 and 0.8 m, below the injection surface in a heterogeneous two-dimensional saturated flow field ({sigma}f2 = 1, {gamma}1 = {gamma}2 = 10 cm, {lambda}dL = 0.2 cm, {lambda}dT = 0.02 cm, {theta}s = 0.5) of two different tracers: a nonlinearly sorbing reactive tracer (kf = 1, nf = 0.67) (left panel), and a strongly nonlinearly sorbing tracer (kf = 0.6, nf = 0.3) (right panel). Solid lines are numerically simulated breakthrough curves in realizations of the conductivity field, dashed lines are the predictions by the equivalent CDE (Eq. [42]) and symbols are predictions by a STM (Eq. [41], [43], [44], [45]).

 

Figure 5
View larger version (18K):

[in a new window]
 
Fig. 5. Effect of (a) the travel time/distance and (b) the local-scale dispersion (Pe = {gamma}1/{lambda}dL) on the spatial correlation in the direction transverse to the mean flow direction of the local concentration, {rho}CC (open circles), and stream tube velocities, {rho}vsvs (solid circles). First-order predictions in a heterogeneous aquifer with {sigma}f2 = 1, {gamma}1 = {gamma}2 = 5 m, {gamma}3 = 1 m, {lambda}dL = 0.1 m, {lambda}dT = 0.01 m are shown.

 

Figure 6
View larger version (74K):

[in a new window]
 
Fig. 6. Transport in a layered soil profile. Top panel represents the spatial distribution of the hydraulic conductivity in the profile [{sigma}f2 = 1, {lambda}dL = 0.2 cm, {lambda}dT = 0.02 cm, {theta}s = 0.5, {gamma}1 = {gamma}2 = 10 cm (top layer), {gamma}1 = 20 cm, {gamma}2 = 2 cm (bottom layer)], middle panel the simulated vertical component of the pore water velocity, vx1, and the bottom panel the concentration profile when the center of mass of the plume is at 50 cm below the injection surface. In the areas between the two top and two bottom horizontal lines, the spatial correlation in the horizontal direction of the concentrations was calculated (see Fig. 7).

 

Figure 7
View larger version (15K):

[in a new window]
 
Fig. 7. Spatial correlations of (a) concentrations, {rho}CC, and (b) stream tube velocities, {rho}vsvs, in the direction perpendicular to the mean flow direction in a layered soil profile (see Fig. 6). Symbols refer to correlations derived from numerical simulations and lines from first-order predictions. Correlations of concentrations were derived at 30 cm below the injection (top layer; solid circles, dashed lines) and at 70 cm below the injection (bottom layer; open circles, solid lines) when the center of mass of the plume is at the boundary of the two layers. Correlations of stream tube velocities were derived at the 40-cm depth (top layer; solid circles, dashed lines) and at the 80-cm depth (bottom layer; open circles, solid lines).

 

Figure 8
View larger version (12K):

[in a new window]
 
Fig. 8. Simulated breakthrough curves of spatially averaged solute fluxes at two depths in a realization of a two-dimensional heterogeneous soil profile. Averaged fluxes are plotted vs. the cumulative amount of leachate. Solid circles represent simulated breakthrough in a steady-state flow field; open circles are the results for a transient flow field simulated for climatic boundary conditions.

 





HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome
Copyright © 2006 by the Soil Science Society of America.