Degradation and Transport of the Chiral Herbicide S-Metolachlor at

Evaluating pesticide degradation and transport in the soil-surface water continuum remains challenging at the catchment scale. Here we investigated th...
1 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF MISSOURI ST LOUIS

Article

Degradation and transport of the chiral herbicide S-metolachlor at the catchment scale: combining observation scales and analytical approaches Marie Lefrancq, Sylvain Payraudeau, Benoît Guyot, Maurice Millet, and Gwenaël Imfeld Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b02297 • Publication Date (Web): 22 Oct 2017 Downloaded from http://pubs.acs.org on October 26, 2017

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

Environmental Science & Technology 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.

Subscriber access provided by UNIV OF MISSOURI ST LOUIS

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.

Environmental Science & Technology 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 32

Environmental Science & Technology

Infiltration Runoff

Soil Degradation Volatilization SM MESA Dissolved

PLOT

SM Particulate

MESA MOXA

SM Dissolved SM Plot export Particulate 8.1% ACS Paragon Plus Environment

MOXA Catchment export 11.2%

Environmental Science & Technology

1

Degradation and transport of the chiral herbicide S-

2

metolachlor at the catchment scale: combining

3

observation scales and analytical approaches

4

Lefrancq Marie.a,b, Payraudeau Sylvain.a*, Guyot Benoit.a, Millet Maurice.c, Imfeld

5

Gwenaël.a

6 7

a

8

University of Strasbourg, CNRS, ENGEES, 1 Rue Blessig, 67084 Strasbourg cedex, France

9

b

Laboratory of Hydrology and Geochemistry of Strasbourg (LHyGeS UMR 7517),

LETG-Angers (UMR CNRS 6554), University of Angers, 2 bd Lavoisier, 49045 Angers,

10

France

11

c

12

Strasbourg, CNRS, 1 rue Blessig, 67084 Strasbourg cedex, France

Atmospheric Physical Chemistry Department (ICPEES UMR 7515), University of

13 14

*Corresponding author: [email protected]

15

Laboratory of Hydrology and Geochemistry of Strasbourg (LHyGeS UMR 7517), Université

16

de Strasbourg, CNRS, ENGEES, 1 Rue Blessig, 67084 Strasbourg cedex, France. phone: + 33

17

(0)3 6885 0437 and fax: +33 (0)3 68 85 04 02

18 19

Manuscript for Environmental Science and Technology

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Environmental Science & Technology

20

Abstract

21

Evaluating pesticide degradation and transport in the soil-surface water continuum remains

22

challenging at the catchment scale. Here we investigated the dissipation of the chiral herbicide

23

S-metolachlor (SM) in soil in relation to its transport in runoff. Analyses of SM,

24

transformation products (TPs, i.e. MESA and MOXA), and enantiomers were combined to

25

determine SM degradation at plot and catchment scales. Assisted by modelling, we found that

26

the main dissipation pathways of SM at the plot scale were degradation (71%), volatilization

27

(5%), leaching (8%) and runoff (3%), while 13% of SM persisted in top-soil. This highlights

28

the relevance of degradation processes. TPs could trace the different discharge contributions:

29

MOXA prevailed in runoff water, whereas MESA was associated with slower flowpaths. At

30

the catchment outlet, 11% of SM applied was exported in dissolved or particulate phases or as

31

TPs (in SM mass equivalent). A single event one week after application exported 96% of SM,

32

which underlined the potential importance of severe rainfall on seasonal SM export.

33

Enantioselective degradation enriched SM in the R-enantiomer over longer periods and may

34

be associated with slower flowpaths. Altogether, combining observation scales and analytical

35

approaches enabled to quantify SM degradation and to identify how degradation controls SM

36

export at the catchment scale.

ACS Paragon Plus Environment

Environmental Science & Technology

37

Introduction

38

Understanding pesticide degradation in agricultural systems is essential as degradation

39

is the only process that can reduce mass transport of pesticides.1 Pesticide degradation is

40

controlled by a complex interplay between physico-chemical characteristics of molecules and

41

environmental conditions, including crop management practices, climatic forcing and water

42

flowpaths.2–4 Degradation processes are often estimated at agricultural plot scales,5,6 whereas

43

water management measures are defined on a larger scale, such as that of the hydrological

44

catchment.7 Degradation is rarely evaluated at the catchment scale8 because degradative

45

processes can hardly be distinguished from non-degradative ones, such as dilution or

46

sorption.1 As a result, relevant quantitative information is mostly lacking to understand factors

47

driving pesticide degradation in relation to transport at the catchment scale.

48 49

Chloroacetanilides belong to the top ten pesticide classes used worldwide.9 They are

50

used for pre-emergent control of annual weeds, mainly on corn, sugar beet and sunflower.

51

Chloroacetanilides are frequently detected in surface and ground water together with their

52

transformation products (TPs), notably ethane sulfonic acid (ESA) and oxanilic acids

53

(OXA).10–14 S-metolachlor (SM) is one of the most used chloroacetanilides and equates to

54

4.2% of the global use of pesticides.1 Similar to approximatively 25% of currently used

55

chemical pesticides, SM is a chiral compound, i.e. with two enantiomers which are non-

56

superimposable mirror images of each other.15 Enantiomers of chiral pesticides have identical

57

physico-chemical properties.16 However, enantiomers may interact selectively with molecular

58

receptors and thus differ in their biological efficacity, toxicity and environmental

59

behaviour.17–20 Enrichment of specific enantiomers in environmental compartments can

60

indicate biodegradation of chiral pesticides.21,22 The combination of different analytical

61

approaches, such as parent compound and TP quantification, and enantiomeric signatures can

62

help to identify and disentangle the main dissipation pathways of SM (degradation vs.

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Environmental Science & Technology

63

export).21,23 To the best of our knowledge, simultaneous evaluation of SM, its enantiomers

64

and major TPs in surface runoff at the catchment scale has not been reported yet.

65 66

Therefore, the purpose of this study was to evaluate the degradation of SM in soil in

67

relation to its export in runoff concomitantly with its major TPs, i.e., metolachlor-ESA

68

(MESA) and metolachlor-OXA (MOXA). SM was used as a proxy for moderately soluble

69

herbicides. Our specific objectives were to (i) identify the dissipation pathways of SM in

70

relation to hydrological conditions, and (ii) to quantify transport of SM and its main TPs at

71

two spatial scales: an experimental plot (77 m²) nested within a catchment (47.8 ha).

72

Measurements on both scales were strengthened with pesticide mass-balance modelling to

73

quantify the contribution of degradation and transport of SM in the soil-surface water

74

continuum. Beyond the quantification of SM, enantiomer fractionation and a mass-balance

75

analysis of MESA and MOXA production in the soil and export in runoff enabled to estimate

76

the contribution of degradation to SM dissipation at the catchment scale.

77 78

Material and methods

79

Study catchment. The study site is a 47.8-ha catchment that is part of the Zorn catchment of

80

750 km², situated in the lower Alsatian hills (Bas-Rhin, France) (Figure 1).24 Briefly, it is

81

characterized by silty-clay soils and a mean annual precipitation of 605 mm. Arable land

82

comprises 88% of the catchment area (of which, in 2012, 68% was corn, 16% winter wheat

83

and 4% sugar beet). Corn and sugar beet are sown between mid-March and end of May.

84

Natural drainage and a 5-cm diameter subsurface pipe draining an unknown part of the study

85

area discharged into a ditch located 102 m upstream of the catchment outlet (Figure 1). This

86

drainage continuously sustained the baseflow observed at the catchment outlet during the

87

study period (0.43 ± 0.15 L-1s-1).

88

ACS Paragon Plus Environment

Environmental Science & Technology

89

To isolate the contribution of subsurface and runoff flow, a 77 m² portion of a sugar

90

beet plot was studied in the upstream area of the catchment (plot is in dark orange in

91

Figure 1). This plot is considered representative of the corn and sugar beet area in the

92

catchment because soil types, slope and crop management techniques are spatially

93

homogeneous within the catchment.24 The plot was hydrologically isolated from runoff and

94

subsurface flow in the plowing layer from neighboring plots with a 60-cm high HDPE shield

95

to a depth of 30 cm below the ground surface, as saturation excess overland flow is unlikely

96

to occur24. The study was carried out from March 12 to July 10 (Day 0 to 120) at the plot

97

scale, and up to August 14, 2012 (Day 155) at the catchment scale.

98 99

SM characteristics and applications. SM is nonionic and moderately soluble in water,

100

whereas MESA and MOXA have a higher polarity and solubility and lower Kd values as

101

reported in the Supporting Information (Table S1).25,26 Racemic metolachlor has been applied

102

since the early 1990s at the study site, and, since 1998-2000, it has been replaced by SM,

103

which is enriched in the herbicidally active S-enantiomer.27,28 A total of 10.95 kg of SM was

104

applied in 9 plots in the catchment during the study period (Figure 1), with the percentage of

105

total SM by formulation of Mercantor Gold® (Day 31, 4.9% and Day 50, 4.9%), Camix® (Day

106

63, 78.9%) and Dual Gold® (Day 62, 11.3%). Only Mercantor Gold® was applied twice (Day

107

31 and Day 50) in the experimental plot. On Day 31, 16 petri dishes were used to quantify SM

108

deposition on the experimental plot (Annex S1), which ranged from 192 to 336 g ha-1 for 288

109

g ha-1 reported by the farmer. This highlights the spatial variability of SM deposition at the

110

plot scale. Drift and rainfall depositions of metolachlor were previously measured during 64

111

days for similar rainfall and higher application rate.29 This may represent 0.7 µm) to mass equivalent

164

load of SM (MELSM) was evaluated as a percentage of total loads of SM (illustrated here for

165

MESA [%]).

166

 MWSM  load (MESA) ×    MWMESA  × 100 %MESA = MELSM

(2)

167

Enantiomer fraction. The enantiomer fraction (EF [-]) was used to describe enantiomer

168

composition of SM:31  =

    +  

(3)

169

The EF of S-metolachlor in the commercial formulations is approximately 88% S isomer and

170

12% R isomer32, and is obtained through a selective and likely identical synthetic

171

manufacturing process by Syngenta.32–34 The EF of the commercial product Mercantor Gold®

172

applied in the studied plot was 0.866 ± 0.008 (n = 8), which is consistent with previous

173

studies for Dual II®, Dual magnum®32,33 or for Mercantor Gold®23, indicating an EF value

174

ranging between 0.872-0.882. Shift in the enantiomer fractions was calculated as follow: ∆ =    −      

175 176

(4)

A positive ∆EF thus indicates S-enantiomer enrichment, and, inversely, a negative one indicates R-enantiomer enrichment.

177 178

Pesticide mass balance model. A parsimonious model35,36 was developed to estimate the

179

different terms of the SM mass balance (i.e. degradation, volatilization, infiltration, runoff

180

loads and remaining mass in top-soil). Rainfall partitioning between infiltration and runoff

181

was determined using the Green-Ampt method37, which is adapted to the studied catchment.24

182

Retention and release of SM in the top-soil was simulated using a mixing layer approach36,38

ACS Paragon Plus Environment

Environmental Science & Technology

183

and assuming linear equilibrium sorption39 (Table S3). Despite a moderate risk of

184

volatilization for SM according to its Henry’s law constant (Table S1), significant

185

volatilization rates have been observed previously40 and were thus implemented in the model

186

(Annex S3). Long range atmospheric transport of metolachlor has never been shown to be

187

significant,41 and was therefore neglected. Between rainfall events, SM degradation and

188

volatilization were estimated by a first order dissipation rate (kb), fitted on the observed SM

189

concentrations in the top-soil. Predictions were compared with the observed dissolved

190

concentrations in runoff water (Figure S1). Degradation led to the formation of TPs and it was

191

assumed that MESA and MOXA were the only TPs for SM (in testing the full range of

192

MESA/MOXA production ratios) and that they have the degradation half-lives and sorption

193

coefficients reported in Table S1. Abiotic degradation (e.g. photolysis and hydrolysis) of SM

194

was considered negligible compared to biodegradation in soil.42 Plant uptake was not

195

explicitely integrated in the model as no correlation could be found between the remaing mass

196

of SM in the top-soil and the sugar beet growth (from 0 to 40 cm during the investigation

197

period) or the surface cover (from 0 to 100%). The confidence intervals (95% CI) of model

198

predictions, including parameter uncertainties for SM sorption and degradation, were obtained

199

from best runs (Nash-Sutcliffe criterion > 0.55, n = 138) using a Monte-Carlo method (1000

200

permutations) (Table S3).

201 202

Results and discussion

203 204

SM dissipation in the top-soil. Before applications, SM concentrations in top-soil (0 to 1 cm)

205

reached 0.006 µg gdw-1 emphasizing the low risk of SM runoff after a year. Following SM

206

applications, SM concentrations in the top-soil gradually decreased from 2.1 to 0.5 µg gdw-1

207

after 89 days (Day 120) (Figure 2B and Table S4). Changes in SM concentrations in the top-

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

Environmental Science & Technology

208

soil was fitted with the model by a dissipation rate (kb) including both degradation and

209

volatilization. Accordingly, kb was evaluated to 0.05 day-1 (DT50 = 14 days), which is in the

210

range of literature values.43,44 Among the tested range of cumulative volatilization losses, only

211

the low volatilization hypothesis, corresponding to 4.5% of SM, was consistent with SM

212

concentrations measured in top-soil (Figure S1C). However, the dissipation rate

213

underestimated the mass of SM remaining in the soil during the second part of the study

214

period (model ‘kb constant' in Figure 2B). Therefore, a gradual decrease of the rate, which

215

reached 0.005 day-1 (i.e. DT50 = 140 days) 30 days after application, was fitted to estimate the

216

SM mass in soil until the end of the study (model ‘2 successive kb’ in Figure 2B). As most of

217

the volatilization occurred during the 24h following the application40, the use of a lower

218

dissipation rate during the second study period suggests that SM sorption onto the soil

219

reduced its availability for degradation over time, thereby increasing its persistence in the

220

soil.45,46 SM reached its maximum concentration in the depth 2-5 cm (Figure 3A). Overall, the

221

SM mass remaining in the plot corresponded to 13% (± 2%) of the applied mass in the top-

222

soil (0-1 cm) after 89 days (Day 120) (or 25% (± 4%) in the first soil meter, Day 122) (Table

223

1).

224 225

The concentrations of MESA and MOXA in top-soil were systematically below the

226

quantification limits, i.e. 1 µg gdw−1. This was confirmed by MESA and MOXA predicted

227

concentrations (Figure S2). SM was enriched in the R-enantiomer before SM application

228

(∆EF = -0.057 ± 0.013) (Figure 2B). This highlights preferential degradation of the S-

229

enantiomer in the top-soil during the winter preceeding the study period. Following

230

applications, the ∆EF values in the top-soil averaged 0.001 in the soil, ranging from -0.005 to

231

0.01 (Figure 2B). These values reflect the signature of freshly applied, non-degraded SM. No

232

significant shift of enantiomer fractions could be observed between the first SM application

ACS Paragon Plus Environment

Environmental Science & Technology

233

(Day 31, April 12) and Day 99 (June 19) (Figure 2B). This suggests either no significant

234

biodegradation of SM during the first 68 days following applications, or no significant

235

enantioselective processes. Then, the R-enantiomer was preferentially degraded in the top-soil

236

(0-5 cm) (∆EF = 0.032 ± 0.004, provided as mean ± σ over the samples). In contrast, the S-

237

enantiomer was preferentially degraded in deeper soil layers (∆EF = -0.028 ± 0.013) (5-

238

100 cm) (Figure 3B). Different shifts of enantiomer fractions with increasing soil depth

239

suggest changes of bacterial community47, and/or to decreasing oxygen concentrations with

240

depth. Shifts of enantiomer fractions were also observed in lab-scale wetlands, where ΔEF of

241

SM was lower in the anoxic zone than in oxic one.21

242 243

Athough the SM pool in soil was depleted, in situ SM dissipation pathways could not

244

be precisely identified because TPs in the soil could not be detected and shifts of enantiomer

245

fractions present opposite trends in surface and deeper layers. A detailed mass balance is thus

246

needed to infer the contribution of SM dissipation pathways in the soil (i.e., degradation vs

247

export).

248 249

SM mass balance: degradation versus dissolved and particulate export. At the

250

experimental plot scale, dissolved concentrations of SM in runoff correlated with the SM

251

mass in top-soil (ps < 0.001) (Figure S1). While successive runoff gradually decreased the SM

252

pool in top-soil, SM remained largely available for off-site transport across the entire study

253

period. Based on measured concentrations and predicted discharges, SM export was estimated

254

to 3.7% of the applied SM mass, with a dissolved and particulate contribution of 97.9% and

255

2.1%, respectively (Table 1). The main predicted SM dissipation pathways in the top-soil

256

(± 95% confidence interval) were degradation (71 ± 4% of applied SM), followed by SM

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32

Environmental Science & Technology

257

leaching via infiltration (8 ± 3%), volatilization (5 ± 4%), runoff (3 ± 1%) and SM remaining

258

in the top-soil (13 ± 2%) at the end of the observation period (Day 120) (Figure S1C).

259 260

At the catchment scale, 3.4% of the applied SM was exported during the study period

261

in both the dissolved and particulate phases of runoff (Table 1). This is an order of magnitude

262

higher than export coefficients previously reported at the event scale (0.16%)4 or for several

263

years (0.072%)48 and similar to yearly exports coefficients (1.0 ± 1.6%) reported for 54 US

264

streams.14,49 Export of SM mainly occurred during nine significant runoff events (> 10 m3)

265

that represented 69% of the total outflow water volume (i.e. 17,419 m3) (Table S5), with 53%

266

for one single event (May 21). The catchment lag time, defined as the time between the

267

rainfall gravity center and the resulting peak discharge, ranged between 10 and 184 min.24

268

The contribution of drainage was 8.6 ± 6.7% during the nine major runoff events (> 10 m3)

269

and 33% during the entire study period (Table S5)24, as evaluated with the constant-slope

270

method for hydrograph separation.50 Discharge was then predominantly controlled by fast

271

flowpaths driven by Hortonian runoff (i.e unsaturated overland flow), with a lower volumetric

272

contribution of slow flowpaths driven by drainage.24 The hydrochemical signature of the drain

273

water (e.g. TSS, TIC, DIC, TOC, DOC, PO43−, Ptotal, NO2− and NH4+) also differed from that

274

of runoff from the catchment (pW < 0.05), suggesting different flowpaths (Table S6 and S7).

275

The TSS export during the study period represented 10.1 t ha-1, i.e. a top-soil loss of

276

0.4 mm over the catchment, which is in the range of regional annual erosion for a corn

277

catchment.51 Thus, 46% of the SM was exported in association with the particulate phase

278

(>0.7 µm), in agreement with previous estimates at the plot scale (20 - 46%)52,53 but higher

279

than export reported at the catchment scale (12%)54. The largest runoff event occurred on May

280

21, one week after applications, with 54 mm in 4 h and a runoff coefficient of 40.8%24, and

281

accounted for, respectively, 92%, 42%, 87% and 96% of TSS, DOC, POC and SM exported

ACS Paragon Plus Environment

Environmental Science & Technology

282

during the study period. TSS (60 g L-1), NO3- (92 mg L-1) and SM (62 µg L-1 and 34 µg g dw -1)

283

concentrations were particularly high during this intense event, whereas other physico-

284

chemical characteristics, especially DOC or POC, remained at background concentrations

285

(Figure S3). This suggests that the extreme event did not present particular characteristics,

286

except for TSS, which can be explained by the high intensity and duration of the rainfall

287

event. The study catchment is frequently prone to significant mudflows,55 and following this

288

event, the soil surface of the arable land was crusted, which enhanced runoff in subsequent

289

rainfall events. Such an erosive event may mobilize not only freshly applied SM, but also a

290

persistent fraction of SM in the soil.56 Without this event, SM export from the catchment

291

would represent only 0.1%, of which 85% of the total loads should have been found in the

292

dissolved phase (< 0.7 µm).

293 294

Rainfall intensity positively correlated with TSS and SM concentrations in both the

295

dissolved and particulate phases (pS < 0.001). The Kd values of SM negatively correlated with

296

TSS concentrations at the catchment scale (pS < 0.001) (Table S8 and Annex S5), which is

297

consistent with lower POC concentrations with increasing TSS concentrations (ps < 0.05).

298

This suggests that organic carbon of eroded materials progressively get depleted, which in

299

turn decreased SM sorption onto TSS.57,58 However, Koc values largely varied over time with

300

a large range of values (7 – 16180 ((L kg-1)/(gC g-1)), Table S8), whereas stable Koc over time

301

was expected if POC concentrations were the only reason for SM sorption. Koc values were

302

about two orders of magnitude lower than those estimated previously.57 This underlines that

303

Koc or Kd values obtained from pesticide concentrations in surface water57,59 can hardly be

304

compared as equilibration sorption time, and to a lesser extent sorbent characteristics (size of

305

the particles, aromaticity and polarity), largely vary between runoff events.

306

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32

Environmental Science & Technology

307

Overall, the complete mass balance of SM enabled to quantify the contribution of

308

degradation versus export at the plot scale. Similar export coefficients of SM at the plot and

309

catchment scales (3.7% and 3.4%) (Table 1) support the idea that SM dissipation mainly

310

occurred at the plot scale, and to a lesser extent during transport from the plot to the

311

catchment outlet. The analysis of MESA and MOXA may thus provide another line of

312

evidence to estimate SM degradation at the catchment scale.

313 314

Production and export of MESA and MOXA from soil. Concentrations of MESA and

315

MOXA were systematically below the quantification limits in soils and TSS. At the

316

experimental plot scale, MESA and MOXA in runoff generally contributed equally or more to

317

the MELSM compared to SM (%MOXA + %MESA = 52.1 ± 44.4%) (Figure 4C). The

318

seasonal MELSM export coefficient reached 8.1% of the applied SM mass, including 3.7% of

319

MOXA and 0.7 % of MESA (in SM mass equivalent). At the catchment scale, the seasonal

320

MELSM export coefficient reached 11.2% of the applied SM, including 5.5% of MOXA and

321

2.4% of MESA. MOXA was frequently detected (88% and 53% detection frequency at the

322

plot and catchment scales, respectively), and more so than MESA (25% and 29%,

323

respectively) (Table S9).

324

At the catchment scale, SM and MOXA concentrations were significantly higher

325

during runoff events compared to low-flow periods mainly controlled by drainage

326

(pW < 0.01), whereas MESA concentrations did not significantly differ (pW > 0.05). The ratios

327

of SM/MOXA/MESA at both scales highly varied over time (Fig 4C and 4D). Before major

328

applications (Day 29 for plot and Day 58 for catchment, as only 10% of SM was applied in

329

upstream areas before), TPs predominantly contributed to MELSM with mainly MOXA at the

330

plot and MESA at the catchment scale. At the plot scale, the two runoff events following

331

applications on Day 36 and 50 exported mainly freshly applied SM. Then from Day 78

ACS Paragon Plus Environment

Environmental Science & Technology

332

onwards, the contribution of MOXA to MELSM gradually increased until the end of the study

333

period (Day 120), which indicates SM degradation. SM degradation was also supported by a

334

shift of the S-enantiomer fraction in top-soil samples from Day 91 (Figure 2B). At the

335

catchment scale, the temporal variability of the ratios SM/MOXA/MESA was less clear (Fig

336

4D), likely due to overlapping and complex hydrological pathways. At Day 8 and at the end

337

of the study period (Day 120-155), most of the exported MELSM seemed to be SM instead of

338

TPs. However, this finding may be an artifact because LOQ was two orders of magnitude

339

lower for SM compared to TPs.

340 341

The high contribution of TPs to the MELSM in runoff water at both scales, although

342

TPs were not detected in the soils or TSS, suggests that MESA and MOXA were mostly

343

present in the soil solution due to their high water solubility and low Koc.12,25,42,60 The higher

344

mass export rate of MOXA compared to MESA at both scales may be due either to a smaller

345

production of MESA in the top-soil or due to more significant off-site transfer of MOXA

346

from top-soil to runoff. Both processes likely coexisted in the top-soil. Fluctuations of soil

347

moisture in top-soil (Figure 4A) with an alternation of oxic and anoxic conditions may favour

348

MOXA production in soils, as previously observed in an artificial wetland during wet/dry

349

cycles.23 In addition, MESA is less sorptive than MOXA60 (Table S1). Therefore, it may

350

preferentially leach from top-soil at the beginning of rainfall and thus be less available, when

351

ponding and runoff occur.25,35,42,61

352

An estimated TP formation ratio of 70/30 (MOXA/MESA) for unsaturated zones61

353

supports the idea that MOXA may be preferentially produced. Similarly, about two moles of

354

MOXA were exported from the catchment for each mole of MESA, which corresponds

355

approximately to the ratio of runoff (67%) to drainage (33%) volumes in the study catchment.

356

However, without the May 21 event, MELSM reached 2.8%, with 0.43% for MOXA and

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Environmental Science & Technology

357

2.27% for MESA. Hence, without the impact of the May 21 event, about four moles of MESA

358

would have been exported for each mole of MOXA giving a ratio of runoff (22%) to drainage

359

(78%). This suggests that both SM and MOXA were associated with faster surface flows,

360

such as Hortonian runoff, as a result of shorter half-life (Table S1), whereas MESA prevailed

361

in slower and longer flowpaths that sustain baseflow. This is in agreement with previous

362

predictions, which showed that MESA was associated with groundwater recharge, whereas,

363

SM stemmed from event-driven runoff

364

the unsaturated zone compared to MOXA.14,61

62

and that MESA was more deeply transported into

365 366

The TP analyses in runoff directly confirmed that at least 7.9% (5.5% for MOXA and

367

2.4% for MESA) of SM was degraded at the catchment scale. In addition, the model

368

predictions underscored that a significant fraction of MESA and MOXA leached (57% of the

369

total mass of applied SM) or was degraded (4%) at the plot scale during the study period.

370

However, the plot-based pesticide model disregards the travel time and the different pathways

371

of SM from the different plots to the catchment outlet. In addition, SM degradation may

372

generate other TPs that are mostly unknown. In this case, enantiomer analyses may strengthen

373

the evaluation of SM degradation at the catchment scale.

374 375

Enantioselective behaviour of SM in runoff. Prior to SM applications, negative ∆EF were

376

observed in runoff water from the catchment (∆EF = -0.37 ± 0.01) (Day 8) and in runoff

377

water from the plot (-0.067 ± 0.009) (Day 29) (Figure 4E and 4F). Hence, SM in discharge

378

water from the catchment outlet contained only 49.1% of the S-enantiomer (i.e. EF=0.49),

379

whereas the applied commercial product contained 86.6% of the S-enantiomer. This is

380

consistent with the finding that the SM pool available for runoff in the top-soil was also

381

depleted in the S-enantiomer (∆EF = -0.057 ± 0.013) (Day 8) (Figure 2B). The pool of SM

ACS Paragon Plus Environment

Environmental Science & Technology

382

prone to runoff may thus partly correspond to the SM fraction degraded in the top-soil where

383

enantioselective degradation occured. However, the ∆EF value of catchment discharge was

384

one order of magnitude lower than that of the plow layer at Day 8 (-0.37 and -0.057

385

respectively). Enantioselective processes may have occurred deeper in the soil as suggested

386

by higher ∆EF values in the top-soil compared to those found in deeper soil layers (Figure 3)

387

and/or during slower flowpaths driven by natural drainage in fall and winter. This also

388

supports previous findings showing that ∆EF of SM negatively correlated with catchment

389

area, suggesting that catchment retention leads to a greater loss of the S-enantiomer.63

390 391

Following SM application, no significant shift of ∆EF could be observed in runoff

392

water from the plot or the catchment (Figure 4E and 4F). Hortonian runoff prevailed during

393

the season and mainly mobilized SM, which was not enantiomerically enriched in the top-soil.

394

Enantiomeric analysis may thus not indicate degradation or export flowpaths of SM in the

395

first month following SM application. Enantioselective microbial degradation of metolachlor

396

has been shown in laboratory conditions27,64 but remains controversial in experimental

397

wetlands.21,23,65,66 Under field conditions, the fate of SM enantiomers remains poorly

398

described. Enantioselective shift was first accounted for an environmental response to the

399

introduction of non-racemic SM in the 2000s31, and enantioselective degradation was then

400

assumed in 2006-200767. In this study, SM concentrations in the different environmental

401

matrices (dissolved and particulate phases in runoff water, sediment and soil samples)

402

positively correlated with ΔEF values (pS < 0.05). This indicates that samples with low SM

403

concentrations were associated with low ΔEF, i.e. depleted in the S-enantiomer compared to

404

the commercial product. In addition, SM in a grab water sample (< 0.7 µm) collected at Day

405

253 was depleted in the S-enantionmer with ∆EF = -0.050 ± 0.009 (data not shown). This

406

confirmed that enantioselective processes, enriching SM in the R-enantiomer, occured over

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

Environmental Science & Technology

407

longer time scales and may be associated with slower flowpaths at the catchment scale.

408

Altogether, the results emphasized a slight increase of ΔEF in the top-soil during the growing

409

season, and negative ΔEF several months after the applications in fall or winter, at the

410

beginning of the new growing season. These opposite trends may be related to a seasonal

411

effect (temperature and nutrient availability) over the year, impacting the diversity of

412

microoganisms associated with SM degradation in the agricultural soil.68 A year-long field

413

survey would be required to confirm the enantioselective trends observed during SM

414

degradation.

415 416

Environmental implications for herbicide transport in agricultural catchments

417

By means of two investigative scales assisted with a modelling approach, we could

418

evaluate SM dissipation pathways at the catchment scale. A combination of analytical

419

approaches enabled to disentangle SM degradation from other non-degradative processes.

420

One significant implication of our results is that export of TSS-bound pesticides should be

421

accounted for to evaluate risks of moderately soluble pesticides transport in loess covered

422

catchments prone to extreme rainfall-runoff events. Export of SM sorbed onto TSS did not

423

prevail at the plot scale (2.1%) but represented 46% of the total load of SM exported by

424

runoff at the catchment scale, mainly due to a single and extreme rainfall-runoff event. The

425

impact of this event also underlines how severe rainfalls may alter seasonal export of SM.

426

Pesticide application during periods prone to muddy floods that occur in the beginning of

427

spring when bare soil prevails should therefore be subject to caution. Approximately 450 km²

428

of crops on sensitive loess surface soil drain into the Zorn River (Figure 1B). This underlines

429

the potential contribution of both dissolved and particulate SM flows from more than 900

430

headwater catchments with features similar to those of the study catchment. Sorption to and

431

transport with TSS may reduce, bioavailability45,69, biodegradation and toxicity70 of SM. In

ACS Paragon Plus Environment

Environmental Science & Technology

432

parallel, the gradual degradation of the S-enantiomer leads to exported SM residues, which are

433

also herbicidally less active over time. Such temporal changes of SM fate may in turn control

434

the exposure level of biota to SM and its enantiomers.

435 436

From a methodological point of view, the combined scales approach (in isolating the

437

runoff component at the plot) and the combined analysis (especially MOXA and MESA

438

ratios) may help to distinghuish flowpaths (fast and slow) contributing to catchment

439

discharge, by analogy to groundwater dating with MESA.12 This potential approach to trace

440

the contribution of different flowpaths can be reused in other agricultural catchments because

441

metolachlor is often detected concomitantly with MESA and MOXA.14

442 443

Multiple lines of evidence for SM degradation were observed. TP export represented

444

more than 50% of the total molar loads of SM exported at both plot and catchment scales, and

445

pronounced enantiomer fractionation of SM could be observed in the field. MESA and

446

MOXA concentration and loading patterns provided direct evidence of SM degradation and

447

preferential water flowpath in the short term. Complementarily, R/S enantiomers analysis may

448

emphasize enantio-selective degradation on medium and longer terms. However, only few

449

quantitative information on persistence of enantiomers and TPs have been reported so far.

450

From a regulatory point of view, there is a growing need not only to improve the

451

interpretation of physico-chemical characteristics of TPs but also to consider pesticide active

452

ingredients manufactured as isomeric mixtures. Chiral analysis of pesticides to evaluate in-

453

situ degradation processes must be done with caution because mechanistic aspects of

454

enantioselectivity of biotic and abiotic reactions have not been fully explored yet. Reference

455

experiments for enantiomer fractionation are currently missing to evaluate the extent of

456

degradation and relate it to TP formation. From a modelling point of view, predicting

ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

Environmental Science & Technology

457

pesticide transport in runoff is burdened by numerous uncertainties as sorption and

458

degradation processes may counterbalance each other.71 Models may incorporate enantiomer

459

fractionation under field conditions to improve the estimation of degradation, and thus assist

460

interpretation of complex data sets, while reducing uncertainties in risk assessments.

461 462

Acknowledgements

463

This work was funded by the European INTERREG IV program Upper Rhine

464

(PhytoRET project C.21). Support from REALISE, the Network of Laboratories in

465

Engineering and Science for the Environment in the Alsace Region (http://realise.unistra.fr) is

466

also gratefully acknowledged. Marie Lefrancq was supported by a fellowship from Région

467

Alsace. The authors thank the farmers, the soil laboratory EOST UMS 830 CNRS and

468

Martine Trautmann, Anne Véronique Auzet, Tristan Meyer, Diogo Reis, Jeanne Dollinger,

469

Thomas Dreidemy, René Boutin, Sophie Gangloff, Marie-Pierre Ottermatte, Eric Pernin and

470

Brian Sweeney for their support in data collection, sampling and analysis. We would like to

471

thank Richard Coupe and Pablo Alvarez-Zaldivar for proof-reading the manuscript and three

472

anonymous reviewers for relevant comments.

473 474

Supporting information

475

Additional information contains the summary of the physico-chemical properties of

476

SM, MESA and MOXA, detailed description of analysis standards and norms, values and

477

sources of the dissipation model parameters, SM concentrations and enantiomeric signature in

478

soil and sediment samples, hydrological characteristics of the runoff events, hydrochemical

479

data of the plot, drain and catchment outlets, detailed calculations for Kd and Koc values, data

480

on SM, MESA and MOXA detection frequencies, predicted concentrations of SM, MESA and

ACS Paragon Plus Environment

Environmental Science & Technology

481

MOXA, detailed procedure for extraction of SM, MESA and MOXA from solid samples and

482

for the SM deposition measurement during the application, detailed calculations for

483

volatilization or rainfall deposition estimation. This material is available free of charge via the

484

Internet at http://pubs.acs.org/.

485

References

486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524

(1)

Fenner, K.; Canonica, S.; Wackett, L. P.; Elsner, M. Evaluating pesticide degradation in the environment: Blind spots and emerging opportunities. Science 2013, 341 (6147), 752–758. (2) Arias-Estévez, M.; López-Periago, E.; Martínez-Carballo, E.; Simal-Gándara, J.; Mejuto, J.-C.; García-Río, L. The mobility and degradation of pesticides in soils and the pollution of groundwater resources. Agric. Ecosyst. Environ. 2008, 123 (4), 247– 260. (3) Chen, Z.; Chen, Y.; Vymazal, J.; Kule, L.; Koželuh, M. Dynamics of chloroacetanilide herbicides in various types of mesocosm wetlands. Sci. Total Environ. 2017, 577, 386– 394. (4) Doppler, T.; Lück, A.; Camenzuli, L.; Krauss, M.; Stamm, C. Critical source areas for herbicides can change location depending on rain events. Agric. Ecosyst. Environ. 2014, 192, 85–94. (5) Dechesne, A.; Badawi, N.; Aamand, J.; Smets, B. F. Fine scale spatial variability of microbial pesticide degradation in soil: Scales, controlling factors, and implications. Front. Microbiol. 2014, 5 (667). (6) Rodriguez-cruz, M. S.; Jones, J.; Bending, G. Field-scale study of the variability in pesticide biodegradation with soil depth and its relationship with soil characteristics. Soil Biol. Biochem. 2006, 38 (9), 2910–2918. (7) Macary, F.; Morin, S.; Probst, J.-L.; Saudubray, F. A multi-scale method to assess pesticide contamination risks in agricultural watersheds. Ecol. Indic. 2014, 36, 624– 639. (8) Ghafoor, A.; Moeys, J.; Stenström, J.; Tranter, G.; Jarvis, N. J. Modeling spatial variation in microbial degradation of pesticides in soil. Environ. Sci. Technol. 2011, 45 (15), 6411–6419. (9) Food and Agriculture Organization of the United Nations. FAOSTAT database collections; Rome, 2013. (10) Gadagbui, B.; Maier, A.; Dourson, M.; Parker, A.; Willis, A.; Christopher, J. P.; Hicks, L.; Ramasamy, S.; Roberts, S. M. Derived Reference Doses (RfDs) for the environmental degradates of the herbicides alachlor and acetochlor: Results of an independent expert panel deliberation. Regul. Toxicol. Pharmacol. 2010, 57 (2–3), 220–234. (11) Kalkhoff, S. J.; Vecchia, A. V.; Capel, P. D.; Meyer, M. T. Eleven-year trend in acetanilide pesticide degradates in the Iowa River, Iowa. J. Environ. Qual. 2012, 41 (5), 1566–1579. (12) Rice, C. P.; McCarty, G. W.; Bialek-Kalinski, K.; Zabetakis, K.; Torrents, A.; Hapeman, C. J. Analysis of metolachlor ethane sulfonic acid (MESA) chirality in groundwater: A tool for dating groundwater movement in agricultural settings. Sci. Total Environ. 2016, 560–561, 36–43.

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572

Environmental Science & Technology

(13) Sidoli, P.; Lassabatere, L.; Angulo-Jaramillo, R.; Baran, N. Experimental and modeling of the unsaturated transports of S-metolachlor and its metabolites in glaciofluvial vadose zone solids. J. Contam. Hydrol. 2016, 190, 1–14. (14) Rose, C. E.; Coupe, R. H.; Capel, P. D.; Webb, R. M. T. Holistic assessment of occurrence and fate of metolachlor within environmental compartments of agricultural watersheds. Sci. Total Environ. 2018, 612, 708–719. (15) Ulrich, E. M.; Morrison, C. N.; Goldsmith, M. R.; Foreman, W. T. Chiral pesticides: Identification, description, and environmental implications. In Reviews of Environmental Contamination and Toxicology Volume 217; Whitacre, D. M., Ed.; Springer US: Boston, MA, 2012; Vol. 217, pp 1–74. (16) Celis, R.; Gamiz, B.; Adelino, M. A.; Hermosin, M. C.; Cornejo, J. Environmental behavior of the enantiomers of the chiral fungicide metalaxyl in Mediterranean agricultural soils. Sci. Total Environ. 2013, 444, 288–297. (17) Gámiz, B.; Facenda, G.; Celis, R. Evidence for the effect of sorption enantioselectivity on the availability of chiral pesticide enantiomers in soil. Environ. Pollut. 2016, 213, 966–973. (18) Jing, X.; Yao, G.; Liu, D.; Qu, H.; Zhou, Q.; Zhou, Z.; Wang, P. Enantioselective toxicity and degradation of chiral herbicide fenoxaprop-ethyl in earthworm Eisenia fetida. Ecol. Indic. 2017, 75, 126–131. (19) Qu, H.; Ma, R.; Liu, D.; Gao, J.; Wang, F.; Zhou, Z.; Wang, P. Environmental behavior of the chiral insecticide fipronil: Enantioselective toxicity, distribution and transformation in aquatic ecosystem. Water Res. 2016, 105, 138–146. (20) Wang, F.; Liu, D.; Qu, H.; Chen, L.; Zhou, Z.; Wang, P. A full evaluation for the enantiomeric impacts of lactofen and its metabolites on aquatic macrophyte Lemna minor. Water Res. 2016, 101, 55–63. (21) Elsayed, O. F.; Maillard, E.; Vuilleumier, S.; Millet, M.; Imfeld, G. Degradation of chloroacetanilide herbicides and bacterial community composition in lab-scale wetlands. Sci. Total Environ. 2015, 520, 222–231. (22) Milosevic, N.; Qiu, S.; Elsner, M.; Einsiedl, F.; Maier, M. P.; Bensch, H. K. V.; Albrechtsen, H. J.; Bjerg, P. L. Combined isotope and enantiomer analysis to assess the fate of phenoxy acids in a heterogeneous geologic setting at an old landfill. Water Res. 2013, 47 (2), 637–649. (23) Maillard, E.; Lange, J.; Schreiber, S.; Dollinger, J.; Herbstritt, B.; Millet, M.; Imfeld, G. Dissipation of hydrological tracers and the herbicide S-metolachlor in batch and continuous-flow wetlands. Chemosphere 2016, 144, 2489–2496. (24) Lefrancq, M.; Van Dijk, P.; Jetten, V.; Schwob, M.; Payraudeau, S. Improving runoff prediction using agronomical information in a cropped, loess covered catchment. Hydrol. Process. 2017, 31 (6), 1408–1423. (25) Krutz, L. J.; Senseman, S. A.; Dozier, M. C.; Hoffman, D. W.; Tierney, D. P. Infiltration and adsorption of dissolved metolachlor, metolachlor oxanilic acid, and metolachlor ethanesulfonic acid by buffalograss (Buchloe dactyloides) filter strips. Weed Sci. 2004, 52 (1), 166–171. (26) McCarty, G. W.; Hapeman, C. J.; Rice, C. P.; Hively, W. D.; McConnell, L. L.; Sadeghi, A. M.; Lang, M. W.; Whitall, D. R.; Bialek, K.; Downey, P. Metolachlor metabolite (MESA) reveals agricultural nitrate-N fate and transport in Choptank River watershed. Sci. Total Environ. 2014, 473–474 (0), 473–482. (27) Ma, Y.; Liu, W. P.; Wen, Y. Z. Enantioselective degradation of rac-metolachlor and Smetolachlor in soil. Pedosphere 2006, 16 (4), 489–494.

ACS Paragon Plus Environment

Environmental Science & Technology

573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621

(28) Xu, D.; Wen, Y.; Wang, K. Effect of chiral differences of metolachlor and its (S)isomer on their toxicity to earthworms. Ecotoxicol. Environ. Saf. 2010, 73 (8), 1925– 1931. (29) Rice, C. P.; Bialek, K.; Hapeman, C. J.; McCarty, G. W. Role of riparian areas in atmospheric pesticide deposition and its potential effect on water quality. JAWRA J. Am. Water Resour. Assoc. 2016, 52 (5), 1109–1120. (30) Elsayed, O. F.; Maillard, E.; Vuilleumier, S.; Imfeld, G. Bacterial communities in batch and continuous-flow wetlands treating the herbicide S-metolachlor. Sci. Total Environ. 2014, 499, 327–335. (31) Buser, H. R.; Poiger, T.; Muller, M. D. Changed enantiomer composition of metolachlor in surface water following the introduction of the enantiomerically enriched product to the market. Environ. Sci. Technol. 2000, 34 (13), 2690–2696. (32) Shaner, D. L.; Brunk, G.; Belles, D.; Westra, P.; Nissen, S. Soil dissipation and biological activity of metolachlor and S-metolachlor in five soils. Pest Manag. Sci. 2006, 62 (7), 617–623. (33) Blaser, H.; Spindler, F. Enantioselective catalysis for agrochemicals. The case histories of (S)–metolachlor, (R)–metalaxyl and clozylacon. Top. Catal. 1997, 4 (3/4), 275–282. (34) Blaser, H. U. The chiral switch of (S)-metolachlor: A personal account of an industrial odyssey in asymmetric catalysis. Adv. Synth. Catal. 2002, 344 (1), 17–31. (35) McGrath, G. S.; Hinz, C.; Sivapalan, M.; Dressel, J.; Putz, T.; Vereecken, H. Identifying a rainfall event threshold triggering herbicide leaching by preferential flow. Water Resour. Res. 2010, 46 (2), W02513. (36) McGrath, G. S.; Hinz, C.; Sivapalan, M. Modeling the effect of rainfall intermittency on the variability of solute persistence at the soil surface. Water Resour. Res. 2008, 44 (9), W09432. (37) Green, H. W.; Ampt, G. A. Studies on soil physics. J. Agric. Sci. 1911, 4 (01), 1–24. (38) Joyce, B. A.; Wallender, W. W.; Mailapalli, D. R. Application of pesticide transport model for simulating diazinon runoff in California’s central valley. J. Hydrol. 2010, 395 (1–2), 79–90. (39) Wauchope, R. D.; Yeh, S.; Linders, J. B. H. J.; Kloskowski, R.; Tanaka, K.; Rubin, B.; Katayama, A.; Kordel, W.; Gerstl, Z.; Lane, M.; et al. Pesticide soil sorption parameters: Theory, measurement, uses, limitations and reliability. Pest Manag. Sci. 2002, 58 (5), 419–445. (40) Gish, T. J.; Prueger, J. H.; Daughtry, C. S. T.; Kustas, W. P.; McKee, L. G.; Russ, A. L.; Hatfield, J. L. Comparison of field-scale herbicide runoff and volatilization losses: An eight-year field investigation. J. Environ. Qual. 2011, 40 (5), 1432. (41) Thurman, E. M.; Cromwell, A. E. Atmospheric transport, deposition, and fate of triazine herbicides and their metabolites in pristine areas at Isle Royale National Park. Environ. Sci. Technol. 2000, 34 (15), 3079–3085. (42) Baran, N.; Gourcy, L. Sorption and mineralization of S-metolachlor and its ionic metabolites in soils and vadose zone solids: Consequences on groundwater quality in an alluvial aquifer (Ain Plain, France). J. Contam. Hydrol. 2013, 154, 20–28. (43) Lizotte, R.; Locke, M.; Bingner, R.; Steinriede, R. W.; Smith, S. Effectiveness of integrated best management practices on mitigation of atrazine and metolachlor in an agricultural lake watershed. Bull. Environ. Contam. Toxicol. 2017, 98 (4), 447–453. (44) PPDB. The Pesticide Properties DataBase (PPDB) developed by the Agriculture & Environment Research Unit (AERU) at the University of Hertfordshire, from the database that originally accompanied the EMA (Environmental Management for Agriculture) software (also developed by AERU), with additional input from the EU–

ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669

Environmental Science & Technology

(45) (46)

(47)

(48)

(49) (50) (51) (52)

(53) (54)

(55)

(56)

(57)

(58)

(59)

(60)

funded FOOTPRINT project (FP6–SSP–022704). http://www.herts.ac.uk/aeru/footprint. 2009. Louchart, X.; Voltz, M. Aging effects on the availability of herbicides to runoff transfer. Environ. Sci. Technol. 2007, 41 (4), 1137–1144. Trigo, C.; Spokas, K. A.; Hall, K. E.; Cox, L.; Koskinen, W. C. Metolachlor sorption and degradation in soil amended with fresh and aged biochars. J. Agric. Food Chem. 2016, 64 (16), 3141–3149. Dennis, P. G.; Miller, A. J.; Hirsch, P. R. Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities? FEMS Microbiol. Ecol. 2010, 72 (3), 313–327. Boithias, L.; Sauvage, S.; Taghavi, L.; Merlina, G.; Probst, J. L.; Perez, J. M. S. Occurrence of metolachlor and trifluralin losses in the Save river agricultural catchment during floods. J. Hazard. Mater. 2011, 196, 210–219. Capel, P. D.; Larson, S. J.; Winterstein, T. A. The behaviour of 39 pesticides in surface waters as a function of scale. Hydrol. Process. 2001, 15 (7), 1251–1269. McCuen, R. H. Hydrologic analysis and design; Upper Saddle river, New Jersey, 2004. Neboit, R. L’homme et l’érosion: l’érosion des sols dans le monde; Association des publications de la faculté des lettres et sciences humaines de Clermont-Ferrand, 1991. Dores, E. F. G. C.; Spadotto, C. A.; Weber, O. L. S.; Carbo, L.; Vecchiato, A. B.; Pinto, A. A. Environmental behaviour of metolachlor and diuron in a tropical soil in the central region of Brazil. Water. Air. Soil Pollut. 2009, 197 (1–4), 175–183. Buttle, J. M. Metolachlor transport in surface runoff. J. Environ. Qual. 1990, 19 (3), 531–538. El Azzi, D.; Probst, J. L.; Teisserenc, R.; Merlina, G.; Baqué, D.; Julien, F.; Payre-Suc, V.; Guiresse, M. Trace element and pesticide dynamics during a flood event in the Save agricultural watershed: Soil-river transfer pathways and controlling factors. Water. Air. Soil Pollut. 2016, 227 (12), 227–442. Heitz, C.; Flinois, G.; Glatron, S. Protection against muddy floods: perception for local actors in Alsace (France) of a protection measure (fascines). Int. Disaster Risk Conf. Davos Suisse 26-30 Août 2012 2012, 1–4. Le Gall, M.; Evrard, O.; Foucher, A.; Laceby, J. P.; Salvador-Blanes, S.; Thil, F.; Dapoigny, A.; Lefèvre, I.; Cerdan, O.; Ayrault, S. Quantifying sediment sources in a lowland agricultural catchment pond using 137 Cs activities and radiogenic 87 Sr/ 86 Sr ratios. Sci. Total Environ. 2016, 566–567, 968–980. Boithias, L.; Sauvage, S.; Merlina, G.; Jean, S.; Probst, J.-L.; Sánchez Pérez, J. M. New insight into pesticide partition coefficient Kd for modelling pesticide fluvial transport: Application to an agricultural catchment in south-western France. Chemosphere 2014, 99 (0), 134–142. Cerro, I.; Sanchez-Perez, J. M.; Ruiz-Romera, E.; Antigüedad, I. Variability of particulate (SS, POC) and dissolved (DOC, NO3) matter during storm events in the Alegria agricultural watershed. Hydrol. Process. 2014, 28 (5), 2855–2867. Taghavi, L.; Merlina, G.; Probst, J.-L. The role of storm flows in concentration of pesticides associated with particulate and dissolved fractions as a threat to aquatic ecosystems - Case study: the agricultural watershed of Save river (Southwest of France). Knowl. Manag. Aquat. Ecosyst. 2011, No. 400, 06. Capel, P. D.; McCarthy, K. A.; Barbash, J. E. National, holistic, watershed-scale approach to understand the sources, transport, and fate of agricultural chemicals. J. Environ. Qual. 2008, 37 (3), 983–993.

ACS Paragon Plus Environment

Environmental Science & Technology

670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704

(61) Bayless, E. R.; Capel, P. D.; Barbash, J. E.; Webb, R. M. T.; Hancock, T. L. C.; Lampe, D. C. Simulated fate and transport of metolachlor in the unsaturated zone, Maryland, USA. J. Environ. Qual. 2008, 37 (3), 1064–1072. (62) Huntscha, S.; Singer, H.; Canonica, S.; Schwarzenbach, R. P.; Fenner, K. Input dynamics and fate in surface water of the herbicide metolachlor and of its highly mobile transformation product metolachlor ESA. Environ. Sci. Technol. 2008, 42 (15), 5507–5513. (63) Kurt-Karakus, P. B.; Teixeira, C.; Small, J.; Muir, D.; Bidleman, T. F. Current-use pesticides in inland lake waters, precipitation, and air from Ontario, Canada. Environ. Toxicol. Chem. 2011, 30 (7), 1539–1548. (64) Mueller, M. D.; Buser, H.-R. Environmental behavior of acetamide pesticide stereoisomers. 2. Stereo- and enantioselective degradation in sewage sludge and soil. Environ. Sci. Technol. 1995, 29 (8), 2031–2037. (65) Eish, M. Y. Z. A.; Wells, M. J. M. Monitoring stereoselective degradation of metolachlor in a constructed wetland: Use of statistically valid enantiomeric and diastereomeric fractions as opposed to ratios. J. Chromatogr. Sci. 2008, 46 (3), 269– 275. (66) Klein, C.; Schneider, R. J.; Meyer, M. T.; Aga, D. S. Enantiomeric separation of metolachlor and its metabolites using LC-MS and CZE. Chemosphere 2006, 62 (10), 1591–1599. (67) Kurt-Karakus, P. B.; Bidleman, T. F.; Muir, D. C. G.; Struger, J.; Sverko, E.; Cagampan, S. J.; Small, J. M.; Jantunen, L. M. Comparison of concentrations and stereoisomer ratios of mecoprop, dichlorprop and metolachlor in Ontario streams, 2006-2007 vs. 2003-2004. Environ. Pollut. 2010, 158 (5), 1842–1849. (68) Birgander, J.; Rousk, J.; Olsson, P. A. Comparison of fertility and seasonal effects on grassland microbial communities. Soil Biol. Biochem. 2014, 76, 80–89. (69) Barbash, J. E. 11.15 - The geochemistry of pesticides. In Treatise on Geochemistry (Second Edition); Holland, H. D., Turekian, K. K., Eds.; Elsevier: Oxford, 2014; pp 535–572. (70) Karuppiah, M.; Liggans, G.; Gupta, G. Effect of river and wetland sediments on toxicity of metolachlor. Ecotoxicol. Environ. Saf. 1997, 36 (2), 180–182. (71) Gassmann, M.; Khodorkovsky, M.; Friedler, E.; Dubowski, Y.; Olsson, O. Uncertainty in the river export modelling of pesticides and transformation products. Environ. Model. Softw. 2014, 51, 35–44.

705

ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

706

Environmental Science & Technology

Figures

707 708

Figure 1. Location of the study catchment in Bas-Rhin, in France (A) within the Zorn

709

catchment (B), and the experimental setup and landuse in 2012 (C).

710 711

Figure 2. Daily rainfall and temperature, weekly top soil moisture and SM application at the

712

plot scale (A), compared with the measured and predicted remaining mass of SM in the top

713

soil (0-1 cm), with constant or successive dissipation rates (kb) model in relation with its

714

enantiomer enrichment (∆EF) (B). Error bars correspond to the analytical uncertainty

715

calculated via error propagation and incorporating accuracy and reproducibility of n = 3

716

measurements. Confidence intervals (95%) were derived from 123 runs (Nash-Sutcliffe

717

criterion > 0.55) using a Monte-Carlo approach (1000 runs).

718 719

Figure 3. SM mass in the soil profile on July 12th (Day 122) normalized by the SM mass

720

applied in relation to the organic matter content (A), and compared with the vertical

721

enantiomer enrichment in the soil profile (B). Error bars correspond to the analytical

722

uncertainty calculated via error propagation and incorporating accuracy and reproducibility of

723

n = 3 measurements.

724 725

Figure 4. Temporal changes of SM use, hydrological conditions and TSS loads at the plot (A)

726

and catchment scales (B), compared with weekly MELSM values (white diamonds: low flow,

727

black triangles: high flow, defined by a runoff event > 10 m3) with relative fractions of MESA

728

and MOXA and of dissolved and particulate SM in runoff at the plot (C), and at the catchment

729

outlets (D) in relation to enantiomer enrichments of SM in soil, runoff ( 0.7

730

µm) and sediment samples at the plot (E) and the catchment scales (F).

ACS Paragon Plus Environment

Environmental Science & Technology

731

Figure 1.

732

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

733

Environmental Science & Technology

Figure 2

734

28

ACS Paragon Plus Environment

Environmental Science & Technology

735

Page 30 of 32

Figure 3.

29

ACS Paragon Plus Environment

Page 31 of 32

736

Environmental Science & Technology

Figure 4.

737

30

ACS Paragon Plus Environment

Environmental Science & Technology

738 739

Table 1. SM mass balance from March 12 to July 10, 2012 (120 days) at the plot scale and up to August 14 (155 days) at the catchment scale. The unit haapplied refers to the SM load normalized by the surface where SM was applied. Unit [ha] [ha]

740

Page 32 of 32

Plot 7.7 10-3 7.7 10-3

Catchment 47.8 12.5

Area Application area SM application and stock in the soil Stock in the soil (top 1cm) before applications [g] 6.9 10-3 ± 1.1 10-3 nc Application in 2012 [g] 4.4 10.95 103 Application in 2012 [g haapplied-1] 571.4 874.6 Stock in the soil (top 1cm) at day 36 nc [g] (% SM applied) 2.4 ± 0.4 (55 ± 9) Stock in the soil (top 1cm) at day 120 nc [g] (% SM applied) 0.6 ± 0.1 (13 ± 2) Stock in the soil (top 1m) at day 122 nc [g] (% SM applied) 1.1 ± 0.2 (25 ± 4) Hydrology Outflow discharge (0-120 days)a [m3 day-1 ha-1] 5.2 ± 28.3 (0 - 270.2) 3.0 ± 18.1 (0 -195.2) a,b Runoff coefficient (0-120 days) [%] 30.9 ± 27.2 (0.7 - 83.1) 7.8 ± 12.6 (0.2 - 40.8) Erosion TSS export (0-120 Days) [t ha-1] 1.3 10.0 SM and TPs export in runoff/dischargeb 20.8 ± 1.7 15.8 ± 1.3 Dissolved export (0-120 days) [g haapplied-1] c -1 0.4 ± 0.1 13.5 ± 2.2 Particulate export (0-120 days) [g haapplied ] 3.7 ± 0.3 3.4 ± 0.4 Total export (dissolved, particulate) (0-120 days) [%] -1 3.9 ± 1.1 20.7 ± 5.7 MESA-dissolved export (0-120 days) [g haapplied ] 21.3 ± 8.2 47.7 ± 18.4 [g haapplied-1] MOXA-dissolved export (0-120 days) -1 46.4 ± 11 97.7 ± 27.5 Total export (dissolved, particulate, TPs) (0-120 days) [g haapplied ] 8.1 ± 1.9 11.2 ± 3.1 Total export (dissolved, particulate, TPs) (0-120 days) [%] -1 nc 98.4 ± 27.6 [g haapplied ] Total export (dissolved, particulate, TPs) (0-155 days) nc 11.2 ± 3.2 Total export (dissolved, particulate, TPs) (0-155 days) [%] a Mean ± σ (min – max), bsimulated plot runoff or measured catchment discharge, c >0.7 µm, nc: non-calculated 31

ACS Paragon Plus Environment