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 Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (3) Citing Articles via Google Scholar Google Scholar Articles by Haws, N. W. Articles by Rao, P. S. C. Search for Related Content PubMed Articles by Haws, N. W. Articles by Rao, P. S. C. GeoRef GeoRef Citation Agricola Articles by Haws, N. W. Articles by Rao, P. S. C. Related Collections Dual Porosity/Permeability Models Nonequilibrium Transport Vadose Zone Processes and Chemical Transport

ORIGINAL RESEARCH

The Effect of Vertically Decreasing Macropore Fractions on Simulations of Non-Equilibrium Solute Transport

Department of Agronomy, Purdue University, West Lafayette, IN 47907-2054
^* Corresponding author (nhaws{at}jhu.edu)

Received 7 January 2004.

The effects of a depth-decreasing macropore fraction and the corresponding decrease in the diffusive mass-transfer coefficient for solute exchange between soil matrix and macropore regions were simulated for a hypothetical, one-dimensional, mobile–immobile, transport domain. These effects are illustrated using an expression for a Damkhöler number, and with numerical simulations with a dual-porosity module within the HYDRUS-2D code. The depth-wise change in the macroporosity was represented with an exponentially decreasing macropore fraction, based on a pore-network model. Five simulations were run, four with various steady infiltration rates (two at water saturation and two at different levels of unsaturation) and one with transient infiltration rates. Deviations between the predicted solute breakthrough curves (BTCs) for simulations using depth-variable parameters and those using spatially averaged parameters were minor at the lowest infiltration rate, indicating little effect of non-equilibrium conditions. However, at higher infiltration rates, the BTC peak height and peak timing varied by as much as a factor of two, indicating significant influence of the variability of the macropore fractions with depth. Using depth-averaged values to represent the variable macropore fraction will closely approximate the BTC of a spatially variable case for both steady-state and transient infiltration; however, the concentration profile within the transport domain will not be accurately predicted. A sensitivity analysis conducted for the saturated hydraulic conductivity, length of the soil ped, and dispersivity parameters shows that changes in the prediction error between the spatially variable and surface-valued simulations are most sensitive to these parameters as the soil approaches water saturation.

Abbreviations: BTC, breakthrough curve

This article has been cited by other articles:

				T. B. Ramos, M. C. Goncalves, J. C. Martins, M. Th. van Genuchten, and F. P. Pires Estimation of Soil Hydraulic Properties from Numerical Inversion of Tension Disk Infiltrometer Data Vadose Zone J., May 26, 2006; 5(2): 684 - 696. [Abstract] [Full Text] [PDF]