Mineralization of the Common Groundwater ... - ACS Publications

Aug 17, 2016 - Sand filters (SFs), commonly used in DWTPs to remove iron and manganese(7, 8) were recently shown to harbor the capacity to biodegrade ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Northern Illinois University

Article

Mineralization of the common groundwater pollutant 2,6dichlorobenzamide (BAM) and its metabolite 2,6-dichlorobenzoic acid (2,6-DCBA) in sand filter units of drinking water treatment plants Johanna Vandermaesen, Benjamin Horemans, Julie Degryse, Jos Boonen, Eddy Walravens, and Dirk Springael Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b01352 • Publication Date (Web): 17 Aug 2016 Downloaded from http://pubs.acs.org on August 18, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

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 19

Environmental Science & Technology

1

Mineralization

of

the

common

groundwater

pollutant

2,6-

2

dichlorobenzamide (BAM) and its metabolite 2,6-dichlorobenzoic acid

3

(2,6-DCBA) in sand filter units of drinking water treatment plants

4

Johanna Vandermaesen†, Benjamin Horemans†, Julie Degryse§, Jos Boonen§, Eddy Walravens§ and Dirk

5

Springael†*

6



7

Belgium

8

§

9

Heverlee, Belgium

KU Leuven, Division of Soil and Water Management, Kasteelpark Arenberg 20 bus 2459, B-3001 Heverlee,

Centraal laboratorium, De Watergroep, Researchpark Haasrode Leuven 1834 - Technologielaan 23, B-3001

10

* Corresponding author: [email protected], tel. +3216321604, fax. +3216321997

11

Abstract

12

The intrinsic capacity to mineralize the groundwater pollutant 2,6-dichlorobenzamide (BAM) and its metabolite

13

2,6-dichlorobenzoic acid (2,6-DCBA) was evaluated in samples from sand filters (SFs) of drinking water

14

treatment plants (DWTPs). Whereas BAM mineralization occurred rarely and only in SFs exposed to BAM,

15

2,6-DCBA mineralization was common in SFs, including those treating uncontaminated water. Nevertheless,

16

SFs treating BAM contaminated water showed the highest 2,6-DCBA mineralization rates. For comparison,

17

2,6-DCBA and BAM mineralization were determined in various topsoil samples. As in SF samples, BAM

18

mineralization was rare, whereas 2,6-DCBA mineralization capacity appeared widespread, with high

19

mineralization rates found especially in forest soils. Multivariate analysis showed that in both SF and soil

20

samples, high 2,6-DCBA mineralization correlated with high organic carbon content. Adding a 2,6-DCBA

21

degradation deficient mutant of the BAM mineralizing Aminobacter sp. MSH1 confirmed that 2,6-DCBA

22

produced from BAM is rapidly mineralized by the endogenous microbial community in SFs showing intrinsic

23

2,6-DCBA mineralization. This study demonstrates that (i) 2,6-DCBA mineralization is widely established in

24

SFs of DWTPs, allowing the mineralization of any 2,6-DCBA produced during BAM degradation and (ii) the

25

first metabolic step in BAM mineralization is rare in microbial communities, rather than its further degradation

26

beyond 2,6-DCBA.

27

ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 19

28

Introduction

29

2,6-dichlorobenzamide (BAM) is a transformation product of dichlobenil (2,6-dichlorobenzonitrile), a herbicide

30

mainly used for weed control in public and private areas.1 BAM has a high water solubility and low Koc and

31

easily leaches to groundwater,1,2 leading to the widespread occurrence of BAM as a groundwater pollutant with

32

concentrations up to 5 µg/L.1,3,4 The presence of BAM in groundwater resources used for drinking water

33

treatment results in the costly closure of abstraction wells5 or implementation of expensive measures in drinking

34

water treatment plants (DWTPs), such as activated carbon filtration and advanced oxidation processes6 to reach

35

the EU drinking water threshold concentration of 0.1 µg/L (80/778/EEC). Sand filters (SFs), commonly used in

36

DWTPs to remove iron and manganese7,8 were recently shown to harbor the capacity to biodegrade organic

37

micropollutants such as pesticides and pharmaceuticals.9–12 Whether SF microbial communities also degrade

38

BAM is however unknown, but not expected based on its recalcitrance in both soil and groundwater.13–17 To

39

biologically remove BAM in DWTPs, operational conditions should thus be optimized or alternative

40

approaches developed such as bioaugmentation, i.e., inoculating BAM mineralizing bacteria in SFs.18

41

Although uncertainties exist about the metabolic pathway of BAM mineralization, it is generally accepted that

42

BAM is first converted to 2,6-dichlorobenzoic acid (2,6-DCBA). Holtze et al.14 identified 2,6-DCBA as an

43

intermediate of dichlobenil biodegradation based on its transient accumulation in soil. This is supported by

44

findings of Simonsen et al.19 for the BAM mineralizing isolates Aminobacter sp. MSH1 and ASI1. In those

45

strains, a constitutively expressed amidase converts BAM at high rate to 2,6-DCBA20 but further degradation

46

proceeds more slowly and 2,6-DCBA was appointed a bottleneck for BAM mineralization.19 2,6-DCBA was

47

suggested as a potential groundwater contaminant based on its persistence and low sorption to soil.21 However,

48

data on the actual presence of 2,6-DCBA in soil or groundwater are scarce since 2,6-DCBA is usually not

49

included in monitoring campaigns. Nevertheless, 2,6-DCBA was detected in 1.1% of 4739 Danish groundwater

50

samples taken between 2012 and 201422 and in soils after treatment with dichlobenil.23,24 Furthermore, no

51

consensus exists on biodegradability of 2,6-DCBA. Holtze et al.15 reported 2,6-DCBA as degradable in some

52

soils, but persistent in others, with half-lives ranging between 24 and 108 days. Most studies on 2,6-DCBA

53

biodegradability focused on disappearance and not mineralization.

54

Even when BAM is degraded by either the endogenous community or by introduced organisms in SFs of

55

DWTPs, it is thus unclear whether potentially produced 2,6-DCBA is removed by endogenous microbial

56

activity. This knowledge is of importance for acceptance and implementation of biological BAM removal

57

approaches in DWTPs. Therefore, in this study, mineralization of BAM and 2,6-DCBA was evaluated in SF

58

samples from DWTPs differing in intake water BAM contamination. For comparison, the BAM and 2,6-DCBA

59

mineralization capacity was determined in topsoil samples from 17 locations, differing in land use and soil type.

60

Finally, BAM mineralization was evaluated in filter sand showing high 2,6-DCBA mineralization, after adding ACS Paragon Plus Environment

Page 3 of 19

Environmental Science & Technology

61

a 2,6-DCBA degradation deficient variant of Aminobacter sp. MSH1, to examine whether 2,6-DCBA produced

62

during biological conversion of BAM was efficiently mineralized by endogenous SF populations.

63

Materials and methods

64

Sand filter and soil samples. SFs were sampled at 11 Belgian DWTPs between 2011 and 2015. Since they

65

were permanently inundated and regularly mixed by backwashing, and thus expected to be physicochemically

66

and biologically homogenous, one sample was taken at a random position. Samples were taken from the top 20

67

cm, transferred to sterile plastic containers and stored for maximally one month at 4°C until analysis. Pesticide

68

mineralization assays performed with the same sample after different storage times did not differ in

69

mineralization kinetics for at least 6 months after sampling. Table 1 shows the SF samples, their origin and

70

BAM contamination in the DWTP intake water. BAM concentrations in the SF influent were unavailable but

71

not expected to differ from intake water concentrations since the treatment steps implemented before sand

72

filtration do not remove BAM. 2,6-DCBA concentrations were unavailable since its measurement is not

73

included in routine water monitoring campaigns. DWTPs Kluizen, De Blankaart, Antwerp, Eeklo, Zele and

74

Snellegem treated water containing BAM and other pesticides. The intake water treated at DWTP Klein-Sinaai

75

contained other pesticides, but no BAM, and the remaining four DWTPs (Saint-Léger, Bree, Haacht and

76

Lommel) treated pristine groundwater in which no pesticide residues were detected. All SFs were rapid SFs

77

except the one from Antwerp, which was a slow SF. Soil samples were taken from the top 20 cm from different

78

locations in Belgium. At several locations (Heverlee, Kessel-Lo, Linkhout, Rotselaar and Schoonderbuken),

79

samples were taken from different areas in close proximity to each other, with different land uses (Table 2). An

80

overview of the implemented treatment steps at all DWTPs, physicochemical and biological characteristics of

81

SF and soil samples and a map with all sampled locations are given in the supporting information.

82

Chemicals, bacterial strains and culture conditions. [Ring-U-14C]-labeled BAM with > 95% purity was

83

purchased from Izotop (Institute of Isotopes Co., Ltd., Budapest, Hungary) and dissolved in acetone.

84

Aminobacter sp. strain M6.100g is a genetically characterized mutant of the BAM mineralizing Aminobacter sp.

85

MSH1.25 Strain M6.100g stoichiometrically converts BAM to 2,6-DCBA but is deficient in further degradation

86

of 2,6-DCBA and hence in mineralization of the aromatic moiety. M6.100g was grown from frozen stocks on

87

R2A.26 After 4 days of incubation at 25°C, a smear of colonies was taken, inoculated in 50 mL R2B26

88

containing 10 mg/L BAM (Sigma-Aldrich, Steinheim, Germany) and incubated for 2 more days at 25°C on an

89

orbital shaker.

90

Production of 14C-2,6-DCBA. 14C-labeled 2,6-DCBA was produced by incubating [ring-U-14C]-labeled BAM

91

with strain M6.100g. M6.100g cells were harvested by centrifugation (4000g, 15 min, 15°C) and inoculated in ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 19

92

50 mL MMO medium27 amended with 1500000 counts per minute (CPM) 14C-BAM (100 µg/L) at a cell density

93

of 108 cells/mL in flasks equipped with a glass vial containing 1 mL 0.5 M NaOH to trap 14CO2. After 2 days of

94

incubation at 20°C, the deficiency of 14CO2 production was confirmed by measuring radioactivity in the NaOH

95

solution as described.28 After centrifugation (9000g, 15 min),

96

supernatant by solid phase extraction (SPE) using Oasis MAX SPE cartridges (Waters, Zellik, Belgium)

97

preconditioned with 1 mL methanol, 1 mL mQ-H2O and 0.5 mL PO4-buffer (pH 7.5). The supernatant was

98

loaded at a flow rate of 0.1 mL/min. Thereafter, the cartridge was washed with 1 mL mQ-H2O and dried under

99

vacuum. Neutral compounds, including residual

14

C-2,6-DCBA was extracted from the

14

C-BAM, were eluted with 1 mL methanol.

14

C-2,6-DCBA

100

was eluted using a solution of 90% methanol, 7% H2O and 3% HNO3. The purity of the 14C-2,6-DCBA extract

101

was evaluated using reverse phase ultrahigh performance liquid chromatography as described,20 but eluting with

102

only an isocratic flow of 15% acetonitrile and 85% mQ water (acidified with H3PO4 to pH 2.5) for 8 min. The

103

chromatogram was compared with a

104

DCBA were observed, i.e., other compounds such as BAM, ortho-chlorobenzamide, ortho-chlorobenzoic acid,

105

benzamide or benzoate were not detected. Radioactivity was measured28 in both the 2,6-DCBA fraction

106

(fraction eluting from 4.5 to 6.5 min run time, retention time of 2,6-DCBA is 5.3 min) and the rest of the eluted

107

sample, showing that the

108

substrate. The identity of the

109

mutant of MSH1 that lacks the ability to convert BAM to 2,6-DCBA but that mineralizes 2,6-DCBA.

110

Microcosm mineralization assays. Mineralization assays were performed in triplicate in 10 mL vials

111

containing 2.5 g of wet SF material/soil as described.28 Taking into account a detection limit for radioactivity of

112

30 CPM per measurement, a total radioactivity of 15000 CPM 14C-BAM or 14C-2,6-DCBA, was used to enable

113

accurate determination of kinetic mineralization parameters such as the maximum mineralization rate, which

114

was generally < 1 %.d-1 in case of BAM. After evaporation of acetone or methanol, pure 14C-BAM or 14C-2,6-

115

DCBA was dissolved in 25 µL mineral MMO medium (spike concentration of 2 mg/L) and added to 2.5 g SF or

116

soil material (wet weight). Four hundred µL sterilized mQ-H2O was added to the soil microcosms to attain

117

water contents similar to these of the SF samples, leading to final BAM or 2,6-DCBA concentrations of 100 to

118

200 µg/L in the SF/soil water phase. Each assay was replicated three times and complemented with an abiotic

119

control (soil or SF sample sterilized by autoclaving three times at 121 °C for 20 min with 24 h periods in

120

between). Microcosms were incubated on an orbital shaker at 20°C for 50 to 100 days during which 14CO2 was

121

trapped using NaOH and

122

plotting the cumulative percentage 14CO2 as a function of incubation time, subtracted by the background activity

123

measured for abiotic controls (of which none showed mineralization, data not shown).

124

Evaluation of mineralization kinetics. The mineralization capacity in every mineralization assay was

125

evaluated based on the extent of mineralization (after 50 days of incubation, pmax) and the maximum linear

12

C-2,6-DCBA standard. No other peaks than that associated with 2,6-

14

C-2,6-DCBA extract had a purity of 94%, i.e., identical to this of the

14

14

14

C-BAM

C-2,6-DCBA was additionally confirmed by showing its mineralization by a

CO2-radioactivity measured.28 Cumulative mineralization curves were obtained by

ACS Paragon Plus Environment

Page 5 of 19

Environmental Science & Technology

126

mineralization rate (kmax) derived as the maximal slope of the mineralization curve using three consecutive data

127

points. To study mineralization kinetics in more detail, mineralization curves were modelled with the Three-

128

half-order model,29 which combines first and zero order mineralization kinetics with linear growth of the

129

degrading organisms: ࡼ = ࢖૚ ቆ૚ − ࢋ

࢚࢘² ି࢑૚࢚ି൬ ൰ ૛ ቇ+

࢑૙࢚

130 131

in which P is the percentage mineralization at time t, p1 is the total extent of mineralization by the first order

132

process, k0 and k1 are the zero and first order rate constants, respectively, and r is the linear growth rate

133

constant. Parameter values were estimated by nonlinear regression analysis, using the lsqnonlin command in

134

Matlab (Mathworks) with the default trust-region-reflective algorithm, at a termination tolerance of 10-14 and

135

allowing maximally 200000 function evaluations and 30000 iterations. Initial parameter estimates were set at

136

100, 5, 5 and 5 for p1, k0, k1 and r, respectively.

137

The Three-half-order model is most applicable to mineralization curves that show a lag phase, which indicates

138

adaptation/growth of the BAM/2,6-DCBA degrading organisms before significant mineralization is detected.

139

The presence of such a lag phase is reflected by r-values > 0. When r approached zero, i.e., when the average

140

value of r for three replicates was less than 0.001, the modelling procedure was repeated after deleting the linear

141

growth term from the Three-half-order model, reducing it to the First-zero-order model.30 This did not affect the

142

estimated parameter values, but improved their accuracy, i.e., their standard errors decreased. For soil samples

143

He1, KL1, Li2 and SL, r-values predicted with the Three-half-order model were greater than 0.001, but the

144

First-zero-order model yielded better prediction of the mineralization curves, based on visual comparison

145

(Figure S2). r-values for these mineralization assays were therefore set to 0 and parameter values were predicted

146

using the First-zero-order model. Mineralization was considered significant when pmax > 5%.

147

Partial least squares regression (PLSR). Multivariate analysis using PLSR was performed to correlate

148

physicochemical and biological characteristics of DWTP intake water/SF samples or soil samples with kinetic

149

parameters of 2,6-DCBA mineralization in two separate models. DWTP Antwerp was excluded since no raw

150

water data were available. PLSR was performed using SOLO software (Eigenvector Research Inc., Manson,

151

WA, USA) as described in the supporting information. As independent variables for the model corresponding to

152

the SF samples (designated as the SF model), the characteristics shown in Table S1 were used, including

153

exposure to BAM and other micropollutants, nutrient concentrations, organic carbon concentrations, oxygen

154

availability, pH, temperature, iron and manganese concentrations and the bacterial 16S rRNA gene copy

155

number. As independent variables for the model corresponding to the soil samples (the soil model), the

156

characteristics shown in Table S2 were used, including soil texture data, pH, nutrient concentrations, organic

157

carbon concentrations, (specific) UV-absorbance as a measure of aromatic compound content31 and the bacterial ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 19

158

16S rRNA gene copy number. Soil moisture content was not included since mineralization assays were

159

performed with adjusted water content (see above). The dependent variables used for both models were the

160

kinetic parameters of 2,6-DCBA mineralization, i.e., pmax, kmax, p1, k0, k1 and r (Table S3). Latent variables were

161

selected based on simultaneously capturing most of the variance in the independent variables and most of the

162

co-variance with all six dependent variables in the same model.

163

Bioaugmentation assay. M6.100g cells were harvested and diluted to 107 cells/mL in 0.01M MgSO4, of which

164

50 µL was added to 14C-BAM mineralization assays performed with SF sample K4 as described above. As non-

165

inoculated controls, 50 µL sterile 0.01M MgSO4 was added. Mineralization was monitored during 33 days.

166

Results

167

BAM mineralization capacity in sand filters. BAM mineralization was evaluated in SF samples from 11

168

DWTPs. Results for the extent of mineralization after 50 days of incubation (pmax) and the maximum linear

169

mineralization rate (kmax) are shown in Table 1. Using a threshold of pmax > 5%, substantial BAM mineralization

170

was recorded for samples from Kluizen, Snellegem and Eeklo. Sample K1 showed a sigmoidal BAM

171

mineralization curve (Figure 1A) with a pmax and kmax of 34% and 1%.d-1, respectively. In samples Sn1 and E1,

172

BAM mineralization proceeded more slowly (kmax = 0.3 – 0.5%.d-1) and reached a lower extent (pmax = 6 – 8%).

173

To assess whether the BAM mineralization capacity observed for DWTP Kluizen, Snellegem and Eeklo was

174

consistent in time, extra SF samples were taken at different time points (samples K2, K3, K4, E2, E3 and Sn2).

175

Mineralization capacity generally remained constant and an increasing trend was even observed, between

176

sampling time point 1 and 2 to 4 and sampling time point 2 and 3 in samples from Kluizen and Eeklo,

177

respectively (Table 1, Figure S3).

178

Comparison of the different kinetic parameter values (Figure S4) provides indications about the degradation

179

processes that lead to mineralization. Kinetic parameters could not accurately be estimated for sample E3 using

180

the Three-half-order model. However, for samples E1, E2, Sn1 and Sn2, r equaled 0, which indicates that the

181

observed slow BAM mineralization did not support growth of the degrading organisms. On the other hand, the

182

sigmoidal pattern of the BAM mineralization curves in samples from Kluizen was reflected by r-values > 0 and

183

the trend of increasing mineralization capacity (with shorter lag phases) for samples from Kluizen taken at

184

different time points was reflected by differences in r, but not by k0 and k1.

185

2,6-DCBA mineralization capacity in sand filters. 2,6-DCBA mineralization was tested in a selection of SF

186

samples, i.e., K4, E3, Sn2, StL, Br, Hch and Lo, including SFs that either treated BAM contaminated or

187

uncontaminated water and SFs that showed BAM or no BAM mineralization. All samples tested positive for

188

2,6-DCBA mineralization (Figure 1B), but mineralization kinetics differed greatly (Table 2, Figure 2A). SFs ACS Paragon Plus Environment

Page 7 of 19

Environmental Science & Technology

189

treating BAM contaminated water (K4, E3, Sn2) all showed high 2,6-DCBA mineralization capacity with pmax-

190

and kmax-values reaching up to 63% and 16%.d-1, respectively. Results for samples treating pristine groundwater

191

were more divergent. Intermediate mineralization capacity was found in samples Hch and Lo, with pmax- and

192

kmax-values of 27 to 35% and 5 to 7%.d-1, respectively; while 2,6-DCBA mineralization in samples Br and StL

193

was very low, with kmax < 0.3%.d-1 and pmax < 6%.

194

Both k0 and k1 behaved in a similar way as kmax for all SF samples (Figure 2A), which indicates that both the

195

first and zero order mineralization rate determined the maximum mineralization rate. For all SF samples except

196

Br and StL, r equaled 0. For samples Br and StL, kinetic parameters estimated using the Three-half-order model

197

were not accurate, but in accordance with those calculated for pmax and kmax.

198

BAM and 2,6-DCBA mineralization capacity in soil samples. The frequent occurrence of 2,6-DCBA

199

mineralization in both BAM exposed and non-exposed SFs aroused our interest whether this capacity is also

200

common to other ecosystems. Therefore, a broad survey on BAM/2,6-DCBA mineralization capacity was

201

performed using 17 soils from various ecosystems. The only soil that showed substantial BAM mineralization

202

was the railway embankment sample from Rotselaar (R3). Mineralization was nevertheless slow, with an

203

average pmax and kmax of 24.5% and 0.72%.d-1, respectively (Table 2,

204

Figure 1A). On the other hand, all soil samples tested positive for 2,6-DCBA mineralization (Figure 1C).

205

Mineralization capacity varied, being highest in forest soils (pmax = 46 to 70%, kmax = 5 to 15%.d-1), followed by

206

grassland soils (pmax ≈ 54%, kmax ≈ 4%.d-1) and railway embankments, agricultural and garden soil (pmax = 23 to

207

45%, kmax = 1 to 5%.d-1) (Figure 2B). The dune sand sample showed lowest 2,6-DCBA mineralization capacity

208

with a pmax < 20% and kmax ≈ 0.5%.d-1. One of the forest soil samples (R1) differed from the other forest

209

samples, with a pmax and kmax of only 34% and 2%.d-1, respectively, hence being comparable to railway

210

embankments and agricultural soils.

211

In contrast to the SF samples, k0-values for 2,6-DCBA mineralization in the soil samples were all in the same

212

range, i.e., 0.07 to 0.15 (Figure 2B) and pmax > p1. The zero order mineralization process thus contributed to the

213

maximum mineralization extent. For mineralization curves without lag phase, the linear growth rate (r) equaled

214

0, as was the case for the forest soils (except R1), four of the five railway embankment soils and the agricultural

215

soil StT. In all other cases, the linear growth term contributed to the kmax-values, since these were higher than

216

expected based on the corresponding k1-values. Kinetic parameters for soil BL could not accurately be

217

estimated with the Three-half-order model.

218

Partial least squares regression (PLSR). PLSR was performed to determine whether 2,6-DCBA

219

mineralization kinetics correlated with characteristics of either the SF or soil samples. In case of SF samples,

220

pmax, p1, kmax and k0 were successfully calibrated to DWTP/SF characteristics as shown by high R²-values (>0.8) ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 19

221

and low residual mean square errors (RMSE) (Table S4). However, based on results for cross validation, none

222

of the kinetic parameters could be predicted using DWTP/SF characteristics (R² < 0.5, high RMSE) (Table S4).

223

To evaluate which characteristics correlated with mineralization, the regression vector was analyzed for p1, pmax,

224

kmax and k1 (Figure 3). For all four parameters, 2,6-DCBA mineralization correlated positively with exposure to

225

BAM and other micropollutants and with the presence of NH4, alternative carbon sources (NPOC), Fe and Mn.

226

A negative correlation was found with the average temperature of the intake water.

227

Only pmax, p1 and kmax were successfully calibrated to soil characteristics but, in contrast to the SF model, the

228

model could be used to predict pmax and p1 (R² > 0.6 for cross validation) (Table S4). Results for the regression

229

vector (Figure 3) were comparable for pmax, p1 and kmax and indicated that the most important soil characteristics

230

determining 2,6-DCBA mineralization were pH (negative correlation), UV-absorbance (A254, positive

231

correlation) and texture. With regards to the latter, the most important fraction was clay, of which the

232

percentage correlated positively with mineralization.

233

A similar multivariate analysis for BAM mineralization could not be performed since the sample population

234

sizes were insufficiently large as only few BAM mineralization curves could accurately be described with the

235

Three-half-order model and only few SF/soil samples showed substantial BAM mineralization.

236

Mineralization of 2,6-DCBA produced by BAM degrading bacteria in bioaugmented sand filters. The

237

results reported above show that in contrast to 2,6-DCBA, BAM mineralization in SFs occurs rarely and, if so,

238

proceeds slowly. A solution to establish BAM mineralization in SFs is bioaugmentation with bacteria such as

239

Aminobacter sp. MSH1, but this might lead to the release of 2,6-DCBA, being a bottleneck for BAM

240

mineralization in MSH1.19 Therefore, we determined the fate of 2,6-DCBA released in the SF after introducing

241

strain M6.100g, an MSH1 mutant which converts BAM to 2,6-DCBA but is deficient in its further

242

mineralization, imposing as such a worst case scenario. M6.100g was added to SF sample K4 in a BAM

243

mineralization assay. The mineralization curve was compared with those obtained for non-inoculated K4

244

samples amended with 14C-BAM or 14C-2,6-DCBA (Figure 4). Adding strain M6.100g increased kmax from 3 to

245

≈ 8%.d-1, but did not reach the 2,6-DCBA mineralization rate (≈ 13%.d-1) (Figure S5). Similar effects were

246

found for k1. pmax after 33 days and k0 were comparable for all three mineralization curves. The estimated r-

247

values of the bioaugmented BAM mineralization assays did not significantly differ from those for BAM

248

mineralization without bioaugmentation, whereas for 2,6-DCBA mineralization without bioaugmentation, r

249

equaled 0.

ACS Paragon Plus Environment

Page 9 of 19

Environmental Science & Technology

250

Discussion

251

Mineralization of BAM and 2,6-DCBA was evaluated in SF samples from DWTPs differing in intake water

252

BAM contamination and, for comparison, in topsoil samples from different locations. Slow to intermediate

253

BAM mineralization was found for three of the five DWTPs which treated BAM contaminated water. The high

254

BAM mineralization capacity found in samples from Kluizen was consistent in time while none of the non-

255

exposed SF samples mineralized BAM. Furthermore, BAM was mineralized in only one soil sample (R3), taken

256

from a railway embankment located in an industrial area. Although we do not have records of pesticide use at

257

the locations of the soil samples, dichlobenil, the parent compound of BAM, has often been used for weed

258

control on railways, hence indicating plausible pre-exposure of soil R3 to BAM. As such, our results for both

259

SF and soil samples are in accordance with previous studies that report BAM degradation/mineralization as rare

260

and only occasionally found in pre-exposed environments.13–17 The occurrence of BAM mineralization might be

261

the result of selective conditions for adaptation of microbiota to degrade and even metabolize BAM. However,

262

exposure to BAM does not assure the occurrence of BAM mineralization since certain BAM exposed SFs (from

263

De Blankaart and Zele) did not show mineralization. In case of the SF samples from Kluizen and of soil sample

264

R3, mineralization curves were sigmoidal with lag phases of several days before significant BAM

265

mineralization was recorded. Although we have no evidence, this might indicate that BAM mineralization is

266

growth linked in those ecosystems. On the other hand, the long lag phases suggest that initial BAM degrading

267

populations are relatively small. Our study provides the first report of BAM mineralization capacity in an

268

environment different from soil and is of particular interest for DWTPs since BAM is a frequent groundwater

269

contaminant.3,4

270

In contrast to BAM, all SF samples mineralized 2,6-DCBA, showing that 2,6-DCBA mineralization is more

271

distributed in SFs than BAM mineralization. 2,6-DCBA mineralization capacity was highest for DWTPs

272

treating BAM contaminated water, with average maximum rates (kmax) up to 16%.d-1. That the SF samples

273

which showed BAM mineralization also showed 2,6-DCBA mineralization suggests that the observed BAM

274

mineralization pathway in the SFs likely proceeds through 2,6-DCBA. The results further indicate that, as

275

suggested for BAM, 2,6-DCBA mineralization depends on exposure to BAM, which was supported by PLSR

276

analysis. This might be due to the potential presence of 2,6-DCBA as a BAM transformation product in the

277

intake water or produced in the SF units, and agrees with Holtze et al.14 who reported that 2,6-DCBA was only

278

degraded in soils pre-exposed to dichlobenil. However, high 2,6-DCBA mineralization rates (≈ 6%.d-1) were

279

also recorded in samples of two SFs receiving pristine groundwater (Hch and Lo). The other two pristine

280

samples (StL2 and Br) showed 2,6-DCBA mineralization as well, albeit slow. Other factors than exposure to

281

BAM/2,6-DCBA might thus be key or contribute to 2,6-DCBA mineralization. That exposure to

282

dichlobenil/BAM is not a prerequisite to establish 2,6-DCBA mineralization is further supported by the results ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 19

283

obtained with the soil samples. These show that environmental prevalence of 2,6-DCBA mineralization is high

284

since all sampled soils mineralized 2,6-DCBA, with highest mineralization found for forest soils (Li2, SL, He1

285

and KL1). Previous reports on 2,6-DCBA biodegradation capacity did not always find degradation but never

286

included soils from forest environments.14,15 On the other hand, Fulthorpe et al.32 found that the capacity to

287

mineralize another chlorinated benzoic acid, i.e., 3-chlorobenzoate (3-CBA), was high and omnipresent in forest

288

and woodland soils. PLSR analysis indicated that 2,6-DCBA mineralization in the soil samples, as for the SF

289

samples, was related with high organic carbon content and in particular aromatic compounds (measured as UV-

290

absorbance). Likewise, Larsson et al.33 reported that several chlorinated aromatic xenobiotics were more easily

291

mineralized by microbial communities from humic-rich lake water, than by communities from a clear-water

292

lake. Fulthorpe et al.,32 found no correlation between the 3-CBA mineralization rate and any of the measured

293

soil characteristics but suggested that the high 3-CBA mineralization capacity in forest soils might be attributed

294

to high levels of organic compounds like aromatics, naturally present as components of plant litter, and

295

organochlorine compounds, for instance produced by forest litter decomposing fungi.32,34–36 Organisms

296

metabolizing these substrates might possess broad substrate enzymes that degrade analogous compounds such

297

as 2,6-DCBA. In case such organisms are highly abundant among the microbial community, this might explain

298

for the observed high 2,6-DCBA mineralization activity. Interestingly, PLSR analysis showed that, also for the

299

SF samples, 2,6-DCBA mineralization correlated with the intake water organic carbon content. Hence high

300

natural microbial activity towards potentially present aromatic compounds might also explain the widely

301

distributed 2,6-DCBA mineralization capacity in SFs. Otherwise, organisms present in microbial communities

302

able to metabolize (chlorinated) aromatic compounds might be prone to more rapid adaptation to

303

metabolize/degrade analogous anthropogenic compounds like 2,6-DCBA when they become available in the

304

ecosystem, as might be the case in SFs exposed to BAM.

305

Most of the 2,6-DCBA mineralization curves for SF and soil samples showed no lag phase, represented by

306

linear growth rates (r) equal to 0. This indicates that the degrading community did not require an initial

307

adaptation/growth phase to establish

308

mineralization kinetics do not provide information whether or not 2,6-DCBA mineralization might support

309

growth, but also do not exclude it. On the other hand, the absence of a lag phase indicates that the initial number

310

of 2,6-DCBA degrading organisms was relatively high. Based on comparison with results for BAM

311

mineralization by Aminobacter sp. MSH1 in SF material, the number of 2,6-DCBA mineralizing cells should at

312

least be 105 cells/g to achieve lag times shorter than one day (Vandermaesen et al., unpublished results).

313

Our observations that BAM and 2,6-DCBA can be mineralized in SF material opens windows for future

314

bioremediation strategies in DWTPs. However, taking into account (i) the relatively low BAM mineralization

315

rates (max. 2.8%.d-1 for sample K4) (ii) the low numbers of BAM degraders, as discussed above and (iii) the

316

high fluxes and short hydraulic retention times operated in rapid SFs, intrinsic BAM removal is expected to be

14

CO2 production from the added

ACS Paragon Plus Environment

14

C-2,6-DCBA. As such, the

Page 11 of 19

Environmental Science & Technology

317

insufficient in actual SFs to reach EU threshold concentrations. In contrast, 2,6-DCBA was readily mineralized

318

in most SF samples and the actual numbers of intrinsic 2,6-DCBA degraders appear relatively high. Assuming

319

that (i) 2,6-DCBA degradation is metabolic, (ii) the ratio of 14CO2 production to 14C incorporated into biomass

320

is constant throughout the experiment and (iii) 2,6-DCBA is completely degraded when mineralization reaches

321

its plateau (thus the remaining

322

(sample E3) corresponds to 100% 2,6-DCBA degradation and the kmax-value of 16%.d-1 (sample E3)

323

corresponds to a degradation rate of 1.6 µg/L.h at the used initial 2,6-DCBA concentration of 150 µg/L. To

324

mineralize a 2,6-DCBA concentration of 0.14 µg/L, corresponding to the maximum intake water BAM

325

concentration found in this study (Table 1), to below the threshold concentration of 0.1 µg/L in a rapid SF with

326

a flow rate of 20 m/h and a filter depth of 1 m, a mineralization rate of 8%.d-1 or 0.8 µg/L.h is then needed,

327

which is lower than found for sample E3. Actual mineralization rates in operational SFs are however expected

328

to be negatively influenced by regular backwashing,37 lower temperatures38,39 and lower pesticide

329

concentrations.16,25,40 Although the extent to which this affects 2,6-DCBA mineralization cannot be predicted

330

based on the currently available data and although the mechanism of mineralization is unsure, the results are

331

promising in order to achieve complete BAM degradation after bioaugmentation of SFs with Aminobacter sp.

332

MSH1, especially at lower flow rates such as those used in slow SFs (0.4 m/h).

333

Bioaugmentation of SFs with BAM degrading bacterial strains such as Aminobacter sp. MSH1 seems a

334

promising strategy for BAM removal in DWTPs.18 A possible drawback could however be the production of

335

2,6-DCBA after the first step in BAM degradation, as its further mineralization has been described to be rate

336

limiting.19,20 The high 2,6-DCBA mineralization found in SF samples from DWTPs treating BAM contaminated

337

water though suggests that any 2,6-DCBA produced during BAM degradation is rapidly degraded by the

338

endogenous community, as shown in our bioaugmentation assay with strain M6.100g. As such, our observations

339

are also of major interest for bioaugmentation approaches for BAM removal in DWTPs using MSH1. In case

340

BAM is mineralized in SFs via 2,6-DCBA, conversion of BAM to 2,6-DCBA appears the rate limiting reaction.

341

This is in agreement with studies on BAM mineralization in environmental soil samples by Clausen et al.,13 but

342

in contrast to what has been described for MSH1, i.e., degradation of 2,6-DCBA as the bottleneck for BAM

343

mineralization.19,20

14

C is also degraded and incorporated into biomass), then a pmax-value of 61%

344

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 19

345

Supporting information available

346

Supporting Information Available: Overview of SF and soil characteristics and locations; details, data used for

347

and performance of PLSR analysis; a comparison of the Three-half-order and First-zero-order model; BAM

348

mineralization curves and kinetics for SF samples taken at different time points and kinetic parameters from the

349

bioaugmentation assay. This information is available free of charge via the Internet at http://pubs.acs.org.

350

Acknowledgements

351

This study was funded by the FP7 project BIOTREAT (EU grant n° 266039), the Inter-University Attraction

352

Pole (IUAP) “µ-manager” of the Belgian Science Policy (BELSPO, P7/25) and the “Fonds Wetenschappelijk

353

Onderzoek” (FWO) post-doctoral fellow grant n° 12Q0215N to B. Horemans. We thank S. Goethals and R.

354

Jenné for assistance in SF sampling and providing DWTP process information, and K. Simoens, D. Grauwels,

355

K. Moors, J. Plevoets, L. Fondu, A. Deckers, A. Vertommen and N. Croonenborghs for assistance in the

356

experimental work.

357

References

358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380

(1)

(2)

(3)

(4)

(5)

(6)

(7)

Björklund, E.; Anskjaer, G. G.; Hansen, M.; Styrishave, B.; Halling-Sorensen, B. Analysis and environmental concentrations of the herbicide dichlobenil and its main metabolite 2,6dichlorobenzamide (BAM): A review. Sci. Total Environ. 2011, 409 (12), 2343–2356 DOI: 10.1016/j.scitotenv.2011.02.008. Clausen, L.; Larsen, F.; Albrechtsen, H. J. Sorption of the herbicide dichlobenil and the metabolite 2,6dichlorobenzamide on soils and aquifer sediments. Environ. Sci. Technol. 2004, 38 (17), 4510–4518 DOI: 10.1016/j.jconhyd.2006.04.004. Törnquist, M.; Kreuger, J.; Adielsson, S. Occurrence of pesticides in Swedish water resources against a background of national risk-reduction programmes - Results from 20 years of monitoring. In XIII Symposium Pesticide Chemistry - Environmental Fate and Human Health; 2007. VMM. Pesticiden in het grondwater in Vlaanderen (in Dutch); Vlaamse Milieumaatschappij: Aalst, Belgium, 2012; https://www.vlaanderen.be/nl/publicaties/detail/pesticiden-in-het-grondwater-invlaanderen. GEUS. Groundwater monitoring 2013 — groundwater status and development 1989–2012 (in Danish); The Geological Survey of Denmark and Greenland, Ministry of Climate, Energy and Building: Copenhagen, Denmark, 2013; http://www.geus.dk/DK/publications/groundwater_monitoring/Sider/1989_2012.aspx. Westerhoff, P.; Yoon, Y.; Snyder, S.; Wert, E. Fate of endocrine-disruptor, pharmaceutical, and personal care product chemicals during simulated drinking water treatment processes. Environ. Sci. Technol. 2005, 39 (17), 6649–6663 DOI: 10.1021/es0484799. Albers, C. N.; Ellegaard-Jensen, L.; Harder, C. B.; Rosendahl, S.; Knudsen, B. E.; Ekelund, F.; Aamand, J. Groundwater chemistry determines the prokaryotic community structure of waterworks sand filters. Environ. Sci. Technol. 2015, 49 (2), 839–846 DOI: 10.1021/es5046452. ACS Paragon Plus Environment

Page 13 of 19

381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432

(8) (9) (10)

(11) (12)

(13)

(14)

(15)

(16)

(17)

(18)

(19)

(20)

(21)

(22)

(23) (24) (25)

Environmental Science & Technology

Proctor, C. R.; Hammes, F. Drinking water microbiology — from measurement to management. Curr. Opin. Biotechnol. 2015, 33 (0), 87–94 DOI: 10.1016/j.copbio.2014.12.014. Hedegaard, M. J.; Albrechtsen, H.-J. Microbial pesticide removal in rapid sand filters for drinking water treatment – Potential and kinetics. Water Res. 2014, 48, 71–81 DOI: 10.1016/j.watres.2013.09.024. Hedegaard, M. J.; Arvin, E.; Corfitzen, C. B.; Albrechtsen, H.-J. Mecoprop (MCPP) removal in fullscale rapid sand filters at a groundwater-based waterworks. Sci. Total Environ. 2014, 499, 257–264 DOI: 10.1016/j.scitotenv.2014.08.052. Zearley, T. L.; Summers, R. S. Removal of trace organic micropollutants by drinking water biological filters. Environ. Sci. Technol. 2012, 46 (17), 9412–9419 DOI: 10.1021/es301428e. Zuehlke, S.; Duennbier, U.; Heberer, T. Investigation of the behavior and metabolism of pharmaceutical residues during purification of contaminated ground water used for drinking water supply. Chemosphere 2007, 69 (11), 1673–1680 DOI: 10.1016/j.chemosphere.2007.06.020. Clausen, L.; Arildskov, N. P.; Larsen, F.; Aamand, J.; Albrechtsen, H. J. Degradation of the herbicide dichlobenil and its metabolite BAM in soils and subsurface sediments. J. Contam. Hydrol. 2007, 89 (34), 157–173 DOI: 10.1016/j.jconhyd.2006.04.004. Holtze, M. S.; Hansen, H. C. B.; Juhler, R. K.; Sørensen, J.; Aamand, J. Microbial degradation pathways of the herbicide dichlobenil in soils with different history of dichlobenil-exposure. Environ. Pollut. 2007, 148 (1), 343–351 DOI: 10.1016/j.envpol.2006.10.028. Holtze, M. S.; Sorensen, S. R.; Sorensen, J.; Aamand, J. Microbial degradation of the benzonitrile herbicides dichlobenil, bromoxynil and ioxynil in soil and subsurface environments–insights into degradation pathways, persistent metabolites and involved degrader organisms. Environ. Pollut. 2008, 154 (2), 155–168 DOI: 10.1016/j.envpol.2007.09.020. Janniche, G. S.; Clausen, L.; Albrechtsen, H. J. Inherent mineralization of 2,6-dichlorobenzamide (BAM) in unsaturated zone and aquifers–effect of initial concentrations and adaptation. Environ. Pollut. 2011, 159 (10), 2801–2807 DOI: 10.1016/j.envpol.2011.05.010. Simonsen, A.; Holtze, M. S.; Sorensen, S. R.; Sorensen, S. J.; Aamand, J. Mineralisation of 2,6dichlorobenzamide (BAM) in dichlobenil-exposed soils and isolation of a BAM-mineralising Aminobacter sp. Environ. Pollut. 2006, 144 (1), 289–295 DOI: 10.1016/j.envpol.2005.11.047. Albers, C. N.; Feld, L.; Ellegaard-Jensen, L.; Aamand, J. Degradation of trace concentrations of the persistent groundwater pollutant 2,6-dichlorobenzamide (BAM) in bioaugmented rapid sand filters. Water Res. 2015, 83, 61–70 DOI: 10.1016/j.watres.2015.06.023. Simonsen, A.; Badawi, N.; Anskjær, G. G.; Albers, C. N.; Sorensen, S. R.; Sorensen, J.; Aamand, J. Intermediate accumulation of metabolites results in a bottleneck for mineralisation of the herbicide metabolite 2,6-dichlorobenzamide (BAM) by Aminobacter spp. Appl. Microbiol. Biotechnol. 2012, 94 (1), 237–245 DOI: 10.1007/s00253-011-3591-x. T’Syen, J.; Tassoni, R.; Hansen, L.; Sorensen, S. J.; Leroy, B.; Sekhar, A.; Wattiez, R.; De Mot, R.; Springael, D. Identification of the amidase BbdA that initiates biodegradation of the groundwater micropollutant 2,6-dichlorobenzamide (BAM) in Aminobacter sp. MSH1. Environ. Sci. Technol. 2015, 49 (19), 11703–11713 DOI: 10.1021/acs.est.5b02309. Fava, L.; Orrù, M. A.; Crobe, A.; Caracciolo, A. B.; Bottoni, P.; Funari, E. Pesticide metabolites as contaminants of groundwater resources: assessment of the leaching potential of endosulfan sulfate, 2,6dichlorobenzoic acid, 3,4-dichloroaniline, 2,4-dichlorophenol and 4-chloro-2-methylphenol. Microchem. J. 2005, 79 (1–2), 207–211 DOI: 10.1016/j.microc.2004.10.009. GEUS. Grundvand, Status og udvikling 1989 – 2014 (in Danish); The Geological Survey of Denmark and Greenland, Ministry of Climate, Energy and Building: Copenhagen, Denmark, 2015; http://www.geus.dk/DK/publications/groundwater_monitoring/Sider/1989_2014.aspx. Sheets, T. J.; Harris, C. I.; Smith, J. W. Persistence of Dichlobenil and SD-7961 in Soil. Weed Sci. 1968, 16 (2), 245–249. Verloop, A. Residue reviews: Residues of pesticides and other contaminants in the total environment; Gunther, F. A., Gunther, J. D., Eds.; Springer New York: New York, NY, 1972; pp 55–103. Sørensen, S. R.; Holtze, M. S.; Simonsen, A.; Aamand, J. Degradation and mineralization of nanomolar concentrations of the herbicide dichlobenil and its persistent metabolite 2,6-dichlorobenzamide by ACS Paragon Plus Environment

433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474

Environmental Science & Technology

(26) (27)

(28)

(29) (30) (31)

(32)

(33) (34)

(35) (36)

(37) (38) (39) (40)

Page 14 of 19

Aminobacter spp. isolated from dichlobenil-treated soils. Appl. Environ. Microbiol. 2007, 73 (2), 399– 406 DOI: 10.1128/AEM.01498-06. Reasoner, D. J.; Geldreich, E. E. A new medium for the enumeration and subculture of bacteria from potable water. Appl. Environ. Microbiol. 1985, 49 (1), 1–7. Dejonghe, W.; Berteloot, E.; Goris, J.; Boon, N.; Crul, K.; Maertens, S.; Höfte, M.; De Vos, P.; Verstraete, W.; Top, E. M. Synergistic degradation of linuron by a bacterial consortium and isolation of a single linuron-degrading Variovorax strain. Appl. Environ. Microbiol. 2003, 69 (3), 1532–1541 DOI: 10.1128/AEM.69.3.1532-1541.2003. Sniegowski, K.; Bers, K.; Ryckeboer, J.; Jaeken, P.; Spanoghe, P.; Springael, D. Robust linuron degradation in on-farm biopurification systems exposed to sequential environmental changes. Appl. Environ. Microbiol. 2011, 77 (18), 6614–6621 DOI: 10.1128/AEM.05108-11. Brunner, W.; Focht, D. D. Deterministic three-half-order kinetic model for microbial degradation of added carbon substrates in soil. Appl Environ Microbiol 1984, 47 (1), 167–172. Fomsgaard, I. S. Modelling the mineralization kinetics for low concentrations of pesticides in surface and subsurface soil. Ecol. Modell. 1997, 102 (2-3), 175–208 DOI: 10.1016/S0304-3800(97)01982-0. Amery, F.; Degryse, F.; Degeling, W.; Smolders, E.; Merckx, R. The copper-mobilizing-potential of dissolved organic matter in soils varies 10-fold depending on soil incubation and extraction procedures. Environ Sci Technol 2007, 41 (7), 2277–2281 DOI: 10.1021/es062166r. Fulthorpe, R. R.; Rhodes, A. N.; Tiedje, J. M. Pristine soils mineralize 3-chlorobenzoate and 2,4dichlorophenoxyacetate via different microbial populations. Appl. Environ. Microbiol. 1996, 62 (4), 1159–1166. Larsson, P.; Okla, L.; Tranvik, L. Microbial degradation of xenobiotic, aromatic pollutants in humic water. Appl. Environ. Microbiol. 1988, 54 (7), 1864–1867. Helfrich, M.; Ludwig, B.; Buurman, P.; Flessa, H. Effect of land use on the composition of soil organic matter in density and aggregate fractions as revealed by solid-state 13C NMR spectroscopy. Geoderma 2006, 136 (1–2), 331–341 DOI: 10.1016/j.geoderma.2006.03.048. Marschner, B.; Kalbitz, K. Controls of bioavailability and biodegradability of dissolved organic matter in soils. Geoderma 2003, 113 (3–4), 211–235 DOI: 10.1016/S0016-7061(02)00362-2. Verhagen, F. J. M.; Swarts, H. J.; Wunberg, J. B. P. A.; Field, J. A. Organohalogen production is a ubiquitous capacity among Basidiomycetes. Chemosphere 1998, 37 (9–12), 2091–2104 DOI: 10.1016/S0045-6535(98)00272-0. Hozalski, R. M.; Bouwer, E. J. Deposition and retention of bacteria in backwashed filters. J. - Am. Water Works Assoc. 1998, 90 (1), 71. Nair, D. R.; Schnoor, J. L. Effect of soil conditions on model parameters and atrazine mineralization rates. Water Res. 1994, 28 (5), 1199–1205 DOI: 10.1016/0043-1354(94)90208-9. del Pilar Castillo, M.; Torstensson, L. Effect of biobed composition, moisture, and temperature on the degradation of pesticides. J. Agric. Food Chem. 2007, 55 (14), 5725–5733 DOI: 10.1021/jf0707637. Lipthay, J. R. de; Sørensen, S. R.; Aamand, J. Effect of herbicide concentration and organic and inorganic nutrient amendment on the mineralization of mecoprop, 2,4-D and 2,4,5-T in soil and aquifer samples. Environ. Pollut. 2007, 148 (1), 83–93 DOI: http://dx.doi.org/10.1016/j.envpol.2006.11.005.

ACS Paragon Plus Environment

Page 15 of 19

475 476 477

Environmental Science & Technology

Table 1. Overview of sand filter samples used in this study with kinetic data for recorded BAM and 2,6-DCBA mineralization. Average values ± standard deviation are given for three replicates.

SF Sample K1 K2 K3 K4 DB AWW E1 E2 E3 Z Si Sn1 Sn2 StL Br Hch Lo

Sampling date 8/11/2011 12/06/2012 3/04/2013 6/05/2014 27/06/2012 13/06/2012 9/11/2011 13/06/2012 27/11/2014 28/06/2012 28/06/2012 28/06/2012 19/11/2014 17/11/2014 4/03/2015 5/05/2015 7/07/2015

Location Kluizen Kluizen Kluizen Kluizen De Blankaart Antwerp Eeklo Eeklo Eeklo Zele Klein-Sinaai Snellegem Snellegem Saint Léger Bree Haacht Lommel

Type of intake water1 SW + GW SW + GW SW + GW SW + GW SW SW GW GW GW GW GW GW GW GW GW GW GW

BAM exposure (µg/L)2 0.09 0.09 0.09 0.09 0.10 unknown 0.02 0.02 0.02 0.06 0.00 0.14 0.14 0.00 0.00 0.00 0.00

BAM mineralization3 pmax (%) kmax (%.d-1) 33.9 ± 3.74 0.98 ± 0.18 62.4 ± 6.89 2.92 ± 0.47 62.3 ± 3.39 3.24 ± 0.13 63.2 ± 11.9 2.79 ± 0.49 1.57 ± 0.62 0.09 ± 0.04 4.44 ± 1.42 0.20 ± 0.05 5.89 ± 0.65 0.33 ± 0.04 6.99 ± 1.27 0.26 ± 0.10 10.6 ± 1.47 0.31 ± 0.07 1.20 ± 0.16 0.22 ± 0.06 2.08 ± 0.68 0.09 ± 0.00 7.82 ± 5.94 0.55 ± 0.41 4.82 ± 0.58 0.21 ± 0.01 1.04 ± 0.04 0.03 ± 0.00 0.40 ± 0.13 0.02 ± 0.00 2.60 ± 0.19 0.10 ± 0.02 4.57 ± 0.35 0.12 ± 0.02

2,6-DCBA mineralization3,4 pmax (%) kmax (%.d-1) n.d. n.d. n.d. n.d. n.d. n.d. 50.8 ± 9.49 11.64 ± 2.25 n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. 61.2 ± 5.31 15.68 ± 5.69 n.d. n.d. n.d. n.d. n.d. n.d. 62.7 ± 0.74 12.9 ± 1.43 5.22 ± 1.69 0.15 ± 0.05 5.45 ± 4.80 0.22 ± 0.24 35.0 ± 3.45 7.09 ± 0.99 26.8 ± 3.30 5.45 ± 1.43

1

SW = surface water, GW = groundwater Average BAM concentration in the DWTP intake water 3 pmax = extent of mineralization after 50 days of incubation, kmax = maximum linear mineralization rate 4 n.d. = not done 2

478 479 480 481

Table 2. Overview of soil samples used in this study with kinetic data for recorded BAM and 2,6-DCBA mineralization. Average values ± standard deviation are given for three replicates.

BAM mineralization1 Soil Sampling sample date Location Land use pmax (%) kmax (%.d-1) BL 25/05/2015 Blankenberge Sand dune 3.54 ± 0.38 0.11 ± 0.02 Ha 23/03/2015 Halen Agriculture: maize, potatoe 2.54 ± 0.19 0.11 ± 0.03 He1 16/03/2015 Heverlee Deciduous forest 1.69 ± 0.47 0.17 ± 0.04 He2 27/05/2015 Heverlee Railway embankment, domestic area 0.65 ± 0.31 0.05 ± 0.03 KL1 2/08/2015 Kessel-Lo Deciduous forest 2.22 ± 0.19 0.13 ± 0.03 KL2 2/08/2015 Kessel-Lo Railway embankment, domestic area 1.37 ± 0.28 0.08 ± 0.01 Li1 15/03/2015 Linkhout Grassland, historical: coniferous forest 0.76 ± 0.11 0.08 ± 0.01 Li2 23/03/2015 Linkhout Mixed forest 0.77 ± 0.09 0.17 ± 0.01 R1 14/03/2015 Rotselaar Coniferous forest 1.54 ± 0.08 0.12 ± 0.01 R2 14/03/2015 Rotselaar Domestic garden, riverbank 1.55 ± 0.30 0.13 ± 0.02 R3 30/05/2015 Rotselaar Railway embankment, industrial area 24.5 ± 2.92 0.72 ± 0.09 SB1 15/03/2015 Schoonderbuken Grassland, historical: agriculture 2.89 ± 0.18 0.13 ± 0.01 SB2 15/03/2015 Schoonderbuken Agriculture: maize, sugar beet, wheat 1.29 ± 0.59 0.07 ± 0.02 Sch 26/06/2015 Schulen Railway embankment, remote area 1.57 ± 0.13 0.12 ± 0.01 SL 14/03/2015 Somme-Leuze Mixed forest 4.45 ± 0.78 0.22 ± 0.03 StT 12/03/2015 Sint-Truiden Agriculture: pear orchard 2.45 ± 0.38 0.18 ± 0.01 Zm 26/06/2015 Zelem Railway embankment, remote area 3.00 ± 0.31 0.11 ± 0.01 1 pmax = extent of mineralization after 50 days of incubation, kmax = maximum linear mineralization rate

482

ACS Paragon Plus Environment

2,6-DCBA mineralization1 pmax (%) kmax (%.d-1) 15.7 ± 6.21 0.50 ± 0.06 28.7 ± 1.69 1.15 ± 0.29 51.2 ± 4.60 7.13 ± 1.21 31.9 ± 3.21 2.29 ± 0.29 54.5 ± 3.28 8.52 ± 0.55 32.7 ± 2.24 2.09 ± 0.34 54.6 ± 4.54 4.81 ± 1.28 63.2 ± 7.55 9.54 ± 3.61 33.9 ± 1.90 2.39 ± 0.37 42.5 ± 2.44 2.04 ± 0.39 32.0 ± 1.23 1.23 ± 0.11 52.2 ± 2.68 3.04 ± 0.42 26.1 ± 3.66 0.99 ± 0.23 36.8 ± 3.65 2.58 ± 0.23 60.8 ± 3.57 12.6 ± 2.45 36.8 ± 1.51 3.86 ± 0.10 35.6 ± 4.90 4.38 ± 0.89

Environmental Science & Technology

Page 16 of 19

483 484 485 486 487 488 489

Figure 1. Cumulative BAM mineralization curves (A) for sand filter samples (full lines) and soil sample R3 (dashed line), 2,6DCBA mineralization curves for sand filter samples (B) from DWTPs treating water contaminated with BAM (black) and pristine groundwater (grey) and 2,6-DCBA mineralization curves in soil samples (C) from forests (green), agricultural fields (red), railway embankments (brown) and grasslands, garden soil and a sand dune (blue). Mean values and standard deviations based on three replicates are shown for 14CO2 production relative to the initial amount of 14C-BAM/2,6-DCBA added (14C0).

ACS Paragon Plus Environment

Page 17 of 19

Environmental Science & Technology

490 491 492 493 494 495 496 497

Figure 2. Kinetic parameters of 2,6-DCBA mineralization in sand filter samples (A) from DWTPs treating water contaminated with BAM (black) and pristine groundwater (grey) and in soil samples (B) from forests (green), agricultural fields (red), railway embankments (brown) and grasslands, garden soil and a sand dune (blue): extent of mineralization after 50 days (pmax), maximum linear mineralization rate (kmax), extent of mineralization by the first order process (p1), first order mineralization rate (k1), zero order mineralization rate (k0) and linear growth rate (r). Error bars indicate 95% confidence intervals for kinetic parameters returned by the Three-half-order model (P1, k0, k1 and r).

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 19

498 499 500 501 502 503 504 505 506

Figure 3. Regression vectors for DWTP/SF characteristics (left) and soil characteristics (right) correlated with the extent of mineralization by the first order process (p1, black), the extent of mineralization after 50 days (pmax, grey) and the maximum linear mineralization rate (kmax, white) for 2,6-DCBA mineralization in soil/SF samples. For the SF model, results are only shown for p1 because of overlap with data points for pmax, kmax and the zero order mineralization rate (k0). The regression vector returns positive values for characteristics that correlate positively with mineralization and negative values for those that correlate negatively. The absolute value relates to the importance of the characteristic in the regression model. A detailed overview of the DWTP/SF and soil characteristics used is provided in Table S1 and Table S2, respectively.

507 508 509 510

Figure 4. Cumulative mineralization curves for BAM (black) and 2,6-DCBA (grey) with (full symbols) and without (open symbols) bioaugmentation with strain Aminobacter sp. M6.100g in sand filter sample K4. Mean values and standard deviations based on three replicates are shown for 14CO2 production relative to the initial amount of 14C-BAM/2,6-DCBA added (14C0).

ACS Paragon Plus Environment

Page 19 of 19

Environmental Science & Technology

1037x580mm (96 x 96 DPI)

ACS Paragon Plus Environment