Geophysical Monitoring of Hydrocarbon-Contaminated Soils

Jul 7, 2016 - Colorado School of Mines, Department of Geophysics, Golden, ... The removal of hydrocarbon was enhanced by a bioelectrochemical system (...
0 downloads 0 Views 1000KB Size
Subscriber access provided by University of Sussex Library

Article

Geophysical Monitoring of Hydrocarbon-Contaminated Soils Remediated with a Bioelectrochemical System Deqiang Mao, Lu Lu, Andre Revil, Yi Zuo, John Hinton, and Zhiyong Jason Ren Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b00535 • Publication Date (Web): 07 Jul 2016 Downloaded from http://pubs.acs.org on July 8, 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 31

1 2

Environmental Science & Technology

Geophysical Monitoring of Hydrocarbon-Contaminated Soils Remediated with a Bioelectrochemical System

3 4

Deqiang Mao1, Lu Lu2, André Revil3*, Yi Zuo4, John Hinton1, and Zhiyong Jason Ren2

5 6

1

Colorado School of Mines, Department of Geophysics, Golden, 80401, CO, USA

7

2

Department of Civil, Environmental, and Architectural Engineering, University of Colorado Boulder, Boulder,

8

Colorado 80309, United States

9

3

ISTerre, CNRS, UMR CNRS 5275, Université de Savoie, 73376 cedex, Le Bourget du Lac, France

10

4

Chevron Energy Technology Company, San Ramon, California 94583, United States

11 12 13 14

* Corresponding Author: André Revil

15

Phone: 33+ (0)479758715

16

Emails: [email protected]

17 18 19 20 21 22 23 24 25

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 31

26

Abstract. Efficient non-invasive techniques are desired for monitoring the remediation process

27

of contaminated soils. We applied the direct current resistivity technique to image conductivity

28

changes in sandbox experiments where two sandy and clayey soils were initially contaminated

29

with diesel hydrocarbon. The experiments were conducted over a 230 day period. The removal

30

of hydrocarbon was enhanced by a bioelectrochemical system (BES) and the electrical potentials

31

of the BES reactors were also monitored during the course of the experiment. We found that the

32

variation in electrical conductivity shown in the tomograms correlate well with diesel removal

33

from the sandy soil, but this is not the case with the clayey soil. The clayey soil is characterized

34

by a larger specific surface area and therefore a larger surface conductivity. In sandy soil, the

35

removal of the diesel and products from degradation leads to an increase in electrical

36

conductivity during the first 69 days. This is expected since diesel is electrically insulating. For

37

both soils, the activity of BES reactors is moderately imaged by the inverted conductivity

38

tomogram of the reactor. An increase in current production by electrochemically active bacteria

39

activity corresponds to an increase in conductivity of the reactor.

40

2 ACS Paragon Plus Environment

Page 3 of 31

41

Environmental Science & Technology

TOC Art

42

43 44 45 46 47 48 49 50 51 52 53 54 55 56

3 ACS Paragon Plus Environment

Environmental Science & Technology

57

Page 4 of 31

INTRODUCTION

58

Hydrocarbon contamination due to accidental spills and tank leakages represents a threat

59

to the environment1,2. Besides natural attenuation of the contaminants3,4, a number of engineering

60

solutions have been applied to enhance the remediation process5,6. Recently, bioelectrochemical

61

systems (BESs) were developed and tested for the remediation of hydrocarbon contaminated

62

soils7-13. During the remediation process, the BES provides an environment for electrochemically

63

active bacteria (EAB) to catalyze the oxidation of organic electron donors, e.g., hydrocarbon

64

contamination, and deliver the electrons to the anode. Produced electrons are collected on the

65

anode and transferred through an external circuit to the cathode, where O2 serves as the unlimited

66

terminal electron acceptor9,14. Previous studies showed BESs could enhance hydrocarbon

67

degradation by up to 241% in saturated soils12. In addition, no external energy is required by this

68

process and in fact a small current is even produced that can be used in turn to power small

69

electronic sensors for onsite monitoring.

70

Even though BESs have demonstrated improved performance in hydrocarbon-

71

contaminated soil remediation and provide current as a proxy for monitoring13, they still relies on

72

traditional total petroleum hydrocarbon (TPH) measurements, which are intrusive and time

73

consuming. The TPH measurement involves soil sampling, blending with volatile and toxic

74

chemicals, plus multiple stages of extractions using vortex and sonication. Another issue is that

75

the spatial coverage of TPH measurements is usually very limited. In order to provide more

76

accurate, real time, and easy remediation monitoring, new methods that are fast and non-

77

intrusive are urgently needed.

78

Electrical resistivity survey is a non-invasive geophysical method consisting in imaging

79

the electrical resistivity or electrical conductivity of porous media15. Electrical conductivity is 4 ACS Paragon Plus Environment

Page 5 of 31

Environmental Science & Technology

80

sensitive to the salinity and temperature, water content, and cation exchange capacity (CEC) or

81

specific surface area (SSp) of soils16. During a resistivity survey, an electrical current is injected

82

through current electrodes, and at the same time, differences of electrical potentials are measured

83

at voltage electrodes17-19. Resistivity survey has been applied to track the movement of saline

84

plume20, 21 and contaminants22. For petroleum hydrocarbon contamination, a number of studies

85

have been done to show the potential of resistivity survey monitoring23-25. Several laboratory

86

experiments have been conducted on soil samples to understand the relation between electrical

87

conductivity and hydrocarbon contamination by manually mixing different contents of

88

hydrocarbon products with uncontaminated soil26-28. For non-wetting oils, e.g. diesel29, the

89

measured conductivity is well correlated with hydrocarbon content, with a lower content

90

corresponding to a higher conductivity. This is true for both fresh oil and biodegraded oil taken

91

from field sites28.

92

In the current study, diesel is chosen as the hydrocarbon contaminant and mixed with

93

either sandy or clayey soil in two tanks for characterizing the BES reactor efficiency in

94

accelerating hydrocarbon degradation. When non-conductive diesel is consumed by degradation

95

and removed from soils, it is expected to increase the electrical conductivity of the material. A

96

time lapse resistivity survey is designed to capture these changes during the tank experiment.

97

Our goal is to establish the relationship between the diesel TPH content and measured electrical

98

conductivity variations. At the same time, besides the monitoring of voltage generation from

99

BES reactors, reactors were also evaluated in the resistivity survey in terms of conductivity

100

variation.

101 102

MATERIALS AND METHODS

5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 31

103

Tank Setup. The two tanks used for the analysis had a dimension of 40 cm × 5.5 cm ×

104

20 cm (Figure 1). The BES reactor was placed on the left side of each tank. A column-type BES

105

reactor was characterized by an outer diameter of 4.5 cm. A total of 32 stainless steel electrodes

106

were fixed to the side wall of the tank. These electrodes were connected with a banana socket

107

panel for resistivity surveys. Two sites at x = 5 cm and x = 39 cm which were 1 cm and 35 cm

108

away from the reactor were used as soil sampling locations for TPH measurements. Three

109

locations A, B and C, at x = 7 cm, 14 cm and 30 cm, were chosen as soil electrical conductivity

110

measurements at the end of the experiment. For each of the three locations, soil samples at three

111

different levels were taken.

112

Two different types of soils, sandy and clayey soils, at saturated conditions were

113

employed to test the effectiveness of the BES reactors. Two soils were prepared by mixing

114

natural soils with pure silica (Type 7030 from Unimin Corporation) sand or pure clay (Pottery

115

clay from Stone Leaf Pottery), respectively. Natural soil is needed for indigenous bacterial

116

population. The natural soil characteristics were reported by Lu et al. [2014]13. After mixing,

117

sandy soil (57.14% sand, 28.93% silt and 13.93% clay in dry weight) was characterized as sandy

118

loam and clayey soil (33.12% sand, 35.46% silt and 31.42% clay in dry weight) was

119

characterized as clay loam according to USDA soil texture classification. The porosity of sandy

120

soil is 0.41, and 0.50 for clayey soil. SSp of sandy soil was measured at 1.68 m2/g, and 4.07 m2/g

121

for clayey soil with 5-point BET method. The CEC values30 of sandy and clayey soils are

122

calculated as 0.56 meq/100g and 1.35 meq/100g by assuming two elementary charges per nm2.

123

Diesel was chosen to represent the hydrocarbon contaminants. Diesel was well mixed with soils

124

before putting into tanks. Soils were saturated with artificial groundwater. Artificial groundwater

6 ACS Paragon Plus Environment

Page 7 of 31

Environmental Science & Technology

125

refers to a solution of 4 mM Na+, 0.2 mM Ca2+, 0.1 mM Mg2+, 2.4 mM Cl-, 2 mM HCO3- , and

126

0.1 mM SO42- in deionized water (conductivity ≈ 0.05 S/m at 22oC, pH ≈ 8.2)31.

127

During the whole experiment and in order to minimize evaporation, both tanks were

128

covered with plastic wraps on the surface excluding the area where the reactors were located.

129

The hollow pipe of the reactor was left open for exposure with air, which provides the terminal

130

electron acceptor in the form of O2. The water loss due to evaporation and leakage to the

131

perforated BES reactors was approximately compensated by adding deionized (DI) water twice a

132

week to maintain saturation. The saturated condition was observed from the moisture condition

133

at the surface of each tank.

134

Bioelectrochemical Reactor Operation. For both reactors, carbon cloth with activated

135

carbon catalyst served as air cathode32-34. For the anode, a carbon felt anode was used. Cation

136

Exchange Membrane (CEM) was employed as a separator between cathode and anode. Both

137

reactors were wrapped around a perforated PVC pipe, and the height of the PVC pipe was 25 cm.

138

Detailed information about how to construct reactors can be found in Lu et al. [2014]13.

139

Continuous current generated from the remediation process was monitored and evaluated

140

by an external circuit with a resistor during the whole experiment period. A 1000 Ω external

141

resistor was utilized to start the experiment when a sustainable and relative large voltage

142

(>100mV) was observed in a few days after the start of the experiment35, and at 162nd days it

143

was replaced with a small resistor of 100 Ω to speed up the remediation process.

144

Control experiment of natural biodegration was run in separate containers which had the

145

dimension of 8 cm × 5 cm. There were no BES reactors and geophysical electrodes installed for

146

the control experiments, and therefore TPH variation only reflects the natural biodegration. We

7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 31

147

took TPH measurements from the control experiment at the same time as the active tank

148

experiment for comparison.

149

Soil Sampling. Soil samples were taken throughout the experiment at two distances away

150

(x = 5 cm and x = 39 cm) from the reactor (Figure 1) for TPH measurements36. At each location,

151

5 g of soil from the top, middle and bottom of the tank were sampled and then mixed for TPH

152

measurements. Seven TPH measurements at day 0, 8, 40, 60, 120, 150 and 200 were conducted

153

at each location during the whole experiment. The initial TPH values for sandy and clayey soil

154

were 6960 mg/kg and 6495 mg/kg, respectively.

155

At the end of the experiment, soil samples were taken at three different locations A, B,

156

and C (Figure 1) for electrical conductivity measurement, at each location, lower (z ≈ 1.5 cm~4

157

cm), middle (z ≈ 6.5 cm~9 cm) and upper (z ≈ 11.5 cm~14 cm) levels soil samples were taken.

158

At each location, the average conductivities from all three levels were compared with results

159

from the resistivity tomograms. Soils conductivity was measured in a 1:5 (w/v) soil: deionized

160

(DI) water mixture37. With this method, the initial electrical conductivities for sandy and clayey

161

soils were 558 µS/m and 1109 µS/m, respectively.

162

Soil Column Conductivity Measurements. Measured conductivity tomograms are

163

affected by both surface conductivity and pore fluid conductivity, and these two contributions

164

cannot be separated in a resistivity survey. Therefore besides the tank experiment, soil column

165

experiments were setup to separate the influence of surface conductivity and diesel TPH content

166

for the two different types of soils. A four-electrode configuration (two current electrodes and

167

two voltage electrodes) was adopted. The soil column had a diameter of 3 cm and a length of 11

168

cm. Two stainless steel plates were used as current electrodes at both ends of the column, and

169

Ag-AgCl potential electrodes were inserted in the middle of column with a spacing of 3 cm. The

8 ACS Paragon Plus Environment

Page 9 of 31

Environmental Science & Technology

170

soil samples were saturated in a vacuum chamber. All the measurements were performed with

171

ZEL-SIP04-V02 impedance meter38. The real part of the complex conductivity at 1Hz was

172

chosen to represent electrical conductivity during a resistivity survey for analysis.

173 174 175

According to Waxman and Smits39, electrical conductivity ߪ (in S/m) of the material at low frequency can be represented as ଵ

ߪ = ߪ௪ + ߪ௦ , ி

(1)

176

where F (dimensionless) is the formation factor related to the connected porosity of the soil by a

177

power law relationship known as Archie’s law, ߪ௪ is the pore water conductivity, and ߪ௦ is the

178

surface conductivity associated to electrical conduction in the electrical double layer coating the

179

surface of the grains30. This surface conductivity is controlled by the SSp or the CEC of the clays.

180

Six electrolytes were obtained by mixing DI water and pure NaCl at the following conductivity

181

(at 298 K) ߪ௪ =0.002 S/m, 0.02 S/m, 0.08 S/m, 0.2 S/m, 0.8 S/m and 2.0 S/m.

182

For the TPH values in the soils, four different ratios were tried, TPH=0.0 g/kg, 7.5 g/kg,

183

13 g/kg ad 32 g/kg, and 0.0 g/kg means no diesel was added to the soil. Diesel was added to soil

184

first and artificial groundwater was then added. With the added diesel, it took a longer time to

185

expel air from the mixture (approximately 3 hours in this study) compared with only water as

186

pore fluid (0.6 hour) in a vacuum chamber. Potential loss of diesel due to evaporation was not

187

accounted. In such a short measurement period, there was also no consideration of the natural

188

degradation of diesel contaminants on the electrical conductivity.

189

Resistance Data Collection and Inversion. Resistivity surveys were conducted eight

190

times throughout the experiment (Figure 2), at day 12, 26, 42, 69, 105, 135, 187 and 230, with a

191

two week interval at the beginning of the test and a 1.5 month interval at the end. For each

192

survey, 480 resistance measurements were taken to cover the whole tank with a multichannel

9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 31

193

instrument ABEM LS and survey duration was 15 min. For the 480 measurements, current

194

electrode pairs AB were chosen as 14 horizontal pairs and 3 vertical pairs, i.e. AB=15 and 21, 2

195

and 7, 21 and 27, and 15 and 2. For each current injection pair AB, potential electrodes pairs

196

were chosen with skip one scheme18. For example, for AB=2 and 15, potentials electrode pairs

197

MN=3 and 5, 4 and 6, 5 and 7, and so on.

198

All the 480 measurements were used for inversion. For each tank, the resistances from all

199

eight surveys were inverted together with the active time constraint (ATC) technique40, 41. This

200

technique improves the time-lapse imaging especially when the survey data is corrupted with

201

strong noise. Due to the sensitivity difference of a resistivity survey with higher sensitivity close

202

to the electrode arrays19, we presented the results in terms of conductivity ratio, which better

203

shows the variations. The result from the first survey (at day 12) was treated as a reference. The

204

conductivity ratio between the other seven surveys and the first one was presented for discussion.

205

In this way, a conductivity ratio greater than 1.0 indicates an increase of electrical conductivity.

206

The inversion was performed with a three dimensional finite element code built in the

207

Comsol and Matlab environments41. Because of the smaller tank thickness in the y direction, we

208

assumed no conductivity variation in the y direction. Therefore the inversion results were

209

presented as two dimensional cross-tomograms.

210

211

RESULTS AND DISCUSSION

212

Voltage Generation. Figure 2 shows the voltage profile for the two tanks. Stable voltage

213

was generated in both tanks. At the beginning the electrical voltage has a rising limb due to

214

acclimation of indigenous microbial electrochemical community, and then the voltage reaches a

215

stable value. A longer delay period is observed in the clayey soil to reach the stable value due to 10 ACS Paragon Plus Environment

Page 11 of 31

Environmental Science & Technology

216

slower mass transfer rate compared with the sandy soil from the restrict passage of solutes in

217

clayey soils42. For both tanks, a stable voltage is reached at around 160 mV.

218

After 162-day experiment, a 100 Ω resistor was used to replace the 1000 Ω to allow

219

faster electron flow35. Both reactors show immediate voltage fluctuation in response to the

220

resistor change and then the voltage decreased gradually with time. The immediate surge is

221

caused by the reestablishment of the microbial activity with the smaller resistor.

222

TPH Removal. The results from both tank and control experiments are shown in Figure

223

2 (Y axis on the right). Corresponding to the initial increase of voltage, there is an increase of

224

TPH close to the reactor at x = 5 cm for both tanks. This is due to the establishment of EAB

225

activity with high adsorption of carbon felt anode attracting more diesel contaminants towards

226

the reactor.

227

Generally, TPH has larger removal rate closer to the reactor (x = 5 cm) and relatively

228

lower rate further away (x = 39 cm) 13. However, the difference between these two distances is

229

very small, especially for sandy soil. This proves the relatively large radius of influence under

230

sandy soil condition. A higher TPH removal rate close to the reactor compared to TPH further

231

away creates a diesel concentration gradient. This concentration gradient could drive diesel

232

contaminant to diffuse towards the reactors. Another driving force is the hydraulic gradient

233

created by the perforated BES reactors. The continuous leaking of liquid into the reactors and

234

compensation from the adding DI water on the surface also drives the movement of diesel

235

contaminant. Since sandy soil has a larger hydraulic conductivity which permits a fast movement,

236

the difference between TPH measurements at the two distances is smaller compared to the clayey

237

soil.

11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 31

238

Sandy soil shows a higher removal rate than clayey soil based on the general trend. Both

239

BES experiments have higher TPH removal rate than control experiments, but less evident for

240

clayey soil. Clay soil at x = 39 cm is almost close to the control test. First two TPH

241

measurements before day 8 are indistinguishable with control experiment. From the

242

measurement at day 40, we start to see the difference. At day 60, the difference becomes clearer.

243

In the control experiment, for sandy soil, TPH decrease to 5200 mg/kg which corresponds to 25%

244

removal rate, and 5550 mg/kg for clayey soil with removal rate 15%. With BES reactor, sandy

245

soil has remove 41% and 35% TPH at x = 5 cm and at x = 39 cm respectively which are 64% and

246

40% higher than natural biodegration. As expected, even with BES reactor, the TPH removal

247

rate of clayey soil is not as good as sandy soil. At day 60, it only has 23% and 17% TPH

248

removed at the two distances which is slightly higher than the control test value 15%, only 53%

249

and 13% increase from the control experiments. At the end of the experiments, sandy soil has

250

TPH values of 4100 mg/kg and 4210 mg/kg at x = 5 cm and 39 cm, which correspond to 41%

251

and 39% removal rates based on the initial TPH content. For the control test, it has a removal

252

rate 34%. For clayey soil, they are 4061 mg/kg at x = 5 cm and 4861 mg/kg at 39 cm which

253

correspond to 35% and 25% removal rates. As for the control test, it only removes 22% of the

254

TPH. From both tests, we can see that natural biodegration from the indigenous microorganisms

255

plays an important role during the TPH removal. Therefore, when analyzing the resistivity

256

change from the experiments, the change reflects the contribution from both BES and natural

257

biodegration.

258

Conductivity of Soil Column Measurements. The measured conductivity at 1Hz with

259

different salinities is shown in Figure3a obtaining the influence of surface conductivity. The

260

measured data was fitted with the relationship in Eq. (1). With this fitting, the obtained formation

12 ACS Paragon Plus Environment

Page 13 of 31

Environmental Science & Technology

261

factors are 4.4 and 2.9 for sandy and clayey soil respectively. Clayey soil indicates a larger

262

surface conductivity ߪ௦ = 0.051 S/m and sandy soil has a relatively smaller surface conductivity

263

ߪ௦ = 0.022 S/m. With both fluid conductivity and surface conductivity affect the surveyed

264

conductivity values, a larger surface conductivity corresponds to the higher clay content and also

265

indicates the smaller influence from the change of pore water conductivity, which makes

266

conductivity variation less evident when removing diesel from the soil.

267

In Figure3b, ߪ଴ represents the electrical conductivity of soil sample without diesel

268

contaminant (0.0 g/kg). Both sandy and clayey soils have a decreasing conductivity trend with an

269

increasing TPH content as indicated in Figure3b. However, clayey soil shows a smaller decrease

270

trend with only a 5% drop at 32 g/kg TPH due to a larger surface conductivity. Conductivity of

271

sandy soil drops by 6% even at 7.2 g/kg TPH, and 22% drop when TPH reaches 32 g/kg.

272

Conductivity Ratio Tomograms. Figure 4 is the conductivity ratio tomograms for both

273

tanks where the sub titles show the survey day after the start of the experiment (12 being the

274

reference). For the analysis of the soil conductivity variation, we exclude the reactor region

275

marked by the dashed rectangles in black, and the conductivity variation of the reactors will be

276

evaluated separately. Since we didn’t have a parallel tank experiment to monitor the change of

277

conductivity with diesel removal under natural biodegration, the observed electrical conductivity

278

variation represents the contribution from both natural and enhanced biodegradations from BES

279

reactors.

280

Compared to clayey soil results in Figure 4b, there is a distinct variation trend for

281

saturated sandy soil (Figure 4a), which indicates the removal of TPH from the tanks. In Figure 4a,

282

in the region closer to the surface, we observed a continuous drop of conductivity indicated by

283

the blue color. There are two reasons for this drop. There is a buoyance effect from the diesel

13 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 31

284

contaminants. Diesel has a density 0.832 kg/L and therefore has a tendency to move up above the

285

pore water phase. A few weeks after the start of the experiment, there was a visual existence of

286

diesel spots on the surface of the sandy tank. Another reason is due to added DI water on the

287

surface to keep the soil saturated. Even though the tank is covered with plastic wrap to minimize

288

evaporation, the tank is not completely sealed and water is still lost from the surface under

289

Colorado’s dry climate. DI water is added to the surface twice a week. The added DI water only

290

has a conductivity of 10-4 S/m and therefore this low conductive pore water reduces the

291

conductivity of the sandy soil. The added DI water is believed to be the main reason for the

292

observed decrease of conductivity on the top layer of the soil.

293

From the tomograms, in the region below the surface, there is an increase of electrical

294

conductivity at first, which is then followed by a decrease trend. However, the variation is not

295

uniformly distributed inside the tanks. This is potentially from the uneven distribution of diesel

296

contaminants inside the tank. In some regions, there is an increase of conductivity by 30% (Day

297

69 in Figure4a) even with only 3 g/kg TPH removal, from approximately 7g/kg to 4g/kg (see

298

Figure2a the second and fourth TPH measurements). This is a much bigger variation than the

299

conductivity measurements shown in Figure 3b. There may exist some higher diesel content

300

spots in the tank, and the organic acids and bicarbonate accumulation from the degradation also

301

contribute to the increase of conductivity which is not accounted in the soil column

302

measurements. As discussed by Allen et al. [2007]43, an increase of culturable hydrocarbon-

303

degrading microbes is often observed in petroleum-contaminated sites. The presence of certain

304

bacterial phyla corresponds to an increase of conductivity. Microbial activity increase from BES

305

and natural biodegration can alter the subsurface electrolytic and interfacial properties of porous

306

media. This alternation can influence the electrical conductivity25, 44.

14 ACS Paragon Plus Environment

Page 15 of 31

Environmental Science & Technology

307

For the clayey soil, there was less variation in the conductivity ratio tomograms (Figure

308

4b). We also observed some diesel spots on the surface of the tank after a few weeks of operation.

309

However, the conductivity decrease trend close to the surface is not as clear as the sandy tank.

310

The DI water moves slower in the clayey soil due to the low hydraulic conductivity, and

311

therefore does not lead to a significant decrease to conductivity as indicated in the sandy soil (see

312

Figure 4a).

313

There is a significant increase of conductivity close to the reactor at 10 cm at the late

314

stage of the experiment. This happens after the change of smaller resistor (162 day) with a

315

concomitant reactor activity speedup. This increase is related to more consumption of diesel

316

contaminants close to the reactor.

317

Average Conductivity and TPH Variations. Because of the heavy soil sampling

318

activity at the two locations x = 5 cm and 39 cm (Figure1), the soil was highly disturbed in these

319

two regions. Therefore, we only consider the area inside the white dashed rectangle for

320

calculation of average conductivity variation (the area is indicated by white dashed rectangles in

321

the Figure 4a and Figure 4b Day 26). Solid lines in Figure 5a and 5b show the average

322

conductivity ratio variation for both tanks for the area inside the white dashed rectangles

323

indicated in Figure 4 (The last survey at 230 days is not used since it is beyond the last TPH

324

measurement). Figure 5a and 5b also have the average TPH ratio from the two locations (x = 5

325

cm and 39 cm) in dashed lines. The average TPH is calculated by taking the average of the TPH

326

measurements at the two locations to represent the TPH variation inside the tank. Since the

327

conductivity survey and TPH measurements are not taken at the same time, the average TPH

328

values at the same time of the conductivity survey are linearly interpolated. All TPH ratios are

329

referred to the interpolated TPH value at 12th day when the first conductivity survey taken.

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 31

330

As shown by the average conductivity, the conductivity increases first and then

331

decreases. The average conductivity ratio increases by about 12% at the fourth survey 69 days

332

after the start of the experiment and then slowly decreases by 5% at the seventh survey (187 days)

333

compared to the initial state. For the average TPH ratio result of sandy soil, there is a significant

334

drop for the first three results in Figure5a, and then followed by a slight increase and decrease

335

variation trend. The TPH decrease trend of the first three points corresponds to an increase of

336

average conductivity. This is consistent with previous research conclusions that a lower TPH will

337

reflect a higher conductivity value26-29.

338

After the first three data points in Figure 5a, we could not find this relation anymore

339

because there is another reason responsible for the observed variation of average electrical

340

conductivity. Indeed, there are two processes occurring at the same time during remediation

341

causing the change of electrical conductivity: (i) the removal of diesel contaminants and (ii) the

342

leakage of the pore fluid into the reactor. Nonconductive diesel is removed from soil through

343

biodegration and leads to organic acids and bicarbonate accumulation. This leads to the increase

344

of soil conductivity. At the same time, because the whole reactor is wrapped around a perforated

345

PVC pipe to maximize the cathode contact with air, pore fluid leaks through the perforated pipe

346

and creates a hydraulic gradient towards the reactor which is enhanced by the added DI water on

347

the surface. As a result of this head gradient, the removal of fluid and icons loss in the fluid from

348

the tank decrease conductivity of the soils. The saturation maintained by adding low conductive

349

DI water on the surface of the tanks could not compensate the decrease of conductivity.

350

The effect of leakage enhances both water flow and solute transport in this experiment.

351

Compared to solute transport under natural hydraulic gradient, e.g. 1%, which only moves a few

352

meters per year45, BES reactor itself creates a favorite situation for the transport of contamination.

16 ACS Paragon Plus Environment

Page 17 of 31

Environmental Science & Technology

353

In the tank experiment, if we assume the water table is approximately on the soil surface, and it

354

decrease to about 5 cm below the surface where it starts the perforated part of the PVC pipe. The

355

average hydraulic gradient can be calculated

356

as infiltration as what naturally happens from predipitation23, and the sandy soil example shows

357

us that it can help the movement of solute, and therefore it will help the remediation process46.

ଷହିସ.ହ ହ

≈14%. The act of adding water can be treated

358

For the sandy soil, the conductivity trend confirms that higher TPH removal rate is more

359

prevailing than the leakage of fluid at early period. At later period, the TPH removal rate slows

360

down as shown by the TPH variation. The decrease of conductivity from leakage is more

361

dominant which causes the decrease of surveyed conductivity. The smaller surface conductivity

362

for the sandy soil (ߪ௦ =0.022S/m) supports more influence of pore water in this study.

363

For the clayey soil in Figure 5b, we only observed a decrease trend in conductivity,

364

reaching 10% at 187th days with respect to the reference. The TPH ratio also continues dropping

365

throughout the experiment. The removal of diesel does not lead to an increase of conductivity at

366

the beginning of experiment compared to the results in Figure 5a for the sandy soil. As

367

confirmed by soil column results, this is due to the large surface conductivity of the clayey soil

368

(ߪ௦ = 0.051 S/m). Leakage of fluid is more important for conductivity variation throughout the

369

whole experiment. The removal of fluid decreases fluid conductivity, and therefore what is

370

shown in Figure 4a is a steady decrease.

371

Lateral Conductivity Variations. At the end of the experiment, we took soil samples

372

from the three locations A, B, and C (Figure 1) to evaluate the spatial variation away from

373

reactors. The conductivity ratio from the last survey at these three locations was also extracted

374

from the inverted tomograms. The extracted areas are indicated by three solid white rectangles in

375

Figure 4. Figure 5c is the comparison between sampled conductivity ratios and inverted values at

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 31

376

8th survey from the two tanks. In Figure 5c, from both soil samples, there is a decrease trend

377

away from the reactor at locations A to C, indicating the loss of ions during the leakage process

378

throughout the duration of the experiment and its influence at the later period of the experiment.

379

This is supported by the tomograms of sandy soil figure4a. The increase of conductivity closer to

380

the reactor (indicated by red color) is always higher than the value further away, especially in the

381

last tomogram at Day 230.

382

The cross plot between inversion results and measured conductivity values shows a linear

383

relationship which means the geophysical survey can capture the conductivity variation in the

384

lateral direction. The relationship is not 1:1. The first reason is because of the uneven diesel

385

distribution in the tank initially, and there is only one initial value taken from each tank. When

386

taking the ratio with respect to the initial condition for all the samples, it introduces some

387

uncertainty. For the numerical inversion results, the ratio is taken with respect to the conductivity

388

in A, B, and C locations from the first survey at 12th day respectively. Another reason is the way

389

we measured electrical conductivity by drying the soil first and then mixing with DI water with

390

1:5 ratio37. With this method, the measured conductivity from soil samples can only be used

391

relatively compared to the in situ measurements from resistivity surveys.

392

Relationship between Electrical Conductivity and Voltage in BES Reactors. There is

393

no soil existing in the region of reactors with a hollow pipe design. Therefore, the reactor

394

conductivity variations reflect the activity on the anode of the reactor itself rather than bulk

395

conductivity of the soil. We analyzed conductivity variation of both reactors separately. During

396

the experiment, we also observed a conductivity variation for the reactors from the inversion

397

results. As shown in Figure 4a, the conductivity ratio for the reactor was closer to 1.0 first, and

398

then became smaller to 0.4 in blue color for the 6th survey Day 135. It then kept at the same

18 ACS Paragon Plus Environment

Page 19 of 31

Environmental Science & Technology

399

range. If we see the voltage curve in Figure 2a, the first survey on the 12th day was taken when

400

its magnitude reaches the peak, and is then followed by a gradual decrease. The conductivity

401

ratio and voltage ratio from the reactors are plotted in Figure 5d. For the voltage ratio, the

402

magnitudes 316 mV and 83 mV for sandy and clayey soils at day 12 was utilized as a reference.

403

As indicated by Figure 5d, it follows a moderate linear trend with a scattered pattern for

404

the clayey soil. The linear trend supports that a higher voltage corresponds to a more conductive

405

reactor in the inverted tomograms. The mechanism behind this linear trend can be related to the

406

activity of EAB. When the voltage of the BES reactor is high, that means the EAB activity is

407

high. More diesel contaminants are consumed, and therefore more electrons are produced. This

408

leads to the creation of a conductive anomaly during the resistivity survey.

409

We have demonstrated a useful non-invasive monitoring technique for diesel removal in

410

a tank experiment using resistivity survey. An increase of conductivity corresponds to a decrease

411

of TPH content in sandy soil at the early period of the experiment, but not for clay soil due to the

412

high surface conductivity, and these tank results are further supported by separate soil column

413

measurements with different TPH contents. The leakage of the fluid through the perforated BES

414

reactors causes a decrease of conductivity temporally indicated by the inverted conductivity

415

tomograms and spatially supported by the soil samples at the end of the experiment, and fluid

416

leakage also acts as a driving force for the movement of diesel contaminant towards the reactors.

417

Besides the traditional voltage monitoring of BES reactors, the activity of the BES reactors is

418

characterized by the magnitude of the conductivity in the inverted tomograms as well. A more

419

active BES reactor corresponds to an increase of conductivity in the conductivity tomograms.

420

Even though these results are only from a lab scale tank experiment, the encouraging results

421

support the usage of geophysical survey during contamination remediation monitoring. Further

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 31

422

large scale field tests are still needed to evaluate this technique under more complicated

423

environment to understand more influence parameters of electrical conductivity.

424 425 426

ACKNOWLEDGEMENTS

427

This work is funded by the Office of Science (BER), US. Department of Energy (awards DE-

428

SC0007118) and Chevron Energy Technology Company (grants CW852844 and K26625). We

429

also thank the specific surface area measurement from Dr. Ryan Richards’ lab of Colorado

430

School of Mines.

431

20 ACS Paragon Plus Environment

Page 21 of 31

Environmental Science & Technology

432

REFERENCES

433

1. Smith, S. E.; Hult, M. F. Hydrogeologic data collected from a crude-oil spill site near Bemidji,

434

Minnesota, 1983-91. 93-496,1993, https://pubs.er.usgs.gov/publication/ofr93496.

435

2. U.S. Environmental Protection Agency, Leaking Underground Storage Tank (LUST) Trust

436

Fund, 2014, http://www.epa.gov/oust/ltffacts.htm.

437

3. Cozzarelli, I. M.; Bekins, B. A.; Baedecker, M. J.; Aiken, G. R.; Eganhouse, R. P.; Tuccillo,

438

M. E., Progression of natural attenuation processes at a crude-oil spill site: I. Geochemical

439

evolution of the plume. J. Contam. Hydrol. 2001, 53, (3–4), 369-385.

440

4. Bekins, B. A.; Cozzarelli, I. M.; Godsy, E. M.; Warren, E.; Essaid, H. I.; Tuccillo, M. E.,

441

Progression of natural attenuation processes at a crude oil spill site: II. Controls on spatial

442

distribution of microbial populations. J. Contam. Hydrol. 2001, 53, (3–4), 387-406.

443

5. Riser-Roberts, E., Remediation of petroleum contaminated soils: biological, physical, and

444

chemical processes. CRC press: Lewis Publishers, NY, 1998.

445

6. Wang, H.; Luo, H.; Fallgren, P. H.; Jin, S.; Ren, Z. J., Bioelectrochemical system platform for

446

sustainable environmental remediation and energy generation. Biotechnol. Adv. 2015, 33, (3–4),

447

317-334.

448

7. Morris, J. M.; Jin, S.; Crimi, B.; Pruden, A., Microbial fuel cell in enhancing anaerobic

449

biodegradation of diesel. Chem. Eng. J. 2009, 146, (2), 161-167.

450

8. Zhang, T.; Gannon, S. M.; Nevin, K. P.; Franks, A. E.; Lovley, D. R., Stimulating the

451

anaerobic degradation of aromatic hydrocarbons in contaminated sediments by providing an

452

electrode as the electron acceptor. Environ. Microbiol. 2010, 12, (4), 1011-1020.

21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 31

453

9. Wang, X.; Cai, Z.; Zhou, Q.; Zhang, Z.; Chen, C., Bioelectrochemical stimulation of

454

petroleum hydrocarbon degradation in saline soil using U-tube microbial fuel cells. Biotechnol.

455

Bioeng. 2012, 109, (2), 426-433.

456

10. Morris, J. M.; Jin, S., Feasibility of using microbial fuel cell technology for bioremediation

457

of hydrocarbons in groundwater. J. Environ. Sci.Health, Part A: Environ. Sci. Eng. 2008, 43, (1),

458

18-23.

459

11. Rakoczy, J.; Feisthauer, S.; Wasmund, K.; Bombach, P.; Neu, T. R.; Vogt, C.; Richnow, H.

460

H., Benzene and sulfide removal from groundwater treated in a microbial fuel cell. Biotechnol.

461

Bioeng. 2013, 110, (12), 3104-3113.

462

12. Lu, L.; Yazdi, H.; Jin, S.; Zuo, Y.; Fallgren, P. H.; Ren, Z. J., Enhanced bioremediation of

463

hydrocarbon-contaminated soil using pilot-scale bioelectrochemical systems. J. Hazard. Mater.

464

2014, 274, (0), 8-15.

465

13. Lu, L.; Huggins, T.; Jin, S.; Zuo, Y.; Ren, Z. J., Microbial Metabolism and Community

466

Structure in Response to Bioelectrochemically Enhanced Remediation of Petroleum

467

Hydrocarbon-Contaminated Soil. Environ. Sci. Technol. 2014, 48, (7), 4021-4029.

468

14. Wang, H.; Ren, Z. J., A comprehensive review of microbial electrochemical systems as a

469

platform technology. Biotechnol. Adv. 2013, 31, (8), 1796-1807.

470

15. Loke, M. H.; Chambers, J. E.; Rucker, D. F.; Kuras, O.; Wilkinson, P. B., Recent

471

developments in the direct-current geoelectrical imaging method. J. Appl. Geophys. 2013, 95,

472

135-156.

473

16. Revil, A.; Cathles, L. M.; Losh, S.; Nunn, J. A., Electrical conductivity in shaly sands with

474

geophysical applications. J. Geophys. Res. B: Solid Earth 1998, 103, (B10), 23925-23936.

22 ACS Paragon Plus Environment

Page 23 of 31

Environmental Science & Technology

475

17. Daily, W.; Ramirez, A., Electrical resistance tomography during in-situ trichloroethylene

476

remediation at the Savannah River Site. J. Appl. Geophys. 1995, 33, (4), 239-249.

477

18. Slater, L.; Binley, A. M.; Daily, W.; Johnson, R., Cross-hole electrical imaging of a

478

controlled saline tracer injection. J. Appl. Geophys. 2000, 44, (2–3), 85-102.

479

19. Kemna, A.; Kulessa, B.; Vereecken, H., Imaging and characterisation of subsurface solute

480

transport using electrical resistivity tomography (ERT) and equivalent transport models. J.

481

Hydrol. 2002, 267, (3–4), 125-146.

482

20. Singha, K.; Gorelick, S. M., Saline tracer visualized with three-dimensional electrical

483

resistivity tomography: Field-scale spatial moment analysis. Water Resour. Res. 2005, 41,

484

W05023.

485

21. Cassiani, G.; Bruno, V.; Villa, A.; Fusi, N.; Binley, A. M., A saline trace test monitored via

486

time-lapse surface electrical resistivity tomography. J. Appl. Geophys. 2006, 59, (3), 244-259.

487

22. Revil, A.; Skold, M.; Karaoulis, M.; Schmutz, M.; Hubbard, S. S.; Mehlhorn, T. L.; Watson,

488

D. B., Hydrogeophysical investigations of the former S-3 ponds contaminant plumes, Oak Ridge

489

Integrated Field Research Challenge site, Tennessee. Geophysics 2013, 78, (4), EN29-EN41.

490

23. Sauck, W. A., A model for the resistivity structure of LNAPL plumes and their environs in

491

sandy sediments. J. Appl. Geophys. 2000, 44, (2–3), 151-165.

492

24. Atekwana, E.; Werkema, D.; Duris, J.; Rossbach, S.; Atekwana, E.; Sauck, W.; Cassidy, D.;

493

Means, J.; Legall, F., In‐situ apparent conductivity measurements and microbial population

494

distribution at a hydrocarbon‐contaminated site. Geophysics 2004, 69, (1), 56-63.

495

25. Atekwana, E. A.; Atekwana, E.; Legall, F. D.; Krishnamurthy, R. V., Biodegradation and

496

mineral weathering controls on bulk electrical conductivity in a shallow hydrocarbon

497

contaminated aquifer. J. Contam. Hydrol. 2005, 80, (3–4), 149-167.

23 ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 31

498

26. Schmutz, M.; Revil, A.; Vaudelet, P.; Batzle, M.; Viñao, P. F.; Werkema, D. D., Influence of

499

oil saturation upon spectral induced polarization of oil-bearing sands. Geophys. J. Int. 2010, 183,

500

(1), 211-224.

501

27. Revil, A.; Schmutz, M.; Batzle, M., Influence of oil wettability upon spectral induced

502

polarization of oil-bearing sands. Geophysics 2011, 76, (5), A31-A36.

503

28. Abdel Aal, G. Z.; Atekwana, E. A., Spectral induced polarization (SIP) response of

504

biodegraded oil in porous media. Geophys. J. Int. 2014, 196, (2), 804-817.

505

29. Schwartz, N.; Huisman, J. A.; Furman, A., The effect of NAPL on the electrical properties of

506

unsaturated porous media. Geophys. J. Int. 2012, 188, (3), 1007-1011.

507

30. Revil, A., Spectral induced polarization of shaly sands: Influence of the electrical double

508

layer. Water Resour. Res. 2012, 48, (2), W02517.

509

31. Hort, R. D.; Revil, A.; Munakata-Marr, J.; Mao, D., Evaluating the potential for quantitative

510

monitoring of in situ chemical oxidation of aqueous-phase TCE using in-phase and quadrature

511

electrical conductivity. Water Resour. Res. 2015, 51, (7), 5239-5259.

512

32. Yazdi, H.; Alzate-Gaviria, L.; Ren, Z. J., Pluggable microbial fuel cell stacks for septic

513

wastewater treatment and electricity production. Bioresource Technol. 2015, 180, 258-263.

514

33. Haeger, A.; Forrestal, C.; Xu, P.; Ren, Z. J., High performance spiral wound microbial fuel

515

cell with hydraulic characterization. Bioresource Technol. 2014, 174, 287-293.

516

34. Zhang, F.; Cheng, S.; Pant, D.; Bogaert, G. V.; Logan, B. E., Power generation using an

517

activated carbon and metal mesh cathode in a microbial fuel cell. Electrochem. Commun. 2009,

518

11, (11), 2177-2179.

24 ACS Paragon Plus Environment

Page 25 of 31

Environmental Science & Technology

519

35. Menicucci, J.; Beyenal, H.; Marsili, E.; Veluchamy; Demir, G.; Lewandowski, Z., Procedure

520

for Determining Maximum Sustainable Power Generated by Microbial Fuel Cells. Environ. Sci.

521

Technol. 2006, 40, (3), 1062-1068.

522

36. U.S. Environmental Protection Agency, Non-Halogenated Organics usingGC/FID, Method

523

8015D. In June 2003.

524

37. Kettler, T.; Doran, J. W.; Gilbert, T., Simplified method for soil particle-size determination

525

to accompany soil-quality analyses. Soil Sci. Soc. Am. J. 2001, 65, (3), 849-852.

526

38. Zimmermann, E.; Kemna, A.; Berwix, J.; Glaas, W.; Münch, H. M.; Huisman, J. A., A high-

527

accuracy impedance spectrometer for measuring sediments with low polarizability. Meas. Sci.

528

Technol. 2008, 19, (10), 105603.

529

39. Waxman, M. H.; Smits, L. J. M., Electrical Conductivities in Oil-Bearing Shaly Sands. Soc.

530

Petr. Eng. J. 1968, 8, (2), 107-122.

531

40. Karaoulis, M. C.; Kim, J. H.; Tsourlos, P. I., 4D active time constrained resistivity inversion.

532

J. Appl. Geophys. 2011, 73, (1), 25-34.

533

41. Mao, D.; Revil, A.; Hort, R. D.; Munakata-Marr, J.; Atekwana, E. A.; Kulessa, B.,

534

Resistivity and self-potential tomography applied to groundwater remediation and contaminant

535

plumes: Sandbox and field experiments. J. Hydrol. 2015, 530, 1-14.

536

42. Malusis, M. A.; Shackelford, C. D.; Olsen, H. W., Flow and transport through clay

537

membrane barriers. Eng. Geol. 2003, 70, (3–4), 235-248.

538

43. Allen, J. P.; Atekwana, E. A.; Atekwana, E. A.; Duris, J. W.; Werkema, D. D.; Rossbach, S.,

539

The Microbial Community Structure in Petroleum-Contaminated Sediments Corresponds to

540

Geophysical Signatures. Appl. Environ. Microb. 2007, 73, (9), 2860-2870.

25 ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 31

541

44. Williams, K. H.; Ntarlagiannis, D.; Slater, L. D.; Dohnalkova, A.; Hubbard, S. S.; Banfield, J.

542

F., Geophysical Imaging of Stimulated Microbial Biomineralization. Environ. Sci. Technol. 2005,

543

39, (19), 7592-7600.

544

45. Panday, S.; Wu, Y. S.; Huyakorn, P. S.; Springer, E. P., A three-dimensional multiphase flow

545

model for assessing NAPL contamination in porous and fractured media, 2. Porous medium

546

simulation examples. J. Contam. Hydrol. 1994, 16, (2), 131-156.

547

46. Li, X.; Wang, X.; Ren, Z. J.; Zhang, Y.; Li, N.; Zhou, Q., Sand amendment enhances

548

bioelectrochemical remediation of petroleum hydrocarbon contaminated soil. Chemosphere 2015,

549

141, 62-70.

550

26 ACS Paragon Plus Environment

Page 27 of 31

Environmental Science & Technology

Figures

551

552

Reactor

Electrodes

TPH locations

EC location A

EC location B

EC location C

12 7

18

1 13

32 26

20

Z (cm)

8

27 3

35

21

) cm Y(

15 3 0

5

25 15

X (cm)

553

554 555 556 557 558

Figure 1. Tank configuration used for the experiment. Current from the BES reactor was monitored by a data acquisition system. A total of 32 resistance electrodes (black cubes) were attached on the side wall of the tank. Two sites (x = 5 cm and 39 cm, light blue) were used as TPH sampling locations during the whole experiment. Finally, 9 soil samples were taken from 3 locations A, B, and C at the end of experiments for electrical conductivity (EC) measurements.

559

27 ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 31

560

561

a) Sandy soil sandbox experiment TPH-5 cm TPH-39 cm Control

300

8000

200 6000 100

0

0

50

100

150

200

TPH (mg/kg)

Voltage(mV)

400

4000

Time(day)

b) Clayey soil sandbox experiment

8000

TPH-5 cm TPH-39 cm

300

Control

200 6000 100

0

0

50

100

150

200

TPH (mg/kg)

Voltage(mV)

400

4000

562

Time(day)

563

Figure 2. Time-series of the voltage generation of BES reactors and TPH degradation during the experiment. a) Sandy soil. b) Clayey soil. Eight resistivity surveys were conducted at day 12, 26, 42, 69, 105, 135, 187 and 230. They are indicated the inverted filled triangles. TPH variation measured for each tank and the results are referred to the Y axis on the right-side of the plot. Seven TPH measurements were taken at day 0, 8, 40, 60, 120, 150 and 200. Control experiments of the two soils were prepared in different containers without BES reactors. It represents the natural biodegration. The vertical solid lines at day 162 indicate the change of resistor from 1000Ω to 100Ω.

564 565 566 567 568 569 570

28 ACS Paragon Plus Environment

Page 29 of 31

Environmental Science & Technology

a) Conductivity data at saturation Conductivity of the soil (S/m)

0.60 0.40 0.20

Sandy soil Clayey soil

0.05

0.01 -3 10

Conductivity ratio (unitless)

574 575 576 577 578 579 580 581

100

Clayey soil 0.95

0.90

0.85

Sandy soil 0.80 0

573

10-1

b) Conductivity ratio vesus TPH

1.00

571 572

10-2

Pore water conductivity (S/m)

10

20

30

TPH (g/kg)

Figure 3: Petrophysical measurements. a) The solid dots show the electrical conductivity at 1Hz with different salinities for sandy and clayey soils. The data were fitted by a linear conductivity equation given by Eq. (1) of the main text. The fit of the data with this equation is materialized by the solid lines. Clayey soil has a bigger surface conductivity (ߪ௦ = 0.051 S m-1) than sandy soil (ߪ௦ = 0.022 S m-1). A bigger surface conductity will lower the influence of pore fluid conductivity created by diesel contaminants removal in this experiment. b) Electrical conductivity ratio σ / ߪ଴ decreases with increasing TPH. ߪ଴ denotes the soil conductivity with no diesel added in the soil sample.Sandy soil has a larger decrease rate than clayey soil due to a smaller surface conductivity as indicated the in Figure 3a.

582 29 ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 31

583

a) Sandy soil sandbox tomograms Day 42

R E A 10 C T 5 O R

Z (cm)

Day 69

R E A C T O R

15

R E A C T O R

Day 105

Z (m)

Day 26

R E A C T O R

0

Day 187

Z (cm)

0

Day 230

Z (m)

R E A C T O R

R E A C T O R

Conductivity ratio Z (m)

Day 135 R E A 10 C T 5 O R 15

0

0

0.3

0.2

0.5

0.7

0.9

1.1

1.2

1.4

x

5

10

15

20

X (cm)

25

30

35

b) Clayey soil sandbox tomograms Day 42

Z (cm)

Day 69

R E A C T O R

R E A C T O R

Day 105

Z (m)

Day 26 R E A 10 C T 5 O R 15

R E A C T O R

0

Z (cm)

0

584

585 586 587 588 589 590 591 592 593

Day 230 R E A C T O R

Conductivity ratio Z (m)

Day 187 R E A C T O R

Z (m)

Day 135 R E A 10 C T 5 O R 15

0

0

0.2

0.3

0.5

0.7

0.9

1.1

1.2

1.4

x

5

10

15

20

X (cm)

25

30

35

Figure 4. Electrical conductivity ratio tomograms with respect to the first tomogram of electrical conductivity taken at day 12. a) Sandy soil. b) Clayey soil. Blue color (1.0) indicates the removal of diesel contaminants. Title of each subfigure denotes the measurement day after the start of the experiment. The dashed rectangle (on the left side) indicates the location of the reactor. The dashed rectangles in white in the first figure show the area used to determine the average electrical conductivity used in Figure 5a and 5b. The solid rectangles (last plot) are three soil sample locations A, B and C at the end of the test. The solid black lines are the locations of surface and borehole electrode shown in Figure 1.

594

30 ACS Paragon Plus Environment

Page 31 of 31

Environmental Science & Technology

595

596

c) Measured vesus inverted conductivity ratios

a) Conductivity and TPH ratio variations: sandy soil 1.2

1.05

Conductivity

0.8

TPH

0.95

TPH ratio

Conductivity ratio

0.9

0.7

0.85 10

40

70

100

130

160

Meausurements ratio

1.15

0.6 190

A

Sand

B C

0.9

A

Clay

B C

0.6

0.3 0.6

0.9

1.2

1.5

Inversion ratio

Time (day)

b) Conductivity and TPH ratio variations: clayey soil

d) Comparison between conductivity and voltage ratios

1.15

2.5

0.8

TPH 0.95

598 599 600 601 602 603 604 605 606 607 608

Sand Clay

1.5

1

0.5

0.85 10

597

0.7

Conductivity

Voltage ratio

2 1.05

TPH ratio

Conductivity ratio

0.9

40

70

100

130

160

0.6 190

0

0

0.5

Time (day)

1

1.5

2

2.5

Conductivity ratio

Figure 5. Electrical conductivity and TPH ratio variations for a) sandy soil and b) clayey soil. Electrical conductivity ratio is calculated from tomograms in the dashed white rectangle area in Figure 4 (except last survey at 230th day), and TPH ratio variation is averaged from the two TPH sampling locations at x = 5 cm and x = 39 cm. For sandy soil, an increase of conductivity corresponds to a decrease of TPH for the first three data points. c) Comparison of the measured conductivity ratio of the sandy and clayey soil samples (measured from tanks) versus the inverted conductivity ratio at the same locations A, B, and C from the inverted tomograms. This linear relation shows that the laboratory measurements are consistent with the tomogram data. d) Linear relationship between the inverted conductivity ratio and the measured voltage ratio for sandy soil (red dots) and clayey soil (green dots) for BES reactors indicated by the dashed rectangles in black in Figure 4.

609

31 ACS Paragon Plus Environment