SSSAJ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
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 HighWire
Right arrow Citing Articles via ISI Web of Science (11)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.
Soil Science Society of America Journal 65:1585-1593 (2001)
© 2001 Soil Science Society of America

DIVISION S-1 - SOIL PHYSICS

Diffusion of Sorbing Organic Chemicals in the Liquid and Gaseous Phases of Repacked Soil

T. Olesena, J. Gamsta, P. Moldrup*,a, T. Komatsub and D. E. Rolstonc

a Dept. of Environmental Engineering, Aalborg Univ., Sohngaardsholmsvej 57, DK-9000 Aalborg, Denmark
b Dep. of Social and Environmental Engineering, Graduate School of Engineering, Hiroshima Univ., 1-4-1 Kagamiyama, Higashi-Hiroshima, 739, Japan
c Soils and Biogeochemistry, Dep. of Land, Air and Water Resources, Univ. of California, Davis, CA 95616

* Corresponding Author (i5pm{at}civil.auc.dk)


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 
Transport models for sorbing organic chemicals in soil require accurate predictions of the diffusion and sorption processes in both the liquid and gaseous phases. In this study, the ability of recently-presented diffusivity models in combination with equilibrium sorption models to predict the effective (i.e., including sorption effects) diffusion coefficient, Deff, as a function of soil-water content, {theta}, is tested for different sorbing organic chemicals in different soils. The water-induced linear reduction (WLR) gas diffusivity model, combined with a two-component (hydrophobic and vapor) equilibrium sorption model, well described short-term (<24 h) Deff({theta}) of two volatile organic compounds, toluene and trichloroethylene (TCE), in Yolo loam (fine-silty, thermic Typic Xerorthent). Further extended with the constant slope impedance factor (CSIF) solute diffusivity model, the resulting two-phase diffusion and two-component equilibrium sorption (DATES) model accurately predicted short-term (48 h) Deff({theta}) of a semivolatile pesticide, lindane (C6H6Cl6), in Gila silt loam (coarse-loamy, thermic Typic Torrifluvents). To test the DATES model for longer-term data, Deff({theta}) of naphthalene (C10H8) in Lerbjerg sandy clay (L5) was measured at different soil-water contents and incubation times in air-tight diffusion cells. Because of sorption kinetics, both the apparent hydrophobic sorption coefficient, KD, derived from the column diffusion experiments and KD derived from batch desorption experiments increased with naphthalene-soil contact time for t > 700 h. Consequently, an increasing KD value with time was required in the DATES model to obtain successful Deff({theta}) predictions. Since chemical-soil contact time and adsorption-desorption history typically vary in each compartment of a soil profile, transport models, including sorption kinetics, will become very complex and DATES with a time-dependent KD relation may represent a useful model alternative.

Abbreviations: CSIF, Constant slope impedance factor • Deff, effective diffusion coefficient • DATES, Diffusion and two-component equilibrium sorption • L5, Lerbjerg sandy clay • TCE, Trichloroethylene • WLR, Water-induced linear reduction


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 
DIFFUSION of sorbing organic chemicals in the soil liquid and gaseous phases affects and potentially controls the spreading and fate of organic pollutants in the soil. Recently, Moldrup et al. (2000a) and Olesen et al. (2001) suggested improved models to predict the diffusion coefficients of nonsorbing tracers in the gaseous and liquid phases, respectively, of repacked soil. For sorbing organic chemicals, however, the effective diffusion rate is determined partly by the diffusion properties and partly by the sorption of the organic chemicals in the soil gaseous and liquid phases.

Hydrophobic sorption onto soil organic matter is described for many chemicals and soil types (e.g., Karickhoff et al., 1979; Abdul and Gibson, 1986; Weber et al., 1992; Pignatello and Xing, 1996), while studies of vapor sorption onto dry soil mineral surfaces are limited to very few different chemicals and in most cases only TCE (e.g., Petersen et al., 1994, 1995, 1996a; Ong and Lion, 1991a,b; Shimuzu et al., 1992, 1994). Consequently, several predictive expressions to estimate hydrophobic adsorption coefficients from soil and chemical characteristics are available (e.g., Briggs, 1973; Karickhoff, 1981; Abdul et al., 1987), while it remains to be verified that recent predictive expressions for vapor sorption of TCE (e.g., Petersen et al., 1995; Poulsen et al., 1998, 2000; Yamaguchi et al., 1999; Chen et al., 2000) can be applied for other chemicals.

A correct description of diffusion and sorption processes in both the soil liquid and gaseous phases is required to accurately predict the effective diffusion of a sorbing organic chemical (Ryan and Cohen, 1990). In this study, a combined diffusion and sorption theory is derived, yielding a two-phase diffusion and two-component equilibrium sorption (DATES) model. The diffusion in the gaseous and liquid phase is predicted using recently presented gas and solute diffusivity models by Moldrup et al. (2000a) and Olesen et al. (2001). Sorption from the liquid (hydrophobic) and gaseous (vapor) phase need to be experimentally determined or estimated from empirical expressions. The DATES model is tested against short term (<=48 h) data from the literature for volatile (TCE and toluene) and semivolatile (lindane) organic chemicals. Finally, the DATES model is used to describe longer-term data for naphthalene (semivolatile), measured in this study, to evaluate possible effects of sorption kinetics.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 
Radiolabeled 14C naphthalene (Sigma Chemical Co., St. Louis, MO) with a specific activity of 1.16 x 109 Bq mmol-1 (31.3 mCi mmol-1) was used. A stock solution containing 14C naphthalene with a specific activity of 1.17 x 109 Bq L-1 (31.5 mCi L-1) was prepared in methanol (128 mg L-1) and a stock solution containing 10 µg cm-3 unlabeled naphthalene in 0.01 M CaCl2 was prepared for the adsorption or desorption, and diffusion experiments. Sieved (<2 mm) L5 was used for all experiments. Soil physical properties for L5 and chemical properties for naphthalene are shown in Tables 1 and 2, respectively. To avoid microbial degradation, the soil and the water used for naphthalene solutions were autoclaved for 2 h prior to each experiment. The autoclaving procedure did not affect the sorption results (not shown).


View this table:
[in this window]
[in a new window]
 
Table 1. Soil physical properties.

 

View this table:
[in this window]
[in a new window]
 
Table 2. Chemical properties.

 
Naphthalene Adsorption
Adsorption experiments were carried out at 10°C using a batch technique. Solution to solid ratios were chosen to ensure ~50% sorption. Five grams of air dried L5 soil (water content ~0.045 g g-1) were weighed into 36-mL glass centrifuge tubes, and 30 mL of naphthalene in 0.01 M CaCl2 solution were added, and subsequently a Teflon cap was used to seal the glass tube. Three replicates were made for each naphthalene concentration (0.1–5 µg cm-3). The vials were gently rotated vertically (50 rpm) for 3 wk. Two 1-mL samples were taken from the supernatant after centrifuging for 10 min at 2665 x g, then transferred to glass scintillation vials, and 10 mL of Packard Ultima Gold XR (Packard Bioscience, Groningen, The Netherlands) scintillation counting liquid added. The scintillation vials were shaken for 1 min and after 24 h analyzed by liquid-scintillation counting (Packard 1600 TR, Packard Instruments, Downers Grove, IL) and counted for 3 min. For each concentration, two controls without soil material, but otherwise treated similarly, were included. It was assumed that the difference between the solution concentration in vials with soil and without soil was because of adsorption only. Sorption to centrifuge tubes and Teflon caps was small (<2%) and was taken into account by assuming that the amount of adsorbed naphthalene to glass and cap was the same whether the vial was with or without soil.

Naphthalene Desorption
The samples with the initial naphthalene concentration of 5.0 µg cm-3 from the adsorption experiment were continued for 12 (2-d) desorption steps. Prior to each desorption step, the samples were centrifuged at 2665 x g for 10 min, and 20 mL of the supernatant plus two 1-mL samples for liquid-scintillation counting was replaced by 22 mL of 0.01 M CaCl2. The amount of naphthalene adsorbed to the soil after each desorption step was calculated as the initial amount of naphthalene minus the amount removed during previous desorption steps and minus the amount present in the liquid phase. At the end of the desorption experiments, the soil was extracted with hexane to determine the mass balance. The samples were gently rotated (50 rpm) for 24h, and two 1-mL samples were removed for liquid-scintillation counting. This extraction procedure was carried out for up to 5 successive extraction steps. The average recovery was 98.2%.

Naphthalene Effective Diffusion
The effective diffusion of naphthalene was measured in L5 at a bulk density of 1.35 g cm-3 at eight different soil-water contents, ranging from 0.05 to 0.34 cm3 cm-3 (from air-dry soil to soil at a soil-water potential of ~-100 cm H2O). Air-dry soil was mixed with water to obtain the desired water content. Subsequently, 14C labeled (466 Bq per g dry soil; 0.0126 µCi per g dry soil) and unlabeled naphthalene were added in hexane to obtain a concentration of ~5 µg cm-3 soil-water which was ~50 µg g-1 dry soil. The hexane was allowed to volatilize off before the soil was packed in air-tight (tested over the range of -0.09 to 1.3 MPa [-0.9–13 bar]) 8-cm aluminum half-cells, 2 cm at a time. Between the packing of each 2-cm layer, the surface was disturbed to obtain good contact between the layers. The half-cells were sealed with aluminum lids, and allowed to equilibrate for 10 to 14 d depending on the soil-water content. For the three highest soil-water contents, the half cells were packed with air-dry soil, and the required amount of water was added with a pipette after packing each 1-cm layer. These half-cells were allowed 2 to 3 wk of equilibration before incubation. After the equilibration period, the source and recipient half-cells (cells without addition of naphthalene) were connected and incubated horizontally at 10°C. The incubation time varied between 5 d and 6 wk depending on the water content, to obtain sufficient spreading of the concentration versus distance profile.

The cells were subsequently sliced in 0.5-cm slices (within 0–4 cm on each side of the interface) and 1.0-cm slices (within 4–8 cm distance from the interface). Between one and three replicate 1-g soil samples were transferred to glass scintillation vials, and 15 mL of Packard Ultima XR Gold scintillation counting liquid (Packard Bioscience, Groningen, The Netherlands) was added. The vials were shaken and then allowed to settle for 24 h, and the samples were then counted with 24-h intervals until maximum counts were obtained. The water content was measured in three samples from each half-cell. At the highest water content (0.34 cm3 cm-3) 3-cm half-cells were used. The cells were packed in 1-cm layers, and sectioned in 0.3- and 0.5-cm slices after incubation.


    THEORY
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 
The gas and solute diffusion coefficients are governed by the length, geometry, and connectivity of the diffusive pathways in the gaseous and liquid phases, and are often predicted from empirical equations for nonsorbing chemicals (e.g., freon-12 and Cl). For predicting diffusion in the soil gaseous phase, a model for gas diffusivity in sieved, repacked soils was proposed by Moldrup et al. (2000a), the so-called water-induced linear reduction (WLR) model,

[1]
where DP,g is the gas diffusion coefficient in soil (cm3 soil air cm-1 soil s-1), D0,g is the gas diffusion coefficient in free air (cm2 air s-1), {Phi} is the soil total porosity (cm3 pores cm-3 soil), and {epsilon} is the air-filled porosity (cm3 soil air cm-3 soil). Equation [1] is used in this study where all data are for repacked soil columns. A macroporosity dependent model for gas diffusivity in undisturbed soil was developed by Moldrup et al. (2000b) and should be used instead of Eq. [1] in the case of undisturbed soil systems. The second term ({epsilon}/{Phi}) in Eq. [1] takes into account the increased tortuosity (reduced diffusivity) in a wet soil compared with a dry soil at the same air-filled porosity. The WLR model, Eq. [1], accurately predicted gas diffusivity as a function of soil air contents in 12 repacked soils with different soil textures and total porosities (Moldrup et al., 2000a).

For predicting diffusion in the soil-liquid phase, a soil-type dependent solute diffusivity model, the so-called CSIF model, was proposed by Olesen et al. (2001),

[2a]

[2b]
where DP,l is the solute diffusion coefficient in soil (cm3 soil water cm-1 soil s-1), D0,l is the solute diffusion coefficient in water (cm2 water s-1), {theta} is the soil-water content (cm3 soil water cm-3 soil), and b is the Campbell (1974) pore-size distribution parameter. The term 0.02b in Eq. [2] describes the threshold soil-water content, {theta}th, where the diffusion ceases because of disconnection of the diffusive pathway, and the constant (1.1) is a factor describing the meandering of the diffusive pathway. Olesen et al. (2001) suggested that the term 0.02b can be replaced by an empirical soil texture and bulk density dependent expression if the soil-water characteristic curve is not known. Equation [2] is assumed valid for both sieved, repacked, and for undisturbed soil since no significant difference in measured diffusivities for repacked and for undisturbed soil was observed by Olesen et al. (2000) and Barraclough and Tinker (1982). Equation [2] well predicted solute diffusivities for 22 repacked soils with different soil textures (Olesen et al., 2001).

The combined diffusion of a sorbing chemical in the soil gaseous and liquid phases, assuming constant {epsilon} and {theta}, can be described by,

[3]

[4]
where F is the chemical flux (µg cm-2 soil s-1); Ctotal is the total chemical concentration (µg cm-3 soil); t is time (s); z is soil depth (cm soil); Deff is the effective diffusion coefficient (cm2 soil s-1), and Rg and Rl are the retardation factors (due to sorption) for the gas and liquid phases, respectively.

Assuming immediate equilibrium among the soil phases, where Henry's constant, KH (cm3 soil water cm-3 soil air), describes the distribution between the liquid and gas phase, and Kd (cm3 soil water g-1 soil) is the liquid/solid partition coefficient (assuming linear sorption/desorption behavior), the retardation factors are,

[5a]

[5b]
where {rho}b is the soil bulk density (g soil cm-3 soil).

At low soil-water contents, i.e., typically below a water content corresponding to three to five molecular layers of water coverage, the adsorption capacity of soils increase drastically because of chemical adsorption from the gas phase onto dry mineral surfaces (Ong and Lion, 1991a,b; Petersen et al., 1994, 1995), and Henry's law and Eq. [5] are not applicable (this is referred to as the non-Henry range in the following). Instead the adsorption can be described by the combined vapor sorption (soil mineral surface–gas phase partitioning) and hydrophobic sorption (soil organic matter–aqueous phase partitioning) coefficient, K'D (cm3 soil air g-1 soil), describing the partitioning of the chemical between the gas phase and the rest of the soil system,

[6]
where Cg is the concentration in the gaseous phase (µg cm-3 soil air), Cl is the concentration in the liquid phase (µg cm-3 soil water), and CS is the sorbed concentration (µg g-1 soil). Including vapor sorption, the retardation factors can be expressed as,

[7a]

[7b]

The partitioning coefficient K'D can be expressed as a function of the water content by,

[8]
where Ksg is the solid/vapor phase partitioning coefficient (cm3 soil air g-1 soil), w is the gravimetric water content (g H2O g-1 dry soil), {phi} is the aqueous activity coefficient [equal to 1], and {rho} is the water density (equal to 1 g cm-3). It should be noted that Eq. [8] is valid only at low chemical vapor pressures, where chemical vapor sorption is linearly related to the vapor pressure (Amali et al., 1994). At high soil-water contents where chemical vapor sorption is negligible (Ksg = 0), Eq. [8] inserted in Eq. [7] yields Eq. [5]. Thus, Eq. [7] is generally valid in both the dry (nonHenry) and wetter (Henry) soil-water content regions.

The retardation factors given by Eq. [7] combined with Eq. [8] can be used to predict K'D if Ksg(w) is known. By rewriting Eq. [8], Petersen et al. (1995) proposed a semi-empirical model that can be used to describe Ksg(w) and K'D,

[9]
where K'D is the value of K'D at zero moisture content, {alpha} describes the decrease in K'D with increasing w in the low water content range, and the last term of the equation describes the hydrophobic sorption in the high water content (Henry) range. In the first (vapor sorption) term, K'D and {alpha} can be obtained either from laboratory experiments (e.g., Petersen et al., 1995, 1996a) or predicted from soil texture and surface area by empirical expressions (e.g., Poulsen et al., 1998, 2000; Yamaguchi et al., 1999).

Equation [4] together with Eq. [1], [2], [7], and [9] represent a two-phase diffusion and two-component equilibrium sorption (DATES) model, that can be used to predict the effective diffusion coefficient, Deff, of a sorbing (hydrophobic and vapor) chemical in sieved, repacked soil.

Assuming constant {epsilon} and {theta}, the flux equation, Eq. [3], combined with the continuity equation yields the governing differential equation for one-dimensional diffusion of a sorbing chemical,

[10]
where the soil chemical concentration Ctotal = Cg{epsilon} + Cl{theta} + CS{rho}b, and Deff is given by Eq. [4].

Predicting Short-term (<=48 h) Diffusion of Volatile and Semivolatile Chemicals
The WLR gas diffusivity prediction model, Eq. [1], in combination with the two- component (hydrophobic and vapor) equilibrium sorption model, Eq. [4], [7] and [9], is first tested against data from Petersen et al. (1994) for TCE and toluene effective diffusion coefficients in Yolo loam soil (Fig. 1) . For these two volatile chemicals, the contribution from liquid phase diffusion on Deff({theta}) is negligible. The effective TCE and toluene diffusion coefficients were measured using a two-chamber diffusion apparatus consisting of two air-filled chambers separated by a packed soil column (7.57-cm diameter and 10- or 20-cm length). The incubation time was <24 h. Soil physical properties for Yolo loam and chemical properties for TCE and toluene are shown in Tables 1 and 2, respectively. From independent batch sorption experiments, Petersen et al. (1994)(1995) found KD = 0.58 cm3 g-1, K'D = 3401 cm3 g-1, and {alpha} = 37.5 for TCE and KD = 0.87 cm3 g-1, K'D = 21678 cm3 g-1, and {alpha} = 34.3 for toluene. Figure 1 shows that the effective diffusion coefficients of Deff({theta}) for TCE and toluene are well described by the combined gas diffusion (WLR) and two-component (hydrophobic and vapor) sorption model, when the independently measured sorption characteristics [KD, K'D, and {alpha}] are used. The average relative prediction error is 15% for both TCE and toluene, excluding the measurement at the lowest water content. At the lowest water content (air-dry soil) there is a relatively large deviation between measured and predicted values. Some of this deviation is because of uncertainty in determining the exact water content for air-dry soil. To illustrate that the WLR part, Eq. [1], of the DATES model is accurate, Fig. 1 shows that the WLR model well predicts the effective diffusion coefficients for freon 12 (nonsorbing).



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 1. Measured freon-12, TCE, and toluene effective diffusion coefficients in Yolo loam, compared with predictions by the WLR model, Eq. [1], in combination with retardation factors calculated from Eq. [7]. For freon-12 Rg = Rl = 1. Data from Petersen et al. (1994).

 
The simultaneous application of both the gas and the solute diffusivity prediction models (Eq. [1] and [2]) in the full DATES model (Eq. [4] plus [1], [2], [7], and [9]) is verified for the semivolatile pesticide lindane, where Deff({theta}) were measured in Gila silt loam by Ehlers et al. (1969) given by Fig. 2 . The effective diffusion coefficients were measured using a half-cell method with a source half-cell containing soil treated with 14C labeled lindane and a recipient half-cell containing untreated soil (the half-cell diameter and length were 1.5- and 4.5-cm, respectively). The incubation time was 48 h. Soil physical properties for Gila silt loam and chemical properties for lindane are shown in Tables 1 and 2, respectively. Since Campbell b was not available for Gila silt loam, the term 0.02b in Eq. [2] is replaced by the soil texture and bulk density dependent expression suggested by Olesen et al. (2001) (yielding 0.077 for Gila silt loam). The value of KD is estimated from the fraction of organic C, foc, and the organic C partitioning coefficient, Koc (Table 1) as KD = foc Koc. Poulsen et al. (1998) found {alpha} = 89 for TCE diffusion in Gila silt loam. As Petersen et al. (1995) and Poulsen et al. (2000) suggested that {alpha} depends only on the soil characteristics, the same value was assumed valid for lindane. However, as noted by Poulsen et al. (2000), more data are needed to verify that {alpha} is indeed independent of the chemical properties. In this case, for lindane sorption onto Gila silt loam K'D was not available (measured). Therefore the model was fitted (least square optimization) to measured data in the non-Henry range yielding a K'D value of 0.7 x 108. Figure 2 shows that the effective diffusivities for lindane at 20°C and 30°C are well described by the DATES model. The diffusion in the liquid phase and the gaseous phase is of equal magnitude at a soil-water content around 0.25 cm3 cm-3, emphasizing that an accurate description of both solute and gas diffusivities is required in the case of semivolatile chemicals. Hydrophobic sorption decreases the effective diffusion by up to approximately two orders of magnitude (as compared with effective diffusion in the gas and liquid phases without taking sorption into account). Vapor sorption further decreases Deff by up to three orders of magnitude in the low water content region.



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 2. Comparison of measured lindane effective diffusion in Gila silt loam and predictions by the two-phase diffusion and two-component equilibrium sorption (DATES) model. The influence of the separate diffusion and sorption processes on Deff is illustrated. Data from Ehlers et al. (1969). (a) 20°C; (b) 30°C.

 
Predicting Naphthalene Diffusion
To test the DATES model for other than the short-term data (<=48h for TCE, toluene, and lindane in the previous section), the DATES model was compared with the data for effective diffusion of naphthalene (semivolatile) in L5 measured in this study. As a preliminary test, Fig. 3 shows diffusivities in L5 for nonsorbing tracers in the liquid and gaseous phases [data from Olesen et al. (1999) and Moldrup et al. (2000a)]. It is obvious that Eq. [1] and [2] provide a good description of the measured gas and solute diffusivities for the L5 soil.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 3. Measured and predicted diffusivities for nonsorbing tracers (O2 and Cl) in Lerbjerg sandy clay (L5).

 
Naphthalene diffusion in L5 was measured in air-tight diffusion cells at eight different soil-water contents ranging from air-dry soil to a soil-water potential of {Psi}~-100 cm H2O (see Materials and Methods). Figure 4 shows examples of measured concentration versus distance profiles. A numerical solution to Eq. [10] (Moldrup et al., 1996) was used to calculate the effective naphthalene diffusion coefficients from the measured concentration profiles. The corresponding calculated concentration profiles are shown as solid lines. For comparison, predicted profiles based on diffusion coefficients of nonsorbing tracers in L5 soil (Fig. 3) and retardation factors calculated from Eq. [5] (using the value of KD = 9.5 cm3 g-1 obtained from the adsorption isotherm, Fig. 5) , are shown as dotted lines. This is obviously not a realistic model approach. In Fig. 4a and 4b, the predictions fail since the water content is low (nonHenry range) and vapor sorption is not taken into account . Petersen et al. (1996b) suggest that for fine to medium textured soils, the soil-water content at wilting point (at {Psi} = -15000 cm H2O) can be used as an indication for the soil-water content below which vapor sorption becomes important. For L5, this corresponds to {theta} = 0.17 cm3 cm-3. The profiles in Fig. 4c and 4d are reasonably well predicted since the water content is higher ({theta} > 0.17 cm3 cm-3) and the incubation time is relatively short. By comparing Fig. 4d and 4e representing cores at the same {theta} but different incubation times, it appears that the poor predictions for the long incubation time (Fig. 4e) must be because of adsorption and desorption kinetics (aging). The same is the case for Fig. 4f. Thus, Fig. 4 implies that both vapor sorption and sorption kinetics (aging) needs to be taken into account to obtain a correct description of chemical mobility.



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 4. Examples of measured concentration versus distance profiles for naphthalene diffusion at different soil-water contents ({theta}) and incubation times (t). The solid lines are model fit to measured data to obtain the Deff value, while the dotted lines represent model predictions using measured gas and solute diffusivities (Fig. 3) and adsorption KD (equal to 9.5 cm3 g-1).

 


View larger version (24K):
[in this window]
[in a new window]
 
Fig. 5. Twenty-one day naphthalene adsorption isotherm for Lerbjerg sandy clay (L5) and subsequent 2 d desorption steps. The adsorption is approximated by a linear isotherm . A linear isotherm is plotted for the first and the last desorption step corresponding to an apparent KD value between 13 cm3 g-1 and 93 cm3 g-1 for contact times (adsorption plus desorption time) between 552 and 1080 h.

 
Figure 5 shows the 21-d naphthalene adsorption isotherm for L5, approximated by a linear isotherm and yielding a KD value of 9.5 cm3 g-1. Twelve desorption steps after the 21-d adsorption are shown, and adsorption or desorption hysteresis is evident. Apparent desorption KD values (KD,app) increased with contact time (adsorption time plus desorption time) from KD,app = 13 cm3 g-1 at t = 552 h (21-d adsorption and 2-d desorption) to KD,app = 93 cm3 g-1 at t = 1080 h (21-d adsorption plus 24-d desorption).

For comparison, apparent KD values can be found from each column diffusion experiment as the KD value in Eq. [5] that combined with the measured diffusivities of nonsorbing tracers (Fig. 3) yields the best model fit of the measured concentration versus distance profile (assuming instantaneous linear sorption equilibrium). In Fig. 6 , the apparent KD values are shown as a function of the contact time (equal to equilibration plus incubation time) for the experiments performed in the Henry range ({theta} > 0.17 cm3 cm-3). It is obvious that the apparent KD is fairly constant for short-term experiments (contact time <700 h), while an increase is seen for long-term experiments. Similarly White et al. (1999) showed that increased aging slowed the initial desorption rate of phenanthrene (C14H10) and also increased the resistant fraction (this would correspond to an increase in the apparent KD if diffusion experiments were carried out). Alexander (1995) suggested that changes in extractability with residence time in soil is because of sequestration in inaccessible microsites within the soil matrix.



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 6. Apparent KD value as a function of contact time (equilibration plus incubation time) for the column diffusion experiments. Only experiments in the Henry range are included ({theta} > 0.17 cm3 cm-3). For comparison, apparent KD values from the desorption batch experiments (contact time equals adsorption plus desorption time) are shown.

 
Two things are noticed in Fig. 6. First, apparent KD values from column and batch experiments are not the same. The reason for this is that the adsorption and desorption history typically is different in each compartment within a soil column. This illustrates the complexity of the problem since a transport and fate model would then need to include the adsorption and desorption history at each point in the soil profile. Second, the increase in KD with time seems fairly linear between 700 and 1400 h of contact time. We cannot be sure if this trend will continue at higher contact times, however, the results imply that a simple linear increase in KD with contact time may be a reasonable approximation.

In Fig. 7 , the measured effective diffusion coefficients are compared with predictions by the DATES model. The apparent KD value derived from the short-term column experiments (17 cm3 g-1, Fig. 5) was used. The K'D was expressed by Eq. [9], where the value of {alpha} = 30 was estimated from the clay content using an empirical expression (based on vapor sorption of TCE) suggested by Yamaguchi et al. (1999). The value of K'D was fitted to 3.5 x 105 to provide a good description of Deff({theta}) in the nonHenry range. In accordance with Fig. 5, Fig. 7 shows that the effective diffusion coefficients are largely overestimated for long-term experiments. Using apparent KD values between 35 and 70 cm3 g-1 gave a fair description of Deff for the long-term experiments.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 7. Effects of apparent KD on DATES predicted effective diffusion coefficients for naphthalene in Lerbjerg sandy clay(L5). The influence of each diffusion and sorption process on Deff is illustrated.

 
It is obvious from Fig. 7, as well as from Fig. 4 through 6, that it is important to gain a more detailed knowledge of sorption kinetics. Furthermore, it is clear that a dynamic description of the sorption and desorption process is required to obtain a conceptually realistic description of the effective diffusion of naphthalene in L5. The introduction of a time-dependent apparent KD relation in the DATES Deff({theta}) model may represent a useful alternative to a more complex, kinetics-based model, and provides a good preliminary evaluation of the effects of sorption kinetics on chemical mobility.

The importance of each of the diffusion and sorption processes are compared in Fig. 7, showing that diffusion in the liquid phase is unimportant as compared with diffusion in the gaseous phase except at water contents very close to water-saturation. This is in contrast to lindane diffusion in Gila silt loam, Fig. 2, where diffusion in the liquid phase was shown to be very important at {theta} > 0.25 cm3 cm-3. Hydrophobic sorption decreases the effective diffusion by two to three orders of magnitude as compared with diffusion of a nonsorbing compound, which is similar to what was seen for lindane (Fig. 2). Vapor sorption further decreases Deff up to around two orders of magnitude. Vapor sorption is important in a wider soil-water content range for L5 as compared with Gila silt loam, because of the difference in soil mineral surface area. This is also evident from the different {alpha} values of 30 and 89 for L5 and Gila soil, respectively.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 
Recently presented gas and solute diffusivity models were combined with hydrophobic and vapor sorption models yielding a two-phase diffusion and two-component equilibrium sorption (DATES) model. The DATES model successfully described short-term (<=48h) effective diffusion data from the literature for volatile (TCE and toluene) and semivolatile (lindane) organic chemicals.

For longer-term experiments with naphthalene diffusion in a sandy clay, the effective diffusion could not be accurately predicted using a constant value of the hydrophobic sorption coefficient KD in DATES. Because of sorption kinetics, both the apparent hydrophobic sorption coefficient, KD, derived from the column diffusion experiments and KD obtained from batch desorption experiments increased with naphthalene-soil contact time. Since sorption history is typically different in each compartment of a soil-profile, results for KD(t) from batch and diffusion column experiments were different. Thus, a both time-dependent (kinetics) and space-dependent (sorption history) based modeling seems required to get a fully accurate description of the effective diffusive transport.

The DATES model is suggested as a model platform for further development of a two-phase, diffusive-reactive chemical transport model, and is useful to evaluate the relative importance of each of the diffusion (gaseous and liquid phase) and sorption (hydrophobic and vapor) processes for different combinations of chemicals and soil types. Using a time-dependent KD term in DATES may provide a good preliminary evaluation of the possible effects of sorption kinetics on chemical mobility.


    ACKNOWLEDGMENTS
 
This study was funded by the Danish Technical Research Council, Research Talent Project entitled: ‘New methods for measuring and predicting liquid and gaseous phase transport properties in undisturbed soils’, grant 5P42ESO4699 from the National Institute of Environmental Health Sciences, NIH, and the U.S. EPA (R819658) Center for Ecological Health Research at University of California–Davis. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official view of the NIEHS, NIH, or EPA. The authors gratefully acknowledge a research grant from the Japanese Ministry of Education, Science, Sports, and Culture (Monbushu International Scientific Research Program: Joint Research No. 12555156).

Received for publication October 16, 2000.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 THEORY
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Soil Sci.Home page
W. S. D. Rocha, J. B. Regitano, and L. R. F. Alleoni
2,4-D Residues in Aggregates of Tropical Soils as a Function of Water Content
Soil Sci. Soc. Am. J., October 27, 2006; 70(6): 2008 - 2016.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
J. Gamst, P. Moldrup, D. E. Rolston, T. Olesen, K. Scow, and T. Komatsu
Comparison of Naphthalene Diffusion and Nonequilibrium Adsorption-Desorption Experiments
Soil Sci. Soc. Am. J., May 1, 2003; 67(3): 765 - 777.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
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 HighWire
Right arrow Citing Articles via ISI Web of Science (11)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Olesen, T.
Right arrow Articles by Rolston, D. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Vadose Zone Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome