HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (3) Citing Articles via Google Scholar Google Scholar Articles by Russo, D. Articles by Laufer, A. Search for Related Content PubMed Articles by Russo, D. Articles by Laufer, A. GeoRef GeoRef Citation Agricola Articles by Russo, D. Articles by Laufer, A. Related Collections Soil Salinity Spatial Variability

ORIGINAL RESEARCH

Numerical Analysis of Transport of Interacting Solutes in a Three-Dimensional Unsaturated Heterogeneous Soil

^a Dep. of Environmental Physics and Irrigation, Institute of Soils, Water and Environmental Sciences, Agricultural Research Organization, The Volcani Center, Bet Dagan 50250, Israel
^b Currently, AMEC Earth & Environmental Ltd., 160 Traders Blvd. East, Suite 110, Mississauga, ON, Canada L4Z 3K7
^* Corresponding author (vwrosd{at}volcani.agri.gov.il)

Received 4 December 2003.

The present study focuses on field-scale flow and transport in three-dimensional, heterogeneous, variably saturated formations, for the case in which the flow is coupled to the transport through the dependence of the hydraulic conductivity and water retentivity on solute concentrations. Numerical simulations of flow and transport of both tracer and mixed Na–Ca solutes were employed to analyze long-term effects of the interactions between the soil solution and the soil matrix on water and solute movement, under transient, nonmonotonic flows. The simulated flows were derived from actual irrigation and weather records, along with water uptake by plant roots. Results of this study suggest that enhanced soil solution–soil matrix interactions, induced by soil alkalinity and dilution of the soil solution, may reduce both the mean and the spatial variability of the hydraulic conductivity, and, concurrently, of the velocity. Consequently, enhanced soil solution–soil matrix interactions may slow down the tracer solute movement, decrease the spreading of the tracer solute about its center of mass (particularly in the vertical direction), diminish the skewing of the tracer solute breakthrough, and decrease both the magnitude of the effective retardation factor and the rate at which it approaches its asymptotic value. In addition, our results suggest that under realistic conditions, the three-dimensionality of the flow domain, the periodicity of the rain or irrigation events, along with the substantial redistribution periods between successive events, and the spatial heterogeneity of the hydraulic properties of the variably saturated formation may compensate in part for the adverse effects of soil alkalinity on flow and transport on the field scale. This last finding has practical implications regarding the use of sewage water for irrigation.

Abbreviations: BTC, breakthrough curve • CDE, convection–dispersion equation • CEC, cation-exchange capacity • CP, control plane • DDL, diffuse double layer • PDF, probability density function • PSD, pore-size distribution • RSF, random space function • SAR, Na adsorption ratio

This article has been cited by other articles:

				D. Russo, J. Zaidel, and A. Laufer Numerical Analysis of Solute Transport from Trickle Sources in a Combined Desert Soil-Imported Soil Flow System Vadose Zone J., January 23, 2008; 7(1): 53 - 66. [Abstract] [Full Text] [PDF]

				G. Crescimanno and P. Garofalo Response to "Comments on 'Application and Evaluation of the SWAP Model for Simulating Water and Solute Transport in a Cracking Clay Soil'" Soil Sci. Soc. Am. J., August 22, 2006; 70(5): 1817 - 1818. [Full Text] [PDF]