Simulating and Explaining Passive Air Sampling Rates for

Jul 9, 2013 - Here we present a method for modeling hourly sampling rates for ... Citation data is made available by participants in Crossref's Cited-...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Glasgow Library

Article

Simulating and explaining passive air sampling rates for semivolatile compounds on polyurethane foam passive samplers Nicholas T Petrich, Scott Spak, Greg R Carmichael, Dingfei Hu, Andres Martinez, and Keri C. Hornbuckle Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es401532q • Publication Date (Web): 09 Jul 2013 Downloaded from http://pubs.acs.org on July 14, 2013

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24

Environmental Science & Technology

1

Simulating and explaining passive air sampling rates for semi-volatile

2

compounds on polyurethane foam passive samplers

3 4

Nicholas T. Petrich1,2, Scott N. Spak1,2,3*, Gregory R. Carmichael1,2,4, Dingfei Hu1,5, Andres

5

Martinez1,5, Keri C. Hornbuckle1,2,5*

6 7 8 9 10 11 12 13 14 15 16

1. Department of Civil & Environmental Engineering, The University of Iowa, Iowa City, Iowa, United States, 52242 2. Center for Global and Regional Environmental Research, The University of Iowa, Iowa City, Iowa, United States, 52242 3. Public Policy Center and School of Urban & Regional Planning, The University of Iowa, Iowa City, Iowa, United States, 52242 4. Department of Chemical & Biochemical Engineering, The University of Iowa, Iowa City, Iowa, United States, 52242 5. IIHR-Hydroscience and Engineering, The University of Iowa, Iowa City, Iowa, United States, 52242

17 18

*Phone: (319) 335-9993 (S.N.S); (319) 384-0789 (K.C.H.). Fax: (319) 335-6801 (S.N.S.); (319)

19

335-5660 (K.C.H.). E-mail: [email protected] (S.N.S.); [email protected]

20

(K.C.H.).

21 22

Keywords: Great Lakes, PCB concentrations, Chicago, persistent organic pollutants

23

1

ACS Paragon Plus Environment

Environmental Science & Technology

24

Abstract

25

Passive air samplers (PAS) including polyurethane foam (PUF) are widely deployed as an

26

inexpensive and practical way to sample semi-volatile pollutants. However, concentration

27

estimates from PAS rely on constant empirical mass transfer rates, which add unquantified

28

uncertainties to concentrations. Here we present a method for modeling hourly sampling rates

29

for semi-volatile compounds from hourly meteorology using first-principle chemistry, physics,

30

and fluid dynamics, calibrated from depuration experiments. This approach quantifies and

31

explains observed effects of meteorology on variability in compound-specific sampling rates and

32

analyte concentrations; simulates nonlinear PUF uptake; and recovers synthetic hourly

33

concentrations at a reference temperature. Sampling rates are evaluated for polychlorinated

34

biphenyl congeners at a network of Harner model samplers in Chicago, Illinois during 2008,

35

finding simulated average sampling rates within analytical uncertainty of those determined from

36

loss of depuration compounds, and confirming quasi-linear uptake. Results indicate hourly, daily

37

and interannual variability in sampling rates, sensitivity to temporal resolution in meteorology,

38

and predictable volatility-based relationships between congeners. We quantify importance of

39

each simulated process to sampling rates and mass transfer and assess uncertainty contributed

40

by advection, molecular diffusion, volatilization, and flow regime within the PAS, finding PAS

41

chamber temperature contributes the greatest variability to total process uncertainty (7.3%).

42

TOC/Abstract Art

43

44 2

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24

45 46

Environmental Science & Technology

Introduction Passive air samplers (PAS) collect atmospheric trace compounds under ambient

47

conditions onto polyurethane foam (PUF), semi-permeable membranes, and polymer resin

48

sampling media.1-8 This versatile, low-cost technique proves useful for establishing spatial

49

gradients in air toxics and emerging contaminants of concern, and in sampling semi-volatile

50

organic compounds (SVOCs) for which real-time observations remain challenging. Many

51

SVOCs, including polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons, are

52

known or potential human carcinogens, toxic, bio-accumulative, and subject to long-range

53

chemical transport.2 Passive air sampling is also routinely used in atmospheric monitoring for

54

mercury, ammonia, aerosols, and nitrogen dioxide.9, 10

55

In passive air sampling, ambient air flows through the PAS chamber and the

56

compound’s mass accumulates onto sampling media throughout the sampler’s deployment

57

period. Without the controlled flow rates of active methods, concentrations must be calculated

58

from estimated passive sampling rates. Traditionally, PAS has been described as the uptake of

59

a chemical onto sampling media (e.g. polyurethane foam7) as a function of sampling rate (Rs)

60

and the fugacity gradient between the ambient air and the PUF disk, with uptake onto the disk

61

defined as 1, 4, 5, 7:

62

 = ∙ ∙(  



  )

(1)



63

where  represents the mass of a gaseous or semi-volatile compound on the PUF (ng), 

64

the empirical constant mass transfer coefficient (m s-1), the surface area of the PUF (m2), 

65

the analyte concentration in the air (ng m-3),  the compound’s concentration on the PUF (ng

66

m-3), and  the PUF:air equilibrium partition coefficient (m3 g-1). The product of  x is RS,

67

which may be determined experimentally from loss of depuration compounds or uptake of native

68

compounds.4 For RS estimation by loss of depuration compounds, the PUF disk is spiked with a

69

known amount of chemical prior to exposure and then the mass transfer coefficient is calculated 3

ACS Paragon Plus Environment

Environmental Science & Technology

70

based on its first order elimination rate during the entire deployment.5, 11 Alternatively, the

71

uptake of native compounds entails deploying a series of clean PUF disks in a setting with a

72

known ambient concentration and deriving a linear uptake curve to calculate RS.5, 12 These two approaches in determining RS lead to a range of issues regarding

73 74

concentration observations gathered using PAS. As with the analyzed mass from a PAS, these

75

techniques result in only one RS value for the entire sampling period, with unquantified

76

uncertainties and unknown variability. Beyond the large quantified and unquantified errors in the

77

analytical approach to calculating RS for individual samples, empirical RS values are often used

78

for different periods and different locations from those in which they were originally determined,

79

leading to even greater uncertainties in resultant analyte air concentrations.

80

Empirical and theoretical studies suggest that much of the uncertainty stems from the

81

effects of meteorology on RS, primarily with a focus on isolating how individual meteorological

82

parameters such as wind speed11 and temperature13 affect observed variability. Klánová et al.

83

found that RS for gas-phase compounds, including PCBs, exhibits a linear correlation to wind

84

speed and inverse correlation to temperature.14 The angle at which the ambient air enters the

85

PAS chamber also affects flow within the chamber, with positive angles of 15° and 30°

86

increasing RS by 30% and 50%, respectively, and a -10° angle leading to a 40% reduction.15

87

While studies have attempted to address these issues through computational fluid dynamic

88

modeling,15, 16 and theoretical models of the passive sampling uptake process have been

89

developed,1 no comprehensive process model to date has directly simulated and quantified the

90

comprehensive effects of atmospheric dynamics on observed passive sampling flow rates. The PCB concentration in air can be experimentally determined from Equation 2:

91

  ∗

92

 

93

where MPCB is the total accumulated mass on the PUF (ng) and t is the total deployment time

94

(d). This approach assumes that  is constant over the entire sampling period. However, daily

(2)

4

ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24

Environmental Science & Technology

95

observations in urban areas consistently invalidate this assumption.17, 18 Hsu et al. showed that

96

variability in ambient air concentration for ∑PCBs in 22 daily samples taken in a 50-day period

97

near Chicago (0.22 to 1.80 ng m-3) using active sampling was larger than the average

98

concentration (1.0 ± 0.5 ng m-3). 17 The high short-term variability in PCB analyte air

99

concentration can be explained by the influences of meteorology on volatilization (temperature)

100

and source location (wind); in that study, concentration was highly correlated with daily time

101

series of temperature (r2 = 0.46) and slightly correlated with wind speed (0.19) and wind

102

direction (0.17). The assumption of negligible loss from the PUF (i.e.,  ≈ 0) throughout the

103

deployment period, another necessary precondition for this reduced-form expression of the

104

fundamental passive sampling relationships, may also be invalid for SVOCs depending on the

105

length of the deployment.





106

In this study, a method to estimate hourly sampling rates, analyte air concentrations

107

(CAir), and PCB uptake on PAS-PUF samplers (CPUF) from first principle chemistry, physics,

108

mass transfer, and computational fluid dynamics is presented. Using data from seven

109

deployments of samplers from our PCB sampling network in Chicago6, 19 flow rates are

110

calibrated with empirical depuration approaches to refine the model, and process contributions

111

to uncertainties are estimated.

112 113

Experimental Methods & Data

114

Modeling Approach This model reconsiders the theoretical passive sampling relationship in Equation 1 to

115 116

incorporate temporal variability in mass transfer, sampling rate, and PUF concentration based

117

on discrete temporal information on ambient air concentrations and meteorology (Equation 3).

118

 

  , !", !#, $, %& ∗ ∗   , %& 

 & & 

(3)

5

ACS Paragon Plus Environment

Environmental Science & Technology

119

where CAIR is a function of both ambient temperature (T) and time (t), CPUF varies in time

120

due to resolved mass transfer, and kv varies in time as a result of meteorological conditions

121

known from prior empirical and theoretical studies to impact passive sampling rates, including

122

wind speed (WS), wind direction (WD), temperature and pressure (P). The dependency of the

123

mass transfer coefficient on these ambient parameters is derived from boundary layer mass

124

transfer theory, where the Sherwood number representing the ratio of convective to diffusive

125

mass transfer is dependent on the flow conditions represented by the Reynolds Number and the

126

ratio of the diffusive momentum to the diffusion of mass.20 From this theory, the sampling rate

127

equation becomes

128

'  ()* ∗ +, ∗ # 1-. ∗ / .-0 ∗ 1 0 ∗ 20-1 3 ∗

129

where fPBL is the fraction of ambient analyte air concentration inside the PUF boundary layer,

130

and the first bracket represents the mass transfer coefficient, where γ is an advective mass

131

transfer coefficient specific to the PAS housing and media (dimensionless), D is the molecular

132

diffusivity for each compound (m s-1), ν is the kinematic viscosity (m2 s-1), VA is the air velocity

133

(m s-1), L is the diameter of the PUF disk (m), β represents the differences between momentum

134

transport and mass transport (dimensionless), and α represents the flow regime (laminar or

135

turbulent) as a function of air velocity. We considered α and β values consistent with flow over a

136

flat plate. Expressions and values for each component in Equation 4 are presented in the SI.

137

(4)

Uptake of an SVOC on the PAS is a function of accumulation of mass on the PUF over

138

time and the mass of the PUF, as indicated in Equation 5:

139

  4-1

140

where dm/dt is the accumulation of SVOC mass onto the PUF for that time period (ng), and

141

mPUF is mass of the PUF (g). To calculate the CPUF to KPUF ratio, KPUF was computed using a

142

logarithmic relationship between the octanol-air partition coefficient KOA and KPUF from Shoeib &

143

Harner7 (Equation 6):



5⁄ 5

(5)

6

ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24

144

Environmental Science & Technology

  100.6366∗log89

&-3.1774

(6)

145

Using this dynamic model, we can predict the meteorology-driven variation in the

146

ambient concentrations of the SVOC during the deployment period. The hourly concentration of

147

SVOCs in the atmosphere at ambient temperature T can be computed by first correcting the

148

measured integrated air concentration to a standard temperature, then modeling the hourly

149

concentrations as a function of hourly meteorology during deployment period. Using Equation 7:

150

  :;< ∗ =

151

where C298 is the PCB concentration at 298K (ng m-3) calculated from the integrated CPUF and

152

analyzed mass, ∆UOA is octanol-air enthalpy for PCBs (J/mol), R is the gas constant (J mol-1-K-

153

1

154

CAir varies predictably with temperature for SVOCs, the Van’t Hoff relationship explains

155

variability between T and a reference temperature (here, 298 K). The ratio of T to TPAS considers

156

volatility of PCBs in the conditions within PAS chamber, as temperature within the chamber

157

differs from ambient air temperature due to solar heating during the day and radiative cooling at

158

night,13 and is consistent with the temperature of the local built environment.13 Thus, the ratio of

159

T/TPAS ensures that the computed CAir reflects ambient air concentrations rather than the

160

concentration in the PAS to which the PUF disk was directly exposed. Further details are

161

presented in the Supporting Information (SI).

162

>89 1 1 @ - B A A298 ?

∗ C

C

9D

(7)

), T298 is the temperature 298 K, and TPAS is the PAS internal chamber temperature (K). Since

The fraction of ambient analyte air concentration present in the air flowing immediately

163

over the PUF disk, fPBL, determines the SVOC mass capable of transfer to the PUF.

164

Computational fluid dynamics modeling found that simulated PCB concentration in the turbulent

165

boundary layer around the PUF disk in a Harner model sampler is lower than that in free flowing

166

air in the PAS chamber, with fPBL during turbulent flow (internal air velocities ≥ 0.5 m s-1) of

167

approximately 0.5. At lower air velocities, laminar flow prevails, with fPBL ≈ 1.16

168 7

ACS Paragon Plus Environment

Environmental Science & Technology

169 170

PAS PUF Observations & Loss of Depuration Compound Experiments Hourly gas phase RS were calculated in MATLAB21 for 156 PCB congeners and co-

171

eluting groups of congeners22 in Chicago during 2008, including two isotope labeled depuration

172

compounds with negligible ambient concentrations, 2,4,4’-trichlorobiphenyl (13C12) (C13-PCB

173

28) and 2,3,3’,5,5’-pentachlorobiphenyl (13C12) (C13-PCB 111), representing low and high

174

molecular mass PCB compounds. A total of 83 PAS PUFs were deployed during the study

175

period in Harner model samplers at 32 locations in Chicago, and sampling rates determined

176

from loss of depuration compounds from seven PUF samples at seven sites.19 The PUF

177

deployments with depuration compounds covered most of 2008, except for an eight week period

178

from mid-September to mid-November, and 11 days in late December and early January.

179

Deployments ranged from 24 to 76 d, with an average of 43 ± 8 d. All samples were extracted

180

and analyzed following methods described by Hu et al. (2008).23 Using methods described by

181

Persoon et al. (2010),6 the depuration compound experiment consisted of spiking clean PUF

182

disks with 40 ng of PCB congeners 28 (13C-PCB 28) and 111 (13C-PCB 111) from a mixture of

183

EPA Method 1668A labeled clean-up standard solution (Cambridge Isotope Laboratories). The

184

PUFs were spiked by placing them in amber jars and then adding the depuration compounds in

185

20 mL acetone. The PUFs were dried in the jars, capped with Teflon lids, sealed in Ziploc bags

186

and stored in a -20°C freezer until deployment. Each PUF was deployed at one of the seven

187

locations. After deployment, PUF disks were collected and analyzed to determine the remaining

188

PCB mass of each congener. The remaining fraction of the spiked depuration compounds was

189

29% and 96%, for C13-PCB28 and C13-PCB11, respectively. We determined kv from

190

difference between the initial spiked tracer mass and final PUF analyzed mass, and then

191

determined RS from the product of kv x As (Table 1).

192 193

Meteorology 8

ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24

194

Environmental Science & Technology

Temporal variability in mass transfer was calculated using hourly modeled and observed

195

meteorology. The Weather Research and Forecasting Model (WRF) version 3.324, 25 was used

196

to simulate meteorology and land surface temperatures at 1.33 km horizontal resolution over

197

Chicago and 12 km over the Upper Midwest and Northeast United States. WRF operates by

198

integrating the atmospheric primitive equations with physics parameterizations, including sub-

199

models simulating land surface processes, boundary layer dynamics, convection, microphysics

200

and radiation.26 A detailed description of the WRF configurations is presented in the SI in the

201

Meteorology Simulation section. Results from WRF for this application met or exceeded

202

community standards for atmospheric chemical transport modeling, as presented in SI under

203

Meteorological Evaluation. We also employed hourly observations from Automated Weather

204

Observing System (AWOS) units at 16 local area airports, obtained from the Research Data

205

Archive at the National Center for Atmospheric Research, Computational and Information

206

Systems Laboratory.27 From both WRF and AWOS observations, we used temperature, three-

207

dimensional wind vectors, surface pressure and relative humidity as inputs to the sampling rate

208

model. We employed 2 m ambient air temperature from WRF and AWOS, and used WRF

209

surface temperature (TSKIN) as a proxy for TPAS, the sampler’s internal temperature.

210

As meteorology simulated by numerical weather models includes errors in accuracy that

211

may influence sampling rate, hourly de-biased meteorology was calculated by subtracting the

212

hourly domain-average bias (the average difference between model and observation at the 16

213

airports) from the WRF simulated meteorology to yield bias-corrected meteorological fields.

214

These bias-corrected fields thus combine the accuracy of hourly observations with the fine

215

spatial features resolved by urban and regional scale modeled meteorology.

216

The impacts of PAS chamber air velocity on kv were calculated from the empirical linear

217

relationship determined from wind tunnel experiments relating ambient wind speed Va to internal

218

air velocity,11 as the Harner design reduces internal wind speed, advection, and in turn RS. We 9

ACS Paragon Plus Environment

Environmental Science & Technology

219

also employed a polynomial regression on empirical results reported by May et al. from wind

220

tunnel studies, where RS depends on the angle between the 3D wind vector and the PAS

221

apparatus.15 We considered wind vectors for swinging tethered samplers as deployed,

222

accounting for sampler mass in an hourly equilibrium angle to the wind. Results differed by less

223

than 1% using this second approach (SI).

224 225

Results and Discussion

226 227

We first evaluate the model using depuration experiments and quantify hourly variability and

228

process contributions and estimate uncertainty at a representative site. We identify variability

229

between compounds, months, and years at the Chicago sampling network, and extend analysis

230

to simulated spatial patterns across Chicago and the Great Lakes. We then evaluate linearity in

231

hourly PUF-air mass transfer and contextualize RS uncertainty.

232 233

Model Evaluation

234

We evaluated modeled RS for the two C13-labeled depuration compounds, C13-PCB 28 and

235

C13-PCB 111. First, values for the advective mass transfer coefficient γ were calibrated to yield

236

an empirical constant relating the observed sampling rates specific to the Harner samplers’

237

design to first principles mass transfer processes. We fitted γ using Equation 4 at seven

238

sampling sites in Chicago by Euler’s method from C-13 PCB 28 depuration compound RS

239

(Table 1) and the numerical model. The median γ C13-PCB 28 value from the seven

240

deployments was 0.076. First principles only account for flow over a foam disk floating in space,

241

with no enclosure, so this reduction is expected, especially since the Harner sampler is

242

designed to modulate flow rate. We then estimated modeled RS for each PUF. Normalized

243

mean difference between modeled and depuration C13-PCB 28 RS using the median γ value

244

was -5% and normalized mean absolute difference was 22%. These differences are consistently 10

ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24

Environmental Science & Technology

245

less than or comparable to the estimated analytical uncertainty of 15-25% for congener-specific

246

PCB depuration compounds for our laboratory.22 Thus, model results fall within the

247

contemporary range of analytical uncertainty even when using a single representative median

248

value for γ. Results were similar for C13-PCB 111 (SI), and indicated that γ scales between

249

compounds by the ratio of log(KOA), with normalized mean difference of -4% and normalized

250

mean absolute difference of 36% using γ values (median = 0.087) scaled from C13-PCB 28.

251

Sample-specific anomalies in fitted γ (Table 1) were investigated further and found to be

252

correlated with both wind speed (R2 = 0.76) and temperature (R2 = 0.46), suggesting additional

253

process contributions to mass transfer not yet resolved by the model. Multiple linear regression

254

on wind speed and temperature (SI) provide γ estimates tailored to each deployment (MLR γ),

255

which further improve estimates of modeled RS, (Table 1).

256 257 258

Table 1. Modeled and depuration PCB 28 RS for each depuration compound sample. RS and γ reported for sample-specific fit, multiple linear regression, and mean values. A median value for γ of 0.076 is suggested from all samples.

Site

Sampling Period

Sample Specific γ

MLR γ

Depuration RS (m3 d-1)

Dewey Dawes Corkery Gresham Chase St. Elizabeth Webster

1/17 – 2/26 2/22 – 4/11 4/18 – 5/16 5/16 – 7/31 6/30 – 8/11 8/11 – 9/19 11/13 – 12/19

0.06 0.08 0.09 0.08 0.19 0.10 0.06

0.07 0.07 0.06 0.11 0.18 0.09 0.06

5.00 5.53 5.94 3.84 6.79 5.67 5.06

Mean

0.09

0.09

5.40

259 260 261 11

ACS Paragon Plus Environment

Average Model RS (m3 d-1) Using Using Median MLR γ γ (0.076) 6.43 6.06 5.55 5.31 5.13 4.28 3.87 5.63 2.66 6.28 4.55 5.52 6.51 4.81 4.96

5.41

Environmental Science & Technology

262

Hourly Compound-Specific Sampling Rates

263

Using this modeling approach, hourly RS can be determined anywhere and at any time for all

264

SVOCs from meteorology, KOA, and molecular weight alone, using empirical mass transfer

265

coefficient γ from prior depuration experiments. The sampling rate RS is not constant during a

266

deployment period, and is highly dependent on local hourly weather. Modeled hourly RS for PCB

267

congener 28 vary by a factor of 5 (Figure 1) over a PAS deployment period due to meteorology.

268

Modeled RS resolves diurnal and synoptic weather cycles in sampling, seen as the small and

269

large amplitude peaks, respectively. Sampling rate was also computed using averaged

270

simulated meteorology during each PUF deployment to determine if average values could be

271

simply scaled from reference conditions based on representative average deployment

272

meteorology. As shown in Figure 1, average meteorology greatly underestimated RS. Average

273

VA and TSkin were approximately 40% and 21% of the hourly variability seen in VA and TSkin

274

during the deployment period. Since sampling rate is an exponential function of both VA and T,

275

the use of average values produces consistently lower advective and diffusive sampling rates.

276

The PCB 28 RS time series at the Dawes site from 2/22/2008 to 4/11/2008 indicated

277

minimal impact of weather model biases on modeled RS. We found a 1.8% average difference

278

between debiased and raw model meteorological data. In contrast, temperature plays a vital

279

role in variability in RS, where using TSkin rather than ambient T as a proxy for TPAS leads to a

280

7.3% average reduction in total RS. A complete sensitivity analysis is presented in SI.

281

12

ACS Paragon Plus Environment

Page 12 of 24

Page 13 of 24

Environmental Science & Technology

282 283 284 285

Figure 1. RS at the Dawes sampling site from 2/22/2008 to 4/11/2008 for PCB congener 28, including the hourly Rs predicted from our model (black); average loss of depuration Rs (red); modeled RS using average meteorology (green); and RS due to molecular diffusion (blue).

286 287 288

Diffusive and Advective Contributions to Sampling Rate Sampling rate was disaggregated into advective and diffusive components to quantify

289

their roles in the net sampling rate. Figure 1 shows hourly total and diffusive sampling, with

290

advective contribution as the difference. The relative contributions vary with wind speed, with

291

higher PAS chamber air velocities leading to a larger contribution by advective mass transfer.

292

The average RS from diffusion for PCB 28 at the Dawes sampling site during February-April

293

2008 was 2.73 ± 1.07 m3 d-1, consistent with average reported indoor RS values of 2.45 for PCB

294

28 determined by loss of depuration compounds.5 Thus, this model might also be used to

13

ACS Paragon Plus Environment

Environmental Science & Technology

295

estimate indoor RS, due primarily to diffusion, using indoor temperatures and air circulation

296

rates.

297 298 299

Sampling Rate Variability among SVOCs We found simulated RS to be congener dependent, although the effect is smaller than

300

variability in meteorology. Rs increases with KOA. There is a consistent 27% range in RS

301

between low and high molecular weight congeners in each deployment (Figure 2) and a 20%

302

range among homolog groups (Table 2). Sampling rate decreases with increasing molecular

303

weight along KOA isopleths, indicating that differences in RS between SVOCs under ambient

304

conditions are due more to variability in KOA (γ) than to differences in molecular diffusion (D).

305

Therefore, accurate estimation of SVOC concentrations from PAS PUF samples requires

306

compound-specific sampling rates rather than the use or averaging of an arbitrary set of

307

depuration compounds, and highlights this practice as a source of uncertainty in contemporary

308

concentration estimates from passive air sampling. By removing the analytical uncertainty

309

inherent to depuration experiments, these findings explain the variability in RS values derived for

310

a range of compounds in prior depuration studies,28-31 and show that the theoretical variability in

311

RS is of similar magnitude to contemporary analytical uncertainties of ~20-25% for PCBs in

312

ambient air 32, 33.

14

ACS Paragon Plus Environment

Page 14 of 24

Page 15 of 24

Environmental Science & Technology

313

Figure 2. RS dependency on Log(KOA) and homolog group (colors, from pink = 1 to black = 10)

314

for each PCB congener during the PUF deployment at the Dawes sampling site from 2/22/2008

315

to 4/11/2008.

316 317

Not only is there congener-specific variability in RS between using regional versus local

318

scale meteorology, there is also seasonal and annual variability, with maxima in winter and

319

summer (~40% greater than boreal spring). There is less than 5% difference in average monthly

320

RS in the warmest months (July and August) compared to the coldest (January and February),

321

with advection and diffusion each contributing 50 ±6% to the total. Differences in monthly

322

average RS between 2008 and 2011 using observed meteorology at Chicago O’Hare Airport

323

varied from 0.97) over 60+ day deployment periods in winter and spring. In the





18

ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24

Environmental Science & Technology

384

cooler months, hourly uptake follows a staircase pattern, with spikes corresponding to large

385

episodic increases in ambient wind speed and passive sampling rate during synoptic frontal

386

passages, and higher linear uptake rates during sustained periods of high winds. While PUF

387

uptake follows this quasi-linear trend for 12-13 weeks throughout the winter, this is not the case

388

during the warmer months. More volatile and low molecular weight (LMW) compounds begin to

389

saturate after 6 weeks of deployment in the summer, when linear uptake gradually slows,

390

followed by loss in concentration from the PUF into the atmosphere. Sampling rates respond to

391

saturation, and are also sensitive to the dynamic impacts of temperature on saturation levels.

392

After the PUF becomes saturated at high temperatures, LMW compounds volatize off the PUF

393

disk until a new PUF concentration equilibrium point is reached. That equilibrium concentration

394

may reach higher values during cold spells (e.g. near 8/15 in Figure 4), and then lower again as

395

the temperature rises. This trend continues if the PUFs are left out longer than 12 weeks.

396

Saturation can lead to a >20% underprediction of quantified mass for LMW compounds.

397

However, this is only the case for volatile LMW compounds, as linear uptake occurs for less

398

volatile congeners throughout long summer deployments (Figure 4).

399

Modeled mass uptake onto the PUF shows that it is likely that previous methods used

400

for determining CAir have underpredicted concentrations for LMW compounds during long

401

deployments in warm conditions if one assumes constant linear uptake throughout as PAS PUF

402

deployment. In turn, the quantified analyzed mass may not be representative of the true mass

403

uptake during PAS PUF deployment period. Saturation can be estimated from KPUF and can be

404

corrected for using temperature, deployment length and kV from depuration compounds, 35, 36

405

but only if a representative range of such compounds is included on each PUF. In the Great

406

Lakes and other areas with high variability in temperature, cloud cover and wind speed during a

407

PUF deployment, even correcting using daily average temperature may not resolve the effects

408

of PAS internal temperatures within a metal housing under full sunlight. However, nonlinear 19

ACS Paragon Plus Environment

Environmental Science & Technology

409

uptake is not a requirement when using the numerical model described here, as loss can be

410

quantified and CAir calculated from the linear uptake period.

411 412

413 414

Figure 4. Simulated PUF concentration (ng/g) at the Dawes sampling site from June 1 to

415

August 31 2008: PCB 11 (blue, left axis) and PCB 153 (green, right axis). For low molecular

416

weight congeners, concentration on the PUF during the first 6 weeks (red) increases linearly

417

(black trendline, R2 = 0.9).

418 419

Remaining Uncertainties in Passive Sampling

420

This study identifies and quantifies uncertainties in the knowledge of passive air

421

sampling processes, and provides a modular platform for evaluating their impacts on sampling

422

rate in the field. We identify temperature in the PAS chamber as the single most important

423

unresolved condition for predicting seasonality in RS. The simulated net process uncertainty of

424