Perfluoroalkyl Acids in Great Lakes Precipitation and Surface Water

Jun 24, 2019 - Perfluoroalkyl Acids in Great Lakes Precipitation and Surface Water .... were greater in Lake Ontario precipitation than in precipitati...
0 downloads 0 Views 1MB Size
Subscriber access provided by BUFFALO STATE

Characterization of Natural and Affected Environments

Perfluoroalkyl Acids in Great Lakes Precipitation and Surface Water (2006 – 2018) Indicate Response to Phase-outs, Regulatory Action, and Variability in Fate and Transport Processes Sarah B Gewurtz, Lisa E. Bradley, Sean Backus, Alice Dove, Daryl McGoldrick, Hayley Hung, and Helena Dryfhout-Clark Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b01337 • Publication Date (Web): 24 Jun 2019 Downloaded from pubs.acs.org on July 18, 2019

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 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 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.

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 35

Environmental Science & Technology

1

Perfluoroalkyl Acids in Great Lakes Precipitation

2

and Surface Water (2006 – 2018) Indicate Response

3

to Phase-outs, Regulatory Action, and Variability in

4

Fate and Transport Processes

5

Sarah B. Gewurtz,† Lisa E. Bradley,*,† Sean Backus,† Alice Dove,† Daryl McGoldrick,† Hayley

6

Hung,‡ Helena Dryfhout-Clark,‡

7 8



Water Quality Monitoring and Surveillance, Environment and Climate Change Canada, 867

9 10 11

Lakeshore Road, Burlington, Ontario L7S 1A1, Canada ‡

Air Quality Processes Research Section, Environment and Climate Change Canada, 4905 Dufferin Street, Toronto, Ontario M3H 5T4, Canada

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 35

12

ABSTRACT: Perfluoroalkyl acids (PFAAs) were determined in precipitation from three

13

locations across the Great Lakes between 2006 and 2018 and compared to surface water.

14

Perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) concentrations generally

15

decreased in precipitation, likely in response to phase-outs/regulatory actions. In comparison,

16

shorter-chained PFAAs, which are not regulated in Canada, did not decrease and

17

perfluorohexanoate and perfluorobutanoate (PFBA) recently increased, which could be due to their

18

use as replacements as the longer-chained PFAAs are phased-out by industry. PFOS and PFOA

19

concentrations were greater in Lake Ontario precipitation than in the more remote locations. In

20

comparison, PFBA concentrations were comparable across locations, suggesting greater

21

atmospheric transport either through its more volatile precursors and/or directly in association with

22

particles/aerosols. In Lake Ontario, comparison of PFAAs in precipitation to surface water

23

provides evidence of sources (e.g., street dust and wastewater effluent) in addition to wet

24

deposition to surface water, whereas wet deposition appears to be dominant in Lakes Huron and

25

Superior. Our results suggest that source control of shorter-chained PFAAs may be slow to be

26

reflected in environmental concentrations due to emissions far from the location of detection and

27

continued volatilization from existing in-use products and waste streams.

ACS Paragon Plus Environment

2

Page 3 of 35

28

Environmental Science & Technology

TABLE OF CONTENTS (TOC)/ABSTRACT ART

29 30

ACS Paragon Plus Environment

3

Environmental Science & Technology

31

Page 4 of 35

INTRODUCTION

32

Concern over perfluoroalkyl acids (PFAAs) has increased in recent years despite phase-outs and

33

regulatory actions ongoing since the early 2000s.1,2,3 Much of this concern is due to the derivation

34

of regulatory criteria and advisory guidance for PFAAs that are lower than previous values1,4,5 as

35

well as the finding that the drinking water supply of approximately six million U.S. residents may

36

have concentrations of perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) above

37

the United States Environmental Protection Agency (USEPA) advisory level of 70 ng/L.2

38

The fate and transport of PFAAs are complex and difficult to predict because, unlike many other

39

persistent organic contaminants, they are confounded by their surfactant-properties and more

40

volatile precursors.6 As such, detailed long-term assessments are critical for determining if

41

environmental concentrations are responding to phase-outs, regulatory actions, and control

42

measures on a vast array of sources. Precipitation and surface water are ideal media for evaluating

43

PFAAs. Despite the fact that PFAAs are not considered volatile under environmental conditions,

44

they are commonly detected in the atmosphere as a result of their association with particles and

45

aerosols and following transformation of their more volatile precursors.7,8,9 Precipitation reflects

46

atmospheric concentrations10 and is an important removal mechanism of PFAAs from the

47

atmosphere through sequestration of particles11 and/or sorption of PFAAs from the atmosphere to

48

water droplets.10,12 Precipitation has been shown to impact PFAA concentrations in surface

49

water,13 groundwater,14 and soil.15 In comparison, surface water can reflect a variety of sources,

50

including point source releases in addition to atmospheric deposition.13,16,17,18,19 Surface water is

51

also relevant in terms of human exposure given its importance with respect to drinking water,

52

which is a major source of PFAAs to humans.20 Therefore, simultaneous assessment of PFAA

53

concentrations in precipitation and surface water can generate hypotheses related to fate and

ACS Paragon Plus Environment

4

Page 5 of 35

Environmental Science & Technology

54

transport and the relative importance of atmospheric deposition versus local anthropogenic sources

55

in influencing surface water concentrations. This information has implications for source control.

56

Chemicals whose fate and transport are largely influenced by atmospheric transport and wet

57

deposition would be more difficult to control as they may have been emitted to the environment

58

far from the location of detection.

59

The first objective of this study was to use a high-resolution dataset to assess the time trends of

60

PFAAs in precipitation at three locations across the Laurentian Great Lakes between 2006 and

61

2018 to test the hypothesis that phase-outs and regulatory measures have resulted in reduced

62

concentrations of PFAAs in the environment.

63

differences in the role of atmospheric deposition as a source to surface water. Precipitation

64

samples were integrated on a monthly basis and compared to open surface water measurements

65

collected three or four times across the Great Lakes between 2008 and 2017. These datasets were

66

generated in support of Environment and Climate Change Canada (ECCC)’s mandate to monitor

67

the quality of the environment for Canadians and priorities under the Chemicals Management Plan.

68

Given that the Great Lakes span 244,106 km2 and range from highly urbanized Lakes Erie and

69

Ontario to the more remote Lake Superior, the datasets also provide a unique opportunity to assess

70

spatial differences in fate and transport and the relative importance of the atmosphere as a potential

71

source. To our knowledge, this is the first study to evaluate high-resolution temporal trends of

72

PFAAs in precipitation over more than ten years.

The second objective was to assess spatial

73 74

MATERIALS AND METHODS

75

Sample Collection. Integrated monthly precipitation samples were collected between April

76

2006 and August 2018 from three locations in the Great Lakes basin: Point Petre, Burnt

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 35

77

Island/Evansville, and Sibley (Supporting Information (SI) Figure SI1).

Point Petre

78

(43.8428, -77.1536) is located in eastern Lake Ontario and is considered a lightly populated rural

79

region, 160 km east of Toronto, ON and approximately 85 km north of Rochester, NY.21 Burnt

80

Island (45.8283, -82.94801)/Evansville (45.8181, -81.3386) is considered a remote site that is

81

located on the southwestern end of Manitoulin Island in northern Lake Huron, about 140 km

82

southeast of Sault Ste. Marie, ON.21 Burnt Island was the station of sample collection between

83

April 2006 and March 2013, after which the site was shut down. In July 2013, the site was

84

relocated to Evansville, 22 km to the east on Manitoulin Island. The Evansville site is not close to

85

significant sources and, due to its close proximity to Burnt Island, the relocation is not expected to

86

influence atmospheric concentrations of PFAAs. Sibley (48.4986, -88.6694) is located on the

87

northwest corner of Lake Superior and is considered a remote location about 90 km east of Thunder

88

Bay, ON and approximately 300 km northeast of Duluth, Minnesota.21 Monthly integrated

89

precipitation samples were collected using a MIC-B wet-only precipitation automated sampler

90

(Meteorological Instruments of Canada, Thornhill, ON). The samples were collected in a sample

91

bucket lined with a polyethylene bag with an opening of 0.0314 m2. After each month, the samples

92

were transferred to 2 L high-density polyethylene (HDPE) bottles for shipment. The excess

93

portion of any sample greater than 2 L was measured in a graduated cylinder and then discarded.

94

The samples were transported to ECCC’s Canada Centre for Inland Waters (CCIW, Burlington,

95

ON, Canada) and stored at -10 ºC until analysis.

96

Surface water samples were collected using a pole sampler from a 1 m depth at open water

97

stations in Lakes Ontario (March/April 2008, 2012, 2013, and 2015), Erie (April/May 2012, 2013,

98

2014, and 2015), Huron/Georgian Bay (May/June 2012, 2014, and 2017), and Superior (May 2008

99

and 2016) (SI Figure SI1 and SI Table SI1). HDPE sample bottles were rinsed twice with sample

ACS Paragon Plus Environment

6

Page 7 of 35

Environmental Science & Technology

100

water before filling. Samples were collected in the spring when there is a high degree of mixing

101

in the water column.22 Following collection, the samples were transported to CCIW and stored at

102

4 °C until analysis.

103

Chemical and Instrumental Analysis. Extraction and instrumental analysis of PFAAs in

104

surface water and precipitation samples collected up until November 2017 were conducted by SGS

105

AXYS Analytical Services Ltd. (SGS AXYS, Victoria, BC, Canada) with SGS AXYS Method

106

MLA-060. Precipitation samples collected from December 2017 onwards were analyzed using

107

SGS AXYS Method MLA-110. In brief, MLA-060 and MLA-110 involve analysis of samples by

108

liquid chromatography-tandem mass spectrometry (LC-MS/MS), following solid-phase extraction

109

(SPE) and selective elution procedures as described in the Supporting Information and detailed

110

previously.23 For both methods, the separation employed a gradient elution with a mobile phase

111

comprised of acetonitrile and aqueous ammonium acetate and calibration curves were based on

112

six to eight calibration levels. PFAA quantification was performed using an isotope dilution

113

internal standard approach by comparing the area of the quantification ion to that of the labeled

114

standard and correcting for response factor. The samples were spiked with surrogate standards

115

prior to extraction and recovery standards prior to analysis with LC-MS/MS (SI Table SI2).

116

The SGS AXYS Methods MLA-060 and MLA-110 are accredited to ISO 17025 standards by

117

the Canadian Association for Laboratory Accreditation Inc. (CALA) and SGS AXYS routinely

118

participates in inter-laboratory comparison studies. Both methods use similar SPE and LC-MS/MS

119

approaches and data from the two methods are considered equivalent. Details on differences and

120

method equivalence are shown in the Supporting Information. A procedural blank and spiked

121

reference sample, consisting of a commercially available brand of bottled drinking water that was

122

tested by lot for PFAA background, were analyzed with each sample batch. PFAAs were not

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 35

123

detected in any of the 74 blank samples analyzed except for trace level detections listed in the

124

Supporting Information. Mean recoveries in native PFAA spiked matrices ranged from 97% to

125

100% (n = 59). 5% of samples within each sample set were analyzed in duplicate. Duplicate

126

samples were required to have a relative percent deviation of less than 40%. Reporting limits were

127

determined by prorating the concentration of the lowest calibration limit for the sample size and

128

extract volume. Concentrations were quantified if the peak met the retention time requirement for

129

a given chemical and the signal to noise ratio was greater than three. Detected concentrations in

130

precipitation that were in the region above the 3:1 signal to noise ratio and below the reporting

131

limit were J-flagged, which denotes estimated values.

132

Data Analysis. Data analysis was performed on PFAAs that were detected (including J-flagged

133

data) in greater than 30% of the samples, namely perfluorobutanoate (PFBA), perfluoropentanoate

134

(PFPeA),

135

perfluorononanoate (PFNA), perfluorodecanoate (PFDA), perfluoroundecanoate (PFUnA), and

136

perfluorododecanoate (PFDoA). Given the relatively low detection frequency ranging from 5%

137

to 65% (30% to 86% including J-flagged data) (SI Table SI3), methods for statistical analysis were

138

selected for their ability to handle censored measurements. We treated J-flagged data as detected

139

measurements in statistical analysis. However, they were distinguished from detected data above

140

the reporting limit in tables and figures. Statistical analysis of precipitation time trends with the

141

Seasonal Kendall trend test and spatial patterns were completed on u-scores, as recommended by

142

Helsel24 for datasets with multiple reporting limits. The u-score was calculated as the sum of the

143

algebraic sign of differences comparing the ith observation to all other observations within the

144

same variable as follows:

145

=∑

perfluorohexanoate

(



(PFHxA),

perfluoroheptanoate

(PFHpA),

PFOA,

) ACS Paragon Plus Environment

8

Page 9 of 35

Environmental Science & Technology

146

where ui is the u-score of the ith observation and xi and xk are the concentration values of the ith

147

and kth observations, respectively. The u-score forms the basis of the well-known Mann-Whitney

148

test and is related to Kendall’s tau and the calculation of Kaplan-Meier percentiles.24 In the

149

calculation of u-scores, when a detected concentration cannot be determined to be higher or lower

150

than a non-detect value, the sign of the difference is considered zero.

151

Previous studies have found that PFAA depositional flux is greatest at the beginning of a rainfall

152

event (first few millimeters of deposition)25 and concentrations may become diluted when the

153

amount of precipitation is high (> 2 mm).10 As such, we completed the statistical analysis for

154

PFAAs in precipitation on both a concentration- and flux-basis, where flux (ng/m2) = PFAA

155

concentration in the sample (ng/L) × volume of precipitation during the sampling period (L) /

156

surface area of the bucket used for sampling (0.0314 m2). The relationships between PFAA

157

concentration u-scores and volume of precipitation and between volume of precipitation and time

158

were evaluated with the Mann-Kendall and Seasonal Kendall trend tests, respectively.26

159

Time trends were analyzed using two methods: the Seasonal Kendall trend test26 on u-scores and

160

the Digital Filtration (DF) technique on concentration values.27,28,29 Only detected and J-flagged

161

data were included in the DF technique, as has been done previously.28 The F-test was used to

162

assess the significance of the coefficient of determination (r2) between the trend derived by the DF

163

model and the measured data. First-order half-lives or doubling times were approximated as

164

ln(2)/k where k was defined as the linear regression slope of the DF-estimated concentrations

165

versus time. It should be noted that PFAAs did not necessarily decline/increase linearly or

166

consistently in the first order manner throughout the entire monitoring period. Therefore, half-

167

lives and doubling times are only used to compare the relative rates of decline/increase in

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 35

168

concentrations between chemicals. Additional technical details of the DF technique are presented

169

in the Supporting Information.

170

Principal component analysis (PCA) was conducted to explore overall spatial patterns. Factor

171

loadings were rotated using the Varimax normalized rotation. Two factors were retained based on

172

the results of the Scree test. PCA assumes a constant rate of change between points in space.24 In

173

order to evaluate if this assumption influenced the results of the PCA, we also applied nonmetric

174

multidimensional scaling (NMDS), which plots distances between points in the same rank order

175

as distances in a resemblance matrix between sites, and makes no assumptions about linear

176

relationships.24 We used Euclidean distance to calculate the resemblance matrix. The initial runs

177

of the PCA and NMDS found that the results were highly influenced by non-detect values of

178

PFPeA, PFUnA, and PFDoA, which were only detected in 51%, 47%, and 32% of the samples,

179

respectively. As such, the PCA and NMDS were re-run and presented with the following PFAAs:

180

PFBA, PFHxA, PFHpA, PFOA, PFNA, and PFDA.

181

Finally, analysis of similarity (ANOSIM), using Euclidean distance to calculate the resemblance

182

matrix, was used to test for significant differences in PFAA patterns between locations and

183

between precipitation and surface water. ANOSIM is a nonparametric multivariate procedure that

184

operates on ranks of resemblances and uses permutation tests to obtain p-values.24

185

Statistical analyses were performed at the α = 0.05 significance level. R Version 3.4.4 was used

186

for all statistical methods except for the DF model application. U-scores were calculated using the

187

R script uscore.r.30 The Mann-Kendall and Seasonal Kendall tests for trends were conducted using

188

the kendallTrendTest and kendallSeasonalTrendTest functions of the EnvStats package.31 The

189

PCA was conducted using the principal function of the psych package.32 The NMDS and

ACS Paragon Plus Environment

10

Page 11 of 35

Environmental Science & Technology

190

ANOSIM tests were applied using the metaMDS and anosim functions, respectively, of the vegan

191

package.33

192 193

RESULTS AND DISCUSSION

194

PFAAs in Precipitation and Surface Water. PFAA concentrations in precipitation and surface

195

water collected in the Great Lakes are presented in SI Tables SI4 and SI5, respectively, and

196

summarized in SI Table SI3. Concentrations of PFAAs in Great Lakes precipitation were

197

comparable to other locations observed in previous studies, particularly in remote and rural

198

locations, although sometimes in heavily urbanized regions as well.9,13,34,35 This is not surprising

199

given that the more volatile precursors of some PFAAs are transported far distances in the

200

atmosphere, leading to their occurrence in remote locations,28 as discussed below. In addition,

201

PFAAs can be emitted directly to air in the form of aerosols or particles, which can also be

202

transported through the atmosphere.7,8,9

203

comparable to previous open water measurements from the Great Lakes.23,36

Surface water concentrations observed here were

204

Canadian Federal Environmental Quality Guidelines (FEQGs) exist for PFOS and the maximum

205

concentrations of PFOS observed in surface water (7.4 ng/L in Lake Ontario on April 21, 2015)

206

and precipitation (14 ng/L at Point Petre in February 2007) were more than two orders of

207

magnitude below the FEQG for surface water of 6,800 ng/L that is intended to protect all forms of

208

aquatic life for indefinite exposure periods.37 Furthermore, Health Canada has published drinking

209

water screening values for the majority of PFAA compounds evaluated in this study that range

210

from 20 ng/L for PFNA to 30,000 ng/L for PFBA38 and the PFAA concentrations observed here

211

were at least 8-fold lower than these values. It should be noted, however, that the maximum PFAA

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 35

212

concentrations observed in surface water in this study were just below the lowest screening and

213

guidance values published by regulatory agencies throughout the world.1

214

Time Trends. Time trends of PFAA concentrations in Great Lakes precipitation between April

215

2006 and August 2018 are presented in SI Figures SI2a, b, and c. Time trends for four selected

216

PFAAs are presented in Figure 1 (PFBA and PFHxA) and Figure 2 (PFOA and PFOS). Although

217

time trends of PFAAs in Great Lakes precipitation vary among chemicals and locations, some

218

overall patterns are evident. For example, PFOA, PFNA, PFDA, and PFOS, for the most part,

219

significantly decreased over time (Seasonal Kendall trend test, p < 0.05, SI Table SI6) with

220

approximate half-lives of less than 15 years as determined through application of the DF technique.

221

PFNA and PFDA at Sibley and PFDA at Burnt Island/Evansville were the only exceptions with

222

longer half-lives of 22, 67, and 20 years, respectively. These overall decreasing trends were likely

223

in response to North-American phase-outs and regulatory actions for PFOS, PFOA, and other long-

224

chained perfluoroalkyl carboxylates (PFCAs) and their precursors ongoing since the early

225

2000s.1,39 Although PFUnA and PFDoA are also long-chained PFCAs, their concentrations did

226

not monotonically decrease between April 2006 and August 2018. PFUnA and PFDoA are less

227

water soluble than PFOA, PFNA, PFDA, and PFOA40 and percent detection was relatively low

228

(46% and 31%, respectively) for these chemicals in Great Lakes precipitation (SI Table SI3),

229

which may be influencing their observed trends. PFBA, PFPeA, PFHxA, and PFHpA did not

230

decrease monotonically at the three Great Lakes sites between April 2006 and August 2018

231

(Seasonal Kendall trend test, p > 0.05), with the exception of PFBA at Burnt Island/Evansville (p

232

= 0.02) and PFHpA at Point Petre (p = 0.023). Furthermore, the DF trend line indicates that

233

concentrations of PFBA and PFHxA appear to be increasing in most recent years (~2010/2014 to

234

2018 at Point Petre and Burnt Island and ~2016 to 2018 at Sibley) (Figure 1). These results may

ACS Paragon Plus Environment

12

Page 13 of 35

Environmental Science & Technology

235

not be surprising considering that these short-chained PFCAs and their precursors have not been

236

regulated in Canada37 and are being used as replacement chemicals for PFOS, PFOA, and the

237

longer-chain PFCAs.41 For example, PFBA is an impurity in perfluorobutane sulfonyl fluoride-

238

based derivatives currently used as a replacement for perfluorooctane sulfonyl fluoride in surface

239

treatment products.41 In addition, in response to the 2010/2015 USEPA voluntary stewardship

240

program, most aqueous film forming foam (AFFF) manufacturers have switched to short-chain

241

(C6) fluorotelomer-based per- and polyfluoroalkyl substances (PFAS).42

242

Time trends of PFAA concentrations in Great Lakes surface water were not statistically

243

evaluated due to the fact that only two (Superior), three (Huron), and four (Erie and Ontario)

244

years/sampling events were evaluated between 2008 and 2017. More than four years of data are

245

typically required to detect the changes in PFAA concentrations observed in the environment.43

246

For PFAAs detected in Lake Ontario (i.e., the shorter chained-PFCAs containing eight or fewer

247

carbon atoms and PFOS), surface water concentrations were within the upper range of those

248

observed in precipitation. A similar pattern was observed in Huron surface water except for PFBA.

249

The concentrations of PFBA in Huron surface water was comparable to precipitation. The

250

frequency of detection of PFAAs was low in Superior surface water (SI Table SI3). Only two

251

chemicals, PFBA and PFOA, with detection frequencies of 36% and 55%, were detected in

252

Superior surface water, at concentrations comparable and slightly elevated compared to

253

precipitation, respectively. These results are discussed below.

254

PFAA concentration u-scores in precipitation decreased significantly with precipitation volume

255

in about 50% of the chemicals/locations (Mann-Kendall, p < 0.05, SI Table SI7). There was no

256

monotonic relationship (Mann-Kendall, p > 0.05) between PFAA concentration u-scores and

257

precipitation volume for the remaining chemicals/locations. These results correspond to previous

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 35

258

studies10,25 and indicate that PFAA concentrations in precipitation are, in some cases, dependent

259

on the volume of precipitation within each sampling month. Precipitation volume did not change

260

significantly with time at Point Petre (Seasonal Kendall trend test, p = 0.60, tau = -0.03, Z = -0.53)

261

or Burnt Island/Evansville (Seasonal Kendall trend test, p = 0.77, tau = 0.019, Z = 0.29) but

262

increased at Sibley (Seasonal Kendall trend test, p = 0.0010, tau = 0.25, Z = 3.3). However, the

263

overall time trends of PFAAs in Great Lakes precipitation were consistent whether evaluated on a

264

concentration- or flux-basis (SI Figure SI2 and SI Table SI6).

265

The DF-derived seasonal trend lines indicated considerable intra-annual variability of PFAAs in

266

precipitation on both a concentration- and flux-basis (Figure 1, Figure 2, and SI Figure SI2).

267

However, no consistent seasonal patterns were evident, which corresponds to previous findings

268

for PFAAs in precipitation collected from two locations in Japan between June 2006 and June

269

200810 and PFOA and PFNA at three northeastern U.S. sites in 1998 and 1999.35

270

Spatial Patterns. The PCA conducted on PFAA concentration (Figure 2) and flux (SI Figure

271

SI3) u-scores indicated that precipitation collected at Point Petre in Lake Ontario grouped

272

separately from Burnt Island/Evansville (Huron) and Sibley (Superior) with relatively high scores

273

for factor 2, which was generally associated with PFOS. PFOA also loaded relatively high on

274

factor 2, with values slightly below PFOS. Application of ANOSIM showed that the pattern of

275

concentrations and fluxes at Point Petre were significantly different than Burnt Island (R = 0.17

276

and 0.15 for concentration and flux, respectively, p = 0.001) and Sibley (R = 0.19 and 0.24 for

277

concentration and flux, respectively, p = 0.001) but there were no significant differences between

278

Burnt Island and Sibley (R = 0.012 and -5.7×10-5, p = 0.078 and 0.48 for concentration and flux,

279

respectively). Point Petre is considered a lightly populated rural region. However, it is relatively

280

close (160 km and 85 km, respectively) to major urban centers of Metropolitan Toronto, ON

ACS Paragon Plus Environment

14

Page 15 of 35

Environmental Science & Technology

281

(population = 6.4 million) and Rochester, NY (population = 208,046). As such, it is likely that

282

relatively elevated PFOS and PFOA concentrations in precipitation at this location results from

283

the influence of human activity, as found in other studies for PFOS and in some cases

284

PFOA.9,44,45,46 PFOS and PFOA are not considered volatile; however, they can be released into

285

the gas phase from point sources indirectly in the form of their more volatile precursors,

286

transported in the atmosphere in the gas phase, and then degraded to PFOS and PFOA.9,47,48,49

287

PFOS and PFOA can also be emitted from point sources to the atmosphere directly in the form of

288

particles and/or aerosols and travel through the air associated with particles/aerosols.8,49,50,51

289

Recent studies show that PFAAs can enrich significantly in sea spray aerosols and be emitted to

290

the atmosphere;11,52,53 although only about 3% of the aerosolized PFOA and PFOS were estimated

291

to be transported and deposited to land areas.52 In comparison, all three locations had similar factor

292

1 scores, which were generally associated with PFBA and other short-chained PFCAs. PFBA,

293

other short-chain PFCAs, and their precursors such as the fluorotelomer alcohols (FTOHs) are

294

more volatile compared to the longer-chained PFAAs, are found at higher concentrations in gas

295

compared to particle phase in air, and are subject to transport in the atmosphere (either through

296

their precursors and/or associated with particles/aerosols).28,40,54 Furthermore, as discussed by

297

Wong et al.,28 atmospheric transformation of certain volatile hydrofluorocarbons and

298

hydrofluoroethers to PFBAs could have occurred, which has been demonstrated in laboratory

299

studies under zero or low NOx conditions.55,56 This provides evidence that the major source of

300

PFBA and other short-chain PFCAs to Great Lakes precipitation is through long-range

301

atmospheric transport, which leads to similar concentrations throughout Great Lakes precipitation.

302

Our results correspond with those of Wong et al.28 who found that air concentrations of PFBA at

303

Alert, located in the Canadian high Arctic (median = 1.7 pg/m3, range: