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