Modeling pesticide volatilization from turf

Pesticide volatilization models are typically based on equilibrium partitioning of the chemical into solid, liquid, and gaseous phases in the soil environment. In turf systems direct vaporization from vegetation surfaces is a more likely source, and it is difficult to apply equilibrium methods to plant material due to the uncertainties of solid-liquid-gas partitioning. An alternative approach is to assume that pesticide volatilization is governed by the same processes that affect water evaporation. A model was developed in which evapotranspiration values, as determined by the Penman equation, were adjusted to chemical vaporization using ratios of water and chemical saturated vapor pressures and latent heats of vaporization. The model also assumes first-order degradation of pesticide on turf vegetation over time. The model was tested by comparisons of predictions with measurements of volatilization for eight pesticides measured during 3 to 7 d in 11 field experiments. Measured volatilization fluxes ranged from 0.1 to 22% of applied chemical. Pesticides were divided into two groups based on saturated vapor pressures and organic C partition coefficients. One pesticide was selected from each group to calibrate the model's volatilization constant for the group, and the remaining pesticides were used for model testing. Testing results indicated that the model provides relatively conservative estimates of pesticide volatilization. Predicted mean losses exceeded 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.     

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.     

EVAPOTRANSPIRATION MEASUREMENTS AND MODELING FOR THREE WETLAND SYSTEMS IN SOUTH FLORIDA1

ABSTRACT: At the Everglades Nutrient Removal project in south Florida, three lysimeters were installed to measure daily evapotranspiration (ET) rates from cattails (Typha domingensis), mixed marsh vegetation, and an open water/algae system. The cattail lysimeter began operation in February 1993. The mixed marsh vegetation lysimeter began operation in January 1994, and the open water lysimeter with occasional algae cover began operation in December 1993. The mean measured ET rate was 3.6 mm, 3.5 mm, and 3.7 mm per day for the cattail, mixed marsh vegetation, and open water/algae system, respectively. High resolution weather data were continuously measured at the site. Six models were applied to estimate daily ET rates of the three systems. The Penman-Monteith equation best estimated ET of cattail and mixed marsh vegetation, and the Penman Combination equation was most suitable for the open water/algae system. Empirical equations based on solar radiation and maximum temperature produced estimates of daily ET from the three systems that are comparable to models that require many more parameters. In cases where limited data is available, the calibrated simple models can be used to estimate ET from wetlands in south Florida.     

Do Energy‐Based PET Models Require More Input Data than Temperature‐Based Models? — An Evaluation at Four Humid FluxNet Sites

It is well established that wet environment potential evapotranspiration (PET) can be reliably estimated using the energy budget at the canopy or land surface. However, in most cases the necessary radiation measurements are not available and, thus, empirical temperature‐based PET models are still widely used, especially in watershed models. Here we question the presumption that empirical PET models require fewer input data than more physically based models. Specifically, we test whether the energy‐budget‐based Priestley‐Taylor (P‐T) model can reliably predict daily PET using primarily air temperature to estimate the radiation fluxes and associated parameters. This method of calculating PET requires only daily minimum and maximum temperature, day of the year, and latitude. We compared PET estimates using directly measured radiation fluxes to PET calculated from temperature‐based radiation estimates at four humid AmeriFlux sites. We found good agreement between P‐T PET calculated from measured radiation fluxes and P‐T PET determined via air temperature. In addition, in three of the four sites, the temperature‐based radiation approximations had a stronger correlation with measured evapotranspiration (ET) during periods of maximal ET than fully empirical Hargreaves, Hamon and Oudin methods. Of the three fully empirical models, the Hargreaves performed the best. Overall, the results suggest that daily PET estimates can be made using a physically based approach even when radiation measurements are unavailable.     

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.     

Assessing Calibration Uncertainty and Automation for Estimating Evapotranspiration from Agricultural Areas Using METRIC

Agricultural irrigation accounts for a large fraction of the total water use in the western United States. The Mapping Evapotranspiration at high Resolution with Internalized Calibration (METRIC) remote sensing energy balance model is being used to estimate historical agricultural water use in western Nevada to evaluate basin‐wide water budgets. Each METRIC evapotranspiration (ET) estimate must be calibrated by a trained user, which requires some iterative time investment and results in variation in ET estimates between users. An automated calibration algorithm for the METRIC model was designed to generate ET estimates comparable to those from trained users by mimicking the manual calibration process. Automated calibration allows for rapid generation of METRIC ET estimates with minimal manual intervention, as well as uncertainty and sensitivity analysis of the model. The variation in ET estimates generated by the automated calibration algorithm was found to be similar to the variation in manual ET estimates. Results indicate that uncertainty was highest for fields with low ET levels and lowest for fields with high ET levels, with a seasonal mean uncertainty of approximately 5% for all fields. In addition, in a blind comparison, automated daily and seasonal ET estimates compared well with flux tower measurement ET data at multiple sites. Automated methods can generate first‐order ET estimates that are similar to time intensive manual efforts with less time investment.     

Using Landsat data to estimate evapotranspiration of winter wheat

An evapotranspiration (ET) model that accurately estimates daily water use and soil moisture on a regional basis is required for many agricultural and hydrological studies. The model should use meterological data that are readily available and crop information that is responsive to the changing vigor of the plants.We evaluated an ET model with a weighing lysimeter and then applied it to winter wheatfields at four Kansas locations. Model inputs are solar radiation, temperature, precipitation, and leaf area index (LAI); included in the outputs are estimates of transpiration, evaporation, and soil moisture. An equation was developed to estimate LAI from Landsat data. Because LAI can be estimated from satellites, the ET model can potentially be used on a regional basis.     

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.     

AUTOMATED METHODS FOR ESTIMATING BASEFLOW AND GROUND WATER RECHARGE FROM STREAMFLOW RECORDS1

ABSTRACT: To quantify and model the natural ground water recharge process, six sites located in the midwest and eastern United States where previous water balance observations had been made were compared to computerized techniques to estimate: (1) base flow and (2) ground water recharge. Results from an existing automated digital filter technique for separating baseflow from daily streamflow records were compared to baseflow estimates made in the six water balance studies. Previous validation of automated baseflow separation techniques consisted only of comparisons with manual techniques. In this study, the automated digital filter technique was found to compare well with measured field estimates yielding a monthly coefficient of determination of 0.86. The recharge algorithm developed in this study is an automated derivation of the Rorabaugh hydrograph recession curve displacement method that utilizes daily streamflow. Comparison of annual recharge from field water balance measurements to those computed with the automated recession curve displacement method had coefficients of determination of 0.76 and predictive efficiencies of 71 percent. Monthly estimates showed more variation and are not advocated for use with this method. These techniques appear to be fast, reproducible methods for estimating baseflow and annual recharge and should be useful in regional modeling efforts and as a quick check on mass balance techniques for shallow water table aquifers.     

A SENSITIVITY AND ERROR ANALYSIS CF PROCEDURES USED FOR ESTIMATING EVAPORATION1

ABSTRACT: Sensitivity and error analyses were used to examine the following objectives: (1) analyze the structure of commonly used evaporation models; (2) provide estimates of the effect of variation in meteorological factors on observed evaporation rates; and (3) estimate the effect of error in measurements of the meteorological factors. The results indicate error in evaporation estimates resulting from measurement error in meteorological factors is probably much less than five percent of the computer evaporation rate. Variation with both time and space of the importance of the different meteorological factors is demonstrated. The sensitivity analysis indicates that the Fractional-Evaporation Equivalent method is structurally inadequate and the Weather Bureau model is more flexible than the Penman model. However, the Penman model appears to provide more realistic estimates of the importance of the various meteorological factors.     

Global Daily Reference Evapotranspiration Modeling and Evaluation1

Abstract: Accurate and reliable evapotranspiration (ET) datasets are crucial in regional water and energy balance studies. Due to the complex instrumentation requirements, actual ET values are generally estimated from reference ET values by adjustment factors using coefficients for water stress and vegetation conditions, commonly referred to as crop coefficients. Until recently, the modeling of reference ET has been solely based on important weather variables collected from weather stations that are generally located in selected agro‐climatic locations. Since 2001, the National Oceanic and Atmospheric Administration’s Global Data Assimilation System (GDAS) has been producing six‐hourly climate parameter datasets that are used to calculate daily reference ET for the whole globe at 1‐degree spatial resolution. The U.S. Geological Survey Center for Earth Resources Observation and Science has been producing daily reference ET (ETo) since 2001, and it has been used on a variety of operational hydrological models for drought and streamflow monitoring all over the world. With the increasing availability of local station‐based reference ET estimates, we evaluated the GDAS‐based reference ET estimates using data from the California Irrigation Management Information System (CIMIS). Daily CIMIS reference ET estimates from 85 stations were compared with GDAS‐based reference ET at different spatial and temporal scales using five‐year daily data from 2002 through 2006. Despite the large difference in spatial scale (point vs. ～100 km grid cell) between the two datasets, the correlations between station‐based ET and GDAS‐ET were very high, exceeding 0.97 on a daily basis to more than 0.99 on time scales of more than 10 days. Both the temporal and spatial correspondences in trend/pattern and magnitudes between the two datasets were satisfactory, suggesting the reliability of using GDAS parameter‐based reference ET for regional water and energy balance studies in many parts of the world. While the study revealed the potential of GDAS ETo for large‐scale hydrological applications, site‐specific use of GDAS ETo in complex hydro‐climatic regions such as coastal areas and rugged terrain may require the application of bias correction and/or disaggregation of the GDAS ETo using downscaling techniques.     

WATER BUDGETS FOR SMALL DIKED MARSHES1

ABSTRACT: This paper reports an analysis of the water budgets of 10 small (5–6 ha) diked areas (cells) within the Delta Marsh in southcentral Manitoba, Canada. The important terms in the water budget equation in this study were precipitation (P), water pumped in (SWI), evapotranspiration (ET), seepage in (GWI) and out (GWO), and change in storage (ΔS). P, SWI, and S were measured directly, and the sum of ET and GWO determined by difference. Estimating ET as 0.7 pan evaporation gave a seepage loss of 2.9 mm/day from the most intensively studied cell. Other methods of estimating ET produced estimates of GWO ranging from 2.4 to 3.8 mm/day. Water budgets for less intensively studied cells indicated seepage loss increased as perimeter available for seepage increased, but not proportionately. Efforts to measure seepage directly or estimate it from measured hydraulic gradients and hydraulic conductivity produced estimates much lower than the estimates from the water budget equation. Hydraulic conductivities were very heterogeneous, reflecting the sorting of water deposited sediments. Comparison of the hydraulic conductivities with seepage estimates from the water budget strongly suggests water movement downward as well as laterally from these diked areas.     

A Satellite Model of Forest Flammability

We describe a model of forest flammability, based on daily satellite observations, for national to regional applications. The model defines forest flammability as the percent moisture content of fuel, in the form of litter of varying sizes on the forest floor. The model uses formulas from the US Forest Service that describe moisture exchange between fuel and the surrounding air and precipitation. The model is driven by estimates of temperature, humidity, and precipitation from the moderate resolution imaging spectrometer and tropical rainfall measuring mission multi-satellite precipitation analysis. We provide model results for the southern Amazon and northern Chaco regions. We evaluate the model in a tropical forest-to-woodland gradient in lowland Bolivia. Results from the model are significantly correlated with those from the same model driven by field climate measurements. This model can be run in a near real-time mode, can be applied to other regions, and can be a cost-effective input to national fire management programs.     

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.     

ESTIMATION OF REGIONAL EVAPOTRANSPIRATION FROM ARID AND SEMI-ARID SURFACES1

ABSTRACT: Evapotranspiration (ET) from oasis and Gobi surfaces in the Heihe River region of China was estimated by Kotoda (1986) and Advection-Aridity (Brutsaert and Stricker, 1979) models. The ET estimates were compared with eddy-correlation ET estimates. The ET estimated by both models are in good agreement with ET obtained by means of eddy-correlation method for the oasis surface while underestimating ET as compared to the eddy-correlation method. For the Gobi surface, the models yielded obviously overestimates of evapotranspiration. In order to estimate evapotranspiration from arid surfaces, the Kotoda model was modified by introducing the surface moisture availability parameter a from Barton (1979). The modified Kotoda model yielded ET estimates that were very close to that from the eddy-correlation method for the Gobi surface. The modified Kotoda model was used to estimate evapotranspiration from the Heihe River watershed, an area with complicated topography and land use, and the results compared with those from a water balance method. A sensitivity analysis of the modified model was performed. The results show that the modified Kotoda model can reflect the relationship between the actual evapotranspiration and the main controlling factors on it for both wet and arid surfaces reasonably. From this study, it can be concluded that the modified Kotoda model is applicable fro the estimation of regional evapotranspiration from areas with complicated topography and land use.     

ESTIMATING AVERAGE MONTHLY LAKE EVAPORATION IN THE NORTHEAST UNITED STATES1

ABSTRACT: A methodology to estimate the average monthly lake evaporation, E(τ), (month τ=1,12) for fresh water bodies located in the northeast United States is presented. The approach combines analysis of at‐site, lake‐specific vertical water temperature profile data and a previously developed regional air temperature based model approximation of the widely accepted modified Penman energy budget estimate of mean monthly potential evaporation, E_p(τ) (mm/day). The paper presents procedures to develop site‐specific estimates of E_p(τ) and to convert water temperature data to average monthly conductive heat flux, G(τ). With monthly estimates of G(τ), the average monthly potential evaporation, E_p(τ), is then convertible to estimates of the average monthly lake evaporation, E(τ). This new method permits a good estimate of site‐specific lake evaporation rates without the data and computational requirements of the Penman energy budget procedure nor the comparatively expensive, time consuming field eddy correlation approach.     

Annual Estimates of Recharge,Quick‐Flow Runoff,and Evapotranspiration for the Contiguous U.S. Using Empirical Regression Equations

This study presents new data‐driven, annual estimates of the division of precipitation into the recharge, quick‐flow runoff, and evapotranspiration (ET) water budget components for 2000‐2013 for the contiguous United States (CONUS). The algorithms used to produce these maps ensure water budget consistency over this broad spatial scale, with contributions from precipitation influx attributed to each component at 800 m resolution. The quick‐flow runoff estimates for the contribution to the rapidly varying portion of the hydrograph are produced using data from 1,434 gaged watersheds, and depend on precipitation, soil saturated hydraulic conductivity, and surficial geology type. Evapotranspiration estimates are produced from a regression using water balance data from 679 gaged watersheds and depend on land cover, temperature, and precipitation. The quick‐flow and ET estimates are combined to calculate recharge as the remainder of precipitation. The ET and recharge estimates are checked against independent field data, and the results show good agreement. Comparisons of recharge estimates with groundwater extraction data show that in 15% of the country, groundwater is being extracted at rates higher than the local recharge. These maps of the internally consistent water budget components of recharge, quick‐flow runoff, and ET, being derived from and tested against data, are expected to provide reliable first‐order estimates of these quantities across the CONUS, even where field measurements are sparse.     

Estimating the Extent of Impervious Surfaces and Turf Grass across Large Regions

The ability of researchers to accurately assess the extent of impervious and pervious developed surfaces, e.g., turf grass, using land‐cover data derived from Landsat satellite imagery in the Chesapeake Bay watershed is limited due to the resolution of the data and systematic discrepancies between developed land‐cover classes, surface mines, forests, and farmlands. Estimates of impervious surface and turf grass area in the Mid‐Atlantic, United States that were based on 2006 Landsat‐derived land‐cover data were substantially lower than estimates based on more authoritative and independent sources. New estimates of impervious surfaces and turf grass area derived using land‐cover data combined with ancillary information on roads, housing units, surface mines, and sampled estimates of road width and residential impervious area were up to 57 and 45% higher than estimates based strictly on land‐cover data. These new estimates closely approximate estimates derived from authoritative and independent sources in developed counties.     

Evaluation of the macroalga,muskgrass, for the phytoremediation of selenium-contaminated agricultural drainage water by microcosms

Previous field studies suggested that the macroalga, muskgrass (Chara canescens Desv. & Lois), plays an important role in the removal of selenium (Se) from agricultural drainage water. This study evaluated the efficiency of Se removal from drainage water by muskgrass-vegetated wetland microcosms, and determined the extent to which muskgrass removed Se through phytoextraction and biovolatilization. Six flow-through wetland microcosms were continuously supplied with drainage water containing an average Se concentration of 22 microg L(-1) over a 24-d experimental period. The Se mass input and outflow and the rate of Se volatilization were monitored daily for each microcosm. Three microcosms containing muskgrass reduced the daily mass Se input in the inflow drainage water by 72.1%; this compared with a reduction of 50.6% of the mass Se input for three unvegetated control microcosms. Selenium accumulated in muskgrass tissues accounted for 1.9% of the total mass Se input in the microcosm, followed by 0.5% via biological volatilization. The low rates of Se volatilization from selenate-supplied muskgrass, which were 10-fold less than from selenite, were probably due to a major rate limitation in the reduction of selenate to organic forms of Se in muskgrass. This conclusion was derived from X-ray absorption spectroscopy speciation analysis, which showed that muskgrass treated with selenite contained 91% of the total Se in organic forms (selenoethers and diselenides), compared with 47% in muskgrass treated with selenate.