Influence of Air Pollution on Inhalation and Dermal Exposure of

Mar 8, 2019 - Related Content: Water-soluble Brown Carbon in Atmospheric Aerosols from Godavari (Nepal), A Regional Representative of South Asia...
0 downloads 0 Views 912KB Size
Subscriber access provided by ECU Libraries

Ecotoxicology and Human Environmental Health

Influence of Air Pollution on Inhalation and Dermal Exposure of Human to Organophosphate Flame Retardants: A Case Study During a Prolonged Haze Episode Zhiguo Cao, Leicheng Zhao, Yacai Zhang, Meihui Ren, Yajie Zhang, Xiaotu Liu, Jianye Jie, Zhiyu Wang, Changhe Li, Mohai Shen, and Qingwei Bu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b07053 • Publication Date (Web): 08 Mar 2019 Downloaded from http://pubs.acs.org on March 8, 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 32

Environmental Science & Technology

1

Influence of Air Pollution on Inhalation and Dermal Exposure of Human to

2

Organophosphate Flame Retardants: A Case Study During a Prolonged Haze

3

Episode

4

Zhiguo Cao,*,†,‡,1 Leicheng Zhao,†,1 Yacai Zhang,† Meihui Ren,† Yajie Zhang,† Xiaotu

5

Liu,‡ Jianye Jie,† Zhiyu Wang,† Changhe Li,† Mohai Shen,† and Qingwei Bu*,§

6 7



8

Environment and Pollution Control, Ministry of Education, Henan Key Laboratory for

9

Environmental Pollution Control, Henan Normal University, Xinxiang 453007, China

School of Environment, Key Laboratory for Yellow River and Huai River Water

10



11

Environment, Tsinghua University, Beijing 100084, China

12

§

13

Technology-Beijing, Beijing 100083, China

14

1 Co-first

Beijing Key Laboratory for Emerging Organic Contaminants Control, School of

School of Chemical & Environmental Engineering, China University of Mining &

authors

1

ACS Paragon Plus Environment

Environmental Science & Technology

15

ABSTRACT

16

The health impact of haze is of great concern, but few studies have explored its

17

influence on human inhalation and dermal exposure to trace pollutants. Size-segregated

18

atmospheric particles (n=72) and forehead wipe samples (n=80) from undergraduates

19

were collected in Xinxiang, China, during a prolonged haze episode and analysed for

20

ten organophosphate flame retardants (OPFRs). ∑TCPP and TCEP were the most

21

abundant OPFR substances in all samples. The arithmetic mean particle-bound and

22

forehead OPFR concentrations under a heavy pollution condition (air quality index

23

(AQI), 350-550) were 41.9 ng/m3 (∑8OPFRs) and 7.4 μg/m2 (∑6OPFRs), respectively,

24

apparently greater than the values observed under a light pollution condition (AQI, 60-

25

90) (19.5 ng/m3 and 3.9 μg/m2, respectively). Meteorological conditions played

26

distinctive roles in affecting the OPFR occurrence in atmospheric particles (statistically

27

significant for TCEP and ∑TCPP) and forehead wipes (excluding TPHP), implying that

28

OPFR exposure through inhalation and dermal absorption was synchronously

29

influenced by air quality, and OPFRs on the forehead may be mainly absorbed from the

30

air. Inhalation contributed dominantly to the total OPFR exposure dose for humans

31

when using the relative absorption method to assess dermal exposure, while according

32

to the permeability coefficient method, dermal exposure was much more significant

33

than inhalation. The results of this study indicate that OPFR exposure should attract

34

particular concern in regions with heavy air pollution.

2

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Environmental Science & Technology

35

1.

INTRODUCTION

36

As a substitute for polybrominated diphenyl ethers (PBDEs), organophosphate flame

37

retardants (OPFRs) have been widely used in a range of products, including textiles,

38

furniture, insulation materials, baby products, floor polishes and electronics.1 OPFRs

39

are physically (rather than chemically) mixed into original materials and are thus readily

40

released into the environment by volatilization and abrasion.2 Some OPFRs,

41

particularly those with chlorinated alkyl groups, have low degradation potential and can

42

persist in the environment.3,

43

environmental matrices, and their occurrence, fate, behaviour and consequent human

44

health risk have caused increasing concern in recent years.5-7

4

OPFRs have been detected extensively in multiple

45

OPFRs tend to bind to atmospheric particles and resist metabolic breakdown and

46

photodegradation.8 One of the significant factors determining atmospheric transport of

47

particle-bound contaminants and the associated human health risk is the aerodynamic

48

diameter of atmospheric particles.9 Human inhalation exposure to particle-bound

49

contaminants is also a function of particle size, and smaller particles are transported

50

deep into the human respiratory system and can even enter the bloodstream.10 In recent

51

years, relevant studies have generally focused on the occurrence of OPFR in the total

52

suspended particles in the air,11 but the particle size distribution patterns of OPFRs in

53

the air remain poorly understood.

54

While dust ingestion, diet and inhalation are commonly recognized as the main

55

OPFR exposure pathways for humans,12 several recent studies have illustrated that

56

dermal absorption is also a significant pathway for OPFR exposure.5, 13 In addition, skin 3

ACS Paragon Plus Environment

Environmental Science & Technology

57

wipes have been found to be an effective method to assess the dermal exposure of

58

humans to contaminants.14-16 Xu et al. examined the magnitude of OPFR uptake from

59

dermal absorption using hand wipe,17 and Liu et al. reported the variation of OPFR

60

dermal absorption from palms, backs of hands, and forearms.18

61

To date, investigations on particle-bound OPFRs have been limited to their spatial

62

and temporal variations.19,

20

63

influencing the contamination characteristics of OPFRs in ambient air and on the human

64

skin surface, which has typically been overlooked. Severe haze pollution has frequently

65

occurred over China in recent years, and its health impact has caused increasing

66

concern.21, 22 With different haze levels, the concentrations of particulate matter in the

67

ambient air and meteorological factors, including the height of the atmospheric

68

boundary layer, wind speed, humidity, and temperature, are expected to vary greatly,

69

with differences of up to orders of magnitude.23, 24 It is unknown whether the variations

70

of these parameters would strongly influence the OPFR distribution in ambient air and

71

on the human skin surface. Furthermore, while the air quality around the world differs

72

considerably, few studies have considered the influence of haze (the metric of haze can

73

be PM2.5, PM10 and air quality index (AQI)) when reporting the contamination

74

characteristics of trace organic pollutants in ambient air, thus the comparability of data

75

from different studies is questionable. Specifically, little is known about the variation

76

of OPFR occurrence in ambient air and on human skin under different haze levels. In

77

addition, it remains to be clarified whether OPFRs can accumulate in the air

78

simultaneously with particulate matter during a haze episode and how haze affects

However, air quality may also be a major factor

4

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Environmental Science & Technology

79

human inhalation and dermal exposure to OPFRs.

80

The process and routes of dermal exposure are complex and are currently not well

81

understood. Dermal exposure routes of OPFRs may vary among skin locations.

82

Comparatively, OPFRs on surface of the exposed skin locations may primarily come

83

from the air (absorption of gaseous phase by skin lipid and the deposition of airborne

84

particles).25 OPFRs on surface of the covered skin locations may come from both the

85

air and clothing. Further, sources of OPFRs on hands are more complex.18 Therefore,

86

forehead is a more suitable skin location to reflect the influence of haze on human

87

dermal exposure to OPFRs. To the best of our knowledge, however, no study has

88

reported OPFR levels on the human forehead.

89

The present study addresses these gaps by characterizing OPFRs in size-segregated

90

particulate samples from ambient air, and on human forehead surfaces under light and

91

heavy haze pollution conditions during a prolonged haze episode. The principal

92

objectives of this study are as follows: (1) determine the contamination characteristics

93

of OPFRs in atmospheric particles and on the foreheads of undergraduates in Xinxiang,

94

China; (2) examine the extent to which haze influences OPFR occurrence and the

95

corresponding mechanism; and (3) assess the human inhalation and dermal exposure

96

patterns to OPFRs.

97 98

2.

99

2.1 Sampling information

100

MATERIALS AND METHODS

Atmospheric particles: Particulate samples were collected during a prolonged haze 5

ACS Paragon Plus Environment

Environmental Science & Technology

101

episode from December 2016 to January 2017. The sampling sites were located on the

102

roofs of four university buildings (generally approximately 15 m high) in the Muye

103

District of Xinxiang, China. Seventy-two particulate samples (8 batches with 9 size

104

fractions in each batch) were obtained through two rounds of sampling processes under

105

a light haze pollution condition (with the AQI ranging from 60 to 90) and a heavy haze

106

pollution condition (AQI, 350 to 550). Each sample was collected on a Whatman quartz

107

fibre filter (preheated for 4 h at 450 °C) with a diameter of 81 mm using an Anderson

108

eight-stage cascade impactor (Tisch Environmental Inc., Cleves, OH, USA), and the

109

flow rate was controlled at 28.3 L/min. The cut off aerodynamic diameters for each

110

stage were > 9.0 μm, 5.8-9.0 μm, 4.7-5.8 μm, 3.3-4.7 μm, 2.1-3.3 μm, 1.1-2.1 μm, 0.7-

111

1.1 μm, 0.4-0.7 μm and < 0.4 μm (backup filter), respectively.26 The sampling duration

112

was 48 h for the heavy pollution condition and 120 h for the light pollution condition.

113

After sampling, the filter samples were wrapped in aluminium foil carefully and stored

114

at -20 °C until analysis. Meteorological data, including AQI, PM10, PM2.5, temperature

115

and wind speed were recorded during the sampling process in the Supporting

116

Information (SI, TableS1).

117

Human forehead wipe samples: All samples were collected under normal condition.

118

All of the sampling protocols were approved by the Research Ethics Committee of

119

Henan Normal University, and all participants gave informed consent prior to providing

120

personal information or samples. Similar to former studies, sterile gauze pad (used as

121

skin wipe) obtained from a medicine shop was used as sampling matrix to collect

122

forehead wipe samples. The forehead surface of a specific participant was wiped thrice 6

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

Environmental Science & Technology

123

using one surface of a sterile gauze pad, and then wiped for three additional times using

124

the other side.18, 27 Eighty forehead wipe samples were obtained through two rounds

125

(with the same participants) of sampling processes simultaneously with the collection

126

of atmospheric particles from 20 males and 20 females. All of the 40 participants were

127

recruited from the undergraduate class at the corresponding campus. We informed the

128

participants to keep similar activity patterns between the two rounds of sampling

129

activities and not to wash faces at least for two hours before sample collection. More

130

detailed information in terms of the sample collection and preparation are provided in

131

the SI (Section S2.1).

132 133

2.2 Sample analysis

134

Ten OPFRs were analysed, including triethyl phosphate (TEP), tri-n-propyl

135

phosphate (TPP), tri-isobutyl phosphate (TIBP), tri-n-butyl phosphate (TNBP), tris-

136

(chloropropyl) phosphate isomers (∑TCPP, mixture of 3 isomers), tris-(2-chloroethyl)

137

phosphate (TCEP), tris-(2-butoxyethyl) phosphate (TBOEP), triphenyl phosphate

138

(TPHP), tris-(1,3-dichloro-2-propyl) phosphate (TDCIPP) and tri-cresyl-phosphate

139

isomers (∑TMPP, mixture of 4 isomers). Detailed information on the physical-

140

chemical properties of these compounds and the instrumental analysis procedure18, 27

141

are provided in the SI (Section S2.2 and Table S2).

142 143 144

2.3 Quality assurance and quality control (QA/QC) The spiked recovery of the sample analysis was determined by spiking 10 ng, 50 ng 7

ACS Paragon Plus Environment

Environmental Science & Technology

145

and 100 ng of ten standard OPFRs into a quartz fibre filter (n=5) and a forehead wipe

146

(n=5), respectively, and analysing using the same method as the samples. One

147

procedural blank and one field blank were included for every batch of 8 particulate

148

samples or 8 forehead wipe samples. The recovery of triaryl phosphate (TAP) and

149

triphenyl phosphate-d15 (TPHP-d15) (surrogate standards) for all samples averaged 89

150

± 19% and 93 ± 25% (mean ± SD), respectively. A value of three times of the deviation

151

of all blank samples was used as the method detection limit (MDL) for each compound

152

detected in blank samples. For compounds undetected in the blanks, a signal of ten

153

times of the noise level was used as the MDL (Table S3). More detailed information

154

and data on QA/QC are provided in SI (Section S2.3).

155 156

2.4 Data analysis

157

Statistical analysis was performed using SPSS statistical software package, version

158

17.0 (SPSS Inc.). In the statistical analysis, the concentration below MDL were treated

159

as the half value of MDL. Shapiro-Wilk test was used to examine whether the data were

160

normally or log-normally distributed. All relevant data were log transformed in case

161

their distributions were highly skewed before bivariate comparisons and the level of

162

significance was set to p = 0.05 for all statistical tests. Statistical significance of OPFR

163

concentrations in particulate and forehead wipe samples between light and heavy haze

164

levels was evaluated by a paired-sample t-test. An independent-sample t-test was

165

conducted between males and females to test the statistical significance of OPFR

166

concentration variation in forehead wipes. A Spearman correlation analysis was used 8

ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

Environmental Science & Technology

167

to explore the correlations between the OPFR concentrations in particulate samples

168

(forehead wipe samples) and meteorological parameters.

169 170

3

RESULTS AND DISCUSSION

171

3.1 Concentrations, profiles and size distribution of particle-bound OPFRs

172

The arithmetic mean concentrations of the sum of ∑TCPP, TCEP, TPHP, TEP,

173

TBOEP, TDCIPP, TIBP and TNBP (designated as ∑8OPFRs) in all nine particle

174

fractions were 19.5 ± 1.0 ng/m3 (ranging from 18.5 ng/m3 to 20.7 ng/m3) and 41.9 ± 4.1

175

ng/m3 (ranging from 38.0 ng/m3 to 47.6 ng/m3) in light and heavy pollution conditions,

176

respectively (Figure 1a and Table S4). The OPFR contamination in Xinxiang was much

177

lower than that in the e-waste recycling zone of Qingyuan (∑8OPFRs, 130 ± 130 ng/m3)

178

and Guangzhou, China (∑8OPFRs, 138 ± 127 ng/m3),19 which was expected and

179

consistent with the fact that e-waste recycling zones and metropolitan areas may

180

represent prominent OPFR sources. The OPFR concentration under the light pollution

181

condition was similar to that in the Baoshan District of Shanghai, China, (∑6OPFRs,

182

19.4 ± 10.3 ng/m3)11 but approximately 10 times greater than that in the Great Lakes,

183

USA (∑12OPFRs, 2.1 ± 0.4 ng/m3).28

184

The mean atmospheric concentration of the chlorinated OPFRs (∑ClOPFRs, sum of

185

∑TCPP, TCEP and TDCIPP) was 12.3 ± 0.6 ng/m3 and 34.2 ± 4.4 ng/m3 in light and

186

heavy pollution conditions, respectively, accounting for 63.3% and 81.7% of the

187

∑8OPFRs (Figure S1a), respectively. Similar results were reported by Yang et al.,26

188

Bergh et al.29 and Ren et al.,11 who found that ∑ClOPFRs accounted for 77% (median) 9

ACS Paragon Plus Environment

Environmental Science & Technology

189

of ∑10OPFRs, 94% of ∑9OPFRs, and 69.0% of ∑6OPFRs in the ambient air,

190

respectively. This might be owing to the fact that chlorinated OPFRs had been used

191

more widely than non-halogen OPFRs commercially.26, 30

192

To the best of our knowledge, this study provides the first comparison of OPFR

193

occurrence in atmospheric size-fractionated particles under different air pollution

194

conditions. Ratios of OPFR concentration under a heavy pollution condition to that

195

under light pollution condition were calculated for individual OPFRs at all four

196

sampling sites. Detailed information on the calculation process is described in the SI

197

(Section S3.1). Notably, the concentration ratios were apparently greater for chlorinated

198

OPFRs (ranging from 0.8 to 3.7) than for non-chlorinated OPFRs (ranging from 0.4 to

199

1.7) (Figure 1b), indicating that chlorinated OPFRs can massively accumulate in

200

particulate matter during heavy haze pollution condition. Table S5 shows the Spearman

201

correlation between particle-bound OPFR concentrations and meteorological

202

parameters, including AQI, PM10 (μg/m3), PM2.5 (μg/m3), wind speed (km/h) and

203

temperature (°C). In summary, TCEP and ∑TCPP in particulate samples were

204

positively (p < 0.01) correlated with AQI, PM10 and PM2.5 levels, but presented

205

significant negative correlation (p < 0.01) with wind speed and temperature.

206

Conclusively, the contamination of TCEP and ∑TCPP in atmospheric particles was

207

more significantly (p < 0.05) affected by haze and the related meteorological parameters

208

(Table S6), apparently different from other OPFR constituents. Furthermore, it is

209

speculated that the influence of haze on OPFR occurrence in atmospheric particles (as

210

well as in forehead wipe samples) is compound-specific because of their different 10

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

Environmental Science & Technology

211

emission sources and pathways,8 degradation potential31 and physical-chemical

212

properties.30 More evidence should be obtained to fully understand the influencing

213

mechanism in future.

214

According to the size distribution patterns of particle-bound OPFRs (Figure 1c), it is

215

clear that OPFRs concentrations in all nine fractions were apparently greater in the

216

heavy haze pollution condition than in the light pollution condition. Generally, OPFRs

217

in both haze pollution conditions are highly enriched in the 0.7-1.1 μm and 1.1-2.1 μm

218

fractions, where ∑ClOPFRs and ∑8OPFRs account for 29.4% and 28.3%, respectively,

219

of the total OPFRs in all fractions in the light pollution condition, and 30.4% and 29.0%,

220

respectively, of the total OPFRs in all fractions in the heavy pollution condition (Figure

221

S2). The abundance of most OPFRs in fine (aerodynamic diameter < 2.1 μm) and coarse

222

(aerodynamic diameter > 2.1 μm) particles under both haze pollution conditions were

223

approximately equal (Figure S3), representing high tendency to bind to respirable

224

particles and potential to be exposed by humans for these OPFRs in atmospheric

225

particles. These results are similar to the distribution patterns of atmospheric OPFRs in

226

Guangzhou, China,19 but different from those of atmospheric PBDEs32 and polycyclic

227

aromatic hydrocarbons (PAHs)33 in Guangzhou, China.

228

Mass median aerodynamic diameter (MMAD) is the diameter at which the mass of

229

contaminants in all particle fractions are divided into two equal parts and is calculated

230

for

231

(http://www.mmadcalculator.com/andersen-impactor-mmad.html).

232

information on MMAD is presented in the SI (Section S3.1). MMAD values of each

individual

OPFRs

with

a

11

ACS Paragon Plus Environment

MMAD

calculator More

detailed

Environmental Science & Technology

233

OPFR in light and heavy pollution conditions are shown in Figure 1d. For light

234

pollution condition, TPHP, TCEP, TDCIPP and TEP had MMADs less than 2.5 μm,

235

while the MMADs of other compounds varied from 2.5 μm to 4 μm, except for TIBP

236

(the MMAD of which varied from 4 μm to 5 μm). For the heavy pollution condition,

237

TPHP, TCEP, TBOEP, TIBP and ∑TCPP had MMADs less than 2.5 μm, while the

238

MMADs of TDCIPP, TEP and TNBP varied from 2.5 μm to 4 μm. While the MMADs

239

of both TPHP and TCEP exhibited no differences (less than 2.5 μm) between light and

240

heavy pollution conditions, those of other OPFRs showed apparent differences between

241

light and heavy pollution conditions, indicating for these compounds the particle size

242

distribution patterns changed under the influence of haze level. Since smaller particles

243

can travel deeper into the human respiratory system and potentially enter the

244

bloodstream,34 TPHP, TCEP, TDCIPP and TEP may lead to high human exposure risks

245

in the light pollution condition, and TPHP, TCEP, TBOEP, TIBP and ∑TCPP may lead

246

to high human exposure risks in the heavy pollution condition.

247

The size distribution of atmospheric OPFRs can be correlated to the physiochemical

248

properties of these compounds. A former study reported a significant positive

249

correlation between the MMADs of OPFRs and their vapour pressures.26 However, no

250

significant (p > 0.05) positive correlation between the MMADs of atmospheric OPFRs

251

and their vapour pressures in either light or heavy pollution conditions was observed in

252

the present study (Figure S4). Similarly, negative correlation between the geometric

253

mean diameters (GMDs, a concept similar to MMADs) of OPFRs and their vapour

254

pressures was discovered in the atmosphere at different heights of an e-waste recycling 12

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32

Environmental Science & Technology

255

zone and in urban Guangzhou, China.19 Consequently, in addition to gas-particle

256

partitioning, which depends on the vapour pressure of OPFRs, there might be other

257

sources of particle-bound OPFRs. As speculated in our previous study, abrasion

258

particles or fibres of flame-retarded materials may also contribute to the occurrence of

259

OPFRs in coarse particles containing OPFRs.35

260

The mean normalized size distribution of atmospheric OPFRs in light and heavy haze

261

pollution conditions is shown in Figure 2, indicating there were differences for

262

individual OPFRs in terms of particle size distribution patterns. For the light pollution

263

condition, TNBP, TIBP, TBOEP and ∑TCPP shared a unimodal distribution pattern,

264

with a peak in the 4.7-5.8 μm fraction. TDCIPP, TCEP and TPHP were characterized

265

by a bimodal distribution, with two peaks in the 0.7-1.1 μm and 4.7-5.8 μm fractions.

266

TEP was characterized by a bimodal distribution, with two peaks in the 2.1-3.3 μm and

267

4.7-5.8 μm fractions. For the heavy pollution condition, TEP, TNBP, TBOEP and

268

TDCIPP were characterized by a unimodal distribution pattern, with a peak in the 4.7-

269

5.8 μm fractions. TPHP, TIBP and TCEP shared a bimodal distribution pattern, with

270

two peaks in the 0.7-1.1 μm and 4.7-5.8 μm fractions. ∑TCPP was characterized by a

271

bimodal distribution, with two peaks in the 2.1-3.3 μm and 4.7-5.8 μm fractions.

272

Overall, the size distribution patterns of OPFRs in atmospheric particles cannot be

273

generalized. These differences in size distribution patterns of OPFRs among different

274

studies can be attributed to the physicochemical properties and source characteristics of

275

these compounds, chemical composition and particle size distribution of the

276

atmospheric particles, as well as meteorological conditions. 13

ACS Paragon Plus Environment

Environmental Science & Technology

277 278

3.2 OPFRs in forehead wipe samples, gender variation and implication for dermal

279

exposure route

280

The arithmetic mean concentrations of the sum of ∑TCPP, TCEP, TPHP, TEP,

281

TBOEP and TDCIPP (designated as ∑6OPFRs) in wipe samples were 3.9 ± 0.8 μg/m2

282

(ranging from 2.9 μg/m2 to 5.5 μg/m2) and 7.4 ± 1.6 μg/m2 (ranging from 5.6 μg/m2 to

283

12.4 μg/m2) in light and heavy pollution conditions, respectively (Figure 3a and Table

284

S7). The ∑6OPFR levels in this study were generally less than those on the hands (palms

285

+ backs of hands: 22.7 μg/m2 for ∑3OPFRs), but greater than those on the forearms (3.2

286

μg/m2 for ∑3OPFRs), of participants in Beijing, China.18 Among the OPFRs, ∑TCPP,

287

TCEP and TPHP were the dominant components in both light and heavy pollution

288

conditions, accounting for 36%, 28% and 19% of the ∑6OPFRs, respectively, in the

289

light pollution condition, and 41%, 26% and 15% of the ∑6OPFRs, respectively, in the

290

heavy pollution condition (Figure S1b). The OPFR profiles on human foreheads in this

291

study were similar to those observed on human hands18 and in atmospheric particles, as

292

reported by former studies,11, 26 with ∑TCPP, TCEP and TPHP as the most abundant

293

analogues, indicating wide usage of these compounds in China. Limited literature

294

reported that the consumption amount of OPFRs in China reached approximately 70000

295

tons in 2007, and the marketing demand was expected to grow with an annual rate of

296

15%.8, 36 Globally, the usage of OPFRs was 500000 tons in 2011 and was expected to

297

reach 680000 tons in 2015.30 With such intensive application, it is suggested that human

298

exposure to OPFRs should cause more concern in the near future. 14

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32

Environmental Science & Technology

299

There were no significant differences for OPFRs (except for TPHP) on foreheads

300

between male and female participants (p > 0.05), but the p-values of TPHP in light (p

301

= 0.04) and heavy (p = 0.03) pollution conditions suggested potential gender-based

302

differences (Table S8). Former studies have investigated OPFR levels on the skin

303

surface without considering the air quality of the sampling sites over time.5, 17, 18 In the

304

present study, the mean concentrations of ∑TCPP, TCEP, TEP, TBOEP and TDCIPP

305

in forehead wipes in heavy pollution condition were all significantly (p < 0.01) greater

306

than those in the light pollution condition, but the p-value (0.06) for TPHP suggested

307

that the influence of haze on dermal exposure to TPHP was not significant (Figure 3b

308

and Table S9). As shown in Table S5, similar with the results for chlorinated OPFRs in

309

particulate samples, most OPFRs (except for TPHP) in forehead wipe samples were

310

positively correlated (p < 0.01) with AQI, PM10 and PM2.5 levels, but wind speed and

311

temperature exhibited negative correlations (p < 0.01) with concentrations of OPFRs

312

except for TPHP. The OPFR profiles in forehead wipe samples were similar with that

313

in particulate samples in the present study, with ∑TCPP and TCEP as the major

314

contributor, accounting for more than 56% and 64% of the total OPFRs in particulate

315

samples and forehead wipe samples, respectively (Figure S1). Nevertheless, TPHP

316

contributed less than 2% to the total OPFRs in atmospheric particles but more than 15%

317

to the total OPFRs in forehead wipes. Nail polish is a specific TPHP source as we

318

know.13 As semi-volatile organic compound, TPHP can volatilize into the ambient air

319

from personal products, and we deduce that there might be higher level of TPHP in the

320

air of females’ dormitories than males’ living environment. Then females’ exposure to 15

ACS Paragon Plus Environment

Environmental Science & Technology

321

TPHP might be relatively higher accordingly.37 Further, females’ fingers may contact

322

with forehead on occasion while tidying their hair. Therefore, the application of

323

personal care products (e.g., nail polishes) might be the possible reason for the

324

significant differences for TPHP between genders. Because human forehead directly

325

contacts with the air, and our results indicated that OPFR levels in atmospheric particles

326

and in forehead wipes were synchronously influenced by air quality, it could be

327

concluded that OPFRs on forehead might mainly origin from the air.

328 329

3.3 OPFR exposure doses and comparison between inhalation and dermal

330

absorption pathways

331

Particle size plays a crucial role in human exposure assessment because inhalation

332

exposure to particle-bound OPFRs is a function of particle size. The data on the particle

333

size distributions of OPFRs were used to calculate OPFR fractions that could be inhaled

334

through the nose or mouth (inhalable fraction, IF, %), penetrate progressively into the

335

lung below the larynx (thoracic fraction, TF, %) and even reach deep into the gas

336

exchange region (respirable fraction, RF, %).38 The inhaled particles generally

337

accumulate in three main regions of the human respiratory tract: head airway (HA),

338

tracheobronchial region (TB) and alveolar region (AR). Deposition concentration (DC,

339

ng/m3) is defined as the OPFR concentration potentially deposited in a specific region

340

of the human respiratory tract. Furthermore, the deposition fraction (DF, %) is the ratio

341

of DCTotal (sum of DCHA, DCTB and DCAR) to CTotal (total concentrations of OPFRs in

342

all size-fractionated particles in the air). The daily inhalation dose (DID, ng/kg BW/day) 16

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Environmental Science & Technology

343

is estimated based on DCTotal (compared with CTotal) and daily inhalation rate (m3/day)

344

of human.19,

345

provided in the SI (Section S3.3).

26, 39

Detailed information on the equations and parameter setting are

346

Overall, the values of inhalable fraction, thoracic fraction, respirable fraction and

347

deposition fraction for individual OPFRs were estimated to be 89.5-93.0%, 83.3-89.7%,

348

50.9-72.6% and 51.9-70.1%, respectively, in light pollution condition, and 89.3-93.7%,

349

84.3-91.3%, 52.0-77.5% and 54.0-69.2%, respectively, in heavy pollution condition

350

(Table S10). The ratios of deposition fraction to inhalable fraction for individual OPFRs

351

ranged from 55.7% to 78.7%, with a mean value of 69.4% in light pollution condition,

352

and from 57.6% to 77.5%, with a mean value of 67.9% in heavy pollution condition

353

(Figure S5).

354

Specifically, the deposition concentrations and fractions of size-fractionated OPFRs

355

in the head airway, tracheobronchial region and alveolar region in light and heavy

356

pollution conditions are shown in Table S11 and S12. The deposition fractions of

357

individual OPFRs in HA, TB and AR of the human respiratory tract were 40.8-59.6%,

358

2.7-3.7% and 7.0-8.7%, respectively, in light pollution condition, and 41.6-58.8%, 3.1-

359

3.6% and 6.9-9.3%, respectively, in heavy pollution condition (Figure 4a and Table

360

S12). The majority of OPFRs were deposited in the head airway (DFHA), averagely

361

accounting for 82.1% and 81.5% of DFTotal (sum of DFHA, DFTB and DFAR) in light and

362

heavy pollution conditions, respectively. Apparently, the deposition fraction of OPFRs

363

in different human respiratory tracts was particle size dependent (Figure 4b and Table

364

S13). With increasing particle size, the deposition concentration of OPFRs increased in 17

ACS Paragon Plus Environment

Environmental Science & Technology

365

the HA, remained steady in the TB and decreased in the AR (Figure 4b). Fine particles

366

(particle size < 2.1 μm) were dominant in the alveolar region (60.1% for light pollution

367

condition and 59.7% for heavy pollution condition), while coarse particles (particle size

368

> 2.1 μm) contributed the most in the head airway (77.7% for the light pollution

369

condition and 77.4% for the heavy pollution condition) (Table S13).

370

For adults, Figure S6a and Table S14 depict the total mean DIDDC values (based on

371

DCTotal) of individual particle-bound OPFRs in light and heavy pollution conditions.

372

Chlorinated OPFRs were the main contributors to the total mean DIDDC in both light

373

and heavy pollution conditions. Specifically, the total mean DIDC values of ∑8OPFRs

374

(based on CTotal) in light and heavy pollution conditions were 4.2 ng/kg BW/day and

375

9.0 ng/kg BW/day, respectively, both 1.6 times greater than DIDDC values (2.6 ng/kg

376

BW/day and 5.6 ng/kg BW/day, respectively) (Figure S6a), suggesting apparent

377

overestimations on daily intake dose with bulk OPFR concentrations. Furthermore,

378

DIDDC and DIDC values in heavy pollution condition were 2.1 and 2.2 times greater,

379

respectively, than those in light pollution condition. Comparatively, DIDDC and DIDC

380

values of individual OPFRs for children were all 2.4 times greater than that for adults

381

in both light and heavy pollution conditions (Table S14).

382

Permeability coefficient (Kp) and relative absorption (RA) methods were conducted

383

to calculate daily dermal absorption doses (DAD) of OPFRs in the present study. The

384

Kp method was developed from Fick’s first law of diffusion,40 and the RA method was

385

developed based on the US EPA Exposure Factors Handbook.41, 42 Detailed information

386

on the equations and parameter setting are provided in the SI (Section S3.3 and Table 18

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

387

Environmental Science & Technology

S15).

388

For adults, the DAD values for human foreheads, according to the Kp and RA

389

methods, in light and heavy pollution conditions are shown in Figure S6a and Table

390

S16. Chlorinated OPFRs were the dominant contributors to the total mean DAD in both

391

light and heavy pollution conditions. The DADKp and DADRA values of ∑6OPFRs in

392

heavy pollution condition were 1.8 and 1.9 times greater than those in light pollution

393

condition, respectively. Assuming the OPFR contamination on foreheads of children

394

and adults are similar, it is deduced that DAD values of individual OPFRs for children

395

might be 2.1 times greater than that for adults (Figure S7a and Table S16) in both light

396

and heavy pollution conditions, either with Kp method or RA method.

397

Furthermore, according to the DID and DAD data calculated above, the contributions

398

of inhalation and dermal absorption (forehead) to human exposure to OPFRs were

399

compared. DAD values based on the Kp method were approximately 1.5 and 1.2 times

400

greater than DIDDC values in light and heavy pollution conditions, respectively,

401

however, DAD values based on the RA method were approximately 9.3 and 10.3 times

402

less than DIDDC values in light and heavy pollution conditions, respectively (Figure

403

S6a).

404 405

3.4 Health risk assessment

406

The non-cancer risk from inhalation and dermal exposure to OPFRs can be evaluated

407

by the hazard quotient (HQ).18, 43 We used the oral reference dose (RfD) of OPFRs to

408

estimate the health risk because RfD data of OPFRs for inhalation and dermal exposure 19

ACS Paragon Plus Environment

Environmental Science & Technology

409

are scarce. Detailed information on the equations and parameters are provided in the SI

410

(Section S3.4 and Table S2).

411

For inhalation exposure of adults, the total mean HQC values of ∑8OPFRs (based on

412

CTotal) in light and heavy pollution conditions were 3.3×10-4 and 9.0×10-4, respectively,

413

which were 1.6 and 1.6 times greater, respectively, than HQDC values (based on

414

deposition concentration, 2.1×10-4 and 5.4×10-4 in light and heavy pollution conditions)

415

(Figure S6b). Furthermore, HQDC and HQC values of ∑8OPFRs for adults in heavy

416

pollution condition were 2.6 and 2.7 times greater, respectively, than those in the light

417

pollution condition. For children, the HQDC and HQC values of individual OPFRs were

418

all 2.4 times greater than that for adults in both light and heavy pollution conditions

419

(Figure S7b and Table S14).

420

For adults’ dermal exposure, the total mean HQKp values for ∑6OPFRs on human

421

forehead (based on the Kp method) in light and heavy pollution conditions were 5.4×10-

422

4

423

calculated by RA method (HQRA: 3.1×10-5 and 5.9×10-5, respectively) (Figure S6b).

424

HQKp and HQRA values in heavy pollution condition were 1.8 and 1.9 times greater,

425

respectively, than those in the light pollution condition. For children, HQKp and HQRA

426

values of individual OPFRs were all 2.1 times greater than that for adults in both light

427

and heavy pollution conditions (Figure S7b and Table S16).

and 9.5×10-4, respectively, which were 17.4 and 16.1 times greater than those

428

Furthermore, for adults, the HQDC of particle-bound OPFRs was compared with the

429

HQRA and HQKp of forehead wipe OPFRs to evaluate the contribution of inhalation and

430

dermal absorption to human health risk (Figure S6c). Specifically, relative to HQRA, 20

ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

Environmental Science & Technology

431

HQDC was the major contributor to total HQ in light pollution (87.0%) and heavy

432

pollution (90.3%) conditions, indicating inhalation to be the dominant exposure

433

pathway. Nevertheless, HQKp contributed substantially to the total HQ relative to HQDC,

434

accounting for 72.6% of the total risk in the light pollution condition and 63.6% of the

435

total risk in the heavy pollution conditions (Figure S6c).

436 437

3.5 Implications for future studies

438

Influence of haze on human inhalation and dermal exposure to OPFRs was firstly

439

synchronously reported in this study. It is clear that air quality is a crucial factor

440

influencing both inhalation and dermal exposure of humans to OPFRs, and OPFR

441

contamination on human forehead might substantially depend on their distribution in

442

the ambient air, especially for chlorinated OPFRs. As values of meteorological

443

parameters can vary significantly (especially in China) under different sampling

444

conditions, it is recommended that air quality should be considered and labelled when

445

investigating organic contamination in the atmosphere and on human skin surfaces.

446

Without considering air quality, the results and conclusions of relevant studies may

447

suffer from prominent contingencies and poor universality. In addition, in regions and

448

countries with severe air pollution, human exposure to contaminants through inhalation

449

and dermal absorption could be of particular concern. The correlation between OPFR

450

pollution in the air and human dermal exposure, as well as the corresponding

451

mechanism are in need to be explored and verified in the near future.

452 21

ACS Paragon Plus Environment

Environmental Science & Technology

453

■ ASSOCIATED CONTENT

454

Supporting Information: The Supporting Information is available free of charge on

455

the ACS Publications website at DOI:

456 457

The SI contains additional detailed information on chemicals, materials, and instrumental analysis and some additional tables and figures, as noted in the text.

458 459

■ AUTHOR INFORMATION

460

Corresponding author

461

*Phone: +86-373-3325971, Fax: +86-373-3325971. E-mail: [email protected]

462

(Zhiguo Cao)

463

* Phone: +86-10-62339298, Fax: +86-10-62339298. E-mail:

464

[email protected] (Qingwei Bu)

465

ORCID

466

Zhiguo Cao: 0000-0002-8580-3368

467

Leicheng Zhao: 0000-0002-9846-5126

468

Xiaotu Liu: 0000-0002-3698-019X

469

Author contributions

470

Zhiguo Cao and Leicheng Zhao contributed equally to the study and should be regarded

471

as joint first authors.

472

Notes

473

The authors declare no competing financial interest.

474

Ethics 22

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

Environmental Science & Technology

475

Ethical approval for this investigation was obtained from the Research Ethics

476

Committee of Henan Normal University.

477 478

ACKNOWLEDGEMENTS

479

This work was supported by National Natural Science Foundation of China (21607038,

480

21806030), China Postdoctoral Science Foundation (2016T90668, 2015M570629),

481

Open Fund of State Key Laboratory of Environmental Chemistry and Ecotoxicology,

482

Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences

483

(KF2015-09, KF2017-01), Science Foundation (2016PL14, 20180572) of Henan

484

Normal University and Key Scientific Research Project Plan of Henan Province

485

(16A610002, 17A610007). We thank the participants for their active involvement and

486

support for this study.

23

ACS Paragon Plus Environment

Environmental Science & Technology

487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529

REFERENCES (1) Li, J.; Tang, J.; Mi, W.; Tian, C.; Emeis, K.-C.; Ebinghaus, R.; Xie, Z. Spatial Distribution and Seasonal Variation of Organophosphate Esters in Air above the Bohai and Yellow Seas, China. Environ Sci Technol. 2018, 52 (1), 89-97. (2) Castro-Jimenez, J.; Gonzalez-Gaya, B.; Pizarro, M.; Casal, P.; Pizarro-Alvarez, C.; Dachs, J. Organophosphate Ester Flame Retardants and Plasticizers in the Global Oceanic Atmosphere. Environ Sci Technol. 2016, 50 (23), 12831-12839. (3) van den Eede, N.; Dirtu, A. C.; Neels, H.; Covaci, A. Analytical developments and preliminary assessment of human exposure to organophosphate flame retardants from indoor dust. Environment International. 2011, 37 (2), 454-461. (4) Liu, Y.; Liggio, J.; Harner, T.; Jantunen, L.; Shoeib, M.; Li, S.-M. Heterogeneous OH Initiated Oxidation: A Possible Explanation for the Persistence of Organophosphate Flame Retardants in Air. Environ Sci Technol. 2014, 48 (2), 1041-1048. (5) Hoffman, K.; Garantziotis, S.; Birnbaum, L. S.; Stapleton, H. M. Monitoring Indoor Exposure to Organophosphate Flame Retardants: Hand Wipes and House Dust. Environmental Health Perspectives. 2015, 123 (2), 160-165. (6) Hoffman, K.; Butt, C. M.; Chen, A.; Limkakeng, A. T., Jr.; Stapleton, H. M. High Exposure to Organophosphate Flame Retardants in Infants: Associations with Baby Products. Environ Sci Technol. 2015, 49 (24), 14554-14559. (7) Salamova, A.; Hermanson, M. H.; Hites, R. A. Organophosphate and halogenated flame retardants in atmospheric particles from a European Arctic site. Environ Sci Technol. 2014, 48 (11), 6133-6140. (8) Wei, G. L.; Li, D. Q.; Zhuo, M. N.; Liao, Y. S.; Xie, Z. Y.; Guo, T. L.; Li, J. J.; Zhang, S. Y.; Liang, Z. Q. Organophosphorus flame retardants and plasticizers: Sources, occurrence, toxicity and human exposure. Environ Pollut. 2015, 196, 29-46. (9) Luo, P.; Ni, H. G.; Bao, L. J.; Li, S. M.; Zeng, E. Y. Size distribution of airborne particle-bound polybrominated diphenyl ethers and its implications for dry and wet deposition. Environ Sci Technol. 2014, 48 (23), 13793-9. (10) Lv, Y.; Li, X.; Xu, T. T.; Cheng, T. T.; Yang, X.; Chen, J. M.; Iinuma, Y.; Herrmann, H. Size distributions of polycyclic aromatic hydrocarbons in urban atmosphere: sorption mechanism and source contributions to respiratory deposition. Atmospheric Chemistry and Physics. 2016, 16 (5), 2971-2983. (11) Ren, G.; Chen, Z.; Feng, J.; Ji, W.; Zhang, J.; Zheng, K.; Yu, Z.; Zeng, X. Organophosphate esters in total suspended particulates of an urban city in East China. Chemosphere. 2016, 164, 75-83. (12) Vykoukalova, M.; Venier, M.; Vojta, S.; Melymuk, L.; Becanova, J.; Romanak, K.; Prokes, R.; Okeme, J. O.; Saini, A.; Diamond, M. L.; Klanova, J. Organophosphate esters flame retardants in the indoor environment. Environment International. 2017, 106, 97-104. (13) Mendelsohn, E.; Hagopian, A.; Hoffman, K.; Butt, C. M.; Lorenzo, A.; Congleton, J.; Webster, T. F.; Stapleton, H. M. Nail polish as a source of exposure to triphenyl phosphate. Environment International. 2016, 86, 45-51. (14) US.EPA. Dermal Exposure Assessment: Principles and Applications. U.S. Environmental Protection Agency, Washington, DC (EPA/600/8-91/011B). 1998. (15) Stapleton, H. M.; Eagle, S.; Sjoedin, A.; Webster, T. F. Serum PBDEs in a North Carolina Toddler Cohort: Associations with Handwipes, House Dust, and Socioeconomic Variables. Environmental Health Perspectives. 2012, 120 (7), 1049-1054. 24

ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

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 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573

Environmental Science & Technology

(16) Hammel, S. C.; Hoffman, K.; Webster, T. F.; Anderson, K. A.; Stapleton, H. M. Measuring Personal Exposure to Organophosphate Flame Retardants Using Silicone Wristbands and Hand Wipes. Environ Sci Technol. 2016, 50 (8), 4483-4491. (17) Xu, F.; Giovanoulis, G.; van Waes, S.; Padilla-Sanchez, J. A.; Papadopoulou, E.; Magner, J.; Haug, L. S.; Neels, H.; Covaci, A. Comprehensive Study of Human External Exposure to Organophosphate Flame Retardants via Air, Dust, and Hand Wipes: The Importance of Sampling and Assessment Strategy. Environ Sci Technol. 2016, 50 (14), 7752-7760. (18) Liu, X.; Yu, G.; Cao, Z.; Wang, B.; Huang, J.; Deng, S.; Wang, Y. Occurrence of organophosphorus flame retardants on skin wipes: Insight into human exposure from dermal absorption. Environment International. 2017, 98, 113-119. (19) Luo, P.; Bao, L.-J.; Guo, Y.; Li, S.-M.; Zeng, E. Y. Size-dependent atmospheric deposition and inhalation exposure of particle-bound organophosphate flame retardants. Journal of Hazardous Materials. 2016, 301, 504-511. (20) Zhou, L.; Hiltscher, M.; Gruber, D.; Puettmann, W. Organophosphate flame retardants (OPFRs) in indoor and outdoor air in the Rhine/Main area, Germany: comparison of concentrations and distribution profiles in different microenvironments. Environmental Science and Pollution Research. 2017, 24 (12), 10992-11005. (21) Han, D.; Wang, Z.; Cheng, J.; Wang, Q.; Chen, X.; Wang, H. Volatile organic compounds (VOCs) during non-haze and haze days in Shanghai: characterization and secondary organic aerosol (SOA) formation. Environmental Science and Pollution Research. 2017, 24 (22), 18619-18629. (22) Zhuang, X.; Wang, Y.; He, H.; Liu, J.; Wang, X.; Zhu, T.; Ge, M.; Zhou, J.; Tang, G.; Ma, J. Haze insights and mitigation in China: An overview. Journal of Environmental Sciences. 2014, 26 (1), 2-12. (23) Tang, G.; Zhang, J.; Zhu, X.; Song, T.; Muenkel, C.; Hu, B.; Schaefer, K.; Liu, Z.; Zhang, J.; Wang, L.; Xin, J.; Suppan, P.; Wang, Y. Mixing layer height and its implications for air pollution over Beijing, China. Atmospheric Chemistry and Physics. 2016, 16 (4), 2459-2475. (24) Zhu, Q.; Liu, G.; Zheng, M.; Zhang, X.; Gao, L.; Sua, G.; Liang, Y. Size distribution and sorption of polychlorinated biphenyls during haze episodes. Atmospheric Environment. 2018, 173, 38-45. (25) Gong, M.; Zhang, Y.; Weschler, C. J. Measurement of Phthalates in Skin Wipes: Estimating Exposure from Dermal Absorption. Environ Sci Technol. 2014, 48 (13), 7428-7435. (26) Yang, F.; Ding, J.; Huang, W.; Xie, W.; Liu, W. Particle size-specific distributions and preliminary exposure assessments of organophosphate flame retardants in office air particulate matter. Environ Sci Technol. 2014, 48 (1), 63-70. (27) Liu, X.; Cao, Z.; Yu, G.; Wu, M.; Li, X.; Zhang, Y.; Wang, B.; Huang, J. Estimation of exposure to organic flame retardants via hand wipe, surface wipe, and dust: Comparability of different assessment strategies. Environ Sci Technol. 2018, 52 (17), 9946-9953. (28) Salamova, A.; Ma, Y.; Venier, M.; Hites, R. A. High Levels of Organophosphate Flame Retardants in the Great Lakes Atmosphere. Environmental Science & Technology Letters. 2014, 1 (1), 8-14. (29) Bergh, C.; Torgrip, R.; Emenius, G.; Ostman, C. Organophosphate and phthalate esters in air and settled dust - a multi-location indoor study. Indoor Air. 2011, 21 (1), 67-76. (30) van der Veen, I.; de Boer, J. Phosphorus flame retardants: Properties, production, environmental occurrence, toxicity and analysis. Chemosphere. 2012, 88 (10), 1119-1153. (31) Wang, T.; Tian, M.; Ding, N.; Yan, X.; Chen, S. J.; Mo, Y. Z.; Yang, W. Q.; Bi, X. H.; Wang, X. M.; Mai, B. X. Semivolatile Organic Compounds (SOCs) in Fine Particulate Matter (PM2.5) during Clear, Fog, and Haze Episodes in Winter in Beijing, China. Environ Sci Technol. 2018, 52 (9), 519925

ACS Paragon Plus Environment

Environmental Science & Technology

574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603

5207. (32) Zhang, B. Z.; Zhang, K.; Li, S. M.; Wong, C. S.; Zeng, E. Y. Size-dependent dry deposition of airborne polybrominated diphenyl ethers in urban Guangzhou, China. Environ Sci Technol. 2012, 46 (13), 7207-14. (33) Zhang, K.; Zhang, B.-Z.; Li, S.-M.; Zhang, L.-M.; Staebler, R.; Zeng, E. Y. Diurnal and seasonal variability in size-dependent atmospheric deposition fluxes of polycyclic aromatic hydrocarbons in an urban center. Atmospheric Environment. 2012, 57, 41-48. (34) Ni, H.-G.; Cao, S.-P.; Chang, W.-J.; Zeng, H. Incidence of polybrominated diphenyl ethers in central air conditioner filter dust from a new office building. Environ Pollut. 2011, 159 (7), 1957-1962. (35) Cao, Z.; Xu, F.; Covaci, A.; Wu, M.; Wang, H.; Yu, G.; Wang, B.; Deng, S.; Huang, J.; Wang, X. Distribution Patterns of Brominated, Chlorinated, and Phosphorus Flame Retardants with Particle Size in Indoor and Outdoor Dust and Implications for Human Exposure. Environ Sci Technol. 2014, 48 (15), 8839-8846. (36) Wang, X.; Liu, J.; Yin, Y. The pollution status and research progress on organophosphate ester flame retardants. Prog. Chem. (in Chinese). 2010, 22, 1983-1992. (37) Hong, X.; Chen, R.; Hou, R.; Yuan, L.; Zha, J. Triphenyl Phosphate (TPHP)-Induced Neurotoxicity in Adult Male Chinese Rare Minnows ( Gobiocypris rarus). Environ Sci Technol. 2018. (38) Hinds, W. C., Aerosol Technology: Properties, Behavior, and Measurement of Airborne Particles, 2nd Edition. Wiley Interscience, New York, 1999. (39) ICRP. International Commission on Radiological Protection; Publication 66: Human respiratory tract model for radiological protection. Ann. ICRP. 1994, 24 (1-3). (40) Weschler, C. J.; Nazaroff, W. W. SVOC exposure indoors: fresh look at dermal pathways. Indoor Air. 2012, 22 (5), 356-377. (41) US.EPA. Exposure Factors Handbook Chapter 7. U.S. Environmental Protection Agency. 2011. (42) Cao, Z.; Zhao, L.; Zhu, G.; Chen, Q.; Yan, G.; Zhang, X.; Wang, S.; Wu, P.; Sun, L.; Shen, M.; Zhang, S. Propositional modification for the USEPA models for human exposure assessment on chemicals in settled dust or soil. Environmental Science and Pollution Research. 2017, 24 (24), 2011320116. (43) US.EPA, United States Environmental Protection Agency, Risk assessment guidance for superfund. Volume I: human health evaluation manual (Part A). In Washington, D.C,. In 1989.

26

ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

Environmental Science & Technology

605

Figure captions

606

Figure 1. Concentrations, size distribution and MMAD (mass median aerodynamic

607

diameter of OPFRs) of OPFRs in light and heavy pollution conditions in particulate

608

samples.

609 610

Figure 2. Mean normalized size distribution of particle-bound OPFR analogues in light

611

and heavy pollution conditions.

612 613

Figure 3. OPFR concentrations and ratios between heavy and light pollution conditions

614

in forehead wipe samples.

615 616

Figure 4. (a): Deposition fraction (DF) of size-fractionated OPFRs in the head airway

617

(HA), tracheobronchial region (TB) and alveolar region (AR) of the human respiratory

618

tract in light and heavy haze pollution conditions; (b): Total mean contribution of size-

619

fractionated OPFR deposition concentrations (DC) in the three areas of human

620

respiratory tracts. ∑ClOPFRs is the sum of TDCIPP, TCEP and ∑TCPP; DFTotal is the

621

sum of DFHA, DFTB and DFAR.

27

ACS Paragon Plus Environment

Environmental Science & Technology

622 623

Figure 1

28

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Environmental Science & Technology

624 625

Figure 2

29

ACS Paragon Plus Environment

Environmental Science & Technology

626 627

Figure 3

30

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

Environmental Science & Technology

628 629

Figure 4

31

ACS Paragon Plus Environment

Environmental Science & Technology

630

TOC Art

631

32

ACS Paragon Plus Environment

Page 32 of 32