Degradation and sorption of the herbicides 2,4-D and quizalofop-P-ethyl and their metabolites in soils from railway tracks

Broad spectrum of activity and low potential for leaching to groundwater are important arguments for the application of the herbicide glyphosate on railway tracks. Nevertheless, certain weeds are insufficiently controlled or develop resistance, and there is also an ongoing controversial discussion about possible carcinogenicity of glyphosate. Alternatives are thus strongly desired. 2,4-D and quizalofop-P-ethyl (QE) are two selective herbicides with a complementary spectrum of activity. When used in agriculture, the compounds and their metabolites exhibit low groundwater contamination potential. Uses on railway tracks may, however, be more critical, since degradation likely is slower and mobility higher than in agricultural soils. In this study, we investigated degradation and sorption of the two active substances as well as five metabolites in three soils collected from railway tracks and in a crushed sand, used for construction works. In these railway materials, the compounds were indeed degraded slower than in agricultural soils (mean half-lives differed by a factor of 1.4–26, depending on the substance) and mobility was higher (mean sorption coefficients differed by a factor of 3–19). Half-lives and sorption coefficients were also estimated by extrapolation of data from agricultural soils, considering the organic carbon content of railway soils and agricultural soils. Estimated values were more conservative than measured values. Based on our experimental data, possible leaching to groundwater is predicted to be highest for 2,4-D and quizalofop-acid, the primary metabolite of QE, moderate for 2,4-dichlorophenol, the primary metabolite of 2,4-D, but low for QE. Secondary and tertiary metabolites were formed in only low quantities. For herbicides, for which no measured parameters are available for railway soils, estimated values may also be a viable alternative for a first tier groundwater assessment.


Background
In many countries, weed control on railway tracks relies on glyphosate [1][2][3][4]. Glyphosate is a broad-spectrum herbicide with systemic properties [5] and, therefore, shows good activity not only against annual, but also most perennial weeds, including weeds with deep roots that are commonly present on railway tracks. A further advantage is that glyphosate and also its metabolite AMPA exhibit low potential for leaching to groundwater [6][7][8].
Despite these advantages, weed control on railway tracks should not depend on glyphosate alone. Some weeds are not sufficiently controlled, for example, horsetails (Equisetum) with their acicular stems and extensive rhizomes [2,9]. Long-term use of glyphosate may also promote development of resistance [4,5]. Moreover, there is an ongoing, controversial discussion about possible carcinogenicity of glyphosate [10][11][12]. Alternatives are thus strongly desired, along with herbicides that complement the activity of glyphosate on railway tracks.
However, active substances with a wide spectrum of activity comparable to glyphosate are rare. In Europe, for example, the broad-spectrum herbicides paraquat or imazapyr are not registered and approval of diquat, glufosinate, or amitrole has not been extended. Therefore, selective herbicides need to be considered for use on railway tracks, although in this case, for sufficient weed control, two or more active substances with a complementary activity are required.
Auxine-type herbicides and ACCase-inhibitors are two prominent groups of selective herbicides. For this study, we selected one compound of each group. The auxin-type herbicide 2,4-D acts against most annual and perennial broad-leaved weeds with some activity against sedges and rushes [13,14]. In particular, it shows activity against horsetails [15,16]. The second compound, quizalofop-P-ethyl, is an ACCase-inhibitor and has a rather narrow spectrum of activity against grasses of the Poaceae family [17]. Altogether, the two compounds complement each other reasonably well. The compounds are systemic in plants and they are applied after emergence of weeds, i.e., they allow for treatment of single plants and do not need to be applied on the whole surface of railway tracks.
In addition to the spectrum of activity, groundwater protection is of particular importance when evaluating herbicides for use on railway tracks. The tracks are designed to efficiently drain rainwater. This ensures their stability and avoids deformation due to freezing water. Soils below the railway ballast usually have low organic carbon contents, which results in weak sorption of most organic compounds. A few decades ago, atrazine and diuron were widely used for weed control on railway tracks, and in several cases, contamination of aquifers with these compounds or their metabolites was linked to uses on railways [1,18]. Glyphosate and AMPA, however, are strongly retained by mineral substrates such as iron oxides and are thus less susceptible to leaching [6,19,20]. To avoid the risk of groundwater contamination, potential glyphosate alternatives must exhibit sufficient sorption and/or fast degradation in soils below the railway ballast.
Based on data for agricultural topsoils, it can be concluded that 2,4-D and quizalofop-P-ethyl pose low risk of groundwater contamination [13,17,21]. However, topsoils are normally removed when constructing railway tracks and for the remaining subsoils, these basic environmental parameters are often not available.
The aim of the present study was, therefore, to determine sorption and degradation rate coefficients for 2,4-D and quizalofop-P-ethyl in a number of railway soils, selected to represent a realistic worst-case regarding both properties. In soils, 2,4-D is transformed to 2,4-dichlorophenol and 2,4-dichloroanisole ( Fig. 1) [13]. Quizalofop-P-ethyl is rapidly cleaved to quizalofop-P-acid that is further transformed to two hydroxylated metabolites, 3-OH-quizalofop-acid and 3-OH-CQO ( Fig. 1) [17]. Sorption and degradation studies were thus also performed with these 5 metabolites. Finally, we investigated whether such substance properties for railway soils may also be estimated from data for agricultural soils.
The internal standard 2,4-DCA-D 3 was prepared by methylation of 2,4-DCP-D 3 using diazomethane [22]. Quizalofop-acid-D 3 was prepared from quizalofop-ethyl-D 3 by alkaline hydrolysis. For that, 2.5 mg of quizalofopethyl-D 3 was dissolved in 0.5 mL methanol and diluted with 10 mL of a sodium hydroxide solution (0.01 M). After 4 h, the solution was acidified with sulfuric acid to pH ≈ 2 (note that the hydrolysis product quizalofop-acid-D 3 partly precipitated in the acidic, aqueous solvent). Thereafter, the solution/suspension was partitioned with 10 mL of dichloromethane and another three times with 5 mL of dichloromethane. The combined organic phase was then passed through a small column filled with anhydrous sodium sulphate to remove residual water. After complete evaporation of dichloromethane, the residues were dissolved in 10 mL of methanol. Quizalofop-acid-D 3 was used as internal standard for quizalofop-acid and 3-OH-quizalofop-acid.

Soil sampling
Soil samples were collected in September 2016 at three locations from 10-to 15-year-old railway tracks in Switzerland. At location Münchenbuchsee (MBS), beigecolored subsoil with some rusty concretions was taken at a depth of 16-55 cm below the railway ballast. At Müntschemier (MR), white-beige subsoil was taken at a depth of 2-15 cm. From site Erlenbach (EB), we collected a 7-cm thick layer of gray, sandy construction material below the ballast. Sampling depths were selected based on earlier construction geologic investigations with the aim to collect soils with high sand contents. Since the early 90s, weed control at these sites has been undertaken as spot treatments with glyphosate, and before that with atrazine. In addition to these three field samples, new crushed sand (CS), typically used in the construction of drainage channels in renovated tracks, was obtained from the Swiss Federal Railways. Some properties of the soils are listed in Table 1 (for convenience, we use the term "soil" for all four materials). The field-moist soils were sieved (2 mm) and stored in plastic bags at 4 °C.

Soil incubation
Incubation experiments were set up 7-9 months after collection of the soils. Separate experiments were conducted with all test substances (2,4-D, quizalofop-Pethyl, and 5 metabolites) in all 4 soils, yielding a total of 28 incubations. Portions of 800 g field-moist soil were spread in a crystallizing dish (width, 23 cm) and the test substance was applied evenly to the soil surface with a polyethylene spray bottle. For that, 800 µg 2,4-D, quizalofop-P-ethyl, or quizalofop-P-acid, or 80 µg 2,4-DCP, 2,4-DCA, 3-OH-quizalofop-acid, or 3-OH-CQO, were dissolved in 5-10 mL water (note that the spike solution of quizalofop-P-ethyl and quizalofop-P-acid contained 40% and 16% methanol, respectively). The resulting spike levels were 1 µg/g for 2,4-D, quizalofop-P-ethyl, and quizalofop-P-acid, and 0.1 µg/g for the other metabolites. After spiking, the soils were thoroughly mixed and filled into 2-L Erlenmeyer flasks that were closed with airpermeable cellulose plugs. The soils were then incubated at 20 °C in the dark under aerobic conditions for up to 121 days. The soil moisture (adjusted 2 weeks before start  Tables 2 and 3) was maintained by regular addition of distilled water. At appropriate time intervals, aliquots of 10 g soil were removed, filled into 20-mL glass vials, and stored in a freezer at − 20 °C until extraction. To samples from incubation studies with 2,4-D, 2,4-DCP, and 2,4-DCA, 10 mL methanol were added immediately after sampling to minimize possible volatilization of the metabolites during storage.

Extraction of 2,4-D and metabolites
All samples of an incubation experiment with a particular compound were processed on the same day. The methanolic soil suspensions were fortified with internal standards (1 µg 2,4-D-D 3 , 2,4-DCP-D 3 , and 2,4-DCA-D 3 in 100 µL methanol). After vigorous shaking (≈ 1 min), the suspensions were centrifuged (≈ 1500g for 5 min; Eppendorf 5804, Hamburg, Germany) and the supernatants were transferred to 40-mL glass vials. A second extraction was performed with 10 mL deionized water, followed by 7 min centrifugation. The combined extracts were then acidified with H 2 SO 4 (1:4 in water) to pH ≈ 2 and partitioned with 10 mL dichloromethane and then a second time with 5 mL dichloromethane. The combined dichloromethane phases were evaporated to ≈ 200 µL on a heating plate under a gentle stream of nitrogen. After addition of 50 µL 0.5% trifluoroacetic acid in methanol, the samples were derivatized using approximately 1 mL diazoethane in methyl tert-butyl ether (for 30 min, preparation described in [23]). Then, the reaction mixture was allowed to evaporate to ≈ 200 µL. After addition of ≈ 100 mg of anhydrous sodium sulphate, the samples were quantitatively transferred with hexane to 1.5-mL clear glass autosampler vials. Once more the solvent volume was reduced to ≈ 200 µL and filled up to 1 mL with hexane.

Extraction of quizalofop-ethyl and metabolites
As for 2,4-D and metabolites, all samples of an incubation experiment were processed on the same day. The soil samples (10 g) were suspended in 10 mL methanol, internal standards (1 µg each of quizalofop-ethyl-D 3 and quizalofop-acid-D 3 in 100 µL methanol and 1 µg 2,3-dihydroxy-quinoxaline in 100 µL methanol/ water (50% v/v)) were added, and the vials were vigorously shaken (≈ 1 min) and then centrifuged (≈ 1500g for 10 min). The supernatants were transferred to 20-mL glass vials. A second extraction was performed with 10 mL aqueous 0.01 M CaCl 2 , followed by 10 min centrifugation. Aliquots of 0.4 mL of the combined extracts were then transferred to 1.5-mL clear glass autosampler vials and diluted with 0.6 mL water.

GC-MS analysis of 2,4-D and metabolites
2,4-D (as ethyl ester), 2,4-DCP (as ethyl ether), 2,4-DCA, and the corresponding internal standards were analysed on a GC-MS/MS system consisting of a PAL autosampler (CTC Analytics, Zwingen, Switzerland), a HP6890N gas chromatograph (Agilent Technologies, Santa Clara, CA), and a Quattro Micro mass spectrometer (Waters Corporation, Milford, MA). GC conditions were as follows: 1 µL split/splitless injection (250 °C, initial 60 s splitless); 30-m column (internal diameter, 320 µm), coated with (5%-phenyl)-methylpolysiloxane with a film thickness of 0.25 µm (BGB5, BGB Analytics, Boeckten, Switzerland); temperature program: 50 °C, 2 min isothermal, 25 °C/min to 90 °C, 5 °C/min to 160 °C, 25 °C/min to 280 °C; constant flow, 4 mL/min helium. The MS was operated in electron impact ionization (70 eV, 180 °C) and selected-ion-monitoring mode. (166), respectively. As the mass differences between target compounds and deuterated internal standards are small, the masses of the 37 Cl 35 Cl-rather than the 35 Cl 2isotopologues were selected for the internal standards to minimize interferences. A typical chromatogram of a soil extract is shown in Additional file 1: Fig. S1 (left). Quantification was based on peak area ratios relative to the internal standards and in reference to suitable standard solutions. Information on recoveries, precision, storage stability, and limits of quantification are given in the Additional file 1.

LC-MS/MS analysis of quizalofop-ethyl and metabolites
Quizalofop-ethyl, the three metabolites, and the corresponding internal standards were analysed with liquid chromatography-tandem mass spectrometry. The instrument was configured with an autosampler (HTS PAL, CTC), a binary HPLC pump (Agilent 1100, with microvacuum degasser), and a triple quadrupole mass spectrometer (API 4000, with turbo ion spray source, Sciex,

Table 2 Panel A: Half-lives (DT 50 ) and DT 90 values of 2,4-D, 2,4-DCP, and 2,4-DCA. Panel B: Half-lives (DT 50 ) and formation fractions of metabolites resulting from sequential fitting, not considering the reversibility of the reaction 2,4-DCP to 2,4-DCA
T: temperature, M: moisture in % of water holding capacity a SFO: single first-order kinetics, DFOP: double first-order in parallel model, FOMC: first-order multi-compartment model. χ 2 values < 10% b Calculated according to [26], for details see Additional file 1: Framingham, MA). LC conditions were as follows: injection via a 100 µL PEEK loop; Gemini NX column (C18, 150 × 2 mm, 5 µm particle size, protected by a 4 × 2 mm pre-column with the same stationary phase, Phenomenex, Torrance, CA); gradient elution with the solvents formic acid (1‰ v/v in water) and methanol (initial conditions, 20% methanol, linear increase to 85% during 6.5 min, linear increase to 100% during 4 min, 3 min isocratic hold, initial conditions re-established within 1 min, followed by an equilibration time of 5 min); flow, 0.2 mL/min. The MS was operated in negative (ion spray voltage, − 4.2 kV, 400 °C) or positive mode (5 kV, 400 °C; only for detection of quizalofop-ethyl; time window, 11-12 min) and the ion transitions listed in Additional file 1: Table S1 were monitored. As for 2,4-D and metabolites, the mass differences between target compounds and deuterated internal standards are small. Therefore, the ion transitions of the 37 Cl-rather than the 35 Cl-isotopologues were selected for the internal standards to minimize interferences.
A representative chromatogram of a soil extract is shown in Additional file 1: Fig. S1 (right). Quantification was based on peak area ratios relative to the internal standards and in reference to matrix-matched standard solutions. For that, untreated soil was extracted as described above and target compounds and internal standards were added to the extracts directly before transfer to the autosampler vial and dilution with water. For quality assurance data, see Additional file 1.

Soil adsorption experiments
Batch adsorption experiments were performed in accordance with OECD guideline 106 at 20 °C in the dark [24]. Separate experiments were performed with all test substances (2,4-D, quizalofop-P-ethyl, and 5 metabolites) in all 4 soils, at 5 concentration levels each and in duplicate, yielding a total of 280 batch adsorption tests. In all tests, the indirect method was used, where only the remaining Table 3  b Calculated according to [26], for details see Additional file 1: concentration in the aqueous phase was measured after achievement of adsorption equilibrium. The water content of the sieved, field-moist soils was accounted for in the calculation of soil:solution ratios. Typically, aliquots of moist soil corresponding to 10 g dry weight were weighed into 40 mL clear glass vials with Teflon-lined screw caps. To these soils, aqueous 0.01 M CaCl 2 solution was added to yield a water volume of exactly 9 mL (including the water content of the soil samples). The vials were capped, mounted on a reciprocal shaker (90 oscillations/min), and agitated for 24 h for preequilibration. Then, 1 mL of spiking solution containing the test substance was added (resulting soil:solution ratio, 1:1) and the vials were agitated for an additional 24 h. The soil slurries were then centrifuged at ≈ 1500g for 10 min and an aliquot of the clear supernatant was transferred to an autosampler vial, where internal standard was added. The solutions were analysed by LC-MS/MS. Preliminary experiments with 2,4-DCA indicated that volatilisation may have affected the solute concentrations in some samples. Therefore, smaller vials (20 mL) were used for the experiments to minimise the air volume above the soil suspension. The same precautionary measures were used for 2,4-DCP.
Further deviations from this procedure included (i) different soil:solution ratios (Additional file 1: Table S3), (ii) different equilibration times (for substances with low stability, Additional file 1: Table S3), and (iii) analysis by GC-MS, after liquid-liquid extraction with dichloromethane and solvent exchange to hexane (only 2,4-DCA). More details on preliminary tests, adsorption isotherms, and LC-MS/MS analysis of 2,4-D and 2,4-DCP are given in the Additional file 1.

Kinetic analysis
Kinetic parameters for the degradation of the individual test substances as well as metabolites formed during the experiments were determined using the software Cake [25]. For the incubated substances, we tested the single first-order model (SFO), the double first-order in parallel model (DFOP, this bi-phasic model assumes two compartments in which the compound is degraded according to first-order kinetics, but with different rate constants [26]), the first-order multi-compartment model (FOMC, this bi-phasic model assumes a continuum of micro-compartments in which the compound is degraded according to first-order kinetics [27]), and the hockey-stick model (HS, this bi-phasic model assumes two sequential firstorder curves with a breakpoint at a certain time [26]), whereas for metabolites formed during the experiments, we always used SFO kinetics. To determine kinetic formation fractions of metabolites, their concentrations were expressed in "parent" equivalents (note that a formation fraction is defined here as the molar fraction of a test substance being transformed to a respective metabolite). Initial concentrations were adjustable, those of the metabolites were set to 0). For fitting, the iteratively reweighted least squares optimizer was selected. Fits were only accepted, when statistically significant parameters (based on the 95% confidence interval) could be determined, with a χ 2 error < 15% [26]. Visual assessment and residuals were further acceptance criteria. Two types of DT 50 values (half-lives) were derived: (i) values resulting from the best-fitting model for optimal description of the degradation curves and (ii) more conservative values describing the slow phase of degradation in case of bi-phasic kinetics. The latter are generally used in leaching models (so-called modeling endpoints) and served here also for a comparison of data for railway soils with agricultural soils.

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][29][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 CO 2 -C/ kg/h ( Table 1) and was thus in the range of values reported for Swedish railway soils (0.01-0.32 mg CO 2 -C/ kg/h) [30], but clearly lower than in agricultural topsoils [31].
Overall, the soils selected represent a realistic worstcase 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 DT 50 and DT 90 values in the respective soils (Table 2, Panel A; DT x is defined as the time required for x% dissipation). The DT 50 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 DT 50 of the entire experimental duration was considered adequate.
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 DT 50 values of 1-66 days ( DT 50 values of 2,4-DCA, the secondary metabolite of 2,4-D, also varied considerably between soils, with very fast elimination in MR (DT 50 , 0.3 days), but slow elimination in CS (101 days). 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, Panel A). DT 90 values were between 2 and 335 days.

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 DT 50 values for groundwater assessments [26]. Furthermore, DT 50 values shall be normalized to a reference temperature of 20 °C and a reference moisture of pF2, and the geometric mean of the normalized DT 50 values shall be used in pesticide leaching models (modeling endpoints) [34]. Following these recommendations, we determined conservative geomean DT 50 values of 115, 12, and 14 d for 2,4-D, 2,4-DCP, and 2,4-DCA (Table 2, Panel A).
These values were then compared with DT 50 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, Panel A), 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 DT 50 values, Table 2, Panel A).
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 Additional file 1) and formation fractions including volatilization may have been higher.
DT 50 values resulting from sequential fitting generally were in good agreement with DT 50 values determined in experiments with direct incubation of the respective compounds ( Table 2, Panel A and Panel B). 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 parametersfor 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 DT 50 and DT 90 values are listed in Table 3 (Panel A).
Degradation of QE was extremely fast in soils MBS, MR, and CS with DT 50 values < 1 day, but clearly slower in soil EB with a DT 50 value of 21 days (Fig. 4, left). The bi-phasic degradation was most pronounced in soils EB and CS with high DT 90 values of > 1000 days and 68 days, respectively. Degradation of the three metabolites was consistently slower. The primary metabolite QA was degraded with a DT 50 value of ≈ 100 days in soils EB and MBS. In CS, the DT 50 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 days, Fig. 4, third panel). Finally, the tertiary metabolite 3-OH-CQO was degraded with DT 50 values of 24-474 days. 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 (Panel A) were again calculated from normalized DT 50 values, considering only the slow phase in case of bi-phasic kinetics (modeling endpoints). These DT 50 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.
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, Panel B). For the transformation of 3-OH-QA to 3-OH-CQO, no 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 c w is the concentration in the aqueous phase and c s the concentration in soil. The Freundlich adsorption coefficients (K F ) and the Freundlich exponents (1/n) were determined from linear regressions of log c w vs log c s . 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 K F 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 (K F , 1.5-66 mL/g) was about two orders of magnitude stronger than adsorption of 2,4-D (K F , 0.04-0.28 mL/g), which may be assumed to be present as a carboxylate anion in the railway soils (pK a value, 3.4) [13]. K F 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 wider confidence intervals (Additional file 1: Tables S6 and S7). For a precise determination of K F values, a higher fraction of adsorbed test compounds would have been desirable.    Geometric mean organic carbon normalized adsorption coefficients (K Foc ) are currently used as input values for groundwater assessments in the context of registration in Europe [34]. In the pesticide leaching models, K Foc values are converted to K F values based on the organic carbon content of the soils (C org ) implemented in the leaching scenarios. The geomean K F 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 K F , adsorption in the railway soils was clearly weaker than in agricultural soils, by a factor of 3-19.
Adsorption was non-linear for all compounds, with mean Freundlich exponents of 0.62-0.80 (Tables 4 and  5), indicating weaker adsorption at higher concentrations. The mean 1/n values were consistently lower in railway soils than in agricultural soils, i.e., non-linearity of adsorption was more pronounced.

Estimation of adsorption and degradation in railway soils from data with agricultural soils
Since adsorption and degradation data are usually not available for railway soils, the question may be posed whether such data can be estimated by extrapolation from agricultural soils. Sorption to organic matter is often considered as the predominant sorption mechanism for organic compounds in soils [40], i.e., it is assumed that K F values are roughly proportional to the organic carbon content of soils (C org ). K F values for railway soils (r) may thus simply be estimated from K F values for agricultural soils (a), correcting for C org : With this formula, we calculated K F (r) values for the seven test substances in this study (2,4-D, QE, and their metabolites), based on the geometric means of K F (a), C org (a), and C org (r) (Additional file 1: Tables S6 and S7). These estimates were generally lower (1.5-7.6×) and thus more conservative than our measured adsorption coefficients, except for 2,4-DCP, where estimated and measured K F values were almost equal. In the railway soils, sorption of the test substances (except 2,4-DCP) to mineral substrates may be important as well and the above formula does not account for their contribution to sorption.
Freundlich exponents (1/n) for agricultural soils were consistently higher than those determined in this study for railway soils (Tables 4 and 5), i.e., the use of 1/n data from agricultural soils would also be more conservative regarding the prediction of leaching in soils. Overall, the K F (r) = K F (a) C org (r) C org (a) .
estimation of sorption in railway soils from data with agricultural soils using the above approach seems to be conservative, in some cases very conservative. Sorption parameters measured specifically in railway soils may thus not be essential for a first tier groundwater assessment. However, it needs to be highlighted that this conclusion is based on only seven test substances. Degradation half-lives in railway soils may be estimated in a similar way by extrapolation of data from agricultural soils, assuming that biological activity is linked to the organic matter content of a soil. This approach was, for example, implemented in the HardSPEC model, a simple tool developed in the UK for estimating surface water and groundwater exposure resulting from herbicides applied to hard surfaces [41]: DT 50 (r) values estimated with this formula (excluding 2,4-D) yielded 1.7-13 × higher values than actually measured (Additional file 1: Tables S4 and S5), i.e., as for adsorption, this approach tends to be conservative as well. For 2,4-D, the estimated mean half-life was ≈ 2 × lower, when considering only the initial slow phase of the degradation curves.
Concerning formation of metabolites, our experiments with railway soils showed substantial differences to agricultural soils. Generally, metabolites were formed in lower amounts. In particular, formation of the secondary and tertiary metabolite of QE (3-OH-QA and 3-OH-CQO) seems to be unimportant in railway soils.

Conclusions
Groundwater protection is of particular importance when applying herbicides on railway tracks. The degradation and adsorption studies presented here allow a first assessment of the potential leaching of 2,4-D, quizalofop-P-ethyl, and their metabolites to groundwater. Based on these data, leaching is predicted to be highest for 2,4-D and quizalofop-acid. Degradation of these two compounds was considerably slower and the adsorption weaker in railway soils than in agricultural soils. For 2,4-DCP, it can be assumed that leaching is moderate. Its adsorption to railway soils was clearly weaker than in agricultural soils, but elimination was comparatively fast. Even less leaching is expected for 2,4-DCA. This metabolite was also formed in relatively low quantities of ≤ 2.6%. Finally, negligible leaching is expected for quizalofop-P-ethyl. Its degradation was fast in railway soils and adsorption was relatively strong. For 3-OH-QA and 3-OH-CQO, a high inherent potential for leaching would DT 50 (r) = DT 50 (a) C org (a) C org (r) .