The effect of the runoff size on the pesticide concentration in runoff water and in FOCUS streams simulated by PRZM and TOXSWA.

Within the European Union the exposure of aquatic organisms to pesticides is assessed by simulations with the so-called FOCUS Surface Water Scenarios. Runoff plays an important role in these scenarios. As little is known about the effect of runoff size on the exposure, we investigated the effect of runoff size on the concentration in the runoff water and in streams simulated with the PRZM and TOXSWA models for two FOCUS runoff scenarios. For weakly sorbing pesticides (KF,oc<100Lkg-1) the pesticide concentration in the runoff water decreased exponentially with increasing daily runoff size. The runoff size hardly affected the pesticide concentration in the runoff water of strongly sorbing pesticides (KF,oc≥1000Lkg-1). For weakly sorbing pesticides the concentration in the FOCUS stream reached a maximum at runoff sizes of about 0.3 to 1mm. The concentration increased rapidly when the runoff size increased from 0 to 0.1mm and gradually decreased when runoff exceeded 1mm. For strongly sorbing pesticides the occurrence of the maximum concentration in the stream is clearly less pronounced and lies approximately between 1 and 20mm runoff. So, this work indicates that preventing small runoff events (e.g. by vegetated buffer strips) reduces exposure concentrations strongly for weakly sorbing pesticides. A simple metamodel was developed for the ratio between the concentrations in the stream and in the runoff water. This model predicted the ratios simulated by TOXSWA very well and it demonstrated that (in addition to runoff size and concentration in runoff) the size of the pesticide-free base flow and pesticide treatment ratio of the catchment determine the stream concentration to a large extent.

Effects of aged sorption on pesticide leaching to groundwater simulated with PEARL.

Leaching to groundwater is an important element of the regulatory risk assessment of pesticides in western countries. Including aged sorption in this assessment is relevant because there is ample evidence of this process and because it leads to a decrease in simulated leaching. This work assesses the likely magnitude of this decrease for four groundwater scenarios used for regulatory purpose in the EU (from the UK, Portugal, Austria and Greece) and for ranges of aged-sorption parameters and substance properties using the PEARL model. Three aged-sorption parameters sets were derived from literature, representing approximately 5th, 50th and 95th percentile cases for the magnitude of the effect of aged sorption on leaching concentrations (called S, M and L, respectively). The selection of these percentile cases was based only on the fNE parameter (i.e. the ratio of the aged sorption and the equilibrium sorption coefficients) because leaching was much more affected by the uncertainty in this parameter than by the uncertainty in the desorption rate coefficient of these sites (kd). For the UK scenario, the annual flux concentration of pesticide leaching at 1m depth decreased by typically a factor of 5, 30 and >1000 for the S, M and L parameter sets, respectively. This decrease by a factor of 30 for the M parameter set appeared to be approximately valid also for the other three scenarios. Decreasing the Freundlich exponent N from 0.9 into 0.7 for the M parameter set, increased this factor of 30 into a factor of typically 1000, considering all four scenarios. The aged-sorption sites were close to their equilibrium conditions during the leaching simulations for two of the four scenarios (for all substances considered and the M parameter set), but this was not the case for the other two scenarios.

An energetics-based honeybee nectar-foraging model used to assess the potential for landscape-level pesticide exposure dilution.

Estimating the exposure of honeybees to pesticides on a landscape scale requires models of their spatial foraging behaviour. For this purpose, we developed a mechanistic, energetics-based model for a single day of nectar foraging in complex landscape mosaics. Net energetic efficiency determined resource patch choice. In one version of the model a single optimal patch was selected each hour. In another version, recruitment of foragers was simulated and several patches could be exploited simultaneously. Resource availability changed during the day due to depletion and/or intrinsic properties of the resource (anthesis). The model accounted for the impact of patch distance and size, resource depletion and replenishment, competition with other nectar foragers, and seasonal and diurnal patterns in availability of nectar-providing crops and wild flowers. From the model we derived simple rules for resource patch selection, e.g., for landscapes with mass-flowering crops only, net energetic efficiency would be proportional to the ratio of the energetic content of the nectar divided by distance to the hive. We also determined maximum distances at which resources like oilseed rape and clover were still energetically attractive. We used the model to assess the potential for pesticide exposure dilution in landscapes of different composition and complexity. Dilution means a lower concentration in nectar arriving at the hive compared to the concentration in nectar at a treated field and can result from foraging effort being diverted away from treated fields. Applying the model for all possible hive locations over a large area, distributions of dilution factors were obtained that were characterised by their 90-percentile value. For an area for which detailed spatial data on crops and off-field semi-natural habitats were available, we tested three landscape management scenarios that were expected to lead to exposure dilution: providing alternative resources than the target crop (oilseed rape) in the form of (i) other untreated crop fields, (ii) flower strips of different widths at field edges (off-crop in-field resources), and (iii) resources on off-field (semi-natural) habitats. For both model versions, significant dilution occurred only when alternative resource patches were equal or more attractive than oilseed rape, nearby and numerous and only in case of flower strips and off-field habitats. On an area-base, flower strips were more than one order of magnitude more effective than off-field habitats, the main reason being that flower strips had an optimal location. The two model versions differed in the predicted number of resource patches exploited over the day, but mainly in landscapes with numerous small resource patches. In landscapes consisting of few large resource patches (crop fields) both versions predicted the use of a small number of patches.

Effect of Harsh or Mild Extraction of Soil on Pesticide Leaching to Groundwater.

Assessment of leaching to groundwater is an important aspect of pesticide risk assessment. The first leaching tier usually consists of simulations with leaching scenarios based on pesticide-soil properties derived from laboratory studies. Because the extractability of pesticide residues in such studies decreases with time, the harshness of the extraction method influences these pesticide-soil properties. This study investigates the effect of using a mild or harsh extraction method on simulated leaching to groundwater with consideration of substances with a range of half-lives and organic matter sorption coefficient values for selected leaching scenarios. The model for linking the concentrations of the mild and the harsh systems was based on laboratory studies with two pesticides and a Dutch sandy soil and was tested against Canadian field studies with atrazine (6-chloro-2-ethyl-4-isopropyl-1,3,5-triazine-2,4-diamine). The degradation rate and the aged-sorption parameters of each "mild" soil-substance system were derived from a hypothetical laboratory incubation study using prescribed parameter values for the corresponding "harsh" soil-substance system. Simulations were performed for three European leaching scenarios (United Kingdom, France, Portugal). For the best-guess parameter set, the leaching concentrations of the harsh system were approximately equal to those of the mild system at leaching concentrations greater than 1 µg L and were at most approximately a factor of two higher than those of the mild systems at mild leaching concentrations between 0.01 and 0.1 µg L. However, an extreme parameter set led to harsh leaching concentrations that were at most approximately 10 times higher than the mild leaching concentrations at levels between 0.01 and 0.1 µg L.

European scenarios for exposure of soil organisms to pesticides.

Standardised exposure scenarios play an important role in European pesticide authorisation procedures (a scenario is a combination of climate, weather and crop data to be used in exposure models). The European Food Safety Authority developed such scenarios for the assessment of exposure of soil organisms to pesticides. Scenarios were needed for both the concentration in total soil and for the concentration in the liquid phase. The goal of the exposure assessment is the 90th percentile of the exposure concentration in the area of agricultural use of a pesticide in each of three regulatory European zones (North, Centre and South). A statistical approach was adopted to find scenarios that are consistent with this exposure goal. Scenario development began with the simulation of the concentration distribution in the entire area of use by means of a simple analytical model. In the subsequent two steps, procedures were applied to account for parameter uncertainty and scenario uncertainty (i.e. the likelihood that a scenario that is derived for one pesticide is not conservative enough for another pesticide). In the final step, the six scenarios were selected by defining their average air temperature, soil organic-matter content and their soil textural class. Organic matter of the selected scenarios decreased in the order North-Centre-South. Because organic matter has a different effect on the concentration in total soil than it has on the concentration in the liquid phase, the concentration in total soil decreased in the order North-Centre-South whereas the concentration in the liquid phase decreased in the opposite order. The concentration differences between the three regulatory zones appeared to be no more than a factor of two. These differences were comparatively small in view of the considerable differences in climate and soil properties between the three zones.

Estimating degradation rates in outdoor stagnant water by inverse modelling with TOXSWA: a case study with prosulfocarb.

BACKGROUND: The regulatory risk assessment of pesticides requires the assessment of exposure of aquatic ecosystems in small surface waters adjacent to agricultural fields. This exposure is predicted using simulation models, for which an important input parameter is the degradation rate in water. In regulatory dossiers, the decline rate in water from outdoor mesocosms is often available, but this rate encompasses more processes than degradation. Therefore, a procedure was designed for estimating the degradation rate in water that was suitable for mesocosm studies with limited datasets, e.g. datasets lacking site-specific sorption coefficients and relevant sediment properties. The procedure, based upon inverse modelling with TOXSWA, has been tested on a dataset for prosulfocarb in stagnant ditches. RESULTS: A degradation half-life in the ditch water of 2.9 days (20 °C) was found. This short half-life was to a significant extent accounted for by biodegradation rather than hydrolysis or photolysis. This half-life was considerably shorter than the half-life in the water layer of two laboratory water-sediment system experiments. CONCLUSION: The estimation procedure resulted in a unique half-life for the degradation rate in water. Such outdoor mesocosm studies seem to be better suited to assessing the degradation rate in water in ditches than the conventional water-sediment studies.

Simulation of movement of pesticides towards drains with a preferential flow version of PEARL.

BACKGROUND: As part of the Dutch authorisation procedure for pesticides, an assessment of the effects on aquatic organisms in surface waters adjacent to agricultural fields is required. The peak concentration is considered to be the most important exposure endpoint for the ecotoxicological effect assessment. Macropore flow is an important driver for the peak concentration, so the leaching model PEARL was extended with a macropore module. The new model has two macropore domains: a bypass domain and an internal catchment domain. The model was tested against data from a field leaching study on a cracking clay soil in the Netherlands. RESULTS: Most parameters of the model could be obtained from site-specific measurements, pedotransfer functions and general soil structural knowledge; only three macropore-flow-related parameters needed calibration. The flow-related macropore parameters could not be calibrated without using the concentration in drain water. Sequential calibration strategies, in which firstly the water flow model and then the pesticide fate model are calibrated, may therefore be less suitable for preferential flow models. CONCLUSION: After calibration, PEARL could simulate well the observed rapid movement towards drains of two pesticides with contrasting sorption and degradation rate properties. The calibrated value for the fraction of the internal catchment domain was high (90%). This means that a large fraction of water entering the macropores infiltrates into the soil matrix, thus reducing the fraction of rapid flow.

Effect of pesticide fate parameters and their uncertainty on the selection of 'worst-case' scenarios of pesticide leaching to groundwater.

BACKGROUND: For the registration of pesticides in the European Union, model simulations for worst-case scenarios are used to demonstrate that leaching concentrations to groundwater do not exceed a critical threshold. A worst-case scenario is a combination of soil and climate properties for which predicted leaching concentrations are higher than a certain percentile of the spatial concentration distribution within a region. The derivation of scenarios is complicated by uncertainty about soil and pesticide fate parameters. As the ranking of climate and soil property combinations according to predicted leaching concentrations is different for different pesticides, the worst-case scenario for one pesticide may misrepresent the worst case for another pesticide, which leads to 'scenario uncertainty'. RESULTS: Pesticide fate parameter uncertainty led to higher concentrations in the higher percentiles of spatial concentration distributions, especially for distributions in smaller and more homogeneous regions. The effect of pesticide fate parameter uncertainty on the spatial concentration distribution was small when compared with the uncertainty of local concentration predictions and with the scenario uncertainty. CONCLUSION: Uncertainty in pesticide fate parameters and scenario uncertainty can be accounted for using higher percentiles of spatial concentration distributions and considering a range of pesticides for the scenario selection.

Pesticide Leaching from Agricultural Fields with Ridges and Furrows.

In the evaluation of the risk of pesticide leaching to groundwater, the soil surface is usually assumed to be level, although important crops like potato are grown on ridges. A fraction of the water from rainfall and sprinkler irrigation may flow along the soil surface from the ridges to the furrows, thus bringing about an extra load of water and pesticide on the furrow soil. A survey of the literature reveals that surface-runoff from ridges to furrows is a well-known phenomenon but that hardly any data are available on the quantities of water and pesticide involved. On the basis of a field experiment with additional sprinkler irrigation, computer simulations were carried out with the Pesticide Emission Assessment at Regional and Local scales model for separate ridge and furrow systems in a humic sandy potato field. Breakthrough curves of bromide ion (as a tracer for water flow) and carbofuran (as example pesticide) were calculated for 1-m depth in the field. Bromide ion leached comparatively fast from the furrow system, while leaching from the ridge system was slower showing a maximum concentration of about half of that for the furrow system. Carbofuran breakthrough from the furrow system began about a month after application and increased steadily to substantial concentrations. Because the transport time of carbofuran in the ridge soil was much longer, no breakthrough occurred in the growing season. The maximum concentration of carbofuran leaching from the ridge-furrow field was computed to be a factor of six times as high as that computed for the corresponding level field. The study shows that the risk of leaching of pesticides via the furrow soil can be substantially higher than that via the corresponding level field soil.

Measurement and computation of movement of bromide ions and carbofuran in ridged humic-sandy soil.

Water flow and pesticide transport in the soil of fields with ridges and furrows may be more complex than in the soil of more level fields. Prior to crop emergence, the tracer bromide ion and the insecticide carbofuran were sprayed on the humic-sandy soil of a potato field with ridges and furrows. Rainfall was supplemented by sprinkler irrigation. The distribution of the substances in the soil profile of the ridges and furrows was measured on three dates in the potato growing season. Separate ridge and furrow systems were simulated by using the pesticide emission assessment at regional and local scales (PEARL) model for pesticide behavior in soil-plant systems. The substances travelled deeper in the furrow soil than in the ridge soil, because of runoff from the ridges to the furrows. At 19 days after application, the peak of the bromide distribution was measured to be in the 0.1-0.2 m layer of the ridges, while it was in the 0.3-0.5 m layer of the furrows. After 65 days, the peak of the carbofuran distribution in the ridge soil was still in the 0.1 m top layer, while the pesticide was rather evenly distributed in the top 0.6 m of the furrow soil. The wide ranges in concentration measured with depth showed that preferential water flow and substance transport occurred in the sandy soil. Part of the bromide ion distribution was measured to move faster in soil than the computed wave. The runoff of water and pesticide from the ridges to the furrows, and the thinner root zone in the furrows, are expected to increase the risk of leaching to groundwater in ridged fields, in comparison with more level fields.

Simulation of pesticide leaching in a cracking clay soil with the PEARL model.

Testing of pesticide leaching models is important to increase confidence in their use in pesticide registration procedures world-wide. The chromatographic PEARL model was tested against the results of a field leaching study on a cracking clay soil with a tracer (bromide), a mobile pesticide (bentazone) and a moderately sorbing, persistent pesticide (imidacloprid). Input parameters for water flow and solute transport were obtained from site-specific measurements and from literature. The model was tested using a stepwise approach in which each sub-model was sequentially and separately tested. Uncalibrated simulations for the water flow resulted in moisture profiles that were too wet. Calibration of the hydraulic relationships resulted in a good description of the moisture profiles. Calibration of the dispersion length was necessary to obtain a good description of bromide leaching. The calibrated dispersion length was 61 cm, which is very long and indicates a large non-uniformity of solute transport. The half-life of bentazone had to be calibrated to obtain a good description of its field persistence. The calibrated half-life was 2.5 times shorter than the half-life derived from the laboratory studies. Concentrations of bentazone in drain water and groundwater were described reasonably well by PEARL. Although measured and simulated persistence of imidacloprid in soil corresponded well, the bulk of the imidacloprid movement was overestimated by PEARL. However, imidacloprid concentrations in drain water were underestimated. In spite of the extensive calibration of water flow and tracer movement, the behaviour of the moderately sorbing pesticide imidacloprid could not be simulated. This indicates that the convection-dispersion equation cannot be used for accurate simulation of pesticide transport in cracking clay soils (even if extremely long dispersion length is used). Comparison of the model results from a poorly sorbed chemical (bentazone) and a moderately sorbed chemical (imidacloprid) were useful in defining the limitations of using a chromatographic flow model to simulate the effects of preferential flow.

Influence of dispersion length on leaching calculated with PEARL, PELMO and PRZM for FOCUS groundwater scenarios.

Harmonisation of the assessment of pesticide leaching to groundwater for EU registration is desirable to minimise confusion in the decision-making process at EU level. Recently, the FOCUS groundwater scenarios have been developed for three chromatographic models (PEARL, PELMO and PRZM) to increase this harmonisation. This study investigates the role of dispersion parameterisation in explaining the cause of the differences in pesticide leaching calculated by these models. PEARL describes dispersion via a physical parameter, ie the dispersion length. PELMO and PRZM simulate dispersion via a numerical procedure which generates an effective dispersion length equal to 0.5 times the thickness of the numerical compartments. The hypothesis was tested that the difference in the dispersion length input parameter (ie 5 cm for PEARL and about 2.5 cm for PELMO and PRZM) is a major cause of the difference in calculated leaching. It was tested whether results of PEARL calculations with a dispersion length of 2.5 cm corresponded much better to results of PELMO or PRZM than results of PEARL calculations with a dispersion length of 5 cm. This was done by calculations for one substance and all nine FOCUS scenarios and by calculations for a range of substances and two FOCUS scenarios (Châteaudun and Sevilla). All calculations were for winter wheat and an application at 1 day before emergence. Both tests showed that reduction of the dispersion length from 5 to 2.5 cm in PEARL led to a much better correspondence between PEARL and either PELMO or PRZM. Hence the hypothesis was supported. It is likely that harmonisation of the dispersion length in the FOCUS groundwater scenarios would reduce the differences in calculated leaching between PEARL and PELMO or PRZM considerably for part of these scenarios.

Preferential flow of bromide, bentazon, and imidacloprid in a Dutch clay soil.

Leaching to ground water and tile drains are important parts of the environmental assessment of pesticides. The aims of the present study were to (i) assess the significance of preferential flow for pesticide leaching under realistic worst-case conditions for Dutch agriculture (soil profile with thick clay layer and high rainfall) and (ii) collect a high-quality data set that is suitable for testing pesticide leaching models. The movement of water, bromide, and the pesticides bentazon [3-isopropyl-1H-2, 1,3-benzothiadiazine-4(3H)-one-2,2-dioxide] and imidacloprid [1-[(6-chloro-3-pyridinyl)-methyl]-N-nitro-2-imidazolidinimine] was monitored in a clay soil for about 1 yr. The 1.2-ha field was located in the central part of the Netherlands (51 degrees 53' N, 5 degrees 43' E). The soil was a Eutric Fluvisol cropped with winter wheat (Triticum aestivum L.). Tile drains were present at a 0.8- to 0.9-m depth and the ground water level fluctuated between a 0.5- and 2-m depth. All chemicals were applied in spring. None of the soil concentration profiles showed bimodal concentration distributions. However, for each substance the highest concentration in drain water was found in the first drainage event after its application, which indicates preferential flow. This preferential flow is probably caused by permanent macropores that were present in the 0.3- to 1.0-m layer. At the time of the first drainage event, the drain water concentration of each substance was about an order of magnitude higher than its ground water concentration. Thus, the flux concentrations in drain water proved to be a more sensitive detector of preferential flow than the resident concentrations in the soil profile and the ground water.