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