Evaluation of the absorption Ångström exponents for trafﬁc and wood burning in the Aethalometer-based source apportionment using radiocarbon measurements of ambient aerosol

. Equivalent black carbon (EBC) measured by a multi-wavelength Aethalometer can be apportioned to traf-ﬁc and wood burning. The method is based on the differences in the dependence of aerosol absorption on the wave-length of light used to investigate the sample, parameterized by the source-speciﬁc absorption Ångström exponent ( α) . While the spectral dependence (deﬁned as α values) of the trafﬁc-related EBC light absorption is low, wood smoke particles feature enhanced light absorption in the blue and near ultraviolet. Source apportionment results using this methodology are hence strongly dependent on the α values assumed for both types of emissions: trafﬁc α TR , and wood burning α WB . Most studies use a single α TR and α WB pair in the Aethalometer model, derived from previous work. However, an accurate determination of the source speciﬁc α values is currently lacking and in some recent publications the applicability of the Aethalometer model was questioned. Here we present an indirect methodology for the determination of α WB and α TR by comparing the source appor-tionment of EBC using the Aethalometer model with 14 C measurements of the EC fraction on 16 to 40 h ﬁlter samples from several locations and campaigns across Switzer-land during 2005–2012, mainly in winter. The data obtained at eight stations with different source characteristics also enabled the evaluation of the performance and the uncertainties of the Aethalometer model in different environments. The best combination of α TR and α WB (0.9 and 1.68, respectively) was obtained by ﬁtting the Aethalometer model outputs (calculated with the absorption coefﬁcients at 470 and 950 nm) against the fossil fraction of EC (EC F / EC) derived from 14 C measurements. Aethalometer and 14 C source apportionment results are well correlated ( r = 0.81

Abstract.Equivalent black carbon (EBC) measured by a multi-wavelength Aethalometer can be apportioned to traffic and wood burning.The method is based on the differences in the dependence of aerosol absorption on the wavelength of light used to investigate the sample, parameterized by the source-specific absorption Ångström exponent (α).While the spectral dependence (defined as α values) of the traffic-related EBC light absorption is low, wood smoke particles feature enhanced light absorption in the blue and near ultraviolet.Source apportionment results using this methodology are hence strongly dependent on the α values assumed for both types of emissions: traffic α TR , and wood burning α WB .Most studies use a single α TR and α WB pair in the Aethalometer model, derived from previous work.However, an accurate determination of the source specific α values is currently lacking and in some recent publications the applicability of the Aethalometer model was questioned.
Here we present an indirect methodology for the determination of α WB and α TR by comparing the source apportionment of EBC using the Aethalometer model with 14 C measurements of the EC fraction on 16 to 40 h filter samples from several locations and campaigns across Switzerland during 2005-2012, mainly in winter.The data obtained at eight stations with different source characteristics also enabled the evaluation of the performance and the uncertainties of the Aethalometer model in different environments.The best combination of α TR and α WB (0.9 and 1.68, respectively) was obtained by fitting the Aethalometer model outputs (calculated with the absorption coefficients at 470 and 950 nm) against the fossil fraction of EC (EC F / EC) derived from 14 C measurements.Aethalometer and 14 C source apportionment results are well correlated (r = 0.81) and the fitting residuals exhibit only a minor positive bias of 1.6 % and an average precision of 9.3 %.This indicates that the Aethalometer model reproduces reasonably well the 14 C re-Published by Copernicus Publications on behalf of the European Geosciences Union.P. Zotter et al.: Evaluation of the absorption Ångström exponents sults for all stations investigated in this study using our best estimate of a single α WB and α TR pair.Combining the EC, 14 C, and Aethalometer measurements further allowed assessing the dependence of the mass absorption cross section (MAC) of EBC on its source.Results indicate no significant difference in MAC at 880 nm between EBC originating from traffic or wood-burning emissions.Using EC F / EC as reference and constant a priori selected α TR values, α WB was also calculated for each individual data point.No clear stationto-station or season-to-season differences in α WB were observed, but α TR and α WB values are interdependent.For example, an increase in α TR by 0.1 results in a decrease in α WB by 0.1.The fitting residuals of different α TR and α WB combinations depend on EC F / EC such that a good agreement cannot be obtained over the entire EC F / EC range using other α pairs.Additional combinations of α TR = 0.8, and 1.0 and α WB = 1.8 and 1.6, respectively, are possible but only for EC F / EC between ∼ 40 and 85 %.Applying α values previously used in the literature such as α WB of ∼ 2 or any α WB in combination with α TR = 1.1 to our data set results in large residuals.Therefore we recommend to use the best α combination as obtained here (α TR = 0.9 and α WB = 1.68) in future studies when no or only limited additional information like 14 C measurements are available.However, these results were obtained for locations impacted by black carbon (BC) mainly from traffic consisting of a modern car fleet and residential wood combustion with well-constrained combustion efficiencies.For regions of the world with different combustion conditions, additional BC sources, or fuels used, further investigations are needed.

Introduction
Recently, the World Health Organization (WHO) reported around 3.7 million premature deaths in 2012 as a result of exposure to ambient air pollution, demonstrating that health risks in areas of low air quality are far greater than previously thought (WHO, 2014).Atmospheric particulate matter (PM) contributes significantly to ambient air pollution and adversely affects human health causing respiratory and cardiopulmonary diseases associated with increased morbidity and mortality (Pope and Dockery, 2006;WHO, 2006).Although PM levels were decreasing in the last decade in Europe and also in Switzerland, legal thresholds are still exceeded (Barmpadimos et al., 2011(Barmpadimos et al., , 2012)).Carbonaceous material (total carbon, TC) is a major fraction of the fine aerosol mass (up to 90 % of the PM mass < 2.5 µm, Gelencsér, 2004;Putaud et al., 2004;Jimenez et al., 2009) and is further classified into the sub-fractions organic carbon (OC) and black carbon (BC) or elemental carbon (EC; Jacobson et al., 2000).BC is the light-absorbing part of carbonaceous material and, compared to other aerosol components, it contributes significantly to global warming due to its optical and radiative properties (Jacobson, 2001(Jacobson, , 2010;;IPCC, 2013).Because of the relatively short atmospheric lifetime of BC, its radiative forcing ends within weeks after emission.Thus, reducing BC emissions may rapidly reduce climate warming (Shindell et al., 2012;Bond et al., 2013 and references therein).Therefore, the identification of different BC sources and their emission strength is crucial for the implementation of effective mitigation strategies.
The emission sources of BC are combustion processes of fossil and non-fossil carbonaceous fuels.In Switzerland, large parts of Europe and other parts of the world, BC mainly originates either from traffic or biomass burning in winter (e.g., Szidat et al., 2007;Favez et al., 2010;Lanz et al., 2010;Piazzalunga et al., 2011;Harrison et al., 2012;Larsen et al., 2012;Crippa et al., 2013;Herich et al., 2014).Different methods exist to quantify carbonaceous aerosol fractions based on light absorption, thermo-optical or laser-induced incandescence measurements.The quantities measured are defined based on the instrument and protocol used, with BC and EC related to optical and thermo-optical as well as chemical measurements, respectively.When BC is obtained by light absorption measurements it is referred to as mass equivalent black carbon (EBC; Petzold et al., 2013).In recent years, the Aethalometer, an online measurement technique of the aerosol light absorption at seven different wavelengths ranging from near-ultraviolet (N-UV) to near-infrared (N-IR), has become widely used, since it is rather inexpensive, portable, easy to operate and suitable for long-term measurements.Furthermore, multi-wavelength Aethalometer data may be used to derive the traffic and the wood-burning contributions to EBC (EBC TR and EBC WB , respectively) taking advantage of the light absorption in the blue and N-UV of aerosols from biomass-combustion likely due to coemitted organics, which is enhanced compared to aerosols from fossil sources (Sandradewi et al., 2008a).The so-called "Aethalometer model" assumes that light-absorbing particles only originate from vehicle and biomass-burning emissions, and uses absorption Ångström exponent (α) values specific to these sources to derive their contributions.Therefore, the source apportionment of EBC using the Aethalometer model is inherently dependent on the a priori assumed absorption Ångström exponents for traffic (α TR ) and biomass burning (α WB ), which are based on a few emission studies.α TR values cluster in a narrow range (0.8-1.1), whereas a large range of α WB values (0.9 to 3.5) is reported (Schnaiter et al., 2003(Schnaiter et al., , 2005;;Kirchstetter et al., 2004;Lewis et al., 2008;Saleh et al., 2013).Some studies also obtained α TR from ambient Aethalometer measurements by investigating the α values calculated from the ambient absorption coefficient (b abs ) values of the total light-absorbing aerosol during periods and locations that were only influenced by traffic emissions (e.g., in summer close to roads, Sandradewi et al., 2008b;Herich et al., 2011).It should be noted, however, that α values depend not only on different emission sources but also on the choice of wavelengths and different calculation procedures used, although deviations due to the latter are expected to be low (Moosmüller et al., 2011).
Another independent and more direct approach than the Aethalometer model to distinguish between modern (wood burning) and fossil (traffic) contributions is the radiocarbon analysis.Radiocarbon ( 14 C) is completely depleted in fossil fuel emissions ( 14 C half-life = 5730 years) and can, therefore, be separated from non-fossil carbon sources, which have a similar 14 C signal as atmospheric carbon dioxide (CO 2 ; Currie, 2000;Szidat, 2009).Measuring 14 C in the EC fraction therefore directly enables the quantification of the wood-burning and fossil sources of EC.However, the 14 C analysis can only be performed on filter samples and is therefore limited in time resolution.Furthermore, such analysis is rather expensive and time consuming.The 14 C measurement in the EC fraction remains additionally challenging in contrast to TC (Szidat et al., 2013), since a clear physical separation between OC and EC is necessary to avoid interferences from OC in the 14 C signal.Nevertheless, recent developments and method adaptations from different groups show more consistent approaches and yield more robust 14 C results (Zhang et al., 2012;Bernardoni et al., 2013;Dusek et al., 2014).Sandradewi et al. (2008a) first employed the Aethalometer model on winter data from a polluted Swiss Alpine valley and used 14 C measurements of the EC fraction to test the assumed α WB and α TR .Based on this work, subsequent studies using the Aethalometer model employed similar α TR (0.9-1.1) and α WB (1.8-2.2) values (e.g., Sandradewi et al., 2008b;Favez et al., 2010;Perron et al., 2010;Herich et al., 2011;Harrison et al., 2012Harrison et al., , 2013;;Crippa et al., 2013;Mohr et al., 2013), without further evaluations of these parameters using external data.Others compared the Aethalometer model outputs to apportionments using specific source tracers (Favez et al., 2010;Herich et al., 2014;Crilley et al., 2015).However, such approaches heavily rely on a priori assumed tracer-to-BC emission ratios, which are highly variable (Schmidl et al., 2008;El Haddad et al., 2011, 2013;Heringa et al., 2011), and as such are not suitable for the evaluation of the α TR and α WB values used in the Aethalometer model.Even though the Aethalometer model is widely used there are also studies (Harrison et al., 2013;Garg et al., 2016) that question the applicability of this model when other and/or additional combustion sources may contribute to the BC burden and combustion efficiencies are less well constrained.
In this study we present an evaluation of the Aethalometer model by comparing its outputs to 14 C results of the EC fraction in order to validate the choice of the absorption Ångström exponents for wood burning (α WB ) and traffic emissions (α TR ).To this end, we use 14 C and Aethalometer data from different campaigns across Switzerland, mostly from the winter season.The data set in this study (n = 101) is significantly larger than previously reported (n = 12 and n = 18 in Sandradewi et al., 2008a, b, respectively).In addi- tion, the data presented here were obtained at eight different stations in various area types with different source characteristics (e.g., urban, suburban, rural, Alpine valley, traffic, background, etc.) thereby enabling the evaluation of the performance and the uncertainties of the Aethalometer model in different environments.

Aerosol sampling
Aerosol sampling presented in this study (see Table 1) was carried out at several stations of the Swiss National (NABEL) and Cantonal air pollution monitoring networks (EMPA, 2013;Cercl'Air, 2012).The stations ZUR, PAY, REI, and SIS are located north of the Alps, whereas MAG, ROV, and MOL are located south of the Alps and MAS is situated in the Rhone valley.The location of these stations in different areas allowed the sampling of a broad range of particles, with different characteristics ranging from urban to rural and from traffic to background.The exact locations of the stations are shown in Fig. 1 and the details and full names of the sites as well as the different campaigns carried out at these stations are listed in Table 1.
Filter sampling was conducted using quartz fiber filters (Pallflex 2500QAT-UP) and high-volume samplers (Digitel DHA-80, Switzerland) at a flow rate of 30 m 3 h −1 .The collection time as well as the size cut of the PM inlet varied between 16 and 40 h as well as PM 1 and PM 10 , respectively, depending on the campaign (see Table 1).After sampling, filters were stored at −20 • C until analysis.Most of the results presented here (n = 69) were obtained on PM 10 filters with a sampling time of 24 h from the 5-year 14 C project Switzerland (Zotter et al., 2014).The samples from this campaign were collected on days with high PM 10 concentrations (almost all of them exceeding the Swiss and EU daily limit of  Herich et al., 2011;EMPA, 2013).The same type of instrument was also placed at SIS in the winters 2010/2011 and 2011/2012 during the 14 C project Switzerland and earlier campaigns in 2005 and 2006.In total 101 samples with parallel 14 C and Aethalometer measurements are available (n = 9, 24, 19, 19, 13, 4, 5, 8 for SIS, ZUR, MAG, PAY, ROV, MOL, REI and MAS, respectively).

Measurement principle
The Aethalometer provides a real-time optical measurement of light-absorbing carbonaceous aerosols at seven wavelengths (λ = 370, 470, 520, 590, 660, 880, and 950 nm;Hansen et al., 1984;Hansen, 2003).It measures the attenuation (ATN) of a light beam transmitted through a filter on which aerosols are continuously collected: where I 0 and I denote the intensity of a light beam through an empty and particle-laden spot of a filter tape, respectively.The change in ATN over a certain time period (t) is proportional to the attenuation coefficient (b ATN ) given a known flow rate (Q) and spot size (A) onto which particles are col-lected: Like all filter-based absorption techniques, the Aethalometer uses integration of the sample on the filter to increase the sensitivity of the measurement.Scattering by the filter fibers enhances absorption of the light by the aerosols collected on the filter tape.As the filter gets loaded by lightabsorbing aerosols and ATN increases, nonlinear loading effects become apparent (Liousse et al., 1993;Petzold et al., 1997;Bond et al., 1999;Park et al., 2010;Drinovec et al., 2015).To compensate for these effects, the algorithm developed by Weingartner et al. (2003) was used to derive the final absorption coefficient (b abs ): where C λ and R(f λ , ATN λ ) are factors to compensate for multiple scattering of the filter fibers and the loading effect, respectively: In Eq. (4) f λ denotes the slope between the linear function R(f λ , ATN λ ) vs. ln(ATN λ ) and allows estimating the instrumental error that occurs when the shadowing effect is disregarded (Weingartner et al., 2003).This approach is routinely applied to the Aethalometer data from the NABEL stations using a single C value of 2.14 for all wavelengths and wavelength-dependent f values (1.155, 1.137, 1.128, 1.116, 1.103, 1.064, and 1.051 for 370, 470, 520, 590, 660, 880, and 950 nm, respectively) as proposed by Weingartner et al. (2003) and Sandradewi et al. (2008c), respectively.
The same values were also used to compensate the data from SIS and the previous campaigns in Switzerland.Several other algorithms for the compensation of the Aethalometer data are available (Collaud Coen et al., 2010 and references therein) and some studies slightly adapted the Weingartner et al. (2003) approach (Sandradewi et al., 2008c;Favez et al., 2010;Mohr et al., 2013;Segura et al., 2014).It should be noted that these different compensation algorithms might yield slightly different b abs (λ).However, the comparison of these approaches or the improvement of the compensation methodology used is beyond the scope of this study.Also the recently developed dual spot Aethalometer (AE33, Drinovec et al., 2015) allows for an improved and time-dependent loading compensation.
The compensated b abs is then converted into a EBC mass using the mass absorption cross section (MAC): Usually nominal MAC values are used, to directly infer EBC mass from the non-compensated b ATN .These MAC values can be calculated from the parameters furnished by the Aethalometer manufacturer (Hansen, 2003) or are provided in the literature (e.g., Bond et al., 2013, and references therein) and include a C value.Here MAC values are obtained empirically by comparing b abs with simultaneous measurements of EC from thermo-optical methods (e.g., Moosmüller et al., 2001;Bond et al., 2013 and references therein), and the EBC concentration is assumed to be identical to the EC concentration.From Eqs.
(3) and ( 5) it is evident that empirically derived MAC values for absorption photometers strongly depend on the assumed C value.Different C values were previously empirically derived from instrumental comparisons and used to determine the absorption coefficient from Aethalometer measurements (e.g., Collaud Coen et al., 2010;Segura et al., 2014;Crilley et al., 2015).The separation of the C value and the MAC is therefore relative to the methods used, and empirically determined MAC values using Aethalometers should always be reported together with the applied C values (C = 2.14 in our case).

Source apportionment using Aethalometer data
The spectral dependence of the absorption is described by the power law b abs (λ) ∼ λ −α (Moosmüller et al., 2011), where α is the absorption Ångström exponent and consequently for a wavelength pair the following relation can be derived: BC is a strong broadband absorber over the entire visible wavelength range (N-UV to N-IR) with only a weak spectral dependence (α for BC ∼ 1).Traffic emissions mainly contain BC and basically no other light-absorbing compounds and consequently α for traffic emissions (α TR ) ∼ 1 (Bond et al., 2013;Kirchstetter et al., 2004;Schnaiter et al., 2003Schnaiter et al., , 2005)).Biomass-burning aerosols, on the other hand, contain additionally to BC a substantial fraction of light-absorbing organic substances which strongly enhance the light absorption in the N-UV and blue part of the spectrum and have no contribution in the N-IR wavelength range resulting in an α for biomass-burning emissions (α WB ) that is larger than α TR .Based on this, Sandradewi et al. (2008a) developed a two-component model to apportion b abs measured with the Aethalometer at different wavelengths into a wood-burning (WB) and a traffic (TR) contribution assuming that the total b abs is only influenced by these two sources: This assumption is valid for Switzerland and other Alpine regions in Europe, especially in winter, where emissions from other sources are negligible.Coal burning is not used in these areas (Eurostat, 2017) and biogenic secondary organic aerosol (SOA) is mostly absorbing in the UV range (Romonosky et al., 2016) not covered by wavelengths used in the Aethalometer (especially given that we recommend the use of the absorption at 470 nm rather than at 370 nm; see Sect.3.2.4).Mineral dust can usually be neglected in this region (contribution to total PM < ∼ 10 %; Gianini et al., 2012), and special events possibly influencing the absorption at Aethalometer wavelengths 470-590 nm can be identified due to a drop of the absorption Ångström exponent clearly below one during such events (Collaud Coen et al., 2004).Using Eqs. ( 6)-( 7) and the measured b abs at two different wavelengths, a traffic and wood-burning contribution can be apportioned using the following equations: The contributions of wood-burning and traffic to total EBC (EBC WB and EBC TR ) are then derived via the corresponding MAC values (MAC WB and MAC TR , respectively): Consequently the ratio EBC TR to total EBC (EBC TOT ) can be derived from the measured ratio b abs (λ 1 ) to b abs (λ 2 ) and Using Eqs. ( 3)-( 6), Eq. ( 13) can be written as λ 2 has to be a wavelength in the N-IR range, where EBC is the only light absorber, whereas λ 1 should be taken from the N-UV range where also organics contribute to the light absorption.In this model α WB and α TR have to be assumed a priori or determined comparing the contributions of EBC TR and EBC WB to other techniques which apportion BC or EC into those two sources (e.g., 14 C measurements).Additional uncertainties may arise from the compensation factors applied to the attenuation coefficients.In this study, a fixed C λ value was used for the multi-scattering correction (Sect.2.2.1) and thus the ratio C λ,1 / C λ,2 becomes unity in Eq. ( 14).This is justified and introduces very little uncertainty, as the wavelength dependence of the f and C values across the range λ = 470-950 nm was reported to be less than 10 and 12 %, respectively, for the Aethalometer model AE31 (Weingartner et al., 2003;Sandradewi et al., 2008c;Segura et al., 2014).If data from other photometer models, which exhibit a wavelength dependence of the C value, are used for the source apportionment, the correct ratio C λ,1 / C λ,2 must be used in Eq. ( 14) to ensure consistency of the Aethalometer model parameters.The loading compensation factor R(f λ , ATN λ ) depends on wavelength, even if f (λ) is independent of the wavelength, since the ATN depends considerably on the wavelength.Nevertheless, uncertainties in the EBC TR to EBC ratio associated with the filter-loading compensation can be kept small by carefully determining the f values, following the approach in Weingartner et al. (2003) or Sandradewi et al. (2008c).The Aethalometer AE33 measures the compensation parameters and therefore the compensation is performed on-line.The precision of this compensation can be checked using the EBC(ATN) or b abs (ATN) analysis (Drinovec et al., 2015).It should be noted that the calculation of the EBC TR to EBC ratio (Eq.13) might not only be sensitive to the choice of compensation parameters but also on the choice of compensation algorithm.However, large uncertainties of the EBC TR to EBC ratio due to the use of different Aethalometer data compensation algorithms are not expected since in Eqs. ( 13) and ( 14) only fractional contributions of b abs (λ) or b ATN (λ) are used.Therefore, only differences in the wavelength dependency of the compensation parameters in different compensation methods would slightly affect the determination of EBC TR / EBC.An investigation of such effects is beyond the scope of this study; however, future work should be carried out exploring possible influences of different compensation methodologies on EBC TR / EBC.A detailed comparison of the different Aethalometer compensation algorithms can be found in Collaud Coen et al. (2010) and only an overall assessment of the methodology used will be discussed below.Sandradewi et al. (2008a) and subsequent studies that used the Aethalometer model utilized the same MAC for traffic (MAC TR ) and wood burning (MAC WB ) emissions at the N-IR wavelength, based on the fact that MAC values for freshly generated EBC were previously found to fall within a relatively narrow range (Bond and Bergstrom, 2006 and references therein).However, MAC values depend on particle size, morphology, and mixing state and thus different values for biomass-burning and traffic emissions may be possible.Therefore, we assess the ratio of MAC TR to MAC WB for our data set in Sects.3.1 and 3.2.1.Sandradewi et al. (2008a) and many other studies used 470 and 950 nm as N-UV and N-IR wavelengths, respectively.However, also other combinations of wavelengths have been used (e.g., 370 and 880, or 470 and 880 nm; see Perron et al., 2010;Herich et al., 2011;and Fuller et al., 2014), especially in studies that performed Aethalometer measurements with the two-wavelength instrument (370 and 880 nm, model AE22, Magee Scientific).Therefore, we will also investigate the sensitivity of the Aethalometer model using different wavelength combinations.

Separation of the carbonaceous particle fractions
Two different methods to isolate EC for the 14 C analysis were used.For the samples from the 14 C project Switzerland, the Swiss_4S protocol was applied for the EC isolation using a Sunset OC / EC analyzer as described by Zhang et al. (2012).This approach is optimized such that biases in the 14 C result of EC due to OC charring or losses of the least refractory EC during the OC removal are minimized.In brief, to minimize positive artifacts from OC charring, water-soluble OC (WSOC) is first eliminated by a water extraction and the remaining water-insoluble OC (WINSOC) is then removed using the Sunset analyzer by a thermal treatment in three steps: (1) 375 • C for 150 s in pure oxygen (O 2 ); (2) 475 • C for 180 s in O 2 ; (3) 450 • C for 180 s followed by 180 s at 650 • C in helium.Finally, in a fourth step, EC is isolated by the combustion of the remaining carbonaceous material at 760 • C for 150 s in O 2 .The evolving CO 2 is separated from interfering gaseous products, cryo-trapped and sealed in glass ampoules for 14 C measurements.By using the Sunset analyzer, which monitors the transmission of light through the filter with a laser during the combustion, the quantification of OC charring and EC losses before the fourth step is achieved.For the samples of the 14 C project Switzerland, charred OC only contributed ∼ 5 % to EC recovered in the fourth step and on average 74 ± 11 % of the EC was recovered for the 14 C measurement.Charring OC of a given thermal step is quantified as the difference of the maximum ATN and the initial ATN normalized to the initial ATN.The EC recovery is related to the loss of EC during the first three steps and is defined as the ratio between the ATN of the laser signal through the filter before step 4 (EC step) and the initial ATN before the thermal treatment before the first step (Zhang et al., 2012).
EC from samples collected during the campaigns in ROV, MOL, MAS, REI, ZUR in 2005 and 2006 was isolated for the 14 C analysis using the THEODORE system and the approach described by Szidat et al. (2004).In brief, after removal of WSOC by water extraction, WINSOC was evaporated during 4 h in a muffle furnace in air at 375 • C. EC was finally combusted in the THEODORE system at 640 • C for 10 min with O 2 .The evolving CO 2 was recovered in the same manner as described above.The EC recovery for these samples was estimated by Zhang et al. (2012) and was on average 60 ± 12 %.The 14 C results of EC were corrected to 100 % EC recovery (see Sect. 2.4.3 below) and results obtained with the THEODORE and the Swiss_4S method were previously found to agree within the uncertainties (see Zhang et al., 2012).

Radiocarbon measurement
The analysis of the 14 C content in the CO 2 from the separated EC fraction collected as described above was carried out with the MIni radio CArbon DAting System, MICADAS (Synal et al., 2007) at the Swiss Federal Institute of Technology (ETH) Zurich and the Laboratory for the Analysis of Radiocarbon with AMS (LARA), University of Bern (Szidat et al., 2014) using a gas ion source (Ruff et al., 2010;Wacker et al., 2013).The results of the 14 C measurement are presented as fraction of modern (f M ) denoting the 14 C / 12 C content of the sample related that of the reference year 1950 (Stuiver and Polach, 1977).The f M values are corrected for δ 13 C fractionation and for 14 C decay between 1950 and the year of measurement (Wacker et al., 2010).The f M measurement uncertainty for the EC samples from the 14 C project Switzerland and ROV, MOL, MAS, REI, ZUR from 2005 and 2006 is ∼ 2 % (Zotter et al., 2014) and ∼ 3 % respectively (Zhang et al., 2012).

Determination of the non-fossil fraction of EC
As shown above (see Sect. 2.2), on average only 74 ± 11 and 60 ± 12 % of the total EC (EC yield) was isolated for the 14 C measurement of the samples from the 14 C project Switzerland and ROV, MOL, MAS, REI, ZUR from 2005 and 2006, respectively.However, Zhang et al. (2012) showed that f M values are lower for lower EC yields suggesting that the EC that is removed before the fourth step (the step in which EC is recovered for the 14 C measurement), is mainly from biomass burning due to its lower thermal stability (Zhang et al., 2012).Therefore, an extrapolation of the measured EC f M values to 100 % EC yield was applied to account for this underestimation of f M (Zhang et al., 2012).This method was applied to all samples discussed here, and the detailed description of the procedure used for the samples from the 14 C project Switzerland and ROV, MOL, MAS, REI, ZUR from 2005 and 2006 can be found in Zotter et al. (2014) and Zhang et al. (2012), respectively.
The f M of contemporary carbon including biogenic sources and biomass burning (f M,bio and f M,WB , respectively) is characterized by values of 1 whereas f M is equal to 0 for fossil sources due to the decay of 14 C with a half-life of 5730 years.Due to the nuclear weapon tests in the 1950s and 1960s, however, the atmospheric 14 C content increased and f M exhibits values > 1 (Levin et al., 2010).Therefore, f M values for EC were converted into non-fossil fractions (f NF,EC ; Szidat et al., 2006) using a reference value.Since biomass burning is the only non-fossil source of EC (neglecting possible small contributions from bio-fuels) this reference value is equal to f M,WB and was estimated using a tree-growth model as described in Mohn et al. (2008) including 10, 20, 40, 70, and 85-year old trees with weight fractions of 0.2, 0.2, 0.4, 0.1, and 0.1, respectively, harvested 3 years before aerosol sampling.Values of 1.140, 1.135, 1.127, 1.123, 1.119, 1.114, and 1.106 were calculated and consequently used to correct the f M values extrapolated to 100 % EC yield from samples collected in 2005, 2006, 2008, 2009, 2010, 2011, and 2012, respectively.The final uncertainties for f NF,EC (∼ 5 and ∼ 6 % for samples from the 14 C project Switzerland and ROV, MOL, MAS, REI, ZUR from 2005 and 2006, respectively) are derived from an error propagation and include all the individual uncertainties of f M (measurement uncertainty, extrapolation to 100 % EC yield) and f M,WB (Zotter et al., 2014).

Elemental carbon measurement
The EC concentrations on samples from the 14 C project Switzerland (see Table 1) were measured using a thermooptical OC / EC analyzer (Model 4L, Sunset Laboratory Inc., USA), equipped with a non-dispersive infrared (NDIR) detector following the thermal-optical transmittance protocol (TOT) EUSAAR2 (Cavalli et al., 2010).EC concentrations from the campaigns in ROV, MOL, MAS, REI, ZUR in 2005 and 2006 (see Table 1) are not included for the MAC calculations, since in earlier campaigns they were not measured or obtained with a different TOT protocol.We assigned a high uncertainty of 25 % for all measured EC concentrations to account for possible differences between different TOT protocols (Schmid et al., 2001) MAC determination is affected by the uncertainty of the EC concentrations whereas the evaluation of the choice of α WB and α TR using the fossil fraction of EC is influenced by the combined uncertainty of the 14 C measurement of EC, the extrapolation of f M,EC to 100 % EC yield and the bomb peak correction which was on average only 5-6 % (see Sect. 2.3).
No EC was detected on blank filters and consequently no blank correction was necessary (see also Zotter et al., 2014).

Additional data
Nitrogen oxides (NO x ) are routinely measured at the NABEL stations ZUR, MAG, and PAY using reference instrumentation with molybdenum converters according to valid European standards (EMPA, 2013).Since no large sources of NO x (e.g., fossil fuel power plants) are present in Switzerland besides traffic, NO x will be used here for the comparison with EBC TR (see Sect. 3.3 below).
Levoglucosan, a thermal degradation product of cellulose and thus a tracer for primary emissions of organic aerosol from biomass burning and often used to estimate OC mass from this source (Gelencsér et al., 2007), was also measured on 52 samples presented in this study.A description of the measurement details can be found in the corresponding references as listed in Table 1.Levoglucosan data are available for most of the samples from winter 2005 and 2006 from ROV, MOL, REI, MAS, and ZUR (n = 27) as well as from the winter 2008/2009 for ZUR, MAG, and PAY (n = 8) from the 14 C project Switzerland (see Zotter et al., 2014).In addition, data from these three stations (n = 17) with parallel Aethalometer measurements available were also taken from Herich et al. (2011).Levoglucosan data will be used here for the comparison with EBC WB (see Sect. 3.3 below).As photochemical degradation of levoglucosan was previously observed under summertime conditions (Kessler et al., 2010;Hennigan et al., 2011), spring and summer levoglucosan data from ZUR are not used here.

MAC determination
MAC values are determined empirically by comparing b abs with EC thermo-optical measurements (see Fig. 2a).b abs at 880 nm and EC are strongly correlated (r = 0.86) and the geometric mean of the MAC at 880 nm was found to be 11.8 m 2 g −1 (9.2-15.1 m 2 g −1 ), similar to values obtained in Herich et al. (2011) for ZUR (10.0 m 2 g −1 ), PAY (13.2 m 2 g −1 ), and MAG (9.9 m 2 g −1 ) for a 2-year data set (note that the MAC values reported in this study as well as that by Herich et al. (2011) both apply for EC mass based on the thermal-optical transmittance protocol EUSAAR2 and absorption coefficients inferred from Aethalometer AE31 data with assuming a C value of 2.14).No systematic yearto-year or station-to-station variations in the MAC values at 880 nm are observed.While the MAC values determined at SIS are lower on average, they remain within the previously reported range, and given the relatively modest number of samples, this observation cannot be generalized.It should be noted that MAC values depend on the aerosol mixing state, size, and morphology (see, e.g., Bond and Bergstrom, 2006), and empirically derived MAC values also depend on the limitations of the measurement techniques used to determine b abs .The results of our study would translate to ∼ 9.7-10.0m 2 g −1 at 637 nm when recalculating our MAC values from a wavelength of 880 to 637 nm with an absorption Ångström exponent of 0.9-1.0 and if a C value of 3.5 instead of 2.14 was assumed.This is in good agreement with the average MAC value of 10.0 m 2 g −1 at 637 nm reported by Zanatta et al. (2016) for nine European background sites, who also used the EUSAAR2 protocol for EC mass and either multi-angle absorption photometers, particle soot absorption photometers, or Aethalometers with assuming C = 3.5 for the absorption coefficient.Deviations from other previously reported MAC values at similar wavelengths (−5-26 m 2 g −1 , Liousse et al., 1993;Bond and Bergstrom, 2006;Genberg et al., 2013) can be due to different methods used to determine EC and the absorption coefficient and/or possible differences in BC size and mixing state.
Only few studies attempted the empirical determination of MAC values for biomass-burning and traffic EBC emissions using ambient measurements (e.g., Laborde et al., 2013;Bond et al., 2013 and references therein).Since the ratio of MAC TR to MAC WB at the N-IR wavelength is needed in the Aethalometer model (see Eq. 13), it is important to assess possible differences between MAC TR and MAC WB .Sandradewi et al. (2008a) and all other studies that applied the Aethalometer model assumed, implicitly or explicitly, a MAC TR to MAC WB ratio of unity at 880 nm.Having an independent measurement for the relative contributions of traffic and wood-burning to total EC from the 14 C measurements allows us to test this assumption by plotting the MAC values at 880 nm against the corresponding relative traffic contribution to EC (EC F / EC) obtained with the 14 C measurements (see Fig. 2b).No correlation between the two parameters was found, indicating that it is justified to simplify the Aethalometer model (Eq.12) and set the ratio of MAC TR to MAC WB at the N-IR wavelength to unity.This is in agreement with Herich et al. (2011) who did not find differences in MAC for the stations ZUR, MAG and PAY for a 2-year data set between summer and winter, where there is a large seasonality in the relative wood-burning contribution.The variability in Fig. 2b is due to day-to-day and station-to-station variability but could to some degree also originate from different size cuts (PM 10 or PM 1 and PM 2.5 ) of the filter samplers and Aethalometer measurements for some campaigns (see Table 1).Alternatively, the ratio of MAC TR to MAC WB at the N-IR wavelength can be used as a third free parameter, besides α TR and α WB , when fitting the Aethalometer model (Eq.12) against a data set of independent EC F / EC mea- surements.We tested this for the data set of this study and obtained a MAC TR to MAC WB ratio of 0.97, which confirms the finding of Fig. 2b.Therefore, in the following we will use a fixed MAC TR to MAC WB ratio of 1 in the Aethalometer model.

Best α TR and α WB pair, and analysis of uncertainties and biases
Independent measurements of the contribution of wood burning and traffic to BC (or EC) are often not available; therefore in most studies a single α TR and α WB pair is usually used in the Aethalometer model, derived from previous work.However, α WB and α TR may be highly variable, depending on the combustion conditions and efficiency, fuel type and aerosol aging (Lack et al., 2013;Saleh et al., 2013Saleh et al., , 2014;;Zhong and Jang, 2014;Sharpless et al., 2014;Kirchstetter et al., 2004;Bond and Bergstrom, 2006 and references therein;Herich et al., 2011;Garg et al., 2016).In this section we use EC F / EC values from 14 C measurements to determine the best combination of α TR and α WB and assess the performance of the Aethalometer model using this single pair of α values.In practice, the best pair of α values is determined by fitting Eq. ( 13) against EC F / EC from the 14 C analyses using the ratio b abs,470 / b abs,950 from the Aethalometer as independent variable (and assuming MAC TR,950 / MAC WB,950 = 1, as justified in Sect.3.1).We use a least-square fitting weighted by the inverse number of data points in EC F / EC bins of 0.1 as most of the data presented in this study fall within a range of EC F / EC = 0.4-0.6.The absorption Ångström exponents α TR and α WB that fit best our data were found to be 0.90 and 1.68, respectively.
The same α values were obtained when MAC TR / MAC WB was included as a third fitting parameter, because the best-fit MAC ratio is 0.97, which is virtually equal to unity (see also Sect.3.1).EBC TR / EBC at 950 nm, derived with the above bestfit Aethalometer model parameters, and EC F / EC are well correlated (r = 0.81; see Fig. 3a) and the fitting residuals ( EBC TR / EBC = EBC TR / EBC-EC F / EC, Fig. 3b) are normally distributed with only a minor positive bias of 1.6 %.We estimate that the precision of the model ( EBC TR / EBC) is on average 9.3 % in our case, using the standard deviation (σ ) of the Gaussian fit of EBC TR / EBC in Fig. 3b.This indicates that the Aethalometer model reproduces reasonably well the 14 C results for all stations investigated in this study using our best estimate of a single α WB and α TR pair.Since this analysis includes data from urban stations as well as from spring and summer this shows that the Aethalometer model also works for other areas than for polluted Alpine valleys in winter.It should be noted that the determination of EBC TR / EBC using the fitted α values cannot be more accurate than the uncertainty of EC F / EC.The estimated EBC TR / EBC is affected by (1) random measurement uncertainties of EC F / EC and b abs,470 and b abs,950 and (2) day-to-day and stationto-station variability in α WB and α TR values.Investigating the effect of a MAC TR,950 / MAC WB,950 different from one (MAC TR / MAC WB = 0.7-1.3) it is evident that there is no large influence on α WB (1.66-1.71),α TR (0.8-0.95), the mean bias (0.2-2.4 %), and EBC TR / EBC (9.4-9.9 %).This further justifies fixing the MAC ratio at unity when applying the Aethalometer model.
Without an alternative method for the source apportionment of EC or BC, the determination of α values and related uncertainties is unattainable.Therefore, we determined the distribution of α WB values for our data set and investigated whether there are other combinations of α TR and α WB that yield similarly acceptable agreement with 14 C data.For this  13) against EC F / EC.The error bars for EC F / EC represent the combined uncertainty of the 14 C measurement of EC, the extrapolation of f M,EC to 100 % EC yield, and the bomb peak correction (see Sect. 2.3).The error bars for EBC TR / EBC denote the standard deviation (σ ) of the Gaussian fit of EBC TR / EBC as obtained in Fig. 3b.purpose, Eq. ( 13) was solved for α WB : This makes it possible to analytically calculate α WB for every single data point, if a fixed α TR is assumed and setting MAC TR to MAC WB to unity.α WB values were calculated for three different α TR values of 0.9, 1.0, and 1.1, which represent the range previously used in the literature.The resulting three α WB distributions are displayed in Fig. 4. It is evident that an increase in α TR by 0.1 results in a concurrent decrease in α WB by 0.1.This covariance between α TR and α WB implies that using combinations of α TR and α WB randomly altered (e.g., ±0.1) from the best α pair could result in high EBC TR / EBC.No clear station-to-station or seasonto-season differences in α WB were observed (see Table 2), though the number of samples from each station inspected here is limited for such analysis.
Investigating the different distributions in Fig. 4 only the range of α combinations resulting in the best agreement between the Aethalometer model and 14 C results of all individual data points can be obtained but it is not possible to determine other single α pairs representative for the entire data set.To investigate the bias in EBC TR / EBC with respect to EC F / EC ( EBC TR / EBC) due to deviations of α TR and α WB ( α TR and α WB , respectively) from the best α pair, Eq. ( 13) was differentiated with respect to both parameters as a function of EBC TR / EBC.From Fig. S1 (in the Supplement) it is evident that EBC TR / EBC is dependent on EBC TR / EBC: for high and low values of the latter, EBC TR / EBC is mainly driven by α TR and α WB , respectively.A α WB of 0.1 yields a max.EBC TR / EBC of 17 % and a α WB of 0.2 already results in a max.
EBC TR / EBC of 33 %.On the other hand, a α TR of 0.2 results in only a max.EBC TR / EBC of 19 %.Exploring EBC TR / EBC for different α combinations (α TR = 0.9-1.1 and α WB = 1.4-2.2) as a function of EC F / EC (see Fig. 5) it is evident that other α pairs exist yielding low EBC TR / EBC but, in contrast to the best α pair (α TR = 0.9 and α WB = 1.68) not over the entire range of EC F / EC found in this study.Especially for EC F / EC < 30 % almost all α combinations, except the best pair, lead to a significant overor underestimation of EBC TR / EBC compared to EC F / EC.Considering the 1σ confidence interval of EBC TR / EBC (minimum of −0.6 % and maximum of 14 %) as acceptable deviation from EC F / EC also combinations of α TR = 0.8 (see Fig. S2) and 1.0 and α WB = 1.8 and 1.6, respectively, are possible but only for a range of EC F / EC between ∼ 40 and ∼ 85 %.The α pair obtained by Sandradewi et al. (2008a) (α TR = 1.1 and α WB = 1.86) who first used the Aethalometer model results in a constant positive bias of EBC TR / EBC compared to EC F / EC and does not even fall within the 3σ confidence interval of EBC TR / EBC (upper range ∼ 30 %).Furthermore, for α TR of 1.1 only a very narrow range of EC F / EC (spanning maximum 20 %) exists resulting in EBC TR / EBC within the 1σ confidence interval.In addition, from Fig. 5 it is also evident that almost no α WB previously used in the literature (1.8-2.2) would yield EBC TR / EBC within the 1σ confidence for our data set indicating that lower values of α WB should be used in the future in the Aethalometer model.

Evaluation of the Aethalometer model against external data
A further evaluation of the source apportionment results of the Aethalometer model was carried out by comparing EBC WB and EBC TR calculated with the best α TR and α WB pair (0.90 and 1.68, respectively) with other markers for traffic and biomass-burning emissions.Figure 6a presents the correlation of NO x , considered to be from traffic emissions, with EBC TR , both averaged to 24 h from the NABEL stations PAY, MAG, and ZUR for the winter seasons 2009-2012, where Aethalometer and NO x measurements were performed continuously for several years (see Sect. 2.5 and Table 1).Good correlations are found (r = 0.76-0.83)and all stations exhibit similar slopes (24.7-30.7 ppb µg −1 m 3 ) and small axis intercepts (Fig. 6a).These slopes are comparable to London (18-28 ppb µg −1 m 3 , Liu et al., 2014), Grenoble (33 ppb µg −1 m 3 , Favez et al., 2010), and several other locations in Switzerland (32 ppb µg −1 m 3 , Zotter et al., 2014).Levoglucosan obtained on filter samples collected during the winter season and EBC WB were also found to be well corre- lated (r = 0.77; see Fig. 6b) with also only a minor intercept.The slope obtained here (1.08) is also similar to other locations (1.0 for several other locations in Switzerland, Zotter et al., 2014;0.76for three sites in Austria, Caseiro et al., 2009;1.12 in the Po-Valley, Gilardoni et al., 2009;Piazzalunga et al., 2011;and 1.68 in Grenoble, Favez et al., 2010).

Comparison of α TR and α WB with literature values
The α TR value obtained here (0.9) is lower than the values used in Sandradewi et al. (2008a) and many other studies (1.0-1.1,Favez et al., 2010;Crippa et al., 2013;Mohr et al., 2013).However, our findings are in agreement with those reported in Herich et al. (2011) showing that ambient α values in ZUR, MAG, and PAY in summer, when hardly any biomass-burning influence is expected, are around ∼ 0.9.Herich et al. (2011) consequently used then this value as α TR in the Aethalometer model.Also Fuller et al. (2014) determined a value below 1 (α TR of 0.96) for London.
The α WB values obtained in this study are consistent with those reported from smog chamber experiments for fresh and aged biomass-burning emissions (1.63 ± 0.32, Saleh et al., 2013), but are significantly lower than the values from Sandradewi et al. (2008a) often used by other source appor-P.Zotter et al.: Evaluation of the absorption Ångström exponents tionment studies, i.e., 1.8-2.2(Sandradewi et al., 2008a, b;Favez et al., 2010;Perron et al., 2010;Herich et al., 2011;Harrison et al., 2012Harrison et al., , 2013;;Crippa et al., 2013;Mohr et al., 2013).Note that Sandradewi et al. (2008a) determined their best pair of α values (α TR = 1.1 and α WB = 1.86 calculated with b abs,470 , b abs,950 and MAC TR / MAC WB = 1) by optimizing the ratio of the total fossil carbonaceous matter (CM F / CM instead of EC TR / EC) obtained from the 14 C measurements (see Sect. 3.4).Furthermore, Sandradewi et al. (2008a) did not account for the slight underestimation of biomass-burning EC as discussed in Zhang et al. (2012).Applying the approach presented in Sect.3.2.2 to the data in Sandradewi et al. (2008a), using their value for α TR and MAC TR / MAC WB , yields α WB of 1.64 and 1.72 with and without extrapolation of EC F / EC to 100 % EC yield, respectively, which is very similar to the values obtained in this study.Meanwhile, applying the pair of α values previously used from Sandradewi et al. (2008a) to determine EBC TR / EBC from our data results in a mean positive bias ( EBC TR / EBC = 18 %), and therefore the use of this pair is not recommended in future studies.
Recently, Garg et al. (2016) investigated ambient α values in India for various biomass-combustion plumes including paddy-and wheat-residue burning, leaf litter, and garbage burning as well as traffic plumes.They found α values down to 1 for flaming biomass-combustion and α > 1 for older vehicles operating with poorly optimized engines and that α was mostly determined by the combustion efficiency.Therefore, if more than two tightly regulated BC sources, with well constrained combustion efficiencies are present, the α values might be different and additional evaluations of the choice of α WB and α TR in the Aethalometer model using a reference method are needed as well.

Sensitivity due to different wavelength combinations
As different pairs of N-UV and N-IR wavelengths (470 and 950, 470 and 880, and 370 and 880 nm; see Perron et al., 2010;Herich et al., 2011 andFuller et al., 2014) have been previously used in literature, we investigated the sensitivity of the Aethalometer model using different wavelength combinations by performing the same analysis as presented in Sect.3.2.1 with different N-UV and N-IR pairs.In contrast to 470 and 950 nm, no physically meaningful values for α TR could be obtained for the other combinations by fitting Eq. ( 13) against EC F / EC (see Table 3).Consequently, α TR was set to 0.9 to infer α WB for the combinations 470 and 880, 370 and 950, and 370 and 880 nm.As shown in Table 3 different α WB values for these wavelength pairs were obtained than for 470 and 950 nm.Especially using 370 nm as the N-UV wavelength yielded a significantly higher α WB (2.09) than using 470 nm (1.68 and 1.75 for 470 and 950 and 470 and 880 nm, respectively).It has been reported that α is wavelength dependent (e.g., Bond and Bergstrom, 2006), and might be more affected by fuel type, combustion, and atmospheric processes in the N-UV than in the visible part of the spectrum (Sandradewi et al., 2008c).Consequently, α WB can be different for different wavelength pairs.However, for all combinations, especially with 370 nm as N-UV wavelength, the mean residuals of EBC TR / EBC compared to EC F / EC, were higher than using the 470 and 950 nm combination (see Table 3a).Next, EBC TR / EBC was calculated with the best pair of α values (α WB = 1.68 and α TR = 0.90 as obtained in Sect.3.2.1)for the different wavelength combinations.It is evident that using 370 nm as N-UV wavelength EBC TR / EBC exhibits an inferior agreement with EC F / EC (see Table 3b).EBC TR / EBC exhibits larger values, there is a significant number of negative points for EBC TR /EBC and the correlations with EC F / EC are weaker.
On the other hand, similar EBC TR / EBC and hardly any negative EBC TR / EBC values are found for the wavelength combination 470 and 880 nm.Uncertainties in the source apportionment results using the Aethalometer model due to the use of different wavelength pairs are usually not considered and often the same α WB and α TR values are used with different wavelength combinations.However, as shown here, the choice of the wavelengths, especially the one in the N-UV range, and α WB are not independent.Since (1) it was previously shown that adsorption of volatile organic compounds on the filter tape of the Aethalometer can occur which possibly influences the absorption measurement with the 370 nm channel (Vecchi et al., 2014), (2) light-absorbing SOA, other absorbing non-BC combustion particles, and atmospheric processing affect lower wavelengths more than higher ones, and (3) our results indicate an inferior agreement of EBC TR / EBC with EC F / EC using 370 nm as N-UV wavelength, we therefore recommend using 470 nm as the N-UV wavelength in the Aethalometer model while the choice between 950 and 880 nm in the N-IR is less critical.

High time resolution data
Since Aethalometers measure with high time resolutions (e.g., model AE31 down to 2 min and the new model AE33 down to 1 s) the investigation of the temporal behavior of EBC TR and EBC WB is possible (see, e.g., Herich et al., 2011).Figures 7 and S3 show the diurnal cycles for the stations MAG, PAY, and ZUR including continuous data from the entire years 2009 to 2012.It is evident that the Aethalometer can also be applied to high time resolution data and the expected temporal behavior of the sources can be resolved.The contribution of EBC WB is high in winter and during the night, with only small differences between weekends and weekdays.Furthermore, EBC TR exhibits a clear traffic peak in the morning during weekdays whereas during weekends this increase is not evident or only small.
We note that EBC WB also follows EBC TR , with an evident increase during morning hours.This increase is statistically Table 3. Evaluation of the Aethalometer model using different wavelength pairs.(a) Calculation of the α values by fitting Eq. ( 13) (MAC TR / MAC WB = 1) against EC F / EC.(b) comparison between EC F / EC and EBC TR / EBC calculated with α WB = 1.68 and α WB = 0.90, representing the best α pair for all data, for different wavelength pairs.µ and σ denote the center and the standard deviation of the Gaussian fit of the residuals of EBC TR / EBC compared to EC F / EC, respectively.

(a)
Calculation  13) against EC F / EC and therefore α TR was set to 0.9 representing the best value for the wavelength pair 470 and 950 nm.larger than our uncertainties (14-18 in winter and 30-75 % in summer).This indicates that there is some false attribution of EBC TR and EBC WB most probably due to the constant a priori assumed pair of α and α TR .By applying different α combinations for ZUR (α WB = 1.68 and α TR = 0.90, α WB = 1.68 and α TR = 1.1, α WB = 1.9 and α TR = 0.90 as well as α WB = 1.9 and α TR = 1.1) this false attribution between EBC TR and EBC WB during the morning peak disappears (see Fig. S4), indicating that a higher α TR would be more representative of fresh traffic emissions in the case of ZUR.Since the evaluation of α combinations presented in this paper is based on longer sampling times and mostly winter data (16 to 40 h; see Table 1), caution should be taken when applying the Aethalometer model with the best α pair found here to high time resolution data, especially for low EBC concentrations and rush hours.Similar studies with higher time resolutions, and for EBC concentrations, like in summer, are necessary for a further evaluation of the Aethalometer model.

Traffic and wood-burning contributions to EC and EBC
The relative traffic contribution as apportioned by the Aethalometer model (EBC TR / EBC) and the 14 C analysis (EC F / EC) of EBC and EC, respectively, is often > 50 % (see Fig. 3a).However, since hardly any EC F / EC values, except results from the summer season, are above 70 % and the average of EC F / EC over all winter samples is 52 ± 17 %, it is evident that also wood-burning emissions account for a large fraction of EC (and thus EBC) during winter in Switzer-land.The traffic contributions for winter samples range from 7 to 82 and 14 to 84 % for EBC and EC, respectively.The lowest values (31 ± 23 and 36 ± 17 % for EBC TR / EBC and EC F / EC, respectively) were found at ROV which is most likely due to a combination of topography (ROV is located in an Alpine valley), local meteorology (often persistent inversions with low mixing heights) and emissions (high local wood-burning influence; Alfarra et al., 2007;Szidat et al., 2007;Lanz et al., 2008;Sandradewi et al., 2008a, c;Herich et al., 2014;Zotter et al., 2014).The samples from ZUR, the largest city of Switzerland, collected during spring and summer clearly show the highest fossil contributions with an average of 81 ± 10 and 80 ± 7 % and the highest value of 92 and 85 % for EBC TR / EBC and EC F / EC, respectively.Investigating the diurnal cycles of EBC WB , it is evident that the concentrations are high in winter, especially in MAG and during night-time, with no or only small differences between weekends and weekdays.Concentrations in summer are lower but non-negligible, with significantly (paired t test, significance level of 0.05) higher concentrations in Zurich compared to the other locations, especially during night-times on weekends.This suggests an additional source of brown carbon in ZUR, likely related to primary emissions enhanced with urban activities, during weekends, and with higher emissions in an increasingly shallower nighttime boundary layer.Contribution of secondary processes to the brown carbon background concentrations observed at all sites cannot be excluded.EBC TR concentrations for weekdays are significantly higher in winter compared to summer and also for weekends in MAG (t test, significance level of  13) against EC F / EC with the additional traffic (NO x ) and wood-burning (levoglucosan) markers: (a) the correlation between EBC TR at 950 nm and NO x averaged to 24 h; and (b) the scatterplot between the EBC WB at 950 nm and levoglucosan.Details of the origin of the data are given in Table 1.0.05).In contrast, in ZUR and PAY, average EBC TR weekend concentrations are very similar.Weekday EBC TR concentrations exhibit a clear morning peak for all stations and seasons, which is less evident on weekends.The lowest concentrations are found in PAY and the highest in MAG in winter.In summer, EBC TR is highest in ZUR days.

Traffic and wood-burning contributions to PM
It has been attempted to also apportion the total carbonaceous material (CM) to wood-burning (CM WB ) and traffic (CM TR ; e.g., Sandradewi et al., 2008a)  In practice, site-specific c 1 , c 2 , and c 3 values may be either fitted using Eq. ( 17) or less commonly fixed based on the knowledge of the OM-to-EC ratios in the primary emissions of interest (most frequently, only c 1 is fixed, e.g., Favez et al., 2010).The two approaches do not necessarily lead to the same result as they are not based on the same conceptual definitions of the organic aerosol fractions.When derived from the multiple linear regression fitting of Eq. ( 17), CM TR and CM WB would represent the fractions that correlate with b abs,TR,950 and b abs,WB,470 , respectively.As SOA production is often very rapid (Huang et al., 2014), CM TR and CM WB are also expected to partially contain not only primary OA, but also SOA produced through the aging of traffic and wood-burning emissions, respectively.Note that correlation is not causation and some of these correlations are a direct consequence of meteorology and planetary boundary layer mixing.This may significantly complicate data interpretation.Using this methodology, Herich et al. (2014) could not precisely quantify the contributions of the different CM sources.They found a standard error of ± 30 % for c 1 , c 2 , and c 3 and a high sensitivity of c 1 and c 2 on the chosen α values for wood burning and traffic emissions.This is one of the few cases where errors related to CM apportionment using Aethalometer data are explicitly estimated.Usually only the sensitivity of c 1 , c 2 , and c 3 on the chosen α values is reported.However, the most profound flaw in the application of Eq. ( 17) using a multiple linear regression is the assumption that c 3 , which is non-absorbing SOA mass mostly, is constant over time.Consequently, assessing how this mass is apportioned among the different sources and model residuals remains elusive and more faithful representation of the complex atmospheric processes would necessitate the use of a robust tracer for SOA.Accordingly, we do not recommend the use of this model in its current state to apportion the CM mass, especially when the SOA fraction is dominant and highly variable.
A more conservative and controlled approach is to fix in the model the values of c 1 and c 2 , based on emission data, and attribute the time-dependent remainder (c 3 ) to SOA.Under these conditions, CM TR and CM WB relate to the primary fraction and the model may better capture the time variability of SOA.Indeed, this approach would entail the precise knowledge of c 1 and c 2 .OM-to-EC ratios in traffic emissions are heavily dependent on the type of fuel used, with lower values reported for diesel exhausts.Accordingly, for a European fleet dominated by diesel cars, El Haddad et al. (2013) report OM-to-EC ratios ranging between 0.25-0.45,whereas in the US, ratios between 0.9 and 1.4 were found (Zhang et al., 2005;Sun et al., 2012;Stroud et al., 2012).As biomass burning is a poorly controlled combustion process, typical OM-to-EC ratios determined at emissions are highly scattered, ranging between 3 and 63 (Schauer et al., 2001;Fine et al., 2001Fine et al., , 2002Fine et al., , 2004a, b;, b;Schmidl et al., 2008).Nevertheless, more useful information may be obtained from examining ambient measurements, where, depending on the approach used to quantify wood smoke, OM-to-EC ratios may range between 3 and 18 (Favez et al., 2010;Herich et al., 2014;Zotter et al., 2014).While these ratios must be selected with extreme caution and the sensitivity of the source apportionment results to this selection must be systematically assessed, additional on-site data (e.g., levoglucosan, 14 C. . . ) may always aid constraining their values.

Conclusions
In this study, we show a comparison of the source apportionment of equivalent black carbon (EBC) using the Aethalometer model with radiocarbon ( 14 C) measurements of elemental carbon (EC).This enables a validation of the choice of the absorption Ångström exponents for wood-burning (α WB ) and traffic (α TR ) emissions which have to be assumed a priori in the Aethalometer model.Data from several campaigns across Switzerland with parallel Aethalometer and 14 C measurements of the EC fraction from eight different stations with different characteristics allow the investigation of the applicability and performance of the Aethalometer model for different locations and conditions.
To obtain the best α pair in the Aethalometer model, outputs (using the 470 and 950 nm channels) were fitted against the fossil fraction of EC (EC F / EC) derived from 14 C measurements resulting in α TR = 0.9 and α WB = 1.68.The source apportionment results from both methods, Aethalometer and 14 C, are well correlated (r = 0.81) and the fitting residuals exhibit only a minor positive bias of 1.6 % and an average precision of 9.3 %, indicating that the Aethalometer model reproduces reasonably well the 14 C results for all stations investigated in this study using our best estimate of a single α WB and α TR pair.We show that the Aethalometer model also works for other areas than for polluted Alpine valleys in winter, since this analysis includes data from urban stations as well as days from spring and summer.
The residuals of the Aethalometer model outputs ( EBC TR / EBC) calculated with other α pairs depend on EC F / EC and a good agreement (within the 1σ confidence interval of EBC TR / EBC) cannot be obtained over the entire EC F / EC range using other α pairs.However, combinations of α TR = 0.8 and 1.0 and α WB = 1.8 and 1.6, respectively, are also possible but only for a range of EC F / EC between ∼ 40 and ∼ 85 %.The α WB values previously used in Aethalometer models (∼ 2) and any combination with α TR = 1.1 yield significant positive biases in the fitting residuals.Therefore we recommend to use the best α combination as obtained here (α TR = 0.9 and α WB = 1.68 for the wavelength pair 470 and 950 nm) in future studies.We also tested the sensitivity of the Aethalometer model due to different pairs of near-UV (N-UV) and near-IR (N-IR) wavelengths (470 and 950, 470 and 880, and 370 and 880 nm).Any combination with 370 nm as N-UV wavelength resulted in larger residuals, a significant number of negative points, and weaker correlations with EC F / EC and, therefore, we recommend to use 470 nm as N-UV wavelength in the Aethalometer model.Using 950 or 880 nm as N-IR wavelengths showed similar results, though the former wavelength performed slightly better in this study.
Having an independent measurement for the relative contributions of traffic and wood-burning to total EC from the 14 C and Aethalometer measurements, also made it possible to assess the dependence of the mass absorption cross section (MAC) of EBC on its source.The results indicate no significant difference in MAC at 880 nm (with a value of 11.8 m 2 g −1 ) between EBC originating from traffic or woodburning emissions.

P. Zotter et al.: Evaluation of the absorption Ångström exponents
Applying the Aethalometer model to apportion total carbonaceous material (CM) it is evident that there are significant uncertainties and model errors (mainly due to assuming constant fitting parameters relating the absorption of traffic, wood-burning and the residuals (comprising non-lightabsorbing CM, mostly secondary organic aerosol) to the separately determined CM).Therefore, in our opinion such a CM apportionment should only be interpreted qualitatively.
The results obtained in this study demonstrate that the evaluation of the choice of α WB and α TR using a reference method is highly valuable and should be performed when applying the Aethalometer model, if possible.In the absence of such reference measurements, however, assuming a single set of α TR and α WB yields acceptable results (i.e., average precision of 9.3 % of EBC TR / EBC compared to EC F / EC in our case) and provides the best estimate of the fossil and non-fossil contributions to EBC as apportioned by the Aethalometer model.Nevertheless, these results were obtained for locations impacted by BC mainly from traffic consisting of a modern car fleet and wood combustion for residential heating in winter with well-constrained combustion efficiencies.Furthermore, mainly winter conditions with only a few summer samples were available.Therefore, additional studies about the performance of the Aethalometer model with respect to seasonality and for sites with different combustion conditions and efficiencies, sources, and fuels used and their temporal evolution are needed to reduce the uncertainties of their choice in the Aethalometer model.

Figure 1 .
Figure 1.Location of the different stations in Switzerland investigated in this study.

Figure 2 .
Figure 2. Distribution of MAC values of EBC at 880 nm (a) and comparison with the fossil fraction of EC (EC F / EC) determined with the 14 C analysis (b).MAC values were determined assuming a C value of 2.14 for the Aethalometer and the EUSAAR-2 thermal optical transmission protocol was used for EC mass measurement.Only data from the 14 C project Switzerland are included, since in earlier campaigns EC concentrations were not determined or measured with the same TOT protocol.

Figure 3 .
Figure 3. (a) Comparison between EBC TR / EBC at 950 nm and EC F / EC and (b) residuals of EBC TR / EBC compared to EC F / EC ( EBC TR / EBC).EBC TR / EBC was calculated using b abs,470 , b abs,950 , MAC TR / MAC WB = 1, and the α values (α WB = 1.68 and α TR = 0.90) obtained by fitting Eq. (13) against EC F / EC.The error bars for EC F / EC represent the combined uncertainty of the 14 C measurement of EC, the extrapolation of f M,EC to 100 % EC yield, and the bomb peak correction (see Sect. 2.3).The error bars for EBC TR / EBC denote the standard deviation (σ ) of the Gaussian fit of EBC TR / EBC as obtained in Fig. 3b.

Figure 5 .
Figure 5. Residuals of EBC TR / EBC compared to EC F / EC ( EBC TR / EBC) as a function of EC F / EC for different combinations of α TR and α WB .Average EBC TR / EBC values for EC F / EC bins of 0.1 are calculated for α WB = 1.4-2.2 and α TR = 0.9 (upper right panel), α TR = 1.0 (lower left panel) and α TR = 1.1 (lower right panel).EBC TR / EBC for α TR = 0.8 can be found in Fig. S2.The number of points per EC F / EC bin is displayed in the upper left panel.The dashed grey line denotes the best α pair (α TR = 0.9 and α WB = 1.68) as obtained in Sect.3.2.1 and the dark and light grey shaded areas mark the 1σ (standard deviation) and 3σ of EBC TR / EBC per EC F / EC bin for this best α pair.The black dashed line in the lower right panel represents the α combination obtained by Sandradewi et al. (2008a) who first used the Aethalometer model.

Figure 6 .
Figure 6.Comparison of the Aethalometer model outputs calculated with the α values (α WB = 1.68 and α TR =0.90) obtained by fitting Eq. (13) against EC F / EC with the additional traffic (NO x ) and wood-burning (levoglucosan) markers: (a) the correlation between EBC TR at 950 nm and NO x averaged to 24 h; and (b) the scatterplot between the EBC WB at 950 nm and levoglucosan.Details of the origin of the data are given in Table1.
according to the following equations: CM = OM + BC, OM. ..organic matter (16) CM = CM TR + CM WB + CM other = c 1 × b abs,TR,950 + c 2 × b abs,WB,470 + c 3 (17) If CM is determined independently, c 1 and c 2 can be obtained by solving Eq. (17), relating the light absorption to the particulate mass of both sources.The intercept c 3 represents a con-stant background concentration of non-absorbing carbonaceous material (CM other ).While Sandradewi et al. (2008a) did not require CM other to achieve mass closure, Favez et al. (2010), Harrison et al. (2013), and Herich et al. (2011) found significant contributions of CM other .

Table 1 .
List of all stations and their classification according to the Swiss Federal Office for the Environment (FOEN) and additional campaign details.

Table 2 .
Ranges and averages of α WB values resulting in a ratio of 1 between EBC TR / EBC (at 950 nm) and EC F / EC for all stations calculated with α TR of 0.9, b abs,470 , b abs,950 , and MAC TR / MAC WB = 1.
of best α values EC vs. EBC TR / EBC with α TR = 0.90 and α WB = 1.68No physically meaningful value for α TR could be obtained by the fitting of Eq. ( *