Anna
Jaeger‡
*ab,
Claudia
Coll‡
*c,
Malte
Posselt
c,
Jonas
Mechelke
de,
Cyrus
Rutere
f,
Andrea
Betterle
dg,
Muhammad
Raza
hi,
Anne
Mehrtens
j,
Karin
Meinikmann
a,
Andrea
Portmann
k,
Tanu
Singh
l,
Phillip J.
Blaen
lm,
Stefan
Krause
l,
Marcus A.
Horn
fn,
Juliane
Hollender
de,
Jonathan P.
Benskin
c,
Anna
Sobek§
c and
Joerg
Lewandowski§
ab
aDepartment Ecohydrology, Leibniz-Institute of Freshwater Ecology and Inland Fisheries, Berlin, Germany
bGeography Department, Humboldt University Berlin, Berlin, Germany. E-mail: anna.jaeger@igb-berlin.de
cDepartment of Environmental Science and Analytical Chemistry (ACES), Stockholm University, Stockholm, Sweden. E-mail: claudia.coll@aces.su.se
dEawag, Swiss Federal Institute of Aquatic Science and Technology, Dübendorf, Switzerland
eInstitute of Biogeochemistry and Pollutant Dynamics, ETH Zürich, Zürich, Switzerland
fDepartment of Ecological Microbiology, University of Bayreuth, Bayreuth, Germany
gDepartment of ICEA and International Center for Hydrology “Dino Tonini”, University of Padova, Padua, Italy
hInstitute of Applied Geosciences, Technical University of Darmstadt, Darmstadt, Germany
iIWW Water Centre, Mülheim an der Ruhr, Germany
jInstitute of Biology and Environmental Sciences, Carl von Ossietzky University Oldenburg, Oldenburg, Germany
kCivil and Environmental Engineering, Colorado School of Mines, Golden, Colorado, USA
lSchool of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, UK
mYorkshire Water, Bradford, UK
nInstitute of Microbiology, Leibniz University of Hannover, Hannover, Germany
First published on 17th October 2019
Enhancing the understanding of the fate of wastewater-derived organic micropollutants in rivers is crucial to improve risk assessment, regulatory decision making and river management. Hyporheic exchange and sediment bacterial diversity are two factors gaining increasing importance as drivers for micropollutant degradation, but are complex to study in field experiments and usually ignored in laboratory tests aimed to estimate environmental half-lives. Flume mesocosms are useful to investigate micropollutant degradation processes, bridging the gap between the field and batch experiments. However, few studies have used flumes in this context. We present a novel experimental setup using 20 recirculating flumes and a response surface model to study the influence of hyporheic exchange and sediment bacterial diversity on half-lives of the anti-epileptic drug carbamazepine (CBZ) and the artificial sweetener acesulfame (ACS). The effect of bedform-induced hyporheic exchange was tested by three treatment levels differing in number of bedforms (0, 3 and 6). Three levels of sediment bacterial diversity were obtained by diluting sediment from the River Erpe in Berlin, Germany, with sand (1:
10, 1
:
1000 and 1
:
100
000). Our results show that ACS half-lives were significantly influenced by sediment dilution and number of bedforms. Half-lives of CBZ were higher than ACS, and were significantly affected only by the sediment dilution variable, and thus by bacterial diversity. Our results show that (1) the flume-setup is a useful tool to study the fate of micropollutants in rivers, and that (2) higher hyporheic exchange and bacterial diversity in the sediment can increase the degradation of micropollutants in rivers.
Environmental significanceContamination of rivers by wastewater-derived organic micropollutants is an emerging problem. Global-scale risk assessment and modeling, as well as regulatory decision making rely on in-depth knowledge of major factors that influence the fate of micropollutants in rivers. Sediment bacterial diversity and hyporheic exchange are parameters that are often neglected when testing persistence of substances. We present a novel experimental setup consisting of 20 recirculating flumes and a response surface model to test the influence of these two parameters on micropollutant half-lives. Our setup proved useful to study the fate of micropollutants in rivers. Our results reveal that the underrated factors, hyporheic exchange and bacterial diversity significantly affect micropollutant half-lives in rivers. |
Identification of key drivers of the attenuation of micropollutants in the environment has been attempted over the last few years in both field4–7 and lab studies.8–12 While the complex interplay of spatially heterogeneous processes and transient boundary conditions in field experiments hinders repeatability and often reduces the validity of general conclusions about mechanisms, the repeatable results obtained in lab experiments often lack diagnostic power for processes in real systems.13 Bottle incubations carried out according to OECD guidelines 308 and 309,14,15 for instance, are commonly used to test aerobic and anaerobic biodegradability of micropollutants in water–sediment systems. Compound specific half-lives resulting from these tests are used in current chemical legislation and for risk and exposure modeling.2,16 However, these stagnant tests do not cover processes specific for river systems. They disregard dynamic hydrological processes, such as hyporheic exchange fluxes and the characteristics of the site-specific bacterial communities.17–19 To this end, systematic tools that are more suitable for simulating conditions in rivers are necessary to study the fate of micropollutants.
Flumes are experimental mesocosms that can bridge the aforementioned gap between controllability of lab experiments and the authenticity of field experiments. Flumes have been proven valuable for investigating drivers of hyporheic exchange under controlled conditions20,21 and are thus suitable to simulate hydraulic conditions in rivers. Flume designs can vary in complexity, from simple recirculating systems to those with advanced flow-control functions which can be used to simulate gaining and losing conditions in the hyporheic zone.22 While more advanced designs are useful for testing complex hypotheses, they are clearly more challenging to replicate due to the need for additional infrastructure. In contrast, simple flume designs offer the opportunity to run experiments in parallel which enables testing of a wide range of scenarios and application of statistical models. Despite the apparent advantages of flume studies, to the best of our knowledge, there have been only a few studies investigating the fate of micropollutants in flumes. When Li et al.23 and Kunkel and Radke24 used a recirculating flume to test degradation of micropollutants, only one flume was available, which is why they repeated the experiment to test the difference between flat streambed and streambed with bedforms23 and two different flow velocities,24 respectively. However, no replicates and no statistical design were implemented.
Biodegradation is a key process for micropollutant transformation in rivers and is controlled by a variety of physical and biological factors. The hyporheic zone in particular represents a hot-spot of microbial activity.25,26 It has been described as a key compartment for nutrient turnover and other biochemical reactions in rivers.27,28 Schaper et al.29 found that reach-scale removal of non-persistent polar organic micropollutants in a South Australian river was controlled by attenuation in the hyporheic zone during the wet season. Hence, enhancing the potential for hyporheic exchange will likely promote micropollutant transformation. While hyporheic flow can occur on many scales,27 small scale exchange and short flow paths are expected to be the most effective in turnover processes.30,31 Small scale hyporheic exchange is determined by sediment characteristics, mainly hydraulic conductivity and sediment morphology.32 Bedforms such as dunes can induce hyporheic exchange causing a so called “pumping effect”. High pressure on the upstream side of a dune forces the surface water into the sediment and low pressure on the downstream slope of the dune causes exfiltration of hyporheic water to the surface water.20 While the extent of mass transport into the hyporheic zone and residence time distributions control the potential for turnover, the biochemical conditions control the quality of turnover. The composition of the hyporheic bacterial community is expected to play a major role in biochemical processes, but there has been little research on their impact on micropollutant degradation in rivers. Recent studies suggest that bacterial diversity in wastewater treatment plants is associated with biodegradation of certain micropollutants.33,34 It was proposed that this positive correlation is observed when the biodegradation mechanism of a certain micropollutant can be performed by only a few bacteria taxa. There is scarce information on the compound-specific catabolic and co-metabolic pathways that are required for the biodegradation of the variety of micropollutants, and whether these are general or rare functions in bacterial communities. Stadler et al.33 used a “dilution-to-extinction” approach to manipulate bacterial communities from activated sludge into different diversity levels and by this means identified taxa potentially driving the biotransformation of individual micropollutants. It was also shown that the change in taxonomic richness was associated with a change in functional richness across the treatment levels. In the absence of detailed knowledge on specific bacteria activity and functional characteristics, taxonomic diversity is a first step to understand the link between bacterial community composition and biodegradation of micropollutants.
In this study we present a novel experimental setup to test the influence of hyporheic flow and bacterial diversity on degradation half-lives of micropollutants. The experiment was based on a central composite face factorial design and used 20 flume mesocosms run in parallel to simulate different river conditions. The hyporheic flow was manipulated with the presence and number of bedforms. The bacterial diversity was controlled with a dilution-to-extinction approach.33 The flumes were inoculated with sediment from the River Erpe in Berlin, Germany, which receives high loads of treated wastewater.6 As the sediment of this river is continuously exposed to micropollutants, its bacterial community has likely developed the capacity to degrade synthetic organic chemicals. A response surface model was used to assess the effects of the two variables (i.e., sediment dilution and number of bedforms) on micropollutant dissipation half-lives (DT50s). The DT50s of the artificial sweetener acesulfame (ACS) and the anti-epileptic drug carbamazepine (CBZ) are discussed to evaluate the performance of the setup. These chemicals were chosen as model compounds due to their widespread occurrence in freshwaters, high average concentrations in River Erpe (μg L−1-range)6 and contrasting behavior in the environment. Although both compounds were previously reported as relatively persistent,35,36 recent research found increasing degradability of ACS within the last decade, likely caused by the adaptation of microbial communities in treatment plants.37 In addition, a study conducted in the hyporheic zone of River Erpe showed that along a vertical flow path of 40 cm, ACS was removed by 78 ± 1% while CBZ was not removed significantly.30 Therefore, we expect generally lower half-lives and a higher influence of the tested variables for ACS than for CBZ. The study aims at: providing a new experimental method to obtain a more accurate understanding of dependencies of micropollutant half-lives on river-specific biological and physical conditions and discussing the performance of the specific experimental setup.
Independent variables | Factors | Coded levels | ||
---|---|---|---|---|
−1 | 0 | 1 | ||
Hyporheic exchange – bedforms (B) | Level name | B0 | B3 | B6 |
Number of bedforms | 0 | 3 | 6 | |
Bacterial diversity – sediment dilution (S) | Level name | S6 | S3 | S1 |
Sediment![]() ![]() |
1![]() ![]() |
1![]() ![]() |
1![]() ![]() |
(1) Bacterial diversity in sediment (S) was achieved by diluting river sediment with commercial sand in different proportions, based on the dilution-to-extinction method. In this method, the less abundant species are “removed” by stepwise dilutions, and this loss in species richness results in lower diversity.33 The low dilution level (S1) had a sediment dilution of 1:
10 and was expected to have the highest level of bacterial diversity, the medium level (S3) was diluted in a 1
:
103 ratio, and the high dilution level (S6) corresponding to the lowest expected bacterial diversity, was set to a 1
:
106 dilution. Bacterial diversity in S1, S3 and S6 levels was investigated through Illumina sequences of the 16S rRNA taxonomic gene (details in Chapter 2.7).
(2) Hyporheic exchange (B) was regulated by forming triangular-shaped stationary bedforms in the flume sediment. Bedforms cause a so-called “pumping” effect, inducing advective flow through the porous streambed by pressure differences between the upstream and downstream side of the bedform.20 Hence, we anticipated, that higher amount of bedforms would lead to higher total exchange flux, i.e. the volume of water exchanged between sediment and surface water per day.21 Consequently, the increasing solute transfer to the hyporheic zone leads to higher contact of solutes to bacterial communities and thus higher potential for micropollutant degradation in general. We aimed at creating three contrasting levels of hyporheic exchange by minimum, medium and maximum number of bedforms feasible within the present setting. Minimum exchange was expected for flat sediment (B0), followed by a medium level (3 bedforms, B3) and a high level (6 bedforms, B6). The amount of sediment was identical in all flumes and the shape of the individual bedforms was the same in the B3 and B6 flumes. The triangular shape was determined by practicality within the setting aiming at uniform shapes across flumes and inducing hyporheic exchange rather than mimicking shapes commonly found in the field. Differences in hyporheic exchange between levels B0, B3 and B6 were investigated through a salt tracer dilution test (performed at the end of the experiment) from which exchange flux, exchange volumes and average residence times were calculated (details in Chapter 2.8).
The response surface model explores non-linear effects of the bacterial diversity (S) and hyporheic exchange (B) on the dissipation half-lives (DT50s) by fitting the responses to a quadratic equation (eqn (1)):
DT50 = β0 + β1S + β2B + β3SB + β4S2 + β5B2 + ε | (1) |
The central composite face design used here is a factorial design consisting of 20 flumes (Fig. 1a): eight flumes with the factorial variable combinations, eight flumes with axial combinations and four replicates of the center-point experiments to validate the response surface model. Central composite designs are commonly used for response surface models because they are easy to expand (e.g. to include more variables) and flexible in terms of choosing the values of each variable at the axial and center-points.
S1 | S3 | S6 | Erpe sediment | Sand | |
---|---|---|---|---|---|
a n.a.: not applicable. | |||||
Sand base | 21.3 kg sand | 23.7 kg sand | 23.7 kg sand | n.a. | n.a. |
Inoculum | 2 L Erpe sediment | Inoculum 1 (2 mL Erpe sediment in 2 L deionised water) | Inoculum 2 (2 mL inoculum 1 in 2 L deionised water) | n.a. | n.a. |
TC [%] | 0.007 ± 0.002 | 0.004 ± 0.003 | 0.003 ± 0.001 | 0.840 ± 0.066 | 0 ± 0 |
Fine gravel [%] | 5 | 5 | 5 | 6 | 5 |
Coarse sand [%] | 6 | 5 | 5 | 13 | 5 |
Medium sand [%] | 82 | 83 | 83 | 68 | 83 |
Fine sand [%] | 6 | 6 | 6 | 12 | 6 |
<0.063 mm [%] | <1 | <1 | <1 | <1 | <1 |
K f at 10 °C [m s−1] | 3.14 × 10−4 ± 4% | 3.37 × 10−4 ± 14% | 3.37 × 10−4 ± 14% | n.a. | n.a. |
Porosity [%] | 35 | 35 | 35 | n.a. | n.a. |
pH in the surface water was measured two times. From day −4 to day 45, the average pH in the flumes rose from 8.1 (±0.1) to 8.5 (±0.3) (Fig. S4†). While there were no significant differences between treatments at day −4 and day 45, the sediment dilution treatment had significantly influenced pH (ANOVA; p < 0.05). Treatment S1 (8.2 ± 0.3) had a significantly lower pH than S3 (8.7 ± 0.1) and S6 (8.6 ± 0.2; Tukey post hoc test, p < 0.05). Dissolved oxygen (Pro 20 DO Instrument, YSI Incorporated, Yellow Springs, OH, USA) in the surface water was measured 4 times during the attenuation phase. Average O2 saturation in all flumes ranged from 101 to 110% between days 28 and 86 (Fig. S5†). Average O2 saturation at days 28, 36, 44 and 86 was significantly influenced by the sediment dilution treatment (ANOVA; p < 0.05). S1 treatments showed significantly lower O2 saturation (103% ± 0.8) than S3 (105% ± 1.7) and S6 (106% ± 1.9; Tukey post hoc test, p < 0.05).
The method for analysis of nutrients (NO3−, NO2−, NH4+, PO43− and total dissolved nitrogen) and dissolved organic carbon (DOC) is detailed in the ESI.† There was little variation in surface water nutrient dynamics between the bedform treatments, but nutrient concentrations were highly impacted by the level of sediment dilution (Fig. S6†). This is why at day 0 (injection of micropollutants), nutrient concentrations differed between sediment dilution levels. Generally the depletion of nitrogen and DOC during pre-incubation (day −12 to day 0) was higher and faster in treatments with lower dilution. Accordingly, after addition of nutrient solution N2 at day 10, removal of NH4+ was especially high in the lowest dilution (S1) (see ESI† for detailed discussion on nutrient dynamics).
Illumina Miseq amplicon sequencing targeting the 16S rRNA gene using the bacteria specific primer pair 341F and 806R was performed by LGC Genomics GmbH (Berlin, Germany) followed by post processing of the raw data as previously published.42 Taxonomy was assigned using the Ribosomal Database Project (RDP) classifier. For comparative diversity index analyses, uneven sequencing depth among the samples was adjusted by rarefying each sample to an even sequencing depth. The sequence data were submitted to NCBI's sequence reads archive (http://www.ncbi.nlm.nih.gov/sra/) under accession no. PRJNA531245.
Fisher-alpha diversity index was calculated at a genus taxonomic level for all the samples and ANOVA was used to compare the diversity between bedform levels (B0, B3, B6) and sediment dilution levels (S1, S3, S6). Calculations were performed in the R software43 using the vegan44 and phyloseq45 packages.
The DT50s for CBZ and ACS were calculated assuming first order kinetics by fitting the timepoint-averaged measured concentrations to an exponential function (eqn (2a) and (2b)). If no dissipation was observed in the initial timepoints, this period was considered as lag-phase and was excluded from the DT50 calculation.
Cx = C0e−kdist | (2a) |
DT50 = −ln(2)/kdis | (2b) |
The goodness of fit of the first-order dissipation assumption was assessed with a one-tailed t-test of the kinetic constants (kdis) to be significantly different from zero. Only DT50s obtained from kinetic constants (kdis) significantly different from zero (p ≤ 0.05), n = 20 for ACS and n = 18 for CBZ, were fitted to a quadratic response surface model (RSM; eqn (1)), using the rsm package (Lenth RV, 2009) in the R software.43
The coded levels (−1, 0 and 1) were used for the sediment dilution (S) and bedform (B) variables (Table 1). The use of the coded variables to fit the RSM is adequate for our specific goal to understand the relative size and effect of the variables, as in the present study we do not aim to predict the DT50s or optimize the attenuation of micropollutants. The model coefficients βx (eqn (1)) were first calculated using ordinary least squares, then tested to be significantly different from zero (two-tailed t-test), and finally an analysis of variance (ANOVA) was used to evaluate the significance of the first (β1 and β2), second order (β4 and β5) and interaction (β3) terms. The adjusted-R2, F-test and model lack-of-fit were calculated to assess goodness of fit and adequacy of the regression.
Although the combination of medium sediment dilution and medium bedform number (S3:
B3) had slightly higher average bacterial diversity than the other medium dilution (S3) samples (S3
:
B0 and S3
:
B6), this difference was not significant in a Tukey post hoc analysis. The overall difference in Fisher's alpha bacterial diversity was less pronounced between medium and high sediment dilution levels (S3 and S6) than between low and medium dilutions (S1 and S3, Fig. 3b). Since Fisher's alpha index gives more weight to the total number of species, and less to the number of individuals in each species,47 the medium and high dilutions (S3 and S6) were similar in taxonomic richness, but the bacterial community composition could still be considerably different. The dispersion of the Fisher's alpha index differed between the sediment dilution levels: the low dilution (S1) had the widest range in values (RSD 18%), followed by the medium dilution (S3; RSD 15%) and finally the high dilution (S6; RSD 13%). Due to the extinction effect, the communities in the low dilution flumes (S1) evolved from the original heterogeneity of river sediment causing a wider spread of diversities, while the communities in the medium and high dilution flumes (S3 and S6) evolved from inocula with similar pre-selected bacterial communities with lower taxonomic richness.
Due to the dilution-to-extinction method, the number of bacterial cells at the start of pre-incubation (day −12) and corresponding number of copies of the 16S rRNA gene, was theoretically lower in the medium and high sediment dilutions (S3 and S6) compared to the low dilution (S1). This was not the case anymore at day 21. Sediment treatments at day 21, obtained an average of 4.7 × 105 copies of the 16S rRNA gene per gram of sediment dry weight after real-time PCR (Fig. S8†) and were not significantly different between the sediment dilution levels (ANOVA; p > 0.05). Further, there was no significant effect of the bedform variable on the copy numbers of 16S rRNA gene (ANOVA; p > 0.05). Thus, substantial bacterial growth in the high dilutions (S3 and S6) was necessary during the pre-incubation phase to reach a similar number of copies to the low dilution (S1) at day 21, suggesting that the pre-incubation was successful. This reveals that the main effect of the sediment dilution was likely caused by differences in bacterial diversity, and not differences in biomass. Still, it should be pointed out that the copy numbers of 16S rRNA genes do not necessarily correlate to bacterial biomass or number of bacteria.48 First, because the number of copies of the 16S rRNA gene in the chromosome is taxon specific (from 1 to more than 10). Second, because the 16S rRNA gene copy number for a specific taxon can differ between growth phases. Our results for bacterial diversity and 16S rRNA gene copies show the status of bacterial communities in the flumes at day 21, which is representative of the period where most micropollutant attenuation occurred, yet it is possible that the communities evolved over the duration of the test and in response to environmental conditions such as sunlight, temperature and addition of nutrients.
The test showed that the sediment dilution-variable significantly influenced the exchange flux (ANOVA; p < 0.05), while the bedform-variable had no significant effect at this late point in time. Neither the bedform, nor the sediment dilution variable influenced the exchange volume or the residence time significantly (Fig. S7†). In contrast to what we found, we expected that an increasing number of bedforms would induce higher hyporheic exchange due to advective pumping20 and this way increase contact of micropollutants to hotspots of favorable turnover conditions. Previous studies have found differences in hyporheic exchange in flumes with flat beds and in flumes with bedforms, containing the same sediment.21
We attribute the lack of significant effect of bedforms on exchange parameters to the fact that the salt-tracer dilution test was conducted at the final point of the experiment (days 78–84), a time by which the flumes had changed from their initial setup conditions (Table 3). Formation of biofilms and algae and settling of fine particulate matter, more pronounced at the end of the test, likely influenced the permeability of the sediment by clogging and affected hyporheic exchange. Bedform heights had dropped by 19% compared to the beginning of the setup due to gradual erosion, potentially reducing exchange flux in B3 and B6 treatments (Table 3). Small ripples had formed in the flumes without bedforms caused by turbulence in the flume curves, which likely increased the exchange flow in the B0 treatments over time. The average flow velocity had decreased by 19% from initial bedform setup (day −3) to day 82 due to decreasing pump performance (Table 3). Differences in flow velocity might also have contributed to the effect of the sediment dilution treatment on the exchange flux at the time of the salt tracer test. By then, S1 treatments had the lowest average flow velocity (6.1 cm s−1), followed by S3 treatments (6.8 cm s−1) and S6 treatments (7.9 cm s−1), although these differences were not significant (ANOVA; p > 0.05). The lower flow velocities in S1 flumes might have been caused by clogging of the pumps due to floating algae, which appeared in the second half of the attenuation phase in S1 flumes (Fig. S10†).
Day of measurement | −3 | 27 | 46 | 82 | Reduction by day 27 | Reduction by day 46 | Reduction by day 82 |
---|---|---|---|---|---|---|---|
a n.m.: not measured. | |||||||
Height of the bedforms [cm] | 8 | 7.8 ± 0.8 | 7.3 ± 0.8 | 6.5 ± 0.9 | 3% | 9% | 19% |
Depth valleys [cm] | 2 | n.m. | 1.94 | 1.84 | 3% | 8% | |
Water level from bottom [cm] | 12 | 11.5 ± 0.4 | 11.7 ± 0.4 | 11.2 ± 0.5 | 4% | 3% | 7% |
Surface water velocity bedforms [cm s−1] | n.m. | 8.2 ± 1.2 | 8.3 ± 1.1 | 6.6 ± 1.8 | −1% | 19% | |
Sediment depth flat [cm] | 3.5 | 3.8 ± 0.3 | 3.7 ± 0.2 | 3.4 ± 0.2 | −9% | −6% | 3% |
Surface water velocity flat [cm s−1] | n.m. | 9.2 ± 1.7 | 10.3 ± 1.4 | 7.4 ± 2.2 | −11% | 20% |
While the salt tracer dilution test showed that hyporheic exchange did not differ between bedform treatments at the final point of the experiment, it can not be ruled out that the bedforms caused differences in hyporheic exchange as anticipated in the beginning of the experiment and most relevant time for degradation (Fig. 2a). Hydrodynamic modelling could be used to calculate bedform-induced hyporheic exchange under initial setup conditions in the different bedform treatments.49,50 However, this was outside the scope of the present study. A follow-up study will cover this endeavor.
ACS and CBZ were ubiquitously found in a field study in River Erpe, the river from which the sediment was collected for the flumes, with surface water concentrations in the μg L−1-range and both compounds were more persistent relative to other micropollutants.6 While DT50s of ACS observed in the River Erpe were in the same order of magnitude (4–30 days) as the DT50s in our flume setup, the in situ DT50s of CBZ were about one order of magnitude lower (4–13 days) compared to in our flumes. Also, Writer et al.52 observed lower DT50s of CBZ of 21.0 ± 4.5 h in a small creek in Colorado and Acuña et al.53 found CBZ DT50s of 4.1 ± 2.4 h in rivers in Spain. On the other hand, a previous flume experiment23 showed infinite DT50s of CBZ. Three in situ studies that tested ACS and CBZ found poor to no degradation of both compounds.5,29,54 The variability in literature values of DT50s for both ACS and CBZ reflect the complex interactions between dissipation processes and environmental conditions that in turn may influence degradation rates.
The RSM of ACS explained roughly 90% of the variance of observed DT50s (adjusted-R2 of 90.2%) and the overall model was significant, meaning that the coefficients together fitted the DT50s better than just the mean (F-statistic: 35.79 on 5 and 14 degrees of freedom, p ≤ 0.05, see Table 4 and Fig. 4c). The ANOVA showed that both first, second order and interaction terms were significant for the RSM. However, the lack of fit of the model was also significant, which means that the difference between the model predictions and the average of measurements for each variable combination (S1:
B0, S3
:
B0, S6
:
B0, etc.) was large compared to the pure error that was expected due to chance. The difference between predicted vs. observed DT50s was higher in S3 and S6 levels, which also have the highest variations in DT50s (Fig. 4b). Therefore, variation in the DT50s was not adequately explained by the variables S and B in the ACS-RSM. The CBZ-RSM was found significant and had an adjusted-R2 of 44.4% (Table 4 and Fig. 4f). Only the linear term of the model was significant for the fit and thus the quadratic terms were not significantly improving the model performance. The overall interpretation of the CBZ-RSM model is that only the sediment dilution variable (S) was a good explanatory variable of the DT50s of CBZ in our flume set-up.
Acesulfame | Carbamazepine | |||
---|---|---|---|---|
Coefficient | p-Value | Coefficient | p-Value | |
a * indicates significance at the 0.05 level. | ||||
RSM parameters | ||||
Intercept – (β0) | 13.95 | 4.30 × 10−10* | 357.72 | 1.78 × 10−4* |
Linear term – S (β1) | −7.55 | 9.05 × 10−9* | −172.78 | 0.002* |
Linear term – B (β2) | −1.99 | 0.007* | −19.18 | 0.671 |
Interaction term – SB (β3) | 2.36 | 0.008* | 3.53 | 0.950 |
Quadratic term – S2 (β4) | −3.70 | 0.002* | −55.58 | 0.440 |
Quadratic term – B2 (β5) | −0.46 | 0.653 | −70.49 | 0.331 |
Adjusted-R2 | 0.902 | 0.444 | ||
Model p-value | 1.71 × 10−7* | 0.029* | ||
![]() |
||||
RSM ANOVA | ||||
First order (β1 and β2) | 2.7 × 10−8* | 0.005* | ||
Second order (β4 and β5) | 0.007* | 0.442 | ||
Interaction (β3) | 0.008* | 0.826 | ||
Lack of fit | 0.048* | 0.364 |
The main removal process of ACS and CBZ in this experiment was expected to be biodegradation. Although CBZ was often found persistent and thus not biodegradable in sediment–water systems,8,13,23,55 in experiments comprising longer incubation times, slow biodegradation has been observed previously.9,56–58 Biodegradation of ACS in WWTPs was only recently reported and has evolved either by independent evolution or global spreading of organisms or genes responsible for the biotransformation pathway.37 Bacteria in sediment of River Erpe have been chronically exposed to ACS via the effluent of the local WWTP and could have developed the capability to degrade ACS similar to the WWTP communities. Neither photolysis nor hydrolysis has been relevant for ACS and CBZ attenuation in previous studies.8,59
Sorption is expected to influence the dissipation from water for compounds with high logKOW as they sorb to hydrophobic fractions of organic matter, or for positively charged compounds due to electrostatic interactions with negatively charged binding sites of clay minerals, biofilm or organic matter.52,60,61 Low sorption potential of ACS is expected, as ACS is negatively charged at ambient pH of around 8 and has a log
KOW of −0.6. Accordingly, previous studies showed low sorption62 and a low retardation coefficient (max. 1.2)30 in sediments. All flume sediments had low organic carbon content (TC < 0.01%) and low content of fine particles (silt + clay < 1%), which generally rendered them poor sorbents. CBZ is neutral, has a log
KOW of 2.77, and has been shown to sorb to organic matter63 and to have relatively high retardation coefficient in sediment (3.6).30 However, a test following OECD guideline 106
64 to examine potential sorption of CBZ to the flume sediments showed no sorption of CBZ to S1, S3 or S6 sediments (see ESI for details, Fig. S11 and S12†). This can be explained mainly by the low carbon content of the flume sediments. Therefore, sorption was not a relevant attenuation process of CBZ or ACS in our experimental setup.
Manipulating the bacterial community without changing chemical and physical background conditions was challenging. The difference in pH (∼0.5) and O2 (∼3%) in the low dilution (S1) compared to the medium and high levels (S3 and S6) was significant, but likely too low to influence the degradation of micropollutants considerably. In contrast, the differences in nutrient concentrations and nutrient depletion over time were high between the different sediment dilution treatments (Fig. S6†). The low sediment dilution level (S1) had the highest depletion of nutrients and lower long-term availability of nutrients for bacteria. Hence, DT50s decreased in the lower sediment dilution although nutrient concentrations were lower. Thus, the deficiencies in nutrient concentrations could have counteracted the positive effect of the high bacterial diversity in S1 which implies that the effect of the sediment dilution on DT50s of ACS and CBZ might have been even higher at equal nutrient availabilities. Therefore, despite some differences in chemical conditions, the significant influence of sediment dilution levels on degradation of ACS and CBZ can be attributed to the differences in bacterial diversity. A positive association between biodegradation rates and biodiversity has been reported for bacteria in WWTPs.33,34 The faster dissipation of ACS and CBZ in sediment with the highest diversity (S1) is consistent with a mechanism referred to as “sampling effect”. Briefly, if the biodegradation of a micropollutant is a rare function performed by few bacteria species, which is likely for slowly degrading compounds, the bacteria capable of degrading these compounds are more likely found in sediment with the highest bacterial diversity. In addition, the “complementarity effect”, meaning that higher diversity leads to higher functionality and more cooperating specialist strains, might increase complexity and efficiency in transformation pathways of particular micropollutants.65 The “sampling” and “complementarity” effects can be applied also to describe the increasing variability of DT50s of ACS in the higher sediment dilutions (S3 and S6): key species necessary for the efficient degradation of ACS are less likely to be present or randomly found in bacteria communities with low diversity. The variation in bacterial communities in the sediment dilutions thus might be responsible for the lack of fit of the RSM of ACS.
The coefficients in the RSM showed that the effect of the bedform variable (B) on the DT50s was smaller than the effect of the sediment dilution (S) and the interaction between the two variables (S and B), was significant and positive for ACS (Table 4, Fig. 4a). Therefore, the combination of high sediment dilution (low bacterial diversity) and absence of bedforms (low hyporheic exchange), had a negative effect on the attenuation of ACS and resulted in longer DT50s (Fig. 4c). The significant interaction between sediment dilution and bedform also implies that the effect of bedform treatment was higher at higher sediment dilution (Fig. 4c), indicating that the strong impact of high diversity on the dissipation rate of ACS in the low dilution treatment might have covered the effect of the bedform variable.
For CBZ, no effect of bedforms on dissipation rates could be shown, because the coefficients in the RSM of the bedform variable (B) and the interaction coefficient were not significant (p > 0.05, Table 4, Fig. 4d). In two of the flumes, one from level B0 and one from B3, it was not possible to calculate DT50s, and they were excluded from the RSM calculations. The absence of these two potentially long DT50s could have additionally masked an effect of the bedform variable. Given the generally slow degradation of CBZ, the effect of the bedform levels might have been too low to be detectable in the present setup.
In the future, the setup we presented can be applied to other parameters that influence degradability, for instance sediment properties, nutrient concentrations or plant cover. For future studies, we suggest a set of preliminary studies or modelling to test the suitability of appropriate treatment levels for all variables tested. For instance, higher effects of and differences in hyporheic exchange can be ensured by preliminary hydrodynamic modeling of bedform morphologies. A second important improvement would be the continuous monitoring of chemical parameters to maintain comparable nutrient and chemical conditions. We used coded variable levels to investigate the suitability of the experimental design, but future studies could enhance the use of this design to quantify the effect of variables on DT50s of micropollutants.
The use of a central composite face design to fit a response surface model provides a robust statistical method to study the response of micropollutant attenuation to certain environmental drivers and highlight bacterial diversity as a disregarded factor in testing persistence of micropollutants. To the best of our knowledge, no other flume study includes a full experimental design with replicates. This study guides a way forward for more elaborate experimental setups that address how environmental processes affect the fate of micropollutants in rivers.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9em00327d |
‡ These authors contributed equally to the publication. |
§ These authors share the last authorship. |
This journal is © The Royal Society of Chemistry 2019 |