Fossil Fuel-Derived Polycyclic Aromatic Hydrocarbons in the Taiwan

PCGC system consisted of an Agilent 6890 GC with a FID and a 7683 Series auto-injector,. 183 combined with a cooled injection system (CIS, Gerstel) an...
0 downloads 0 Views 758KB Size
Subscriber access provided by University of Winnipeg Library

Environmental Processes

Fossil Fuel-Derived Polycyclic Aromatic Hydrocarbons in the Taiwan Strait, China, and Fluxes across the Air-Water Interface Miaolei Ya, Li Xu, Yu-Ling Wu, Yongyu Li, Songhe Zhao, and Xin-Hong Wang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01331 • Publication Date (Web): 01 Jun 2018 Downloaded from http://pubs.acs.org on June 1, 2018

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 25

Environmental Science & Technology

1

Fossil Fuel-Derived Polycyclic Aromatic Hydrocarbons in the Taiwan Strait, China,

2

and Fluxes across the Air-Water Interface

3

Miaolei Ya1, 2, Li Xu2 *, Yuling Wu1, Yongyu Li1, Songhe Zhao1 and Xinhong Wang1 *.

4

1

5

Ecology, Xiamen University, Xiamen 361102, China.

6

2

7

and Geophysics, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts

8

02543, United States.

State Key Laboratory of Marine Environmental Science, College of the Environment & National Ocean Sciences Accelerator Mass Spectrometry Facility, Department of Geology

9 10

ABSTRACT

11

Based on the application of compound-specific radiocarbon analysis (CSRA) and air-

12

water exchange models, the contributions of fossil fuel and biomass burning derived

13

polycyclic aromatic hydrocarbons (PAHs) as well as their air-water transport were

14

elucidated. The results showed that fossil fuel-derived PAHs (an average contribution of

15

89%) present the net volatilization process at the air-water interface of the Taiwan Strait in

16

summer. Net volatile fluxes of the dominant fluorene and phenanthrene (>58% of the total

17

PAHs) were 27±2.8 μg m−2∙day−1, significantly higher than the dry deposition fluxes

18

(average 0.43 μg m−2∙day−1). The Δ14C contents of selected PAHs (fluorene, phenanthrene

19

plus anthracene, fluoranthene, and pyrene) determined by CSRA in the dissolved seawater

20

ranged from −997±4‰ to −873±6‰, indicating 89−100% (95±4%) of PAHs were

21

contributed by fossil fuels. The South China Sea warm current originating from the

22

southwest China in summer (98%) and the Min-Zhe coastal current originating from the

23

north China in winter (97%) input more fossil fuel PAHs than the Jiulong River estuary

24

(90%) and Xiamen harbor water (93%). The more radioactive decayed 14C of fluoranthene

25

(a 4- ring PAH) than phenanthrene and anthracene (3- ring PAHs) represented a greater

26

fossil fuel contribution to the former in dissolved seawater.

27

28 1

ACS Paragon Plus Environment

Environmental Science & Technology

29

Page 2 of 25

1. INTRODUCTION

30

The origin of polycyclic aromatic hydrocarbons (PAHs) and air-water interface

31

transport are important environmental processes affecting the exist and fate of PAHs in

32

marine environment.1, 2 Generally, sources of PAHs to the coastal ocean include, but are

33

not limited to, the incomplete combustion and pyrolysis of fossil fuels (e.g., coal, petroleum)

34

or modern biomass (e.g., wood, straw and grass and other C3 and C4 plants) and the direct

35

release of petroleum products.3, 4 The research on the origin of PAHs is to find out which

36

type of human activity is the most important way for PAHs to import into the offshore

37

ocean. After the land-emitting and ocean-based anthropogenic PAHs are transmitted to

38

marine environment, they are transported across regional seas under the action of coastal

39

water masses.5-7 Synchronously, the dynamic exchanges (atmospheric deposition and free

40

exchange) of PAHs occur at the air-water interface.8, 9 Air-water transport of PAHs is

41

proved to be a critical intermediate process that affects their migration to the global ocean

42

and burial into the sedimentary environment.10-12 PAH species vary in their volatilization

43

/deposition processes at the air-water interface.13 Typically, the air-water interface is a

44

source of low molecular weight (LMW) PAHs and a sink of high molecular weight (HMW)

45

PAHs in the coastal ocean.14,

46

radiocarbon analysis (CSRA) and air-water exchange models to quantitatively study the

47

sources (fossil fuel vs. biomass burning) of PAHs in surface seawater and interface

48

transport (deposition or volatilization), respectively, in the coastal water masses of the

49

Taiwan Strait (TWS), connecting the South China Sea (SCS) and the East China Sea (ECS).

50

Traditionally, the source identifications of PAHs were widely studied using the isomer

51

diagnostic ratios and receptor models such as principal component analysis (PCA), and

52

positive matrix factorization, etc.1, 16-18 However, variations in combustion conditions and

53

environmental degradation processes occurring during the transport of PAH isomers from

54

their sources to receptors 19 undermines the application of diagnostic ratios and statistical

55

methods as reliable source apportionment tools.20 Therefore, the advent of CSRA based on

56

the separation of individual compounds by preparative capillary gas chromatography

57

(PCGC) augments isotope composition heterogeneity by the molecular level

15

Therefore, we applied advanced compound-specific

21, 22

14

C-

58

Accelerator Mass Spectrometry (AMS) measurement.

59

deconvoluting PAH sources in complex, heterogeneous matrices (soil, sedimentary and

60

atmospheric environments).23-26 These studies have capitalized on the difference in

61

radiocarbon contents between fossil fuel sources of PAHs, which contain no significant 2

ACS Paragon Plus Environment

This has proven useful for

Page 3 of 25

Environmental Science & Technology

62

radiocarbon (Δ14C = −1000‰; Δ14C is expressed as a deviation from a known “oxalic acid”

63

standard, which normalizes the radiocarbon content of a sample to the same δ13C (−25‰)

64

and time point (1950)27, 28), and biomass burning sources, which contain modern levels of

65

radiocarbon (such as wood burning with a Δ14C of +225‰ or leaves and annual grass

66

burning with a Δ14C of +70‰).29 Using a radiocarbon isotopic mass balance approach, the

67

relative contributions of these sources to PAHs presented in environmental matrixes can be

68

apportioned.30 However, limited by the complexity and high demands of ultra-microscale

69

(ca. 5.0−25 μg C) radiocarbon measurements of individual pure compounds,24, 31 CSRA

70

method was used only in the identification of PAHs sources in the atmospheric, soil and

71

sediment environments.23-26, 29, 30, 32 Up to now, this method has not been applied to the

72

source identification of PAHs in the aqueous medium, limited by the large volume

73

enrichment of trace PAHs (e.g., phenanthrene (Phen): 21−71 ng L−1 in our previous study33,

74

34

75

sedimentation and burial of PAHs.35, 36 Therefore, it is a vital and worthy research topic to

76

extend the application of CSRA to assess the source characteristic of PAHs in the aquatic

77

environment.

) in seawater. However, coastal seawater plays an important role in global transport,

78

In this study, we focused on the CSRA application in four typical water masses in the

79

TWS,37 which represent different PAH origins and also have differing effects on the

80

transport of anthropogenic PAHs. In detail, Jiulong River diluted water (JRDW), which has

81

the largest runoff discharge in summer, likely contains PAHs emitted from industrial and

82

agricultural activities (e.g., sewage discharges and surface runoff) in surrounding rural and

83

urban areas along the its drainage path (Figure 1).33, 38 Xiamen Western Harbor (XMWH),

84

the biggest cargo terminal in the western TWS, can be assumed to contain fossil fuel-

85

derived PAHs from the release and combustion of petroleum in the port (e.g., shipping

86

activities) and from urban vehicle-derived emissions.34 Min-Zhe coastal current (MZCC),

87

driven by the northeast monsoon in winter, could transport anthropogenic PAHs from the

88

Yangtze River and Hangzhou Bay along the coast of the ECS.39 The South China Sea warm

89

current (SCSWC), originating from the Southeast Asian coasts, represents the long-range

90

transport of PAHs driven by the southwest monsoon in summer.5, 40 At present study, we

91

determine the source (fossil fuel or biomass burning) and their contributions of PAHs by

92

the measurement of the 14C of individual PAHs, based on CSRA in the above-mentioned

93

water masses. Additionally, we also reveal the transport characteristic of fossil fuel versus

94

biomass burning derived PAHs at the air-water interface by examining transport behaviors

95

(volatilization or deposition and their fluxes) of PAHs in the TWS. This will not only 3

ACS Paragon Plus Environment

Environmental Science & Technology

96

provide scientific advice on the reduction of PAHs, but also have important implications

97

for studying the migration and burial of PAHs in the vertical direction of the offshore

98

seawater.

99

2. MATERIALS AND METHODS

100

2.1. Sampling of PAHs.

101

Sampling for Concentration Analysis During the summer (June and July) of 2010,

102

13 surface seawater (sites T1 to T13; ~1 m depth) and 3 aerosol samples (sites TW01,

103

TW02, and TW03: TW01 covered site T1-T6 area; TW02 covered site T6-T10 area; TW03

104

covered site T10-end area) were taken from the R/V Yanping 2 in the central TWS (solid

105

black spots in Figure 1). Seawater sampling procedures have been described in our previous

106

studies (see section 1 of the Supporting Information, SI),5, 41 glass-fiber filters (GF/Fs) and

107

solid-phase extraction (SPE) cartridges were used to collect suspended particulate matters

108

(SPM, particulate phase) and dissolved phase samples. Synchronously, aerosol sampling

109

was carried out on the forecastle of the ship (upwind of the chimney, about 5 m high) using

110

large flow air samplers (Thermo Electron Corporation) while sailing against the wind to

111

avoid the contamination from the ship, such as chimney. Quartz filters (O.D. 100 mm,

112

Advantec Corporation, pre-burned at 600 °C for 6 hours) and polyurethane foams (PUF,

113

O.D. 65 mm, length 40 mm; Soxhlet extracted with dichloromethane for 24 h) were used

114

to synchronously collect total suspended particulates (TSP, aerosol phase) and gas phase

115

(PUF phase) samples at a flow speed of 360 L min−1. Finally, GF/Fs and SPE cartridges

116

from seawater sampling and TSP and PUF from air sampling were hermetically stored at

117

−20 °C until further analysis.

118

4

ACS Paragon Plus Environment

Page 4 of 25

Page 5 of 25

Environmental Science & Technology

119 120

Figure 1 Sampling areas in the TWS. The roman numerals from I to IV show the bulk seawater

121

sampling locations where the Δ14C values of selected PAHs were measured. MZCC: Min-Zhe

122

coastal current; SCSWC: South China Sea warm current; JRDW: Jiulong River diluted water;

123

XMWH: Xiamen Western Harbor.

124 125

Bulk Sampling for 14C Analysis To measure the Δ14C and δ13C values (relative to

126

VPDB: Vienna Pee Dee Belemnite) of selected PAHs dissolved in TWS seawater, bulk

127

seawater samples (5000–20000 L) were collected from the previously described four

128

regional water masses, including the MZCC (I), SCSWC (II), JRDW (III) and XMWH (IV)

129

(Figure 1). After initial filtration by GF/Fs (Whatman, O.D. 142 mm, 0.75 μm), dissolved

130

organic matter (DOM, 1500 mL of XAD-2 resin) were used to extract the DOM to collect enough

136

target dissolved PAHs from the seawater. Prior to use, XAD-2 resin was regenerated

137

following standard procedures (see section 2 of the SI) and ultrasonically extracted with 5

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 25

138

ethyl alcohol (all solvents from TEDIA Company, USA). The XAD-2 resin was then

139

transferred into the glass columns mentioned above. Ethyl alcohol was removed from the

140

gaps under a gentle stream of N2 (99.99%). Then, the XAD-2 columns were pre-extracted

141

with 100 mL acetone, 100 mL dichloromethane (DCM) and 100 mL n-hexane at a slow

142

flow speed (drop-by-drop) to remove contaminants in the matrix. Finally, plenty of

143

ultrapure water (18.25 MΩ·cm, Millipore Company) was used to elute the residual solvents.

144

After field sampling, used XAD-2 columns were stored at 4 °C until further pretreatment

145

for CSRA steps.

146

2.2. Concentration Analysis.

147

The analytical procedures of PAHs concentrations for the GF/F and SPE seawater

148

samples have been described in more detail in our previous research.33,

149

pretreatment procedure for atmospheric TSP samples was the same as for the GF/F samples.

150

For the atmospheric PUF samples, Soxhlet extraction was carried out for 24 h after spiking

151

with five deuterated PAH surrogates. After rotary evaporation, the concentrated extract was

152

cleaned using pretreated alumina/silica gel chromatography. Finally, the eluate was blown

153

gently under a stream of N2 gas to 100 μL and stored at –20 °C before instrumental analysis.

154

The detailed sample pretreatment methods and instrumental analysis by Agilent 6890

155

Series gas chromatography–Agilent 5973 mass spectrometry are described in our previous

156

papers.5,

157

designated by the United States Environmental Protection Agency. Their abbreviations are

158

shown in section 3 of the SI.

159

2.3. 14C Preparation and Analysis.

41

34, 38

The

The target compounds including 15 priority PAHs (excluding naphthalene)

160

Based on a previous field investigation of individual PAH concentrations (e.g.,

161

phenanthrene (Phen): 21−71 ng L−1) in the seawater of the TWS,33, 34 seawater sampling

162

for several thousand liters is required to isolate adequate quantities of PAHs for 14C-AMS

163

measurement.

164

For the bulk seawater samples, residual seawater was drained from the XAD-2 column

165

under a gentle stream of N2 and then discarded. Each XAD-2 column was eluted with 100

166

mL methanol followed by 200 mL DCM at a slow flow rate (drop by drop). Residual

167

solvents were removed from the XAD-2 under a gentle stream of N2. The mixed extracts

168

were transferred into a 1 L separatory funnel, to which 150 mL of ultrapure water was

169

added. After sufficient mixing and subsequent phase separation, 100 mL DCM was added

6

ACS Paragon Plus Environment

Page 7 of 25

Environmental Science & Technology

170

3x for liquid/liquid extractions. DCM extracts were concentrated to near dryness and the

171

solvent was exchanged with 5 mL n-hexane. After being concentrated to ~1 mL, extracts

172

were fractionated in a glass column (i.d.15 mm) packed with 10 g of activated silica (0.063–

173

0.20 mm mesh, activated at 450 °C for 4 h and deactivated at 130 °C for 12 h). The elution

174

solvents used are as follows: fraction 1, 25 mL of hexane followed by 10 mL of DCM/n-

175

hexane (v: v, 1:1); fraction 2, 40 mL of DCM/n-hexane. Fraction 2, containing the PAHs,

176

was concentrated and filtered through a syringe-tip filter (Millipore 4 mm Millex-FH, pore

177

size 0.2 µm). Finally, all extracts were combined and concentrated to ~1 mL. The PAH

178

fraction was cleaned by partitioning with n-pentane/dimethylformamide, dried with

179

anhydrous sodium sulfate,42 and reduced to the proper volume under a gentle stream of N2.

180

PAHs were quantified using gas chromatography with a flame ionization detector (GC-

181

FID).

182

After purification, PCGC was used to isolate and prepare individual PAHs.21, 29, 43 The

183

PCGC system consisted of an Agilent 6890 GC with a FID and a 7683 Series auto-injector,

184

combined with a cooled injection system (CIS, Gerstel) and a preparative fraction collector

185

(PFC, Gerstel) operated at 320 °C. Details of the PCGC method used for PAH isolation are

186

listed in section 4 of the SI. Briefly, individual PAHs were isolated using two GC columns

187

(30 m Agilent HP5 and 30 m Agilent DB17) in sequence and collected by the PFC 2x.

188

Approximately 1% of the column eluate was diverted to the FID detector and the remaining

189

99% of the sample was collected into glass U-tubes by the PFC. The individual PAHs

190

including fluorene (Flu), Phen plus anthracene (An), fluoranthene (Fluo), and pyrene (Py)

191

were collected into the independent glass U-tubes. Phen and An cannot be separated by

192

DB17 chromatographic column, so we combined these two isomers together as one

193

component for

194

acenaphthylene (Ace), acenaphthene (Acen), benzo(a)anthracene (BaA), chrysene (Chry),

195

and 5- and 6- ring PAHs (see section 3 of the SI) in the dissolved phase were too low to be

196

performed for radiocarbon analysis. Trapped individual PAHs were transferred to 4 mL

197

vials by rinsing with n-hexane 4x. Finally, the samples were passed through a silica gel

198

column (4.0 cm × 0.5 cm i.d.) to remove column bleed from the PCGC. A small aliquot

199

from each trap was taken to check the purity of the PAHs by GC-FID.

14

C analysis. The abundances of the other detected PAHs such as

200

Carbon isotope analysis was performed at the National Ocean Sciences Accelerator

201

Mass Spectrometry (NOSAMS) facility at Woods Hole Oceanographic Institute. After the

202

catalytic oxidation of the target PAHs, the generated CO2 was purified and quantified in

203

the vacuum glass line. Finally, the CO2 was reduced to graphite and Δ14C measurements 7

ACS Paragon Plus Environment

Environmental Science & Technology

204

were made by AMS.44 Detailed procedures are shown elsewhere.29 The Δ14C values based

205

on the sampling year were compared to the NBS oxalic acid I (NIST SRM 4990) standard

206

after normalizing the radiocarbon content of a sample to the same δ13C (−25‰) and time

207

point (1950).27, 45

208

2.4. Quality Assurance and Quality Control.

209

We strictly performed analytic procedures to ensure quality control. Field blanks of

210

GF/Fs, quartz filters, SPE cartridges and PUF were taken to represent potential

211

contamination during handling on the ship. Concentrations of target PAHs in the field blank

212

samples were ~2 orders of magnitude less than in the seawater samples, which means we

213

can largely ignore the PAH background from the ship and laboratory. For the water samples,

214

the method recovery of the 15 PAHs ranged from 74% to 126% in the dissolved phase (SPE

215

cartridges) and from 60% to 120% in the particulate phase (baked GF/Fs). For the

216

atmospheric samples, the method recovery ranged from 78% to 105% in the aerosol

217

samples (baked quartz filters) and from 72% to 113% in the gas samples (clean PUF). The

218

method detection limits (MDLs) were calculated as the mean plus three times the standard

219

deviation of the field blanks. For the water samples, MDLs of the 15 PAHs ranged from

220

0.9–4.1 pg L−1 in the dissolved phase and from 0.8–3.5 pg L−1 in the particulate phase

221

(considering an average sampling volume of 8 L). For the atmospheric samples, MDLs

222

ranged from 0.7–2.4×102 pg m−3 in the aerosol phase and from 1.8–4.3×102 pg m−3 in the

223

gas phase (considering the sampled volume of 400 m3).

224

Prior to bulk seawater sampling, the XAD-2 resin was purified, in order, with acetone,

225

DCM and n-hexane to remove contaminants in the matrix. The mass of PAHs remaining in

226

the 100 mL of purified XAD-2 resin ranged from n. d. (not detected) to 25 ng. This showed

227

that PAHs in the matrix were 2 orders of magnitude less abundant than in the field samples

228

(>17 µg). PAHs in the solvents used in the pretreatment procedure were not detected and

229

other impurities in the solvent or introduced during sample handling were removed during

230

the PCGC step. Solvent from the final column bleed cleanup (less than 4 mL) was blew to

231

dryness in a baked quartz tube before combustion. Therefore, to the trapped target PAHs,

232

the pretreatment, PCGC collection and column bleed removal steps could also not cause

233

additional carbon contamination. Studies have shown that the retention time deviation of

234

the PCGC collection process is the main cause of isotope fractionation. Therefore, we

235

strictly controlled the deviation of retention time during collection. And our previous work

236

has shown no isotope fractionation of

14

C during the PCGC collection step.29, 46 Finally, 8

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25

Environmental Science & Technology

237

multiple PCGC blanks were combined to a concentration of 0.15 μmol C, and the fraction

238

modern value of the combined blank was 0.29.

239

3. RESULTS AND DISCUSSION

240

3.1. PAHs in Water and Air of the TWS.

241

The sum of the 15 PAH concentrations (Σ15PAHs) in the central TWS (sites T1 to T13,

242

n=13) ranged from 53 to 79 (62±8.3) ng L−1 in the dissolved seawater (Table S1 and Figure

243

S1), and from 2.5 to 5.6 (4.0±1.1) ng L−1 in the particulate phase in seawater (Table S2).

244

As the dominant phase of PAHs (>91%, Figure S2), dissolved PAH concentrations showed

245

an even distribution in the central TWS (Figures 2a and S3). The composition of PAHs also

246

showed no significant difference among the different stations (p>0.05). Phen (28±4.0 ng

247

L−1) followed by Flu (20±2.9 ng L−1) displayed the highest relative abundance with average

248

contributions of 46% and 32% to Σ15PAHs across the TWS, respectively. Phen and Flu, the

249

predominant components in ultrafine particles of petroleum combustion in the engines,47

250

might derived from the emissions of vehicles and ships in surrounding cities and cruise

251

routes in the TWS. The abundances of these two compounds were higher than found in the

252

open southern Indian Ocean (Phen: 1.7 ng L−1; Flu: 1.4 ng L−1),6 close to the levels in the

253

ECS (Phen: 19 ng L−1; Flu: 32 ng L−1) and the northern SCS (Phen: 14 ng L−1; Flu: 28 ng

254

L−1),5 and significantly lower than those found in the coast of the Bohai Sea in China (Flu:

255

132 ng L−1 and Phen: 449 ng L−1).48 Along the west coast of the TWS, Σ15PAHs in the

256

representative estuarine (site JRDW, 145 ng L−1), harbor (site XMWH, 98 ng L−1) and

257

coastal currents (site MZCC, 101 ng L−1) were significantly higher than those in the central

258

TWS (62±8.3 ng L−1). Σ15PAHs values in the southern TWS (site SCSWC) were lower due

259

to weaker anthropogenic inputs and dilution by the adjacent open seawater (Figures 2a, S1,

260

and S2). The cluster analysis showed PAHs in the different water masses and the central

261

TWS were distributed in the different groups (Figure S3), which also reflected the

262

significant regional difference of the PAHs composition.

263

As in the surface seawater, no significant differences were found in the total or

264

individual concentrations of PAHs in the aerosol samples over the central TWS (sites

265

TW01 to TW03) (Figures 2b, S1, and S2). Σ15PAHs values ranged from 23−31 and 0.2−0.4

266

ng m−3 in the gas and aerosol phases, respectively. In the dominant gas phase, Phen (16±3.4

267

ng m−3) displayed the highest abundance, with an average contribution of 59% to Σ15PAHs.

268

Discharge from the urban agglomerations on both sides of the TWS and volatilization from

9

ACS Paragon Plus Environment

Environmental Science & Technology

269

seawater resulted in 2−5x higher Σ15PAHs values than in the ECS (8 ng m−3)14 and the

270

northwestern Pacific Ocean (6 ng m−3).49 Along the west coast of the TWS, the Σ15PAHs

271

value in the gas phase of the JRDW (47 ng m−3) was ~2x higher than over the central TWS

272

due to significantly higher Fluo and Py inputs from the surrounding cities (Figure 2b).

273

Compared with the central TWS, we found a significantly lower Phen concentration (2.3

274

ng m−3) while the higher Ace, Acen and Flu concentrations in the southern TWS (Figure

275

2b), which was similar to h the open SCS in our previous study.5 The different composition

276

of PAHs in the central TWS was probably owing to the source differences and the effect of

277

long-range transport of the air in summer (Figure S4).

278

279 280

Figure 2. Relative levels of PAHs in the surface dissolved seawater and the gas phase of the air in

281

the TWS.

282 283

Some studies have suggested that soil is an important natural/biogenic source of

284

PAHs.50 Enhanced surface runoff and soil erosion in summer are the crucial input pathways

285

of these PAHs to the coastal ocean.5, 39 In our study area, the southern TWS is influenced

286

by the combined inputs of Pearl River diluted water and the SCSWC (Figure 1).51, 52

287

Therefore, the northeastward flowing SCSWC (e.g., site SCSWC) and surface runoff

288

inputs (e.g., JRDW) could be seen as two crucial sources of PAHs in the central TWS.

289

These two end-members, with distinctly different concentrations and compositions of

290

PAHs, adequately mix together, driven by the confluence of the water masses and 10

ACS Paragon Plus Environment

Page 10 of 25

Page 11 of 25

Environmental Science & Technology

291

atmospheric diffusion.53 Accompanied by photo- and/or bio- degradation and the air-water

292

transport of PAHs in this mixing process, PAHs displayed uniform characteristics

293

(concentrations, compositions and phase partitioning) in the central TWS.

294

According to backtrack air mass trajectories, the air of the central TWS was affected

295

by the long-range transport of air from the eastern SCS, Philippine islands, and the western

296

Pacific (Figure S4). During this transport, photo-degradation, atmospheric deposition, air-

297

water exchange of PAHs and mixing processes with terrigenous air masses resulted in

298

significant variations in PAH concentrations and compositions in the atmosphere of the

299

central TWS compared with these surrounding area (Figure 2). According to their forward

300

trajectories, PAHs in the atmosphere of the central TWS would be transported to South

301

Korea and Japan in the coming days (about 120 hours) after mixing with atmospheric PAHs

302

from the northern Chinese mainland (Figure S4). PAHs in surface seawater of the TWS

303

showed the same transport routes as PAHs in the air because of the wind-driven ocean

304

currents along the Chinese coast in summer.52

305

3.2. PAHs Transport at the Air-Water Interface.

306

The fugacity gradient of semi-volatile chemicals mathematically describes the

307

potential migration direction in which chemicals diffuse or are transported at the air-water

308

interface.9, 14, 54, 55 The fugacity ratio (fa/fw) and the fluxes of PAHs (Fa/w, ng m−2 day−1) and

309

their uncertainties were used to discuss the transport behaviors (volatilization or deposition

310

and their fluxes). The detailed description and calculations are listed in sections 5 and 6 of

311

the SI.

312

In the TWS, uncertainties of the fugacity gradient in summer were ~52% (calculated

313

in section 5 of the SI), controlled by the gas phase and dissolved concentrations of PAHs,

314

ambient temperatures and Henry’s law constant.7, 56 Therefore, 0.59 < fa/fw < 1.7 (standard

315

deviation in ln(fa/fw)=0.52) were considered to not significantly differ from phase

316

equilibrium. 3- ring PAHs showed significant volatilization across the air-water interface

317

of the TWS, which includes the dominant components of the PAH pools (Flu and Phen,

318

averaging 58% of total PAHs in the gas phase and 67% in dissolved water) (Figure 3a).

319

Fluo and Py distributions (averaging 25% of total PAHs in the gas phase and 11% in

320

dissolved water) highlighted the depositional process in the relatively open central TWS

321

and SCSWC. However, Fluo and Py distributions in the JRDW highlighted the exchange

322

equilibrium, because the local PAH emissions into the water were from the Jiulong River,

323

harbors and coastal cities. The opposite transport directions of the isomeric BaA 11

ACS Paragon Plus Environment

Environmental Science & Technology

324

(volatilization) and Chry (deposition) in the TWS are a function of the differences of their

325

sources and degradation rates in air and surface water.57

326

Overall, the transfer velocities (Kol) of the individual PAHs (Table S3) in the TWS

327

averaged 0.73±0.04 m day−1, a function of wind speed and temperature.54 Determined Kol

328

values are comparable to previous reports from Chinese coasts.14 The relative uncertainties

329

of transport fluxes of individual PAHs ranged from 44% to 142%, calculated in section 6

330

of the SI, a function of Kol error propagation, compound concentrations, ambient

331

temperatures and Henry’s law constant.8, 56, 58 In the central TWS (sites TW01-TW03),

332

volatilization fluxes of the air-water exchange (Fa/w) of the dominant PAHs: Flu and Phen,

333

were as high as 14±0.9 and 12±2.7 μg m−2 day−1, respectively (Figure 3b and Table S4).

334

These fluxes are comparable to the ECS,14 lower than in the Bohai and Yellow Seas,59 but

335

higher than at the southern tip of Taiwan Island.8 This means surface seawater was the most

336

important secondary source of 3- ring PAHs to the atmosphere of the TWS. These PAHs

337

were mainly derived from riverine inputs,33 adjacent ocean inputs,5 sediment release by

338

upwelling41 and offshore oil pollution.60, 61 This further confirmed the volatile behavior of

339

LMW PAHs in the marine ocean is an important driving factor of the cross-sea transport

340

of anthropogenic PAHs.13, 14

341

In the JRDW, the volatile fluxes of LMW PAHs (Ace, Acen, and Flu, 28% of the total

342

PAHs) were 1.5−2x higher than those in the central TWS (Figure 3b and Table S4). The

343

volatile fluxes of Py, BaA, and Chry were caused by their greater concentrations in

344

seawater due to the inputs of PAHs from the Jiulong River and the surrounding cities. Fluo,

345

Py and Chry (4- ring PAHs) showed the statistical exchange equilibrium (0.59 < fa/fw < 1.7)

346

in the JRDW, but the Py and Chry still had the volatile trends with the higher fluxes (Figure

347

3a and 3b). The higher volatile flux of BaA than that in the central TWS (Figure 3b) could

348

be attributed to the higher fugacity gradient of BaA between the JRDW and the surrounding

349

urban atmosphere. In the SCSWC, the volatile fluxes of the dominant PAHs: Flu and Phen,

350

were 1.5−2x lower than in the central TWS due to the continual air-water interaction of

351

PAHs during their long-range transport from the SCS (Figures 1 and S4). In other words,

352

controlled by the surrounding input, migration and transformation characteristics of PAHs

353

in the different regions, fugacity gradient of individual PAHs in the air and water ultimately

354

determined their transport behaviors (volatilization or deposition and their fluxes) of PAHs

355

in the TWS.

356

12

ACS Paragon Plus Environment

Page 12 of 25

Page 13 of 25

Environmental Science & Technology

357 358

Figure 3. Air-water interface transport of selected PAHs. (a) Fugacity ratios, the

359

horizontal gray area (standard deviation in ln(fa/fw)= 0.52) were considered to not

360

significantly differ from phase equilibrium; (b) Air-water exchange fluxes (Fa/w, μg

361

m−2 day−1) and their relative uncertainties (top x-axis, %); (c) Dry deposition fluxes

362

(Fd, μg m−2 day−1); (d) Total fluxes (Fa/w+d, μg m−2 day−1).

363 364

It has been an effort to measure the dry deposition fluxes (Fd, μg m−2 day−1) of PAHs

365

in the open ocean.62, 63 The detailed description and calculations are listed in section 7 of

366

the SI. In the TWS, the estimated dry deposition fluxes (Fd) of 3- plus 4- ring PAHs in the

367

JRDW (56 μg m−2 day−1) and SCSWC (23 μg m−2 day−1) were one order of magnitude

368

higher than those in the central TWS (1.3±0.5 μg m−2 day−1) (Figure 3c). In the JRDW, that

369

was because of the local higher particulate PAH discharge in the atmosphere; but in the

370

SCSWC, that could be due to the higher loadings of PAHs on the fine particulates with

371

larger specific surface areas in the long-range transport of SCSWC.64 In addition, because

372

of discontinuity and much shorter durations of rainfall, wet deposition is considered to be

373

at least ten times less contribution to PAHs in the air-water interface than dry deposition.65

374

Therefore, examining the combined impacts of free air-water exchange and dry deposition

375

(Fa/w+d), the transport fluxes of PAHs at the air-water interface in the central TWS showed

376

significant differences with the JRDW and SCSWC (Figure 3d). In the central TWS, the

13

ACS Paragon Plus Environment

Environmental Science & Technology

377

transport of 3- and 4- ring PAHs across the air-water interface were dominated by

378

volatilization and deposition processes, respectively, based on the Fa/w values of these PAHs

379

(Figures 3c and 3d). However, Fd of Phen and An exceeded the volatile Fa/w values in the

380

SCSWC and JRDW. In the JRDW, there were higher depositional fluxes (Fa/w+d) of Phen

381

and Fluo, but other 3- and 4- ring PAHs showed primarily volatilization (Figure 3d). In the

382

SCSWC, Phen, An and 4- ring PAH fluxes were dominated by deposition at the air-water

383

interface and the LMW PAHs (Ace, Acen, and Flu) showed net volatilization (Figure 3d).

384

Similar transport characteristics of PAHs at the air-water interface have been reported in

385

many other sea areas.2, 14

386

3.3. Fossil Fuel vs. Biomass Burning Contributions of PAHs by CSRA.

387

Based on the available concentrated carbon weight (Table S5), carbon isotopic (Δ14C

388

and δ13C) measurements of four selected individual and isomeric PAHs (Flu, Phen plus An,

389

Fluo and Py) dissolved in the seawater from four water masses around the TWS were made

390

by advanced CSRA technology. The total amount of compound mass for Δ14C analysis

391

ranged from 1.8 to 5.0 μmol C, and the purities of the isolated samples were in the range

392

of 92 to 99% (Figure S5). Limited by the sample sizes of the purified PAHs, only δ13C

393

measurements of Phen plus An (−29±0.6‰) in the XMWH were obtained. The result

394

showed a similar δ13C composition as atmospheric particulate PAHs in Croatia and Greece

395

(−29 and −28‰),66 but slightly lighter than the δ13C values of atmospheric PAHs (−29 to

396

−24‰),24, 29 PAHs in central European forest soils (−25 to −23‰)26 and sedimentary PAHs

397

(−27 to −25‰).30 Although the δ13C signature has been attempted to be used to identify the

398

sources of PAHs,67-69 its application is extremely difficult because of the overlapping end-

399

member values of various PAH sources.68, 70, 71

400

In the four water masses around the TWS, the measured Δ14C contents of selected

401

PAHs ranged from −997±4‰ to −873±6‰ (Table S5 and Figure 4), which indicated

402

decayed carbon contributions, in terms of 14C content. The depleted values were explained

403

by the emission of PAHs originating from fossil contributions, including the pyrolysis of

404

fossil fuels and the release of petroleum products.3, 4 The measured Δ14C values of PAHs

405

around the TWS were comparable to those in the atmosphere of North Birmingham in the

406

southeastern US (−990 to −911‰),29 background areas of Croatia and Greece (−941 and

407

−888‰),66 oil sands in Northern Alberta, Canada (−962 to −849‰)32 and forest soils in

408

central Europe (−942 to −819‰).26 In contrast, Δ14C values of selected PAHs were more

409

depleted than those reported in sediment cores from Siskiwit (−783 to −388‰),72 14

ACS Paragon Plus Environment

Page 14 of 25

Page 15 of 25

Environmental Science & Technology

410

atmospheric PAHs in Sweden (−388 and −381‰),30 residential areas of suburban Tokyo

411

(−787 to −514‰)25 and the Western Balkans (−573 to −288‰),24 which all revealed a

412

greater contribution from biomass burning. The regional differences in Δ14C values implied

413

the significant regional differences in the contribution of biomass burning and fossil fuel

414

to PAHs. To a certain extent, this also indirectly reflected the regional differences in energy

415

structure.

416

At the air-water interface of the TWS, in consideration of the net volatilization of Flu,

417

Phen and An (3- ring PAHs) and the net deposition of Fluo and Py (4- ring PAHs) (Figure

418

3d), the depleted Δ14C values of Flu, Phen and An in the water would lead to a gradual

419

increase of estimated fossil fuel contributions to 3- ring PAHs in the atmosphere. Likewise,

420

the depleted Δ14C values of Fluo and Py could also be influenced by the atmospheric

421

deposition of PAHs derived from fossil fuels or the combustion of fossil fuels.

422

423 424

Figure 4. Δ14C (‰) and ffossil fuel (%) of PAHs determined by CSRA in the water masses of the

425

TWS

426 427

Δ14C values of individual PAHs in the different water masses, in increasing order,

428

were SCSWC, MZCC, XMWH, and JRDW (Figure 4). This indicates, from lower to higher

429

values, a gradually decreasing proportion of fossil fuel contributions to the PAH pool. In

430

the SCSWC, the ocean-based anthropogenic discharge of PAHs, including potential ship

431

emissions, oil spills and produced water discharge from offshore oil platforms could

432

explain the most depleted Δ14C values. In the MZCC, the coastal southwestward current

433

could carry Δ14C depleted PAHs derived from domestic coal burning for indoor heating

434

from the coast of the ECS during winter.5 In the XMWH, the local ports and the 15

ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 25

435

surrounding cities discharged the bulk of fossil fuel derived PAHs, and the latter also

436

exported a small portion of the biomass burning derived PAHs. In the JRDW, river erosion

437

carries PAHs derived from the countryside along its path, which likely contains more PAHs

438

derived from biomass burning and causes the higher Δ14C values of PAHs than those found

439

in the other water masses. Finally, comparing the Δ14C values of each individual PAH,

440

Δ14C values of Phen plus An (3- ring PAHs) were significantly higher those of Fluo (4-

441

ring PAHs) (Figure 4). This was likely due to the mixing effects of different input pathways

442

(atmospheric deposition, particle suspension, water mass transport, etc.) which loaded

443

different proportions of biomass burning versus fossil fuel derived PAHs.

444

We estimated fractional contributions of biomass burning (i.e., fbiomass) versus fossil

445

fuel (i.e., ffossil=1−fbiomass) derived sources to individual PAHs using a two-source mixing

446

method according to the mass balance equation of

447

Δ14Cfossil (1 − fbiomass),24, 25, 29 where, Δ14CPAH is the measured

448

fraction. The radiocarbon free Δ14C value of the fossil fuel end-member (i.e., Δ14Cfossil fuel)

449

is −1000‰.29 Since there is not one specific biomass burning source to the TWS, we

450

adapted a Δ14C value of 152‰ as the modern end-member in this study.29 Therefore, the

451

calculated percent contributions of biomass burning for the individual PAHs reported here

452

were all very low in the different water masses around the TWS, with a mean of 10±0.5%

453

in JRDW, 7.5±0.7% in XMWH, 2.8±0.8% in MZCC and 1.7±0.6% in SCSWC (Table S5

454

and Figure 4). These values are very different compared with total annual global

455

atmospheric emissions, which are around 60% derived from residential/commercial

456

biomass burning.73 This discrepancy could be explained by the dominant anthropogenic

457

input sources of PAHs from the ocean or surrounding regional cities in the TWS. Primarily

458

fossil fuel derived PAHs have been found in many aquatic and atmospheric environments

459

in China.18, 74, 75 Finally, we could convert the biomass burning and fossil fuel contributions

460

to the corresponding concentrations of PAHs. That result showed the contributions of the

461

leakage and combustion of fossil fuels to Phen and An were the main factor of PAHs

462

pollution in the seawater of the Western TWS in summer (Figure 5). Detailly, in the MZCC,

463

XMWH, JRDW and SCSWC, the concentrations of fossil fuel derived Flu, Phen, An and

464

Fluo were 35, 76, 68 and 18 ng L−1, respectively; and the concentrations of biomass burning

465

derived Flu, Phen, An and Fluo were 1.2, 6.8, 8.1 and 0.4 ng L−1, respectively (Table S5).

14

C: Δ14CPAH =Δ14Cbiomass fbiomass +

466

16

ACS Paragon Plus Environment

14

C content of each PAH

Page 17 of 25

Environmental Science & Technology

467 468

Figure 5. The calculated weight of fossil fuel and biomass burning derived individual PAHs in

469

the dissolved phase of seawater around the TWS.

470 471

Principal component analysis-multiple linear regression (PCA-MLR) (Figure S6) and

472

the diagnostic ratios (Figure S7) that have been used widely to identify the sources of PAHs

473

were compared in this study. For the PCA-MLR analysis in this study (see section 8 of the

474

SI), the result showed that the first principal component, with high loadings of the 3- ring

475

PAHs and Fluo and Py, contributed 92% of the PAH variability. The other principal

476

components, with high loadings of HMW PAHs, contributed to 8.0% of the observed PAHs

477

variability. However, we could not definitively assign the principal components to fossil

478

fuel or biomass burning derived sources because of the lack of a clear source spectrum of

479

PAHs. In addition, selective degradation of PAHs in the environment from their sources to

480

sinks is substantial enough to lead serious error of source identification of PAHs by

481

diagnostic ratios.19 Therefore, our work suggests CSRA is the most advantageous method

482

for identifying the sources of PAHs in the environment.

483

Next, we will further apply CSRA to the study of individual PAHs in diverse

484

environmental media (such as atmospheric particulate phase, atmospheric gas phase,

485

seawater particle phase and sediments). If that, we can investigate the migration and

486

transport behavior of PAHs in the air-water-sediment system in the coast of TWS based on 17

ACS Paragon Plus Environment

Environmental Science & Technology

487

the changes of 13C and 14C compositions of PAHs.

488

ASSOCIATED CONTENT

489

Supporting Information

490

Details on sample pretreatment methods, calculations of air-water fugacity ratios, air-

491

water fluxes and their uncertainty, datasets of PAH concentrations, the backward and

492

forward trajectories of air masses, and other auxiliary figures and tables are given in the

493

supporting information. This information is available free of charge via the Internet at

494

http://pubs.acs.org.

495

AUTHOR INFORMATION

496

Corresponding Authors

497

*Xinhong Wang: Phone: +86 592 2187857; email: [email protected].

498

*Li Xu: Phone: +1 508 289 3673; email: [email protected].

499

Notes

500

The authors declare no competing financial interest.

501

ACKNOWLEDGMENTS

502

This work was supported by the National Natural Science Foundation of China (NSFC)

503

Project (41276066). We thank Meihui Lin from the Fujian Marine Forecasts for providing

504

sampling support in the Min-Zhe coastal current. We also thank the crew of the R/V Hai

505

Yang 2 and R/V Yan Ping 2 for sampling in the Taiwan Strait and Xiamen Western Harbor,

506

respectively.

507

REFERENCES

508 509 510 511 512 513 514 515 516

1. Yunker, M. B.; Macdonald, R. W.; Vingarzan, R.; Mitchell, R. H.; Goyette, D.; Sylvestre, S., PAHs in the Fraser River basin: a critical appraisal of PAH ratios as indicators of PAH source and composition. Org. Geochem. 2002, 33 (4), 489-515; DOI 10.1016/S0146-6380(02)00002-5. 2. Nizzetto, L.; Lohmann, R.; Gioia, R.; Jahnke, A.; Temme, C.; Dachs, J.; Herckes, P.; Guardo, A. D.; Jones, K. C., PAHs in air and seawater along a North-South Atlantic transect: trends, processes and possible sources. Environ. Sci. Technol. 2008, 42 (5), 1580-1585; DOI 10.1021/es0717414. 3. Manzetti, S., Polycyclic aromatic hydrocarbons in the environment: environmental 18

ACS Paragon Plus Environment

Page 18 of 25

Page 19 of 25

517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555

Environmental Science & Technology

fate and transformation. Polycycl. Aromat. Comp. 2013, 33 (4), 311-330; DOI 10.1080/10406638.2013.781042 4. Khairy, M. A.; Lohmann, R., Field validation of polyethylene passive air samplers for parent and alkylated PAHs in Alexandria, Egypt. Environ. Sci. Technol. 2012, 46 (7), 39903998; DOI 10.1021/es300012u. 5. Ya, M. L.; Wang, X. H.; Wu, Y. L.; Li, Y. Y.; Yan, J. M.; Fang, C.; Zhao, Y. Y.; Qian, R. R.; Lin, X. L., Seasonal variation of terrigenous polycyclic aromatic hydrocarbons along the marginal seas of china: input, phase partitioning, and ocean-current transport. Environ. Sci. Technol. 2017, 51 (16), 9072-9079; DOI 10.1021/acs.est.7b02755. 6. Cai, M. G.; Liu, M. Y.; Hong, Q. Q.; Lin, J.; Huang, P.; Hong, J. J.; Wang, J.; Zhao, W. L.; Chen, M.; Cai, M. H., Fate of polycyclic aromatic hydrocarbons in seawater from the Western Pacific to the Southern Ocean (17.5°N to 69.2°S) and their inventories on the Antarctic Shelf. Environ. Sci. Technol. 2016, 50 (17), 9161-9168; DOI 10.1021/acs.est.6b02766. 7. Lohmann, R.; Klanova, J.; Pribylova, P.; Liskova, H.; Yonis, S.; Bollinger, K., PAHs on a west-to-east transect across the tropical Atlantic Ocean. Environ. Sci. Technol. 2013, 47 (6), 2570-2578; DOI 10.1021/es304764e. 8. Cheng, J. O.; Ko, F. C.; Lee, C. L.; Fang, M. D., Air-water exchange fluxes of polycyclic aromatic hydrocarbons in the tropical coast, Taiwan. Chemosphere 2013, 90 (10), 2614-2622; DOI 10.1016/j.chemosphere.2012.11.020. 9. Lammel, G.; Audy, O.; Besis, A.; Efstathiou, C.; Eleftheriadis, K.; Kohoutek, J.; Kukučka, P.; Mulder, M. D.; Přibylová, P.; Prokeš, R., Air and seawater pollution and airsea gas exchange of persistent toxic substances in the Aegean Sea: spatial trends of PAHs, PCBs, OCPs and PBDEs. Environ. Sci. Pollut. R. 2015, 22 (15), 11301-11313; DOI 10.1007/s11356-015-4363-4. 10. Huang, Y. M.; Zhang, R. J.; Li, K. C.; Cheng, Z. N.; Zhong, G. C.; Zhang, G.; Li, J., Experimental study on the role of sedimentation and degradation processes on atmospheric deposition of persistent organic pollutants in a subtropical water column. Environ. Sci. Technol. 2017, 51 (8), 4424-4433; DOI 10.1021/acs.est.7b00568. 11. Castro-Jiménez, J.; Berrojalbiz, N.; Wollgast, J.; Dachs, J., Polycyclic aromatic hydrocarbons (PAHs) in the Mediterranean Sea: atmospheric occurrence, deposition and decoupling with settling fluxes in the water column. Environ. Pollut. 2012, 166, 40-47; DOI 10.1016/j.envpol.2012.03.003. 12. Dachs, J.; Eisenreich, S. J.; Hoff, R. M., Influence of eutrophication on air-water exchange, vertical fluxes, and phytoplankton concentrations of persistent organic pollutants. Environ. Sci. Technol. 2000, 34 (6), 1095-1102; DOI 10.1021/es990759e. 13. González-Gaya, B.; Fernández-Pinos, M. C.; Morales, L.; Méjanelle, L.; Abad, E.; Piña, B.; Duarte, C. M.; Jiménez, B.; Dachs, J., High atmosphere-ocean exchange of semivolatile aromatic hydrocarbons. Nat. Geosci. 2016, 9 (6), 438-442; DOI 19

ACS Paragon Plus Environment

Environmental Science & Technology

556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594

10.1038/ngeo2714. 14. Wu, Z. L.; Lin, T.; Li, Z. X.; Jiang, Y. Q.; Li, Y. Y.; Yao, X. H.; Gao, H. W.; Guo, Z. G., Air-sea exchange and gas-particle partitioning of polycyclic aromatic hydrocarbons over the northwestern Pacific Ocean: role of East Asian continental outflow. Environ. Pollut. 2017, 230, 444-452; DOI 10.1016/j.envpol.2017.06.079. 15. Qin, N.; He, W.; Kong, X. Z.; Liu, W. X.; He, Q. S.; Yang, B.; Ouyang, H. L.; Wang, Q. M.; Xu, F. L., Atmospheric partitioning and the air-water exchange of polycyclic aromatic hydrocarbons in a large shallow Chinese lake (Lake Chaohu). Chemosphere 2013, 93 (9), 1685-1693; DOI 10.1016/j.chemosphere.2013.05.038. 16. Larsen, R. K.; Baker, J. E., Source apportionment of polycyclic aromatic hydrocarbons in the urban atmosphere: a comparison of three methods. Environ. Sci. Technol. 2003, 37 (9), 1873-1881; DOI 10.1021/es0206184. 17. Yang, B.; Zhou, L. L.; Xue, N. D.; Li, F. S.; Li, Y. W.; Vogt, R. D.; Cong, X.; Yan, Y. Z.; Liu, B., Source apportionment of polycyclic aromatic hydrocarbons in soils of Huanghuai Plain, China: comparison of three receptor models. Sci. Total Environ. 2013, 443, 31-39; DOI 10.1016/j.scitotenv.2012.10.094. 18. Zhang, Y.; Guo, C. S.; Xu, J.; Tian, Y. Z.; Shi, G. L.; Feng, Y. C., Potential source contributions and risk assessment of PAHs in sediments from Taihu Lake, China: comparison of three receptor models. Water Res. 2012, 46 (9), 3065-3073; DOI 10.1016/j.watres.2012.03.006. 19. Zhang, X. L.; Tao, S.; Liu, W. X.; Yang, Y.; Zuo, Q.; Liu, S. Z., Source diagnostics of polycyclic aromatic hydrocarbons based on species ratios: a multimedia approach. Environ. Sci. Technol. 2005, 39 (23), 9109-9114; DOI 10.1021/es0513741. 20. Katsoyiannis, A.; Sweetman, A. J.; Jones, K. C., PAH molecular diagnostic ratios applied to atmospheric sources: a critical evaluation using two decades of source inventory and air concentration data from the UK. Environ. Sci. Technol. 2011, 45 (20), 8897-8906; DOI 10.1021/es202277u. 21. Eglinton, T. I.; Aluwihare, L. I.; Bauer, J. E.; Druffel, E. R.; McNichol, A. P., Gas chromatographic isolation of individual compounds from complex matrices for radiocarbon dating. Anal. Chem. 1996, 68 (5), 904-912; DOI 10.1021/ac9508513. 22. Pearson, A.; McNichol, A. P.; Benitez-Nelson, B. C.; Hayes, J. M.; Eglinton, T. I., Origins of lipid biomarkers in Santa Monica Basin surface sediment: a case study using compound-specific Δ14C analysis. Geochim. Cosmochim. Ac. 2001, 65 (18), 3123-3137; DOI 10.1016/S0016-7037(01)00657-3. 23. Sheesley, R. J.; Kruså, M.; Krecl, P.; Johansson, C.; Gustafsson, Ö., Source apportionment of elevated wintertime PAHs by compound-specific radiocarbon analysis. Atmos. Chem. Phys. 2009, 9 (10), 3347-3356; DOI 10.5194/acp-9-3347-2009. 24. Zencak, Z.; Klanova, J.; Holoubek, I.; Gustafsson, Ö., Source apportionment of atmospheric PAHs in the Western Balkans by natural abundance radiocarbon analysis. 20

ACS Paragon Plus Environment

Page 20 of 25

Page 21 of 25

595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633

Environmental Science & Technology

Environ. Sci. Technol. 2007, 41 (11), 3850-3855; DOI 10.1021/es0628957. 25. Kumata, H.; Uchida, M.; Sakuma, E.; Uchida, T.; Fujiwara, K.; Tsuzuki, M.; Yoneda, M.; Shibata, Y., Compound class specific 14C analysis of polycyclic aromatic hydrocarbons associated with PM10 and PM1.1 aerosols from residential areas of suburban Tokyo. Environ. Sci. Technol. 2006, 40 (11), 3474-3480; DOI 10.1021/es052407f. 26. Bosch, C.; Andersson, A.; Kruså, M.; Bandh, C.; Hovorková, I.; Klánová, J.; Knowles, T. D.; Pancost, R. D.; Evershed, R. P.; Gustafsson, O., Source apportionment of polycyclic aromatic hydrocarbons in central european soils with compound-specific triple isotopes (δ13C, Δ14C, and δ2H). Environ. Sci. Technol. 2015, 49 (13), 7657-7665; DOI 10.1021/acs.est.5b01190. 27. Stuiver, M.; Polach, H. A., Discussion reporting of 14C data. Radiocarbon 1977, 19 (03), 355-363; DOI 10.1017/S0033822200003672. 28. McNichol, A. P.; Aluwihare, L. I., The power of radiocarbon in biogeochemical studies of the marine carbon cycle: insights from studies of dissolved and particulate organic carbon (DOC and POC). Chem. Rev. 2007, 107 (2), 443-466; DOI 10.1021/cr050374g. 29. Xu, L.; Zheng, M.; Ding, X.; Edgerton, E. S.; Reddy, C. M., Modern and fossil contributions to polycyclic aromatic hydrocarbons in PM2.5 from North Birmingham, Alabama in the southeastern US. Environ. Sci. Technol. 2012, 46 (3), 1422-1429; DOI 10.1021/es2043189. 30. Mandalakis, M.; Gustafsson, Ö.; Reddy, C. M.; Xu, L., Radiocarbon apportionment of fossil versus biofuel combustion sources of polycyclic aromatic hydrocarbons in the Stockholm metropolitan area. Environ. Sci. Technol. 2004, 38 (20), 5344-5349; DOI 10.1021/es049088x. 31. Shah, S. R.; Pearson, A., Ultra-Microscale (5-25 μg C) analysis of individual lipids by 14 C AMS: assessment and correction for sample processing blanks. Radiocarbon 2007, 49 (1), 69-82; DOI 10.1017/S0033822200041904. 32. Jautzy, J. J.; Ahad, J. M. E.; Hall, R. I.; Wiklund, J. A.; Wolfe, B. B.; Gobeil, C.; Savard, M. M., Source apportionment of background PAHs in the Peace-Athabasca Delta (Alberta, Canada) using molecular level radiocarbon analysis. Environ. Sci. Technol. 2015, 49 (15), 9056-9063; DOI 10.1021/acs.est.5b01490. 33. Wu, Y. L.; Wang, X. H.; Li, Y. Y.; Hong, H. S., Occurrence of polycyclic aromatic hydrocarbons (PAHs) in seawater from the Western Taiwan Strait, China. Mar. Pollut. Bull. 2011, 63 (5), 459-463; DOI 10.1016/j.marpolbul.2011.03.008. 34. Ya, M. L.; Wang, X. H.; Wu, Y. L.; Ye, C. X.; Li, Y. Y., Enrichment and partitioning of polycyclic aromatic hydrocarbons in the sea surface microlayer and subsurface water along the coast of Xiamen Island, China. Mar. Pollut. Bull. 2014, 78 (1), 110-117; DOI 10.1016/j.marpolbul.2013.10.053. 35. Berrojalbiz, N.; Dachs, J.; Ojeda, M. J.; Valle, M. C.; Castro‐Jiménez, J.; Wollgast, J.; Ghiani, M.; Hanke, G.; Zaldivar, J. M., Biogeochemical and physical controls on 21

ACS Paragon Plus Environment

Environmental Science & Technology

634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672

concentrations of polycyclic aromatic hydrocarbons in water and plankton of the Mediterranean and Black Seas. Global. Biogeochem. Cy. 2011, 25 (4), GB4003; DOI 10.1029/2010GB003775. 36. Lang, C.; Tao, S.; Liu, W. X.; Zhang, Y. X.; Simonich, S. M., Atmospheric transport and outflow of polycyclic aromatic hydrocarbons from China. Environ. Sci. Technol. 2008, 42 (14), 5196-5201; DOI 10.1021/es800453n. 37. Hu, J. Y.; Kawamura, H.; Li, C. Y.; Hong, H. S.; Jiang, Y. W., Review on current and seawater volume transport through the Taiwan Strait. J. Oceanogr. 2010, 66 (5), 591-610; DOI 10.1007/s10872-010-0049-1. 38. Wu, Y. L.; Wang, X. H.; Li, Y. Y.; Ya, M. L.; Luo, H.; Hong, H. S., Polybrominated diphenyl ethers, organochlorine pesticides, and polycyclic aromatic hydrocarbons in water from the Jiulong River Estuary, China: levels, distributions, influencing factors, and risk assessment. Environ. Sci. Pollut. R. 2017, 24 (10), 8933-8945; DOI 10.1007/s11356-0154782-2. 39. Lin, T.; Hu, L. M.; Guo, Z. G.; Zhang, G.; Yang, Z. S., Deposition fluxes and fate of polycyclic aromatic hydrocarbons in the Yangtze River estuarine-inner shelf in the East China Sea. Global. Biogeochem. Cy. 2013, 27 (1), 77-87; DOI 10.1029/2012GB004317. 40. Liu, L. Y.; Wang, J. Z.; Wei, G. L.; Guan, Y. F.; Zeng, E. Y., Polycyclic aromatic hydrocarbons (PAHs) in continental shelf sediment of China: implications for anthropogenic influences on coastal marine environment. Environ. Pollut. 2012, 167, 155162; DOI 10.1016/j.envpol.2012.03.038. 41. Ya, M. L.; Wu, Y. L.; Li, Y. Y.; Wang, X. H., Transport of terrigenous polycyclic aromatic hydrocarbons affected by the coastal upwelling in the northwestern coast of South China Sea. Environ. Pollut. 2017, 229, 60-68; DOI 10.1016/j.envpol.2017.05.054. 42. Mandalakis, M.; Zebühr, Y.; Gustafsson, Ö., Efficient isolation of polyaromatic fraction from aliphatic compounds in complex extracts using dimethylformamide-pentane partitionings. J. Chromatogr. A 2004, 1041 (1), 111-117; DOI 10.1016/j.chroma.2004.04.038. 43. Zhang, X. Y.; Zhao, L.; Wang, Y. X.; Xu, Y. P.; Zhou, L. P., Optimization of programmed-temperature vaporization injection preparative capillary GC for compound specific radiocarbon analysis. J. Sep. Sci. 2013, 36 (13), 2136-2144; DOI 10.1002/jssc.201300088. 44. Pearson, A.; Mcnichol, A. P.; Schneider, R. J.; Von Reden, K. F.; Zheng, Y., Microscale AMS 14C measurement at NOSAMs. Radiocarbon 1997, 40 (1), 61-75; DOI 10.1017/S0033822200017902. 45. Von Reden, K. F.; McNichol, A. P.; Pearson, A.; Schneider, R. J., 14C AMS measurements of < 100 μg samples with a high-current system. Radiocarbon 1997, 40 (1), 247-253; DOI 10.1017/S0033822200018117. 46. Ball, G. I.; Xu, L.; McNichol, A. P.; Aluwihare, L. I., A two-dimensional, heart-cutting 22

ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25

673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711

Environmental Science & Technology

preparative gas chromatograph facilitates highly resolved single-compound isolations with utility towards compound-specific natural abundance radiocarbon (14C) analyses. J. Chromatogr. A 2012, 1220, 122-131; DOI 10.1016/j.chroma.2011.11.021. 47. Lu, T.; Huang, Z.; Cheung, C. S.; Ma, J., Size distribution of EC, OC and particlephase PAHs emissions from a diesel engine fueled with three fuels. Sci. Total Environ. 2012, 438, 33-41; DOI 10.1016/j.scitotenv.2012.08.026. 48. Zheng, B. H.; Wang, L. P.; Lei, K.; Nan, B. X., Distribution and ecological risk assessment of polycyclic aromatic hydrocarbons in water, suspended particulate matter and sediment from Daliao River estuary and the adjacent area, China. Chemosphere 2016, 149, 91-100; DOI 10.1016/j.chemosphere.2016.01.039. 49. Ma, Y. X.; Xie, Z. Y.; Yang, H. Z.; Möller, A.; Halsall, C.; Cai, M. H.; Sturm, R.; Ebinghaus, R., Deposition of polycyclic aromatic hydrocarbons in the North Pacific and the Arctic. J. Geophys. Res.: Atmos. 2013, 118 (11), 5822-5829; DOI 10.1002/jgrd.50473. 50. Cabrerizo, A.; Dachs, J.; Moeckel, C.; Ojeda, M. J.; Caballero, G.; Barceló, D.; Jones, K. C., Ubiquitous net volatilization of polycyclic aromatic hydrocarbons from soils and parameters influencing their soi-air partitioning. Environ. Sci. Technol. 2011, 45 (11), 47404747; DOI 10.1021/es104131f. 51. Wang, D. X.; Shu, Y. Q.; Xue, H. J.; Hu, J. Y.; Chen, J.; Zhuang, W.; Zu, T. T.; Xu, J. D., Relative contributions of local wind and topography to the coastal upwelling intensity in the northern South China Sea. J. Geophys. Res.: Oceans 2014, 119 (4), 2550-2567; DOI 10.1002/2013JC009172. 52. He, X. Q.; Xu, D. F.; Bai, Y.; Pan, D. L.; Chen, C. T. A.; Chen, X. Y.; Gong, F., Eddyentrained Pearl River plume into the oligotrophic basin of the South China Sea. Cont. Shelf Res. 2016, 124, 117-124; DOI 10.1016/j.csr.2016.06.003. 53. Inomata, Y.; Kajino, M.; Sato, K.; Kurokawa, J.; Tang, N.; Ohara, T.; Hayakawa, K.; Ueda, H., Source-receptor relationship analysis of the atmospheric deposition of PAHs subject to long-range transport in Northeast Asia. Environ. Sci. Technol. 2017, 51 (14), 7972-7981; DOI 10.1021/acs.est.7b00776. 54. Ruge, Z.; Muir, D.; Helm, P.; Lohmann, R., Concentrations, trends, and air-water exchange of PAHs and PBDEs derived from passive samplers in Lake superior in 2011. Environ. Sci. Technol. 2015, 49 (23), 13777-13786; DOI 10.1021/acs.est.5b02611. 55. Tidwell, L. G.; Allan, S. E.; O’Connell, S. G.; Hobbie, K. A.; Smith, B. W.; Anderson, K. A., Polycyclic aromatic hydrocarbon (PAH) and oxygenated PAH (OPAH) air-water exchange during the deepwater horizon oil spill. Environ. Sci. Technol. 2014, 49 (1), 141149; DOI 10.1021/es503827y. 56. Liu, Y.; Wang, S. Y.; McDonough, C. A.; Khairy, M.; Muir, D.; Lohmann, R., Estimation of uncertainty in air-water exchange flux and gross volatilization loss of PCBs: a case study based on passive sampling in the lower Great Lakes. Environ. Sci. Technol. 2016, 50 (20), 10894-10902; DOI 10.1021/acs.est.6b02891. 23

ACS Paragon Plus Environment

Environmental Science & Technology

712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750

57. Ghosal, D.; Ghosh, S.; Dutta, T. K.; Ahn, Y., Current state of knowledge in microbial degradation of polycyclic aromatic hydrocarbons (PAHs): a review. Front. Microbiol. 2016, 7, 1-27; DOI 10.3389/fmicb.2016.01369 58. Fang, M. D.; Lee, C. L.; Jiang, J. J.; Ko, F. C.; Baker, J. E., Diffusive exchange of PAHs across the air-water interface of the Kaohsiung Harbor lagoon, Taiwan. J. Environ. Manage. 2012, 110, 179-187; DOI 10.1016/j.jenvman.2012.06.001. 59. Chen, Y. J.; Lin, T.; Tang, J. H.; Xie, Z. Y.; Tian, C. G.; Li, J.; Zhang, G., Exchange of polycyclic aromatic hydrocarbons across the air-water interface in the Bohai and Yellow Seas. Atmos. Environ. 2016, 141, 153-160; DOI 10.1016/j.atmosenv.2016.06.039. 60. Harman, C.; Langford, K.; Sundt, R. C.; Brooks, S., Measurement of naphthenic acids in the receiving waters around an offshore oil platform by passive sampling. Environ. Toxicol. Chem. 2014, 33 (9), 1946-1949; DOI 10.1002/etc.2651. 61. Liu, B.; Chen, B.; Zhang, B. Y.; Jing, L.; Zhang, H.; Lee, K., Photocatalytic degradation of polycyclic aromatic hydrocarbons in offshore produced water: effects of water matrix. J. Environ. Eng. 2016, 142 (11), 04016054; DOI 10.1109/mnano.2017.2708818. 62. Del Vento, S.; Dachs, J., Atmospheric occurrence and deposition of polycyclic aromatic hydrocarbons in the northeast tropical and subtropical Atlantic Ocean. Environ. Sci. Technol. 2007, 41 (16), 5608-5613; DOI 10.1021/es0707660. 63. González-Gaya, B. N.; Zúñiga-Rival, J.; Ojeda, M. J.; Jiménez, B. A. Dachs, J., Field measurements of the atmospheric dry deposition fluxes and velocities of polycyclic aromatic hydrocarbons to the global oceans. Environ. Sci. Technol. 2014, 48 (10), 55835592; DOI 10.1021/es500846p. 64. Yang, M. M.; Wang, Y.; Liu, Q.; Ding, A. J.; Li, Y. H., The influence of sandstorms and long-range transport on polycyclic aromatic hydrocarbons (PAHs) in PM2.5 in the highaltitude atmosphere of southern China. Atmosphere 2015, 6 (11), 1633-1651; DOI 10.3390/atmos6111633 65. Eagar, J. D.; Ervens, B.; Herckes, P., Impact of partitioning and oxidative processing of PAH in fogs and clouds on atmospheric lifetimes of PAH. Atmos. Environ. 2017, 160, 132-141; DOI 10.1016/j.atmosenv.2017.04.016. 66. Mandalakis, M.; Gustafsson, Ö.; Alsberg, T.; Egebäck, A. L.; Reddy, C. M.; Xu, L.; Klanova, J.; Holoubek, I.; Stephanou, E. G., Contribution of biomass burning to atmospheric polycyclic aromatic hydrocarbons at three European background sites. Environ. Sci. Technol. 2005, 39 (9), 2976-2982; DOI 10.1021/es048184v. 67. Wilcke, W.; Krauss, M.; Amelung, W., Carbon isotope signature of polycyclic aromatic hydrocarbons (PAHs):  evidence for different sources in tropical and temperate environments? Environ. Sci. Technol. 2002, 36 (16), 3530-3535; DOI 10.1021/es020032h. 68. Okuda, T.; Kumata, H.; Zakaria, M. P.; Naraoka, H.; Ishiwatari, R.; Takada, H., Source identification of Malaysian atmospheric polycyclic aromatic hydrocarbons nearby forest 24

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25

751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780

Environmental Science & Technology

fires using molecular and isotopic compositions. Atmos. Environ. 2002, 36 (4), 611-618; DOI 10.1016/S1352-2310(01)00506-4. 69. Wang, R. W.; Yousaf, B.; Sun, R. Y.; Zhang, H.; Zhang, J. M.; Liu, G. J., Emission characterization and δ13C values of parent PAHs and nitro-PAHs in size-segregated particulate matters from coal-fired power plants. J. Hazard. Mater. 2016, 318, 487-496; DOI 10.1016/j.jhazmat.2016.07.030. 70. Okuda, T.; Kumata, H.; Naraoka, H.; Takada, H., Origin of atmospheric polycyclic aromatic hydrocarbons (PAHs) in Chinese cities solved by compound-specific stable carbon isotopic analyses. Org. Geochem. 2002, 33 (12), 1737-1745; DOI 10.1016/S01466380(02)00180-8. 71. Ahad, J. M.; Jautzy, J. J.; Cumming, B. F.; Das, B.; Laird, K. R.; Sanei, H., Sources of polycyclic aromatic hydrocarbons (PAHs) to northwestern Saskatchewan lakes east of the Athabasca oil sands. Org. Geochem. 2015, 80, 35-45; DOI 10.1016/j.orggeochem.2015.01.001. 72. Slater, G. F.; Benson, A. A.; Marvin, C.; Muir, D., PAH fluxes to Siskiwit revisted: trends in fluxes and sources of pyrogenic PAH and perylene constrained via radiocarbon analysis. Environ. Sci. Technol. 2013, 47 (10), 5066-5073; DOI 10.1021/es400272z. 73. Shen, H. Z.; Huang, Y.; Wang, R.; Zhu, D.; Li, W.; Shen, G. F.; Wang, B.; Zhang, Y. Y.; Chen, Y. C.; Lu, Y., Global atmospheric emissions of polycyclic aromatic hydrocarbons from 1960 to 2008 and future predictions. Environ. Sci. Technol. 2013, 47 (12), 6415-6424; DOI 10.1021/es400857z. 74. Zhang, J.; Liu, G.; Wang, R.; Huang, H., Polycyclic aromatic hydrocarbons in the water-SPM-sediment system from the middle reaches of Huai River, China: distribution, partitioning, origin tracing and ecological risk assessment. Environ. Pollut. 2017, 230, 6171; DOI 10.1016/j.envpol.2017.06.012. 75. Liu, Y.; Chen, L.; Huang, Q. H.; Li, W. Y.; Tang, Y. J.; Zhao, J. F., Source apportionment of polycyclic aromatic hydrocarbons (PAHs) in surface sediments of the Huangpu River, Shanghai, China. Sci. Total Environ. 2009, 407 (8), 2931-2938; DOI 10.1016/j.scitotenv.2008.12.046.

25

ACS Paragon Plus Environment