|
|
||||||||
a Biological and Environ. Eng., Riley-Robb Hall, Cornell Univ., Ithaca, NY 14853
b Dep. of Entomology, Fernald Hall, Univ. of Massachusetts, Amherst, MA 01003-2410
* Corresponding author (dah13{at}cornell.edu)
Received for publication May 11, 2001.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Field measurements of pesticide volatilization are time-consuming, expensive, and are impractical for the full range of chemicals, weather, and site conditions encountered in practice. Environmental and health assessments of pesticide volatilization will likely be based on mathematical modeling. Available volatilization models are typically based on equilibrium partitioning of solid, liquid, and gaseous chemical phases in the soil environment (Carsel et al., 1998; Jury et al., 1990; Wang et al., 2000; Yates, 1993). However, soil may be a minor source of gaseous losses in turf systems due to the dense protective layer of foliage and thatch. Direct vaporization of the chemical from vegetation, which is neglected in most models, is a more likely source. It is difficult to extend the equilibrium modeling approach to plant material due to the uncertainties of solidliquidgas partitioning. At least one model, PRZM-3 (Carsel et al., 1998), includes volatilization from foliage, but only as a constant daily sink.
Nonpartitioning approaches to volatilization may be more feasible. Hill and Schaalje (1985) developed a two-compartment model of pesticide dissipation in soil. Losses from compartment one are surface losses, such as volatilization, and redistribution to compartment two, where biochemical degradation takes place. These three transformations are assumed first-order, with empirical rate constants. Weed et al. (1999) created a modified two-compartment model by rationalizing the rate constants for surface losses and redistribution. Surface losses were treated explicitly as volatilization, with rates related to water evaporation, and redistribution was assumed to be a physical movement of the chemical into the soil due to wash off by rain.
Determination of chemical volatilization from water evaporation rates makes sense because of the similarities between the two processes (Taylor and Spencer, 1990). Factors such as air temperature, solar radiation, and wind movement could be expected to have comparable effects on vaporization of both water and chemicals. Moreover, models for estimating evaporation are readily available and evaporation values from such models may be converted to pesticide vaporization by scaling factors that reflect the differences in chemical properties of the water and pesticide. Weed et al. (1999) based their scaling on the ratio of vapor densities, and fitted the resulting relationship to alachlor [2-chloro-N-(2,6-diethylphenyl)-N-(methoxymethyl)acetamide] volatilization from soil.
This paper describes an evaporation-based approach for estimating pesticide volatilization from turfgrass. Since volatilization is from vegetation rather than soil, water losses are based on an evapotranspiration model. Scaling factors are determined from vapor pressures and heats of vaporization. The model was tested by comparison of predictions with field data for eight pesticides.
| METHODS |
|---|
|
|
|---|
![]() | [1] |
![]() | [2] |
In this equation ETt = evapotranspiration during hour t (mm); psct, pswt = saturated vapor pressures of the chemical and water, respectively, during hour t (kPa); and
ct,
wt = latent heats of vaporization of the chemical and water, respectively, during hour t (J g-1).
The first ratio in Eq. [2] adjusts the evapotranspiration value for the differences in chemical and water vapor pressures. The second ratio reflects the differences in energy requirements for vaporizing water and the pesticide. The overall pesticide mass balance must account for other losses, such as bio- or photochemical degradation. Assuming these losses are first-order, or exponential, we have
![]() | [3] |
= degradation rate of the pesticide on vegetation surfaces (h-1). Equation [3] neglects other forms of chemical loss including runoff, leaching, and removal of clippings. As a result, the model may overestimate the pesticide available for leaching.
Evapotranspiration
The Penman equation, as described by Jensen et al. (1990), provides an estimate of evapotranspiration for well-watered short grass:

[4]
where,
t = slope of saturation vapor pressure curve during hour t (kPa °C-1),
t = psychrometric constant during hour t (kPa °C-1), Rnt = net radiant energy available at surface during hour t (kJ m-2 h-1), Gt = net sensible heat flux from the surface to soil during hour t (kJ m-2 h-1), ut = mean wind velocity during hour t (m s-1), and pwt = actual water vapor pressure during hour t (kPa).
Parameters in Eq. [4] are given in the following equations, as described by Jensen et al. (1990):
![]() | [5] |
![]()
![]() | [7] |
![]() | [8] |
![]() | [9] |
The net radiant energy term in Eq. [4] is given by incoming solar radiation minus reflected radiation and net thermal radiation. Only the first of these is generally known, and net energy is usually estimated from regression equations. Jensen et al. (1990) provide a number of these equations determined for various covers and locations, including the following that was obtained for grass in Minnesota:
![]() | [10] |
Pesticide Relationships
Vapor pressures of volatile pesticides vary significantly with temperatures, as determined from Grain (1982):

where psc0 = vapor pressure (kPa) at absolute temperature Ta0 (K),
c0 = latent heat of vaporization (J/g) at Ta0, M = molecular weight of the chemical, R = gas constant (8.32 J/mol-K),
Zb = compressibility factor at boiling point (dimensionless), Tat = absolute temperature during hour t (K), and m is a constant.
The compressibility factor is assumed to be 0.97 in the examples given in Grain (1982). Heat of vaporization is approximated from the ideal gas law, and for the above dimensions,
![]() | [12] |
ct, is also estimated from Eq. [12], with Tat substituted for Ta0. The constant m in Eq. [11] is 0.19 for liquid chemicals, and may be either 0.36, 0.8, or 1.19 for solids, depending on boiling point. However, for the range of field temperatures of interest to pesticide volatilization (280310 K), Eq. [11] is relatively insensitive to m, and a value of m = 0.8 can be used for solid chemicals with minimal loss of accuracy.
Field Studies
The volatilization model was tested using data from field turf experiments conducted at the University of Massachusetts Turfgrass Research Center in South Deerfield, MA. Experimental design and sampling methods were as described by Murphy et al. (1996a)(1996b). The 0.2-ha plots had well-established creeping bentgrass maintained at 13 mm. Thatch thickness ranged from 10 to 15 mm. The soil was a Hadley silt loam (coarse-silty, mixed, superactive, nonacid, mesic Typic Udifluvents). Irrigation was applied as necessary to prevent drought stress. Testing data included the measured concentrations of volatile residues following application of eight pesticides in 11 experiments conducted in the growing seasons of 1995, 1996, and 1997, as described in Table 1. Plots were mown on Monday, Wednesday, and Friday and applications occurred either immediately after mowing or on Tuesday or Thursday. Ethoprop and isofenphos were the most extensively studied chemicals, with applications in 7 and 6, respectively, of the 11 experiments. The remaining pesticides (bendiocarb, carbaryl, chlorpyrifos, diazinon, isazofos, and trichlorfon) were included in four of the experiments.
|
For the 7-d experiments, an initial sample was taken for 1 h after application, and two or three additional sampling intervals of 2 to 4 h covered the remainder of the first day. On Days 2 and 3, three sampling intervals of 4 h each were carried out on Days 2 and 3, and one sample each was taken on Days 5 and 7. The shorter, 2- and 3-d experiments were similar, but omitted sampling on the latter days. Measured concentrations generally fell to very low levels after the second day following application. The TheoreticalProfileShape method (Murphy et al., 1996a, 1996b; Wilson et al., 1982) was used to estimate volatilization mass flux.
All pesticide applications were made as sprays and immediately were watered in by 6.5 mm postapplication irrigation (Murphy et al., 1996a, 1996b). Applications were made to the 0.2-ha circular turfgrass plots in a manner and at a rate of application commonly used for the selected pesticides on golf courses and according to the manufacturer's instructions. Formulated pesticides were applied in 51.1 L (13.5 gal) of water to each plot. Tank mixtures consisted of 333 mL (11.25 fl oz) of diazinon (diazinon, Giba Geigy) [O,O-diethyl O-(2-isopropyl-6-methyl-4-pyrimidinyl) phosphorothioate], 555 mL (18.75 fl oz) of Dursban Pro (chlorpyrifos, Dow AgroSciences) [O,O-diethyl O-(3,5,6-trichloro-2-pyridyl) phosphorothioate], 333 mL (11.25 fl oz) of Mocap (ethoprop, Rhone-Poulenc) (O-ethyl S,S-dipropyl phosphorodithioate), 333 mL (11.25 fl oz) of Oftonol 2 (isofenphos, Bayer) [1-methylethyl 2-((ethoxy((1-methylethyl)amino)phosphinothioyl)oxy) benzoate], 416 mL (14.1 fl oz) of Proxol (trichlorfon, Nor-Am) [dimethyl (2,2,2-trichloro-1-hydroxyethyl)phosphonate], 167 mL (5.6 fl oz) of Sevin (carbaryl, Lesco Inc.) (1-naphthyl-N-methylcarbamate), 167 mL (5.6 fl oz) of Triumph (isazofos, Ciba Geigy) [O,O-diethyl O-(5-chloro-1-(1-methylethyl)-1H-1,2,4-triazol-3-yl) phosphorothioate], and 213 g (0.47 lb) of Turcam (bendiocarb, NOR-AM) (2,2-dimethyl-1,3-benzoldioxol-4-yl methylcarbamate) per 57 L (15 gal) of water.
Sample Collection and Residue Analyses
Volatile pesticide residues were collected onto XAD-4 resin during high-volume air sampling and extracted with acetone for instrumental analysis (Murphy et al., 1996a, 1996b). Extraction efficiencies ranged from 89 to 99% recovery. Chlorpyrifos, diazinon, ethoprop, isofenphos, and trichlorfon were analyzed simultaneously with isazofos using the method described by Murphy et al. (1996a). Bendiocarb and carbaryl were analyzed using the postcolumn derivatization/HPLC method (Environ. Monit. Systems Lab., 1989). Detection limits for all analytes were 0.02 µg/m3.
Model Testing
Properties of the eight pesticides are listed in Table 2. Equations [1] and [2] were used to predict the volatilization mass flux (Vt) for each pesticide during each sampling interval. Remaining pesticide available for the next interval was determined from Eq. [3]. First-order degradation rates (
) are not generally available for turf vegetation, so the soil half-lives for aerobic biodegradation given in Table 2 were used to determine
. There is limited evidence that soil and turf decay rates are comparable (Cisar and Snyder, 1996), but other studies have found substantial differences (Carroll et al., 2000; Horst et al., 1996). Pesticide saturated vapor pressures and heats of vaporization were calculated by Eq. [11] and [12] using vapor pressure data and molecular weights from Table 2.
|
It was necessary to develop methods to estimate the volatilization constant (k) in Eq. [1]. If vaporization of the chemical is completely analogous to water evaporation, then it would seem that k = 1. However, pesticide volatilization will be inhibited by adsorption by plant material, incorporation into plant tissue, and drying of sprayed chemical. These processes are not modeled explicitly; rather, their effects are approximated by the volatilization constant.
Since each pesticide in the field studies was included in several experiments, a unique volatilization constant could be determined for each pesticide by model calibration. However, this would limit subsequent applications of the model to the eight chemicals measured in the field studies. A more general approach is to group the pesticides into chemicals with similar characteristics, and determine a single rate constant for each group. These constants would presumably be applicable to other chemicals that could be placed in these groups.
Clark et al. (2000) classified pesticides according to health hazards associated with volatilization. The most hazardous were chemicals with vapor pressures exceeding 10-5 mm Hg (<1.3 x 10-6 kPa). This vapor pressure was used to divide the eight chemicals in this study into the two groups shown in Table 2. The exception was chlorpyrifos, which at 2.7 x 10-6 kPa would just exceed the group 2 boundary. However, chlorpyrifos is very strongly adsorbed to organic matter, as indicated by Koc = 9930 cm3 g-1, and this property could presumably have effects similar to a low vapor pressure. The volatilization constant k was estimated by selecting one pesticide from each group and adjusting the constant to produce a match of modeled and measured volatilization fluxes. The remaining three pesticides were used to test the model results using the calibrated constant.
Total measured and modeled volatile pesticide losses were compared for each experiment. These totals are the sum of the measured or modeled losses for each sampling interval in the experiment. Since these intervals were not continuous, and regular intervals were sometimes omitted due to equipment or human failures, actual volatile losses may have been larger than the measured values.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
Other Losses
The model's tendency to overestimate volatilization may be due to its neglect of pesticide losses due to clipping removal, runoff, and leaching. None of these losses were measured directly, but their approximate magnitudes can estimated. Assuming that the postapplication watering by 6.5 mm of irrigation distributed the chemical uniformly through the foliage and thatch, clippings losses are not likely to be significant because mowing would remove no more than a small fraction of the total vegetation. For example, Cisar and Snyder (1996) found that total removals of sprayed pesticides over several months in grass clippings were generally <1% of applications. However, this conclusion may not be appropriate for granular pesticides. Cisar and Snyder (1996) reported a total clipping removal of nearly 8% of applied granular chlorpyrifos over a 3-mo period.
Estimates of pesticide losses from runoff and leaching are complicated by lack of water balance data. Irrigation applications and rainfall were not recorded and neither runoff nor percolation from the plots were measured. However, runoff would be a relatively rare occurrence from these plots. The Hadley silt loam belongs to soil hydrologic group B, and the associated runoff curve number for heavily thatched short grass is only 55 (Haith, 2001). Leaching of the pesticide from the vegetation to the soil is more likely, and can be estimated from Haith (2001) as:
![]() | [13] |
Equation [13] was used to estimate the pesticide leaching from the 6.5 mm of irrigation immediately following application (Rt = 6.5). Runoff was assumed to be zero, Koc values were taken from Table 2, and OC was determined from typical values given in Haith (2001) (1120 and 109 kg ha-1 per mm thickness of thatch and grass, respectively, or OC = 15420 kg/ha). The results of these computations are given in Table 5. It is apparent that leaching losses for at least three of the pesticides (ethoprop, isazofos, trichlorfon) could be significant, particularly if additional water inputs occurred early in any of the experimental periods.
|
| CONCLUSIONS |
|---|
|
|
|---|
Testing results indicated that the model is a relatively conservative approach for predicting pesticide volatilization. Predicted mean losses exceeding observations by 20%, and the model explained 67% of the observed variation in volatilization fluxes. The model was most accurate for those chemicals that exhibited the largest volatilization losses (ethoprop, diazinon, and chlorpyrifos). Conversely, the model underpredicted losses for two of the three chemical that had relatively small losses (carbaryl and isofenphos).
The determination of volatilization constant based on the two pesticide groups was largely successful, since model predictions approximated observations for all pesticides in a group when the group rate constant was used. It should be reasonable to apply the model, without further calibration, to additional pesticides in these two groups. However, placement of chemicals in the appropriate group may not always be obvious. Based on its 25°C saturated vapor pressure (2.7 x 10-6 kPa), chlorpyrifos would appear to fit into Group 1. However, if the calibrated rate constant for Group 1 (k = 130 mm-1) had been used for modeling chlorpyrifos, estimated volatilization would have been only 30% of observed values. Chlorpyrifos was placed in Group 2 based on its high adsorption potential (Koc = 9930), but the appropriate Koc limits for Group 2 are unknown. Alternatively, since the vapor pressure limits proposed in Clark et al. (2000) are relatively arbitrary, the group limits could be changed to say, 3 to 4 x 10-6 kPa, thus placing chlorpyrifos unambiguously in Group 2.
All pesticides in this study were applied as sprays, and the model may not exhibit similar accuracy with granular chemicals, particularly in the absence of postapplication irrigation. In such cases, pesticide may be less available for volatilization, and the model would likely overestimate volatilization losses.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. R. Walden and D. A. Haith Estimating Turf Pesticide Volatilization from Simple Evapotranspiration Models J. Environ. Qual., May 1, 2003; 32(3): 1138 - 1143. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Vadose Zone Journal | Journal of Plant Registrations | ||||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal |