Spatial Variation of Aerosol Chemical Composition and Organic

Aug 3, 2015 - Positive matrix factorization (PMF) applied to organic aerosol (OA) data yielded 5 factors: hydrocarbon-like OA (HOA), cooking OA (COA),...
1 downloads 6 Views 2MB Size
Subscriber access provided by The Univ of Iowa Libraries

Article

Spatial Variation of Aerosol Chemical Composition and Organic Components Identified by Positive Matrix Factorization in the Barcelona Region Claudia Mohr, Peter DeCarlo, Maarten F. Heringa, Roberto Chirico, Rene Richter, Monica Crippa, Xavier Querol, Urs Baltensperger, and Andre Prevot Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b02149 • Publication Date (Web): 03 Aug 2015 Downloaded from http://pubs.acs.org on August 11, 2015

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 30

Environmental Science & Technology

1

Spatial Variation of Aerosol Chemical Composition and Organic

2

Components Identified by Positive Matrix Factorization in the Barcelona

3

Region

4 5

Claudia Mohr1,2,*, Peter. F. DeCarlo1,3, Maarten. F. Heringa1, Roberto Chirico1,4, René Richter1, Monica Crippa1,5, Xavier Querol6, Urs Baltensperger1, and André S. H. Prévôt1

6 7

1

8 9

2

10 11

3

12 13

4

5

14 15 16 17 18 19 20 21 22

6

Laboratory of Atmospheric Chemistry, Paul Scherrer Institute (PSI), Villigen, Switzerland

now at: Institute for Meteorology and Climate Research, Karlsruhe Institute of Technology, Karlsruhe, Germany

now at: Department of Civil, Architectural, and Environmental Engineering and Department of Chemistry, Drexel University, Philadelphia, PA, USA now at: Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA), UTAPRAD-DIM, Frascati, Italy

now at: European Commission Joint Research Centre (JRC), Institute for Environment and Sustainability, Ispra, Italy

Institute for Environmental Assessment and Water Research (IDAEA-CSIC), Barcelona, Spain *Corresponding author: Claudia Mohr, Karlsruhe Institute of Technology, Hermannvon-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany Phone: +49 -721 / 608 - 2 2696 Fax: +49 -721 / 608 - 2 4332 [email protected]

23

ACS Paragon Plus Environment

1

Environmental Science & Technology

24

Page 2 of 30

Abstract

25

The spatial distribution of PM1 components in the Barcelona metropolitan area

26

was investigated using on-road mobile measurements of atmospheric particle- and gas-

27

phase compounds during the DAURE campaign in March 2009. Positive matrix

28

factorization (PMF) applied to organic aerosol (OA) data yielded 5 factors: hydrocarbon-

29

like OA (HOA), cooking OA (COA), biomass burning OA (BBOA), and low volatility

30

and semi-volatile oxygenated OA (LV-OOA and SV-OOA). The area under investigation

31

(~ 500 km2) was divided into six zones (city center, harbor, industrial area, pre-coastal

32

depression, 2 mountain ranges) for measurements and data analysis. Mean zonal OA

33

concentrations are 4.9 - 9.5 µg m-3. The area is heavily impacted by local primary

34

emissions (HOA 14 – 38%, COA 10 – 18%, BBOA 10 – 12% of OA); concentrations of

35

traffic-related components, especially black carbon, are biased high due to the on-road

36

nature of the measurements. The formation of secondary OA adds more than half of the

37

OA burden outside the city center (SV-OOA 14 – 40%, LV-OOA 17 – 42 % of OA). A

38

case study of one measurement drive from the shore to the pre-coastal mountain range

39

furthest downwind of the city center indicates the importance of non-fossil over

40

anthropogenic secondary OA based on OA/CO.

41

ACS Paragon Plus Environment

2

Page 3 of 30

Environmental Science & Technology

42

Introduction

43

Adverse health effects of atmospheric aerosols are of special concern in highly populated

44

areas 1. The identification and quantification of their sources, basis for the design of

45

mitigation strategies, remains challenging due to their complex nature. The organic

46

fraction of atmospheric aerosols (organic aerosol, OA), often the major component in

47

submicron particulate matter (PM1) 2, is composed of 1000s of compounds

48

chemical reactions within the particles

49

condensation, evaporation, gas-phase oxidation, and recondensation processes 6-8. Recent

50

advances in the apportionment of OA sources have been achieved through the

51

combination of aerosol mass spectrometer (AMS) measurements

52

analytical model positive matrix factorization (PMF)

53

has led to the identification of diverse OA components: Hydrocarbon-like OA (HOA,

54

related to fresh traffic emissions), biomass burning OA (BBOA, from domestic wood

55

burning or wildfire emissions), cooking OA, and different types of oxygenated OA

56

(OOA, related to secondary OA), i.e., low volatility OOA (LV-OOA), characterized as a

57

regional, heavily aged OA, and the less aged semi-volatile OOA (SV-OOA) 11, 12 11, 13-16.

58

Most of these studies were based on stationary measurements during a certain period of

59

time. However, the evolution of pollutants in the plume of a major source region such as

60

an urban center or the characterization of the sources and processes of atmospheric

61

pollutants of a whole region are best assessed by mobile measurements, ground based or

62

using aircraft 17-25.

4, 5

3

subject to

and the surrounding gas phase via

10

9

and the factor

. The application of this method

63

The densely populated metropolitan area of Barcelona, Spain (5.4 million

64

inhabitants) is subject to very high air pollution levels due to the combination of its

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 30

65

Mediterranean climate, abundant biogenic and anthropogenic emissions, orography and

66

atmospheric dynamics 26, 27. Limited to the SE by the Mediterranean Sea and to the NW-

67

SW by the Catalan coastal ranges, Barcelona experiences a daily cycle of nocturnal

68

offshore flows, transporting air masses from the surrounding valleys (e. g. the densely

69

populated Vallés) to the city center, and diurnal sea breeze, advecting polluted coastal air

70

masses inland 28. The Llobregat river valley SW of the city center contains both extensive

71

industrial and agricultural areas.

72

This paper investigates the sources, processes and evolution of PM1 in the

73

Barcelona metropolitan area during late winter (when frequent cyclonic episodes further

74

enhance pre-existing high PM1 concentrations

75

fraction. Mobile measurements were performed in March 2009 within the international

76

DAURE campaign (Determination of the sources of atmospheric Aerosols in Urban and

77

Rural Environments in the western Mediterranean)

78

measurements of various particle- and gas-phase parameters at an urban background site

79

investigated and quantified PM sources such as road traffic, construction-demolition

80

works, shipping emissions, and photochemical processes

81

burning

82

carbon

83

results from the mobile measurements complement these findings by adding information

84

on the evolution of the Barcelona city plume and the emission sources in the regions

85

surrounding the city. Since measurements were carried out on-road, levels and

86

composition of PM1 are influenced by close traffic emissions from circulating vehicles.

31

32

or cooking

21

29

), with a special focus on the organic

26

30

. Concurrent stationary

. Contributions of biomass

emissions to OA, the fossil and non-fossil fractions of organic

, and the particulate trace metal content

33

were assessed, among others. The

87

ACS Paragon Plus Environment

4

Page 5 of 30

88

Environmental Science & Technology

Experimental Section

89

Mobile measurements. 16 days of mobile measurements (61 sampling hours)

90

were performed in the Barcelona metropolitan area with the Paul Scherrer Institute

91

mobile laboratory, an IVECO Daily van equipped with several continuous online gas-

92

and particle-phase instruments 34, between March 1 and 26, 2009. Figure 1 shows a map

93

of the Barcelona metropolitan area and the driving route. It was divided into 6 different

94

zones for both the measurements and the data analysis: City center (zone 1), harbor (zone

95

2), the industrialized Llobregat valley (zone 3), the coastal mountain range (zone 4), the

96

pre-coastal mountain range (zone 6) and the heavily populated pre-coastal depression

97

Vallés (zone 5). Table 1 summarizes the zones, driving dates and times. For a discussion

98

of zone identification see SI.

99

Instrumentation. Table S1 in the Supporting Information (SI) lists the

100

instruments deployed in the mobile laboratory and gives information on the parameters

101

measured, the detection limits and the time resolution of the data acquisition. Details are

102

provided elsewhere 34; only deviations from this setup are provided here. The focus of the

103

present analysis is on size-resolved chemical composition data of non-refractory PM1

104

measured by the Aerodyne high-resolution time-of-flight aerosol mass spectrometer

105

(AMS). AMS data were acquired in the higher sensitivity V-Mode 35 at a time resolution

106

of 6 s. The PM1 equivalent black carbon fraction (EBC) was derived from aerosol light

107

absorption measurements at 670 nm by a multi-angle absorption photometer (MAAP

108

5012, Thermo), time resolution 1 s. An ozone monitor (Dual Beam, Model 205, 2B

109

Technologies, Inc.) added for this campaign measured O3 concentrations based on the

110

absorption of UV light at 254 nm. NO2 concentrations are available from March 10, 2009

ACS Paragon Plus Environment

35

5

Environmental Science & Technology

Page 6 of 30

34

111

onward (when the chromium trioxide converter used to measure NOx

ceased to

112

function; earlier data, where a separation of NO and NO2 is not possible, are excluded

113

from the analysis). All gases were measured with a time resolution of 1 s.

114

Data Analysis. AMS data were processed using the software SQUIRREL v.1.51B

115

and PIKA v1.1.10B within Igor Pro 6.22A (Wavemetrics). The collection efficiency (CE)

116

was determined based on a) the nitrate fraction fNO3 of the species measured by AMS

117

and b) comparison of the AMS when the mobile laboratory was measuring at the urban

118

background site (location shown in Figure 1) with a scanning mobility particle sizer

119

(SMPS, TSI 3071 type differential mobility analyzer 37 custom built by PSI, condensation

120

particle counter CPC 1022, TSI), a second AMS, and PM1 measured by an optical

121

particle counter (1107, Grimm Labortechnik GmbH & Co. KG) and corrected with PM1

122

sample loadings collected on filters 21. This resulted in a CE = 0.7 for fNO3 ≤ 0.25, CE = 1

123

for fNO3 ≥ 0.78, and a linear increase between those values for 0.25 ≤ fNO3 ≤ 0.78 (Figure

124

S1 in the SI). The CO2+ ion signal was corrected for elevated gas-phase contributions by

125

a time-dependent factor based on gaseous CO2 measurements.

36

126

PMF was applied to organic high resolution (HR) mass spectral data (time

127

resolution 6 s). The input matrices, comprising the complete mobile data set, were

128

prepared as described previously34. All data are in local standard time (UTC + 1) and at

129

local pressure and temperature conditions (10 - 25 ° C and 943 - 1048 hPa).

130

Mobile measurements are usually carried out during short periods of time only

131

(few hours per day during few weeks per year). Analyses of the spatial distribution of

132

parameters based on mobile measurements thus always bear the risk of being overly

ACS Paragon Plus Environment

6

Page 7 of 30

Environmental Science & Technology

133

influenced by meteorological conditions and diurnal patterns 23. This is especially true for

134

the present study where the mean concentration values for the different locations are

135

based on aggregate data from several measurement drives spread out across several

136

weeks. To reduce the influence of the temporal component on the spatial distribution, the

137

time series of each parameter X measured in the mobile laboratory during a single drive

138

(see Table 1) was divided by a normalization factor N based on the concentration of X

139

measured simultaneously at the urban background station. The normalization factor NDy,

140

of species X for an individual drive, Dy, is given as:

141

𝑁𝐷𝑦 =

142

with 𝑋_𝐵𝐺𝐷𝑦 denoting the background station mean value of X for the time interval of Dy

143

and 𝑋_𝐵𝐺 the background station campaign mean value of X. Individual factors for

144

particle components instead of a single factor for all components based on total PM were

145

used since different meteorological conditions (e. g. relative humidity or radiation) can

146

influence individual particle components differently. For PMF factors the Org

147

normalization factor was used (normalization was done after PMF analysis).

148

Normalization factors per drives and species are presented in the SI, as well as a detailed

149

discussion on the influence of the normalization on zonal mean values, and limitations of

150

the method employed. We ask the reader to keep in mind that the mass concentrations

151

reported in Figures 3 and 4 are normalized data. In addition, the data were binned into

152

different time of day intervals to give an indication of the influence of the diurnal pattern

153

of parameters on the spatial variation. Table 1 specifies the time of day bins and the

154

prevailing atmospheric scenarios as defined by Pandolfi et al.

𝑋_𝐵𝐺𝐷𝑦

(1),

𝑋_𝐵𝐺

ACS Paragon Plus Environment

26

which influenced

7

Environmental Science & Technology

Page 8 of 30

155

meteorological conditions, strength of the sea breeze cycle, mixing/residual layer height

156

and thus pollutant concentrations in the region. Briefly, scenario A is characterized by the

157

most severe pollution episodes, recirculation of air masses, the pre-coastal mountain

158

range lying within the mixing and residual layers, and a weakened sea breeze cycle;

159

scenario C by North Atlantic advection flushing the planetary boundary layer; T denotes

160

the transition periods with strong sea breeze cycles and a decoupling of pollutant

161

concentrations between the city center and the pre-coastal mountain range.

162 163

Results and Discussion

164

PMF results. Five OA factors were identified in the mobile dataset, similar to the

165

parallel stationary AMS data from the urban background site 21. These factors are shown

166

in Figure 2, and are identified as the following: Hydrocarbon-like OA (HOA) from traffic

167

(mostly diesel) emissions; biomass burning OA (BBOA) mostly from open agricultural

168

fires

169

secondary OA: Low-volatility OOA (LV-OOA), and semi-volatile OOA (SV-OOA). The

170

R2 values of the correlations of the stationary and the mobile factor profiles, and of the

171

time series of the PMF factors and ancillary data are shown in Table S3 in the SI, as well

172

as the factors’ atomic oxygen to carbon (O:C) and organic matter (OM) to organic carbon

173

(OC) ratios for both solutions. The noisy characteristics of the 6-s mobile data and

174

different magnitudes of the peaks resulted in high scatter in the correlations, thus 5-min

175

averages were used for most of the time series correlations (numbers in italics). PMF was

176

run in the robust mode and the solution was analyzed as outlined in the SI from Mohr et

177

al.

21

31

; cooking OA (COA); and two different oxygenated OA (OOA), related to

, and the 5-factor solution was chosen. The most central rotation (fpeak = 0) and a

ACS Paragon Plus Environment

8

Page 9 of 30

Environmental Science & Technology

178

pseudorandom start value SEED = 23 were selected based on the correlations displayed

179

in Table S3, SI. Figures S4 and S5 in the SI show the mass fractions of the 5-factorial

180

solution as a function of fpeak and SEED values tested, respectively. The ratio of Q

181

(parameter to be minimized by the PMF algorithm

182

freedom of the fitted data is) 1.18 for five factors, indicating a good error estimation12, 38.

183

The 4-factor solution was discarded due to one factor being a combination of COA and

184

BBOA (Figure S6, SI); the 6-factor solution yielded a “residual” factor (Figure S7, SI).

10

) over Qexp (expected Q, degree of

185

Mobile and stationary HOA, COA and LV-OOA factor mass spectra are

186

correlated with R2 values of 0.99, 0.92, and 0.96, respectively (Table S3). The O:C ratios

187

of the mobile and stationary HOA factors are similarly low (0.01 and 0.03, respectively).

188

For the COA, the stationary factor exhibits a higher O:C ratio (0.21) than the mobile

189

factor (0.1), indicating that measurements on-road were closer to the cooking sources

190

with lower atmospheric residence time and lower exposure to oxidants. The mobile LV-

191

OOA factor has a higher O:C ratio (0.98) than the stationary LV-OOA factor (0.75),

192

likely due to mobile measurements including more aged and processed air masses

193

downwind of the city center. The BBOA spectra (R2 = 0.77) exhibit differences in the

194

signals at m/z 44 (and 28), which are higher in the mobile spectrum. This is also reflected

195

in the O:C ratios (0.31 and 0.24 for the mobile and stationary factor profile, respectively).

196

The SV-OOA spectra show the least similarities (R2 = 0.69), mostly due to the signal of

197

the CHO ion (R2 = 0.86 for the linear regression fit excluding CHO). Again, the O:C ratio

198

of the mobile factor is higher than that of the stationary factor (0.54 and 0.32,

199

respectively).

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 30

200

Correlation of EBC with HOA (R2 = 0.63) is better than with BBOA (R2 = 0.24)

201

and is likely due in part to sample bias from on-road sampling and variability in emission

202

of BC from biomass combustion sources 39 (Table S3). This indicates that the majority of

203

EBC sampled originated from traffic related fossil fuel combustions (supported by the R2

204

= 0.71 of the correlation of HOA and CO, and R2 = 0.52 for HOA and NO2). These

205

values are comparable to correlation factors found for measurements of traffic-related

206

pollutants in Queens, New York with a soot particle AMS (R2 = 0.5 and 0.6 for HOA

207

(including refractory BC) and BC or NOx, respectively)

208

the time series of LV-OOA and SO4 are well correlated due to their regional

209

characteristics

210

shipping emissions which occur spatially and temporally separated from urban emissions

211

locally influence the SO4 concentrations 41, 42. LV-OOA exhibits a temporal evolution (R2

212

= 0.40) similar to e.g. the organic fraction of m/z 44 (f44), a surrogate for the degree of

213

oxidation and thus photochemical age of the air mass

214

NO3 (R2 = 0.74) due to its volatile nature, as observed in other datasets e. g. 11, 21.

11

40

. Usually in ambient datasets

. This is not the case for the present data (R2 = 0.30). We speculate that

43

. SV-OOA correlates well with

215

Regional overview. Figure 3 gives an overview of the normalized mean

216

concentrations of PM1 components measured on-road per zone and time interval. For a

217

comparison of normalized to measured mean values see Table S2. Corrected particle

218

number and gas phase concentrations are given in Table S4 in the SI. Apparent from

219

Figure 3 is the relatively uniform distribution of PM1 component concentrations across

220

the Barcelona metropolitan region. EBC is the major component in all zones, with the

221

highest zonal mean value measured in the heavily populated Vallés region (zone 5)

222

followed by the city center (zone 1). These EBC values are biased high due to the

ACS Paragon Plus Environment

10

Page 11 of 30

Environmental Science & Technology

223

influence of surrounding vehicles when measuring on-road. This influence is also visible

224

in the large differences of EBC mean zonal values between different times of day

225

(horizontal bars). In addition, urban EBC background concentrations might be more

226

reduced compared to on-road concentrations (due to tailpipe emissions being a major

227

source) than for AMS species. This effect is even more distinct during daytime when

228

traffic activity is high. Mobile measurements were only performed during the day. The

229

normalization to relatively lower background data will thus bias EBC values high relative

230

to AMS species and make the direct comparison difficult. Consequently, zonal mean

231

EBC values should only be interpreted qualitatively for concentration levels in different

232

zones, but at the same time are indicative of high emissions in the entire region.

233

The inorganic PM1 components NO3, NH4, SO4, and Chloride (Chl) show little

234

variability in concentration and fractional contribution across the Barcelona metropolitan

235

area. An exception are the industrial zones of the Llobregat valley, where Chl shows

236

significantly higher concentrations, most likely due to HCl emissions from industries

237

(not discernible in Figure 3, compare SI, Figure S8). The ammonium nitrate/ammonium

238

sulfate ratio is < 1 except in zone 3 and zone 6 (3.6 and 4.4 µg m-3 of NO3, respectively,

239

compared to ~2.5 µg m-3 in the other zones), which can be attributed to higher NOx

240

emissions (in the industrial area) and/or traffic activity while measuring. Shipping

241

emissions in front of the Barcelona coast are very high due to the high traffic density of

242

ships traveling from all across the Mediterranean towards the Gibraltar straight and

243

impacting coastal air quality 41. The measured SO4 concentrations can thus to some extent

244

be related to oxidation of SO2 or direct emissions of H2SO4 from diesel fuel combustion

245

in ship engines

45

44

. The diurnal differences of the inorganic PM1 component

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 30

246

concentrations are relatively small and indicate that they are less influenced by very local

247

emissions but more part of the regional pollution problem. Zone 4 shows slight decreases

248

of all components in the afternoon compared to morning hours. With the geographical

249

orientation of the Llobregat valley (perpendicular to the shoreline, not downwind of the

250

city center), the sea breeze might have a slight cleansing effect on this region.

251

OA exhibits the second highest concentrations of PM1 components in all zones.

252

Total OA concentrations are similar across all zones (zonal mean values between 4.9

253

(zone 4, t4) and 9.5 (zone 1, t1) µg m-3) and relatively high compared to concentrations

254

measured across Europe

255

PMF vary more considerably. Figure 4 shows the spatial variation of OA and the

256

corrected relative contributions of the PMF factor loadings to OA per zone and time

257

interval (for a comparison of normalized to measured zonal mean values of PMF factors

258

see Table S2). The highest organic concentrations are measured in the city center during

259

t1. The contribution to OA from OOA (sum of SV- and LV-OOA) is only 34%, HOA

260

dominates the primary components with 38% (3.3 µg m-3) followed by COA with 16%

261

(1.4 µg m-3). Similar HOA fractions were found for roadside measurements with a mobile

262

laboratory in New York City 46. In the other zones, HOA makes up between 14% (1.1 µg

263

m-3, zone 4, t5) and 29% (2.4 µg m-3, zone 4, t1). Note again that these values are highly

264

influenced by traffic conditions during on-road measurements. For all zones, BBOA

265

makes up between 10 and 12% of OA (range of absolute values: 0.6 – 1.1 µg m-3) and is

266

attributed to a mix of domestic solid fuel combustion

267

(observed especially in zone 3, where the highest absolute concentration was measured).

268

Relative COA contributions to OA are highest in zone 4 (14 – 18%, 0.7 – 1.3 µg m-3),

15

. The contributions of the different components as found by

31

ACS Paragon Plus Environment

and open agricultural fires

12

Page 13 of 30

Environmental Science & Technology

269

10% in the remote zone 6 (0.7 µg m-3) and between 10 and 15% (0.7 – 0.9 µg m-3) in

270

zones 2, 3, and 5. Even though relative contributions of COA to OA exhibit similar

271

values in the entire region, absolute values decrease from city center to zone 6 by a factor

272

of 2. There is a slight increase in relative HOA and COA contributions in the afternoon

273

hours, likely influenced by the diurnal wind pattern. In the afternoon, sea breezes advect

274

air masses from the city center 26 to the zones north of the city. These air masses contain

275

relatively more HOA and COA leading to a slight decrease in fractional SV-OOA

276

contribution, formed in-situ during the morning hours by photochemical reactions of

277

gaseous precursors. The relative contribution of the more regional LV-OOA stays

278

constant during this meteorological pattern. Exceptions are zone 4 (coastal range) which

279

is less impacted by local emissions and often above the city’s mixing layer 26, and zone 2

280

(harbor), where the fractional contribution of LV-OOA to OA increases in the afternoon

281

due to the sea breeze advecting aged air masses.

282

In all zones except zone 1, OOA mostly dominates (48 – 63%), as observed in

283

other locations (both urban and rural) in the northern hemisphere 2. Biogenic or, more

284

generally, non-fossil emissions might contribute substantially to OOA, especially

285

downwind of the city center 32. The highest OOA contribution (40% SV-OOA, 2.7 µg m-

286

3

287

We attribute the high SV-OOA level to newly formed OA during hours with high

288

photochemical activity from gaseous precursors emitted in the city center, which can be

289

upwind of the harbor at this time of day 28. The high LV-OOA contribution to OA (32%,

290

2.2 µg m-3) may be partly due to aging of SV-OOA and advection of processed air

291

masses by the sea breeze. In addition, the two afternoon drives in zone 2 were both

, and 23% LV-OOA, 1.6 µg m-3) was observed in zone 2 (harbor) in the late morning.

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 30

292

during meteorological episode A, characterized by the most severe pollution episodes and

293

high LV-OOA levels due to recirculation of air masses. The LV-OOA concentrations

294

might thus be biased high despite normalization to urban background values. Zone 6

295

(Montserrat), a forested, mountainous, and sparsely populated region, exhibits the

296

second-highest contribution of OOA. In contrast to zone 2, t1, LV-OOA (42%, 2.9 µg m-

297

3

298

general trend in increasing contributions from (LV-)OOA, indicative of formation and

299

aging of OA in the plume as it evolves in the atmosphere downwind of the city. This will

300

be investigated in the next section.

) is higher than SV-OOA (21%, 1.5 µg m-3). From the city center to zone 6 there is a

301

OOA formation in the urban plume. The formation of secondary OA from

302

anthropogenic precursors in the plume downwind of a city can be investigated by the

303

ratio of OA to CO

304

value) is used as an inert tracer of urban emissions including aerosol precursors to

305

account for dilution effects (evaporation upon dilution is assumed not to affect this

306

relationship 24, 49). We use the portion of the measurement drive on March 20 from zone 1

307

(city center) to zone 6 (Montserrat, pre-coastal mountain range) to estimate

308

anthropogenic secondary OOA formation downwind of the Barcelona city center. The

309

meteorological conditions were governed by a transition period with low mixing layer

310

height and a strong sea breeze effect

311

well after the sea breeze was established, at the shore upwind of the city center and

312

reached Montserrat, with 711 m. a. s. l. the highest point in the measurement route, at

313

16:05 local time (see portion of the route highlighted in red, Figure 1). Figure 5 shows the

314

measured OA plotted vs CO concentrations. Note that no CO background was subtracted

22, 47-49

. Generally, ΔCO (the CO concentrations minus a background

26

. The measurement started at ~14:00 local time,

ACS Paragon Plus Environment

14

Page 15 of 30

Environmental Science & Technology

315

here (issues with CO background subtractions were also reported by e. g. Kleinman et al.

316

47

317

well below the urban background mean of 110 ± 15 ppbv. In addition, the OA, but more

318

so the CO concentrations measured on-road are heavily impacted by the emissions from

319

traffic close by and not necessarily representative of plume mean values (see Figure S9a

320

in the SI). OA/CO ratios can also vary depending on the influence of primary sources

321

such as biomass burning. The numbers reported here should thus be regarded as a case

322

study only and are prone to high levels of uncertainties. Two subsets of data, marked by

323

the stars in Figure 1 and representative of a near-source (blue) and distant-from-source

324

(orange) OA/CO ratio, respectively, were linearly fit (Figure 5). The slope for the city

325

center data points (near-source) equals ~10 µg m-3 ppmv-1, comparable to values found

326

for e. g. Mexico City 47. Fitting the data points measured in zone 6 at Montserrat yields a

327

slope of ~85 µg m-3 ppmv-1, comparable to slopes found for aged (> 1 day) air masses in

328

other studies (pink slope in Figure 5 e. g. 47-49). Since the aging time of this case study is

329

much less than a day (~2.5 hours), the secondary OA mass formed downwind of the

330

Barcelona city center might not be purely of anthropogenic origin. In fact, Minguillón et

331

al.

332

coastal mountain range, ~72% of secondary OA carbon was non-fossil and from

333

biogenic, biomass burning, and urban non-fossil sources (e. g. cooking). Ignoring that

334

urban non-fossil sources also count as anthropogenic precursors, and taking only the

335

fossil 28% of secondary OA into account, yields a slope of OA/CO of ~ 24 µg m-3 ppmv-1

336

(black slope in Figure 5). Kleinman et al.

337

photochemical lifetimes – a rough comparison shows that a slope of 24 µg m-3 ppmv-1

) due to the measurement route crossing two mountain ranges with CO concentrations

32

found that at Montseny, a measurement station during DAURE on the same pre-

32

calculated OA/CO ratios for different

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 30

338

corresponds to ~0.1 unit of photochemical age which in turn corresponds to ~2.5 hours,

339

in good agreement with the time difference of 2 hours between the start and end point of

340

this measurement drive. The calculated secondary OA mass yield of 75 µg m-3 ppmv-1

341

(difference between ~85 µg m-3 ppmv-1 at Montseny and ~10 µg m-3 ppmv-1 in the city

342

center) was compared to the mass increase in SV- and LV-OOA (normalized to the fit

343

CO plume mean values scaled to 1, Figure S9) between start and end point of the

344

measurement drive and showed good agreement as well (Figure S9b, SI).

345

Overall, the spatial distribution of on-road measured PM1 components shows that

346

during March 2009 the entire Barcelona metropolitan area was heavily impacted by local

347

emissions from traffic, cooking activities and, to a lesser extent, biomass burning.

348

Industrial areas, the Barcelona airport, harbor, an extensive highway network and

349

agricultural areas surround the city center. During the measurements several open fires

350

were observed by the drivers (zone 3). The municipality of Barcelona has one of the

351

highest car densities per km2 in Europe 50, and (primary) road traffic emission effects also

352

include high concentrations of NO2 (48 ppbv), CO (556 ppbv), and particle number

353

concentrations (115000 cm-3). Note that since measurements were done on road, the

354

concentrations of compounds related to traffic emissions will be enhanced compared to

355

background concentrations, but may be indicative of near-road concentrations. The

356

formation of secondary aerosol, mostly OA, favored by high precursor emissions of both

357

anthropogenic and biogenic origin, and high photochemical activity in the Mediterranean

358

climate adds more than half of the OA burden outside the city center. We conclude that

359

the Barcelona metropolitan area is an important source of regional PM1, compared to e. g.

360

Paris, which has a very low impact on the air quality of its surroundings

ACS Paragon Plus Environment

51

. In addition,

16

Page 17 of 30

Environmental Science & Technology

361

the diurnal cycle of mountain wind and sea breeze can lead to recirculation of processed

362

air masses mixing with new emissions, possibly causing slight changes in the relative

363

contributions of the various components as a function of time of day. However, since this

364

affects both the city center and the surrounding areas, and since the whole region is

365

heavily populated and has abundant primary emissions sources, the overall concentrations

366

and composition of PM1 are relatively similar across the whole area.

367

368

Acknowledgement

369

We thank the organizers of the DAURE project and the staff at IDAEA for their support.

370

We acknowledge the IMBALANCE project of the Competence Center Environment and

371

Sustainability (CCES) and the EU-FP7 project EUCAARI for financial support and

372

“Accion Complementaria DAURE” from the Spanish Ministry of Science and Innovation

373

(CGL2007-30502-E/CLI) for infrastructure support. P. F. DeCarlo is grateful for the

374

postdoctoral support from the US-NSF (IRFP# 0701013).

375 376

Supporting Information Available

377

This information is available free of charge via the Internet at http://pubs.acs.org.

378

ACS Paragon Plus Environment

17

Environmental Science & Technology

379 380 381 382 383 384 385 386

Page 18 of 30

Table 1: Zones the Barcelona measurement route was divided into, their driving times, and dates, the prevailing atmospheric scenarios as defined by Pandolfi et al. 26, and the time intervals the data were binned to for the investigation of the influence of the time of day on the spatial variation. A single drive corresponds to a single time interval in a zone on a particular day, e. g. 10:30 – 12:45 in zone 1 on 19.03. The times for AM drives are printed in bold for clarity. *For a short description of the atmospheric scenarios see main text

387 388 389 390

ACS Paragon Plus Environment

18

Page 19 of 30

391 392 393 394 395 396 397 398 399

Environmental Science & Technology

Figure 1: Route driven for the mobile measurements, and zones defined for the analysis. The part highlighted in red shows the traverse driven on 20 March. The two stars with the same color bracket the subparts of the route used for the calculation of the OA/CO slopes in the two different regions, respectively. The red dot depicts the location of the urban background station during the DAURE campaign 21.

ACS Paragon Plus Environment

19

Environmental Science & Technology

400 401 402 403 404 405

Page 20 of 30

Figure 2: Factor profiles of the 5-factor high resolution PMF solution found for the organic data matrix measured by AMS. Peaks are colored according to their elemental composition.

ACS Paragon Plus Environment

20

Page 21 of 30

Environmental Science & Technology

406

407 408 409 410 411 412

Figure 3: Average normalized concentrations for PM1 components per zone. Horizontal lines give the average concentrations for the different time intervals. The error bars represent the standard error of the mean.

ACS Paragon Plus Environment

21

Environmental Science & Technology

Page 22 of 30

413

414 415 416 417 418 419 420

Figure 4: Spatial variation of the normalized OA (Org) concentration and average relative composition of OA per zone and time interval. The standard errors of the mean of the reported percentage values are below 1% except for HOA, where they are up to 3%. The higher uncertainty for HOA is due to the influence of surrounding traffic on-road.

ACS Paragon Plus Environment

22

Page 23 of 30

Environmental Science & Technology

421

422 423 424 425 426 427 428 429 430 431

Figure 5: OA (Org, µg m-3) vs CO (ppbv) for the portion of the drive highlighted in Figure 1 on March 20, colored by local time. The subsets of the data used to calculate the slopes using orthogonal distance regression are indicated by the stars in Figure 1 and color coded accordingly. The black curve gives the slope for purely fossil secondary OA based on the findings by Minguillón et al. 32. a and b denotes slopes and intercepts ( only where data was fit, see text).

ACS Paragon Plus Environment

23

Environmental Science & Technology

Page 24 of 30

432

References

433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474

1. Gurjar, B. R.; Nagpure, A. S.; Singh, T. P. H. H. T. E. A. q. i. m., Air quality in megacities. In Encyclopedia of Earth, Cutler J. Cleveland, Environmental Information Coalition, National Council for Science and the Environment: Washington, D.C., 2010. 2. Zhang, Q.; Jimenez, J. L.; Canagaratna, M. R.; Allan, J. D.; Coe, H.; Ulbrich, I.; Alfarra, M. R.; Takami, A.; Middlebrook, A. M.; Sun, Y. L.; Dzepina, K.; Dunlea, E.; Docherty, K.; DeCarlo, P. F.; Salcedo, D.; Onasch, T.; Jayne, J. T.; Miyoshi, T.; Shimono, A.; Hatakeyama, S.; Takegawa, N.; Kondo, Y.; Schneider, J.; Drewnick, F.; Borrmann, S.; Weimer, S.; Demerjian, K.; Williams, P.; Bower, K.; Bahreini, R.; Cottrell, L.; Griffin, R. J.; Rautiainen, J.; Sun, J. Y.; Zhang, Y. M.; Worsnop, D. R., Ubiquity and dominance of oxygenated species in organic aerosols in anthropogenically-influenced Northern Hemisphere midlatitudes. Geophys. Res. Lett. 2007, 34, (13), L13801. 3. Goldstein, A. H.; Galbally, I. E., Known and unexplored organic constituents in the earth's atmosphere. Environ. Sci. Technol. 2007, 41, (5), 1514-1521. 4. Kalberer, M.; Paulsen, D.; Sax, M.; Steinbacher, M.; Dommen, J.; Prevot, A. S. H.; Fisseha, R.; Weingartner, E.; Frankevich, V.; Zenobi, R.; Baltensperger, U., Identification of Polymers as Major Components of Atmospheric Organic Aerosols. Science 2004, 303, (5664), 1659-1662. 5. Shiraiwa, M.; Yee, L. D.; Schilling, K. A.; Loza, C. L.; Craven, J. S.; Zuend, A.; Ziemann, P. J.; Seinfeld, J. H., Size distribution dynamics reveal particle-phase chemistry in organic aerosol formation. Proceedings of the National Academy of Sciences 2013, 110, (29), 11746–11750. 6. Hallquist, M.; Wenger, J. C.; Baltensperger, U.; Rudich, Y.; Simpson, D.; Claeys, M.; Dommen, J.; Donahue, N. M.; George, C.; Goldstein, A. H.; Hamilton, J. F.; Herrmann, H.; Hoffmann, T.; Iinuma, Y.; Jang, M.; Jenkin, M. E.; Jimenez, J. L.; Kiendler-Scharr, A.; Maenhaut, W.; McFiggans, G.; Mentel, T. F.; Monod, A.; Prevot, A. S. H.; Seinfeld, J. H.; Surratt, J. D.; Szmigielski, R.; Wildt, J., The formation, properties and impact of secondary organic aerosol: current and emerging issues. Atmos. Chem. Phys. 2009, 9, (14), 5155-5236. 7. Donahue, N. M.; Kroll, J. H.; Pandis, S. N.; Robinson, A. L., A two-dimensional volatility basis set - Part 2: Diagnostics of organic-aerosol evolution. Atmos. Chem. Phys. Discuss. 2011, 11, (9), 24883-24931. 8. Jimenez, J. L.; Canagaratna, M. R.; Donahue, N. M.; Prevot, A. S. H.; Zhang, Q.; Kroll, J. H.; DeCarlo, P. F.; Allan, J. D.; Coe, H.; Ng, N. L.; Aiken, A. C.; Docherty, K. S.; Ulbrich, I. M.; Grieshop, A. P.; Robinson, A. L.; Duplissy, J.; Smith, J. D.; Wilson, K. R.; Lanz, V. A.; Hueglin, C.; Sun, Y. L.; Tian, J.; Laaksonen, A.; Raatikainen, T.; Rautiainen, J.; Vaattovaara, P.; Ehn, M.; Kulmala, M.; Tomlinson, J. M.; Collins, D. R.; Cubison, M. J.; Dunlea, E. J.; Huffman, J. A.; Onasch, T. B.; Alfarra, M. R.; Williams, P. I.; Bower, K.; Kondo, Y.; Schneider, J.; Drewnick, F.; Borrmann, S.; Weimer, S.; Demerjian, K.; Salcedo, D.; Cottrell, L.; Griffin, R.; Takami, A.; Miyoshi, T.; Hatakeyama, S.; Shimono, A.; Sun, J. Y.; Zhang, Y. M.; Dzepina, K.; Kimmel, J. R.; Sueper, D.; Jayne, J. T.; Herndon, S. C.; Trimborn, A. M.; Williams, L. R.; Wood, E. C.;

ACS Paragon Plus Environment

24

Page 25 of 30

475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520

Environmental Science & Technology

Middlebrook, A. M.; Kolb, C. E.; Baltensperger, U.; Worsnop, D. R., Evolution of organic aerosols in the atmosphere. Science 2009, 326, (5959), 1525-1529. 9. Canagaratna, M. R.; Jayne, J. T.; Jimenez, J. L.; Allan, J. D.; Alfarra, M. R.; Zhang, Q.; Onasch, T. B.; Drewnick, F.; Coe, H.; Middlebrook, A.; Delia, A.; Williams, L. R.; Trimborn, A. M.; Northway, M. J.; DeCarlo, P. F.; Kolb, C. E.; Davidovits, P.; Worsnop, D. R., Chemical and microphysical characterization of ambient aerosols with the Aerodyne aerosol mass spectrometer. Mass Spectrom. Rev. 2007, 26, (2), 185-222. 10. Paatero, P.; Tapper, U., Positive matrix factorization - a nonnegative factor model with optimal utilization of error-estimates of data values. Environmetrics 1994, 5, (2), 111-126. 11. Lanz, V. A.; Alfarra, M. R.; Baltensperger, U.; Buchmann, B.; Hueglin, C.; Prevot, A. S. H., Source apportionment of submicron organic aerosols at an urban site by factor analytical modelling of aerosol mass spectra. Atmos. Chem. Phys. 2007, 7, (6), 1503-1522. 12. Ulbrich, I. M.; Canagaratna, M. R.; Zhang, Q.; Worsnop, D. R.; Jimenez, J. L., Interpretation of organic components from positive matrix factorization of aerosol mass spectrometric data. Atmos. Chem. Phys. 2009, 9, (9), 2891-2918. 13. Zhang, Q.; Worsnop, D. R.; Canagaratna, M. R.; Jimenez, J. L., Hydrocarbon-like and oxygenated organic aerosols in Pittsburgh: insights into sources and processes of organic aerosols. Atmos. Chem. Phys. 2005, 5, 3289-3311. 14. Allan, J. D.; Williams, P. I.; Morgan, W. T.; Martin, C. L.; Flynn, M. J.; Lee, J.; Nemitz, E.; Phillips, G. J.; Gallagher, M. W.; Coe, H., Contributions from transport, solid fuel burning and cooking to primary organic aerosols in two UK cities. Atmos. Chem. Phys. 2010, 10, (2), 647-668. 15. Crippa, M.; Canonaco, F.; Lanz, V. A.; Äijälä, M.; Allan, J. D.; Carbone, S.; Capes, G.; Ceburnis, D.; Dall'Osto, M.; Day, D. A.; DeCarlo, P. F.; Ehn, M.; Eriksson, A.; Freney, E.; Hildebrandt Ruiz, L.; Hillamo, R.; Jimenez, J. L.; Junninen, H.; KiendlerScharr, A.; Kortelainen, A. M.; Kulmala, M.; Laaksonen, A.; Mensah, A. A.; Mohr, C.; Nemitz, E.; O'Dowd, C.; Ovadnevaite, J.; Pandis, S. N.; Petäjä, T.; Poulain, L.; Saarikoski, S.; Sellegri, K.; Swietlicki, E.; Tiitta, P.; Worsnop, D. R.; Baltensperger, U.; Prévôt, A. S. H., Organic aerosol components derived from 25 AMS datasets across Europe using a consistent ME-2 based source apportionment strategy. Atmos. Chem. Phys. 2014, 14, 6159-6176. 16. Lanz, V. A.; Prévôt, A. S. H.; Alfarra, M. R.; Weimer, S.; Mohr, C.; DeCarlo, P. F.; Gianini, M. F. D.; Hueglin, C.; Schneider, J.; Favez, O.; D'Anna, B.; George, C.; Baltensperger, U., Characterization of aerosol chemical composition with aerosol mass spectrometry in Central Europe: an overview. Atmos. Chem. Phys. 2010, 10, (21), 1045310471. 17. Bukowiecki, N.; Dommen, J.; Prevot, A. S. H.; Weingartner, E.; Baltensperger, U., Fine and ultrafine particles in the Zurich (Switzerland) area measured with a mobile laboratory: an assessment of the seasonal and regional variation throughout a year. Atmos. Chem. Phys. 2003, 3, 1477-1494. 18. Weimer, S.; Mohr, C.; Richter, R.; Keller, J.; Mohr, M.; Prévôt, A. S. H.; Baltensperger, U., Mobile measurements of aerosol number and volume size distributions in an Alpine valley: Influence of traffic versus wood burning. Atmos. Environ. 2009, 43, (3), 624-630.

ACS Paragon Plus Environment

25

Environmental Science & Technology

521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565

Page 26 of 30

19. Bahreini, R.; Jimenez, J. L.; Wang, J.; Flagan, R. C.; Seinfeld, J. H.; Jayne, J. T.; Worsnop, D. R., Aircraft-based aerosol size and composition measurements during ACEAsia using an Aerodyne aerosol mass spectrometer. J. Geophys. Res-Atmos. 2003, 108, (D23), D238645. 20. Schneider, J.; Hings, S. S.; Hock, B. N.; Weimer, S.; Borrmann, S.; Fiebig, M.; Petzold, A.; Busen, R.; Karcher, B., Aircraft-based operation of an aerosol mass spectrometer: Measurements of tropospheric aerosol composition. J. Aerosol. Sci. 2006, 37, (7), 839-857. 21. Mohr, C.; DeCarlo, P. F.; Heringa, M. F.; Chirico, R.; Slowik, J. G.; Richter, R.; Reche, C.; Alastuey, A.; Querol, X.; Seco, R.; Peñuelas, J.; Jimenez, J. L.; Crippa, M.; Zimmermann, R.; Baltensperger, U.; Prevot, A. S. H., Identification and quantification of organic aerosol from cooking and other sources in Barcelona using aerosol mass spectrometer data. Atmos. Chem. Phys. 2012, 12, (4), 1649 - 1665. 22. DeCarlo, P. F.; Ulbrich, I. M.; Crounse, J.; de Foy, B.; Dunlea, E. J.; Aiken, A. C.; Knapp, D.; Weinheimer, A. J.; Campos, T.; Wennberg, P. O.; Jimenez, J. L., Investigation of the sources and processing of organic aerosol over the Central Mexican Plateau from aircraft measurements during MILAGRO. Atmos. Chem. Phys. 2010, 10, (12), 5257-5280. 23. von der Weiden-Reinmüller, S. L.; Drewnick, F.; Crippa, M.; Prévôt, A. S. H.; Meleux, F.; Baltensperger, U.; Beekmann, M.; Borrmann, S., Application of mobile aerosol and trace gas measurements for the investigation of megacity air pollution emissions: the Paris metropolitan area. Atmos. Meas. Tech. 2014, 7, (1), 279-299. 24. DeCarlo, P. F.; Dunlea, E. J.; Kimmel, J. R.; Aiken, A. C.; Sueper, D.; Crounse, J.; Wennberg, P. O.; Emmons, L.; Shinozuka, Y.; Clarke, A.; Zhou, J.; Tomlinson, J.; Collins, D. R.; Knapp, D.; Weinheimer, A. J.; Montzka, D. D.; Campos, T.; Jimenez, J. L., Fast airborne aerosol size and chemistry measurements above Mexico City and Central Mexico during the MILAGRO campaign. Atmos. Chem. Phys. 2008, 8, (14), 4027-4048. 25. Canagaratna, M. R.; Onasch, T. B.; Wood, E. C.; Herndon, S. C.; Jayne, J. T.; Cross, E. S.; Miake-Lye, R. C.; Kolb, C. E.; Worsnop, D. R., Evolution of Vehicle Exhaust Particles in the Atmosphere. J. Air Waste Manage. 2010, 60, (10), 1192-1203. 26. Pandolfi, M.; Querol, X.; Alastuey, A.; Jimenez, J. L.; Jorba, O.; Day, D.; Ortega, A.; Cubison, M. J.; Comerón, A.; Sicard, M.; Mohr, C.; Prévôt, A. S. H.; Minguillón, M. C.; Pey, J.; Baldasano, J. M.; Burkhart, J. F.; Seco, R.; Peñuelas, J.; van Drooge, B. L.; Artiñano, B.; Di Marco, C.; Nemitz, E.; Schallhart, S.; Metzger, A.; Hansel, A.; Lorente, J.; Ng, S.; Jayne, J.; Szidat, S., Effects of Sources and Meteorology on Particulate Matter in the Western Mediterranean Basin: An overview of the DAURE campaign J. Geophys. Res-Atmos. 2014, 119, (8), 4978-5010. 27. Pérez, N.; Pey, J.; Castillo, S.; Viana, M.; Alastuey, A.; Querol, X., Interpretation of the variability of levels of regional background aerosols in the Western Mediterranean. Sci. Total Environ. 2008, 407, (1), 527-540. 28. Jorba, O.; Pandolfi, M.; Spada, M.; Baldasano, J. M.; Pey, J.; Alastuey, A.; Arnold, D.; Sicard, M.; Artiñano, B.; Revuelta, M. A.; Querol, X., The DAURE field campaign: meteorological overview. Atmos. Chem. Phys. Discuss. 2011, 11, (2), 49535001.

ACS Paragon Plus Environment

26

Page 27 of 30

566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610

Environmental Science & Technology

29. Pérez, N.; Pey, J.; Querol, X.; Alastuey, A.; López, J. M.; Viana, M., Partitioning of major and trace components in PM10-PM2.5-PM1 at an urban site in Southern Europe. Atmos. Environ. 2008, 42, (8), 1677-1691. 30. Reche, C.; Viana, M.; Moreno, T.; Querol, X.; Alastuey, A.; Pey, J.; Pandolfi, M.; Prévôt, A.; Mohr, C.; Richard, A.; Artiñano, B.; Gomez-Moreno, F. J.; Cots, N., Peculiarities in atmospheric particle number and size-resolved speciation in an urban area in the western Mediterranean: Results from the DAURE campaign. Atmos. Environ. 2011, 45, (30), 5282-5293. 31. Reche, C.; Viana, M.; Amato, F.; Alastuey, A.; Moreno, T.; Hillamo, R.; Teinilä, K.; Saarnio, K.; Seco, R.; Peñuelas, J.; Mohr, C.; Prévôt, A. S. H.; Querol, X., Biomass burning contributions to urban aerosols in a coastal Mediterranean City. Sci. Total Environ. 2012, 427–428, (0), 175-190. 32. Minguillón, M. C.; Perron, N.; Querol, X.; Szidat, S.; Fahrni, S.; Alastuey, A.; Jimenez, J. L.; Mohr, C.; Ortega, A.; Day, D. A.; Lanz, V. A.; Wacker, L.; Reche, C.; Cusack, M.; Amato, F.; Kiss, G.; Hoffer, A.; Decesari, S.; Moretti, F.; Hillamo, R.; Teinilä, K.; Seco, R.; Peñuelas, J.; Metzger, A.; Schallhart, S.; Müller, M.; Hansel, A.; Burkhart, J.; Baltensperger, U.; Prevot, A. S. H., Fossil versus contemporary sources of fine elemental and organic carbonaceous particulate matter during the DAURE campaign in Northeast Spain. Atmos. Chem. Phys. 2011, 11, (23), 12067 - 12084. 33. Moreno, T.; Querol, X.; Alastuey, A.; Reche, C.; Cusack, M.; Amato, F.; Pandolfi, M.; Pey, J.; Richard, A.; Prévôt, A. S. H.; Furger, M.; Gibbons, W., Variations in time and space of trace metal aerosol concentrations in urban areas and their surroundings. Atmos. Chem. Phys. 2011, 11, (17), 9415-9430. 34. Mohr, C.; Richter, R.; DeCarlo, P. F.; Prévôt, A. S. H.; Baltensperger, U., Spatial variation of chemical composition and sources of submicron aerosol in Zurich during wintertime using mobile aerosol mass spectrometer data. Atmos. Chem. Phys. 2011, 11, (15), 7465-7482. 35. DeCarlo, P. F.; Kimmel, J. R.; Trimborn, A.; Northway, M. J.; Jayne, J. T.; Aiken, A. C.; Gonin, M.; Fuhrer, K.; Horvath, T.; Docherty, K. S.; Worsnop, D. R.; Jimenez, J. L., Field-deployable, high-resolution, time-of-flight aerosol mass spectrometer. Anal. Chem. 2006, 78, (24), 8281-8289. 36. Middlebrook, A. M.; Bahreini, R.; Jimenez, J. L.; Canagaratna, M. R., Evaluation of composition-dependent collection efficiencies for the Aerodyne aerosol mass spectrometer using field data. Aerosol. Sci. Technol. 2012, 46, (3), 258–271. 37. Wiedensohler, A.; Birmili, W.; Nowak, A.; Sonntag, A.; Weinhold, K.; Merkel, M.; Wehner, B.; Tuch, T.; Pfeifer, S.; Fiebig, M.; Fjäraa, A. M.; Asmi, E.; Sellegri, K.; Depuy, R.; Venzac, H.; Villani, P.; Laj, P.; Aalto, P.; Ogren, J. A.; Swietlicki, E.; Williams, P.; Roldin, P.; Quincey, P.; Hüglin, C.; Fierz-Schmidhauser, R.; Gysel, M.; Weingartner, E.; Riccobono, F.; Santos, S.; Grüning, C.; Faloon, K.; Beddows, D.; Harrison, R.; Monahan, C.; Jennings, S. G.; O'Dowd, C. D.; Marinoni, A.; Horn, H. G.; Keck, L.; Jiang, J.; Scheckman, J.; McMurry, P. H.; Deng, Z.; Zhao, C. S.; Moerman, M.; Henzing, B.; de Leeuw, G.; Löschau, G.; Bastian, S., Mobility particle size spectrometers: harmonization of technical standards and data structure to facilitate high quality long-term observations of atmospheric particle number size distributions. Atmos. Meas. Tech. 2012, 5, (3), 657-685.

ACS Paragon Plus Environment

27

Environmental Science & Technology

611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656

Page 28 of 30

38. Paatero, P.; Hopke, P. K.; Song, X. H.; Ramadan, Z., Understanding and controlling rotations in factor analytic models. Chemometr. Intell. Lab. 2002, 60, (1-2), 253-264. 39. Weimer, S.; Alfarra, M. R.; Schreiber, D.; Mohr, M.; Prévôt, A. S. H.; Baltensperger, U., Organic aerosol mass spectral signatures from wood-burning emissions: Influence of burning conditions and wood type. J. Geophys. Res. 2008, 113, (D10), D10304. 40. Massoli, P.; Fortner, E. C.; Canagaratna, M. R.; Williams, L. R.; Zhang, Q.; Sun, Y.; Schwab, J. J.; Trimborn, A.; Onasch, T. B.; Demerjian, K. L.; Kolb, C. E.; Worsnop, D. R.; Jayne, J. T., Pollution Gradients and Chemical Characterization of Particulate Matter from Vehicular Traffic near Major Roadways: Results from the 2009 Queens College Air Quality Study in NYC. Aerosol. Sci. Technol. 2012, 46, (11), 1201-1218. 41. Viana, M.; Hammingh, P.; Colette, A.; Querol, X.; Degraeuwe, B.; Vlieger, I. d.; van Aardenne, J., Impact of maritime transport emissions on coastal air quality in Europe. Atmos. Environ. 2014, 90, (0), 96-105. 42. Minguillón, M. C.; Ripoll, A.; Pérez, N.; Prévôt, A. S. H.; Canonaco, F.; Querol, X.; Alastuey, A., Chemical characterization of submicron regional background aerosols in the Western Mediterranean using an Aerosol Chemical Speciation Monitor. Atmos. Chem. Phys. Discuss. 2015, 15, (1), 965-1000. 43. Ng, N. L.; Canagaratna, M. R.; Zhang, Q.; Jimenez, J. L.; Tian, J.; Ulbrich, I. M.; Kroll, J. H.; Docherty, K. S.; Chhabra, P. S.; Bahreini, R.; Murphy, S. M.; Seinfeld, J. H.; Hildebrandt, L.; Donahue, N. M.; DeCarlo, P. F.; Lanz, V. A.; Prevot, A. S. H.; Dinar, E.; Rudich, Y.; Worsnop, D. R., Organic aerosol components observed in northern hemispheric datasets from aerosol mass spectrometry. Atmos. Chem. Phys. 2010, 10, (10), 4625-4641. 44. Querol, X.; Alastuey, A.; Rodriguez, S.; Plana, F.; Ruiz, C. R.; Cots, N.; Massagué, G.; Puig, O., PM10 and PM2.5 source apportionment in the Barcelona Metropolitan area, Catalonia, Spain. Atmos. Environ. 2001, 35, (36), 6407-6419. 45. Pandolfi, M.; Gonzalez-Castanedo, Y.; Alastuey, A.; de la Rosa, J.; Mantilla, E.; de la Campa, A.; Querol, X.; Pey, J.; Amato, F.; Moreno, T., Source apportionment of PM10 and PM2.5 at multiple sites in the strait of Gibraltar by PMF: impact of shipping emissions. Environ. Sci. Pollut. R. 2011, 18, (2), 260-269. 46. Sun, Y. L.; Zhang, Q.; Schwab, J. J.; Chen, W. N.; Bae, M. S.; Hung, H. M.; Lin, Y. C.; Ng, N. L.; Jayne, J.; Massoli, P.; Williams, L. R.; Demerjian, K. L., Characterization of near-highway submicron aerosols in New York City with a highresolution aerosol mass spectrometer. Atmos. Chem. Phys. 2012, 12, (4), 2215-2227. 47. Kleinman, L. I.; Springston, S. R.; Daum, P. H.; Lee, Y. N.; Nunnermacker, L. J.; Senum, G. I.; Wang, J.; Weinstein-Lloyd, J.; Alexander, M. L.; Hubbe, J.; Ortega, J.; Canagaratna, M. R.; Jayne, J., The time evolution of aerosol composition over the Mexico City plateau. Atmos. Chem. Phys. 2008, 8, (6), 1559-1575. 48. Aiken, A. C.; Salcedo, D.; Cubison, M. J.; Huffman, J. A.; DeCarlo, P. F.; Ulbrich, I. M.; Docherty, K. S.; Sueper, D.; Kimmel, J. R.; Worsnop, D. R.; Trimborn, A.; Northway, M.; Stone, E. A.; Schauer, J. J.; Volkamer, R. M.; Fortner, E.; de Foy, B.; Wang, J.; Laskin, A.; Shutthanandan, V.; Zheng, J.; Zhang, R.; Gaffney, J.; Marley, N. A.; Paredes-Miranda, G.; Arnott, W. P.; Molina, L. T.; Sosa, G.; Jimenez, J. L., Mexico

ACS Paragon Plus Environment

28

Page 29 of 30

657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674

Environmental Science & Technology

City aerosol analysis during MILAGRO using high resolution aerosol mass spectrometry at the urban supersite (T0) - Part 1: Fine particle composition and organic source apportionment. Atmos. Chem. Phys. 2009, 9, (17), 6633-6653. 49. Freney, E. J.; Sellegri, K.; Canonaco, F.; Colomb, A.; Borbon, A.; Michoud, V.; Doussin, J. F.; Crumeyrolle, S.; Amarouche, N.; Pichon, J. M.; Bourianne, T.; Gomes, L.; Prevot, A. S. H.; Beekmann, M.; Schwarzenböeck, A., Characterizing the impact of urban emissions on regional aerosol particles: airborne measurements during the MEGAPOLI experiment. Atmos. Chem. Phys. 2014, 14, (3), 1397-1412. 50. Ajuntament de Barcelona Dades bàsiques 2006; Direcció de serveis de mobilitat: Barcelona, 2007. 51. Crippa, M.; DeCarlo, P. F.; Slowik, J. G.; Mohr, C.; Heringa, M. F.; Chirico, R.; Poulain, L.; Freutel, F.; Sciare, J.; Cozic, J.; Di Marco, C. F.; Elsasser, M.; Nicolas, J. B.; Marchand, N.; Abidi, E.; Wiedensohler, A.; Drewnick, F.; Schneider, J.; Borrmann, S.; Nemitz, E.; Zimmermann, R.; Jaffrezo, J. L.; Prévôt, A. S. H.; Baltensperger, U., Wintertime aerosol chemical composition and source apportionment of the organic fraction in the metropolitan area of Paris. Atmos. Chem. Phys. 2013, 13, (2), 961-981.

ACS Paragon Plus Environment

29

Environmental Science & Technology

200x143mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 30 of 30