Subscriber access provided by NEW YORK MED COLL
Article
Drift from the use of handheld knapsack pesticide sprayers in Boyacá (Colombian Andes) Glenda García-Santos, Giuseppe Feola, David Nuyttens, and Jaime Diaz J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.5b03772 • Publication Date (Web): 19 Oct 2015 Downloaded from http://pubs.acs.org on October 25, 2015
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 34
Journal of Agricultural and Food Chemistry
1
Drift from the use of handheld knapsack pesticide sprayers in Boyacá
2
(Colombian Andes)
3 4
Glenda García-Santos1,2 *, Giuseppe Feola1,3, David Nuyttens 4, Jaime Diaz5
5 6
1
Department of Geography, University of Zurich, Switzerland
7
2
Department of Geography, Alpen-Adria-University Klagenfurt, Austria (current
8
address)
9
3
Department of Geography and Environmental Science, University of Reading,
10
United Kingdom
11
4
12
Agricultural and Fisheries Research (ILVO), Belgium
13
5
14
Boyacá, Colombia
Agricultural Engineering, Technology and Food Science Unit, Institute for
Department of Sanitary and Environmental Engineering, Universidad de
15 16
* Corresponding author at: Department of Geography, Alpen-Adria-University
17
Klagenfurt, Universitaetstrasse 65-67, 9020, Austria. E-mail address:
18
[email protected] (G. García-Santos)
19 20
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 34
21 22
Abstract
23
Offsite pesticide losses in tropical mountainous regions have been little studied.
24
One example is measuring pesticide drift soil deposition, which can support
25
pesticide risk assessment for surface water, soil, bystanders, off target plants
26
and fauna. This is considered a serious gap, given the evidence of pesticide-
27
related poisoning in those regions. Empirical data of drift deposition of a
28
pesticide surrogate, Uranine tracer, within one of the highest potato producing
29
regions in Colombia, characterized by small plots and mountain orography, is
30
presented. High drift values encountered in our study reflect the actual spray
31
conditions using handled knapsack sprayers. Comparison between measured
32
and predicted drift values using three existing empirical equations showed
33
important underestimation. However, after their optimization based on
34
measured drift information, the equations showed a strong predictive power for
35
this study area and the study conditions. The most suitable curve to assess
36
mean relative drift was the IMAG calculator after optimization.
37 38
Keywords: drift curve, knapsack sprayer, “Manual Técnico Andino”, developing
39
countries, mountain region, Andean region, tropical soils, potato production
2 ACS Paragon Plus Environment
Page 3 of 34
Journal of Agricultural and Food Chemistry
40 41
Introduction
42
Tropical mountainous regions in developing countries are often neglected in
43
research and policy, but represent key areas to be considered, if a sustainable
44
agricultural and rural development is to be promoted.1 These mountain
45
ecosystems are fragile sources of ecosystem services to communities living in
46
the lowlands (e.g. water resources) and are often directly threatened by human
47
activities, e.g. land use competition between forest and agriculture.2-3 Most
48
importantly, despite often being geographically marginal, tropical mountainous
49
regions are relevant both in terms of contribution to the agricultural production
50
and for sustaining a significant part of the population’s livelihoods, despite low
51
productivity.4 This is observed, for example, in the Andean countries (Peru,
52
Ecuador, Bolivia and Colombia). In the highlands of Colombia, for instance, the
53
mountainous department of Boyacá contributes to about 26% of the national
54
potato,5 which relies mainly on smallholders (more than the 95% of the
55
workforce), who are tenants of more than the 56% of the potato cultivated land,
56
and supply 45% of the regional production .5
57
Due to their relative marginality and low degree of mechanization, agricultural
58
research has only had a limited impact on the communities and crops in those
59
regions, which has also been recognised in the context of the so called Green
60
Revolution.6 One of the issues yet to be properly addressed is that of measuring
61
pesticide drift in typical field conditions using handheld application techniques,
62
knapsack sprayers, and developing accurate pesticide drift prediction tools,
63
which can support pesticide risk assessment. A recent review on drift mass
64
balance in pesticide application omitted the handheld techniques and focused
3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 34
65
on boom sprayers, possibly due to the lack of knapsack sprayer studies. In this
66
regard, few drift studies based on knapsack technique were found by the
67
authors .8-17 This is considered a serious gap, given the evidence of significant
68
presence of pesticide-related environmental and health risks e.g. in Ecuador,18
69
Bolivia19 and Colombia .14, 20-23
70
Pesticide risk after boom sprayers is commonly evaluated by empirical or
71
probabilistic approaches.24 Several empirical curves or models have been
72
developed for this technique in, among others, Germany,25 The Netherlands,26
73
and Belgium.27,
74
regulatory bodies in the respective countries.
75
The Ganzelmeier et al.25 and Holterman and van de Zande26 equations
76
estimate drift deposition as a function of the distance from the field border and
77
crop type, whereas Nuyttens et al.28,
78
equation for a reference spray boom application including drift distance as well
79
as meteorological conditions (wind speed, temperature and humidity).
80
In Colombia, drift is included in the assessment of risk to terrestrial ecosystems,
81
specifically to earthworms, as part of the control and registration process
82
applicable to agrochemicals. Drift is evaluated through table-based values
83
(worst case scenario) given by the Andean technical handbook/manual,
84
“Manual técnico Andino”, in function of distance from the applied field border.30
85
The origin of the table-based values is not quoted and is not known by the
86
authors. However, values seem to be derived by the Ganzelmeier et al.
87
equation 25.
88
One of the major advantages of empirical models compared with mechanistic
89
models is that they are easy to use and limited data are needed,31 which makes
28, 29
Most of these curves were designed to be used by
29
developed a non-linear drift prediction
4 ACS Paragon Plus Environment
Page 5 of 34
Journal of Agricultural and Food Chemistry
90
them especially appealing for pesticide risk assessment purposes in developing
91
countries, where low data availability is often an issue. A major disadvantage of
92
empirical models is that they produce single drift values and the evaluation of
93
probability is usually insufficiently considered.24 In addition, evidence exists that
94
empirical models do not perform well when used in contexts different from those
95
they were developed in.32 It is therefore sensible to verify whether existing table-
96
based values and curves are appropriate for use elsewhere.
97
In the Andean region, pesticides are mainly applied using knapsack sprayers,
98
and orography and micro-meteorology are complex. 90% of the farmers are
99
smallholders, who mainly apply pesticides using knapsack sprayers with old
100
and sometimes obstructed nozzles .16 In this region, there is poor maintenance
101
of equipment, no monitoring of meteorological conditions and no respect for
102
buffer zones (observed by the authors).
103
It is expected that the use of existing drift curves, designed under European
104
conditions, might produce an important discrepancy between the actual and the
105
expected pesticide drift deposition. In addition, the risk for rural residents might
106
also be inaccurately estimated .16,
107
assessment based on these drift models might be seriously flawed. Local scale
108
studies will allow the creation of datasets from which the future distribution
109
functions can be interpolated
110
pesticide risk more realistically.
111
The framework of this study is the project entitled “Reducing human health and
112
environmental risks from pesticide use: Integrating decision-making with spatial
113
risk assessment models: The case of Vereda la Hoya, Colombia”.34 Under the
114
umbrella of the project our working hypothesis was that there would be
24
33
As a consequence, pesticide risk
and used to assess environmental and human
5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 34
115
differences between measured and predicted drift depositions using existing
116
drift curves due to important differences in environmental and meteorological
117
conditions encountered in the studied area (agricultural and anthropogenic
118
effects kept constant) and differences in spray application practices. The
119
objectives of this study were i) to present for the first time an extended drift
120
dataset in marginal potato production located in the second biggest potato
121
production region of Colombia and ii) to assess the suitability of the existing
122
empirical drift worst case scenario curve suggested for the Andean region and
123
of the empirical mean drift curves used in Europe for risk assessment.
124
Materials and methods
125
Study area and sprayed plot. The study area is within the potato crop
126
producing region of Boyacá, the second highest potato production area after
127
Nariño in Colombia.60 The experimental catchment and plot selected to
128
measure off-target drift deposition on downwind soil is in the district La Hoya
129
community (Tunja, Boyacá). La Hoya basin covers an area of 840 ha in the
130
altitude range between 2600 and ~3000 m a.s.l. and mainly features small
131
agricultural plots ( 100 m) and outside of the sprayed plot.
152
153
Spraying technique. For all the 25 spray replications, a farmer from La Hoya
154
sprayed the plot using the conventional application technique in low-
155
mechanised potato production, i.e. a knapsack sprayer equipped with one high
156
discharge hollow cone nozzle at a height above the canopy between 0.2 – 0.35
157
m, with a solution of a green fluorescent tracer dye (Uranine) (further details in
158
the next section). The nozzle used in the study was selected from among the
159
most used in potato production by 10 local farmers. Droplet size spectrum
160
information, missing in García-Santos et al.,35 was measured by the Institute for
161
Agricultural and Fisheries Research (ILVO) in Belgium, using a Phase Doppler
162
Particle Analyser (PDPA) laser-based measuring set-up and protocol as
163
described by Nuyttens et al.38, 39 at 2.75 bar and at a distance of 0.40 m below
164
the nozzle. Other relevant measured sprayer related parameters as described
37
, 16 HAPs were
7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 34
165
in García-Santos et al.35 were distance of nozzle to edge of crop (~1 m), liquid
166
pressure (2.75 bar), average spray volume (530 L ha-1), application speed of
167
the farmer (100 min ha-1), tracer concentration (7 mg in a 20 L tank), number of
168
20 L tanks per trial (n= 1), application time (15 minutes), and nozzle diameter (1
169
mm; no manipulation was observed).
170
Sprayed solution. Uranine (sodium fluorescein, C20H10Na2O5) was used to
171
quantify drift by dry deposition on the horizontal and vertical collectors as
172
surrogate of agrochemicals. Though there might be an effect of formulation on
173
pesticide drift, this effect is limited compared with the effect of application
174
technique and environmental conditions. It is therefore generally accepted to
175
use artificial tracers instead of pesticides for spray drift studies. Uranine is
176
generally accepted as a tracer and was used before by many others.
177
tracer is highly soluble. If released into the atmosphere, Uranine has an
178
estimated vapour pressure of 4.0 10-14 mm Hg at 25 oC, which indicates
179
Uranine exist solely in the particulate-phase in the atmosphere and may be
180
removed by wet and dry deposition. Dry deposition on HAPs of Uranine is not
181
expected to volatilize based upon its low vapour pressure (Hazardous
182
Substances Data Bank, HSDB). This tracer is inexpensive, toxicologically safe
183
(HSDB) and has an extremely low detection limit of ~0.005 µg L-1. Uranine was
184
mixed with 20 L of water in a tank. Samples of the solution in the tank were
185
collected before spraying, to measure the initial tracer concentration per trial.
186
The HAPs were collected from the field and the vertical mesh immediately after
187
the trial, dried in an oven, placed in plastic bags and stored in a dark place until
188
their analysis. Details of the tracer extraction protocol, recovery (99%) and drift
189
calculation are described in detail in García-Santos et al.35
35, 42
This
8 ACS Paragon Plus Environment
Page 9 of 34
Journal of Agricultural and Food Chemistry
190
Meteorological conditions and border effect. Drift experiments were
191
conducted in 2008 and 2009 out of the rainy season, when farmers are more
192
likely to spray agrochemicals. The area has prevailing south east winds with an
193
average wind speed of 1.8 ± 1.39 m s-1 and a maximum of 7.6 m s-1 (1 year
194
meteorological
195
meteorological station, Davis Vantage Pro-2, was installed at 10 m distance
196
from the experimental plot to measure the meteorological variables during the
197
trials. Sensors such as a thermo-hygrometer (± 3% accuracy or error), a wind
198
speed sensor (± 5% accuracy), a rain gauge (± 4% accuracy) and a weather
199
vane (± 7o accuracy) were located in a mast 3 m above ground. The variables
200
considered in the study were ambient temperature (T, oC), relative humidity
201
(RH, %), dew point (Td, oC), vapour pressure deficit (VPD, kPa), wind speed (u,
202
m s-1) and wind direction (wd, o). All variables were registered every 1 minute
203
and averaged per trial (average trial duration 12 minutes). Relative drift
204
measurements at 7 and 11 m distances from the corner of the plot were studied
205
statistically (Kolmogorov-Smirnov Z test) to evaluate the existence of a possible
206
border effect on the results.
207
Spray drift deposition curves. We selected and tested three accepted drift
208
curve types. One of these equations is derived from the table base drift
209
information to assess risk to earthworms, suggested by the Andean technical
210
manual (“Manual técnico Andino”) and used for the Andean region (no
211
information on field datasets were found). The other two are used in different
212
regions of Europe (more details below). Mean values and 90th percentile values
213
for each drift distance and wind speed from our dataset (25 trials), were
214
compared to the predicted values. The goodness of fit of measured data was
measurements;
data
from
this
study).
An
automatic
9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 34
215
studied through the coefficient of determination (r2), predicting efficiency (EF) as
216
the Nash–Sutcliffe model efficiency coefficient43 and the root mean square error
217
(RMSE). EF and r2 close to 1 and low RMSE express a good agreement
218
between measured and predicted data.
219
The first equation tested predicts a worst case scenario of drift as a function of
220
distance and is used to estimate terrestrial risk for earthworms in a tiered
221
approach. We obtained this equation from the table-based values in the Andean
222
technical handbook/manual for pesticide registration and control in the Andean
223
region (Annex 7 E- Table 1 in Gaceta Oficial30). The origin of the table-based
224
values or reference to previous studies are not mentioned in the official
225
document. We observed that the derived equation from the table follows the
226
same mathematical expression as that developed by Ganzelmeier et al.25 in
227
Germany i.e. a single exponential function
228
%D = a x b
229
where %D is relative drift expressed as the percentage of the applied dose and
230
x (m) is the distance from the border to a point downwind from the treated field,
231
and a and b are 3.7839 and -1.1219, respectively. Note that the equation in
232
Ganzelmeier et al.25 had the following coefficients 2.7593 (a) and -0.9778 (b).
233
The second equation tested, IMAG drift calculator (v 1.1), was developed in the
234
Netherlands26 for conventional boom sprayers in function of distance as the
235
sum of two exponential functions,
236
%D = a e-xb + c e-xd
237
where x (m) is the downwind distance from the last nozzle to field edge and the
238
values of a (%), b (m-1), c (%) and d (m-1) depend on the crop and application
(1)
(2)
10 ACS Paragon Plus Environment
Page 11 of 34
Journal of Agricultural and Food Chemistry
239
technique. In case of an extensive crop (0.5 m height) and conventional
240
standard horizontal boom sprayers (1 m height), a, b, c, d are 114, 1.74, 1.29
241
and 0.139, respectively. For bare soil ( 0.5 m, so it might be a combination of other factors, like for
411
example how well the lance is directed downwards (not investigated here).
412
Some recommendations to decrease drift while using a knapsack sprayer point
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 34
413
out to lower the sprayer tip and use of drift shield attached to the lance of the
414
knapsack sprayer,17 not observed in the study area.
415
We found that meteorological variables alone explained 45% of the measured
416
drift variance. This low prediction power is explained by the complexity of drift
417
as measured under field conditions, also observed by Snelder et al..13 Though
418
high correlations were found by Nuyttens et al.28,
419
this is explained because of the use of wind tunnels. Thereby, drift potential is
420
found to increase with wind speed e.g. double spray drift deposition is found
421
with an increase from 3 m s-1 to 5 m s-1.53, 54 In our study and despite the low
422
prediction power of wind speed, very high wind speed increased the amount of
423
drift (Fig. 3). The studied wind conditions do not follow recommendations by
424
good agricultural practices55 (1 – 1.8 m s-1) to minimize the potential for drift, but
425
rather represent higher wind velocities, actual conditions, in which farmers
426
eventually or accidentally might spray. This might happen often in marginal
427
mountainous regions.36 Low relative humidity and higher temperatures increase
428
the potential of drift.44,
429
humidity data do not represent a sufficiently wide range to make a sound
430
analysis. Nevertheless, we observed that trials with a higher temperature
431
combined with low relative humidity tended to yield higher drift.
432
To conclude, factors like the use of high pressure cone nozzles producing a
433
high volume of small droplets, high distance between the nozzle and crop
434
canopy and high wind speed were found to affect drift positively. Other factors
435
like farmers’ mode of spraying were not considered in our study i.e. only one
436
farmer sprayed, but it is contemplated the use of different farmers to add
437
variability to the presented results as also observed in rice production studies.13
54, 56
29, 44
and Holterman et al.,53
However, our measured temperature and relative
18 ACS Paragon Plus Environment
Page 19 of 34
Journal of Agricultural and Food Chemistry
438
Performance of the drift curves. Drift deposition on soil was estimated
439
according to the table-based information provided by the Andean technical
440
manual (“Manual técnico Andino”). The worst drift scenario at 1 m distance from
441
the field edge underestimated the measured data (90th percentile) by 74%
442
(derived from Table 2, Andean technical manual). These results need to be
443
treated with caution because table-based drift values might correspond to boom
444
spray methods instead of knapsack spraying. Information on spray technique,
445
original data, geographical information of the study, theoretical background of
446
the table-based drift values are not mentioned in the Andean technical manual.
447
As mentioned before, we observed that the table information follows the
448
Ganzelmeier et al.25 equation. In this respect, in the Catania plain (Italy)57 as
449
well as in Adana-Turkey58 a low performance of the equation as applied to
450
agricultural areas was found. However, its prediction power significantly
451
improved after optimization of the coefficients i.e. overestimation of measured
452
drift of 3%. An optimized Andean curve as in Ganzelmeier et al.25, may be used
453
to assess risk characterisation using conservative data (realistic worst case
454
condition) under good agricultural practices (low wind speed, < 2.5 m s-1) for the
455
highlands of marginal potato production in Boyacá using knapsack sprayers.
456
The IMAG calculator and Belgian curve44 had low performance. Low
457
performance of the IMAG calculator was also found by others.58 One could
458
explain this because they were designed for boom sprayers. Another source of
459
error might be due to the limited trial number (n = 25), although, this number
460
greatly varies from 2,59 1235 to hundreds.25 We found that optimized curves by
461
IMAG calculator and by Nuyttens et al.44 performed well. The best prediction
19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 34
462
came from the optimized IMAG calculator in function of distance-only with an
463
overestimation of measured drift of 9% at the first meter from the border edge.
464
We found that the inclusion of meteorological variables did not result in a higher
465
predictive power of the measured drift values. This might also reflect the lower
466
correlations encountered before. However, the curve in function of wind speed
467
and other meteorological conditions by Nuyttens et al.44 performed well and
468
therefore it may be used to assess actual drift deposition under different wind
469
speeds up to 6.5 m s-1 (measured in our study) and drift risks.
470
The results showed that the existing empirical equations derived from boom
471
sprayer technique were insufficient to predict handled spray drift for tropical
472
mountainous environments. Existing equations generally underestimated the
473
spray drift and hence pesticide risk would also be underestimated. However, the
474
existing equations could be optimised with application of relevant newly
475
collected data from the local region i.e. the optimized IMAG calculator was
476
suitable to predict drift deposition and the optimized curve derived from the
477
Andean technical manual is suggested to assess the highest tier for risk
478
assessment. In that regard, these equations with revised coefficients may have
479
greater predictive capacity for broader tropical mountainous regions, although
480
this would need to be confirmed with further study.
481
Acknowledgement. This research was funded by the Swiss National Science
482
Foundation (SNSF) through the project number 110807. We thank Prof. C.R.
483
Binder and Prof. H. Barnard for their revision of a previous version of the
484
manuscript; the students D. Scheiben, S. Karrer, P. Muñoz, F. Pui and the
485
farmer who collaborated during the drift trials. Writing of the manuscript was
486
possible thanks to A. Santos-Martin and M. Boudet-Garcia.
20 ACS Paragon Plus Environment
Page 21 of 34
Journal of Agricultural and Food Chemistry
487
Supporting information description.
488
Table 1 includes the Kolmogorov-Smirnov Z test of the difference of relative drift
489
expressed as a percentage of the applied dose (%) at 7 m and 11 m distance
490
from the field corner. Exact significance: *0.01 level.
491
References
492
(1)
Parthenon Publishing, 1997, New York and Carnforth.
493 494
Messerli, B.; Ives, J.D. Eds. Mountains of the world: a global priority.
(2)
Etter, A.; Villa, L.A. Andean Forests and Farming Systems in part of the
495
Eastern Cordillera (Colombia). Mountain Research and Development.
496
2000, 20, 236-245
497
(3)
Etter, A.; McAlpine, C.; Wilson, K.; Phinn, S.; Possingham, H. Regional
498
patterns of agricultural land use and deforestation in Colombia. Agriculture,
499
Ecosystems and Environment. 2006, 114, 369-386
500
(4)
Andes. Rimisp-Latin American Center for Rural Development. 2007.
501 502
Aldana, U. The importance of agriculture in isolated areas in the Peruvian
(5)
Gobernación de Boyacá Direccion de Desarrollo Agropecuario. Acuerdo
503
de competitividad de la cadena agroalimentaria de la papa, para el
504
Departamento de Boyacá. 2004, Tunja, Colombia.
505
(6)
Century. Cornell University Press, Ithaca, New York, USA, 1998.
506 507
(7)
510
Jensen, P. K.; Olesen, M. H. Spray mass balance in pesticide applications. A review. Crop protection. 2014, 61, 23-31.
508 509
Conway, G. The Doubly Green Revolution: Food for All in the Twenty-First
(8)
FAO. Guidelines on Good Practice for Ground Application of Pesticides. Ed. Food and Agriculture Organization of the United Nations. 2001, Rome.
21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
511
(9)
Page 22 of 34
Ohayo-Mitoko, G. J. A. Occupational pesticide exposure among Kenyan
512
agricultural workers. An epidemiological and public health perspective.
513
PhD dissertation. 1997, Wageningen University. The Netherlands.
514
(10) Sibanda, T.; Dobson, H. M.; Cooper, J. F.; Manyangarirwa, W.; Chiimba,
515
W. Pest management challenges for smallholder vegetable farmers in
516
Zimbabwe. Crop Protection. 2000. 19: 807-815.
517
(11) Miller, A.; Bellinder, R. Herbicide application using a knapsack sprayer.
518
New Delhi, India: Facilitation Unit, Rice-Wheat Consortium for the Indo-
519
Gangetic Plains. 2001.
520 521 522
(12) Bateman, R. The use of narrow-angle cone nozzles to spray cocoa pods and other slender biological target. Crop Protection. 2004, 23, 989–999. (13) Snelder, D.J.; Masipiquen, M.D.; de Snoo, GR. Risk assessment of
523
pesticide usage
by
smallholder
farmers
in
524
(Philippines). Crop Protection. 2008, 27, 747-762.
the Cagayan
Valley
525
(14) García-Santos, G.; Keller-Forrer, K. Avoidance behaviour of Eisenia Fetida
526
to Carbofuran, Chlorpyrifos, Mancozeb and Metamidophos in the
527
highlands of Colombia. Chemosphere. 2011, 84, 651-656.
528
(15) Juraske, R.; Vivas C.S.; Velasquez A.E.; García-Santos, G.; Moreno M.B.;
529
Gomez J.D.; Binder C.R.; Hellweg, S.; Dallos J.A. Pesticide uptake in
530
potatoes: model and field experiments. Environ Sci Technol. 2011, 45:651-
531
657.
532
(16) Lesmes, C.F.; García-Santos, G.; Leuenberger, F.; Nuyttens, D.; Binder,
533
C.R. Dermal Exposure Assessment of Pesticide Use: The Case of
534
Sprayers in Potato Farms in the Colombian Highlands. Science of the
535
Total Environment. 2012, 430, 202–208.
22 ACS Paragon Plus Environment
Page 23 of 34
Journal of Agricultural and Food Chemistry
536
(17) Franke, A.C.; Kempenaar, C.; Holterman H.J.; van de Zande J.C. Spray
537
Drift from Knapsack Sprayers. Plant Research International B.V.,
538
Wageningen. 2010, 658, 26.
539
(18) Crissman, C.C.; Cole, D.C.; Carpio, F. Pesticide Use and Farm Worker
540
Health in Ecuadorian Potato Production. American Journal of Agricultural
541
Economics. 1994, 76, 593-597.
542
(19) Jors, E.; Morant, R.; Aguilar, G.; Huici, O.; Lander, F.; Baelum, J.;
543
Konradsen, F. Occupational pesticide intoxications among farmers in
544
Bolivia: a cross-sectional study. Environmental Health. 2006, 5, 10.
545
(20) Cardenas, O.; Silva, E.; Morales, L.; Ortiz, J. Estudio epidemiológico de
546
exposición a plaguicidas organofosforados y carbamatos en siete
547
departamentos colombianos, 1998-2001. Biomedica. 2005, 25: 170-180.
548
(21) Ospina, J.M.; Manrique, F.G.; Ariza, N.E. Salud, ambiente y trabajo en
549
poblaciones vulnerables: los cultivadores de papa en el centro de Boyacá.
550
Revista de la Facultad Nacional de Salud Pública. 2008, 26: 143-152.
551
(22) Feola, G.; Binder, C.R. Why don’t
pesticide applicators
protect
552
themselves? Exploring the use of Personal Protective Equipment among
553
Colombian smallholders. International Journal of Occupational and
554
Environmental Health. 2010, 16: 11-23.
555
(23) Olivero-Verbel, J. Colombia: Environmental health issues. In Encyclopedia
556
of Environmental Health, edited by J.O. Nriagu, Elsevier, Burlington. 2011,
557
740-754.
558 559
(24) Balderacchi, M. and Trevisan, M. Comments on pesticide risk assessment by the revision of Directive EU 91/414, 2010.
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 24 of 34
560
(25) Ganzelmeier H.; Rautmann, D.; Spangenberg, R.; Streloke, M. Studies on
561
the Spray Drift of Plant Protection Products. Mitteilungen aus der
562
Biologischen Bundesanstalt für Land- und Forstwirtschaft Berlin-Dahlem,
563
Heft 305 Blackewll Wissenschafts-Verlag GMBH 1995, Berlin/Wein.
564 565 566
(26) Holterman, H.J.; van de Zande, J.C. IMAG Draft report, IMAG Drift Calculator v 1.1, User Manual. 2003. http://idc.holsoft.nl. (27) Nuyttens, D. Drift from field crop sprayers: The influence of spray
567
application
technology
determined
using
indirect
and
direct
drift
568
assessment means. Katholieke Universiteit Leuven. PhD thesis. 2007, 293
569
pp. ISBN 978-90-8826-039-1
570
(28) Nuyttens D.; De Schampheleire M.; Baetens K.; Sonck B. The influence of
571
operator controlled variables on spray drift from field crop sprayers.
572
Transactions of the ASABE. 2007, 50:1129-1140.
573
(29) Nuyttens D.; De Schampheleire M.; Baetens K.; Brusselman E.; Dekeyser
574
D.; Verboven P.; Drift from Field Crop Sprayers using an Integrated
575
Approach: Results of a 5 Year Study. Transactions of the ASABE. 2011,
576
54: 403-408
577
(30) Gaceta Oficial. Manual técnico andino para el registro y control de
578
plaguicidas químicos de uso agrícola. Ed. Secretaría General de la
579
Comunidad Andina, Perú. 2002, 169.
580
(31) McKinlay, R.; Plant, J. A.; Bell, J.N.B.; Voulvoulis, N. Calculating human
581
exposure to endocrine disrupting pesticides via agricultural and non-
582
agricultural exposure routes. Science of the Total Environment. 2008, 398,
583
1-12.
24 ACS Paragon Plus Environment
Page 25 of 34
Journal of Agricultural and Food Chemistry
584 585
(32) Meli, S.M.; Renda, A.; Nicelli, M.; Capri, E. Studies on pesticides spray drift in a Mediterranean citrus area. Agronomie. 2003, 23, 667-672.
586
(33) Lubin, J.; Giglierano, J.S.; Colt, C.; Wolter, N.; Bekiroglu, D. Proximity to
587
crops and residential exposure to agricultural herbicides in Iowa, Environ
588
Health Perspect. 2006, 114, 893–897.
589
(34) Binder, C.R., Garcia-Santos, G.,.Andreoli, R., Feola, G., Wittensoeldner,
590
M., Yang, J., Diaz, J. Simulating human and environmental exposure from
591
handheld knapsack pesticide application: WetSpa-Pest, Can integrative,
592
spatially explicit modeling approach. This issue.
593
(35) García-Santos, G.; Scheiben, D.; Binder, C.R. The weight method: a new
594
screening method for estimating pesticide deposition from knapsack
595
sprayers in developing countries. Chemosphere. 2011, 82, 1571-1577.
596
(36) Feola, G.; Binder, C.R. Identifying and investigating pesticide application
597
types to promote a more sustainable pesticide use. The case of
598
smallholders in Boyacá, Colombia. Crop Protection. 2010, 29: 612-622.
599
(37) De Schampheleire, M.; Baetens, K.; Nuyttens, D.; Spanoghe P. Spray drift
600
measurements to evaluate Belgian drift mitigation measures in field crops.
601
Crop Protection. 2008, 27: 577-589.
602
(38) Nuyttens D.; Baetens K.; De Schampheleire M.; Sonck B. Effect of nozzle
603
type, size and pressure on spray droplet characteristics. Biosystems
604
Engineering. 2007, 97: 333-345.
605
(39) Nuyttens D.; De Schampheleire M.; Verboven P.; Brusselman E.;
606
Dekeyser D. Droplet size-velocity characteristics of agricultural sprays.
607
Transactions of the ASABE. 2009, 52: 1471-1480.
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 34
608
(40) Southcombe, E.S.E.; Miller, P.C.H.; Ganzelmeier, H.; van de Zande, J.C.;
609
Miralles, A.; Hewitt, A.J. The international (BCPC) spray classification
610
system including a drift potential factor. The Brighton Crop Protection
611
Conference – Weeds. 1997, 371-380.
612 613
(41) Lefebvre,
A.H.
Atomization
and
sprays.
Hemisphere
Publishing
Corporation, 1989, New York.
614
(42) Gil, Y.; Sinfort, C. Emission of pesticides to the air during sprayer
615
application: a bibliographic review. Atmospheric Environment. 2005, 39,
616
5183–5193.
617
(43) Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual
618
models: Part 1. A discussion of principles, J. Hydrol. 1970, 10(3), 282-290.
619
(44) Nuyttens, D.; De Schampheleire, M.; Steurbaut, W.; Baetens, K.;
620
Verboven, P.; Nicolai, B.; Ramon, H.; Sonck, B. Experimental study of
621
factors influencing the risk of drift from field sprayers, Part 1:
622
Meteorological conditions. Aspects of Applied Biology, International
623
Advances in Pesticide Application. 2006, 77: 321-329.
624
(45) Vischetti, C.; Cardinali, A.; Monaci, E.; Nicelli, M.; Ferrari, F.; Trevisan, M.;
625
Capri, E. Measures to reduce pesticide spray drift in a small aquatic
626
ecosystem in vineyard estate. Science of the Total Environment. 2008,
627
389, 497-502.
628
(46) Mosquera-Vivas, C.S.; García-Santos, G.; Celis-Ossa, R. E.; Hellweg, S.
629
H.; Juraske, R.; González-Murillo, C. A.; Guerrero-Dallos, J. A. Adsorption,
630
desorption and degradation behaviour of the pesticides pyrimethanil,
631
dimethomorph and chlorpyrifos through a tropical Colombian agricultural
632
soil profile. Not published
26 ACS Paragon Plus Environment
Page 27 of 34
Journal of Agricultural and Food Chemistry
633
(47) Mosquera-Vivas, C. S.; Martinez-Cordon, M. J.; García-Santos, G.;
634
Guerrero-Dallos, J. A. Adsorption-desorption studies of tebuconazole in
635
tropical Colombian agricultural soils. Not published
636
(48) De Snoo, G.R.; de Wit, P.J. Pesticide drift from knapsack sprayers to
637
ditches and ditch banks. Proceedings of the Brighton Crop Protection
638
Conference - Weeds, Farnham UK: British Crop Protection Council. 1993,
639
879-884.
640
(49) De Snoo, G.R.; de Wit, P.J. Buffer zones for reducing pesticide drift to
641
ditches and risk to aquatic organisms. Ecotoxicology and Environmental
642
Safety. 1998, 41: 112-118.
643 644
(50) Gilbert, A. J. Local environmental risk assesment for pesticides (LERAP) in the UK. Aspects of Applied Biology. 2000, 57, 83-90.
645
(51) Baetens, K.; Ho, Q. T.; Nuyttens, D.; De Schampheleire, M.; Melese
646
Endalew, A.; Hertog, M.; Nicolaï, B., Ramon, H., Verboven, P. A Validated
647
2-D diffusion-advection model for prediction of drift from ground boom
648
sprayers. Atmospheric environment. 2009, 43: 1674-1682.
649
(52) Heidary, M. Al.; Douzals, J.P.; Sinfort, C.; Vallet, A. Influence of spray
650
characteristics on potential spray drift of field crop sprayers: A literature
651
review. Crop protection. 2014, 1- 11.
652
(53) Holterman, H.J., Van De Zande, J.C., Porskamp, H.A.J., Huijsmans, J.F.M.
653
Modelling spray drift from boom sprayers. Computers and electronics in
654
agriculture 1997, 19 (1), 1-22.
655
(54) Thistle, H.W.; Teske, M.E.; Reardon, R.C. Weather effects on drift
656
meteorological factors and spray drift: An overview. Proceedings of the
27 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 28 of 34
657
North American Conference on Pesticide Spray Drift Management March
658
29 - April 1, Portland, Maine, 1998, 64-74.
659
(55) MAVDT. Guias ambientales para el subsector de plaguicidas. 2003, 104 p.
660 661
(56) Luo, Y.; Miller, DR.; Yang X.; McManus, MA.; Krider, HM. Characteristics
662
of evaporation from water based bacterial pesticide droplet. Trans. Am.
663
Soc. Agric. Eng. 1994, 37, 1473-1479.
664
(57) Padovani, L.; Capri, E.; Trevisan, M. Landscape-level approach to assess
665
aquatic exposure via spray drift for pesticides: a case study in a
666
Mediterranean area. Environmental Science and Technology. 2004, 38,
667
3239-3246.
668
(58) Yarpuz-Bozdogan, N.; Bozdogan, A. M. Comparison of field and model
669
percentage drift using different types of hydraulic nozzles in pesticides
670
applications. Int. J. Sci. Tech. 2009, 6, 191 – 1996.
671
(59) Wolters, A.; Linnemann, V.; van de Zande, J. C.; Vereecken, H. Field
672
experiment on spray drift: Deposition and airborne drift during application
673
to a winter wheat crop. Science of The Total Environment. 2008, 405, 269-
674
277.
675
(60) Ministerio de Agricultura y Desarrollo Rural. La cadena de la papa en
676
Colombia. Una mirada global de su estructura y dinámica 1991-2005,
677
Ministerio de Agricultura y Desarrollo Rural, Bogotá, 2006, Colombia.
678
28 ACS Paragon Plus Environment
Page 29 of 34
Journal of Agricultural and Food Chemistry
679 Figure captions Figure 1. Experimental set up to measure drift deposition on soil from edge border up to 20 m from the sprayed field (right site) and in one vertical transect at 1 and 2 m height at 8 m distance from the field border (□: Horizontal High Absorbent Paper (HAP), ∆: Vertical HAPs). Figure 2. Relative drift deposition expressed as the percentage of the applied dose from the field border at 0 m up to 20 m distance after 25 trials (n= 134 samples) and mean relative drift calculated with the mean of the points (wind speed between 1 and 6.5 m s-1). Figure 3. Mean relative drift horizontal deposition from the field border at 0 m up to 20 m distance at different wind speed. Vertical lines denote one standard deviation from the mean. Inside figure: Airborne relative drift (%) at 1 and 2 m height from the ground and 8 m distance from the treated field as a function of wind speed (m s-1) measured at the nozzle height. Vertical lines represent one standard deviation from the mean.
29 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 30 of 34
Tables Table 1. Left: Measured mean relative drift for wind speed 1 -2.5 m s-1 and 3 – 6.5 m s-1 from the field edge to 20 m distance, the 90th percentile drift values and the mean values of the meteorological conditions during the spray trials (T: temperature; Td: dew point; RH: relative humidity (%); u: wind speed; VPD: vapour pressure deficit). Right: Calculated mean drift after the Andean curve, the IMAG calculator and the curve by Nuyttens et al. (2006). Table on the left shows drift values using the original equations and table on the right shows drift values after optimization of the coefficients. Below part shows performance of the curves as r2 (coefficient of determination), RMSE (root mean squared error) and EF (efficiency) and curve coefficients.
VPD (kPa)
1 - 2,5 m s -1 61,0 ± 22,2 11,3 -10,9 ± 4,5 6,8 ± 4,0 4,0 ± 2,3 2,7 ± 2,1 2,1 ± 2,2 1,6 ± 0,9 0,7 ± 0,8 0,6 ± 0,7 0,4 ± 0,2 0,2 ± 0,4 0,2 ± 0,4 0,2 ± 0,2 14,4 ± 0,2 13,1 ± 0,7 91,9 ± 4,8 2,5 ± 0,6 0,1 ± 0,1
82,6 11,3 14,6 10,8 5,8 4,7 4,3 2,4 1,5 1,2 0,6 0,6 0,6 0,4
3 - 6,5 m s -1 73,3 ± 26,9 38,2 ± 15,5 15,6 ± 9,0 9,4 ± 5,1 6,5 ± 4,0 4,8 ± 3,5 2,9 ± 2,1 3,7 ± 1,6 1,3 ± 1,1 2,1 ± 0,5 1,1 ± 0,8 0,6 ± 0,4 0,4 ± 0,3 0,5 ± 0,6 15,0 ± 0,1 12,3 -84,1 ± 1,3 5,0 ± 0,39 0,30
--
Relative drift (%) 90th precentile
Calculated drift (original curve coefficients)
Calculated drift (optimized curve coefficients)
Andean curve (90th IMAG percentile) (potato)
Andean curve (90th percentile)
IMAG
f(x)
f(x)
--15,0 8,6 6,2 4,9 4,1 3,4 2,8 2,5 2,2 1,8 1,6 1,4
--11,8 7,6 4,8 3,1 2,0 1,0 0,5 0,3 0,1 0,0 0,0 0,0
0,97 0,93 2,0
1,00 0,98 0,3
0,97 0,93 1,1
15 -0,8
29,00 -6,82 18,35 -0,44
-1,06 -1210,96 -71,60 3,74 12,57 68,77
IMAG Nuyttens et al. (bare soil) (2006)
f(x)
f(x)
f(x)
--3,8 1,7 1,1 0,8 0,6 0,5 0,4 0,3 0,3 0,2 0,2 0,1
--21,1 4,5 1,5 0,8 0,7 0,5 0,4 0,3 0,3 0,2 0,1 0,1
--6,9 2,4 1,3 1,0 0,8 0,7 0,5 0,4 0,4 0,2 0,2 0,1
2
0,92 -0,13 10,9
0,80 -0,04 0,05
0,91 0,60 4,5
0,98 0,13 6,2
a b c d e d
3,7839 -1,1219
114,00 1,74 1,29 0,139
25,00 1,50 1,54 0,133
-1,03 48,89 1,08 0,5 -0,45 -1,41
100,0 55,5 22,2 16,3 12,1 7,3 5,9 4,8 3,0 2,5 1,8 1,0 0,7 1,1 r EF RSME
f(x , u, T, Td, RH) 2,5 m s -1 5 m s -1 ----3,1 5,6 1,5 2,8 1,0 1,8 0,8 1,4 0,6 1,1 0,5 0,8 0,37 0,66 0,31 0,55 0,27 0,48 0,21 0,37 0,17 0,30 0,14 0,26 0,97 0,16 9,5
Optimized coefficients
0 0,5 1 2 3 4 5 6,5 8 9,5 11 14 17 20 o T ( C) o Td ( C) RH (%) -1 u (m s )
Measured drift Relative drift (%)
Performance
Meteorological variables
x (distance, m)
Relative drift (%) 90th precentile
Coefficients
Measured drift Relative drift (%)
Nuyttens et al. (2006) f(x , u, T, Td, RH) 2,5 m s -1 5 m s -1 ----9,6 18,9 4,6 9,1 3,0 5,9 2,2 4,3 1,7 3,4 1,3 2,6 1,1 2,1 0,9 1,7 0,8 1,5 0,6 1,1 0,5 0,9 0,4 0,8 0,96 0,93 1,0
30
ACS Paragon Plus Environment
Page 31 of 34
Journal of Agricultural and Food Chemistry
Figure graphics Figure 1. Experimental set up to measure drift deposition on soil from edge border up to 20 m from the sprayed field (right site) and in one vertical transect at 1 and 2 m height at 8 m distance from the field border (□: Horizontal High Absorbent Paper (HAP), ∆: Vertical HAPs).
31 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 32 of 34
Figure 2. Relative drift deposition expressed as the percentage of the applied dose from the field border at 0 m up to 20 m distance after 25 trials (n= 134 samples) and mean relative drift calculated with the mean of the points (wind speed between 1 and 6.5 m s-1).
32 ACS Paragon Plus Environment
Page 33 of 34
Journal of Agricultural and Food Chemistry
Figure 3. Mean relative drift horizontal deposition from the field border at 0 m up to 20 m distance at different wind speed. Vertical lines denote one standard deviation from the mean. Inside figure: Airborne relative drift (%) at 1 and 2 m height from the ground and 8 m distance from the treated field as a function of wind speed (m s-1) measured at the nozzle height. Vertical lines represent one standard deviation from the mean.
33 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 34 of 34
TOC Graphic
34 ACS Paragon Plus Environment