Estimating turf pesticide volatilization from simple evapotranspiration models

A previously developed model by Haith et al. (2002) related pesticide volatilization from turf to evapotranspiration (ET) by scaling factors determined from vapor pressures and heats of vaporization. Although the model provided volatilization estimates that compared well with field measurements, it relied on the Penman ET equation, requiring hourly temperature, wind speed, and solar radiation data, none of which are routinely available at field sites. The current study determined that the volatilization model works equally well with a simpler ET equation requiring only daily temperatures. Three daily temperature-based ET models were evaluated as vehicles for estimating pesticide volatilization from turf: Hamon, Hargreaves-Samani, and a modified Priestley-Taylor. When compared with field volatilization measurements for eight pesticides, volatilization estimates produced from the Hargreaves-Samani model most closely approximated both the field observations and the previous estimates based on the more data-intensive Penman model. Mean estimated volatilization exceeded mean observations by 15% and the coefficient of variation (R2) between estimates and observations was 0.65. The comparable values based on Penman ET were 17% and 0.63, respectively.     

TurfPQ, a pesticide runoff model for turf

Environmental assessments of golf courses and other turf systems must often rely on mathematical modeling. However, in the case of pesticide runoff, successful modeling applications are rare. Available models were developed for agricultural applications and have seen very limited testing for turf. TurfPQ is a pesticide runoff model developed exclusively for turf. The model is based on a curve number calculation for runoff volume and linear partitioning of pesticide into adsorbed and dissolved components during a precipitation or irrigation event. Calibration is optional, so the model can be applied, using default parameter values, to situations where runoff and chemical loss data are unavailable. TurfPQ was tested with default parameter values for 52 pesticide runoff events involving six pesticides measured in plot studies in four states. The model typically produced conservative overpredictions of pesticide runoff, particularly with strongly adsorbed pesticides. Mean predicted pesticide runoff was 2.9% [corrected] of application, compared with an observed mean of 2.1%. TurfPQ captured the dynamics of the pesticide runoff events well with R2 = 0.65 [corrected]. Sensitivity analyses indicated that prediction errors could be reduced by better estimates of adsorption parameters and runoff curve numbers. However, even with default parameters, TurfPQ predictions are at least as accurate as those produced by more complex models.     

Comparison of field-scale herbicide runoff and volatilization losses: an eight-year field investigation

An 8-yr study was conducted to better understand factors influencing year-to-year variability in field-scale herbicide volatilization and surface runoff losses. The 21-ha research site is located at the USDA-ARS Beltsville Agricultural Research Center in Beltsville, MD. Site location, herbicide formulations, and agricultural management practices remained unchanged throughout the duration of the study. Metolachlor [2-chloro--(2-ethyl-6-methylphenyl)--(2-methoxy-1-methylethyl) acetamide] and atrazine [6-chloro--ethyl--(1-methylethyl)-1,3,5-triazine-2,4-diamine] were coapplied as a surface broadcast spray. Herbicide runoff was monitored from a month before application through harvest. A flux gradient technique was used to compute volatilization fluxes for the first 5 d after application using herbicide concentration profiles and turbulent fluxes of heat and water vapor as determined from eddy covariance measurements. Results demonstrated that volatilization losses for these two herbicides were significantly greater than runoff losses ( < 0.007), even though both have relatively low vapor pressures. The largest annual runoff loss for metolachlor never exceeded 2.5%, whereas atrazine runoff never exceeded 3% of that applied. On the other hand, herbicide cumulative volatilization losses after 5 d ranged from about 5 to 63% of that applied for metolachlor and about 2 to 12% of that applied for atrazine. Additionally, daytime herbicide volatilization losses were significantly greater than nighttime vapor losses ( < 0.05). This research confirmed that vapor losses for some commonly used herbicides frequently exceeds runoff losses and herbicide vapor losses on the same site and with the same management practices can vary significantly year to year depending on local environmental conditions.     

Predicting and measuring environmental concentration of pesticides in air after soil application

Pesticides can volatilize into the atmosphere, which affects the air quality. The ability to predict pesticide volatilization is an essential tool for human risk and environmental assessment. Even though there are several mathematical models to assess and predict the fate of pesticides in different compartments of the environment, there is no reliable model to predict volatilization. The objectives of this study were to evaluate pesticide volatilization under agricultural conditions using malathion [ O,O-dimethyl-S-(1,2-dicarbethoxyethyl)-dithiophosphate], ethoprophos (O-ethyl S,S-dipropylphosphorodithioate), and procymidone [N-(3,5-dichlorophenyl)-1,2-dimethylcyclopropane-1,2-dicarboximide] as test compounds and to evaluate the ability of the Pesticide Leaching Model (PELMO) to calculate the predicted environmental concentrations of pesticides in air under field conditions. The volatilization rate of procymidone, malathion, and ethoprophos was determined in a field study during two different periods (December 1998 and September 1999) using the Theoretical Profile Shape (TPS) method. The experiments were performed on bare silty soil in the Bologna province, Italy. Residues in the air were continuously monitored for 2 to 3 wk after the pesticide applications. The amount of pesticide volatilized was 16, 5, and 11% in December 1998 and 41, 23, and 19% in September 1999 for procymidone, malathion, and ethoprophos, respectively. In both these experiments, the PELMO simulations of the concentration of ethoprophos and procymidone were in good agreement with the measured data (factor +/- 1.1 on average). The volatilization of malathion was underestimated by a factor of 30 on average. These results suggest that volatilization described by PELMO may be reliable for volatile substances, but PELMO may underpredict volatilization for less-volatile substances.     

Risk assessment of pesticide runoff from turf

The TurfPQ model was used to simulate the runoff of 15 pesticides commonly applied to creeping bentgrass (Agrostis stolonifera L.) fairways and greens on golf courses in the northeastern USA. Simulations produced 100-yr daily records of water runoff, pesticide runoff, and pesticide concentration in runoff for three locations: Boston, MA, Philadelphia, PA, and Rochester, NY. Results were summarized as annual and monthly means and annual maximum daily loads (AMDLs) corresponding to 10- and 20-yr return periods. Mean annual pesticide runoff loads did not exceed 3% of annual applications for any pesticide or site, and most losses were substantially less than 1% of application. However, annual or monthly mean concentrations of chlorothalonil, iprodione, and PCNB in fairway runoff often exceeded concentrations that result in 50% mortality of the affected species (LC50) for aquatic organisms. Concentrations of azoxystrobin, bensulide, cyfluthrin, and trichlorfon in extreme (1 in 10 yr or 1 in 20 yr) events often approached or exceeded LC50 levels. Concentrations of halofenozide, mancozeb, MCPP, oxadiazon, propiconazole, thiophanate-methyl, triadimefon, and trinexapac-ethyl were well below LC50 levels, and turf runoff of these chemicals does not appear to be hazardous to aquatic life in surface waters.     

Role of Riparian Areas in Atmospheric Pesticide Deposition and Its Potential Effect on Water Quality

Riparian buffers are known to mitigate hydrologic losses of nutrients and other contaminants as they exit agricultural fields. The vegetation of riparian buffers can also trap atmospheric contaminants, and these pollutants can subsequently be delivered via rain to the riparian buffer floor. These processes, however, are poorly understood especially for pesticide residues. Therefore, we conducted a four‐year study examining stemflow and throughfall to a riparian buffer which was adjacent a cultured Zea mays field treated with atrazine and metolachlor. Stemflow is rain contacting the tree canopy traveling down smaller to larger branches and down the tree trunk, whereas throughfall is rain that may or may not contact leaves and branches and reaches the earth. Stemflow concentrations of the herbicides were larger than throughfall concentrations and accounted for 5‐15% of the atrazine and 6‐66% of the metolachlor depositional fluxes under the canopy. Larger depositional fluxes were measured when leaves were more fully emerged and temperatures and humidity were elevated. Rain collected outside the riparian buffer on the field side and on the back side revealed the trees trapped the herbicide residues. Herbicide loading to the riparian buffer stream was found to be linked to tree canopy deposition and subsequent washoff during rain events. These results indicate that in agricultural areas canopy washoff can be an important source of pesticides to surface waters.     

Sensitivity and first-step uncertainty analyses for the preferential flow model MACRO

Sensitivity analyses for the preferential flow model MACRO were carried out using one-at-a-time and Monte Carlo sampling approaches. Four different scenarios were generated by simulating leaching to depth of two hypothetical pesticides in a sandy loam and a more structured clay loam soil. Sensitivity of the model was assessed using the predictions for accumulated water percolated at a 1-m depth and accumulated pesticide losses in percolation. Results for simulated percolation were similar for the two soils. Predictions of water volumes percolated were found to be only marginally affected by changes in input parameters and the most influential parameter was the water content defining the boundary between micropores and macropores in this dual-porosity model. In contrast, predictions of pesticide losses were found to be dependent on the scenarios considered and to be significantly affected by variations in input parameters. In most scenarios, predictions for pesticide losses by MACRO were most influenced by parameters related to sorption and degradation. Under specific circumstances, pesticide losses can be largely affected by changes in hydrological properties of the soil. Since parameters were varied within ranges that approximated their uncertainty, a first-step assessment of uncertainty for the predictions of pesticide losses was possible. Large uncertainties in the predictions were reported, although these are likely to have been overestimated by considering a large number of input parameters in the exercise. It appears desirable that a probabilistic framework accounting for uncertainty is integrated into the estimation of pesticide exposure for regulatory purposes.     

Turfgrass thatch effects on pesticide leaching: a laboratory and modeling study

Process-based models are frequently used to assess the water quality impacts of turfgrass management emanating from proposed or existing golf courses. Thatch complicates the prediction of pesticide transport because surface-applied pesticides must pass through an organic-rich layer before entering the soil. This study was conducted to (i) compare the use of a linear equilibrium model (LEM) and two-site nonequilibrium (2SNE) model to predict pesticide transport through soil and thatch + soil columns, and (ii) evaluate thatch effects on pesticide transport through soil columns with a volume-averaging approach. Pesticide breakthrough curves were obtained for soil and thatch + soil columns from a 1 cm h(-1) flux applied one day after applying triclopyr (3,5,6-trichloro-2-pyridinyloxyacetic acid) and carbaryl (1-napthyl-methyl carbamate). Pesticide and bromide transport parameters indicated that nonequilibrium processes were affecting pesticide transport. Columns containing zoysiagrass (Zoysia japonica Steud.) thatch had lower triclopyr and carbaryl leaching losses than did soil-only columns, although total reductions attributable to thatch did not exceed 15% of the applied pesticide. When laboratory-based retardation factors were used, the 2SNE model explained 88 to 93% of the variability for triclopyr and 70 to 94% of the variability for carbaryl. Laboratory-based retardation factors performed well in a 2SNE model to predict the peak concentration and tailing behavior of triclopyr and carbaryl with a volume-averaging approach. These results suggest that separate representation of the thatch layer in process-based models is not a prerequisite to obtain reasonable estimates of pesticide transport under steady state flow conditions.     

Mapping ground water vulnerability to pesticide leaching with a process-based metamodel of EuroPEARL

To support EU policy, indicators of pesticide leaching at the European level are required. For this reason, a metamodel of the spatially distributed European pesticide leaching model EuroPEARL was developed. EuroPEARL considers transient flow and solute transport and assumes Freundlich adsorption, first-order degradation and passive plant uptake of pesticides. Physical parameters are depth dependent while (bio)-chemical parameters are depth, temperature, and moisture dependent. The metamodel is based on an analytical expression that describes the mass fraction of pesticide leached. The metamodel ignores vertical parameter variations and assumes steady flow. The calibration dataset was generated with EuroPEARL and consisted of approximately 60,000 simulations done for 56 pesticides with different half-lives and partitioning coefficients. The target variable was the 80th percentile of the annual average leaching concentration at 1-m depth from a time series of 20 yr. The metamodel explains over 90% of the variation of the original model with only four independent spatial attributes. These parameters are available in European soil and climate databases, so that the calibrated metamodel could be applied to generate maps of the predicted leaching concentration in the European Union. Maps generated with the metamodel showed a good similarity with the maps obtained with EuroPEARL, which was confirmed by means of quantitative performance indicators.     

GROUNDWATER CONTAMINATION BY NEMATICIDES: INFLUENCE OF RECHARGE TIMING UNDER PINEAPPLE CROP1

ABSTRACT: Toxic organic compounds, such as DBCP, EDB, and c TCP, that are associated with pineapple cultivation in Hawaii have been discovered in drinking water wells on Oahu. In order to reach and contaminate the Pearl Harbor aquifer, pesticides must be transported quickly downward away from the soil surface prior to complete volatilization, degradation, or adsorption of residuals. This paper assesses the role of pesticide application timing relative to subsequent rainfall-induced recharge events in determining the amount and extent of chemical leaching from the soil. A water balance model for a pineapple crop is developed to estimate the time series of recharge from two fields for which soil contamination profiles are available. In general, the amounts of DBCP, EDB, and TCP found in the soil profiles of the two fields are consistent with expectations of leaching based on an analysis of the recharge time series. The results indicate that recharge during and immediately following the application of pesticides is important in determining whether groundwater contamination will result.     

Estimating Watershed Level Nonagricultural Pesticide Use From Golf Courses Using Geospatial Methods1

Abstract: Limited information exists on pesticide use for nonagricultural purposes, making it difficult to estimate pesticide loadings from nonagricultural sources to surface water and to conduct environmental risk assessments. A method was developed to estimate the amount of pesticide use on recreational turf grasses, specifically golf course turf grasses, for watersheds located throughout the conterminous United States (U.S.). The approach estimates pesticide use: (1) based on the area of recreational turf grasses (used as a surrogate for turf associated with golf courses) within the watershed, which was derived from maps of land cover, and (2) from data on the location and average treatable area of golf courses. The area of golf course turf grasses determined from these two methods was used to calculate the percentage of each watershed planted in golf course turf grass (percent crop area, or PCA). Turf‐grass PCAs derived from the two methods were used with recommended application rates provided on pesticide labels to estimate total pesticide use on recreational turf within 1,606 watersheds associated with surface‐water sources of drinking water. These pesticide use estimates made from label rates and PCAs were compared to use estimates from industry sales data on the amount of each pesticide sold for use within the watershed. The PCAs derived from the land‐cover data had an average value of 0.4% of a watershed with minimum of 0.01% and a maximum of 9.8%, whereas the PCA values that are based on the number of golf courses in a watershed had an average of 0.3% of a watershed with a minimum of <0.01% and a maximum of 14.2%. Both the land‐cover method and the number of golf courses method produced similar PCA distributions, suggesting that either technique may be used to provide a PCA estimate for recreational turf. The average and maximum PCAs generally correlated to watershed size, with the highest PCAs estimated for small watersheds. Using watershed specific PCAs, combined with label rates, resulted in greater than two orders of magnitude over‐estimation of the pesticide use compared to estimates from sales data.     

Ammonia volatilization from surface-applied poultry litter under conservation tillage management practices

Land application of poultry litter can provide essential plant nutrients for crop production, but ammonia (NH(3)) volatilization from the litter can be detrimental to the environment. A multiseason study was conducted to quantify NH(3) volatilization rates from surface-applied poultry litter under no-till and paraplowed conservation tillage managements. Litter was applied to supply 90 to 140 kg N ha(-1). Evaluation of NH(3) volatilization was determined using gas concentrations and the flux-gradient gas transport technique using the momentum balance transport coefficient. Ammonia fluxes ranged from 3.3 to 24% of the total N applied during the winter and summer, respectively. Ammonia volatilization was rapid immediately after litter application and stopped within 7 to 8 d. Precipitation of 17 mm essentially halted volatilization, probably by transporting litter N into the soil matrix. Application of poultry to conservation-tilled cropland immediately before rainfall events would reduce N losses to the atmosphere but could also increase NO(3) leaching and runoff to streams and rivers.     

MODEL FOR PREDICTING TRANSPORT OF PESTICIDES AND THEIR METABOLITES IN THE UNSATURATED ZONE1

ABSTRACT: A finite element model based on Galerkin's upstream weighted residual technique was developed to predict the simultaneous convective-dispersion transport and transformations of pesticides and their metabolites in the unsaturated zone. Transformations of the parent compound and its metabolites were assumed to be first-order reactions for oxidation and hydrolysis, while adsorption of the pesticide species (parent compound and metabolites) to the soil components was assumed to be represented by a linear equilibrium (Freundlich type) isotherm. Volatilization and plant root uptake of pesticides in the solution phase were neglected in the analysis. The proposed model was used to simulate the transport and transformation of aldicarb and its metabolites, aldicarb sulfoxide and aldicarb sulfone, in the soil profile. Several examples are used to demonstrate the accuracy, validity, and applicability of the proposed model. Simulated results indicate that the proposed model can potentially be used to estimate the mass flux of water, and pesticide and pesticide metabolite concentrations in the subsurface environment. However, further verification of the model against actual field data is needed to fully demonstrate the model's potential.     

Selection of plants for optimization of vegetative filter strips treating runoff from turfgrass

Runoff from turf environments, such as golf courses, is of increasing concern due to the associated chemical contamination of lakes, reservoirs, rivers, and ground water. Pesticide runoff due to fungicides, herbicides, and insecticides used to maintain golf courses in acceptable playing condition is a particular concern. One possible approach to mitigate such contamination is through the implementation of effective vegetative filter strips (VFS) on golf courses and other recreational turf environments. The objective of the current study was to screen ten aesthetically acceptable plant species for their ability to remove four commonly-used and degradable pesticides: chlorpyrifos (CP), chlorothalonil (CT), pendimethalin (PE), and propiconazole (PR) from soil in a greenhouse setting, thus providing invaluable information as to the species composition that would be most efficacious for use in VFS surrounding turf environments. Our results revealed that blue flag iris (Iris versicolor) (76% CP, 94% CT, 48% PE, and 33% PR were lost from soil after 3 mo of plant growth), eastern gama grass (Tripsacum dactyloides) (47% CP, 95% CT, 17% PE, and 22% PR were lost from soil after 3 mo of plant growth), and big blue stem (Andropogon gerardii) (52% CP, 91% CT, 19% PE, and 30% PR were lost from soil after 3 mo of plant growth) were excellent candidates for the optimization of VFS as buffer zones abutting turf environments. Blue flag iris was most effective at removing selected pesticides from soil and had the highest aesthetic value of the plants tested.     

Point- and nonpoint-source pesticide contamination in the Zwester Ohm catchment, Germany

Reducing pesticide loads in surface waters implies identifying the pathways responsible for the pollution. The current study documents the pesticide contamination of the river Zwester Ohm, a 4917-ha catchment in Germany with 41% of the land used for crop production. Discharges and concentrations of 19 pesticides were measured continuously at three locations for 15 mo. The load detected at the outlet of the catchment amounted to 9048 g a.i. The losses represent 0.22% of the pesticides applied by the farmers. The contamination showed a seasonal pattern following the pesticide application times. The wastewater treatment plant system (WWTPS) in the catchment (two wastewater treatment plants [WWTP], 14 combined sewer overflows (CSO), four CSO tanks) emits during dry weather periods purified sewage and during storm events sewage mixed with stormwater runoff into the river. The contribution by the WWTPS to the pesticide load was defined as point-source pollution (PSP). The load was dominated by PSP with at least 77% of the total pollution. No significant interdependencies between intrinsic properties of the pesticides, hydrometeorological factors, and the loads occurring in the stream could be found. Therefore, it is not possible to predict PSP for other catchments based on the results from this study. Whereas 65% of the total load entered the river via the WWTP, a portion of 12% was attributed to the CSO. The study points out that the influence of CSO on PSP should be taken into account in future catchment studies in areas with comparable agricultural structure.     

Partitioning and persistence of trichlorfon and chlorpyrifos in a creeping bentgrass putting green

Golf course putting greens typically receive high pesticide applications to meet high quality demands. Research on pesticide fate in turf ecosystems is important to better understand the potential impact of pesticide use on the environment and human health. This research was conducted to evaluate the environmental fate of two commonly used insecticides--trichlorfon (dimethyl 2,2,2-trichloro-1-hydroxyethylphosphonate) and chlorpyrifos (O,O-diethyl O-3,5,6-trichloro-2-pyridylphosphorothioate)--in a creeping bentgrass (Agrostis palustris Huds.) putting green under customary field management practices at the University of California-Riverside Turf Research Facility during 1996 and 1997. The two insecticides were chosen because of their difference in water solubility, persistence, adsorption, and vapor pressure. Volatilization, clipping removal, and soil residues of the insecticides were quantified and leaching was monitored using lysimeters installed in putting green plots. Results showed trichlorfon volatilization, clipping removal, and leaching loss was insignificant (in the range of 0.0001-0.06% of applied mass) both in 1996 and 1997. No significant difference in clipping removal of trichlorfon and chlorpyrifos was observed in both years (0.06 and 0.05% of applied mass for trichlorfon and 0.15 and 0.19% of applied mass for chlorpyrifos, respectively, in 1996 and 1997), but significantly lower cumulative leaching and lower soil concentration was observed in 1997 than in 1996. Volatilization loss of chlorpyrifos was not significantly different between 1996 (2.05%) and 1997 (2.71%). Volatilization loss of trichlorfon in 1996 (0.01%) was significantly higher than in 1997 (0.008%). This study demonstrated the fraction of applied insecticides leaving the turf putting greens was minimal.     

Mapping and Modeling the Biogeochemical Cycling of Turf Grasses in the United States

Turf grasses are ubiquitous in the urban landscape of the United States and are often associated with various types of environmental impacts, especially on water resources, yet there have been limited efforts to quantify their total surface and ecosystem functioning, such as their total impact on the continental water budget and potential net ecosystem exchange (NEE). In this study, relating turf grass area to an estimate of fractional impervious surface area, it was calculated that potentially 163,800 km² (± 35,850 km²) of land are cultivated with turf grasses in the continental United States, an area three times larger than that of any irrigated crop. Using the Biome-BGC ecosystem process model, the growth of warm-season and cool-season turf grasses was modeled at a number of sites across the 48 conterminous states under different management scenarios, simulating potential carbon and water fluxes as if the entire turf surface was to be managed like a well-maintained lawn. The results indicate that well-watered and fertilized turf grasses act as a carbon sink. The potential NEE that could derive from the total surface potentially under turf (up to 17 Tg C/yr with the simulated scenarios) would require up to 695 to 900 liters of water per person per day, depending on the modeled water irrigation practices, suggesting that outdoor water conservation practices such as xeriscaping and irrigation with recycled waste-water may need to be extended as many municipalities continue to face increasing pressures on freshwater.     

Runoff transport of pyrethroids from a residential lawn in central California

An irrigation runoff study on a residential lawn was conducted in California, northeast of Sacramento, during the summer and fall of 2008 to investigate the contribution of turf uses of pyrethroids to residues in Californian urban creek sediments. This study examined how over irrigation (i.e., irrigation that produces runoff) in the summer season may transport recently applied pyrethroids. The study included liquid and granular applications of both bifenthrin [(2-methyl-3-phenyl-phenyl) methyl 3-(2-chloro-3,3,3-trifluoro-prop-1-enyl)-2,2-dimethyl-cyclopropane-1-carboxylate] and beta-cyfluthrin [Cyano(4-fluoro-3-phenoxyphenyl)methyl 3-(2,2-dichloroethenyl)-2,2-dimethyl-cyclopropanecarboxylate]. Generally, runoff did not occur at irrigation rates of 2.03 cm/h (0.8 in/h) but did occur when the irrigation rates were increased to about 3.81 cm/h (1.5 in/h), generating chemical losses in the first runoff event of up to 0.58 and 0.08% of applied for beta-cyfluthrin and bifenthrin, respectively. Chemical runoff losses dropped significantly between over-irrigation events with the third over-irrigation event chemical runoff losses representing 0.026 and 0.015% of applied for beta-cyfluthrin and bifenthrin, respectively. Runoff losses were generally less for liquid formulations than granular formulations but within a factor of three. Additionally, the study included a simulated winter rainstorm 8 wk after application. The low runoff losses from turf seen in this study suggest that other sources could be contributing to observed residues in urban streams. Other sources could include pyrethroids ending up on impervious surfaces, such as concrete driveways from off-target applications to turf, spills, and other poor handling practices, or pyrechroids applied directly to impervious surfaces for insect control.