Within-City Spatial Variations in Multiple Measures of PM2.5 Oxidative

2 days ago - Few studies have characterized within-city spatial variations in the oxidative potential of fine particulate air pollution (PM2.5). In th...
1 downloads 0 Views 2MB Size
Subscriber access provided by University of Massachusetts Amherst Libraries

Ecotoxicology and Human Environmental Health

Within-City Spatial Variations in Multiple Measures of PM2.5 Oxidative Potential in Toronto, Canada Scott Weichenthal, Maryam Shekarrizfard, Alison Traub, Ryan Kulka, Kenan Al-Rijleh, Sabreena Anowar, Greg J. Evans, and Marianne Hatzopoulou Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b05543 • Publication Date (Web): 08 Feb 2019 Downloaded from http://pubs.acs.org on February 10, 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 50

Environmental Science & Technology

1 2 3 4 5

Within-City Spatial Variations in Multiple Measures of PM2.5 Oxidative

6

Potential in Toronto, Canada

7 8

Scott Weichenthal*1,2, Maryam Shekarrizfard3, Alison Traub4, Ryan Kulka2, Kenan Al-

9

Rijleh3, Sabreena Anowar3, Greg Evans4, Marianne Hatzopoulou3

10 11 12 13 14

1McGill

University, Department of Epidemiology, Biostatistics and Occupational Health,

15

Montreal, Quebec, H3A 1A2, Canada

16

2 Air

17

3University

18

Canada

19

4University

20

Toronto, Ontario, M5S 3E5, Canada

Health Science Division, Health Canada, Ottawa, Ontario, K1A 0K9, Canada of Toronto, Department of Civil Engineering, Toronto, Ontario, M5S 1A4, of Toronto, Department of Chemical Engineering and Applied Chemistry,

21 22 23

*Corresponding Author

24

Scott Weichenthal

25

Faculty of Medicine

26

Department of Epidemiology, Biostatistics, and Occupational Health

27

McGill University

28

1020 Pins Ave. West

29

Montreal, QC H3A 1A2, Canada

ACS Paragon Plus Environment

1

Environmental Science & Technology

30

Email: [email protected]

31

Tel: (514) 398-1584

Page 2 of 50

32 33 34 35 36 37 38 39 40 41 42 43 44

Abstract

45

Few studies have characterized within-city spatial variations in the oxidative

46

potential of fine particulate air pollution (PM2.5). In this study, we evaluated multiple

47

measures of PM2.5 oxidative potential across Toronto, Canada (2016-2017) including

48

glutathione/ascorbate-related oxidative potential (OPGSH and OPAA) and dithiothreitol

49

depletion (OPDTT). Integrated 2-week samples were collected from 67 sites in summer

50

and 42 sites in winter. Multivariable linear models were developed to predict OP based

51

on various land use/traffic factors and PM2.5 metals and black carbon were also

52

examined. All three measures of PM2.5 oxidative potential varied substantially across

ACS Paragon Plus Environment

2

Page 3 of 50

Environmental Science & Technology

53

Toronto. OPAA and OPDTT were primarily associated with traffic-related components of

54

PM2.5 (i.e. Fe, Cu, and black carbon) whereas OPGSH was not a strong marker for traffic

55

during either season. During summer, multivariable models performed best for OPAA

56

(R2CV=0.48) followed by OPDTT (R2CV=0.32) and OPGSH (R2CV=0.22). During winter,

57

model performance was best for OPDTT (R2CV=0.55) followed by OPGSH (R2CV=0.50) and

58

OPAA (R2CV=0.23). Model parameters varied between seasons and between-season

59

differences in PM2.5 mass concentrations were weakly/moderately correlated with

60

seasonal differences in OP. Our findings highlight substantial within-city variations in

61

PM2.5 oxidative potential. More detailed information is needed on local sources of air

62

pollution to improve model performance.

63 64 65 66 67 68

1.

Introduction

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 50

69 70

Outdoor fine particulate air pollution (< 2.5 um, PM2.5) is an important contributor

71

to the overall global burden of disease and is responsible for millions of deaths around

72

the world each year.1 Traditionally, PM2.5 exposures have been measured as particle

73

mass concentrations (i.e. ug/m3); however, this metric is limited in that it treats all

74

particles as equally toxic and does not consider spatial or temporal differences in

75

particle composition or “biological activity”. As oxidative stress and the production of

76

reactive oxygen species are thought to play an important role in explaining air pollution

77

health effects,2-5 more recent studies have also examined particle oxidative potential

78

(OP) as a complementary measure to traditional mass-based exposure measurements.

79

Numerous assays are available to estimate particle oxidative potential which

80

provides an integrated measure of the capacity of PM to oxidize target molecules.3,6 In

81

general, acellular assays are most common and include those based on dithiothreitol

82

consumption (OPDTT), electron spin resonance (OPESR), or the depletion of the

83

antioxidants ascorbate (OPAA) and glutathione (OPGSH) in a synthetic respiratory tract

84

lining fluid.7-14 A recent study examined temporal variations in particle OP over a 1-year

ACS Paragon Plus Environment

4

Page 5 of 50

Environmental Science & Technology

85

period in France and noted higher OP levels during the winter months for the OPDTT,

86

OPAA, and OPGSH assays with little seasonal difference reported for OPESR.15 Moreover,

87

OPESR was predominantly associated with transition metal components whereas the

88

other assays were correlated with both organic and inorganic species reflecting sources

89

including biomass burning, brake/tire wear, and traffic/fossil fuel combustion.15 Likewise,

90

other studies of particle oxidative potential and chemical constituents suggest that

91

organic and inorganic particle components both contribute to the oxidative potential of

92

airborne particles. 9-10,12,16-17

93

In addition to temporal variations, spatial variations in particle oxidative potential

94

have also been reported.7, 12,14,17-20 In Canada, between-city differences in OPGSH have

95

been shown to modify associations between long-term PM2.5 exposures and lung cancer

96

mortality21 as well as short-term associations between PM2.5 and emergency room visits

97

for myocardial infarction22 and respiratory outcomes.14 To date, land use regression

98

models for within-city spatial contrasts in particle oxidative potential have been limited to

99

European cities. In particular, these studies suggest that traffic proximity/road length,

100

harbor proximity, population density, green space, emissions from brake and tire wear,

ACS Paragon Plus Environment

5

Environmental Science & Technology

101

and industrial land use may be important predictors of within-city spatial contrasts in

102

particle oxidative potential.12,17-20

103

In this study, we examined within-city spatial variations in three measures of

104

PM2.5 oxidative potential in Toronto, Canada including OPAA, OPGSH, and OPDTT. We

105

also collected detailed information on PM2.5 metals and black carbon content. To our

106

knowledge, this is the first study to attempt to model within-city spatial variations in

107

PM2.5 oxidative potential in North America.

Page 6 of 50

108 109 110 111 112 113 114 115 116 117 118 119 120

2.0

Methods

2.1

Spatial Monitoring Study

121 122

ACS Paragon Plus Environment

6

Page 7 of 50

123

Environmental Science & Technology

The details of our spatial monitoring campaign have been described previously.23

124

Briefly, two PM2.5 monitoring campaigns were conducted in Toronto, Canada during

125

summer 2016 and winter 2017. In total, 67 were monitored in summer and 42 sites were

126

monitored in winter; 29 sites were monitored during both seasons. Daily mean

127

temperatures in summer ranged from 18.3oC to 27.2oC whereas winter temperatures

128

ranged from -8.8oC to 7.5oC. Total precipitation during the summer (38.2 mm) and

129

winter (39.7 mm) sampling periods were similar.

130

Integrated 2-week samples were collected at each site using Teflon filters with

131

cascade impactors at a flow-rate of 5 liters/minute. During each season, all PM2.5

132

samples were collected simultaneously using pre-set pump timers and separate

133

sampling kits at each monitoring site. This eliminated the need to adjust for temporal

134

differences in PM2.5 mass/components between sites as all sites were monitored over

135

exactly the same time-period. PM2.5 mass concentrations were determined

136

gravimetrically prior to metals analysis by X-ray fluorescence (XRF) according to EPA

137

Method IO-3.3 in Compendium of Methods for the Determination of Metals in Ambient

138

Particulate Matter (EPA 625/R-96/010a). All PM2.5 samples were also analyzed using a

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 50

139

SootScan Model OT21 Transmissometer at 880 nm (BC880nm) and 370 nm (BC370nm).

140

Measurements at 880 nm are interpreted as black carbon whereas measurements at

141

370 nm respond more strongly to biomass PM from sources such as

142

residential/commercial wood burning. These two measures were highly correlated in our

143

data (r > 0.93) and thus we simply refer to black carbon.

144 145 146 147

2.2

PM2.5 Oxidative Potential Analyses PM2.5 was extracted from filters into HPLC grade methanol (Caledon) through

148

vortexing (1,800 rpm) followed by sonication, for 10 minutes each, and with rinse with

149

1.0 mL of methanol. The extraction efficiency, which was typically above 75%, was

150

determined from the pre and post extraction filter weight. The methanol was evaporated

151

at 37oC using a 34-position nitrogen evaporator (N-EVAP, Organomation). PM2.5 was

152

then resuspended in ultrapure water (Caledon) containing 5% HPLC methanol to a

153

concentration of 75 μg/mL. The re-suspended PM2.5 samples were analyzed in triplicate

154

using three acellular in vitro assays: the ascorbate (AA) assay, the glutathione (GSH)

ACS Paragon Plus Environment

8

Page 9 of 50

Environmental Science & Technology

155

assay, and the dithiothreitol (DTT) assay. Ascorbate and glutathione oxidative potential

156

(OPAA and OPGSH) were assessed using a methodology similar to that described

157

previously.14 Briefly, PM2.5 samples were incubated at a concentration of 75 μg/mL for 4

158

hours at 37oC with synthetic respiratory tract lining fluid (RTLF) containing equimolar

159

concentrations of AA, GSH, and uric acid in an ultraviolet-visible plate reader (Molecular

160

Devices, Spectra Max 190) alongside positive controls (1.0 μM Cu(NO3)2, 0.02% H2O2)

161

and blanks. Ascorbate depletion was calculated as the percent difference between the

162

initial and final absorbance values measured at 260 nm and blank-corrected. After the

163

4-hour incubation period, GSH depletion was measured using the glutathione-reductase

164

enzyme recycling assay.24 Percent depletion was calculated by comparing GSH

165

concentrations measured in samples and positive controls to those measured in blanks.

166

Depletion was typically under 40% so a linear instead of first order measure of depletion

167

was adequate.

168

Dithiothreitol oxidative potential (OPDTT) was assessed using a methodology

169

similar to that described in Cho et al.25 Briefly, re-suspended PM2.5 samples were

170

incubated with 100 μM dithiothreitol (DTT) in a 96 well plate alongside positive controls

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 50

171

(5.0 μM Cu(NO3)2), blanks, and DTT standards (containing 0 – 100μM DTT) for 35

172

minutes at 37oC, with constant shaking. After 5, 15, 25, and 35 minutes, remaining DTT

173

was measured by adding 1.0 mM 5,5’-dithiobis(2-nitrobenzoic acid) (DTNB) to each well

174

and measuring absorbance at 412 nm. Samples were initially analyzed at a

175

concentration of 75 μg/mL; however, if DTT depletion exceeded 25% after 35 minutes,

176

the sample was re-analyzed at a lower concentration (37.5 μM or 50 μM) as depletion

177

rates can become non-linear beyond 25% depletion.

178

OPAA, OPGSH, and OPDTT were first normalized by PM2.5 mass based on the

179

quantity used in each assay (i.e. % depletion/µg for OPAA and OPGSH and pmol/min/µg

180

for OPDTT); these values were then converted to units of % depletion/m3 (for OPAA and

181

OPGSH) or pmol/min/m3 (for OPDTT) based on PM2.5 mass concentrations at each

182

location in units of µg/m3. To facilitate comparisons with other studies, descriptive data

183

for OPAA and OPGSH are also presented in units of pmol/min/m3 by multiplying

184

%Depletion/m3 values by 9.3 (i.e. 1% depletion/m3 = 9.3 pmol/min/m3 based on the 4-

185

hour incubation period, 200 µM of AA (or GSH) initially present, and the 0.18 ml of PM

186

solution added to the 0.2 ml well). Negative values were observed for OPGSH in 13

ACS Paragon Plus Environment

10

Page 11 of 50

Environmental Science & Technology

187

summer samples and 5 winter samples: these values were replaced with half the

188

analytical detection limit (0.077%/g or 0.72 pmol/min/ug).

189

2.3

190

Parameters for land use regression models Various land-use and built environment characteristics were derived for each

191

monitoring site using ArcMap 10.4.1 (ESRI, Redlands, CA). Land use composition was

192

computed within circular buffers of 100m, 200m, 300m, 500m, 700m and 1000m, and

193

included specific categories for residential, commercial, governmental/institutional,

194

resource/industrial, parks, open area, water, as well as building footprints (DMTI Spatial

195

Inc. Database 2014). Moreover, the length of bus routes, highways, major roads, and

196

rail lines were derived in each buffer (City of Toronto Open Data portal 2016). Other

197

variables computed included population density (Toronto Neighbourhood and

198

Demographics data of 2013), mean and maximum building height, number of bus stops,

199

number of road intersections, and number of trees (City of Toronto Open Data portal

200

2016). Average hourly traffic volume between 6am and 7pm was estimated using a

201

traffic assignment model and an origin-destination (OD) matrix of trips extracted from

202

the Transportation Tomorrow Survey (TTS) (year 2011) for the Greater Toronto and

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 50

203

Hamilton Area (GTHA). Distances between each sampling location and the closest rail

204

line, major road, highway, Pearson airport, lakeshore, and city center were also

205

estimated. Industrial emissions of CO, NOx, total VOCs, total metals, and PM were

206

evaluated as surrogate measures of local industrial combustion activities. Emissions

207

values were assigned to each location using the nearest facility with emissions

208

(tonnes/year) weighted by 1/distance to the closest facility. Emissions data were

209

obtained from the National Pollutant Release Inventory. Finally, the locations of

210

restaurants and gas stations in Toronto were extracted using the Google Places

211

function in R and counts of these were assigned to each monitoring site with the circular

212

buffers described above.

213

2.4

Statistical Analyses

214

Models were developed to predict within-city spatial variations in OPAA (%

215

depletion/m3), OPGSH (% depletion/m3), and OPDTT (pmol/min/m3); separate models were

216

developed for the summer and winter seasons. Annual models are not presented as

217

only 29 sites were monitored during both the summer and winter campaigns.

ACS Paragon Plus Environment

12

Page 13 of 50

Environmental Science & Technology

218

Alternatively, data collected at sites during both the summer and winter months were

219

used to evaluate between-season differences in OP across Toronto.

220

Our model building procedure followed several steps: 1) First, we examined

221

single variable linear regression models and associated scatter plots with linear and

222

non-linear lines of fit to identify potential non-linear relationships (quadratic terms were

223

evaluated for non-linear relationships); 2) Variables that were associated with the

224

outcome (i.e. 95% CI excluded the null) were retained for potential inclusion in the final

225

model; 3) Spearman’s correlations were determined for the candidate predictors

226

identified in step 2 and highly correlated variables (r > 0.7) were removed (retaining the

227

best predictor of each correlated pair); 4) All parameters retained after step 3 were

228

eligible for inclusion in the final model; 5) “Non-significant” parameters were only

229

removed from multivariable models if doing so improved model fit (i.e. decreased the

230

root mean square error); 6) Finally, interaction terms were evaluated between model

231

covariates and were retained if doing so improved model fit. As the primary purpose of

232

statistical modeling was prediction, we decided to include interaction terms to account

233

for the fact that a given land use parameter may be more or less strongly associated

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 50

234

with PM2.5 oxidative potential given the presence of a second land use characteristic.

235

For example, road length and industrial land use may each be related to PM2.5 oxidative

236

potential but road length may be more strongly associated with OP in areas with higher

237

industrial land use because of the types of vehicles that tend to travel in industrial areas

238

(this would be indicated by a positive coefficient for the interaction term between these

239

variables). A leave-one-out cross validation procedure was used to evaluate final model

240 241

performance using the loocv procedure in Stata version 15 (Statacorp, College Station,

242

Texas, USA). As the primary purpose of modelling was for prediction, we did not impose

243

pre-conceived rules on the direction of effects for parameters included in the model. All

244

OPAA, OPGSH, and OPDTT values were log transformed to normalize distributions. All

245

covariates were standardized (by subtracting the mean and dividing by 2 standard

246

deviations) prior to model development.

247

2.5

Mapping the Land Use Regression Models

248

Maps of exposure surfaces were generated by first dividing the city of Toronto

249

into grid cells of 100m by 100m using ArcMap (ArcMap 10.4.1; ESRI, Redlands, CA).

ACS Paragon Plus Environment

14

Page 15 of 50

Environmental Science & Technology

250

Next, the final set of predictors for land use regression models were computed for the

251

mid-point of each grid cell. Finally, the predicted values for each mid-point were

252

calculated and associated with the corresponding grid cell for mapping.

253

3.

254

Results In total, PM2.5 samples were collected from 67 sites in summer and 42 sites in

255

winter; 29 sites were monitored during both seasons. Descriptive data for PM2.5

256

oxidative potential and black carbon are shown in Table 1 and correlations between

257

oxidative potential metrics, PM2.5 metals, and black carbon are shown in Figure 1. In

258

general, correlations between OPAA, OPGSH, and OPDTT were low to moderate during

259

summer (0.27 < r < 0.77) and winter (0.34 < r < 0.58). During summer, OPAA and OPDTT

260

were moderately correlated with traffic-related pollutants including black carbon (OPAA

261

(r=0.63); OPDTT (r=0.51)) and PM2.5 metals including Fe (OPAA (r=0.75); OPDTT (r=0.62)),

262

Cu (OPAA (r=0.66); OPDTT (r=0.63)), and Ba (OPAA (r=0.64); OPDTT (r=0.65)) which are

263

markers of non-tailpipe emissions such as brake/tire wear. OPGSH was weakly

264

correlated with these traffic-related components (0.12 < r < 0.24) during summer and

265

was moderately correlated with Mn and Zn; however, correlations for Mn and Zn were

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 50

266

influenced by a single outlying value and these correlations decreased (r < 0.2) when

267

this point was removed. Correlations between OPDTT and traffic related components (i.e.

268

Fe, Cu, Ba, and BC) decreased during winter (0.18 < r < 0.26) but OPAA remained

269

moderately correlated with these components (0.39 < r < 0.51). OPGSH was weakly

270

correlated with traffic-related components during winter (0.21 < r < 0.30).

271

On average, between-season differences in PM2.5 oxidative potential were small

272

in magnitude for sites monitored during both the summer and winter campaigns (Figure

273

2). Specifically, OPDTT (mean difference: 10.5 pmol/min/m3, 95% CI: 3.70, 17.3) and

274

OPAA (mean difference: 0.227 % Depletion/m3, 95% CI: -0.0728, 0.526) tended to be

275

greater during summer whereas little difference was observed for OPGSH (mean

276

difference: 0.0737 % Depletion/m3, 95% CI: -0.765, 0.912). However, overall averages

277

masked larger between-season differences for individual sites as illustrated in Figure 3.

278

In particular, the magnitudes of between-season differences at individual sites were

279

greater for OPGSH than for OPAA (panels A and B) and only 4 sites monitored during both

280

seasons had higher OPDTT during winter (panel C). Moreover, patterns of seasonal

281

differences for OP were often different than the pattern for PM2.5 mass concentrations;

ACS Paragon Plus Environment

16

Page 17 of 50

Environmental Science & Technology

282

this was particularly true for OPGSH (Figure 3). On average, PM2.5 concentrations were

283

0.992 g/m3 (95% CI: 0.627, 1.36) higher during summer, with only two sites having

284

higher mean concentrations winter (Figure 3). In simple linear models, between-season

285

differences in PM2.5 mass concentrations were not correlated with seasonal changes in

286

OPGSH (R2=0.03) or OPDTT (R2=0.07) but were moderately correlated with between-

287

season differences in OPAA (R2=0.37). Of the specific elements examined, between-

288

season differences in Se (a marker for fossil fuel combustion)26 explained the largest

289

proportion of between-season variability in OPAA (R2=0.26) whereas between-season

290

differences in traffic-related components including Fe, Cu, Ba, and black carbon

291

explained less of this variation (0.10 < R2 < 0.17). Seasonal differences in vanadium

292

(R2=0.13) and black carbon (R2=0.14) explained a small proportion of between-season

293

variability in OPGSH (R2=0.13) whereas, Cu, and Ba explained a similar proportion of

294

between-season variations in OPDTT (0.10 < R2 < 0.12).

295

To explore between-season differences in OP further (i.e. summer – winter), we

296

examined multivariable linear regression models for relationships between the various

297

land use/traffic parameters described above and seasonal differences in OP at each

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 50

298

site (Supplemental Table S1). For OPGSH, bus lines and trees had positive coefficients

299

in these models whereas coefficients for higher traffic counts were negative. For OPAA,

300

park space had a positive coefficient whereas coefficients for rail and street cars had

301

negative coefficients in models for seasonal differences. Finally, models for seasonal

302

differences in OPDTT had positive coefficients for traffic counts and open spaces and a

303

negative coefficient for rail.

304

OPGSH, OPAA, and OPDTT each varied substantially between monitoring sites

305

during both the summer and winter campaigns as illustrated in Figure 4. During

306

summer, multivariable linear models performed best for OPAA followed by OPDTT and

307

OPGSH (Figure 5). Models for OPAA and OPDTT during summer included parameters for

308

traffic, commercial/industrial/residential land use, and industrial emissions. Specifically,

309

length of highways (700m) (R2=0.10), major roads (300m) (R2=0.063), CO emissions

310

(R2=0.084), and government land use (1000m) (R2=0.067) were the strongest predictors

311

of summer OPAA whereas traffic (200m) (R2=0.18), CO emissions (R2=0.12), park space

312

(500m) (R2=0.090), residential land use (1000m) (R2=0.082), and industrial land use

313

(1000m) (R2=0.05) were the strongest predictors for OPDTT. Few important predictors

ACS Paragon Plus Environment

18

Page 19 of 50

Environmental Science & Technology

314

were identified for OPGSH during summer (i.e. CO emissions (R2=0.17), distance to

315

highways (R2=0.05), and restaurants within 1000m (R2=0.09) and this model performed

316

poorly.

317

During winter the multivariable model for OPGSH explained the majority of the

318

spatial variability for this parameter across Toronto and included terms for length of

319

major roads, park space, open space, and interactions between these variables (Figure

320

6). Of these, the inverse association with park space (500m) (R2=0.18) explained the

321

largest proportion of spatial variability in OPGSH during winter followed by major roads

322

(100m) (R2=0.11) and open space (200m) (R2=0.090). The winter model for OPDTT also

323

explained the majority of spatial variability for this parameter and included terms for

324

distance to shore, open space, road junctions, industrial emissions for metals and CO,

325

and interaction terms between these factors. The strongest predictor of spatial

326

variations in OPDTT during winter was open space (300m) (R2=0.18) followed by

327

industrial emissions for metals (R2=0.14), CO emissions (R2=0.11), distance to shore

328

(R2=0.11), and road junctions (100m) (R2=0.041). The winter OPAA model included

329

terms for traffic counts, distance to the airport, commercial land use, distance to rail,

ACS Paragon Plus Environment

19

Environmental Science & Technology

Page 20 of 50

330

buildings, and interaction terms between distance to rail, buildings, and commercial land

331

use but explained little spatial variability in this parameter. The strongest predictors of

332

spatial variations in OPAA during winter was distance to the airport (R2=0.17), followed

333

by commercial land use (1000m) (R2=0.14), traffic (500m) (R2=0.12), length of rail

334

(500m) (R2=0.12), and buildings (500m) (R2=0.040). Predicted spatial variations in

335

summer OPAA and OPDTT across Toronto are shown in Figures 7. These surfaces

336

highlight potential hot-spots for PM2.5 oxidative potential across Toronto, particularly the

337

very high traffic region around the airport.

338

Finally, one additional set of analyses were conducted for OPAA (summer and

339

winter) and OPDTT (summer) to examine predictors of within-city spatial variations not

340

related to transition metals from brake/tire wear. Specifically, we extracted residuals

341

from linear regression models between OPAA/OPDTT and Fe (a marker for brake/tire

342

wear and highly correlated (r > 0.86) with both Cu and Ba) and evaluated predictors of

343

spatial variations in these residuals. The number of trees within 100m (= 0.140, 95%

344

CI: 0.0603, 0.220; R2=0.16) and distance to highways (= -0.0841, 95% CI: -0.162, -

345

0.00577, 0.220; R2=0.068) were the only variables associated with spatial variations in

ACS Paragon Plus Environment

20

Page 21 of 50

Environmental Science & Technology

346

residuals for OPAA during summer whereas park space within 500m was associated

347

with spatial variations in residuals for OPAA during winter (= -0.276, 95% CI: -0.433, -

348

0.120; R2=0.25). Park space within 500m was the only predictor of spatial variation in

349

residuals for OPDTT during summer (= 0.118, 95% CI: 0.00786, 0.229; R2=0.068).

350

4.

351

Discussion We conducted large-scale monitoring campaigns for multiple measures of PM2.5

352

oxidative potential during the summer and winter months in Toronto, Canada and noted

353

several interesting findings. First, OPAA and OPDTT were each correlated with markers of

354

traffic-related air pollution during summer whereas weaker correlations were observed

355

during winter, particularly for OPDTT. OPGSH was not a strong marker of traffic-related air

356

pollution in either season. Interestingly, the summer and winter models for all three OP

357

measures generally contained different predictors. Specifically, no parameters

358

overlapped between seasons for OPGSH, one parameter overlapped for OPDTT (CO

359

emissions), and three parameters overlapped for OPAA (commercial land use, traffic

360

counts, distance to rail). As land use/traffic data are only surrogate measures of the

361

underlying sources/components of interest, it appears that the extent to which these

ACS Paragon Plus Environment

21

Environmental Science & Technology

362

parameters reflect sources/components related to PM2.5 oxidative potential may differ

363

between seasons.

364

Page 22 of 50

In general, the predictive performance of multivariable models for OPAA, OPDTT,

365

and OPGSH was modest and model performance varied between seasons. Indeed, it

366

seems that traditional GIS parameters may not perform as well in predicting spatial

367

variations in PM2.5 oxidative potential as they do for other air pollutants such as NO2.

368

Interestingly, when we examined predictors of spatial variations in OPAA and OPDTT not

369

attributable to transition metals associated with brake/tire wear (i.e. the residual models

370

explained above) few important predictors were identified. This suggests that many of

371

the land use/traffic/built environment parameters included in final models for OPAA and

372

OPDTT were capturing spatial variations in transition metals associated with brake/tire

373

wear. For OPAA, the inverse association with distance to highways in the residual model

374

suggests that other factors (e.g. organic compounds) emitted near highways also

375

contribute to particle oxidative potential independent of transition metals. Moreover, the

376

positive association between OPAA and the number of trees within 100m in the summer

377

residual model is also interesting. One explanation is that trees are a marker for some

ACS Paragon Plus Environment

22

Page 23 of 50

Environmental Science & Technology

378

unmeasured source of PM2.5 with elevated oxidative potential; however, the fact that we

379

did not observe this relationship in the winter OPAA residual model argues against this

380

hypothesis. A second possibility is that biogenic volatile organic compounds released by

381

trees during the summer months undergo photochemical reactions that produce

382

localized increases in secondary organic aerosols that contribute to particle oxidative

383

potential.27 While we cannot reliably confirm or refute either of these hypotheses in the

384

present study, future studies should explore this possibility further. More generally,

385

additional effort is needed to identify and geocode local sources of air pollution

386

(primarily non-traffic sources) that may contribute to PM2.5 oxidative potential,

387

particularly for organic compounds as these also play an important role in predicting OP

388

levels.15 Specifically, detailed inventories of air pollution sources related to biomass

389

burning are not currently available and these sources are particularly relevant as they

390

tend to be local in nature and vary seasonally (i.e. residential heating). Moreover,

391

existing evidence suggests that biomass particles contribute to particle oxidative

392

potential and that OPAA depletion by wood smoke particles is not inhibited by co-

393

incubation with a metal chelator.31,32 Therefore, biomass burning may play an important

ACS Paragon Plus Environment

23

Environmental Science & Technology

394

role in determining within-city spatial variations in PM2.5 oxidative potential and should

395

be explored more thoroughly in future investigations.

Page 24 of 50

396

To date, only a small number of studies have examined within-city spatial

397

variations in OPGSH, OPAA, and OPDTT and none of these studies have been conducted

398

in North America. Yanosky et al.18 examined spatial variations in OPGSH using data from

399

34 PM10 monitoring sites across London, UK and reported a positive association with

400

brake/tire wear emissions and an inverse association with NOx emissions from heavy-

401

goods vehicles. Gulliver et al.20 developed land use regression models for OPAA and

402

OPGSH in five European cities and also measured PM2.5 metals and black carbon (PM2.5

403

absorbance). As in this study, they reported stronger correlations between OPAA and

404

traffic-related components (i.e. Fe, Cu, and PM2.5 absorbance) than for OPGSH. Cross-

405

validation R2 values ranged from 0.31-0.82 for OPAA and 0.05-0.38 for OPGSH with all

406

models containing a small number of predictors related to traffic, green/natural space,

407

and industrial land use. Other studies of within-city spatial variations in OPAA and OPGSH

408

were not identified but three additional studies have examined spatial variations in

409

OPDTT. Specifically, Yang et al.12 reported moderate correlations between OPDTT and

ACS Paragon Plus Environment

24

Page 25 of 50

Environmental Science & Technology

410

traffic-related components of PM2.5 (i.e. Fe, Cu, and PM2.5 absorbance) across 40 sites

411

in the Netherlands/Belgium and a model including terms for traffic, natural space, and

412

regional OPDTT explained the majority of spatial variations across the study area

413

(R2=0.60). Jedynska et al.19 also reported land use regression models for OPDTT across

414

ten European cities but model performance was poor (median R2=0.33) and the authors

415

reported a high number of samples below the limit of detection for OPDTT. Recently

416

,Bates et al.28 described large-scale variations in OPDTT in the United States and found

417

that emissions from vehicles and biomass burning were important predictors of OPDTT.

418

Our findings related to traffic components of PM2.5 and OPDTT in Toronto are consistent

419

with these results. In an earlier study, these authors also reported higher OPDTT in

420

Atlanta during winter32 whereas we generally observed higher values during summer.

421

Reasons for this difference are not entirely clear but could be explained in part by lower

422

rates of particle resuspension in Toronto during the winter months owing to snow.

423

Alternatively, seasonal results between the two studies may not be directly comparable

424

given the much colder winter temperatures in Toronto and different sampling media

425

used between the two studies. Nevertheless, we did note a small number of sites in

ACS Paragon Plus Environment

25

Environmental Science & Technology

426

Toronto with higher OPDTT during winter and thus local sources may play an important

427

role in determining seasonal patterns. As only five sites were monitored in the Atlanta

428

study,32 it is not clear if this seasonal pattern holds across Atlanta area or is limited to

429

the specific monitoring locations.

430

Page 26 of 50

Collectively, existing studies of spatial variations in OPAA, OPGSH, and OPDTT

431

highlight the difficulty of modeling these parameters and provide little consistent

432

evidence related to specific sources of interest beyond fossil fuel/biomass combustion

433

and transition metals associated with brake and tire wear. As noted here and

434

elsewhere20, future studies need to consider sources not currently captured in traditional

435

GIS databases if we are to improve these models for use in population-based studies.

436

This may also help to explain seasonal differences in OP on a local scale as our

437

observation of different model parameters between the summer and winter models

438

suggests that current land use parameters are imperfect measures of local sources that

439

ultimately predict OP. While the models presented above may be useful in future

440

epidemiological studies, a more conservative approach would be to simply use the

441

measurement data collected at each location and apply a buffer around each site to

ACS Paragon Plus Environment

26

Page 27 of 50

Environmental Science & Technology

442

define the study population (e.g. 1km). Indeed, as none of the OP parameters examined

443

had strong models for both seasons, their application in population-based studies may

444

underestimate the true association between particle oxidative potential and chronic

445

health outcomes owing to imprecision in model estimates.

446

This study had a number of important advantages including simultaneous

447

monitoring at all monitoring sites, detailed information on multiple measures of PM2.5

448

oxidative potential, and information on particle components including PM2.5-metals and

449

black carbon. However, it is important to note several limitations. First, only 29 sites had

450

information for both the summer and winter campaigns and thus we could not develop

451

models for spatial variations in annual average PM2.5 oxidative potential across Toronto.

452

It is possible that annual average models may perform better than season-specific

453

models and we will explore this possibility in future studies. As noted above, we

454

observed substantial heterogeneity in model performance between seasons and

455

different predictors were often selected for the summer and winter models for all three

456

measures of oxidative potential. This suggests that the current land use/traffic

457

parameters provide imperfect estimates of spatial variations in air pollution sources

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 50

458

related to OP and further work is needed to document and geocode local sources of air

459

pollution that may contribute to within-city spatial variations in particle oxidative

460

potential.

461

In summary, we conducted a large-scale spatial monitoring study in Toronto,

462

Canada to evaluate within-city variations in multiple measures of PM2.5 oxidative

463

potential. Our findings suggest that PM2.5 oxidative potential varies substantially within

464

cities and that systematic seasonal differences may exist for PM2.5 oxidative potential on

465

a local scale. Moreover, our results highlight the fact that seasonal trends in PM2.5 mass

466

concentrations may not correspond to seasonal changes in particle oxidative potential.

467

Additional work is needed to identify the particle components/sources that explain this

468

discrepancy as well as the geographic parameters than can be used to predict

469

exposures for large population-based studies.

470

Supporting Information

471

Multivariable linear models for between-season differences in OPGSH (%Depletion/m3),

472

OPAA (%Depletion/m3), and OPDTT (pmol/min/m3) in Toronto, Canada (Table S1);

473

Multivariable models for OPAA (%Depletion/m3) in Toronto, Canada (Table S2);

ACS Paragon Plus Environment

28

Page 29 of 50

Environmental Science & Technology

474

Multivariable linear models for spatial variations in OPGSH (%Depletion/m3) in Toronto,

475

Canada (Table S3); Multivariable linear models for spatial variations in OPDTT

476

(pmol/min/m3) in Toronto, Canada (Table S4).

477 478 479 480 481 482 483 484 485 486 487 488 489

ACS Paragon Plus Environment

29

Environmental Science & Technology

490 491

Page 30 of 50

References 1. Landrigan, P.J.; Fuller, R.; Acosta, N.J.R.; Adeyi, O.; Arnold, R.; Basu N.;

492

Baldé, A.; Bertollini,R.; Bose-O’Reilly, S.; Boufford, J.; Breysse, P.; Chiles, T.;

493

Mahidol, C.; Coll-Seck, A.;Cropper, ML.;Fobil, J.;Fuster, V.; Greenstone, M.;

494

Haines, A.; Hanrahan, D.; Hunter, D.; Khare, M.;Krupnick, A.; Lanphear,

495

B.;Lohani, B.;Martin, K.;Mathiasen, K.; McTeer, M.; Murray, C.; Ndahimananjara,

496

J.;Perera, F.;Potočnik, J.;Preker, A.; Ramesh, J.; Rockström, J.; Salinas, C.;

497

Samson, L.; Sandilya, K.; Sly, P.; Smith, K.; Steiner, A.; Stewart, R.; Suk,W.; van

498

Schayck, O.;Yadama, G.; Yumkella, K.; Zhong, M. The Lancet commission on

499

pollution and health. Lancet 2018, 391, 462-512.

500 501 502

2. Kelly, F.J. Oxidative stress: its role in air pollution and adverse health effects. Occup.

Environ. Med. 2003, 60(8), 612–616.

503

3. Borm, P.J.A.; Kelly, F.; Kunzli, N.; Schins, R.F.P.; Donaldson, K. Oxidant generation

504

by particulate matter: from biologically effective dose to a promising, novel

505

metric. Occup. Environ. Med. 2007, 64, 73-74.

ACS Paragon Plus Environment

30

Page 31 of 50

506

Environmental Science & Technology

4. Li, N.; Hao, M.; Phalen, R.F.; Hinds, W.C.; Nel, A.E. Particulate air pollutants and

507

asthma: a paradigm for the role of oxidative stress in PM-induced adverse health

508

effects. Clin. Immunol. 2003, 109(3), 250–265.

509 510 511

5. Weichenthal, S.; Godri-Pollitt, K.; Villeneuve, P.J. PM2.5, oxidant defence and cardiorespiratory health: a review. Environ Health 2013, 12, 40. 6. Ayres, J.G.; Borm, P.; Cassee, F.R.; Castranova, V.; Donaldson, K.; Ghio, A.;

512

Harrison, R.M.; Hider, R.; Kelly, F.; Kooter, I.M.; Marano, R.L.; Mudway, I.; Nel,

513

A.; Sioutas, C.; Smith, S.; Baeza-Squiban, A.; Cho, A.; Duggan, S.; Froines, J.

514

Evaluating the toxicity of airborne particulate matter and nanoparticles by

515

measuring oxidative stress potential-a workshop report and concensus

516

statement. Inhal, Toxicol. 2008. 20(1), 75–99.

517

7. Kunzli, N,; Mudway, I.S.; Gotschi, T.; Shi, T.; Kelly, F.J.; Cook, S.; Burney, P.;

518

Foresberg, B.; Gauderman, J.W.; Hazenkamp, M.E.; Heinrich, J.; Jarvis, D.;

519

Norback, D.; Payo-Losa, F.; Poli, A.; Sunyer, J.; Borm, P.J.A. Comparison of

520

oxidative properties, light absorbance, and total and elemental mass

ACS Paragon Plus Environment

31

Environmental Science & Technology

Page 32 of 50

521

concentration of ambient PM2.5 collected at 20 European Sites. Environ. Health

522

Perspect. 2006, 114, 684-690.

523

8. Godri, K.J.; Duggan, S.T.; Fuller, G.W.; Baker, T.; Green, D.; Kelly, F.J.; Mudway, I.S.

524

Particulate matter oxidative potential from waste transfer station activity. Environ.

525

Health Perspect. 2010, 118, 493-498.

526

9. Boogaard, H,; Janssen, N.A.; Fischer, P.H.; Kos, G.P.; Weijers, E.P.; Cassee, F.R.;

527

van der Zee, S.C.; de Hartog, J.J.; Brunekreef, B.; Hoek, G. Contrasts in

528

oxidative potential and other particulate matter characteristics collected near

529

major streets and background locations. Environ. Health Perspect. 2012, 120,

530

185-191.

531

10. Janssen, N.A.; Yang, A.; Strak, M.; Steenhof, M.; Hellack, B.; Gerlofs-Nijland, M.E.;

532

Kuhlbusch, T.; Kelly, F.; Harrison, R.; Brunekreef, B.; Hoek, G.; Cassee, F.

533

Oxidative potential of particulate matter collected at sites with different source

534

characteristics. Sci. Total Environ. 2014, 472, 572-581.

535 536

11. Bates, J.T.; Weber, R.J.; Abrams, J.; Verma, V.; Fang, T.; Klein, M.; Strickland, M.J.; Ebelt Sarnat, S.; Chang, H.H.; Mulholland, J.A.; Tolbert, P.E.; Russell, A.G.

ACS Paragon Plus Environment

32

Page 33 of 50

Environmental Science & Technology

537

Reactive oxygen species generation linked to sources of atmospheric particulate

538

matter and cardiorespiratory effects. Environ. Sci. Technol. 2015, 49, 13605-

539

13612.

540

12. Yang A, Wang M, Beelen R, Dons E, Leseman DL, Brunekreef B, Cassee FR,

541

Janssen NA, Hoek G. Spatial variation and land use regression modeling of the

542

oxidative potential of fine particles. Environ. Health Perspect. 2015a, 123, 1187-

543

1192.

544

13. Maikawa, C.L.; Weichenthal, S.; Wheeler, A.; Dobbin, N.A.; Smargiassi, A.; Evans,

545

G.; Liu, L.; Goldberg, M.S.; Pollitt, K.J. Particulate oxidative burden as a predictor

546

of exhaled nitric oxide in children with asthma. Environ. Health Perspect. 2016,

547

124, 1616-1622.

548

14. Weichenthal, S.; Lavigne, E.; Evans, G.J.; Godri Pollitt, K.J.; Burnett, R.T. Fine

549

particulate matter and emergency room visits for respiratory illness: effect

550

modification by oxidative potential. Am. J. Respir. Crit. Care Med. 2016a, 194,

551

577-586.

ACS Paragon Plus Environment

33

Environmental Science & Technology

552

Page 34 of 50

15. Calas, A.; Uzu, G.; Kelly, F.J.; Houdier, S.; Martins, J.M.F.; Thomas, F.; Molton, F.;

553

Charron, A.; Dunster, C.; Oliete, A.; Jacob, V.; Besombes, J.L.; Chevrier, F.;

554

Jaffrezo, J.L. Comparison between five acellular oxidative potential measurement

555

assays performed with detailed chemistry on PM10 samples from the city of

556

Chamonix (France). Atmos. Chem. Phys. Discuss. 2018, 18, 7863-7875.

557

16. Yang, A.; Jedynska, A.; Hellack, B.; Kooter, I.; Hoek, G.; Brunekreef, B.; Kuhlbusch,

558

T.A.J.; Cassee, F.R.; Janssen, N.A.H. Measurement of the oxidative potential of

559

PM2.5 and its constituents: the effect of extraction solvent and filter type. Atmos.

560

Environ. 2014, 83, 35-42.

561

17. Yang, A.; Hellack, B.; Leseman, D.; Brunekreef, B.; Kuhlbusch, T.A.J.; Cassee,

562

F.R.; Hoek, G.; Janssen, N.A.H. Temporal and spatial variation of the metal-

563

related oxidative potential of PM2.5 and its relation to PM2.5 mass and elemental

564

composition. Atmos. Environ. 2015c, 102, 62-69.

565

18. Yanosky, J.D.; Tonne, C.C.; Beevers, S.D.; Wilkinson, P.; Kelly, F.J. Modeling

566

exposures to the oxidative potential of PM10. Environ Sci Technol 2012, 46,

567

7612-7620.

ACS Paragon Plus Environment

34

Page 35 of 50

568

Environmental Science & Technology

19. Jedynska, A.; Hoek, G.; Wang, M.; Yang, A.; Eeftens, M.; Cyrys, J.; Keuken, M.;

569

Ampe, C.; Beelen, R.; Cesaroni, G.; Forastiere, F.; Cirach, M.; de Hoogh, K.; De

570

Nazelle, A.; Nystad, W.; Makarem Akhlaghi, H.; Declercq, C.; Stempfelet, M.;

571

Eriksen, K.T.; Dimakopoulou, K.; Lanki, T.; Meliefste, K.; Nieuwenhuijsen, M.; Yli-

572

Tuomi, T.; Raaschou-Nielsen, O.; Janssen, N.A.H.; Brunekreef, B.; Kooter, I.M.

573

Spatial variations and development of land use regression models of oxidative

574

potential in ten European study areas. Atmos. Environ. 2017, 150, 240-32.

575

20. Gulliver, J.; Morley, D.; Dunster, C.; McCrea, A.; van Nunen, E.; Tsai, M.Y.; Probst-

576

Hensch, N.; Eeftens, M.; Imboden, M.; Ducret-Stich, R.; Naccarati, A.; Galassi,

577

C.; Ranzi, A.; Nieuwenhuijsen, M.; Curto, A.; Donaire-Gonzalez, D.; Cirach, M.;

578

Vermeulen, R.; Vineis, P.; Hoek, G.; Kelly, F.J. Land use regression models for

579

the oxidative potential of fine particles (PM2.5) in five European areas. Environ.

580

Res. 2018, 160, 247-255.

581 582

21. Weichenthal, S.; Crouse, D.L.; Pinault, L.; Godri-Pollitt, K.; Lavigne, E.; Evans, G.; van Donkelaar, A.; Martin, R.V.; Burnett, R.T. Oxidative burden of fine particulate

ACS Paragon Plus Environment

35

Environmental Science & Technology

583

air pollution and risk of cause-specific mortality in the Canadian Census Health

584

and Environment Cohort (CanCHEC). Environ. Res. 2016b; 146: 92-99.

585

Page 36 of 50

22. Weichenthal, S.; Lavigne, E.; Evans, G.; Pollitt, K.; Burnett, R.T. Ambient PM2.5 and

586

risk of emergency room visits for myocardial infarction: impact of regional PM2.5

587

oxidative potential: a case-control study. Environ. Health 2016c, 15, 46.

588

23. Weichenthal, S.; Shekarrizfard, M.; Kulka, R.; Lakey, P.S.J.; Al-Rijleh, K.; Anowar,

589

S.; Shiraiwa, M.; Hatzopoulou, M. Spatial variations in the estimated production

590

of reactive oxygen species in the epithelial lung lining fluid by PM2.5 iron and

591

copper. Environ. Epidemiol. 2018, 2, e020.

592

24. Baker, M.A.; Cerniglia, G.J.; Zaman, A. Microtiter Plate Assay for the Measurement

593

of Glutathione and Glutathione Disulfide in Large Numbers of Biological

594

Samples. Anal. Biochem. 1990, 190, 360–365.

595

25. Cho, A.K.; Sioutas, C.; Miguel, A.H.; Kumagai, Y.; Schmitz, D.A.; Singh, M.;

596

Eiguren-Fernandez, A.; Froines, J.R. Redox Activity of Airborne Particulate

597

Matter at Different Sites in the Los Angeles Basin. Environ. Res. 2005, 99, 40–

598

47.

ACS Paragon Plus Environment

36

Page 37 of 50

599

Environmental Science & Technology

26. Bell, M.L.; Dominici, F.; Ebisu, K.; Zeger, S.L.; Samet, J.M. Spatial and temporal

600

variations in PM2.5 chemical composition in the United States for Health Effects

601

Studies. Environ. Health Perspect. 2007, 115, 989-995.

602

27. Di Carlo, P.; Brune, W.H.; Martinez, M.; Harder, H.; Lesher, R.; Ren, X.; Thornberry,

603

T.; Carroll, M.A.; Young, V.; Shepson, P.B.; Riemer, D.; Apel, E.; Campbell, C.

604

Missing OH reactivity in a forest: evidence for unknown reactive biogenic VOCs.

605

Science 2004, 304, 722-725.

606

28. Bates, J.T.; Weber, R.J.; Verma, V.; Fang, T.; Ivey, C.; Liu, C.; Sarnat, S.E.; Chang,

607

H.H.; Mulholland, J.A.; Russell, A. Source impact modeling of spatiotemporal

608

trends in PM2.5 oxidative potential across the eastern United States. Atmos.

609

Environ. 2018, 193, 158-167.

610

30. Wang, Y.; Hopke, P.K.; Rattigan, O.V.; Xia, X.; Chalupa, D.C.; Utell, M.J.

611

Characterization of residential wood combustion particles using the two-wavelength

612

aethalometer. Environ Sci Technol 2011; 45: 7387-7393.

613

31. Kurmi, O.P.; Dunster, C.; Ayres, J.G.; Kelly, F.J. Oxidative potential of smoke from

614

burning wood and mixed biomass fuels. Free. Rad. Res. 2013, 47, 829-835.

ACS Paragon Plus Environment

37

Environmental Science & Technology

Page 38 of 50

615

32. Verma, V., Fang, T., Guo, H., King, L., Bates, J. T., Peltier, R. E., Edgerton, E.,

616

Russell, A. G., and Weber, R. J. Reactive oxygen species associated with water-soluble

617

PM2.5 in the southeastern United States: spatiotemporal trends and source

618

apportionment. Atmos Chem Phys 2014, 14, 12915-12930.

ACS Paragon Plus Environment

38

Page 39 of 50

Environmental Science & Technology

Table 1. Descriptive Statistics for PM2.5, black carbon, and PM2.5 oxidative potential in Toronto, Canada Pollutant

n

Mean

SD

Min

5th

25th

50th

75th

95th

Max

PM2.5 (g/m3)

67

6.41

0.78

4.85

5.16

5.90

6.32

6.84

7.82

8.53

OPAA (%Depletion/m3)

66

2.00

0.54

1.03

1.30

1.71

1.96

2.10

2.83

4.72

OPAA (pmol/min/m3)

66

18.6

5.02

9.58

12.1

15.9

18.2

19.5

26.3

43.9

OPGSH (%Depletion/m3)

66

2.29

2.12

0.0840

0.441

0.698

1.74

3.08

5.99

13.8

OPGSH (pmol/min/m3)

66

21.3

19.7

0.781

4.10

6.49

16.2

28.6

55.7

128

OPDTT (pmol/min/m3)

66

36.7

11.9

20.5

23.4

29.5

35.0

40.1

62.3

88.9

BC880nm (ng/m3)

67

1965

534

653

1194

1652

1933

2168

2753

4187

BC370nm (ng/m3)

67

1079

306

289

585

918

1060

1223

1545

2227

PM2.5 (g/m3)

42

5.32

0.88

4.01

4.34

4.80

5.13

5.86

6.38

8.80

OPAA (%Depletion/m3)

40

1.71

0.532

0.844

0.972

1.34

1.60

1.99

2.70

3.26

OPAA (pmol/min/m3)

40

15.9

4.95

7.85

9.04

12.5

14.9

18.5

25.1

30.3

OPGSH (%Depletion/m3)

40

2.66

2.46

0.336

0.376

0.759

2.14

3.58

8.01

11.7

OPGSH (pmol/min/m3)

40

24.7

22.9

3.12

3.50

7.06

19.9

33.3

74.5

109

OPDTT (pmol/min/m3)

40

27.2

13.2

7.59

10.0

18.8

26.0

31.7

51.5

80.8

BC880nm (ng/m3)

42

1197

319

519

766

1069

1157

1290

1593

2638

BC370nm (ng/m3)

42

1024

218

477

709

882

1035

1105

1339

1778

29

5.77

0.76

4.73

4.73

5.33

5.65

6.20

7.49

8.06

Summer

Winter

Annual* PM2.5 (g/m3)

ACS Paragon Plus Environment

39

Environmental Science & Technology

Page 40 of 50

OPAA (%Depletion/m3)

26

1.81

0.459

1.32

1.33

1.56

1.75

1.93

2.40

3.50

OPAA (pmol/min/m3)

26

16.8

4.27

12.3

12.4

14.5

16.3

17.9

22.3

32.5

OPGSH (%Depletion/m3)

26

2.08

1.03

0.497

0.733

1.30

1.87

2.54

4.08

4.23

OPGSH (pmol/min/m3)

26

19.3

9.58

4.62

6.82

12.1

17.4

23.6

37.9

39.3

OPDTT (pmol/min/m3)

26

30.6

10.7

16.9

17.7

25.7

28.5

34.1

46.3

70.1

BC880nm (ng/m3)

29

1497

450

940

961

1233

1470

1609

1919

3412

BC370nm (ng/m3)

29

1000

247

695

700

837

993

1081

1250

2002

*Annual data reflects those sites with both summer and winter data

ACS Paragon Plus Environment

40

Page 41 of 50

Environmental Science & Technology

Figure 1. Spatial correlations between PM2.5 oxidative potential (OPGSH, OPAA, OPDTT), PM2.5 metals, BC880nm, and BC370nm during summer (A) and winter (B) A

B

ACS Paragon Plus Environment

41

Page 42 of 50

As Rb Pb Zn Sb Cr Co Ni Cs Zr Fe Ba Cu Cl Na Mn Sr Ti Br BC.880 BC.370 K PM2.5 S P Mg Ca Al Si Se V OP.DTT OP.GSH OP.AA

Environmental Science & Technology

As Rb Pb Zn Sb Cr Co Ni Cs Zr Fe Ba Cu Cl Na Mn Sr Ti Br BC.880 BC.370 K PM2.5 S P Mg Ca Al Si Se V OP.DTT OP.GSH OP.AA

0.09 −0.11 0.17 −0.28 0.21

1 0.09

1

0.1

0.15

0.1

−0.12 −0.12 −0.15 −0.19 −0.16 −0.07 0.15 −0.08 −0.14 −0.17 −0.27 −0.19 −0.15 0.17 −0.29 −0.35 −0.16 −0.22 −0.13 −0.15 −0.11 −0.21 −0.28 −0.09 −0.14 −0.19

0.56

0.04 −0.03 0.15

0.01

0.08

0.1

0.16

0.21

0.1

0.09

0.37

0.19

0.18

0.17

0.11

0.2

0.36

0.19

0.25

0.26 −0.01 −0.21 −0.2 −0.13 −0.02

0.69

0.05

0.05

0.14

0.21

0.4

0.19

0.1

0.08

0.13

0.08

0.07

0.01

0.03 −0.12 0.05

0.15

0.08

0.1

0.03

−0.11

0.1

1

0.7

0.17

0.28

0.7

1

0

0

1

−0.28 −0.19 −0.11

−0.11 −0.01 0.33

0.15 −0.01 0.04 −0.22

0.21

0.04

0.4

0.02

0.07

0.16

0.11

0.07 −0.07 −0.01 −0.04

−0.22 −0.01 −0.15 0.02 −0.08 −0.11 −0.04 −0.06 0.19 1

0

0

0.07

−0.1 −0.17 −0.32 −0.18 −0.1

0.33

0.63

0.61

0.6

0.64

0.62

0.42

0.58

0.51

0.55

0.55

0.05

0.43

0.34

0.38

0.2

0.31

0.27

0.39

0.36

0.31

0.27

0.15

0.08

0.08

0.33

0.25 0.32

1

0.69

0.64

0.55

0.73

0.69

0.72

0.62

0.61

0.69

0.68

0.68

0.32

0.63

0.54

0.24

0.39

0.33

0.22

0.6

0.61

0.61

0.53

0.06

0.13 −0.04 0.18

0.69

1

0.43

0.48

0.61

0.49

0.57

0.54

0.64

0.62

0.57

0.57

0.38

0.55

0.46

0.3

0.29

0.29

0.35

0.53

0.53

0.56

0.5

0.01

0.15 −0.09 −0.04 0.12

0.05

0.63

0.64

0.43

1

0.85

0.83

0.92

0.89

0.75

0.7

0.71

0.79

0.88

0.39

0.73

0.61

0.22

0.43

0.38

0.32

0.55

0.52

0.45

0.4

0.14

0.08

0.06 −0.12 −0.03 0.02 −0.08 0.61

0.55

0.48

0.85

1

0.91

0.96

0.93

0.69

0.64

0.78

0.86

0.89

0.28

0.75

0.59

0.23

0.4

0.36

0.4

0.6

0.54

0.5

0.46

0.25

0.07 −0.15 0.15

0.16 −0.11

0.6

0.73

0.61

0.83

0.91

1

0.95

0.97

0.83

0.74

0.93

0.94

0.95

0.44

0.84

0.68

0.4

0.55

0.48

0.5

0.81

0.76

0.75

0.7

0.2

0.11 −0.19 0.01

0.05 −0.04 0.64

0.69

0.49

0.92

0.96

0.95

1

0.98

0.77

0.7

0.83

0.89

0.93

0.35

0.8

0.64

0.27

0.48

0.43

0.4

0.66

0.62

0.57

0.52

0.1

0.14

0.56

0.15 −0.12 0.04

0.4

−0.01 0.49

0.02

0.15

0.36

0.55

0.28

0.3

0.27

0.46

0.36

0.21

0.22

0.5

0.23

0.29

0.24

0.3

0.51 0.49

0.07 −0.16 0.08

0.14 −0.06 0.62

0.72

0.57

0.89

0.93

0.97

0.98

1

0.82

0.76

0.87

0.92

0.96

0.34

0.81

0.66

0.29

0.46

0.4

0.42

0.7

0.67

0.65

0.58

0.16

0.29

0.18

0.21

−0.07 −0.07

0.1

0.21

0.19

0.42

0.62

0.54

0.75

0.69

0.83

0.77

0.82

1

0.89

0.8

0.83

0.84

0.37

0.76

0.63

0.41

0.49

0.28

0.5

0.66

0.78

0.74

0.66

0.02

0.25

0.07

0.08

0.4

−0.01 0.15

0.16

0.4

0.12

0.58

0.61

0.64

0.7

0.64

0.74

0.7

0.76

0.89

1

0.74

0.73

0.74

0.25

0.63

0.53

0.35

0.31

0.13

0.43

0.48

0.62

0.57

0.47

0.04

0.13 −0.06 0.01

0.24

0.19 −0.11 0.51

0.69

0.62

0.71

0.78

0.93

0.83

0.87

0.8

0.74

1

0.89

0.89

0.47

0.81

0.66

0.5

0.63

0.46

0.58

0.83

0.84

0.82

0.8

0.21

0.37

0.17

0.16

0.5

0.1

−0.04 −0.08 0.21

0.55

0.68

0.57

0.79

0.86

0.94

0.89

0.92

0.83

0.73

0.89

1

0.94

0.44

0.86

0.73

0.42

0.51

0.45

0.57

0.8

0.78

0.77

0.72

0.19

0.31

0.27

0.24

0.52

0.01 −0.17 0.09

0.08 −0.04 0.55

0.68

0.57

0.88

0.89

0.95

0.93

0.96

0.84

0.74

0.89

0.94

1

0.37

0.84

0.7

0.38

0.53

0.44

0.5

0.74

0.78

0.72

0.69

0.1

0.26

0.22

0.21

0.55

−0.07 −0.27 0.37

0.13

0.05

0.32

0.38

0.39

0.28

0.44

0.35

0.34

0.37

0.25

0.47

0.44

0.37

1

0.57

0.55

0.31

0.4

0.53

0.41

0.53

0.37

0.49

0.47

0.19

0.21

0.02

0.24

0.21

−0.03 −0.19 0.19

0.08 −0.01 0.43

0.63

0.55

0.73

0.75

0.84

0.8

0.81

0.76

0.63

0.81

0.86

0.84

0.57

1

0.95

0.51

0.51

0.39

0.5

0.74

0.71

0.71

0.69

0.27

0.36

0.26

0.24

0.39

−0.01 −0.15 0.18

0.07

0.34

0.54

0.46

0.61

0.59

0.68

0.64

0.66

0.63

0.53

0.66

0.73

0.7

0.55

0.95

1

0.45

0.38

0.29

0.42

0.6

0.57

0.6

0.58

0.26

0.3

0.26

0.28

0.34

0.01 −0.14 0.38

0.24

0.3

0.22

0.23

0.4

0.27

0.29

0.41

0.35

0.5

0.42

0.38

0.31

0.51

0.45

1

0.39

0.45

0.43

0.56

0.62

0.57

0.69 −0.07 0.14

0.03

0.06

0.06

0.03

0

−0.14

0.1

0.01

0.8

0.12 −0.11 0.01 −0.04 0.07 −0.01 0.04 −0.14 0.07 −0.16 −0.04 −0.21 −0.06 −0.12 −0.16 0.03 −0.16 0.07 −0.04 −0.05

0.49

0.69 −0.15 0.33

0.23 −0.03 0.33

1

0.14

0.14

0.06

0.01 −0.07 −0.03 −0.01

−0.07 0.15 −0.13 −0.05 −0.11 −0.17 −0.17 −0.43 −0.13 −0.13 0.34

0.23

0.28 −0.19 0.15 −0.03

0.6

0.4

0.2

0

0

0.17

0.17

0.07

0.04

0.29

0.43

0.4

0.55

0.48

0.46

0.49

0.31

0.63

0.51

0.53

0.4

0.51

0.38

0.39

1

0.45

0.48

0.57

0.61

0.58

0.59

0.25

0.35

0.27

0.21

0.56

0.33

0.29

0.38

0.36

0.48

0.43

0.4

0.28

0.13

0.46

0.45

0.44

0.53

0.39

0.29

0.45

0.45

1

0.48

0.57

0.5

0.49

0.53

0

0.23

0.07

0.44

0.31

0.2

0.05 −0.04 0.27

0.22

0.35

0.32

0.4

0.5

0.4

0.42

0.5

0.43

0.58

0.57

0.5

0.41

0.5

0.42

0.43

0.48

0.48

1

0.6

0.66

0.61

0.64

0.08

0.35

0.05

0.19

0.33

−0.05 −0.22 0.36

0.15 −0.21 0.39

0.6

0.53

0.55

0.6

0.81

0.66

0.7

0.66

0.48

0.83

0.8

0.74

0.53

0.74

0.6

0.56

0.57

0.57

0.6

1

0.86

0.84

0.86

0.19

0.3

0.17

0.2

0.48

−0.11 −0.13 0.19

0.08 −0.06 0.36

0.61

0.53

0.52

0.54

0.76

0.62

0.67

0.78

0.62

0.84

0.78

0.78

0.37

0.71

0.57

0.62

0.61

0.5

0.66

0.86

1

0.92

0.94

0.07

0.34

0.09

0.06

0.38

−0.17 −0.15 0.25

0.1

−0.12 0.31

0.61

0.56

0.45

0.5

0.75

0.57

0.65

0.74

0.57

0.82

0.77

0.72

0.49

0.71

0.6

0.57

0.58

0.49

0.61

0.84

0.92

1

0.96

0.14

0.45

0.14

0.1

0.4

−0.17 −0.11 0.26

0.03 −0.16 0.27

0.53

0.5

0.4

0.46

0.7

0.52

0.58

0.66

0.47

0.8

0.72

0.69

0.47

0.69

0.58

0.69

0.59

0.53

0.64

0.86

0.94

0.96

1

0.08

0.41

0.16

0.06

0.36

0.15

0.06

0.01

0.14

0.25

0.2

0.23

0.16

0.02

0.04

0.21

0.19

0.1

0.19

0.27

0.26 −0.07 0.25

0

0.08

0.19

0.07

0.14

0.08

1

0.38

0.32

0.04

0.16

−0.13 −0.28 −0.21 −0.17 −0.16 0.08

0.13

0.15

0.08

0.28

0.36

0.29

0.29

0.25

0.13

0.37

0.31

0.26

0.21

0.36

0.3

0.14

0.35

0.23

0.35

0.3

0.34

0.45

0.41

0.38

1

0.08 −0.04 −0.09 0.15

0.3

0.21

0.24

0.18

0.07 −0.06 0.17

0.27

0.22

0.02

0.26

0.26

0.03

0.27

0.07

0.05

0.17

0.09

0.14

0.16

0.32

0.26

−0.07 −0.29 0.11

0.2

0.39

0.15 −0.35 0.07 −0.12 −0.16 0.31 −0.13 −0.16

−0.43 −0.21 −0.01 −0.1

0.07

0.03

−0.13 −0.09 −0.2 −0.32 0.07

0.26 −0.02 0.15 1

0.34 −0.14 −0.13 −0.18 −0.04 0.33

0.18 −0.04 0.36

0.27

0.22

0.3

0.21

0.08

0.01

0.16

0.24

0.21

0.24

0.24

0.28

0.06

0.21

0.44

0.19

0.2

0.06

0.1

0.06

0.04 −0.02 0.34

0.14 −0.19 −0.02 −0.1 −0.05 0.25

0.32

0.46

0.5

0.51

0.49

0.4

0.24

0.5

0.52

0.55

0.21

0.39

0.34

0.06

0.56

0.31

0.33

0.48

0.38

0.4

0.36

0.16

0.12

0.55

0.15

0.57

0.34

0.57

1

0.58

0.58

1

−0.2

−0.4

−0.6

−0.8

−1

Figure 2. Distributions of

OPDTT (panel

A) and

OPGSH

and

OPAA

(panel B) in Toronto,

Canada for sites with summer and winter data (n=26) A

OP−DTT: Winter

OP−DTT: Summer

OP−DTT: Annual 10

20

30

40

50

60

70

80

3 pmol min m

B

ACS Paragon Plus Environment

42

Page 43 of 50

Environmental Science & Technology

ACS Paragon Plus Environment

43

Environmental Science & Technology

Page 44 of 50

Figure 3. Between-season differences in OPGSH (panel A), OPAA (panel B), OPDTT (panel C), and PM2.5 (panel D) by site. Sites with positive values had higher concentrations during summer; sites with negative values had higher concentrations during winter.

ACS Paragon Plus Environment

44

Page 45 of 50

Environmental Science & Technology

Figure 4. Between-site variations in OPGSH, OPAA, and OPDTT during summer (left) and winter (right) in Toronto, Canada. Green bars indicate sites above the mean value (Z=0) and red bars indicate sites below the mean.

ACS Paragon Plus Environment

45

Environmental Science & Technology

Summer

Page 46 of 50

Winter

ACS Paragon Plus Environment

46

Page 47 of 50

Environmental Science & Technology

Figure 5. Multivariable models for OPAA (% Depletion/m3), OPGSH (% Depletion/m3), and OPDTT (pmol/min/m3) during summer

ACS Paragon Plus Environment

47

Environmental Science & Technology

Page 48 of 50

Figure 6. Multivariable models for OPAA (% Depletion/m3), OPGSH (% Depletion/m3), and OPDTT (pmol/min/m3) during winter

ACS Paragon Plus Environment

48

Page 49 of 50

Environmental Science & Technology

Figure 7. Predicted spatial distributions of OPAA (%Depletion/m3) and OPDTT (pmol/min/m3) during summer in Toronto, Canada (ArcMap 10.4.1; ESRI, Redlands, CA) Legend

A)

Summer OPAA

B)

Summer OPDTT

TorontoGrid_OPnew ln(OPAA ln(OPAA Summer) < 0.04 Summer) 0.04 - 0.11 0.12 - 0.19 0.20 - 0.26 0.27 - 0.34 0.35 - 0.42 0.43 - 0.50 0.51 - 0.57 0.58 - 0.65 0.66 - 0.73 0.74 - 0.81 0.82 - 0.89 0.90 - 0.96 0.97 - 1.04 > 1.04

ln(OPDTT Legend

ACS Paragon Plus Environment

Tor_OPnew summer Summer) ln_op_dtt_per_min_per_m3_S < 3.16 3.16 - 3.23 3.24 - 3.30 3.31 - 3.38 3.39 - 3.45 3.46 - 3.53 3.54 - 3.61 3.61 - 3.68 49 3.69 - 3.75 3.76 - 3.83 3.84 - 3.91

Environmental Science & Technology

Page 50 of 50

TOC ART

OP−DTT: Winter

OP−DTT: Summer

OP−DTT: Annual 10

20

30

40

50

60

70

80

pmol min m3

ACS Paragon Plus Environment

50