Spatial Variation of Aerosol Chemical Composition and Organic

Aug 3, 2015 - The spatial distribution of PM1 components in the Barcelona metropolitan area was investigated using on-road mobile measurements of ...
0 downloads 3 Views 8MB Size
Article pubs.acs.org/est

Spatial Variation of Aerosol Chemical Composition and Organic Components Identified by Positive Matrix Factorization in the Barcelona Region Claudia Mohr,*,†,§ Peter F. DeCarlo,†,∥ Maarten F. Heringa,† Roberto Chirico,†,⊥ René Richter,† Monica Crippa,†,# Xavier Querol,‡ Urs Baltensperger,† and André S. H. Prévôt† †

Laboratory of Atmospheric Chemistry, Paul Scherrer Institute (PSI), Villigen 5232, Switzerland Institute for Environmental Assessment and Water Research (IDAEA-CSIC), Barcelona 08034, Spain



S Supporting Information *

ABSTRACT: The spatial distribution of PM1 components in the Barcelona metropolitan area was investigated using on-road mobile measurements of atmospheric particle- and gas-phase compounds during the DAURE campaign in March 2009. Positive matrix factorization (PMF) applied to organic aerosol (OA) data yielded 5 factors: hydrocarbon-like OA (HOA), cooking OA (COA), biomass burning OA (BBOA), and low volatility and semivolatile oxygenated OA (LV-OOA and SV-OOA). The area under investigation (∼500 km2) was divided into six zones (city center, harbor, industrial area, precoastal depression, 2 mountain ranges) for measurements and data analysis. Mean zonal OA concentrations are 4.9−9.5 μg m−3. The area is heavily impacted by local primary emissions (HOA 14−38%, COA 10−18%, BBOA 10−12% of OA); concentrations of traffic-related components, especially black carbon, are biased high due to the on-road nature of the measurements. The formation of secondary OA adds more than half of the OA burden outside the city center (SV-OOA 14−40%, LV-OOA 17−42% of OA). A case study of one measurement drive from the shore to the precoastal mountain range furthest downwind of the city center indicates the importance of nonfossil over anthropogenic secondary OA based on OA/CO.



INTRODUCTION Adverse health effects of atmospheric aerosols are of special concern in highly populated areas.1 The identification and quantification of their sources, the basis for the design of mitigation strategies, remains challenging due to their complex nature. The organic fraction of atmospheric aerosols (organic aerosol, OA), often the major component in submicron particulate matter (PM1),2 is composed of 1000s of compounds3 subject to chemical reactions within the particles4,5 and the surrounding gas phase via condensation, evaporation, gas-phase oxidation, and recondensation processes.6−8 Recent advances in the apportionment of OA sources have been achieved through the combination of aerosol mass spectrometer (AMS) measurements9 and the factor analytical model positive matrix factorization (PMF).10 The application of this method has led to the identification of diverse OA components: Hydrocarbonlike OA (HOA, related to fresh traffic emissions), biomass burning OA (BBOA, from domestic wood burning or wildfire emissions), cooking OA, and different types of oxygenated OA (OOA, related to secondary OA), i.e., low volatility OOA (LVOOA), characterized as a regional, heavily aged OA, and the less aged semivolatile OOA (SV-OOA)11,.11,−16 Most of these studies were based on stationary measurements during a certain period of time. However, the evolution of pollutants in the plume © 2015 American Chemical Society

of a major source region such as an urban center or the characterization of the sources and processes of atmospheric pollutants of a whole region are best assessed by mobile measurements, ground based or using aircraft.17−25 The densely populated metropolitan area of Barcelona, Spain (5.4 million inhabitants) is subject to very high air pollution levels due to the combination of its Mediterranean climate, abundant biogenic and anthropogenic emissions, orography, and atmospheric dynamics. 26,27 Limited to the SE by the Mediterranean Sea and to the NW-SW by the Catalan coastal ranges, Barcelona experiences a daily cycle of nocturnal offshore flows, transporting air masses from the surrounding valleys (e.g., the densely populated Vallés) to the city center, and diurnal sea breeze, advecting polluted coastal air masses inland.28 The Llobregat river valley SW of the city center contains both extensive industrial and agricultural areas. This paper investigates the sources, processes, and evolution of PM1 in the Barcelona metropolitan area during late winter (when frequent cyclonic episodes further enhance pre-existing high Received: Revised: Accepted: Published: 10421

April 29, 2015 July 30, 2015 August 3, 2015 August 3, 2015 DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

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



PM1 concentrations29), with a special focus on the organic fraction. Mobile measurements were performed in March 2009 within the international DAURE campaign (Determination of the sources of atmospheric Aerosols in Urban and Rural Environments in the western Mediterranean).26 Concurrent stationary measurements of various particle- and gas-phase parameters at an urban background site investigated and quantified PM sources such as road traffic, constructiondemolition works, shipping emissions, and photochemical processes.30 Contributions of biomass burning31 or cooking21 emissions to OA, the fossil and nonfossil fractions of organic carbon,32 and the particulate trace metal content33 were assessed, among others. The results from the mobile measurements complement these findings by adding information on the evolution of the Barcelona city plume and the emission sources in the regions surrounding the city. Since measurements were carried out on-road, levels and composition of PM1 are influenced by close traffic emissions from circulating vehicles.

EXPERIMENTAL SECTION Mobile Measurements. Sixteen days of mobile measurements (61 sampling hours) were performed in the Barcelona metropolitan area with the Paul Scherrer Institute mobile laboratory, an IVECO Daily van equipped with several continuous online gas- and particle-phase instruments,34 between March 1 and 26, 2009. Figure 1 shows a map of the Barcelona metropolitan area and the driving route. It was divided into 6 different zones for both the measurements and the data analysis: City center (zone 1), harbor (zone 2), the industrialized Llobregat valley (zone 3), the coastal mountain range (zone 4), the precoastal mountain range (zone 6) and the heavily populated precoastal depression Vallés (zone 5). Table 1 summarizes the zones, driving dates, and times. For a discussion of zone identification, see SI. Instrumentation. Table S1 lists the instruments deployed in the mobile laboratory and gives information on the parameters measured, the detection limits, and the time resolution of the data acquisition. Details are provided elsewhere;34 only 10422

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology Table 1. Zones of Barcelona Measurement Route driving times zone

name

date (2009)

drive 1

1

City center

01.03 04.03 06.03 10.03 25.03 26.03 01.03 18.03 19.03 11.03 18.03 23.03 04.03 09.03 10.03 12.03 20.03 23.03 25.03 26.03. 12.03 20.03 23.03 26.03 12.03 20.03 23.03 26.03

08:45−11:15 10:15−12:30 09:45−12:00

2

3

4

Harbor

Industrial valley Llobregat Coastal mountain range

5

Precostal depression Vallés

6

Precoastal mountain range Montserrat

drive 2

15:45−20:15 15:00−18:30 07:15−11:45 11:15−13:00 10:30−12:45 08:45−12:00 09:30−12:15 11:45−12:30 11:15−11:45 17:15−17:45 14:45−15:00 14:30−15:00 14:00−14:15 16:15−16:30 08:30−09:00 15:00−15:30 15:00−15:30 14:15−14:45 15:00−15:30 15:30−16:30 15:30−16:30 14:45−15:45 15:30−16:15

13:00−16:00 14:45−17:00 14:15−17:00

19:45−20:00 17:15−17:45 17:15−18:00 16:30−16:45 18:00−18:15 11:15−11:30 16:30−17:00 16:30−17:00 15:45−16:15 16:15−16:45

atmospheric scenarioa

time intervals for data analysisb

T A C A T T T A A A A T A A A A T T T T A T T T A T T T

t1: 07:00−12:45 t2: 14:45−20:30

t1: 10:15−13:15 t2: 12:45−17:15 t1: 08:30−12:30 t2: 14:00−17:15 t1: 08:15−09:15 t2:11:00−12:45 t3: 13:45−15:15 t4: 16:00−18:30 t5:19:30−20:15

t1: 14:00−15:45 t2: 15:30−17:15

t1: 14:45−16:30

a The prevailing atmospheric scenarios as defined by Pandolfi et al.26 For a short description of the atmospheric scenarios see main text. bTime 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.

(1107, Grimm Labortechnik GmbH & Co. KG) and corrected with PM1 sample loadings collected on filters.21 This resulted in a CE = 0.7 for f NO3 ≤ 0.25, CE = 1 for f NO3 ≥ 0.78, and a linear increase between those values for 0.25 ≤ f NO3 ≤ 0.78 (Figure S1). The CO2+ ion signal was corrected for elevated gas-phase contributions by a time-dependent factor based on gaseous CO2 measurements. PMF was applied to organic high resolution (HR) mass spectral data (time resolution 6 s). The input matrices, comprising the complete mobile data set, were prepared as previously described.34 All data are in local standard time (UTC + 1) and at local pressure and temperature conditions (10−25 °C and 943−1048 hPa). Mobile measurements are usually carried out during short periods of time only (few hours per day during few weeks per year). Analyses of the spatial distribution of parameters based on mobile measurements thus always bear the risk of being overly influenced by meteorological conditions and diurnal patterns.23 This is especially true for the present study, where the mean concentration values for the different locations are based on aggregate data from several measurement drives spread out across several weeks. To reduce the influence of the temporal component on the spatial distribution, the time series of each parameter X measured in the mobile laboratory during a single drive (see Table 1) was divided by a normalization factor N based on the concentration of X measured simultaneously at the urban

deviations from this setup are provided here. The focus of the present analysis is on size-resolved chemical composition data of nonrefractory PM1 measured by the Aerodyne high-resolution time-of-flight aerosol mass spectrometer35 (AMS). AMS data were acquired in the higher sensitivity V-Mode35 at a time resolution of 6 s. The PM1 equivalent black carbon fraction (EBC) was derived from aerosol light absorption measurements at 670 nm by a multiangle absorption photometer (MAAP 5012, Thermo), time resolution 1 s. An ozone monitor (Dual Beam, Model 205, 2B Technologies, Inc.) added for this campaign measured O3 concentrations based on the absorption of UV light at 254 nm. NO2 concentrations are available from March 10, 2009 onward (when the chromium trioxide converter used to measure NOx34 ceased to function; earlier data, where a separation of NO and NO2 is not possible, are excluded from the analysis). All gases were measured with a time resolution of 1 s. Data Analysis. AMS data were processed using the software SQUIRREL v.1.51B and PIKA v1.1.10B within Igor Pro 6.22A (Wavemetrics). The collection efficiency (CE) was determined based on (a) the nitrate fraction f NO3 of the species measured by AMS36 and (b) comparison of the AMS when the mobile laboratory was measuring at the urban background site (location shown in Figure 1) with a scanning mobility particle sizer (SMPS, TSI 3071 type differential mobility analyzer37 custom built by PSI, condensation particle counter CPC 1022, TSI), a second AMS, and PM1 measured by an optical particle counter 10423

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology

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.

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.

background station. The normalization factor NDy, of species X for an individual drive, Dy, is given as follows: ⎧ X _BGD y ⎫ ⎬ ⎨ND y = X _BG ⎭ ⎩ ⎪







with X _BGDy denoting the background station mean value of X for the time interval of Dy and X _BG the background station campaign mean value of X. Individual factors for particle components instead of a single factor for all components based on total PM were used since different meteorological conditions

(1) 10424

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology

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.

as the factors’ atomic oxygen to carbon (O:C) and organic matter (OM) to organic carbon (OC) ratios for both solutions. The noisy characteristics of the 6-s mobile data and different magnitudes of the peaks resulted in high scatter in the correlations, thus 5 min averages were used for most of the time series correlations (numbers in italics). PMF was run in the robust mode and the solution was analyzed as outlined in the SI from Mohr et al.,21 and the 5-factor solution was chosen. The most central rotation (f peak = 0) and a pseudorandom start value SEED = 23 were selected based on the correlations displayed in Table S3. Figures S4 and S5 show the mass fractions of the 5-factorial solution as a function of fpeak and SEED values tested, respectively. The ratio of Q (parameter to be minimized by the PMF algorithm10) over Qexp (expected Q, degree of freedom of the fitted data is) 1.18 for five factors, indicating a good error estimation.12,38 The 4-factor solution was discarded due to one factor being a combination of COA and BBOA (Figure S6); the 6-factor solution yielded a “residual” factor (Figure S7). Mobile and stationary HOA, COA, and LV-OOA factor mass spectra are correlated with R2 values of 0.99, 0.92, and 0.96, respectively (Table S3). The O:C ratios of the mobile and stationary HOA factors are similarly low (0.01 and 0.03, respectively). For the COA, the stationary factor exhibits a higher O:C ratio (0.21) than the mobile factor (0.1), indicating that measurements on-road were closer to the cooking sources with lower atmospheric residence time and lower exposure to oxidants. The mobile LV-OOA factor has a higher O:C ratio (0.98) than the stationary LV-OOA factor (0.75), likely due to mobile measurements including more aged and processed air masses downwind of the city center. The BBOA spectra (R2 = 0.77) exhibit differences in the signals at m/z 44 (and 28), which are higher in the mobile spectrum. This is also reflected in the O:C ratios (0.31 and 0.24 for the mobile and stationary factor profile, respectively). The SV-OOA spectra show the least similarities (R2 = 0.69), mostly due to the signal of the CHO ion

(e.g., relative humidity or radiation) can influence individual particle components differently. For PMF factors (Figure 2), the Org normalization factor was used (normalization was done after PMF analysis). Normalization factors per drives and species are presented in the SI, as well as a detailed discussion on the influence of the normalization on zonal mean values, and limitations of the method employed. We ask the reader to keep in mind that the mass concentrations reported in Figures 3 and 4 are normalized data. In addition, the data were binned into different time of day intervals to give an indication of the influence of the diurnal pattern of parameters on the spatial variation. Table 1 specifies the time of day bins and the prevailing atmospheric scenarios as defined by Pandolfi et al.26 which influenced meteorological conditions, strength of the sea breeze cycle, mixing/residual layer height, and thus pollutant concentrations in the region. Briefly, scenario A is characterized by the most severe pollution episodes, recirculation of air masses, the precoastal mountain range lying within the mixing and residual layers, and a weakened sea breeze cycle; scenario C by North Atlantic advection flushing the planetary boundary layer; T denotes the transition periods with strong sea breeze cycles and a decoupling of pollutant concentrations between the city center and the precoastal mountain range.



RESULTS AND DISCUSSION PMF Results. Five OA factors were identified in the mobile data set, similar to the parallel stationary AMS data from the urban background site.21 These factors are shown in Figure 2, and are identified as the following: Hydrocarbon-like OA (HOA) from traffic (mostly diesel) emissions; biomass burning OA (BBOA) mostly from open agricultural fires;31 cooking OA (COA); and two different oxygenated OA (OOA), related to secondary OA: Low-volatility OOA (LV-OOA), and semivolatile OOA (SV-OOA). The R2 values of the correlations of the stationary and the mobile factor profiles, and of the time series of the PMF factors and ancillary data are shown in Table S3, as well 10425

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology (R2 = 0.86 for the linear regression fit excluding CHO). Again, the O:C ratio of the mobile factor is higher than that of the stationary factor (0.54 and 0.32, respectively). Correlation of EBC with HOA (R2 = 0.63) is better than with BBOA (R2 = 0.24) and is likely due in part to sample bias from on-road sampling and variability in emission of BC from biomass combustion sources39 (Table S3). This indicates that the majority of EBC sampled originated from traffic related fossil fuel combustions (supported by the R2 = 0.71 of the correlation of HOA and CO, and R2 = 0.52 for HOA and NO2). These values are comparable to correlation factors found for measurements of traffic-related pollutants in Queens, New York with a soot particle AMS (R2 = 0.5 and 0.6 for HOA (including refractory BC) and BC or NOx, respectively).40 Usually in ambient data sets the time series of LV-OOA and SO4 are well correlated due to their regional characteristics.11 This is not the case for the present data (R2 = 0.30). We speculate that shipping emissions which occur spatially and temporally separated from urban emissions locally influence the SO4 concentrations.41,42 LV-OOA exhibits a temporal evolution (R2 = 0.40) similar to, e.g., the organic fraction of m/z 44 (f44), a surrogate for the degree of oxidation and thus photochemical age of the air mass.43 SV-OOA correlates well with NO3 (R2 = 0.74) due to its volatile nature, as observed in other data sets, for example, refs 11 and 21. Regional Overview. Figure 3 gives an overview of the normalized mean concentrations of PM1 components measured on-road per zone and time interval. For a comparison of normalized to measured mean values, see Table S2. Corrected particle number and gas phase concentrations are given in Table S4. Apparent from Figure 3 is the relatively uniform distribution of PM1 component concentrations across the Barcelona metropolitan region. EBC is the major component in all zones, with the highest zonal mean value measured in the heavily populated Vallés region (zone 5) followed by the city center (zone 1). These EBC values are biased high due to the influence of surrounding vehicles when measuring on-road. This influence is also visible in the large differences of EBC mean zonal values between different times of day (horizontal bars). In addition, urban EBC background concentrations might be more reduced compared to on-road concentrations (due to tailpipe emissions being a major source) than for AMS species. This effect is even more distinct during daytime when traffic activity is high. Mobile measurements were only performed during the day. The normalization to relatively lower background data will thus bias EBC values high relative to AMS species and make the direct comparison difficult. Consequently, zonal mean EBC values should only be interpreted qualitatively for concentration levels in different zones, but at the same time are indicative of high emissions in the entire region. The inorganic PM1 components NO3, NH4, SO4, and Chloride (Chl) show little variability in concentration and fractional contribution across the Barcelona metropolitan area. An exception are the industrial zones of the Llobregat valley, where Chl shows significantly higher concentrations, most likely due to HCl emissions from industries44 (not discernible in Figure 3, compare Figure S8). The ammonium nitrate/ammonium sulfate ratio is 1 day) air masses in other studies (pink slope in Figure 5, for example, refs 47−49). Since the aging time of this case study is much less than a day (∼2.5 h), the secondary OA mass formed downwind of the Barcelona city center might not be purely of anthropogenic origin. In fact, Minguillón et al.32 found that at Montseny, a measurement station during DAURE on the same precoastal mountain range, ∼ 72% of secondary OA carbon was nonfossil and from biogenic, biomass burning, and urban nonfossil sources (e.g., cooking). Ignoring that urban nonfossil sources also count as anthropogenic precursors, and taking only the fossil 28% of secondary OA into account, yields a slope of OA/CO of ∼24 μg m−3 ppmv−1 (black slope in Figure 5). Kleinman et al.32 calculated OA/CO ratios for different photochemical lifetimesa rough comparison shows that a slope of 24 μg m−3 ppmv−1 corresponds to ∼0.1 unit of photochemical age which in turn corresponds to ∼2.5 h, in good agreement with the time difference of 2 h between the start and end point of this measurement drive. The calculated secondary OA mass yield of 75 μg m−3 ppmv−1 (difference between ∼85 μg m−3 ppmv−1 at Montseny and ∼10 μg m−3 ppmv−1 in the city center) was compared to the mass increase in SV- and LV-OOA (normalized to the fit CO plume mean values scaled to 1, Figure S9) between start and end point of the measurement drive and showed good agreement as well (Figure S9b). Overall, the spatial distribution of on-road measured PM1 components shows that during March 2009, the entire Barcelona metropolitan area was heavily impacted by local emissions from traffic, cooking activities, and, to a lesser extent, biomass burning. Industrial areas, the Barcelona airport, harbor, an extensive highway network and agricultural areas surround the city center. During the measurements, several open fires were observed by the drivers (zone 3). The municipality of Barcelona has one of the highest car densities per km2 in Europe,50 and (primary) road traffic emission effects also include high concentrations of NO2 (48 ppbv), CO (556 ppbv), and particle number concentrations (115 000 cm−3). Note that since measurements were done on road, the concentrations of compounds related to traffic emissions will be enhanced compared to background concentrations, but may be indicative of near-road concentrations. The formation of secondary aerosol, mostly OA, favored by high precursor emissions of both anthropogenic and biogenic origin, and high photochemical activity in the Mediterranean climate adds more than half of the OA burden outside the city center. We conclude that the Barcelona metropolitan area is an important source of regional PM1, compared to, e.g., Paris, which has a very

photochemical activity from gaseous precursors emitted in the city center, which can be upwind of the harbor at this time of day.28 The high LV-OOA contribution to OA (32%, 2.2 μg m−3) may be partly due to aging of SV-OOA and advection of processed air masses by the sea breeze. In addition, the two afternoon drives in zone 2 were both during meteorological episode A, characterized by the most severe pollution episodes and high LV-OOA levels due to recirculation of air masses. The LV-OOA concentrations might thus be biased high despite normalization to urban background values. Zone 6 (Montserrat), a forested, mountainous, and sparsely populated region, exhibits the second-highest contribution of OOA. In contrast to zone 2, t1, LV-OOA (42%, 2.9 μg m−3) is higher than SV-OOA (21%, 1.5 μg m−3). From the city center to zone 6, there is a general trend in increasing contributions from (LV-)OOA, indicative of formation and aging of OA in the plume, as it evolves in the atmosphere downwind of the city. This will be investigated in the next section. OOA Formation in the Urban Plume. The formation of secondary OA from anthropogenic precursors in the plume downwind of a city can be investigated by the ratio of OA to CO.22,47−49 Generally, ΔCO (the CO concentrations minus a background value) is used as an inert tracer of urban emissions including aerosol precursors to account for dilution effects (evaporation upon dilution is assumed not to affect this relationship24,49). We use the portion of the measurement drive on March 20 from zone 1 (city center) to zone 6 (Montserrat, precoastal mountain range) to estimate anthropogenic secondary OOA formation downwind of the Barcelona city center. The meteorological conditions were governed by a transition period with low mixing layer height and a strong sea breeze effect.26 The measurement started at ∼14:00 local time, well after the sea breeze was established, at the shore upwind of the city center and reached Montserrat, with 711 m. a. s. l. the highest point in the measurement route, at 16:05 local time (see portion of the route highlighted in red, Figure 1). Figure 5 shows the measured OA plotted vs CO concentrations. Note that no CO background was subtracted here (issues with CO back-

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 Parameters a and b denote slopes and intercepts (only where data were fit, see text). 10427

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology low impact on the air quality of its surroundings.51 In addition, the diurnal cycle of mountain wind and sea breeze can lead to recirculation of processed air masses mixing with new emissions, possibly causing slight changes in the relative contributions of the various components as a function of time of day. However, since this affects both the city center and the surrounding areas, and since the whole region is heavily populated and has abundant primary emissions sources, the overall concentrations and composition of PM1 are relatively similar across the whole area.



(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. Proc. Natl. Acad. Sci. U. S. A. 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 twodimensional 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.; 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.; Kiendler-Scharr, 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ä,

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.est.5b02149. Experimental section zone identification; instruments; normalization of mobile data; results PMF results; regional overview; and additional references (PDF)



AUTHOR INFORMATION

Corresponding Author

*Phone: +49−721/608−2 2696; fax: +49-721/608-2 4332; email: [email protected] (C.M.). Present Addresses §

Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology, Karlsruhe, Germany. ∥ Department of Civil, Architectural, and Environmental Engineering and Department of Chemistry, Drexel University, Philadelphia, PA, U.S.A. ⊥ Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA), UTAPRAD-DIM, Frascati, Italy. # European Commission Joint Research Centre (JRC), Institute for Environment and Sustainability, Ispra, Italy. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We thank the organizers of the DAURE project and the staff at IDAEA for their support. We acknowledge the IMBALANCE project of the Competence Center Environment and Sustainability (CCES) and the EU-FP7 project EUCAARI for financial support and “Accion Complementaria DAURE” from the Spanish Ministry of Science and Innovation (CGL2007-30502E/CLI) for infrastructure support. P.F.D. is grateful for the postdoctoral support from the US-NSF (IRFP# 0701013).



REFERENCES

(1) Gurjar, B. R.; Nagpure, A. S.; Singh, T. P. Air quality in megacities. In Encyclopedia of Earth, Cutler J. Cleveland, Environmental Information Coalition; Hanson, H., Eds.; 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. 10428

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology 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), 10453−10471. (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. (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 ACE-Asia using an Aerodyne aerosol mass spectrometer. J. Geophys. Res. 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. Assoc. 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), 4953−5001. (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.; GomezMoreno, 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. (38) Paatero, P.; Hopke, P. K.; Song, X. H.; Ramadan, Z. Understanding and controlling rotations in factor analytic models. Chemom. Intell. Lab. Syst. 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 10429

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430

Article

Environmental Science & Technology 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. Res. 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 high-resolution 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 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.

10430

DOI: 10.1021/acs.est.5b02149 Environ. Sci. Technol. 2015, 49, 10421−10430