Drift from the Use of Hand-Held Knapsack Pesticide Sprayers in

Offsite pesticide losses in tropical mountainous regions have been little studied. One example is measuring pesticide drift soil deposition, which can...
1 downloads 6 Views 573KB Size
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