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 Google Scholar Google Scholar Articles by Zhou, L. Articles by Selim, H. M. Search for Related Content PubMed Articles by Zhou, L. Articles by Selim, H. M. GeoRef GeoRef Citation Agricola Articles by Zhou, L. Articles by Selim, H. M. Related Collections Soil Kinetics Structure and Properties Soil Physics

DIVISION S-1—SOIL PHYSICS

Solute Transport in Layered Soils

Nonlinear and Kinetic Reactivity

Agronomy Dep., Sturgis Hall, Louisiana Agric. Exp. Station, LSU Agricultural Center, Baton Rouge, LA 70803-2110

* Corresponding author (mselim{at}agctr.lsu.edu) or (xp2469{at}unix1.sncc.lsu.edu)

In this study, solute transport in multilayered soils where steady flow is dominant was investigated. We considered nonreactive as well as reactive solutes in layered soils with emphasis on nonlinear solute reactivity with the soil matrix. For individual soil layers, solute-retention mechanisms considered were linear, nonlinear, Langmuir, first-, second-, and nth-order kinetics, and irreversible reactions. The convective-dispersive equation (CDE) for reactive solutes was solved using the finite difference method. First-type and a combination of first- and third-type boundary conditions (BCs) for the interface between soil layers were tested. Unlike the first-type BC, the combined first- and third-type BC always achieved a good solute mass balance in multilayered soils. Physical and chemical properties of each soil layer were assumed to differ significantly from one another. For all retention mechanisms used, our simulation results indicated that solute breakthrough curves (BTCs) were similar, regardless of the layering sequence in a soil profile. This finding is consistent with an earlier finding (Selim et al., 1977), where a linear adsorption mechanism was dominant and contrary to that of Bosma and van der Zee (1992) for nonlinear adsorption. Experimental results based on miscible displacements from soil columns of a two-layer system [sand over Sharkey (very-fine, smectitic, thermic Chromic Epiaquerts) clay and Sharkey clay over sand] support our simulation results. Specifically, BTCs for pulse inputs for tritium, as well as for the Ca–Mg system, support the above conclusion. All tritium and Ca–Mg BTCs were well predicted with our multilayered model where independently derived solute physical and retention parameters were implemented.

Abbreviations: BC, boundary condition • BTC, breakthrough curve • CDE, convective-dispersive equation • CEC, cation-exchange capacity