Published online 5 January 2006
Published in J Environ Qual 35:268-276 (2006)
DOI: 10.2134/jeq2005.0140
© 2006 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Organic Compounds in the Environment
Sorption Kinetics and Equilibria of Organic Pesticides in Carbonatic Soils from South Florida
P. Nkedi-Kizzaa,*,
D. Shindea,
M. R. Savabib,
Y. Ouyangc and
L. Nievesb
a Soil and Water Science Department, University of Florida, Gainesville, FL 32611
b USDA-ARS, Subtropical Horticulture Research Station, Miami, FL 33158
c St. Johns River Water Management District, P.O. Box 1429, Palatka, FL 32178-1429
* Corresponding author (Kizza{at}ifas.ufl.edu)
Received for publication April 26, 2005.
 |
ABSTRACT
|
|---|
A batch reactor was used to determine sorption kinetic parameters (k2, F, and K*) and the equilibrium sorption coefficient (K). The two-site nonequilibrium (TSNE) batch sorption kinetics model was used to calculate the kinetic parameters. Two probe organic pesticides, atrazine [2-chloro-4-ethylamino-6-isopropylamino-s-triazine] and diuron [3-(3,4-dichlorophenyl)-1,1-dimethylurea] were studied using three carbonatic soils from South Florida (Chekika, Perrine, and Krome), one noncarbonatic soil from Iowa (Webster), and one organic soil (Lauderhill) from South Florida. Carbonatic soils contained more than 600 g kg1 CaCO3. Sorption is initially very fast up to 3 h and then slowly reaches equilibrium. All soilchemical combinations reached sorption equilibrium after about 24 h and all sorption isotherms were linear. The sorption kinetics data were well described by the TSNE model for all soil-chemical combinations except for the marl soil data (PerrineAtrazine), which were better described by the one-site nonequilibrium (OSNE) model. Diuron, with higher K, undergoes slower sorption kinetics than atrazine. The Lauderhill soil containing organic carbon (OC) of 450 g kg1 exhibited slowest sorption kinetics for both pesticides. An inverse relationship between k2 and K was observed for atrazine and diuron separately in Chekika, Webster, and Lauderhill soils but not in Perrine and Krome soils. The sorption kinetic parameters were used to distinguish the sorption behavior between atrazine and diuron and to identify differences between soils. Normalizing the sorption coefficient (K) to OC showed that atrazine and diuron had Koc values in carbonatic soils that were a third of reported literature values for noncarbonatic soils. Using existing literature Koc values in solute transport models will most likely underestimate the mobility of atrazine, diuron, and other neutral organic chemicals in carbonatic soils.
Abbreviations: OC, organic carbon OSNE, one-site nonequilibrium TSNE, two-site nonequilibrium
 |
INTRODUCTION
|
|---|
SURFACE AND GROUND WATER contamination resulting from application of agricultural chemicals to soils is a national environmental concern. Assessment of pesticide sorption during transport is a prerequisite for minimizing their potential mobility in the vadose zone (Selim, 2004). In South Florida, approximately 13 230Mg of pesticides are used each year, mainly in agriculture (Miles and Pfeuffer, 1997). In a nationwide survey of pesticide loading to coastal regions, Miami-Dade County in South Florida ranked second in agricultural application (Pait et al., 1992). The Miami-Dade County agricultural area is of particular importance due to its proximity to the Everglades National Park (Downing et al., 2004). Of particular concern is the shallow soil depth (1560 cm) of major agricultural soils that provides only a thin protective layer of soil to the underlying high water table. Very little data are available in the literature for the sorption kinetics and equilibria parameters of pesticides in carbonatic soils. Sorption kinetics and equilibria parameters play an important role in determining the fate and transport of pesticides in the environment. Application practices for pesticides, combined with the high annual rainfall in South Florida, could result in significant amounts of pesticides reaching nontarget sites.
Atrazine is one of the most widely detected herbicides in surface and ground water in the United States. It is usually applied to the soil surface and its fate in soil depends on the interaction between its molecules and soil components (Ben-Hur et al., 2003). Diuron is among the top 10 pesticides frequently detected in shallow ground water and is characterized as a potential threat to ground water resources due to its moderate water solubility, low tendency to sorb to soils, and relatively long half-life (Field et al., 2003). The drinking water maximum concentration level (MCL) for atrazine is 3 µg L1 (USEPA, 1998). Although no MCL is established for diuron in the United States, the health advisory limit is 10 µg L1 (USEPA, 1987).
Two of the major factors known to influence sorption of pesticides are soil properties and the molecular characteristics. Sorption of neutral organic pesticides, such as atrazine and diuron, depends primarily on soil organic carbon (OC) content (Madhun et al., 1986). Sorption determines whether the pesticide will persist, be transported, and pollute the underlying ground water (Rao and Hornsby, 1983). Strongly adsorbed and persistent pesticides that have large (Koc) values are likely to remain near the soil surface. In contrast, weakly adsorbed but persistent pesticides (small Koc) may be readily leached through the soil and are more likely to contaminate ground water. Atrazine and diuron are moderately persistent with half-life of 60 d for atrazine and 90 d for diuron and Koc values are 100 for atrazine and 400 for diuron (Nkedi-Kizza et al., 1985; Wauchope et al., 1992).
Experimental investigations of sorption processes are often conducted in batch systems rather than dynamic systems, to distinguish the sorption process from hydrodynamic dispersive effects (Pedit and Miller, 1994). Batch experiments are also used to obtain data for rates of intermediate and slow sorption (Wu and Gschwend, 1986). Batch experiments utilize static conditions in which forward and reverse reactions are allowed to proceed until an equilibrium condition is established (Harter, 1991). Modeling efforts focusing on describing the retention mechanisms of pesticides in the soil environment have generally been based on results from batch sorption experiments. Batch data based on 24-h adsorption are commonly regarded as a measure of equilibrium-type retention. Kinetic data, which are measured infrequently, have the advantage of taking into account possible time-dependent reactions for adsorption, release, or desorption. Nonequilibrium conditions may be caused by the heterogeneity of sorption sites and slow diffusion to sites within the soil matrix, or organic matter (Brusseau et al., 1991; Selim, 2004).
A variety of chemical and physical nonequilibrium processes affect sorption and transport of pesticides in the subsurface (Nielsen et al., 1986; Aharoni and Sparks, 1991). Chemical nonequilibrium models consider adsorption on some of the adsorption sites to be instantaneous, while adsorption on the remaining sites is governed by first-order kinetics (Selim et al., 1976; Cameron and Klute, 1977; Van Genuchten, 1981; Nkedi-Kizza et al., 1984; Ma and Selim, 1997).
The main objective of this study was to determine sorption kinetics from batch experiments for neutral organic chemicals (i.e., atrazine and diuron) and to ascertain sorption characteristics in carbonatic and noncarbonatic soils. Such information is necessary for understanding soil pesticide mobility through runoff and leaching and for preventing potential contamination of water resources.
 |
MATERIALS AND METHODS
|
|---|
Soils
Carbonatic Soils
The Chekika gravelly loam (loamy-skeletal, carbonatic, hyperthermic Lithic Udorthents) series consists of very shallow, somewhat poorly drained, moderately permeable soils over limestone. Soils formed by scarification of outcrops of oolitic limestone, and the marly sediments that partially cover the limestone and fill the many cavities or solution holes (USDA, 1996). These soils are adjacent to the Miami Ridge. Slopes are 0 to 2%. Chekika soils are transitional between well drained soils of the Miami Ridge and the poorly drained marl soils of Sawgrass Marsh in the Everglades. The water table is between 30 and 90 cm and is always within the limestone bedrock. The soils are cultivated to tomatoes, beans, corn, malanga, and limes as the principal crops. Geographically associated soils include Krome and Perrine.
The Krome (loamy-skeletal, carbonatic, hyperthermic Lithic Udorthents) series consists of very shallow, moderately well drained, moderately permeable soils over limestone. Soils formed by scarification of outcrops of oolitic limestone, and the loamy residuum that partially covers the limestone and fills the many cavities or solution holes (USDA, 1996). These soils are on broad areas of the Miami Ridge. Slopes are predominantly 0 to 2% but range to 5%. The water table is between 1 and 1.5 m and is also always within the limestone bedrock. The soils are cultivated and tomatoes, beans, avocadoes, and limes are the main crops.
The Perrine marl (coarse-silty, carbonatic, hyperthermic Typic Fluvaquents) series consists of moderately deep, poorly drained, moderately slowly to moderately permeable soils in lowlands along the Atlantic Coast of Peninsular Florida. The soils formed in calcareous silty and loamy sediments of marine or freshwater origin over limestone (USDA, 1996). Slopes are less than 1%. The areas that have been drained, or have controlled water regimes, are used for growing vegetable crops, ornamental shrubs, and trees. Perrine soils are closely associated with Chekika, Krome, and Lauderhill.
Noncarbonatic Soils
The Webster (fine-loamy, mixed, superactive, mesic Typic Endoaquolls) series consists of very deep, poorly drained, moderately permeable soils formed in glacial till or local alluvium derived from till on uplands (USDA, 1996; Konen et al., 2003; Reuter and Bellb, 2003). Slope ranges from 0 to 3%. These soils are largely cultivated and cropped intensively with corn and soybeans. Small grain and hay are other major crops (Iowa State University Cooperative Extension Service, 1998). The soil was sampled from Iowa.
The Lauderhill (euic, hyperthermic Lithic Haplosaprists) series consists of very poorly drained soils that are 0.30 to 0.90 m thick over limestone. Lauderhill soils formed in organic deposits of freshwater marshes (USDA, 1996). Lauderhill soils are in freshwater marshes within peninsular Florida. Lauderhill soils are used for growing sugarcane, corn, sod, and improved pasture.
The five soils were selected for this investigation based on contrasting OC and CaCO3 contents. Organic carbon content was determined by the method of Walkley and Black (1934). The thermogravimetric (TG) method (Nestler et al., 2003) was used to determine CaCO3. Particle size distribution was measured by the pipette method (USDA, 1992). Sterile soils for sorption experiments were prepared by autoclaving once for 1 h at 120°C. Important physical and chemical properties of the soils that can influence atrazine and diuron sorption are listed in Table 1.
Pesticides
Atrazine and diuron were selected as probe pesticides because they are widely used in South Florida on carbonatic soils. We also have a large database on the two herbicides from our previous studies on sorption and transport in noncarbonatic soils. The selected properties of the two herbicides (Table 2) show a close resemblance in their molecular weight, total surface area, and hydrocarbonaceous surface area. As such we expected the herbicides to behave similarly during sorption in soils.
Two-Site Nonequilibrium Batch Sorption Kinetics Model
The two-site nonequilibrium (TSNE) batch sorption kinetics model (see Appendix) was employed in this study to characterize sorption kinetics of atrazine and diuron in the five soils. Among the nonequilibrium processes influencing sorption is the possible division of sorption sites into two types. The two-site sorption concept presumes that sorption or exchange sites in soils can be classified into two fractions: Type 1, where sorption is assumed to be instantaneous, and Type 2, where sorption is considered time-dependent. In Fig. 9, the solid line represents reversible reaction (equilibrium sites) and the adsorbed fraction controlled by first-order sorption kinetics. The parameters and the governing equations for the TSNE model used to calculate batch sorption kinetics parameters (k2, F, and K*) are presented in detail in the Appendix.
Two equations were used to calculate the relative solution concentration (C/Cin) and the apparent sorption coefficient (K*) as a function of time:
 | [1] |
The three parameters (k2, F, and K) are determined from the data using nonlinear curve fitting procedures on the C/Cin vs. time data (Eq. [1]). Thus a series of experiments with different Cin are normalized to one line by plotting C/Cin vs. time:
 | [2] |
The apparent sorption coefficient (K*) was calculated as a function of time using parameters C/Cin (from Eq. [1]), and V and M that were measured before the sorption kinetics experiment was initiated.
Procedures for Sorption Experiments
Sorption isotherms were measured with the batch equilibration method (Nkedi-Kizza et al., 1985) at different times between 0.5 to 50 h. A solution of 10 mL of atrazine or diuron of initial concentration of 4, 8, 14, and 18 mg L1 in 0.01 M CaCl2, as a supporting electrolyte, was added to 10 g of autoclaved soils (in triplicates) in Teflon-lined centrifuge tubes. For the Lauderhill soil 0.5 g of soil were used. The tubes were shaken on a reciprocating shaker for different times at ambient temperature (22 ± 1°C). The tubes were then centrifuged at 6720 x g for 15 min and atrazine and diuron in clear supernatant solutions were analyzed by reverse phase high performance liquid chromatography (HPLC) with UV detection, using a Dionex (Sunnyvale, CA) 500 IC/HPLC system and Dionex-Peaknet 5.1 computing integrator. Atrazine and diuron were analyzed at 220 nm with absorbency set at 0.1 AU for all samples. The flow rate was 1 mL per minute. A C-18 column was used with a mobile phase of methanol, acetonitrile, and water (25:40:35). The amount of solute retained in the soil solution normalized to the initial solution concentration (C/Cin) was plotted versus time and the kinetic parameters were calculated using the TSNE model. The sorption coefficient (K) for each chemicalsoil combination was calculated from an isotherm obtained by plotting S vs. C for each of the four initial solution concentrations measured at or after 24 h. All isotherms were described by the linear Freundlich isotherm (S = KC). The obtained K values were used as initial estimates in the TSNE model.
 |
RESULTS AND DISCUSSION
|
|---|
Equilibrium Sorption Coefficient (K)
The data in Table 1 show that the soils of Miami-Dade County (Krome, Chekika, and Perrine) are highly carbonatic with over 600 g kg1 CaCO3 content. Perrine has the highest CaCO3 content (880 g kg1) among the three soils. Considering that both atrazine and diuron are neutral at the pH of the experiments, any observed differences in sorption behavior can be attributed to differences in OC and CaCO3 content in these soils. All sorption equilibrium isotherms were linear and equilibrium was attained in about 24 h. An example of the isotherms is presented in Fig. 1 . The sorption coefficient (K) values calculated from the isotherms and those calculated using the TSNE model are presented in Table 3. For the two noncarbonatic soils, Webster and Lauderhill, the Koc values calculated by normalizing the equilibrium K with soil OC content were close to literature reported values for atrazine (Koc = 100) and for diuron (Koc = 400). However, in carbonatic soils atrazine and diuron had Koc values of about one-third of literature Koc values for noncarbonatic soils. For example comparing Chekika soil (carbonatic) and Webster soil (noncarbonatic) that have similar OC contents, we observe that sorption coefficients were more than twice in Webster over Chekika for both atrazine and diuron (Table 3). This seems to suggest that the nature and makeup of OC somehow affects the sorption efficiency of OC in carbonatic soils. Although Perrine soil has nearly three times the clay content (Table 1) of the Krome and Chekika soils, clay content did not seem to have an effect on the sorption coefficient of atrazine and diuron at pH 7.
Sorption Kinetics Parameters
Time-dependent sorption of atrazine and diuron by different soils with initial concentration of 8 mg L1 are shown in Fig. 2
5.
To avoid crowding of the figures with many data points we use the initial concentration of 8 mg L1 to discuss the data because data analysis in Table 4 indicated that initial solution concentration did not affect the magnitude of the parameters (k2, F, and K) calculated using Eq. [1]. However, the kinetic parameters presented in Table 3 were obtained from all four initial solution concentrations for each soilchemical combination.
View this table:
[in this window]
[in a new window]
|
Table 4. Effect of atrazine initial solution concentration on sorption kinetic parameters in Lauderhill soil, with standard error in parentheses.
|
|
Data in Fig. 2 to 5 show that diuron sorption is more time-dependent in all soils than atrazine. Atrazine in most soils approaches equilibrium in less than 10 h, but diuron approaches equilibrium in about 15 h. This trend is further supported by the smaller k2 values for diuron compared to atrazine (Table 3). The TSNE model described all sorption kinetics data well except for atrazine in Perrine soil. The sorption data for atrazine in Perrine soil show that the TSNE model is not valid since F = 0 (Table 3 and Fig. 4). The atrazine data in Fig. 4, in fact, can be described by the one-site nonequilibrium (OSNE) model. The OSNE model is a special case of the TSNE model, where all sorption sites are assumed to be time-dependent (Selim et al., 1976; Van Genuchten, 1981). Nzengung et al. (1997) used the OSNE model to describe sorption kinetics of several organic chemicals by organoclays. The TSNE model, however, described the data well for diuron sorption in Perrine soil (Fig. 4 and Table 3). Note that Perrine and Krome soils have similar OC content but the TSNE model described atrazine sorption kinetics data well in Krome soil. Perrine is the only marl soil among the carbonatic soils used in this study. The Lauderhill soil, which has the highest OC content, exhibited the slowest sorption kinetics (smallest k2 values) for both pesticides. In Florida, Lauderhill soil is found in association with Chekika, Krome, and Perrine.
k2K Relationship
For soils that do not contain expanding clay minerals, sorption kinetics of neutral organic chemicals have been attributed to soil organic matter. Since organic matter is the predominant component of soils that adsorbs neutral organic chemicals, intra-organic matter diffusion has been proposed as the rate limiting mechanism for sorption kinetics (Khan, 1973; Brusseau et al., 1991). The fist-order desorption rate coefficient (k2) is a lumped parameter which combines several processes that control mass transfer of sorbates into the organic matter matrix. These can include molecular size, molecular weight, hydrocarbonaceous surface area, total molecular surface area, diffusion path length, and makeup of organic matter that determines resistance to diffusion (Brusseau et al., 1991). The diffusion path length should increase with increase in organic matter content of the soil, which in turn should decrease the value of k2. The sorption coefficient K should increase as the organic matter content of the soil increases.
The sorption parameters k2 and K have been reported to be inversely related for sorption of neutral organic chemicals in soils and sediments (Karickhoff, 1980; Karickhoff and Morris, 1985; Rao et al., 1979; Brusseau et al., 1991). The data in Table 2 show that atrazine and diuron have many similar values of parameters that can affect the magnitude of k2. However, the correlation lines between k2 and K show that atrazine and diuron behaved differently in these five soils (Fig. 6
). Atrazine and diuron seem to behave similarly in only the Lauderhill muck. For Webster and Chekika soils diuron and atrazine data were described by different correlations. In addition, for two carbonatic soils (Perrine and Krome) the k2K relationship for either atrazine or diuron does not describe the sorption kinetics data. In Florida Chekika soils are transitional between well drained soils (Krome) and the poorly drained marl soils (Perrine). These data seem to imply that the makeup of organic matter in carbonatic soils is different from organic matter in noncarbonatic soils. This trend in sorption kinetics data can be used to further explain the low Koc values that were calculated for atrazine and diuron in carbonatic soils.

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 6. Relationship between k2 and K. The top line is for atrazine and the bottom line is for diuron. Two soils (Krome and Perrine) were excluded while calculating the correlation lines for atrazine and diuron.
|
|
Apparent Sorption Coefficient (K*)
Data in Fig. 7
and 8
show examples of the apparent sorption coefficient (K*) as a function of time for atrazine and diuron in Lauderhill and Chekika soils. The sorption coefficients asymptotically approached equilibrium K values in about 24 h. The TSNE model (Eq. [2]) describes the apparent sorption coefficient data reasonably well. Chen and Wagenet (1997) observed similar increase of atrazine adsorption in Niagara silt loam from 1 to 237 h. Their data indicated that adsorption did not reach equilibrium even after 237 h of incubation. In this study, all soilchemical combinations during the kinetics experiments reached equilibrium in 24 h or less. These data clearly show that sorption isotherms can be generated at any time, but the correct isotherm is the one that is measured at equilibrium. Sorption coefficient values calculated from isotherms that are not measured at equilibrium will underestimate sorption.
 |
CONCLUSIONS
|
|---|
All three carbonatic soils (Krome, Chekika, and Perrine) adsorbed diuron and atrazine much less than the noncarbonatic soil (Webster) and the organic soil (Lauderhill) as evident from their Koc values. This might be due to differences in the composition of organic matter in these soils. The TSNE sorption model described the sorption kinetics data well for all soilchemical combinations except for Perrineatrazine that seems to be described by the OSNE model. Perrine is the only soil characterized as a marl soil among the three carbonatic soils and contains the largest amount of CaCO3. Diuron exhibited slower sorption kinetics (smaller k2 values) than atrazine in all soils. The difference in sorption kinetics between atrazine and diuron increased with increase in soil OC content. The batch kinetics experiments were used to differentiate the behavior of two similar organic pesticides in five soils. The data were also used to point to the possible differences in organic matter makeup between carbonatic and noncarbonatic soils.
The Koc values for atrazine and diuron obtained in this study for carbonatic soils were a third of the values reported literature for noncarbonatic soils. This implies that use of literature Koc values for predicting fate and transport of similar pesticides in carbonatic soils would result in underestimation of pesticide travel times to ground water.
As a whole, the batch sorption kinetics experiments proved to be a useful tool in distinguishing the behavior of two organic pesticides in carbonatic and noncarbonatic soils and in identifying a potential difference in organic matter makeup between carbonatic and noncarbonatic soils. Further research is, therefore, warranted to investigate the origin and nature of organic matter of these carbonatic soils.
 |
APPENDIX
|
|---|
Figure 9
illustrates the two-site nonequilibrium (TSNE) batch sorption kinetics model. For the TSNE model linear sorption is assumed. All the solute is initially in solution.
Terms include:
Term
|
Definition
|
Unit
|
| A |
total mass of solute added to soil, always constant = VCin |
mg |
| C |
solute concentration in the water |
mg L1 |
| Cin |
initial, added solution concentration |
mg L1 |
| C0 |
solution concentration after partitioning to the instantaneous site at t = 0 |
mg L1 |
C |
solution concentration at equilibrium |
mg L1 |
| S |
total sorbed solute concentration = S1 + S2 |
mg kg1 |
| S1 |
sorbed solute concentration on the soil of the equilibrium region |
mg kg1 |
| S1 |
F KC at all times |
mg kg1 |
| S2 |
sorbed solute concentration on the soil of the kinetic region |
mg kg1 |
| S20 |
sorbed solute concentration on the soil of the kinetic region at (t = 0) |
mg kg1 |
| t |
time |
h |
| F |
fraction of the total sorption in the Type 1 region when the system is in equilibrium |
|
| K |
linear sorption partition coefficient at equilibrium |
L kg1 |
| K* |
apparent linear sorption partition coefficient before equilibrium |
L kg1 |
| V |
volume of water |
L |
| M |
mass of soil |
kg |
| k1 |
kinetic rate coefficient for sorption to the Type 2 sites |
1/h |
| k2 |
kinetic rate coefficient for desorption from the Type 2 sites |
1/h |
| R |
retardation factor |
|
ß
|
fraction of retardation for the instantaneous region
|
|
|
Define:
The Solution for C versus Time
First at equilibrium (time =
) total sorption is given by:
 | [A1] |
 | [A2] |
Solute mass transfer from the Type 1 sites to the Type 2 sites is given by:
 | [A3] |
which implies at time =
:
 | [A4] |
The total solute mass, A (mg), added to the soil is constant and is given by:
 | [A5] |
 | [A6] |
At any given time, S2 is given in terms of C by:
 | [A7] |
And at equilibrium, time =
, where S2 = (1 F)KC, the equilibrium solution concentration is given by:
 | [A8] |
The basic equations for this system are:
 | [A9] |
 | [A10] |
By substituting for k1 from Eq. [A4] and for S2 from Eq. [A7], Eq. [A9] becomes:
 | [A11] |
The solution to Eq. [A11] in terms of C is:
 | [A12] |
where ß is given by:
 | [A13] |
By using Eq. [A8] for C
, Eq. [A12] can be rearranged to:
 | [A14] |
Sorption Experiment
Water with solute at concentration Cin is added to soil containing no sorbed solute. Thus solute mass, A = VCin. The solute partitions instantly to sites of Type 1 such that:
 | [A15] |
 | [A16] |
 | [A17] |
 | [A18] |
With these substitutions into Eq. [A14] the solution becomes:
 | [A19] |
The three parameters (F, k2, and K) are determined from the data using nonlinear curve fitting procedures on the C/Cin vs. time data. Thus a series of experiments with different Cin is normalized to one line by plotting C/Cin vs. time.
Apparent Sorption Coefficient (K*)
The term S vs. time is given by:
 | [A20] |
 | [A21] |
Using C/Cin obtained in Eq. [A19], the apparent sorption coefficient (K*) in Eq. [A21] can be calculated as a function of time in Eq. [A22]:
 | [A22] |
 |
REFERENCES
|
|---|
- Aharoni, C., and D.L. Sparks. 1991. Rates of soil chemical processesA theoretical treatment. p. 281302. In D.L. Sparks and D.L. Suarez (ed.) Rates of soil chemical processes. Spec. Publ. 27. SSSA, Madison, WI.
- Ben-Hur, M., J. Letey, W.J. Farmer, C.F. Williams, and S.D. Nelson. 2003. Soluble and solid organic matter effects on atrazine adsorption in cultivated soils. Soil Sci. Soc. Am. J. 67:11401146.[Abstract/Free Full Text]
- Brusseau, M.L., R.E. Jessup, and P.S. Rao. 1991. Nonequilibrium sorption of organic chemicals: Elucidation of rate-limiting processes. Environ. Sci. Technol. 25:134142.[CrossRef]
- Cameron, D.A., and A. Klute. 1977. Convective-dispersive solute transport with a combined equilibrium and kinetic adsorption model. Water Resour. Res. 19:718724.
- Chen, W., and R.J. Wagenet. 1997. Description of atrazine in soil with heterogeneous nonequilibrium sorption. Soil Sci. Soc. Am. J. 61:360371.[Abstract/Free Full Text]
- Downing, H.F., M.E. Delorenzo, M.H. Fulton, G.I. Scott, C.J. Madden, and J.R. Kucklick. 2004. Effects of agricultural pesticides atrazine, chlorothalonil, and endosulfan on South Florida microbial assemblages. Ecotoxicology 13:245260.[Medline]
- Field, J.A., R.L. Reed, T.E. Sawyer, S.M. Griffith, and J.R. Wigington. 2003. Diuron occurrence and distribution in soil and surface and ground water associated with grass seed production. J. Environ. Qual. 32:171179.[Abstract/Free Full Text]
- Harter, R.D. 1991. Kinetic of sorption/desorption processes in soil. p. 135149. In D.L. Sparks and D.L. Suarez (ed.) Rates of soil chemical processes. Spec. Publ. 27. SSSA, Madison, WI.
- Iowa State University Cooperative Extension Service. 1998. Iowa crop performance test. Bull. Pm-660-2-98. Iowa State Univ., Ames.
- Karickhoff, S.W. 1980. Sorption kinetics of hydrophobic pollutants in natural sediments. p. 193205. In R.A. Baker (ed.) Contaminants and sediments. Vol. 2. Ann Arbor Sci. Publ., Ann Arbor, MI.
- Karickhoff, S.W., and K.R. Morris. 1985. Sorption dynamics of hydrophobic pollutants in sediment suspensions. Environ. Toxicol. Chem. 4:469474.
- Khan, S.U. 1973. Equilibrium and kinetic studies of the adsorption of 2,4-D and picloram on humic acid. Can. J. Soil Sci. 53:429434.
- Konen, M.E., C.L. Burrasb, and J.A. Sandorb. 2003. Organic carbon, texture, and quantitative color measurement relationships for cultivated soils in north central Iowa. Soil Sci. Soc. Am. J. 67:18231830.[Abstract/Free Full Text]
- Ma, L., and H.M. Selim. 1997. Evaluation of nonequilibrium models for predicting atrazine transport in soils. Soil Sci. Soc. Am. J. 61:12991307.[Abstract/Free Full Text]
- Madhun, Y.A., V.H. Freed, J.L. Young, and S.C. Fang. 1986. Sorption of bromacilm chlortoluron, and diuron by soils. Soil Sci. Soc. Am. J. 50:14671471.[Abstract/Free Full Text]
- Miles, C.J., and R.J. Pfeuffer. 1997. Pesticides in canals of South Florida. Arch. Environ. Contam. Toxicol. 32:337345.[Medline]
- Nestler, K., D. Dietrich, K. Witke, R. Rößler, and G. Marx. 2003. Thermogravimetric and Raman spectroscopic investigations on different coals in comparison to dispersed anthracite found in permineralized tree fern Psaronius sp. J. Mol. Struct. 661662:357362.
- Nielsen, D.R., M.Th. van Genuchten, and J.W. Biggar. 1986. Water flow and solute transport processes in the unsaturated zone. Water Resour. Res. 22(Supplement):89108.
- Nkedi-Kizza, P., J.W. Biggar, H.M. Selim, P.J. Wierenga, M.Th. van Genuchten, J.M. Davidson, and D.R. Nielsen. 1984. On the equivalence of two conceptual models for describing ion exchange during transport through an aggregated Oxisol. Water Resour. Res. 20:11231130.
- Nkedi-Kizza, P., P.S.C. Rao, and A.G. Hornsby. 1985. Influence of organic cosolvents on sorption of hydrophonic organic chemicals by soils. Environ. Sci. Technol. 19:975979.[CrossRef]
- Nzengung, V.A., P. Nkedi-Kizza, E.A. Voudrias, and R. Jessup. 1997. Organic cosolvent effects on sorption kinetics of hydrophobic organic chemicals by organoclays. Environ. Sci. Technol. 31:14701475.[CrossRef]
- Pait, A.S., A.E. De Souza, and D.R.G. Farrow. 1992. Agricultural pesticide use in coastal areas: A national summary. NOAA/NOS/ORCA, Silver Spring, MD.
- Pedit, J.A., and C.T. Miller. 1994. Heterogeneous Sorption processes in subsurface systems: Development and application of models to batch systems. Center of Multiphase Research News, Vol. 1, p. 24.
- Rao, P.S.C., J.M. Davidson, R.E. Jessup, and H.M. Selim. 1979. Evaluation of conceptual models for describing nonequilibrium adsorption-desorption of pesticides during steady-flow in soils. Soil Sci. Soc. Am. J. 43:2228.[Abstract/Free Full Text]
- Rao, P.S.C., and A.G. Hornsby. 1983. Behavior of pesticides in soils and water. Soil Sci. Fact Sheet SL40. Univ. of Florida, Gainesville.
- Reuter, R.J., and J.C. Bellb. 2003. Hillslope hydrology and soil morphology for a wetland basin in south-central Minnesota. Soil Sci. Soc. Am. J. 67:365372.[Abstract/Free Full Text]
- Selim, H.M. 2004. Modeling kinetic retention of atrazine and metribuzin in soil. Soil Sci. 169:2534.[CrossRef]
- Selim, H.M., J.M. Davidson, and R.S. Mansell. 1976. Evaluation of a two-site adsorption-desorption model for describing solute transport in soils. p. 444448. In Proc. Summer Computer Simulation Conf., Washington, DC. 1214 July 1976. Simulation Councils, La Jolla, CA.
- USDA. 1992. Soil survey laboratory methods manual. USDA, Washington, DC.
- USDA. 1996. Florida online soil survey manuscriptsDade County. Available at http://soils.usda.gov/survey/online_surveys/florida/dade/index.html (verified 20 Sept. 2005). USDA, Washington, DC.
- USEPA. 1987. Diuron health advisory. USEPA Office of Drinking Water, Washington, DC.
- USEPA. 1998. Safe drinking water information system. EPA 816-F-98-006. USEPA Office of Drinking Water, Washington, DC.
- Van Genuchten, M.Th. 1981. Non-equilibrium transport parameters from miscible displacement experiments. Res. Rep. 119. U.S. Salinity Lab., Riverside, CA.
- Walkley, A., and I.A. Black. 1934. An examination of the digestion method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 37:2938.
- Wauchope, R.D., T.M. Buttler, A.G. Hornsby, P.W. Augustijn-Beckers, and J.P. Burt. 1992. The SCS/ARS/CES pesticide properties database for environmental decision-making. Rev. Environ. Contam. Toxicol. 123:1155.[ISI][Medline]
- Wu, S., and P.M. Gschwend. 1986. Sorption kinetics of hydrophobic organic compounds to natural sediments and soils. Environ. Sci. Technol. 20:717725.[CrossRef][ISI]
This article has been cited by other articles:

|
 |

|
 |
 
A. Mermoud, J. M. F. Martins, D. Zhang, and A. C. Favre
Small-Scale Spatial Variability of Atrazine and Dinoseb Adsorption Parameters in an Alluvial Soil
J. Environ. Qual.,
August 8, 2008;
37(5):
1929 - 1936.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. L. Milfont, J. M. F. Martins, A. C. D. Antonino, E. R. Gouveia, A. M. Netto, V. Guine, H. Mas, and M. B. G. dos Santos Freire
Reactivity of the Plant Growth Regulator Paclobutrazol (Cultar) with Two Tropical Soils of the Northeast Semiarid Region of Brazil
J. Environ. Qual.,
January 4, 2008;
37(1):
90 - 97.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
T. L. Potter, D. D. Bosch, H. Joo, B. Schaffer, and R. Munoz-Carpena
Summer Cover Crops Reduce Atrazine Leaching to Shallow Groundwater in Southern Florida
J. Environ. Qual.,
July 17, 2007;
36(5):
1301 - 1309.
[Abstract]
[Full Text]
[PDF]
|
 |
|