Worst-case character of the railway soils with regard to sorption and degradation
The soils used in this study exhibited high sand (61-84%) and low clay (5-15%) contents and, consequently, low water holding capacities (0.2-0.3 g/g, Table 1), properties which lead to rapid vertical transport of water. Similar soil textures were also assumed in three scenarios that were developed for authorization of herbicides on railway tracks in Germany [28]. These scenarios were implemented in the computer model PELMO that simulates possible leaching of pesticides in soil.
The pH-values of our railway soils were high (7.5-7.7) so that weak acids (e.g., quizalofop-P-acid, 2,4-D, or 2,4-DCP) were present in their more mobile, anionic form. Sorption is therefore expected to be weaker than in acidic soils. Furthermore, the organic carbon contents of <0.06-0.31% (Table 1) were among the lowest values reported for subsoils below railway tracks (0.06-2.3% for the fine material) [18, 28-30]. Low organic carbon and low clay contents are expected to result in weak sorption of most organic compounds.
Finally, the basal respiration, a measure for the microbial activity of soils, varied from 0.07 to 0.11 mg CO2-C kg-1h-1 (Table 1) and was thus in the range of values reported for Swedish railway soils (0.01-0.32 mg CO2-C kg-1h-1) [30], but clearly lower than in agricultural topsoils [31].
Overall, the soils selected represent a realistic worst-case regarding sorption and degradation of organic compounds for the use on railway tracks. However, the soil properties intentionally deviate from those recommended to study degradation of pesticides in agricultural soils, e.g., according to OECD guideline 307 on aerobic and anaerobic transformation in soil, particularly in terms of organic matter content and range of textures and pH values [32].
Degradation of 2,4-D and its metabolites in railway soils
Degradation of 2,4-D differed considerably in the four soils (Fig. 2, left). In soil MR and the crushed sand (CS), a pronounced lag-phase was observed with no or very slow initial degradation, followed by a sharp decline in concentration after about 100 and 50 d, respectively. An initial phase with slower degradation was also found in soil MBS (and, to a certain degree, in EB).
Bi-phasic degradation of pesticides in soil incubation studies normally goes along with a deceleration over time [26]. Occasionally, also an increase of the degradation rate is observed [33]. Microorganisms may require a certain adaptation time to cope with contaminants such as 2,4-D. When transferring the soils from the railway tracks into the laboratory environment, the conditions for microorganisms are changed. In particular, the soil structure is disturbed by sieving, with a possible negative impact on microorganisms. After some time, however, microorganisms may recover and eventually be capable of degrading 2,4-D. Also note that the soils had been stored for 7-9 months prior to start of the incubation experiments which may have affected certain microorganisms.
On the other hand, the observed sharp decline may also be the result of another laboratory artifact. Microorganisms, which were originally not present in the soil, may have been introduced during collection, sieving, sampling, or moisturizing, which may have altered the degradation behavior. It is thus not clear whether such a pronounced increase in the degradation rate of 2,4-D would also occur under natural conditions in soils below railway tracks.
In the kinetic analysis, we considered the slow initial phase as representative for natural conditions, which resulted in conservative estimates of the DT50 and DT90 values in the respective soils (Table 2, top; DTx is defined as the time required for x% dissipation). The DT50 values in this initial phase ranged from 39 d in MBS to 203 d in CS, while in MR, no degradation was observed during the initial 84 d (>1000 d). For soil EB, even though visually bi-phasic, an SFO DT50 of the entire experimental duration was considered adequate.
Table 2
Half-lives (DT50) and DT90 values of 2,4-D, 2,4-DCP, and 2,4-DCA.
|
pH
|
Corg
|
T [°C]/
|
DT50/DT90 [d], best fitb
|
|
(CaCl2)
|
[%]
|
M [%]a
|
2,4-D
|
2,4-DCP
|
2,4-DCA
|
soils from railway tracks
|
EB
|
7.7
|
0.15
|
20/35
|
35/115 SFO
|
27/129 DFOP
|
82/271 SFO
|
MBS
|
7.6
|
0.31
|
20/60
|
39/131 SFOd
|
1.0/7.0 FOMC
|
9.2/31 SFO
|
MR
|
7.6
|
≈0.04
|
20/72
|
>1000/>1000 SFOd
|
4.1/18 DFOP
|
0.3/2.1 DFOP
|
CS
|
7.5
|
<0.06
|
20/65
|
203/675 SFOd
|
66/519 FOMC
|
101/335 SFO
|
|
|
|
DT50 [d], modeling endpointsc
|
geom. mean, 20°C, pF 2
|
≈0.1
|
|
115
|
12
|
14
|
agricultural soils [13]
|
min-max
|
6.2-7.8e
|
0.9-3.7
|
20/50
|
1.6-95
|
3.2-6.2
|
11-16
|
geom. mean, 20°C, pF 2
|
1.5/1.8/1.8f
|
|
4.4
|
7.0
|
10
|
Half-lives (DT50) and formation fractions of metabolites resulting from sequential fitting, not considering the reversibility of the reaction 2,4-DCP to 2,4-DCA.
|
sequence of fitting
|
DT50 [d]
2,4-Dg
|
formation fraction [%]h
|
DT50 [d]
2,4-DCPg
|
formation fraction [%]h
|
DT50 [d]
2,4-DCAg
|
EB
|
2,4-D→2,4-DCP→2,4-DCA→sink
|
35 SFO
|
34
|
35 SFO
|
32
|
33 SFO
|
|
2,4-DCP→2,4-DCA→sink
|
|
|
28 DFOP
|
12
|
103 SFO
|
MBS
|
2,4-D→2,4-DCP→sink
|
40 SFOd
|
16
|
2.6 SFOd
|
|
|
|
2,4-DCP→2,4-DCA→sink
|
|
|
1.0 FOMC
|
6.4
|
16 SFO
|
MR
|
2,4-DCP→2,4-DCA→sink
2,4-DCA→2,4-DCP→sink
|
|
|
4.1 DFOP
4.8 SFO
|
2.2
79i
|
11 SFO
0.3 DFOP
|
CS
|
2,4-DCP→2,4-DCA→sink
|
|
|
66 FOMC
|
2.5
|
101 SFO
|
aT: temperature, M: moisture in % of water holding capacity
bSFO: single first-order kinetics, DFOP: double first-order in parallel model, FOMC: first-order multi-compartment model. χ2 values <10%.
cCalculated according to [26], for details see supplementary information, Table S4 and kinetic fits
dOnly the slow initial phase was considered.
eMeasured in H2O
fIn degradation studies with 2,4-D/2,4-DCP/2,4-DCA
gBest fit. χ2 values <12%, except for 2,4-DCA in CS (20%).
hVolatilization may have contributed to overall dissipation (i.e. actual formation may have been higher).
iBack reaction from 2,4-DCA to 2,4-DCP
Dissipation of 2,4-DCP, the primary metabolite of 2,4-D, when applied as test substance, was consistently faster, particularly in soils MBS and MR (Fig. 2, second panel). For 2,4-DCP, a bi-phasic decline was observed as well, but with the more typical deceleration over time. The dissipation curves were well described with the bi-phasic DFOP or FOMC model, resulting in DT50 values of 1-66 d (Table 2, top). DT90 values ranged from 7 to 519 d, indicating the clearly slower degradation at later time points.
DT50 values of 2,4-DCA, the secondary metabolite of 2,4-D, also varied considerably between soils, with very fast elimination in MR (DT50, 0.3 d), but slow elimination in CS (101 d). 2,4-DCA concentrations could be adequately described by the SFO model, except for soil MR, where a better fit was obtained with the DFOP model (Fig. 2, right; Table 2, top). DT90 values were between 2 and 335 d.
Comparison with agricultural soils
The European “FOrum for the Coordination of pesticide fate models and their USe” (FOCUS) recommends that, in case of bi-phasic degradation, only the slow phase should be considered when deriving DT50 values for groundwater assessments [26]. Furthermore, DT50 values shall be normalized to a reference temperature of 20°C and a reference moisture of pF2, and the geometric mean of the normalized DT50 values shall be used in pesticide leaching models (modeling endpoints) [34]. Following these recommendations, we determined conservative geomean DT50 values of 115, 12, and 14 d for 2,4-D, 2,4-DCP, and 2,4-DCA (Table 2, top).
These values were then compared with DT50 values typically found in agricultural soils. For 2,4-D, numerous studies are available in public literature [21]. However, for a robust comparison, we used kinetic endpoints from studies that were performed for registration of 2,4-D in Europe [13] as these studies were evaluated according to the above cited recommendations of FOCUS. Degradation of 2,4-D was considerably slower in our railway soils than in these agricultural soils (on average by a factor of 26, Table 2, top), which was expected based on the fact that the microbial activity (approximated by the measured basal respiration) was at least one order of magnitude lower in the railway soils (Table 1) than in agricultural soils. Nonetheless, 2,4-D showed a remarkably slow degradation in our soils. To what extent this may be the consequence of the experimental conditions in the laboratory, cannot be answered here (see the above discussion on the lag-phase). Elimination of the metabolites 2,4-DCP and 2,4-DCA was, however, only 1.7 and 1.4× slower in the railway soils than in agricultural soils, respectively (again based on geomean DT50 values, Table 2, top).
Formation of 2,4-D metabolites
Highest amounts of 2,4-D metabolites were observed in soil EB, where 2,4-DCP reached 14% after 49 d and 2,4-DCA 2.6% after 120 d (Fig. 3, left; formation expressed in “parent” equivalents). In the other soils, metabolites were found in much lower amounts of ≈1% (2,4-DCP) and <0.5% (2,4-DCA). In incubation experiments with 2,4-DCP, the secondary metabolite 2,4-DCA was observed in amounts of up to 7% in EB (Fig. 3, second panel), 5% in MBS, and 1% in MR and CS. However, the methylation of 2,4-DCP to 2,4-DCA was reversible, particularly in soils MR and CS, with considerable formation of 2,4-DCP in incubation experiments with 2,4-DCA (up to 55 and 22%, respectively, Fig. 3, right).
Kinetic parameters for formation and further degradation of metabolites were fitted assuming sequential reactions (2,4-D → 2,4-DCP → 2,4-DCA, 2,4-DCP → 2,4-DCA, or 2,4-DCA → 2,4-DCP). Statistically significant DT50 values and formation fractions could be determined in a number of experiments and corresponding results are given in Table 2 (bottom). Some fits of experiments with notable metabolite formation are depicted in Fig. 3.
For the transformation from 2,4-D to 2,4-DCP, kinetic formation fractions of 34% and 16% were obtained for soils EB and MBS, respectively. The subsequent transformation from 2,4-DCP to 2,4-DCA occurred with highly differing formation fractions of 2-32%, again with the highest formation in soil EB. For the back reaction from 2,4-DCA to 2,4-DCP, only one statistically significant formation fraction could be derived, which was 79% in soil MR. Note that 2,4-DCP and 2,4-DCA formed from the respective precursors may have volatilized to some extent (see supplementary information) and formation fractions including volatilization may have been higher.
DT50 values resulting from sequential fitting generally were in good agreement with DT50 values determined in experiments with direct incubation of the respective compounds (Table 2, top and bottom). Differences were primarily observed when a bi-phasic model was used in direct incubation experiments and SFO kinetics in experiments with a precursor compound. Such differences are commonly observed, even if the same kinetic model is used for fitting. Half-lives may indeed be different, whether a compound is directly spiked to a soil or formed in soil from a precursor [35, 36].
In the kinetic analysis of the experiments with 2,4-D and its metabolites, the reversibility of the reaction from 2,4-DCP to 2,4-DCA should, in principle, be considered. However, in many experiments, formation of at least one of the two metabolites was too low for such a kinetic analysis. In other experiments, the kinetic analysis did not result in statistically significant fitting parameters – for example, in incubation experiments with 2,4-DCP and 2,4-DCA in soil MR. Note that the model Aquasim [37] was used to consider the interconversion of the two metabolites, with a similar approach as in [35, 38].
Degradation of quizalofop-P-ethyl and its metabolites in railway soils
Degradation curves of quizalofop-P-ethyl (QE), quizalofop-P-acid (QA), 3-OH-QA, and 3-OH-CQO in the four railway soils are shown in Fig. 4. Most experiments were better fitted by a bi-phasic model, except for 3-OH-QA where SFO fits were acceptable. The resulting DT50 and DT90 values are listed in Table 3 (top).
Table 3
Half-lives (DT50) and DT90 values of QE, QA, 3-OH-QA, and 3-OH-CQO.
|
pH
|
Corg
|
T [°C]/
|
DT50/DT90 [d], best fitb
|
|
(CaCl2)
|
[%]
|
M [%]a
|
QE
|
QA
|
3-OH-QA
|
3-OH-CQO
|
soils from railway tracks
|
|
EB
|
7.7
|
0.15
|
20/49
|
21/>1000 FOMC
|
96/414 DFOP
|
37/122 SFO
|
335/>1000 FOMC
|
MBS
|
7.6
|
0.31
|
20/54
|
0.22/1.4 DFOP
|
113/488 DFOP
|
35/116 SFO
|
24/>1000 FOMC
|
MR
|
7.6
|
≈0.04
|
20/67
|
0.21/0.98 FOMC
|
>1000/>1000 SFO
|
630/>1000 SFO
|
445/>1000 SFO
|
CS
|
7.5
|
<0.06
|
20/55
|
0.69/68 FOMC
|
371/>1000 SFO
|
189/628 SFO
|
474/>1000 FOMC
|
|
|
|
DT50 [d], modeling endpointsc
|
geom. mean, 20°C, pF 2
|
≈0.1
|
|
4.0
|
276
|
101
|
246
|
agricultural soils [17]
|
|
min-max
|
5.0-8.2d
|
0.8-4.6
|
20-22/40-70
|
0.3-1.1
|
7-182
|
7-69
|
42-258
|
geom. mean, 20°C, pF 2
|
2.8/1.9/1.5/2.1e
|
|
0.4
|
24
|
18
|
63
|
Half-lives (DT50) and formation fractions of metabolites resulting from sequential fitting.
|
sequence of fitting
|
DT50 [d]
QEf
|
formation fraction [%]
|
DT50 [d] QAf
|
formation fraction [%]
|
DT50 [d] 3-OH-QAf
|
EB
|
QE→QA→3-OH-QA→sink
|
28 FOMC
|
100
|
92 SFO
|
0.5
|
180 SFO
|
MBS
|
QE→QA→3-OH-QA→sink
|
0.20 DFOP
|
100
|
105 SFO
|
3.0
|
91 SFO
|
MR
|
QE→QA→sink
|
0.21 DFOP
|
98
|
968 SFO
|
|
|
CS
|
QE→QA→sink
|
1.1 FOMC
|
100
|
746 SFO
|
|
|
aT: temperature, M: moisture in % of water holding capacity
bSFO: single first-order kinetics, DFOP: double first-order in parallel model, FOMC: first-order multi-compartment model. χ2 values <10%.
cCalculated according to [26], for details see supplementary information, Table S5 and kinetic fits
dMeasured in H2O
eIn degradation studies with QE/QA/3-OH-QA/3-OH-CQO
fBest fit. χ2 values ≤12%, except for 3-OH-QA in MBS (16%).
Degradation of QE was extremely fast in soils MBS, MR, and CS with DT50 values <1 d, but clearly slower in soil EB with a DT50 value of 21 d (Fig. 4, left). The bi-phasic degradation was most pronounced in soils EB and CS with high DT90 values of >1000 d and 68 d, respectively. Degradation of the three metabolites was consistently slower. The primary metabolite QA was degraded with a DT50 value of ≈100 d in soils EB and MBS. In CS, the DT50 value was ≈1 year, and in MR, no degradation was observed (Fig. 4, second panel). Degradation curves of the secondary metabolite 3-OH-QA were qualitatively similar to QA, but half-lives were 2-3× lower (35-630 d, Fig. 4, third panel). Finally, the tertiary metabolite 3-OH-CQO was degraded with DT50 values of 24-474 d. In comparison to the other metabolites, a slower degradation was observed in soil EB (Fig. 4, right).
As for 2,4-D and its metabolites, no clear correlation could be identified between soil parameters and rate of degradation. Degradation of all QE metabolites was slowest in soil MR (Fig. 4). Degradation of the parent compound itself was fastest in this soil, likely because of abiotic hydrolysis of QE to QA. Note that the water content was highest in soil MR. The geometric means listed in Table 3 (top) were again calculated from normalized DT50 values, considering only the slow phase in case of bi-phasic kinetics (modeling endpoints). These DT50 values were 4-12× higher than those reported for agricultural soils [17].
Formation of quizalofop-P-ethyl metabolites
Quizalofop-P-ethyl is a so-called pro-herbicide, which is, due to its lipophilicity, more readily taken up through the cuticle of leaves. As with other “FOP” herbicides [23], the compound is rapidly transformed in plants to the corresponding acid, which is the herbicidally active substance.
In the railway soils, quizalofop-P-ethyl was also quantitatively transformed to quizalofop-acid (Fig. 5, left). However, the secondary and tertiary metabolites, 3-OH-QA and 3-OH-CQO, were formed in much lower quantities (≤1%, shown for soil EB in Fig. 5, second panel), also in experiments with QA (third panel). Only in experiments with incubation of 3-OH-QA, 3-OH-CQO reached amounts of up to 29% after 4 months (Fig. 5, right).
Sequential fitting confirmed the quantitative transformation of QE to QA (formation fractions, 98-100%). For the subsequent transformation of QA to 3-OH-QA, formation fractions of 0.5% and 3.0% were determined in soil EB and MBS, respectively (Table 3, bottom). For the transformation of 3-OH-QA to 3-OH-CQO, no statistically significant formation fraction could be determined, but in soil EB, it must have been >29% (maximum formation at the end of incubation).
In degradation studies with agricultural soils, formation fractions of 70-100%, 32-76%, and 100%, respectively, were determined for the reaction sequence QE→QA→3-OH-QA→3-OH-CQO [17]. Hydroxylation of QA to 3-OH-QA thus seems to be clearly less important in railway soils than in agricultural soils.
Adsorption to railway soils
Adsorption experiments were evaluated with the Freundlich model [39],
where cw is the concentration in the aqueous phase and cs the concentration in soil. The Freundlich adsorption coefficients (KF) and the Freundlich exponents (1/n) were determined from linear regressions of log cw vs log cs. These so-called Freundlich isotherms for 2,4-D, QE, and their metabolites in the four railway soils are shown in Figs. 6 and 7, and corresponding KF and 1/n values are listed in Tables 4 and 5.
The compounds exhibited strongest adsorption to the construction material from site EB and weakest adsorption to subsoil from MR, except for 3-OH-CQO, which showed strongest adsorption to subsoil MBS and weakest adsorption to the crushed sand (CS) (Figs. 6 and 7). Adsorption of QE (KF, 1.5-66 mL/g) was about two orders of magnitude stronger than adsorption of 2,4-D (KF, 0.04-0.28 mL/g), which may be assumed to be present as a carboxylate anion in the railway soils (pKa value, 3.4) [13]. KF values of the metabolites were in between those of the two parent compounds. The weakest adsorption of 2,4-D was found in soil MR, where the adsorbed fraction represented only 3-18%, depending on the concentration level. In such cases, determination of the concentration in soil is less accurate, resulting in higher confidence intervals (Tables S6 and S7). For a precise determination of KF values, a higher fraction of adsorbed test compounds would have been desirable.
Geometric mean organic carbon normalized adsorption coefficients (KFoc) are currently used as input values for groundwater assessments in the context of registration in Europe [34]. In the pesticide leaching models, KFoc values are converted to KF values based on the organic carbon content of the soils (Corg) implemented in the leaching scenarios. The geomean KF values obtained in our study on railway soils can be compared with typical values for agricultural soils (Tables 4 and 5). For this comparison, we only relied on studies accepted in the context of registration [13, 17]. Based on geomean KF, adsorption in the railway soils was clearly weaker than in agricultural soils, by a factor of 3-19.
Table 4
Freundlich adsorption coefficients (KF) and exponents (1/n) of 2,4-D, 2,4-DCP, and 2,4-DCA.
|
pH
(CaCl2)
|
Corg [%]
|
KF [mL/g]
2,4-D
|
1/n
|
KF [mL/g]
2,4-DCP
|
1/n
|
KF [mL/g]
2,4-DCA
|
1/n
|
soils from railway tracks
|
EB
|
7.7
|
0.15
|
0.28
|
0.69
|
3.32
|
0.83
|
10
|
0.58
|
MBS
|
7.6
|
0.31
|
0.14
|
0.75
|
0.96
|
0.75
|
1.4
|
0.77
|
MR
|
7.6
|
≈0.04
|
0.04
|
0.62
|
0.08
|
0.63
|
0.43
|
0.47
|
CS
|
7.5
|
<0.06
|
0.07
|
0.61
|
0.14
|
0.88
|
0.93
|
0.65
|
geom. mean
arithm. mean
|
|
≈0.1
|
0.10
|
0.67
|
0.43
|
0.77
|
1.5
|
0.62
|
agricultural soils [13]a
|
min-max
|
5.0-7.5
|
0.9-4.4
|
0.19-0.83
|
0.78-0.90
|
3-25
|
0.80-0.94
|
10-27
|
0.85-0.95
|
geom. mean
arithm. mean
|
|
1.8
|
0.40
|
0.83
|
8.3
|
0.88
|
18
|
0.92
|
aSoils M800-M822 considered
Confidence intervals are given in Table S6
Table 5
Freundlich adsorption coefficients (KF) and exponents (1/n) of QE, QA, 3-OH-QA, and 3-OH-CQO.
|
pH
(CaCl2)
|
Corg [%]
|
KF [mL/g]
QE
|
1/n
|
KF [mL/g]
QA
|
1/n
|
KF [mL/g]
3-OH-QA
|
1/n
|
KF [mL/g]
3-OH-CQO
|
1/n
|
soils from railway tracks
|
|
|
|
EB
|
7.7
|
0.15
|
66
|
0.76
|
5.4
|
0.81
|
2.1
|
0.78
|
2.0
|
0.74
|
MBS
|
7.6
|
0.31
|
47
|
0.97
|
0.99
|
0.79
|
0.62
|
0.77
|
3.1
|
0.74
|
MR
|
7.6
|
≈0.04
|
1.5
|
0.77
|
0.17
|
0.76
|
0.16
|
0.84
|
0.52
|
0.65
|
CS
|
7.5
|
<0.06
|
6.5
|
0.70
|
0.96
|
0.64
|
0.83
|
0.60
|
0.35
|
0.48
|
geom. mean
arithm. mean
|
|
≈0.1
|
13
|
0.80
|
0.96
|
0.75
|
0.65
|
0.75
|
1.0
|
0.65
|
agricultural soils [17]
|
|
|
|
min-max
|
4.3-8.1
|
0.06-5.9
|
15-99
|
0.83-0.88
|
0.19-40
|
0.69-0.89
|
0.8-10
|
0.80-1.07
|
5.5-22
|
0.59-0.80
|
geom. mean
arithm. mean
|
|
2.3/1.0/2.0/2.0a
|
39
|
0.86
|
5.2
|
0.81
|
2.6
|
0.93
|
9.9
|
0.66
|
aIn adsorption studies with QE/QA/3-OH-QA/3-OH-CQO
Confidence intervals are given in Table S7