Organic Chemical Characterization and Mass Balance of a

Nov 14, 2017 - †Department of Civil, Environmental, and Architectural Engineering and ‡Center for Environmental Mass Spectrometry, University of C...
0 downloads 9 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Organic Chemical Characterization and Mass Balance of a Hydraulically Fractured Well: From Fracturing Fluid to Produced Water Over 405 Days James S. Rosenblum, E. Michael Thurman, Imma Ferrer, George R. Aiken, and Karl G. Linden Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03362 • Publication Date (Web): 14 Nov 2017 Downloaded from http://pubs.acs.org on November 17, 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 41

Environmental Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

1

Organic Chemical Characterization and Mass Balance of a Hydraulically Fractured

2

Well: From Fracturing Fluid to Produced Water Over 405 Days

3 4

James Rosenbluma, E. Michael Thurmanb, Imma Ferrerb, George Aikenc, and Karl G.

5

Lindena*

6 7

a

8

Colorado, UCB 428, Boulder, CO 80309, United States

9

b

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

Center for Environmental Mass Spectrometry, University of Colorado, UCB 428,

10

Boulder, CO 80309, United States

11

c

U.S. Geological Survey, Boulder, Colorado, 80309, United States

12 13

*Corresponding author: [email protected]

14 15

Abstract

16

A long-term field study (405 days) of a hydraulically fractured well from the Niobrara

17

Formation in the Denver-Julesburg Basin was completed. Characterization of organic

18

chemicals used in hydraulic fracturing and their changes through time, from the pre-

19

injected fracturing fluid to the produced water, was conducted. The characterization

20

consisted of a mass balance by dissolved organic carbon (DOC), volatile organic analysis

21

by gas chromatography/mass spectrometry, and nonvolatile organic analysis by liquid

22

chromatography/mass spectrometry. DOC decreased from 1500 mg/L in initial flowback

23

to 200 mg/L in the final produced water. Only ~11% of the injected DOC returned by the

1

ACS Paragon Plus Environment

Page 2 of 41

Page 3 of 41

Environmental Science & Technology

24

end of the study, with this 11% representing a maximum fraction returned since the

25

formation itself contributes DOC. Furthermore, the majority of returning DOC was of the

26

hydrophilic fraction (60 - 85%). Volatile organic compound analysis revealed substantial

27

concentrations of individual BTEX compounds (0.1 – 11 mg/l) over the 405-day study.

28

Nonvolatile organic compounds identified were polyethylene glycols (PEGs),

29

polypropylene glycols (PPG), linear alkyl-ethoxylates, and triisopropanolamine (TIPA).

30

The distribution of PEGs, PPGs, and TIPA, and their ubiquitous presence in our samples

31

and the literature, illustrate their potential as organic tracers for treatment operations or in

32

the event of an environmental spill.

33 34

1. Introduction

35

Unconventional oil and gas development and hydraulic fracturing has grown

36

steadily in recent years, with the advent of directional drilling and fracturing fluid

37

mixtures that permit economic recovery of shale resources. This is demonstrated by the

38

major U.S. shale basins (lower 48 states) accounting for 92% of the domestic oil

39

production growth and all of the domestic natural gas production growth from 2011 to

40

2014.1 The extraction of shale resources takes place by drilling both vertical and

41

horizontal wells, followed by formation perforation, and the injection of fracturing fluids.

42

Fracturing fluids are primarily water (~90%) and proppant (~9%), along with various

43

chemical additives (~1%) that enable operators to pump and deliver proppant through the

44

borehole to elongate and prop-open fissures, creating surface area and continued access

45

to oil and gas resources.2,3 Although the extraction of shale resources has had significant

46

economic impacts, there are numerous environmental and public health concerns, 2

ACS Paragon Plus Environment

Environmental Science & Technology

47

especially as it relates to water use and chemical additives associated with hydraulic

48

fracturing.4–7

49 50

Water plays a key role in unconventional oil and gas development, from drilling

51

to hydraulic fracturing. Hydraulic fracturing uses anywhere from 7,000,000 to

52

50,000,000 L of water for a single well, depending on geological characteristics,

53

horizontal length, operator, and type of fracturing method (i.e. slick water, gel, or hybrid-

54

frac).8–10 The chemicals associated with hydraulic fracturing are also of environmental

55

and public health interest, especially considering that 1% can amount to 500,000 L of

56

chemicals used per well (vol:vol). The fracturing fluids are injected into the formation

57

and return to the surface along with the hydrocarbons over the life of the well, generating

58

a complex wastewater referred to as “flowback” in the initial stages of the well, and

59

produced water during the oil and gas production phase (collectively referenced herein as

60

“produced water”). Produced waters are a complex mixture of organic and inorganic

61

constituents, that originate either from the fracturing fluid (>1,000 chemicals potentially

62

used)

63

radioactive material [NORM]),13 or as a result of subsurface reactions between the

64

two.12,14 It has been estimated that over 3x1012 liters of produced water are generated

65

annually in the U.S. alone, creating a substantial waste stream for the oil and gas

66

industry.15

9,11,12

, the shale formation (i.e. hydrocarbons, salinity, metals, naturally occurring

67 68

There are environmental concerns related to produced waters and their ability to

69

reach the environment through well casing failures, surface spills, and inadequately

3

ACS Paragon Plus Environment

Page 4 of 41

Page 5 of 41

Environmental Science & Technology

70

treated produced water.16–18 These concerns also arise from the most common

71

management practice, disposal of water through deep-well injection (95% of wells),

72

which has recently been linked to seismic activity19 and the introduction of endocrine

73

disrupting chemicals to the surrounding environment.20 The potential environmental

74

impacts of produced water can be countered by their potential benefits, particularly for

75

water stressed regions, if treated and managed properly (e.g. irrigation of crops).

76

Therefore, the need to improve our understanding of produced water constituents,

77

specifically their organic composition, for environmental release scenarios or reuse

78

situations is vitally important.

79 80

Numerous publications have presented data on the inorganic constituents found in

81

produced water samples,21–25 while a paucity of literature exists on the organic

82

constituents.12,26,27 Review articles have identified and prioritized organic chemicals that

83

have been reported in fracturing fluid mixtures9,11,28 and that could return in produced

84

water.5,7,12 However, there are only a few studies to date that target organic chemicals in

85

produced water samples, with several focused solely on non-polar compounds (i.e.

86

hydrocarbons),14,29 and the few studies that do include the polar fraction have analyzed a

87

limited number of samples.3,30,31 Furthermore, the majority of this research has been

88

conducted on dry formations (e.g. gas producing), with little research as whole on wet

89

formations (e.g. oil and gas producing).27

90 91

The objectives of this study were, first, to characterize and perform a mass

92

balance of the organic constituents found in the fracturing fluid and the returned

4

ACS Paragon Plus Environment

Environmental Science & Technology

93

produced water over a one-year period from a single well using non-specific organic

94

analyses (DOC, hydrophobicity, and fluorescence). The second objective was to identify

95

and to follow over time specific volatile and nonvolatile organic compounds. Thirdly, and

96

finally, we assessed the possible environmental hazards and removal mechanisms of

97

identified organic compounds during different stages of production. As such, this work

98

presents a unique comprehensive long-term study of the presence of organic compounds

99

used and returned in a hydraulically fractured well.

100 101

2. Methods

102

2.1. Well Characteristics and Water Samples

103

The focus of our study was a well located in Weld County, CO, which produces

104

both oil and gas from the Niobrara Formation at a depth of around 2,000 meters. The well

105

was completed using a gel-based hydraulic fracturing fluid with around 11,000,000 L of

106

water. The fracturing fluid included 10 different proprietary chemical mixtures for a total

107

of around 87,000 liters/Kg of chemicals injected (SI Table 1; FracFocus and operator’s

108

well report). Water volume returned was from the separator’s water meter, which data

109

logs the flow of water over time.

110 111

The water samples collected from this individual well included the water used to

112

make the fracturing fluid (groundwater), the pre-injected fracturing fluid prior to

113

injection, and nine produced water samples over time. The first produced water sample

114

was taken during the first day (day 1) of the well flowing back, after a 30-day shut-in

115

period. Samples were collected on days 4, 7, 15, 22, 80, 130, 220, and 405 after flowing

5

ACS Paragon Plus Environment

Page 6 of 41

Page 7 of 41

Environmental Science & Technology

116

back the well. The samples taken on days 1 and 4 were drawn from a tank connected

117

directly to the well head (emptied every 4 to 6 hours), and all other produced water

118

samples were taken directly from the onsite gas/oil/water separator. All samples, except

119

for the fracturing fluid (wide mouth glass jar) were collected in pre-combusted (550°C

120

for 4 hr) 1-L amber bottles with no headspace. Samples were immediately transported

121

back to the lab and stored at 4°C prior to analysis (< 5 days). Further details about the

122

well, sampling procedure, and fracturing fluid can be found in Rosenblum et al. (2017).32

123 124

2.2. Dissolved Organic Matter Characterization, Fractionation, and Mass Balance

125

The characterization and fractionation of dissolved organic carbon (DOC) was

126

conducted at the U.S. Geological Survey laboratory (Boulder, CO). The fractionation was

127

done using XAD resins, according to Aiken et al. (1992),33 and included HPI as the

128

hydrophilic fraction (not retained by XAD resins), TPIA as transphilic acids (retained by

129

XAD-4), and HPOA as hydrophobic organic acids (retained on XAD-8). DOC analysis

130

(non-purgeable organic carbon) of unfractionated water and XAD fractions was

131

conducted using a TOC-Vcph analyzer (Shimadzu Corp., Japan). In addition to these

132

analyses, ultraviolet light absorbance was measured using an Agilent 8453 photodiode

133

array UV-Vis spectrophotometer, with specific ultraviolet absorbance (SUVA254)

134

calculated (UV absorbance in cm-1 divided by TOC in mg/L), and a spectrofluorometer

135

(Horiba FluoroMax-3, NJ; Raman units [RU]) was used to generate an excitation-

136

emission matrix; all samples were filtered through a 0.45-µm prior to analysis.34

137

6

ACS Paragon Plus Environment

Environmental Science & Technology

138

The mass balance was performed, and normalized to the mg of carbon in the

139

fracturing fluid, using Equation 1, where V is volume (L) and C is DOC concentration

140

(mg/L). The index corresponds to the order the samples were taken (i.e., i=2, day 4;

141

i=3…i=10, day 405). Values used and calculated are presented in SI Table 2.  :          =

    

  



"

$    

+   



#

+

∑+ , " &, − , )  $    

 +  ) * " 

142 143

2.3. Anion Determination

144

The Laboratory for Environmental and Geological Studies (LEGS) at the

145

University of Colorado Boulder performed the chloride analyses using ion

146

chromatography (Dionex Series 4500I). Samples were diluted by x 10, x 100, and x 1000,

147

with the average value being calculated from these three dilutions, if the levels were

148

above the analytical method detection limit.30

149 150

2.4. Volatile and Polar Organic Compound Analysis

151

Volatile organic compounds (VOCs) were determined from headspace-free 40-

152

mL volatile organic analysis (VOA) vials. VOCs were determined by purge and trap (OI

153

Analytical- 4660) gas chromatography mass spectrometry Agilent 6890 and 7890 GC-

154

MS) by a NELAC/A2LA accredited lab using EPA method 8260C (Origins Laboratory

155

Inc., Denver CO). Polar organic compounds were determined using ultrahigh pressure

156

liquid chromatography (1290 Agilent UHPLC; Agilent Technologies, CA) with an

157

ultrahigh definition quadrupole time-of-flight mass spectrometer (qTOF-HRMS; 6540

158

Agilent qTOF; Agilent Technologies, CA), with a 2-5% error of signal intensity for direct

7

ACS Paragon Plus Environment

Page 8 of 41

Page 9 of 41

Environmental Science & Technology

159

injection. All samples were filtered through a 0.2-µm filter (Acrodisc), and 10 µL were

160

injected onto a reverse phase C8 analytical column (Zorbax Eclipse XDB-C8). The

161

mobile phases were H2O with 0.1% formic acid and acetonitrile. The chromatographic

162

method was: 5 minutes (10% H2O with 0.1% formic acid/90% acetonitrile), followed by

163

linear gradient to 100% acetonitrile after 30 minutes, with a flow rate of 0.6 mL/min. Ions

164

were generated using electrospray Jet Stream Technology, operating in positive ion

165

mode. The data were processed using MassHunter software (version 6.1). This method

166

and details about the filters used during sample preparation are further described in

167

Thurman et al. (2014)3. Triisopropanolamine (98%) was purchased from Acros Organics.

168 169

3. Results and Discussion

170

3.1. Dissolved Organic Carbon (DOC) and Mass Balance

171

The

non-specific

organic

characterization

(DOC,

hydrophobicity,

and

172

fluorescence) of the groundwater, fracturing fluid, and produced water samples are

173

presented in this section. Figure 1A shows the concentrations of DOC and chloride from

174

the beginning to the end of the study. DOC concentrations ranged from ~1500 mg/L in

175

day 1 flowback water to ~200 mg/L in the final produced water (SI Table 3), which is

176

approximately a ten-fold decrease in total DOC. The pre-injected fracturing fluid had a

177

DOC of ~2200 mg/L, while the groundwater used to mix the fracturing fluid was ~5

178

mg/L DOC (SI Table 3); thus, the addition of fracturing fluid additives increased the

179

DOC by several orders of magnitude. SI Table 1 shows the list of ingredients and their

180

volumes from the FracFocus report and Well Report for this well, including 37,000 L of

181

guar slurry, 11,000 L of non-emulsifiers (proprietary surfactants), and additional organic

8

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 41

182

ingredients that likely contributed to the measured DOC. Chloride concentrations rapidly

183

increased from day 1 to day 22, and then appeared to reach a pseudo steady-state over the

184

remainder of the study (Figure 1A and SI Table 3).32 Figure 1B shows chloride versus

185

DOC concentrations over time.

186 187

B

2000

12000

1500

9000

ClDOC

1000

Cl- (mg/L)

DOC (mg/L)

A

e Tim

1 Day

0

0 0

188

6000 3000

500

400 Day

100

200

300

1500

400

1000

500

0

DOC (mg/L)

Time (Days)

189

Figure 1: Concentration of chloride and DOC in the produced water samples over time

190

(A), and Cl- versus DOC over time (B). Note, the y-axis for chloride represents both plots

191 192

The observed decline in DOC concentrations with time is consistent with other

193

studies that investigated DOC in wells over time.35–37 Hayes, who investigated several

194

wells with 3 to 5 samples over 90 days, had similar findings regarding DOC over time for

195

the majority of the wells studied,35 whereas Cluff et al. had nearly the same trend for the

196

three wells they investigated.36 Kim et al., who investigated two Denver-Julesburg Basin

197

(DJ-Basin) samples over 200 days, showed similar total organic carbon (TOC) results to

198

those seen here; although they had roughly a 50% decrease in TOC over the 200 days

199

compared to our 70% decrease of DOC over 220 days.37 Furthermore, research that has

9

ACS Paragon Plus Environment

Page 11 of 41

Environmental Science & Technology

200

denoted the age of a produced water (time since flowback) have generally shown

201

declining levels of DOC, chemical oxygen demand (COD), or TOC over time.38–40

202 203

A mass balance was performed for the DOC measured in the produced water

204

samples to determine the percentage of DOC returning over the course of the study

205

(Figure 2). Approximately 10.5% of the injected DOC returned by day 220, and only an

206

additional 0.6% returned between days 220 and 405 (11.1%). This 11.1% represents a

207

maximum fraction returned, because the formation itself contributes DOC, which is not

208

discernable from the fracturing fluids DOC. Therefore, these results indicate that a small

209

percentage of DOC from the injected fracturing fluids returns to the surface over 405

210

days of operation (< 11.1%). This 11.1% does not consider the volatile organic

211

compounds (VOCs), because they are removed during non-purgeable organic carbon

212

analysis, so the TOC recovered is likely slightly higher than the 11.1% of the DOC

213

observed. This low level of returning DOC is not seen with the returning water volumes,

214

which had around 34% of the volume injected return by day 405 (Figure 2); more than

215

three times that of the DOC.

216 217

These DOC results could be solely from dilution and mixing with subsurface

218

brine, however it is likely the returning fluid is influenced by more than just dilution and

219

mixing, which is often thought to be a key contributor to the total dissolved solids (TDS)

220

of returning produced waters.32 A combination of effects, such as sorption,

221

biodegradation, subsurface reactions, or imbibition are likely playing a role. Sorption is

222

likely due to the organic-rich mud, shale, and hydrocarbons present within the borehole

10

ACS Paragon Plus Environment

Environmental Science & Technology

223

that could sequester DOC from the produced water. Biodegradation is also possible since

224

there were high levels of DOC that could be labile (i.e. guar).41 A significant diversity of

225

microorganisms have been shown to return in produced water,36,42 even in the presence of

226

biocides.43 Imbibition in particular is important to consider, even for oil and gas bearing

227

formations,44 since it has been proposed that the majority of injected fluids are retained

228

(e.g. imbibed) in shale formations and that the returning waters are primarily natural

229

occurring formation brine.45 Thus, the 34% of injected water that returned (Table SI 2;

230

4% increase from Day 220 to Day 405), which again represents a maximum fraction

231

returned because much of the returning volume is likely influenced by the formation

232

itself, suggests that a substantial volume injected was likely imbibed in the formation.

233

Furthermore, Kim et al (2016)37 performed a mass balance of returning DJ-Basin

234

produced water over time, and saw significant variation in their frac fluid additive

235

recoveries, with recoveries of aluminum at 3%, zirconium at 9%, and potassium at

236

33%.37 Therefore, the data and literature suggest that imbibition, along with

237

mixing/dilution, could both be playing significant roles, while sorption and

238

biodegradation could also be contributing, and that it is additive-specific. Furthermore,

239

this 11.1% returned DOC could be composed of injected, formation, and subsurface

240

reaction based DOC, which make up the ~2,300 kg of DOC returned by day 405,

241

including chemicals of environmental concern from either the fracturing fluid (i.e.

242

biocides) or formation (i.e. benzene), validating the need to characterize and identify this

243

DOC.

11

ACS Paragon Plus Environment

Page 12 of 41

Environmental Science & Technology

% DOC and Water Returned

Page 13 of 41

100

DOC

80

Water

60 40 20 0 0

244

100

200

300

400

Time (Days)

245

Figure 2: Percent of DOC returned based on a mass balance of ~24,000 Kg of DOC

246

injected. Percent of water returned is also plotted and was calculated using the volume

247

returned on that day relative to the volume injected (11,000,000 L).

248 249

3.1.1. Dissolved Organic Carbon Fractionation and Characterization

250

The fractionation of produced waters DOC using XAD resins is a valuable

251

method for determining the distribution of organic matter class in oil and gas

252

wastewater.30 The fractionation results for the produced waters over time indicate the

253

majority of DOC for all of the samples and fracturing fluid was composed of low

254

molecular weight hydrophilic compounds (HPI; Figure 3a). This is in agreement with

255

Lester et al. (2015)30 who investigated a single DJ-Basin flowback sample.30

256

The day 1 sample had the highest percentage of hydrophilic acid (HPI; Figure

257

3A). The percentage of HPI steadily decreased over time, similar to that of the DOC

258

(Figure 3A & B). HPI has been linked to simple organic acids (e.g. acetic acid), which

259

can be found at high levels in produced waters.30 HPI and hydrophilic neutral compounds

260

are also likely responsible for much of the polar organic constituents (i.e. guar), and they

261

have been reported to consist of proteins and carbohydrates.46,47 Carbohydrates are of

262

particular interest, since they are generally used as gelling agents (i.e., guar) during

12

ACS Paragon Plus Environment

Environmental Science & Technology

263

hydraulic fracturing and were used in the well, studied herein (SI Table 1). This result,

264

then suggests that the gelling agents, which would be found in the HPI fraction, could be

265

responsible for a significant portion of the DOC in these samples. This hypothesis is

266

supported by the concentration of guar slurry injected (SI Table 1) and the trends shown

267

in Figures 1 and 3B, which illustrate that DOC and HPI follow nearly identical trends

268

over time. However, it should be noted that the HPI could also consist of acetone, acetic

269

acid, and other hydrophilic DOC compounds that are known to be present in produced

270

waters, either from the formation or the injected fluid.29,30

271 272

The other fractions, transphilic acids (TPIA) and hydrophobic organic acids

273

(HPOA), increase relative to the total measured DOC over time (Figure 1 and Figure 3A).

274

TPIA is thought to consist of intermediately polar compounds and has been linked to

275

aromatic constituents with acidic functional groups, which have been found in produced

276

water.14,48 HPOA is generally thought to consist of the less polar compounds and has

277

been associated with humic material, along with phenols and cresols,49 which have been

278

measured in DJ-Basin produced waters.30 Since these fractions increase relative to the

279

total DOC over time (Figure 3A), it is likely that these fractions are mostly from the

280

formation water and not the injected fluid. However, it is important to note that these

281

fractions could also be forming in the subsurface from the degradation of injected

282

compounds. This is illustrated in Figure 3B, which shows TPIA and HPOA

283

concentrations decrease by 43% and 48%, respectively, while the HPI fraction decreased

284

by 89%.

285

13

ACS Paragon Plus Environment

Page 14 of 41

Page 15 of 41

Environmental Science & Technology

B HPI TPIA

80

HPOA

60 40 20

80

1000

60 40

500

20 0

0 100

200

300

400

Time (Days)

Time (Days)

286

1500

HPI TPIA HPOA

0

80 13 0 22 0 40 5

22

7

15

4

1

0

100

HPI DOC mg/L

Percent of Total DOC

100

TPIA or HPOA DOC (mg/L)

A

287

Figure 3: DOC fractionation of produced waters over time by XAD resins 4 and 8 (HPI,

288

hydrophilic fraction; TPIA, transphilic acids; HPOA, hydrophobic organic acids).

289 290

UV absorbance at 254nm (SUVA254) and fluorescence were measured on all of

291

the samples.50 Results are presented in SI Table 3 and SI Figure 1. The results for UV

292

absorbance at 254nm show decreasing absorbance over time, while SUVA254 increases

293

slightly over the course of the study. As UV absorbance at 254nm is generally linked to

294

aromaticity, it is likely that SUVA254-rich formation-based aromatic constituents (e.g.,

295

benzene) not solely linked to the original fracturing fluid constituents were introduced.

296

That is, the UV 254nm values do not decrease proportionally to the DOC over time

297

(SUVA254), suggesting that a portion of the DOC that absorbs highly at 254nm could be

298

from the formation. Excitation emission matrices were generated from the fluorescence

299

data (SI Figure 1). The discussion of these results is presented in the SI, along with SI

300

Figure 1.

301 302

3.2. Characterization and environmental implications of Volatile Compounds

14

ACS Paragon Plus Environment

Environmental Science & Technology

303

The determination of specific volatile compounds and their environmental

304

implications are presented in this section. Sixty-five volatile organic compounds (VOCs)

305

were determined in all of the samples from the pre-injected fracturing fluid to produced

306

water recovered on day 405. Selected results are presented in Figure 4, which shows

307

VOC concentrations over time (concentrations of VOCs and their reporting limits are

308

presented in SI Table 4). The groundwater, which was stored on site, contained three

309

VOCs: 2-hexanone (312 µg/L), acetone (44.7 µg/L), and benzene (8.08 µg/L). Their

310

presence could be the result of chemicals (i.e. biocides, preservatives) added to the

311

storage container, produced water previously stored in the tank, the groundwater, or from

312

these volatile chemicals being transported (i.e., air) to the open storage tank during

313

drilling operations.

314 315

The pre-injected fracturing fluid contained several VOCs including 2-butanone

316

(10.7 mg/L), acetone (230 mg/L), bromomethane (0.07 mg/L), and chloromethane (0.03

317

mg/L). Acetone has been detected in flowback and produced water, as well as ambient air

318

around oil and gas wells, and is thought to be a common solvent used in cleaning

319

purposes.30,35,51 2-butanone has also been previously reported in produced water, and is

320

thought to be used as a friction reducer or industrial solvent.35 Chloromethane has been

321

described as a fracturing fluid additive, but has yet to be detected in produced water.9,52

322

Neither chloro- or bromomethane were detected in any produced water samples,

323

suggesting that they either transform downhole or sorb to the formation and particles in

324

the borehole.14

325

15

ACS Paragon Plus Environment

Page 16 of 41

Page 17 of 41

326

Environmental Science & Technology

3.2.1. BTEX and Substituted Benzenes

327

Benzene, toluene, ethylbenzene, and xylene (BTEX) were found in the mg/L

328

range in every produced water sample analyzed (Figure 4A), and are frequently detected

329

in produced waters and oil and gas impacted waters.14,29,30,35 Ethylbenzene and xylene

330

concentrations appear to decrease overtime, after day 22, while toluene decreases until

331

day 220 then increases on day 405. Benzene concentrations varied over the course of the

332

study with no particular trend. Benzene concentrations ranged from 6-12 mg/L, with the

333

lowest concentration being detected on day 220, while the second highest concentration

334

was detected on day 405 (Figure 4A). The high concentrations of ethylbenzene and

335

xylene in the early samples that peaked at day 22, suggest that BTEX compounds could

336

have been constituents of fracturing fluid additives (e.g. components of petroleum

337

distillates), even though they were not detected in the pre-injected fracturing fluid

338

sample. Nonetheless, the high concentrations observed for BTEX over the course of the

339

study suggest that substantial amounts are introduced from the formation, probably from

340

the oil in this formation, which would be anticipated for water in contact with an oil-

341

bearing shale. These concentrations are an order of magnitude lower than that recently

342

found in the Permian Basin produced waters (average concentration = 107 mg/L, n=8)29,

343

whereas other oil-bearing formations have shown high variations in reported benzene

344

levels from 0.026 to 4 mg/L.53–56

345 346

Benzene concentrations measured for all of the produced waters in the time series

347

studied are of environmental concern in the event of a spill/release, or if the produced

348

waters are to be used for applications other than recycling for future hydraulic

16

ACS Paragon Plus Environment

Environmental Science & Technology

349

fracturing.29,57,58 These environmental concerns are the result of benzene concentrations

350

being orders of magnitude higher than both the maximum contaminant level (MCL;

351

0.005 mg/L) and discharge permit levels used in produced water treatment facilities

352

(Pennsylvania [PA], U.S.).17 Toluene was also present in concentrations greater than its

353

MCL (1.0 mg/L) and PA discharge permits for all of the produced waters in the time

354

series studied.17 Ethylbenzene was greater than its MCL (0.7 mg/L) for days 7, 22, and

355

80, while xylenes were never greater than its MCL (10 mg/L total xylenes).

356

17

ACS Paragon Plus Environment

Page 18 of 41

Page 19 of 41

Environmental Science & Technology

Xylenes, total Toluene Ethylbenzene Benzene

A 15

mg/L

10

5

0 0

100

200

300

400

Time (Days)

B 1,2,4-Trimethylbenzene 1,3,5-Trimethylbenzene Isopropylbenzene sec-Butylbenzene

2.5

mg/L

2.0 1.5 1.0 0.5 0.0 0

100

200

300

400

Time (Days)

C 2.5

2-Hexanone 2-Butanone 4-Methyl-2-pentanone

mg/L

2.0 1.5 1.0 0.5 0.0 0

100

200

300

400

Time (Days)

357 358

Figure 4: Concentrations of BTEX (A), select methylbenzenes (B), and select ketones (C)

359

in produced water samples over time (day 1 to day 405).

360 361

Substituted benzenes were found at or near mg/L concentrations during the early

362

stages of the well (day 1 to day 80), and were above the limit of detection until day 220 18

ACS Paragon Plus Environment

Environmental Science & Technology

363

(Figure 4B). These four substituted benzenes all have the same trend, with concentrations

364

increasing until day 22 and then decreasing until day 200 and being below the detection

365

limit by day 405, suggesting they are likely associated with fracturing fluid additives.

366

1,2,4-trimethylbenzene, 1,3,5-trimethylbenzene, and isopropylbenzene have been

367

described as fracturing fluid additives, and have also been detected at the µg/L range in

368

produced water from both the Barnett and Marcellus shale.9,35,59 1,2,4-trimethylbenzene

369

has been described as a surfactant, corrosion inhibitor, and friction reducer and has also

370

been reported in surface water near treatment facilities that receive oil and gas fluids in

371

Pennsylvania.17 Both 1,2,4-trimethylbenzene and 1,3,5-trimethylbenzene were also

372

recently reported in monitoring wells near a hydraulically fractured well in Wyoming.60

373

Substituted benzenes have been shown to be irritants (i.e. eye, skin, and respiratory)61,

374

although no MCLs or oral reference doses exist for these four compounds; however a

375

recent study did quantify their potential toxicity using a quantitative structure-activity

376

relationship model.62 Therefore, the environmental impacts associated with substituted

377

benzenes are not well understood and require further research.

378 379

3.2.2. Ketones

380

Numerous ketones were at the mg/L level for portions of the study. Acetone was

381

the ketone detected at the highest concentration on day one (22 mg/L) and was detected

382

over the course of the study with concentrations ranging from 0.6 mg/L on day 7 to 9

383

mg/L on day 405. The levels of acetone over the course of the study are likely due to its

384

high concentration in the fracturing fluid (230 mg/L), meaning it would take a substantial

385

amount of time to dilute (i.e., formation water), sorb, or biodegrade acetone from the

19

ACS Paragon Plus Environment

Page 20 of 41

Page 21 of 41

Environmental Science & Technology

386

returning produced water. Also the polar nature of acetone would support its presence in

387

the water fraction for an extended period of time, whereas non-polar or mid-polar

388

compounds would partition into the oil. Acetone has been detected in numerous produced

389

waters, been found to form during biodegradation of fracturing fluid additives, and has

390

been reported as a fracturing fluid additive.9,30,35,41,52 Acetone has been described as “less

391

toxic” than other industrial solvents,61 and although it does not have an MCL, it has an

392

oral reference dose (RfD) of 0.9 mg/kg/day, which is more than three times that of the

393

total xylenes RfD (0.2 mg/kg/day). Total xylenes has an MCL of 10 mg/L.

394 395

The other ketones detected were 2-hexanone, 2-butanone, and 4-methyl-2-

396

pentanone (Figure 4C). 2-Hexanone was detected at its highest concentration on day 1 (1

397

mg/L) and was last detected on day 130 (0.1 mg/L), suggesting that it could have been

398

formed in the subsurface or from an unknown source. 2-hexanone has not previously

399

been reported in produced water or as a chemical additive used in hydraulic fracturing,

400

thus supporting the notion that it was formed during a subsurface reaction, potentially due

401

to the use of tert-butyl peroxide, which is a strong oxidant. 2-hexanone was historically

402

used in paint, paint thinner, and various chemical processes.63 However, because 2-

403

hexanone was determined to have negative health effects as a neurotoxin, it now has

404

restricted uses in the U.S.63 It has an RfD of 0.005 mg/kg/day, which is the same RfD

405

value as benzene and dichloromethane.64 Thus, 2-hexanone is a chemical of concern

406

because it was found in the produced water samples for at least 130 days.

407

20

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 41

408

2-butanone was detected at 10 mg/L in the fracturing fluid, but was not detected

409

until day 80 in the produced waters samples. This is different than acetone, which was

410

detected at high levels in the fracturing fluid and was also detected in the day 1 sample.

411

These unexpected results for 2-butanone could reflect subsurface reactions. That is, 2-

412

butanone could be yielding 2-hexanone in the early stages of the well, because 2-

413

butanone was detected in the fracturing fluid and appears to increase over time while 2-

414

hexanone is decreasing (Figure 3c). However, the addition of an ethyl group to 2-

415

butanone seems unlikely, even in the presence of a strong oxidant (i.e. tert-butyl

416

peroxide). Nonetheless, 2-butanone appears to be a fracturing fluid additive, even with

417

its unanticipated levels over the course of the study because it was detected in the

418

fracturing fluid, reported as a fracturing fluid additive,9 and has also been detected in

419

produced water and waters impacted by produced water.30,60

420

pentanone was at or near mg/L levels from day 1 until day 22, but was not detected by

421

day 80. This suggests that it was likely associated with the fracturing fluid and not the

422

formation. 4-methyl-2-pentanone has been described as a fracturing fluid agent, and is a

423

common solvent used in industrial processes.9 2-butanone (RfD of 0.6 mg/kg/day) and 4-

424

methyl-2-pentanone (reference concentration for inhalation exposure of 3 mg/m3) have

425

been shown to effect developmental and musculoskeletal systems.64 This could support

426

monitoring for both chemicals in the event of a spill, particularly 4-methyl-2-pentanone

427

in the early stages (i.e. first few weeks) of the well and 2-butanone at the later stages (i.e.,

428

first year).

429 430

3.3. Semi-Polar and Non-Polar Organic Compounds

21

ACS Paragon Plus Environment

Finally, 4-methyl-2-

Page 23 of 41

Environmental Science & Technology

431

The analysis for specific polar organic compounds and their environmental

432

implications are presented in this section. The UHPLC-Q-TOF-HRMS total ion

433

chromatograms (TIC) for the pre-injected fracturing fluid and five produced water

434

samples are shown in Figure 5. The peaks visible in the TIC are dominated by three

435

classes of homologues, which have been previously identified in produced waters (Figure

436

5): polyethylene glycols (PEGs)3, polypropylene glycols (PPGs)65, and linear alkyl

437

ethoxylates (LAEs).3 These three homologous compounds were identified in the

438

fracturing fluid and the produced water samples, and generally decreased in magnitude as

439

the produced waters age increased (Figure 5). During hydraulic fracturing, surfactants are

440

added to aide in the pumping of fracturing fluids and the movement of sand and oily

441

substances.30 The surfactants are thought to be potential fingerprinting compounds, since

442

they are identified so frequently in produced waters and waters impacted by them.3,17,65,66

443

Their substantial signal in both the fracturing fluid and produced waters, along with their

444

presence over the course of the study supports this notion, specifically for PEGs and

445

PPGs. It should also be noted that PEGs were detected in the groundwater samples from

446

3.5 to 11.5 min (SI Figure 2), and could have been the result of biocide addition in the

447

tank, since they were added to prevent microbial growth during storage; although they

448

could have also been from flowback and produced water that were mixed or previously

449

stored in the holding tank.

450

22

ACS Paragon Plus Environment

Environmental Science & Technology

451 452

Figure 5: Total ion chromatogram (TIC) of five produced water samples over 220 days

453

and the pre-injected fracturing fluid (insert above other TICs). Note the inserted

454

fracturing fluid TIC is compressed to fit the figure and is not to scale on the y-axis

455

(relative abundance). Polyethylene glycols (PEGs), polypropylene glycols (PPGs), and

456

linear alkyl ethoxylates (LAEs) are denoted in the figure, along with triisopropanolamine

457

(TIPA). The x- and y-axes represents the acquisition time and the relative size of the

458

peaks (counts), respectively. Samples were analyzed as received (< 48 hours of

459

sampling). It is important to note, that there are numerous other compounds present in

460

these samples and in the regions identified (as PEGs, PPGs, LAEs, and TIPA), and that

461

these identifications are shown to illustrate the major ions present and not the total

462

compounds present.

463 464

The PEGs were detected from 4.4 min to 14.0 min (more hydrophilic zone). The

465

PEGs showed increasing ethylene oxide (EO) units, from PEG-EO4 to PEG-EO16 with

23

ACS Paragon Plus Environment

Page 24 of 41

Page 25 of 41

Environmental Science & Technology

466

accurate masses of 217.1050 m/z (Na+ adduct; retention time 4.4 min; 2.0 ppm error) and

467

740.4645 m/z (NH4+ adduct; retention time 12.4 min; 1.0 ppm error), respectively

468

(sample, day 1). The average mass difference between PEGs was 44.0262 mass units,

469

representing the addition of one ethylene oxide group [-CH2-CH2-O-]. The next class of

470

compounds was PPGs, which were detected from 14.6 min to 22.2 min and show

471

increasing propylene oxide units from PPG-PO4 to PPG-PO12 with accurate masses of

472

273.1673 m/z (Na+ adduct; retention time 14.6 min; 0.2 ppm error) and 732.5470 m/z

473

(NH4+ adduct; retention time 22.2 min; 0.5 ppm error), respectively (sample, day 1). The

474

PPGs had an average mass difference of 58.0419 mass units that is equal to a propylene

475

oxide group [-CH2CH(CH3)-O-].

476 477

The LAEs were detected from 22 min to 26.8 min and show four different classes

478

of LAEs that have varying carbon chain lengths (C9, C10, C11, and C12). This work

479

describes the first putative identification of LAEs in produced water with carbon chain

480

lengths of 9, 10, and 11 (SI Table 5), whereas the 12 carbon chain length LAEs have

481

previously been described.3 The LAEs had an average mass difference of 44.0262 mass

482

units that is equal to an ethylene oxide group [-CH2-CH2-O-]. The LAEs included over 25

483

unique putatively identified compounds, from LAE9-EO13 to LAE12-EO9, which had

484

accurate masses of 734.5266 m/z (NH4 adduct; retention time 22 min; 0.8 ppm error) and

485

600.4688 m/z (NH4 adduct; retention time 26.8; 1.2 ppm error), respectively (samples,

486

day 1 and day 15). It should be noted that individual standards are not available for PEGs,

487

PPGs, and LAEs, with only homologue mixtures available at this time, with these

488

compounds being identified and semi-quantitated.

24

ACS Paragon Plus Environment

Environmental Science & Technology

489 490

PEGs, PPGs, and LAEs are common products in personal care products,

491

medicines, foods, and industrial processes.65 For instance, PEGs are used in laxatives,

492

PPGs are used in ice cream, and LAEs are often found in detergents.65 Their use in

493

commercial products that humans frequently ingest or are in contact with suggests they

494

are not acutely toxic; however, the extent of polymerization is crucial when considering

495

toxicity.65 PEGs of high molecular weight have been described as inert and non-toxic,67

496

while the monomer ethylene glycol has known human toxicity.68 PPGs have been shown

497

to have a similar relationship to that of PEGs, with polymers or high molecular weight

498

PPGs being considered safe,69 while the monomer propylene glycol is considered toxic.70

499

Similarly, LAEs are considered harmless as a polymer.71,72 However, the monomer

500

ethylene glycol can be formed, because the LAEs studied here can have up to 13

501

repeating ethylene oxide groups that are attached to a carbon chain (SI Table 5).68

502

Therefore, consideration of whether the monomer is being formed is a key research

503

question, especially considering the use of strong oxidants as breakers, high temperatures

504

and pressures in the well, and biological processes in the subsurface that could lead to

505

monomer formation. It should be noted that the monomers discussed here are not

506

detectable in our UHPLC-Q-TOF-HRMS, and should be investigated in future research

507

of produced waters.

508 509

The compound triisopropanolamine (TIPA) was found in every sample from the

510

fracturing fluid to day 405. The peak can be seen in Figure 5 around 2.8 min. TIPA was

511

identified by its accurate mass of 192.1594 m/z, corresponding to the protonated molecule

25

ACS Paragon Plus Environment

Page 26 of 41

Page 27 of 41

Environmental Science & Technology

512

and elemental formula C9H21NO3. TIPA was further confirmed by retention time,

513

accurate mass, and by performing MS-MS on a pure standard that validated the putative

514

identification, with major fragments at 174.1494 m/z and 156.1383 m/z, (SI Figure 3).

515

The resulting concentration appears to be in the mg/L range for all of the samples

516

analyzed, with relative intensities comparable to the 1 mg/L standard.

517 518

TIPA is a tertiary amine, which is typically used as an emulsifier, stabilizer, and

519

surfactant during industrial processes. It is also frequently used in cements as a hardening

520

accelerator.73 TIPA has been reported as a fracturing fluid additive11,52 but had yet to be

521

reported in a produced water sample. TIPA is an industrial agent that has been described

522

as a moderate irritant for skin and a severe irritant for eyes,61 but has been shown to be

523

non-toxic or genotoxic in in-vitro bioassays.74

524 525

3.3.2. Surfactant Levels Over Time

526

Figure 5 and 6 illustrate the substantial changes in concentration of PEGs, PPGs,

527

LAEs, and TIPA in the produced water samples over time. All three surfactants

528

decreased in line with their hydrophobic nature (octanol-water partition coefficient; Kow),

529

i.e., the most hydrophobic LAEs decreased most rapidly followed by the PPGs, and then

530

PEGs.75 Because a quantitative measure would be challenging due to matrix-induced

531

ionization effects and specific response factors for these types of compounds,17 relative

532

abundance was used as a way to compare these compounds levels over time (i.e., semi-

533

quantitative).

534

26

ACS Paragon Plus Environment

Environmental Science & Technology

535 536

Figure 6: Relative abundance (concentration at time, t / original concentration; Ct/Co) of

537

the highest detected surfactants from each class (PPG-PO6, PEG-EO6, and LAE11-EO9),

538

and TIPA over the course of the study from the day 1 sample to day 405.

539 540

Hydrophobicity, as expected, is a key indicator for what injected chemicals are

541

present in the produced waters over time. This is demonstrated by the identified

542

compounds’ Kow values and their elution time through the UHPLC, because the least

543

hydrophobic compounds elute from the column first, while the most hydrophobic

544

compounds elute last. The PEGs are the least hydrophobic surfactants with an estimated

545

log Kow of -2.35 earliest elution time and the greatest relative abundance (Figure 5) over

546

the course of the study. The PPGs had the second highest concentrations among the

547

samples analyzed, have an estimated log Kow of -0.25, and eluted after the PEGs. In

548

contrast, LAEs, for which the concentrations decreased the fastest, have an estimated log

27

ACS Paragon Plus Environment

Page 28 of 41

Page 29 of 41

Environmental Science & Technology

549

Kow of 2.4,5 the highest of the surfactants studied, and also eluted the latest. These results

550

suggest the LAEs are likely partitioning into the oil phase or sorbing to the formation or

551

particles in the wellbore (LAE estimated log Koc= 1.45 [carbon-normalized sorption

552

coefficient]), which could explain the rapid removal from the water phase of the LAEs

553

compared to PPGs and PEGs. Conversely, the PEGs and PPGs, both with lower Kow

554

values, indicating that they will remain in the aqueous phase rather than partitioning into

555

the shale or oil phase, were detected in the produced water throughout the study period

556

(Figure 5). This hypothesis is further supported by TIPA, which remained at the highest

557

concentrations over time and had an estimated log Kow of -1.25 (Figure 5). Although its

558

estimated Kow is slightly higher than the PEGs, TIPA’s hydrophilic nature is illustrated

559

by its early elution time in the TIC (2.8 min) compared to PEG-EO6 (9.5 min) or PEG-

560

EO4 (4.4 min), which further supports the role hydrophobicity plays on returning

561

chemicals.

562 563

These results show how surfactants change over time, in terms of what class of

564

compounds might be present and how their concentrations change. The results also

565

support the importance of a chemical’s hydrophobic nature, and how that dictates when

566

the chemical will return in an oil-bearing formation. This is important when considering

567

the chemicals to target as tracers in the event of a spill or water reuse scenario,

568

particularly for an oil and gas containing formation. That is, the water and oil are

569

constantly mixing in the borehole until they reach the separator, providing substantial

570

time for compounds to partition from the water to the oil. These findings coupled with the

571

presence of PPGs, PEGs, and TIPA over the course of the study, illustrate three indicator

28

ACS Paragon Plus Environment

Environmental Science & Technology

572

chemicals that are unique to the water phase and that were found in the fracturing fluid

573

used (Figure 5). TIPA in particular could be of importance, due to its hydrophilicity, that

574

would result in TIPA being in produced waters for extended periods of time, along with

575

TIPAs non-ubiquitous nature compared to many surfactants (i.e. PEGs).

576 577

Overall, this study provides a unique organic characterization that fills a

578

substantial knowledge gap; that much of the injected DOC remains in the subsurface (on

579

a mass balance basis), and that the returning DOC is primarily hydrophilic based on the

580

return of the hydrophilic compounds, either volatile or non-volatile, over the 405-day

581

study period. Furthermore, this work newly-identified persistent compounds that remain

582

for long periods of time, and that have been found in produced waters around North

583

America or effluent discharges from a produced-water treatment facility.3,17,30,65,66 Future

584

research should continue to characterize the organic composition of produced waters (e.g.

585

in general and over time), with particular attention on hydrophilic organic compounds

586

and the suspended solids fraction, while work is also needed to characterize chemical or

587

biological degradation products of these specific constituents, and expound on the

588

physical mechanisms that impact the returning compounds.

589 590

Supplemental Information

591

Supplemental information includes tables and figures for the wells FracFocus report,

592

mass balance values, DOC and non-specific characterization, VOCs analyzed and

593

detected, Kendrick masses for LAEs, excitation-emission matrixes, the total ion

594

chromatogram of the groundwater used to make the fracturing fluid, and MS-MS spectra

29

ACS Paragon Plus Environment

Page 30 of 41

Page 31 of 41

Environmental Science & Technology

595

for Triisopropanolamine. An Excel file includes all data generated by the USGS

596

presented in SI Table 3 and Figure 3. At the time of publication, data other than

597

presented in this paper are not publically available from the University of Colorado.

598 599

Acknowledgements

600

This work was supported by the AirWaterGas Sustainability Research Network, and was

601

funded by the National Science Foundation under Grant No. CBET-1240584. Opinions,

602

findings, recommendations, and conclusions conveyed in this paper are those of the

603

author(s) and do not necessarily reflect the views of the National Science Foundation.

604

Any use of trade, firm, or product names is for descriptive purposes only and does not

605

imply endorsement by the U.S. Government. The researchers would also like to

606

acknowledge Dr. Jessica Rodgers for her review of the manuscript, Dr. Avner Vengosh

607

for his insight on our data set, and our industry collaborators that permitted access and

608

took the time to help us collect the samples. The authors dedicate this paper to the

609

memory of George Aiken, colleague, mentor, and friend.

610 611 612 613

References:

614

2. Vidic, R. D., Brantley, S. L., Vandenbossche, J. M., Yoxtheimer, D. & Abad, J. D.

1. Annual Energy Outlook 2015 with projections to 2040. 154 (2015).

615

Impact of Shale Gas Development on Regional Water Quality. Science 340,

616

1235009 (2013).

30

ACS Paragon Plus Environment

Environmental Science & Technology

617

3. Thurman, E. M., Ferrer, I., Blotevogel, J. & Borch, T. Analysis of Hydraulic

618

Fracturing Flowback and Produced Waters Using Accurate Mass: Identification

619

of Ethoxylated Surfactants. Anal. Chem. 86, 9653–9661 (2014).

620

4. Adgate, J. L., Goldstein, B. D. & McKenzie, L. M. Potential Public Health Hazards,

621

Exposures and Health Effects from Unconventional Natural Gas Development.

622

Environ. Sci. Technol. 48, 8307–8320 (2014).

623

5. Rogers, J. D., Burke, T. L., Osborn, S. G. & Ryan, J. N. A Framework for Identifying

624

Organic Compounds of Concern in Hydraulic Fracturing Fluids Based on Their

625

Mobility and Persistence in Groundwater. Environ. Sci. Technol. Lett. 2, 158–164

626

(2015).

627

6. Yost, E. E., Stanek, J., DeWoskin, R. S. & Burgoon, L. D. Overview of Chronic Oral

628

Toxicity Values for Chemicals Present in Hydraulic Fracturing Fluids, Flowback,

629

and Produced Waters. Environ. Sci. Technol. 50, 4788–4797 (2016).

630

7. Stringfellow, W. T., Domen, J. K., Camarillo, M. K., Sandelin, W. L. & Borglin, S.

631

Physical, chemical, and biological characteristics of compounds used in

632

hydraulic fracturing. J. Hazard. Mater. 275, 37–54 (2014).

633

8. Torres, L., Yadav, O. P. & Khan, E. A review on risk assessment techniques for

634

hydraulic fracturing water and produced water management implemented in

635

onshore unconventional oil and gas production. Sci. Total Environ. 539, 478–493

636

(2016).

637

9. US EPA National Center for Environmental Assessment, I. O. Hydraulic

638

Fracturing for Oil and Gas: Impacts from the Hydraulic Fracturing Water Cycle

639

on Drinking Water Resources in the United States (Final Report). Available at:

31

ACS Paragon Plus Environment

Page 32 of 41

Page 33 of 41

Environmental Science & Technology

640

https://cfpub.epa.gov/ncea/hfstudy/recordisplay.cfm?deid=332990. (Accessed:

641

17th March 2017)

642

10. Gregory, K. B., Vidic, R. D. & Dzombak, D. A. Water Management Challenges

643

Associated with the Production of Shale Gas by Hydraulic Fracturing. Elements 7,

644

181–186 (2011).

645 646 647

11. Waxman, H., Markey, E. J. & Degette, D. Chemicals used in hydraulic fracturing. (2011). 12. Elsner, M. & Hoelzer, K. Quantitative Survey and Structural Classification of

648

Hydraulic Fracturing Chemicals Reported in Unconventional Gas Production.

649

Env. Sci Technol (2016). doi:10.1021/acs.est.5b02818

650

13. Orem, W. et al. Organic substances in produced and formation water from

651

unconventional natural gas extraction in coal and shale. Int. J. Coal Geol. 126, 20–

652

31 (2014).

653

14. Hoelzer, K. et al. Indications of Transformation Products from Hydraulic

654

Fracturing Additives in Shale-Gas Wastewater. Env. Sci Technol (2016).

655

doi:10.1021/acs.est.6b00430

656 657

15. Clark, C. E. & Veil, J. A. Produced Water Volumes and Management Practices in the United States. (Argonne National Laboratory (ANL), 2009).

658

16. Vengosh, A., Jackson, R. B., Warner, N., Darrah, T. H. & Kondash, A. A Critical

659

Review of the Risks to Water Resources from Unconventional Shale Gas

660

Development and Hydraulic Fracturing in the United States. Environ. Sci.

661

Technol. 48, 8334–8348 (2014).

32

ACS Paragon Plus Environment

Environmental Science & Technology

662

17. Getzinger, G. J. et al. Natural Gas Residual Fluids: Sources, Endpoints, and

663

Organic Chemical Composition after Centralized Waste Treatment in

664

Pennsylvania. Environ. Sci. Technol. 49, 8347–8355 (2015).

665

18. Osborn, S. G., Vengosh, A., Warner, N. R. & Jackson, R. B. Methane contamination

666

of drinking water accompanying gas-well drilling and hydraulic fracturing. Proc

667

Natl Acad Sci U A 108, 8172–6 (2011).

668

19. Ellsworth, W. L. Injection-Induced Earthquakes. Science 341, 142–+ (2013).

669

20. Kassotis, C. D. et al. Endocrine disrupting activities of surface water associated

670

with a West Virginia oil and gas industry wastewater disposal site. Sci. Total

671

Environ. 557–558, 901–910 (2016).

672

21. Warner, N. R. et al. Geochemical and isotopic variations in shallow groundwater

673

in areas of the Fayetteville Shale development, north-central Arkansas. Appl.

674

Geochem. 35, 207–220 (2013).

675

22. Barbot, E., Vidic, N. S., Gregory, K. B. & Vidic, R. D. Spatial and Temporal

676

Correlation of Water Quality Parameters of Produced Waters from Devonian-

677

Age Shale following Hydraulic Fracturing. Environ. Sci. Technol. 47, 2562–2569

678

(2013).

679

23. Haluszczak, L. O., Rose, A. W. & Kump, L. R. Geochemical evaluation of flowback

680

brine from Marcellus gas wells in Pennsylvania, USA. Appl. Geochem. 28, 55–61

681

(2013).

682

24. Balashov, V. N., Engelder, T., Gu, X., Fantle, M. S. & Brantley, S. L. A model

683

describing flowback chemistry changes with time after Marcellus Shale

684

hydraulic fracturing. Aapg Bull. 99, 143–154 (2015).

33

ACS Paragon Plus Environment

Page 34 of 41

Page 35 of 41

685

Environmental Science & Technology

25. Ziemkiewicz, P. F. & He, Y. T. Evolution of water chemistry during Marcellus

686

Shale gas development: A case study in West Virginia. Chemosphere 134, 224–

687

231 (2015).

688

26. Ferrer, I. & Thurman, E. M. Chemical constituents and analytical approaches for

689

hydraulic fracturing waters. Trends Environ. Anal. Chem. 5, 18–25 (2015).

690

27. Luek, J. L. & Gonsior, M. Organic compounds in hydraulic fracturing fluids and

691 692

wastewaters: A review. Water Res. 123, 536–548 (2017). 28. Elliott, E. G., Ettinger, A. S., Leaderer, B. P., Bracken, M. B. & Deziel, N. C. A

693

systematic evaluation of chemicals in hydraulic-fracturing fluids and wastewater

694

for reproductive and developmental toxicity. J. Expo. Sci. Environ. Epidemiol. 27,

695

90–99 (2017).

696 697 698

29. Khan, N. A. et al. Volatile-organic molecular characterization of shale-oil produced water from the Permian Basin. Chemosphere 148, 126–36 (2016). 30. Lester, Y. et al. Characterization of hydraulic fracturing flowback water in

699

Colorado: Implications for water treatment. Sci. Total Environ. 512, 637–644

700

(2015).

701 702 703

31. Ferrer, I. & Thurman, E. M. Analysis of hydraulic fracturing additives by LC/QTOF-MS. Anal Bioanal Chem 407, 6417–28 (2015). 32. Rosenblum, J. et al. Temporal characterization of flowback and produced water

704

quality from a hydraulically fractured oil and gas well. Sci. Total Environ. 596–

705

597, 369–377 (2017).

34

ACS Paragon Plus Environment

Environmental Science & Technology

706

33. Aiken, G. R., McKnight, D. M., Thorn, K. A. & Thurman, E. M. Isolation of

707

hydrophilic organic acids from water using nonionic macroporous resins. Org.

708

Geochem. 18, 567–573 (1992).

709

34. Aquatic Organic Matter Fluorescence. (Cambridge University Press, 2014).

710

35. Hayes, T. Sampling and analysis of water streams associated with the

711

development of Marcellus shale gas. Marcellus Shale Initiat. Publ. Database 10,

712

(2009).

713

36. Cluff, M. A., Hartsock, A., MacRae, J. D., Carter, K. & Mouser, P. J. Temporal

714

Changes in Microbial Ecology and Geochemistry in Produced Water from

715

Hydraulically Fractured Marcellus Shale Gas Wells. Environ. Sci. Technol. 48,

716

6508–6517 (2014).

717

37. Kim, S., Omur-Ozbek, P., Dhanasekar, A., Prior, A. & Carlson, K. Temporal analysis

718

of flowback and produced water composition from shale oil and gas operations:

719

Impact of frac fluid characteristics. J. Pet. Sci. Eng. 147, 202–210 (2016).

720

38. Esmaeilirad, N., Carlson, K. & Ozbek, P. O. Influence of softening sequencing on

721

electrocoagulation treatment of produced water. J. Hazard. Mater. 283, 721–729

722

(2015).

723

39. Rosenblum, J. S., Sitterley, K. A., Thurman, E. M., Ferrer, I. & Linden, K. G.

724

Hydraulic Fracturing Wastewater Treatment by Coagulation-Adsorption for

725

Removal of Organic Compounds and Turbidity. J. Environ. Chem. Eng. (2016).

726

doi:10.1016/j.jece.2016.03.013

35

ACS Paragon Plus Environment

Page 36 of 41

Page 37 of 41

727

Environmental Science & Technology

40. Riley, S. M., Oliveira, J. M. S., Regnery, J. & Cath, T. Y. Hybrid Membrane Bio-

728

systems for Sustainable Treatment of Oil & Gas Produced Water and Frac

729

Flowback. Sep. Purif. Technol. (2016). doi:10.1016/j.seppur.2016.07.008

730

41. Kekacs, D., Drollette, B. D., Brooker, M., Plata, D. L. & Mouser, P. J. Aerobic

731

biodegradation of organic compounds in hydraulic fracturing fluids.

732

Biodegradation 26, 271–287 (2015).

733

42. Murali Mohan, A. et al. Microbial community changes in hydraulic fracturing

734

fluids and produced water from shale gas extraction. Env. Sci Technol 47,

735

13141–50 (2013).

736

43. Vikram, A., Lipus, D. & Bibby, K. Produced Water Exposure Alters Bacterial

737

Response to Biocides. Environ. Sci. Technol. 48, 13001–13009 (2014).

738

44. Birdsell, D. T., Rajaram, H., Dempsey, D. & Viswanathan, H. S. Hydraulic

739

fracturing fluid migration in the subsurface: A review and expanded modeling

740

results. Water Resour. Res. 51, 7159–7188 (2015).

741

45. Kondash, A. J., Albright, E. & Vengosh, A. Quantity of flowback and produced

742

waters from unconventional oil and gas exploration. Sci. Total Environ. 574,

743

314–321 (2017).

744

46. Jiang, J., Zhao, Q., Wei, L. & Wang, K. Extracellular biological organic matters in

745

microbial fuel cell using sewage sludge as fuel. Water Res. 44, 2163–2170

746

(2010).

747

47. Wei, L.-L. et al. Fractional, biodegradable and spectral characteristics of

748

extracted and fractionated sludge extracellular polymeric substances. Water Res.

749

46, 4387–4396 (2012).

36

ACS Paragon Plus Environment

Environmental Science & Technology

750

48. Wang, D., Hu, Q., Li, M., Wang, C. & Ji, M. Evaluating the removal of organic

751

fraction of commingled chemical industrial wastewater by activated sludge

752

process augmented with powdered activated carbon. Arab. J. Chem. 9, S1951–

753

S1961 (2016).

754

49. Croue, J. P., Benedetti, M. F., Violleau, D. & Leenheer, J. A. Characterization and

755

copper binding of humic and nonhumic organic matter isolated from the South

756

Platte River: Evidence for the presence of nitrogenous binding site. Environ. Sci.

757

Technol. 37, 328–336 (2003).

758

50. Weishaar, J. L. et al. Evaluation of specific ultraviolet absorbance as an indicator

759

of the chemical composition and reactivity of dissolved organic carbon. Environ.

760

Sci. Technol. 37, 4702–4708 (2003).

761

51. Colborn, T., Schultz, K., Herrick, L. & Kwiatkowski, C. An Exploratory Study of Air

762

Quality Near Natural Gas Operations. Hum. Ecol. Risk Assess. Int. J. 20, 86–105

763

(2013).

764 765

52. FracFocus. Groundwater Protection Council. Interstate Oil & Gas Conservation Commission. FracFocus Chemical Disclosure Registry.

766

53. Utvik, T. I. R. Chemical characterisation of produced water from four offshore oil

767

production platforms in the North Sea. Chemosphere 39, 2593–2606 (1999).

768

54. Horner, J. E., Castle, J. W. & Rodgers, J. H. A risk assessment approach to

769

identifying constituents in oilfield produced water for treatment prior to

770

beneficial use. Ecotoxicol. Environ. Saf. 74, 989–999 (2011).

37

ACS Paragon Plus Environment

Page 38 of 41

Page 39 of 41

Environmental Science & Technology

771

55. Sirivedhin, T. & Dallbauman, L. Organic matrix in produced water from the

772

Osage-Skiatook petroleum environmental research site, Osage county,

773

Oklahoma. Chemosphere 57, 463–9 (2004).

774 775

56. Dorea, H. S. et al. Analysis of BTEX, PAHs and metals in the oilfield produced water in the State of Sergipe, Brazil. Microchem. J. 85, 234–238 (2007).

776

57. Veil, J. A., Puder, M. G., Elcock, D. & Redweik, J., Robert J. A White Paper Describing

777

Produced Water from Production of Crude Oil, Natural Gas, and Coal Bed Methane.

778

(2004).

779

58. Xu, P., Drewes, J. E. & Heil, D. Beneficial use of co-produced water through

780

membrane treatment: technical-economic assessment. Desalination 225, 139–

781

155 (2008).

782 783

59. Hayes, T. & Severin, B. Characterization of flowback waters from the Marcellus and the Barnett shale regions. (2012).

784

60. DiGiulio, D. C. & Jackson, R. B. Impact to Underground Sources of Drinking Water

785

and Domestic Wells from Production Well Stimulation and Completion Practices

786

in the Pavillion, Wyoming, Field. Env. Sci Technol (2016).

787

doi:10.1021/acs.est.5b04970

788 789

61. TOXNET. Available at: https://toxnet.nlm.nih.gov/cgibin/sis/search2/f?./temp/~K2N5Lu:3. (Accessed: 17th March 2017)

790

62. Yost, E. E., Stanek, J., DeWoskin, R. S. & Burgoon, L. D. Estimating the Potential

791

Toxicity of Chemicals Associated with Hydraulic Fracturing Operations Using

792

Quantitative Structure–Activity Relationship Modeling. Environ. Sci. Technol. 50,

793

7732–7742 (2016).

38

ACS Paragon Plus Environment

Environmental Science & Technology

794

63. Public Health Statement: 2-Hexanone. (1992).

795

64. US EPA, O. Integrated Risk Information System. Available at:

796

https://www.epa.gov/iris. (Accessed: 23rd March 2017)

797

65. Thurman, E. M., Ferrer, I., Rosenblum, J., Linden, K. & Ryan, J. N. Identification of

798

polypropylene glycols and polyethylene glycol carboxylates in flowback and

799

produced water from hydraulic fracturing. J. Hazard. Mater. (2017).

800

doi:10.1016/j.jhazmat.2016.02.041

801 802 803

66. He, Y. et al. Chemical and toxicological characterizations of hydraulic fracturing flowback and produced water. Water Res. 114, 78–87 (2017). 67. Chiang, E. T. et al. Protective effects of high-molecular weight Polyethylene

804

Glycol (PEG) in human lung endothelial cell barrier regulation: Role of actin

805

cytoskeletal rearrangement. Microvasc. Res. 77, 174–186 (2009).

806

68. Kassotis, C. D., Tillitt, D. E., Davis, J. W., Hormann, A. M. & Nagel, S. C. Estrogen

807

and androgen receptor activities of hydraulic fracturing chemicals and surface

808

and ground water in a drilling-dense region. Endocrinology 155, 897–907

809

(2014).

810 811 812 813

69. Lanigan, R. S. Final report on the safety assessment of PPG-11 and PPG-15 Stearyl Ethers. Int. J. Toxicol. 20, 53–59 (2001). 70. Glover, M. L. & Reed, M. D. Propylene glycol: The safe diluent that continues to cause harm. Pharmacotherapy 16, 690–693 (1996).

814

71. Alcohol Ethoxylates. (2007).

815

72. Belanger, S. E. et al. Aquatic risk assessment of alcohol ethoxylates in North

816

America and Europe. Ecotoxicol. Environ. Saf. 64, 85–99 (2006).

39

ACS Paragon Plus Environment

Page 40 of 41

Page 41 of 41

Environmental Science & Technology

817

73. Chikh, N., Cheikh-Zouaoui, M., Aggoun, S. & Duval, R. Effects of calcium nitrate

818

and triisopropanolamine on the setting and strength evolution of Portland

819

cement pastes. Mater. Struct. 41, 31–36 (2008).

820

74. Zeiger, E. B., Anderson B., Haworth S., Lawlor T. & Speck W. Salmonella

821

mutagenicity tests: III. Results from the testing of 255 chemicals. Enivron

822

Mutagen 1–110 (1987).

823 824

75. Metcalf and Eddy. Wastewater engineering: treatment and resource recovery. 5, (McGraw-Hill, 2014).

825

40

ACS Paragon Plus Environment