Biological Limitations of Dechlorination of cis-Dichloroethene during

Dec 13, 2017 - School of Public Health, Department of Environmental and Occupational Health, Texas A&M University, College Station, Texas 77843, Unite...
0 downloads 7 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Biological limitations of dechlorination of cisdichloroethene during transport in porous media Itza Mendoza-Sanchez, Robin L Autenrieth, Thomas J McDonald, and Jeffrey A. Cunningham Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04426 • Publication Date (Web): 13 Dec 2017 Downloaded from http://pubs.acs.org on December 14, 2017

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

Environmental Science & Technology

Biological limitations of dechlorination of cis-dichloroethene during transport in porous media Itza Mendoza-Sancheza,b,*, Robin L. Autenriethb, Thomas J. McDonalda, and Jeffrey Cunninghamc a

Texas A&M University, School of Public Health, Department of Environmental and

Occupational Health, College Station, Texas 77843, USA b

Texas A&M University, Department of Civil Engineering, College Station, Texas 77843, USA

c

University of South Florida, Department of Civil and Environmental Engineering, Tampa,

Florida 33620, USA



AUTHOR INFORMATION

Corresponding Author *Address: Texas A&M University, School of Public Health, Department of Environmental and Occupational Health, 1266 TAMU College Station, TX 77843-1266; e-mail: [email protected]; telephone: +1-979-436-9329; fax: +1-979-436-9590.

word count in text (excluding References and Figure Captions):

4427

equivalent word count from 1 table, 2 small figures, 3 large figures:

2700

overall word count:

7127

1 ACS Paragon Plus Environment

Environmental Science & Technology

1

ABSTRACT: We applied a mathematical model to data from experimental column studies to

2

understand the dynamics of successful and unsuccessful reductive dechlorination of chlorinated

3

ethenes in groundwater under different flow conditions. In laboratory column experiments

4

(reported previously), it was observed that complete dechlorination of cis-dichloroethene to ethene

5

was sustained at high flow velocity (0.51 m/d), but that dechlorination failed at medium or low flow

6

velocity (0.080 or 0.036 m/d). The mathematical model applied here accounts for transport of

7

chlorinated ethenes in flowing groundwater, mass transfer of chlorinated ethenes between mobile

8

groundwater and stationary biofilms, and diffusion and biodegradation within the biofilms. Monod

9

kinetics with competitive inhibition are used to describe biodegradation. Nearly all parameters

10

needed to solve the model are estimated independently from batch and non-reactive transport

11

experiments. Comparing the model predictions to the experimental results permits the evaluation of

12

three hypothesized biological limitations: insufficient supply of electron donor, decay of

13

dechlorinators’ biomass, and reduction in bacterial metabolism rates. Any of these three limitations

14

are able to adequately describe observed experimental data, but insufficient supply of electron

15

donor is the most plausible explanation for failure of dechlorination. Therefore, an important

16

conclusion of this investigation is that insufficient hydrogen production occurs if groundwater flow

17

is too slow to provide adequate flux of electron donor. Model simulations were in good agreement

18

with experimental results for both successful and unsuccessful dechlorination, suggesting the model

19

is a valid tool for describing transport and reductive dechlorination. An implication of our findings

20

is that in engineered or natural bioremediation of chloroethene-contaminated groundwater, not only

21

must the proper dechlorinating organisms be present, but also proper groundwater flow conditions

22

must be maintained or else dechlorination may fail.

2 ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

Environmental Science & Technology

23

INTRODUCTION

24

Bioremediation has extensively been studied for transforming perchloroethene (PCE),

25

trichloroethene (TCE), and dichloroethene (DCE) to non-toxic ethene in chloroethene-

26

contaminated groundwater. Under anaerobic conditions, such as are present in many

27

groundwater environments, biodegradation of PCE and TCE occurs through sequential reductive

28

dehalogenation in which the parent compounds are reduced to cis-dichloroethene (cis-DCE),

29

followed by vinyl chloride (VC), and finally ethene. The success of bioremediation in a

30

chloroethene-contaminated aquifer relies on the achievement of complete dechlorination to

31

ethene, because VC is the most toxic of the family of chlorinated ethenes.

32

Although complete or near-complete bioremediation of chloroethene-impacted

33

groundwater has been achieved in many instances, it has also been observed that complete

34

dechlorination to ethene is often not achieved.1–5 This sometimes leads to the accumulation of

35

partial dechlorination products, particularly cis-DCE and VC.3,6,7 Failure of complete

36

dechlorination can be caused by several factors. For instance, although a number of bacteria

37

have been found capable of partially reducing PCE or TCE, only some isolates of the species

38

Dehalococcoides mccartyi are able to fully reduce PCE and TCE to ethene.8–10 Therefore, partial

39

or incomplete dehalogenation in groundwater is sometimes attributed to the absence of D.

40

mccartyi in the aquifer.11,12 However, other factors can also lead to incomplete dechlorination.

41

For instance, it has also been found that dechlorination rates may slow down if there is

42

insufficient electron donor (hydrogen) or nutrients in the system.13–16 At any given contaminated

43

site, it is vital to select and design a remediation approach that will overcome any critical factors

44

that might limit the effectiveness of biodegradation. However, to ensure this, it is first necessary

45

to have an appropriate and sophisticated understanding of the different biological factors that

46

could possibly limit dechlorination.

3 ACS Paragon Plus Environment

Environmental Science & Technology

47

Our research group previously reported17 the results of laboratory column experiments in

48

which near-complete dechlorination of cis-DCE was observed at high flow velocity (0.51 m/d),

49

but not at medium (0.080 m/d) or low (0.036 m/d) flow velocity. To explain this observation, we

50

hypothesized that flow rate controls the flux of electron donor through the column. Because the

51

columns were operated at fixed influent concentrations of yeast extract (the electron donor), the

52

flux of electron donor into the columns was directly proportional to the velocity. Hence, the

53

failure of low- and medium-velocity columns to sustain dechlorination could have been because

54

the flux of electron donor was inadequate for that process. However, that hypothesis has not yet

55

been tested, and therefore the cause of the observed velocity-dependent dechlorination is still

56

unknown.

57

If an appropriate mathematical model is available, the column data of Mendoza-Sanchez

58

et al.17 present an opportunity to quantitatively test different hypotheses for why the success of

59

dechlorination in the columns was velocity-dependent. Our research group previously developed

60

mathematical models and algorithms18,19 for simulating groundwater flow with dechlorination,

61

and these models may be appropriate for describing the laboratory column data. However, until

62

now, the model has not yet been applied to the column data, or indeed to any experimental data

63

set, and thus it is not yet known if the model is able to achieve its intended purpose of describing

64

transport with reductive dechlorination. Application of the mathematical model to the

65

experimental data might therefore address two research needs simultaneously:

66

verification/validation of the model, and elucidation of the critical mechanisms at work in the

67

laboratory experiments.

68

Therefore, the objectives of this paper are: (1) to apply the mathematical model of

69

Mendoza-Sanchez and Cunningham18 to the experimental data of Mendoza-Sanchez et al.17;

70

(2) to estimate all model parameters a priori via literature values or results of independent

71

experiments; (3) to validate or verify the capabilities of the model via comparison of the model 4 ACS Paragon Plus Environment

Page 4 of 29

Page 5 of 29

Environmental Science & Technology

72

predictions to the experimental results; (4) to identify, via the comparison of model predictions

73

to experimental results, what biological mechanisms are capable of explaining the observed

74

effect of flow velocity on dechlorination; and (5) to assess which of the candidate mechanisms

75

are most plausible in the context of the experiments under consideration. Through these

76

objectives, we aim to achieve an overall goal of improving our understanding of what factors can

77

limit complete reductive dechlorination in chloroethene-contaminated groundwater.

78 79 80

MATERIALS AND METHODS Column and Batch Experiments. Details of the column experiments have been

81

presented previously17 and a schematic representation of the experimental set-up is shown in

82

Figure 1. Briefly, 60-cm glass columns were filled with homogeneous glass beads and

83

inoculated with the commercially available KB-1® culture (SiREM, Guelph, Ontario) which is

84

known to be capable of complete dechlorination of PCE to ethene.20 Synthetic anaerobic

85

groundwater containing a known and constant concentration of cis-DCE (30 µM, or 2.9 mg/L)

86

was pumped through the columns, as shown in Figure 1. The anaerobic medium contained yeast

87

extract to act as electron donor to support production of hydrogen for dechlorination. Three

88

different pore velocities, termed slow (0.036 m/d), medium (0.080 m/d), and fast (0.51 m/d),

89

were tested in duplicate. These velocities correspond to residence times of 16.7 d, 7.5 d, and 1.2

90

d, respectively, for one pore volume of fluid. The duration of the experiments were 41 d, 45 d,

91

and 47 d, respectively, corresponding to 2.5, 6, and 40 pore volumes treated for the slow-,

92

medium-, and fast-flow columns. Concentrations of cis-DCE, VC, and ethene were monitored

93

over time and space by sampling at 5 different sampling ports, spaced every 10 cm along the

5 ACS Paragon Plus Environment

Environmental Science & Technology

94

column (see Figure 1), and analyzed by gas chromatography with flame ionization detection

95

(GC/FID).

96

In addition to these column experiments, we also conducted three non-reactive abiotic

97

control experiments with columns that were not inoculated with KB-1®, and two reactive

98

controls with columns inoculated with KB-1® but fed with solution that did not contain any cis-

99

DCE. These control experiments verified that disappearance of cis-DCE and production of VC

100

and ethene are biologically mediated.

101

Prior to conducting the column experiments, batch experiments were conducted inside an

102

anaerobic chamber (10% H2, 10% CO2, 80% N2). Batch microcosms, each containing 60 mL of

103

anaerobic medium, were inoculated with KB-1® culture and spiked with cis-DCE to obtain an

104

initial concentration of 30 µM. Abiotic controls were also conducted with batch bottles that

105

contained cis-DCE but no KB-1®, and bioreactive controls were conducted with batch bottles

106

that contained KB-1® but no cis-DCE. Periodically, samples were taken from the batch reactors

107

and analyzed via GC/FID for concentrations of cis-DCE, VC, and ethene. These batch

108

experiments were used to verify complete biologically mediated dechlorination of cis-DCE, and

109

to estimate biodegradation kinetic rate parameters. Results from the batch experiments are

110

presented in the Supporting Information (SI) and have not been presented previously.

111

Mathematical Model for Transport and Reaction. In a previous paper18 we

112

presented a mathematical model that accounts for three important processes: transport of

113

dissolved cis-DCE in the mobile fluid by advection and dispersion, mass transfer of chemical

114

species between the mobile fluid and stationary biofilms, and diffusion and non-linear Monod

115

kinetics (for both bacterial growth and biodegradation) within the biofilms. The Monod-kinetic

116

model, adopted here, was presented by Cupples et al.21 and describes reductive dehalogenation

6 ACS Paragon Plus Environment

Page 6 of 29

Page 7 of 29

Environmental Science & Technology

117

from cis-DCE to ethene under substrate-limiting conditions and competitive inhibition between

118

cis-DCE and VC. The equations describing these processes are coupled, non-linear, partial

119

differential equations, for which efficient solution algorithms were also developed.18,19 A non-

120

dimensional version of the system of equations, with definitions of all variables, is presented in

121

the SI.

122

Estimation of Model Parameters. An objective of this work is to identify biological

123

mechanisms capable of explaining the observed effect of flow velocity on dechlorination by

124

comparing model predictions to experimental results. Application of the mathematical model

125

requires estimation of a number of parameters (e.g., flow velocity, dispersion coefficients,

126

Monod kinetic parameters). In this work, all model parameters for the column experiments

127

were estimated independently from values in literature or from separate batch and non-reactive

128

transport experiments. Therefore, estimation of parameters via “fitting” the model to

129

experimental column results was not required. In the mathematical model, there are 32 required

130

parameters for conducting the simulations. Of this total, 6 parameters were directly measured, 2

131

were estimated from batch experiments, 4 were estimated from non-reactive transport

132

experiments, 10 were estimated from the results of Cupples et al.,21 and 10 were estimated from

133

other papers in the literature. Details are provided in the following paragraphs and in the SI.

134

Biological Monod-kinetic parameters were defined by three means: (a) estimated or

135

directly measured from the experimental set-up; (b) obtained from values reported by Cupples et

136

al.;21 and/or (c) obtained from values that best reproduced separate batch reactive experimental

137

results. The last method was used to estimate only two parameters: an initial hydrogen

138

concentration, and a hydrogen production rate (i.e., rate at which H2 is produced from the

139

electron donor). To the best of our knowledge, experimental or reference values are not

7 ACS Paragon Plus Environment

Environmental Science & Technology

140

available for these two parameters in the literature. To verify that the set of estimated Monod-

141

kinetic parameters from Cupples et al.21 are appropriate for our experimental setting, we applied

142

our mathematical biodegradation model to the batch experiments and verified that model

143

simulations matched experimental results.

144

Following estimation of the biological kinetic parameters, hydrodynamic dispersion

145

coefficients for the transport model were obtained by fitting advective-dispersive simulations to

146

experimental results of non-reactive column controls. Finally, other model parameters such as

147

mass-transfer coefficients, molecular diffusion coefficients, biofilm diffusion coefficients, and

148

biofilm thickness were estimated from values or correlations available in the literature.

149

Additional details are provided in the SI.

150

Application of the Model to Experimental Data. As discussed elsewhere,17

151

dechlorination of cis-DCE to ethene in the column experiments was successful under fast flow

152

conditions, but not under medium or slow flow conditions. To explain this behavior, we now

153

apply the mathematical model to all three sets of column data. For the fast-flow experiments, the

154

model is applied with all model parameters estimated independently as described above, i.e.,

155

there are no adjustable or “fitting” parameters employed. We refer to these model simulations as

156

“baseline” simulations.

157

For the medium- and slow-flow experiments, some modification to the model is

158

necessary to explain the failure of the dechlorination. As will be seen subsequently,

159

dechlorination began to fail after about 10 d (0.6 pore volumes) for the slow-flow columns, and

160

after about 15 d (2 pore volumes) for the medium-flow columns. A possible explanation for this

161

observation is that the electron donor originally present in the columns was either consumed or

162

flushed out of the columns within 10–15 d, and thereafter, the hydrogen production rate

8 ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

Environmental Science & Technology

163

decreased because hydrogen-producing bacteria lacked sufficient electron donor. Therefore, for

164

the medium- and slow-flow columns, we modify the mathematical model to reduce the hydrogen

165

production by a factor of 10 after some critical time (2 pore volumes for medium-flow columns

166

and 0.6 pore volumes for slow-flow columns). In these “hydrogen-limited” simulations, no other

167

model parameters are changed; the only change is the reduction of hydrogen production after the

168

observed critical time, presumably due to insufficient supply of electron donor.

169

Then, to explore if other possible factors could adequately describe the failure of

170

dechlorination, we also run two other sets of simulations in the medium- and slow-flow columns.

171

In “decay-limited” simulations, we increase the bacterial death rate coefficient by a factor of 10

172

after the observed critical time. In “metabolism-limited” simulations, we decrease the cis-DCE

173

biodegradation rate coefficients (Monod kinetic parameters) by a factor of 10 after the observed

174

critical time. Table 1 summarizes the scenarios considered and the associated model parameter

175

modifications for the unsuccessful dechlorination experiments. Comparison of all model

176

simulations to experimental data enables an evaluation of whether any of these possible

177

explanations can be rejected as invalid.

178 179

RESULTS AND DISCUSSION

180

The simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for three

181

sampling ports (ports 2, 3, and 4, as indicated in Figure 1) of the fast-, medium-, and slow-flow

182

column experiments are depicted in Figures 2, 3, and 4, respectively. Data points represent

183

arithmetic averages of concentrations measured in two duplicate columns, and error bars indicate

184

the two actual concentrations measured in the individual columns. As seen from Figure 1, port 2

185

is the closest to the column inlet of these three ports, and therefore data at port 2 (bottom-row

9 ACS Paragon Plus Environment

Environmental Science & Technology

186

panels in the Figures 2, 3 and 4) represent earlier travel time than the data from port 3 (middle-

187

row panels) or port 4 (top-row panels).

188

Figure 2, 3, and 4 show that dechlorination of cis-DCE to ethene was successful under

189

fast-flow conditions, but not under medium- or slow-flow conditions. To translate these results

190

to in situ conditions, the slow (3.6 cm/d) and medium (8.0 cm/d) flow rates are both in the range

191

of natural gradient conditions, and the fast flow rate (51 cm/d) might be observed under forced-

192

gradient conditions, e.g., if pumping is occurring during an active site remediation.

193

Modeling Successful Dechlorination. As seen in Figure 2, cis-DCE was efficiently

194

converted to ethene in the fast-flow column experiments. Although cis-DCE initially broke

195

through at close to its injected concentration of 30 µM (Figure 2, left-column panels), it was

196

subsequently biodegraded, accompanied by a transitory rise in the concentration of the reactive

197

intermediate VC (Figure 2, middle-column panels), and eventually stoichiometric formation of

198

ethene (Figure 2, right-column panels). To simulate experimental observations of the fast-flow

199

columns, we have applied the mathematical model as described above. All model parameters

200

have been estimated independently, not “fit” to the experimental data, and parameter values are

201

given in Table 1. As seen from Figure 2, the model successfully describes all important behavior

202

of cis-DCE, VC, and ethene, suggesting that the mathematical model developed previously18 is

203

properly accounting for the most important physical and biological processes occurring in the

204

system. In Figure 2, we have labeled the model simulation as the “baseline simulation” because

205

it is the model simulation obtained with all model parameters estimated a priori and with no

206

modifications made to the model.

207

Modeling Unsuccessful Dechlorination. As described above, successful

208

dechlorination was not sustained in the slow- and medium-flow column experiments, so we

10 ACS Paragon Plus Environment

Page 10 of 29

Page 11 of 29

Environmental Science & Technology

209

modified the “baseline” model to account for different factors that could possibly be limiting

210

dechlorination (Table 1).

211

Figure 3 shows a comparison of the “baseline” and “hydrogen-limited” model

212

simulations to the experimental data for the slow-flow columns. The hypothesis of the

213

hydrogen-limited model simulation is that production of hydrogen decreased substantially after

214

the initial electron donor was consumed or flushed out of the column. Examination of Figure 3

215

shows that both models successfully describe the observations initially: concentrations of cis-

216

DCE and VC are low, and observed concentrations of ethene are high, indicating successful

217

complete dechlorination at early times. However, following this initial successful period, the

218

observed cis-DCE concentration rises towards its influent concentration, with some

219

transformation to VC but essentially no formation of ethene. The baseline model does not

220

describe this failure of dechlorination, which begins after about 0.6 pore volumes. However, the

221

modified (hydrogen-limited) model is able to capture all important features of the experimental

222

results. This suggests that the proposed hypothesis, i.e., decreased hydrogen production after

223

some critical time (in this case 0.6 pore volumes), is a viable hypothesis to explain unsuccessful

224

dechlorination in the slow-flow columns.

225

However, the success of the hydrogen-limited model does not rule out the possibility that

226

some other mechanism could possibly be responsible for the failure of sustained dechlorination

227

in the slow- and medium-flow columns. In particular, we have offered two alternate hypotheses,

228

the “decay-limited” and “metabolism-limited” scenarios, as described above and summarized in

229

Table 1. Figure 4 shows a comparison of the hydrogen-limited, decay-limited, and metabolism-

230

limited model simulations to the experimental data for the medium-flow columns. Examination

231

of Figure 4 shows that all three modified models are able to qualitatively describe the observed

11 ACS Paragon Plus Environment

Environmental Science & Technology

232

behavior of the system. As with the slow-flow column (Figure 3), the dechlorination was

233

initially successful, with cis-DCE converted to ethene, but after some time, dechlorination failed,

234

and cis-DCE broke through the column at close to its influent concentration, with partial

235

conversion to VC, and only minimal conversion to ethene. The hydrogen-limited simulation and

236

the decay-limited simulation give nearly identical results, so that neither hypothesis can be

237

considered “better” based on agreement with the experimental data. The metabolism-limited

238

simulation gives somewhat different results, especially at late times, but still cannot be ruled out

239

based on agreement (or lack thereof) between simulation predictions and experimental data.

240

It is also worth noting that none of the three modified simulations have been “optimized”

241

to fit the experimental data. In each of the three modified simulations, one of the relevant rate

242

parameters was increased or decreased by a factor of 10 after some critical time, as summarized

243

in Table 1. It is possible that altering the rate by a different factor (other than 10) might lead to

244

improved agreement between model simulation and experimental results. However, throughout

245

this analysis, we have estimated all model parameters independently of the experimental results.

246

Therefore, since none of the three models have been “optimized,” we do not believe it is possible

247

to conclude that any model is better or worse than the others based on agreement with the

248

experimental data. The salient point is that all three models are able to reasonably describe the

249

observed experimental results, so to conclude which of the hypothesized rate limitations is most

250

reasonable, additional analysis is required, as described subsequently.

251

Assessment of the Mathematical Model. The mathematical model for transport and

252

biodegradation previously presented by Mendoza-Sanchez and Cunningham18 had not heretofore

253

been tested via comparison to experimental data; thus, one of the objectives of this paper was to

254

assess the model. The “baseline” model simulations of Figure 2 were performed with all model

12 ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29

Environmental Science & Technology

255

parameters estimated independently of the column experiments, i.e., the model was not “fit” to

256

the data. Despite this, agreement between the model and the data (Figure 2) is excellent, i.e., the

257

model is able to describe all key behaviors of the cis-DCE, VC, and ethene. We conclude that

258

the proposed model is “valid” for describing one-dimensional transport and reaction of

259

chlorinated ethenes under the proper conditions. It may therefore be useful for such applications

260

as predicting the fractional conversion of cis-DCE to VC and ethene at a contaminated site,

261

predicting breakthrough curves or concentration profiles at a contaminated site, or determining

262

how much travel time is required to achieve a specified remediation goal. The model could also

263

be useful for identifying particular causes of biological limitations in other laboratory or field

264

experiments where dechlorination was observed not to proceed to completion.

265

However, we here note five important areas where the model could potentially be

266

improved in the future. (1) The model accounts for cis-DCE, VC, and ethene, but not the

267

commonly encountered perchloroethene or trichloroethene; the model could be extended to

268

include these parent compounds. (2) The model could be extended to account for different D.

269

mccartyi subpopulations, which has been recommended for environments under electron donor

270

limiting conditions.22 (3) The model is a “biofilm” model, i.e., it accounts for mass transfer

271

between the mobile fluid and stationary biofilms, and it assumes that the biotransformation

272

reactions occur within the biofilms. Cunningham and Mendoza-Sanchez have previously

273

considered23 the conditions under which such a “biofilm” model can be replaced by a simpler

274

model that treats biotransformation reactions as occurring in the aqueous (mobile) phase rather

275

than within stationary biofilms. Modifying the model in this fashion would presumably enable

276

the model to run much faster, and under certain conditions23 would not result in any loss in

277

accuracy. (4) The model is currently capable only of describing one-dimensional transport, such

13 ACS Paragon Plus Environment

Environmental Science & Technology

278

as transport within laboratory columns or sufficiently homogeneous field sites. Most field sites

279

would require the model to be modified to include heterogeneous and multi-dimensional flow

280

fields. (5) As discussed previously, the “baseline” model is appropriate only when sustained

281

biodegradation is successful, and when dechlorination fails, the baseline parameters of the model

282

must be modified. Although the model results with the modified parameters were observed to be

283

successful (Figures 3 and 4), the model does not have a built-in method for a priori estimation of

284

when dechlorination will succeed or fail, nor how to modify the model parameters properly for

285

conditions of unsuccessful dechlorination. Developing these methods would likely require

286

additional experimental work as well as modifications to the model.

287

In particular, modeling efforts could be improved by measuring concentrations of

288

hydrogen gas and dechlorinating bacteria over space and time. This would help in a variety of

289

ways. First, it would provide additional measurements to which model predictions could be

290

compared, thereby further testing or verifying if the model is properly representing actual

291

conditions. Second, quantification of the reactive biomass in the system would be useful to

292

distinguish among different possible biological limitations, i.e., hydrogen, metabolism, or decay

293

limitation. Third, the model currently uses a very rough estimate for the hydrogen production

294

rate, and refinement of that estimate would be useful for ensuring the accuracy of the model.

295

Assessment of Proposed Mechanisms for Unsuccessful Dechlorination. As seen

296

from Figure 4, experimental data for unsuccessful dechlorination can be adequately described

297

with any of three proposed rate-limited model simulations (Table 1). However, we conclude that

298

the most plausible hypothesis is that hydrogen production is limited by insufficient supply of

299

electron donor, as originally proposed by Mendoza-Sanchez et al.17 This conclusion is based on

300

two factors.

14 ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29

301

Environmental Science & Technology

First, in addition to concentration histories such as shown in Figures 2–4, we can examine

302

concentration profiles, i.e., graphs of the measured (or predicted) concentrations as functions of

303

position within the columns. Figure 5 shows simulated and experimental concentration profiles

304

for cis-DCE, VC, and ethene in the slow-flow columns after 1.4 pore volumes and after 2.2 pore

305

volumes. Agreement between the hydrogen-limited simulation and the experimental data is

306

excellent for all three chemicals. As discussed above, none of the modified model simulations

307

have been optimized, and therefore we can not truly compare model performance among the

308

three candidate models; nevertheless, the excellent performance of the hydrogen-limited

309

simulation lends confidence that the model is properly capturing the chemistry and biology of the

310

experiments.

311

Second, the hydrogen-limited simulation appears to have a firmer mechanistic foundation

312

than the other two proposed rate-limited scenarios. It seems plausible that, under low-flow

313

conditions, the original supply of electron donor may have been depleted within a certain period

314

of time, after which hydrogen production decreased due to insufficient supply of electron donor.

315

In contrast, it seems unlikely that the bacterial population(s) would begin to die at a faster rate

316

after a certain period of time, or that the bacterial population’s metabolic processes would slow

317

after a certain period of time. In the absence of a clear rationale for how or why these could

318

occur, we adopt as our working hypothesis that insufficient groundwater flow can limit the flux

319

of electron donor and thereby limit the production of hydrogen needed for sustained

320

dechlorination.

321

Implications for Groundwater Remediation. Based on the preceding analysis, we

322

conclude that reductive dechlorination fails if the supply of electron donor is insufficient to

323

sustain production of hydrogen. Notably, however, in the column experiments of Mendoza-

15 ACS Paragon Plus Environment

Environmental Science & Technology

324

Sanchez et al.,17 the provided concentration of electron donor (yeast extract) was identical in

325

fast-, medium-, and slow-flow columns. The supply of electron donor at a contaminated field site

326

is not controlled only by, say, the concentration in a solution that is injected or amended at the

327

site; it is also controlled by the velocity of the groundwater, which is likely to be distributed

328

heterogeneously at the site. This implies that dechlorination activity may also be distributed

329

heterogeneously, even if the necessary organisms are present throughout the site. Clean-up of a

330

contaminated site must therefore include consideration of the effect of heterogeneous velocity

331

distribution on the flux of delivered chemicals.

332 333



334

Supporting Information

335

Mathematical model description; details on the estimation of transport and biodegradation

336

kinetic parameters; definition of non-dimensional numbers used in the mathematical model;

337

definition of parameters in the non-dimensional numbers; values of Monod kinetic parameters

338

for the baseline case that best describe microbial dechlorination behavior; values of transport

339

and kinetic parameters for the baseline simulation; values of non-dimensional numbers in the

340

three different simulations; simulated and experimental biodegradation results for different

341

cases of initial hydrogen concentration; simulated and experimental biodegradation results for

342

different cases of hydrogen production rates; simulated and experimental breakthrough curves

343

of cis-DCE for non-reactive control columns compared to advective-dispersive simulations.

344

This material is available free of charge via the Internet (PDF).

ASSOCIATED CONTENT

345

16 ACS Paragon Plus Environment

Page 16 of 29

Page 17 of 29

Environmental Science & Technology

346



347

Corresponding Author

348

*I. Mendoza-Sanchez. Address: Texas A&M University, School of Public Health, Department of

349

Environmental and Occupational Health, 1266 TAMU, College Station, TX 77843-1266. E-mail:

350

[email protected]. Tel.: +1-979-436-9329. Fax: +1-979-436-9590.

351

Notes

352

The authors declare no competing financial interest.

AUTHOR INFORMATION

353 354



355

Financial support for this research was provided to Itza Mendoza-Sanchez from the following

356

sources: Consejo Nacional de Ciencia y Tecnología (CONACYT), México; and the United

357

States Geological Survey (USGS) through the Texas Water Resources Institute (TWRI). Any

358

opinions, findings, conclusions, or recommendations expressed in this paper are those of the

359

authors and do not necessarily reflect the views of CONACYT, USGS, or TWRI.

ACKNOWLEDGEMENTS

17 ACS Paragon Plus Environment

Environmental Science & Technology

360



361

(1)

REFERENCES Hendrickson, E. R.; Payne, J. A.; Young, R. M.; Starr, M. G.; Perry, M. P.; Fahnestock,

362

S.; Ellis, D. E.; Ebersole, R. C. Molecular analysis of Dehalococcoides 16s ribosomal

363

DNA from chloroethene-contaminated sites throughout North America and Europe. Appl.

364

Environ. Microbiol. 2002, 68 (2), 485–495.

365

(2)

Liang, Y.; Cook, L. J.; Mattes, T. E. Temporal abundance and activity trends of vinyl

366

chloride (VC) degrading bacteria in a dilute VC plume at Naval Air Station Oceana.

367

Environ. Sci. Pollut. Res. 2017, 24, 13760–13774.

368

(3)

Atashgahi, S.; Lu, Y.; Zheng, Y.; Saccenti, E.; Suarez-Diez, M.; Ramiro-Garcia, J.;

369

Eisenmann, H.; Elsner, M.; Stams, A. J. M.; Springael, D.; et al. Geochemical and

370

microbial community determinants of reductive dechlorination at a site biostimulated with

371

glycerol. Environ. Microbiol. 2017, 19 (3), 968–981.

372

(4)

Adetutu, E. M.; Gundry, T. D.; Patil, S. S.; Golneshin, A.; Adigun, J.; Bhaskarla, V.;

373

Aleer, S.; Shahsavari, E.; Ross, E.; Ball, A. S. Exploiting the intrinsic microbial

374

degradative potential for field-based in situ dechlorination of trichloroethene contaminated

375

groundwater. J. Hazard. Mater. 2015, 300, 48–57.

376

(5)

Mundle, S. O. C.; Johnson, T.; Lacrampe-Couloume, G.; Perez-de-Mora, A.; Duhamel,

377

M.; Edwards, E. A.; Mcmaster, M. L.; Cox, E.; Revesz, K.; Lollar, B. S. Monitoring

378

biodegradation of ethene and bioremediation of chlorinated ethenes at a contaminated site

379

using compound-specific isotope analysis (CSIA). Environ. Sci. Technol 2012, 46, 1731–

380

1738.

18 ACS Paragon Plus Environment

Page 18 of 29

Page 19 of 29

381

Environmental Science & Technology

(6)

Mattes, T. E.; Jin, Y. O.; Livermore, J.; Pearl, M.; Liu, X. Abundance and activity of vinyl

382

chloride (VC)-oxidizing bacteria in a dilute groundwater VC plume biostimulated with

383

oxygen and ethene. Environ. Biotechnol. 2015, 99, 9267–9276.

384

(7)

Zhang, S.; Hou, Z.; Du, X.; Li, D.; Lu, X. Assessment of biostimulation and

385

bioaugmentation for removing chlorinated volatile organic compounds from groundwater

386

at a former manufacture plant. Biodegradation 2016, 27, 223–236.

387

(8)

Maymó-Gatell, X.; Chien, Y.; Gossett, J. M.; Zinder, S. H. Isolation of a bacterium that

388

reductively dechlorinates tetrachloroethene to ethene. Science. 1997, 276 (5318), 1568–

389

1571.

390

(9)

He, J.; Ritalahti, K. M.; Aiello, M. R.; Löffler, F. E. Complete detoxification of vinyl

391

chloride by an anaerobic enrichment culture and identification of the reductively

392

dechlorinating population as a Dehalococcoides species. Appl. Environ. Microbiol. 2003,

393

69 (2), 996–1003.

394

(10)

Löffler, F. E.; Yan, J.; Ritalahti, K. M.; Adrian, L.; Edwards, E.A., Konstantinidis, K.;

395

Müller, J.A.; Fullerton, H.; Zinder, S.H.; Spormann,A.M. Dehalococcoides mccartyi gen.

396

nov., sp. nov., obligately organohalide-respiring anaerobic bacteria relevant to halogen

397

cycling and bioremediation, belong to a novel bacterial class, Dehalococcoidia classis

398

nov., order Dehalococcoidales ord. nov. and family Dehalococcoidaceae fam. nov.,

399

within the phylum Chloroflexi. Int. J. Syst. Evol. Microbiol. 2013, 63, 625-635.

400

(11)

Schaefer, C. E.; Condee, C. W.; Vainberg, S.; Steffan, R. J. Bioaugmentation for

401

chlorinated ethenes using Dehalococcoides sp.: Comparison between batch and column

402

experiments. Chemosphere 2009, 75 (2), 141–148.

403

(12)

Jugder, B.-E.; Ertan, H.; Bohl, S.; Lee, M.; Marquis, C. P.; Manefield, M. Organohalide

19 ACS Paragon Plus Environment

Environmental Science & Technology

404

respiring bacteria and reductive dehalogenases: key tools in organohalide bioremediation.

405

Front. Microbiol. 2016, 7, 1–12.

406

(13) Azizian, M. F.; Behrens, S.; Sabalowsky, A.; Dolan, M. E.; Spormann, A. M.; Semprini,

407

L. Continuous-flow column study of reductive dehalogenation of PCE upon

408

bioaugmentation with the evanite enrichment culture. J. Contam. Hydrol. 2008, 100 (1–2),

409

11–21.

410

(14)

Berggren, D. R. V; Marshall, I. P. G.; Azizian, M. F.; Spormann, A. M.; Semprini, L.

411

Effects of sulfate reduction on the bacterial community and kinetic parameters of a

412

dechlorinating culture under chemostat growth conditions. Environ. Sci. Technol. 2013, 47

413

(4), 1879–1886.

414

(15)

Azizian, M. F.; Marshall, I. P. G.; Behrens, S.; Spormann, A. M.; Semprini, L.

415

Comparison of lactate, formate, and propionate as hydrogen donors for the reductive

416

dehalogenation of trichloroethene in a continuous-flow column. J. Contam. Hydrol. 2010,

417

113 (1–4), 77–92.

418

(16)

Miura, T.; Yamazoe, A.; Masako, I.; Ohji, S.; Hosoyama, A.; Takahata, Y.; Fujita, N. The

419

impact of injections of different nutrients on the bacterial community and its

420

dechlorination activity in chloroethene-contaminated groundwater. Microbes Env. 2015,

421

30 (2), 164–171.

422

(17)

Mendoza-Sanchez, I.; Autenrieth, R.; McDonald, T.; Cunningham, J. Effect of pore

423

velocity on biodegradation of cis-dichloroethene (DCE) in column experiments.

424

Biodegradation 2010, 21 (3), 365–377.

425 426

(18)

Mendoza-Sanchez, I.; Cunningham, J. Efficient algorithms for modeling the transport and biodegradation of chlorinated ethenes in groundwater. Transp. Porous Media 2012, 92

20 ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29

Environmental Science & Technology

427 428

(1), 165–185. (19) Mendoza-Sanchez, I.; Cunningham, J. Efficient algorithm for modeling transport in

429

porous media with mass exchange between mobile fluid and reactive stationary media.

430

Transp. Porous Media 2007, 68 (3), 285–300.

431

(20)

cultures dominated by Dehalococcoides. Fems Microbiol. Ecol. 2006, 58, 538–549.

432 433

Duhamel, M.; Edwards, E. A. Microbial composition of chlorinated ethene-degrading

(21)

Cupples, A. M.; Spormann, A. M.; McCarty, P. L. Vinyl chloride and cis-dichloroethene

434

dechlorination kinetics and microorganism growth under substrate limiting conditions.

435

Environ. Sci. Technol. 2004, 51, 1635–1642.

436

(22)

Mayer-Blackwell, K.; Azizian, M.F.; Green, J.K.; Spormann, A.M.; Semprini, L. Survival

437

of vinyl chloride respiring Dehalococcoides mccartyi under long-term electron donor

438

limitation. Environ. Sci. Technol. 2017, 38 (4), 1102–1107.

439 440

(23)

Cunningham, J. A.; Mendoza-Sanchez, I. Equivalence of two models for biodegradation during contaminant transport in groundwater. Water Resour. Res. 2006, 42 (2), W02416.

441 442

21 ACS Paragon Plus Environment

Environmental Science & Technology

443

Table 1. Summary of modeling scenarios applied to experimental data Scenario baseline simulation

Description of scenario all model parameters estimated from literature values, batch experiments, and non-reactive control experiments

hydrogen-limited simulation

rate of hydrogen production from electron donor reduced by a factor of 10 as compared to baseline simulation

decay-limited simulation

bacterial decay rate coefficient increased by a factor of 10 as compared to baseline simulation

metabolism-limited simulation

maximum bacterial growth rate coefficient and maximum specific substrate utilization rate coefficient both decreased by a factor of 10 as compared to baseline simulation

444 445 446

22 ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29

447

Environmental Science & Technology



FIGURE CAPTIONS

448 449

Figure 1. Schematic representation of the experimental set-up (after Mendoza-Sanchez et

450

al.17). The arrow inside the column indicates the direction of the flow.

451 452

Figure 2. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

453

three sampling ports (port 2, 3, and 4) of the fast-flow column experiments.

454 455

Figure 3. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

456

three sampling ports (port 2, 3, and 4) of the slow-flow column experiments.

457 458

Figure 4. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

459

three sampling ports (port 2, 3, and 4) of the medium-flow column experiments.

460 461

Figure 5. Simulated and experimental concentration profiles for cis-DCE, VC, and ethene in

462

the slow-flow columns after 1.4 and 2.2 pore volumes.

463 464 465

23 ACS Paragon Plus Environment

Environmental Science & Technology

466 467 468 469 470 471 472

Figure 1. Schematic representation of the experimental set-up (after Mendoza-Sanchez et al.17). The arrow inside the column indicates the direction of the flow.

24 ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29

Environmental Science & Technology

473 474 475 476

Figure 2. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

477

three sampling ports (port 2, 3, and 4) of the fast-flow column experiments.

478

25 ACS Paragon Plus Environment

Environmental Science & Technology

479 480 481 482

Figure 3. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

483

three sampling ports (port 2, 3, and 4) of the slow-flow column experiments.

484

26 ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

Environmental Science & Technology

485

486 487 488 489

Figure 4. Simulated and experimental breakthrough curves of cis-DCE, VC, and ethene for

490

three sampling ports (port 2, 3, and 4) of the medium-flow column experiments.

491

27 ACS Paragon Plus Environment

Environmental Science & Technology

492 493 494 495

Figure 5. Simulated and experimental concentration profiles for cis-DCE, VC, and ethene in

496

the slow-flow columns after 1.41 and 2.2 pore volumes.

497

28 ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29

Environmental Science & Technology

TOC art 84x47mm (300 x 300 DPI)

ACS Paragon Plus Environment