Evaluation of equivalent black carbon source apportionment using observations from Switzerland between 2008 and 2018

. Black carbon (BC) or soot is a constituent of particulate matter (PM) which is relevant for negative human health and climate effects, and despite the lack of direct legal limits, it is recognised as an important atmospheric pollutant to monitor, understand, and control. Aethalometers are instruments which continuously monitor BC by measuring absorption at a number of distinct wavelengths. If collocated elemental carbon (EC) observations are used to transform these values into BC mass, by convention, the result is named equivalent black carbon (EBC). BC emitted by different combustion processes has different optical absorption characteristics, and this can be used to apportion EBC mass into trafﬁc (EBC TR ) and woodburning (EBC WB ) components with a data processing technique known as the aethalometer model. The aethalometer model was applied to six EBC monitoring sites across Switzerland (using data between 2008 and 2018) and was evaluated by investigating diurnal cycles, model coefﬁcients, and ambient temperature dependence of the two EBC components. For one monitoring

Abstract. Black carbon (BC) or soot is a constituent of particulate matter (PM) which is relevant for negative human health and climate effects, and despite the lack of direct legal limits, it is recognised as an important atmospheric pollutant to monitor, understand, and control. Aethalometers are instruments which continuously monitor BC by measuring absorption at a number of distinct wavelengths. If collocated elemental carbon (EC) observations are used to transform these values into BC mass, by convention, the result is named equivalent black carbon (EBC). BC emitted by different combustion processes has different optical absorption characteristics, and this can be used to apportion EBC mass into traffic (EBC TR ) and woodburning (EBC WB ) components with a data processing technique known as the aethalometer model. The aethalometer model was applied to six EBC monitoring sites across Switzerland (using data between 2008 and 2018) and was evaluated by investigating diurnal cycles, model coefficients, and ambient temperature dependence of the two EBC components. For one monitoring site, San Vittore, the aethalometer model failed to produce plausible outputs. The reason for this failure was likely due to a high load of freshly emitted wood smoke during the winter which should be thought of as a third distinct emission source. After model evaluation, the trend analysis indicated that EBC TR concentrations at the remaining five locations significantly decreased between 2008 and 2018. EBC WB also demonstrated significant decreases in most monitoring locations but not at a monitoring site south of the Alps with a high PM load sourced from biomass burning. Ratios of EBC and particulate matter with a diameter of less than 2.5 µm (PM 2.5 ) suggested that EBC contributes 6 %-14 % of the PM 2.5 mass in Switzerland. The aethalometer model is a useful data analysis procedure but can fail under certain conditions; thus, careful evaluation is required to ensure the method is robust and suitable in other locations.

Black carbon
Atmospheric particulate matter (PM) has a variety of components, one of which is black carbon (BC), which is commonly referred to as soot. BC is strongly light absorbing (hence the name), is generally found in the fine PM fraction (PM with a diameter less than 2.5 µm; PM 2.5 ) and is generated by the incomplete combustion of fuels (Hansen et al., 1984;Vignati et al., 2010). Biomass burning, internal combustion engines (especially those which are fuelled by diesel), and industrial processes can all be BC emission sources (Bond et al., 2004;Jacobson, 2001). There are no known BC generation or degradation processes in the atmosphere, and because the only removal mechanisms for BC are wet and dry deposition, BC has a lifetime of several days to weeks in the atmosphere. These attributes make BC a reliable tracer for combustion processes.
BC is an important atmospheric pollutant to consider in its own sense however (Anenberg et al., 2012). BC is believed to have a positive (warming) climate forcing effect due to the absorption of radiation and the reduction of albedo, especially when deposited on snow and ice (Ramanathan and Carmichael, 2008;Weinhold, 2012;Bond et al., 2013). BC is harmful to human health, and there is evidence that BC is a particularly potent PM component when considering deleterious health effects due to the aerosol's large surfacearea-to-mass ratio, tendency to be enriched with harmful organic compounds, and ability to penetrate deep into the lungs (Janssen et al., 2011;De Prins et al., 2014;Laeremans et al., 2018). The International Agency for Research on Cancer (IARC) has also classed soot as a Group 1 carcinogen, indicating there is clear evidence that soot causes cancer in humans (International Agency for Research on Cancer, 2019). It is difficult however to tease apart the health effects of BC alone from PM as a whole (Jacobson, 2001;European Environment Agency, 2013). These features have made BC an important pollutant to monitor, understand, manage, and control despite not currently having direct legal limits imposed in Europe and elsewhere (Reche et al., 2011;European Environment Agency, 2016).

Aethalometers and the aethalometer model
BC is measured, sampled, or monitored in a variety of ways, but the most widespread method is the use of aethalometers (Hansen et al., 1984). Aethalometers are described in depth elsewhere; however, briefly they are optical instruments which sample air continuously and deposit the PM onto filter material. The loaded filter is illuminated; the light attenuation is measured as an optical absorption; and the increase in attenuation over time is logged and used to calculate BC concentrations with empirical coefficients. The quartz filter material used by aethalometers can reach saturation, so the instruments employ a filter-tape system where the filter material is advanced if saturation occurs or when a threshold of elapsed time has passed. Continuously depositing aerosol onto a filter causes artefacts known as filter-scattering and filter-shadowing effects, but there are standard algorithms to compensate for these effects, some of which are applied on board by the instrument as part of the measurement cycle (Magee Scientific, 2016).
Modern aethalometers used in ambient air quality monitoring networks are multi-wavelength models measuring absorption between the near-ultraviolet (UV) and the nearinfrared (IR) range, at distinct wavelengths within this range. BC sourced from different families of combustion processes demonstrate different absorption features. Most notably, woodburning particles, i.e. wood smoke, generally contain a rich organic component which is very effective at absorbing light in the UV range, while diesel soot shows a weaker spectral dependence (Kirchstetter et al., 2004) (Fig. 1).
The spectral dependence of the light absorption can be described by b abs ∝ λ −α , where b abs is the absorption, λ is the wavelength, and α is the Ångström exponent (Moosmüller et al., 2011;Harrison et al., 2013). The α value is the exponential slope of a regression model for absorption as a function of wavelength; α for woodburning-sourced BC results in higher values due to the increased absorption in the UV Figure 1. Demonstration of different aethalometer absorption dependence on wavelengths for two monitoring sites in Switzerland with distinct aerosol characteristics. Data have been filtered to a single observation (hourly) to show the dependence on the dominating sources, and the 470, 880, and 950 nm channels are highlighted and explained in text. Magadino-Cadenazzo has a stronger absorption in the UV region due to woodburning emissions. region when compared to BC emitted by vehicular sources (Fig. 1).
The different spectral dependences between these sources can be leveraged with a data processing technique called the aethalometer model, first reported by Sandradewi et al. (2008a) but used many times since in many locations; for examples, see Sandradewi et al. (2008b), Herich et al. (2011), and Fuller et al. (2014. In routine aethalometer use, only the 880 nm absorption wavelength is used with a fixed coefficient (called a mass absorption coefficient discussed in Sect. 2.2), which is the strict definition of BC. The aethalometer model uses UV and IR measurements to represent the two sources, and the wavelengths used are usually 470 and 950 nm for the UV and IR channels respectively (also shown in Fig. 1; Zotter et al., 2017). UV and IR α values are also required for the source apportionment and need to be carefully chosen on a case-by-case basis however. The aethalometer model simply uses information which is generally discarded in most monitoring applications, and the barrier for entry to the method is low.

Black carbon nomenclature note
The nomenclature and vocabulary used for BC is diverse and can lead to confusion (Andreae and Gelencsér, 2006;Petzold et al., 2013). Here, we refer to BC determined by opti-cal measurement techniques, i.e. by aethalometers and transformed with the use of collocated elemental carbon (EC) observations as equivalent black carbon (EBC). This definition has been recommended by the Global Atmospheric Watch (Tarasova, 2012). When EBC is apportioned into its traffic and woodburning components, the TR and WB subscript notation is used respectively, i.e. EBC TR and EBC WB . If total EBC needs to be used for clarity against the EBC components, EBC TOT is used. This notation is consistent with Zotter et al. (2017), which this work somewhat extends. When discussing emissions of soot, BC is still used however because EBC is only an operational definition and does not refer to the pollutant.

Objectives
The primary objective of this work is to apply the aethalometer model data transformation technique to apportion longterm (up to 10 years) EBC time series from six Swiss monitoring sites into two components -a traffic and a woodburning component (EBC TR and EBC WB ) -and evaluate the method's performance. The suitability and limitations of the aethalometer model as a technique to employ before trend analysis will be discussed. A data set is presented where the aethalometer model fails to apportion the two EBC components correctly. This is done to demonstrate the features which will be present if the aethalometer model is inadequate for a particular application, and it outlines that although the aethalometer model is a useful technique, it is not a panacea. The second objective is to expose the EBC TR and EBC WB components to a trend analysis and document the features observed in the Swiss time series.

Data
Absorption observations measured by aethalometers for six monitoring sites in Switzerland were analysed (Table 1; Fig. 2). The monitoring sites were either part of the federal monitoring network (National Air Pollution Monitoring Network; NABEL) or networks run by the Swiss Cantonal authorities (states; Federal Office for the Environment, 2014). Sites were classified as either urban traffic, urban background, rural, or rural mountain according to their surrounds. Four of the sites were located on the Swiss plateau (where the majority of the human population is located), but Rigi-Seebodenalp, despite being located on the plateau, is at altitude (1031 m; Fig. 2). Two monitoring sites, Magadino-Cadenazzo and San Vittore, are located south of the Alps in valleys where residential woodburning is much more common than on the Swiss plateau due to its rural nature Szidat et al., 2007).
The absorption observations were measured by two models of multi-wavelength aethalometers manufactured by Magee Scientific. The two models were the AE31 and AE33, which measure adsorption at seven wavelengths: 370, 470, 520, 590, 660, 880, and 950 nm. The AE31 is an earlier generation instrument discontinued in 2016 (Magee Scientific, 2017), while the AE33 aethalometer replaced the AE31 and is currently manufactured (Magee Scientific, 2019). These two aethalometer models share the same measurement principle; however the newer AE33 uses a "dual-spot" approach to allow for a superior method to compensate for filterloading effects (Drinovec et al., 2015). The differences between the AE31 and AE33 technologies mean that the algorithms which compensate or correct for filter-shadowing effects and filter-loading effects are different. Notably, data from AE31 must be compensated for such effects with a postprocessing procedure, while the AE33's algorithms are conducted on board as part of the measurement cycle. For the full description of the compensation procedures, see Weingartner et al. (2003) and Drinovec et al. (2015). Quartz filters were used by the aethalometers, and therefore, a multiple light-scattering factor (denoted as C) of 2.14 was used (Weingartner et al., 2003). There were 13 distinct aethalometers used in the monitoring network which were irregularly cycled among the monitoring sites as part of the networks' maintenance procedures (Table 2). All aethalometers were operated with PM 2.5 sample inlets.
Data from the aethalometers for five of the six Swiss EBC monitoring sites were queried from the NABEL monitoring network's database, while the remaining site's observations, at San Vittore, were provided directly by the Graubünden canton's environmental department (Amt für Natur und Umwelt). All absorption observations had been compensated for the filter-loading and shadowing effects with the instrument model's respective algorithms before this analysis was undertaken as part of the monitoring networks' operations (Weingartner et al., 2003;Drinovec et al., 2015). Generally, the observations were stored as hourly means, but for the data which were at a higher resolution (10 and 30 min means), observations were aggregated to create a consistent hourly time series which was used for analysis.
The standard procedure for absorption monitoring data was for 10 and 30 min means to be calculated from 1 min observations which were logged on site directly from the instruments. All aggregations required data capture thresholds of 60 % for a valid summary to be calculated. Additional quality control and assurance procedures were undertaken quarterly including the cleaning of the inlet, leak testing, and cleaning of the analytical zone of the instruments. The responses of the AE33 instruments were also checked with optical reference filters regularly. The data were ratified on a monthly basis and compared across the different measurement sites along with other air pollutants, and suspicious measurements were invalidated. The raw instrument outputs are archived but are not routinely used in the data processing within the NABEL monitoring network. However, the raw data can be consulted in case of questions concerning data quality.  The method employed to transform absorption observations to EBC requires elemental carbon (EC) concentrations. EC was determined by the standard EN16909 thermaloptical transmission (TOT) method using the EUSAAR2 temperature protocol (European Committee for Standardization, 2017). Unlike the continuous observations which the aethalometers provided, EC concentrations were available as daily samples and ranged from being sampled every 4 to 12 d. When using aethalometer and EC data together, the aethalometer observations were aggregated (as arithmetic means) to daily resolution (midnight to midnight) to ensure that the observations spanned the same time period and duration. Only days with both EC and absorption observations were used for these comparisons because interpolation of the less frequent EC data was not attempted. The EC samples were also collected with PM 2.5 inlets. PM 2.5 observations for the monitoring sites were required to calcu-late the EBC/PM 2.5 ratios. These data were accessed with the saqgetr R package, which gives access to the European Commission's AirBase and Air Quality e-Reporting (AQER) repositories in a convenient way (European Environment Agency, 2014, 2019; Grange, 2019b). Only daily validated data were kept for analysis with most of the observations being sourced from high-volume samplers.
The analysis presented covered a period between March 2008 and the end of 2018, but the start date of the aethalometer operation varied among the different monitoring sites (Tables 1, 2). All absorption observations, site metadata, and instrument locations data are publicly accessible in a persistent data repository .

Source apportionment
The absorption data from the aethalometers were apportioned into two components: the traffic fraction (EBC TR ) and  Table A1), but any given MAC value is a function of PM size and morphology which will change depending on the PM source (Zotter et al., 2017). For the trend analysis reported here, the MAC values for absorption at the 950 nm wavelength (used for the IR input to the aethalometer model) were calculated with rolling simple least-squares linear regression models between absorption and the mass of EC in the PM 2.5 fraction. A least-squares estimator was used because EC was determined with a reference method and is accepted without uncertainty. The window was 180 d (with centre alignment), and this logic was implemented with the zoo R package (Zeileis and Grothendieck, 2005).
Once the Ångström exponents were chosen and the MAC values were calculated, the aethalometer model was applied to the absorption observations at hourly resolution, and the R function used for this transformation is available (Grange, 2019a). Based on the absorption at 470 and 950 nm (UV and IR spectra respectively), this resulted in three EBC mass variables for each valid absorption observation: EBC WB , EBC TR , and EBC TOT . A flow diagram which represents the data processing steps can be found in Fig. A1.

Trend tests
Formal trend tests were performed on the EBC components with the Theil-Sen slope estimator, a robust, non-parametric estimator provided by the openair R package (Carslaw and Ropkins, 2012). Before the trend was tested, the observations were aggregated to a monthly resolution and deseasonalised with loess models to extract the trend component (R Core Team, 2019). Autocorrelation was taken into account, resulting in conservative slope estimations, and all trend tests were conducted at the 0.05 significance level.
3 Results and discussion 3.1 Model coefficient evaluation

Mass absorption cross-section (MAC) coefficients
The MAC values which were empirically derived using least-squares regression models (using EC and absorption at 950 nm), with a 180 d rolling window showing substantial variability (Fig. 3). This suggests that the use of fixed or static MAC values is questionable for trend analysis applications. The sites' MAC values at 950 nm had an average of 11.3 ± 2.9 m 2 g −1 during the analysis period despite the irregular cycling of different instruments among the sites (Table 2). For perspective, compare this range with the factory MAC values displayed in Table A1. A seasonal pattern was present in the MAC values where summertime absorption measures were higher than those in wintertime, which represents a change of average aerosol characteristics during the year, which is to be expected. The range of the seasonality was as low as 0.4 m 2 g −1 at Zürich-Kaserne but as high as 2.3 m 2 g −1 at Magadino-Cadenazzo (Fig. A2). The seasonal variation was also accompanied by intra-instrumental variation.
The MAC values for Rigi-Seebodenalp were atypical and showed a substantial decrease between 2013 and 2016 (Fig. 3). During this period, the same aethalometer was operating at the monitoring location (Table 2; Fig. A3). The MAC value is an empirical coefficient which represents aerosol composition, and it is very unlikely that Rigi-Seebodenalp's aerosol characteristics altered dramatically between 2013 and 2016; therefore the observed decrease in the MAC value was almost certainly an instrument artefact -perhaps a slow degradation in sensitivity. Indeed, the decrease in the MAC value stopped when the instrument was replaced in late 2015 and remained low for the rest of the monitoring period (Fig. A3). This instrument-derived feature was only uncovered due to the presence of collocated EC observations, but due to a lack of additional data, we do not wish to speculate further on what was the cause for this decline in the MAC coefficient, and this may warrant further investigation elsewhere.
The rolling least-squares regression procedure presented here ensures that the changes in MAC over time are compensated for before EBC mass is calculated with the aethalometer model. The changes or drift in the MAC values, especially for Rigi-Seebodenalp, are not propagated into the EBC masses using this method. Therefore, this approach is recommended over using fixed MAC values if EC time series are available and trend analysis is to be undertaken with several years of aethalometer observations.

Ångström exponents (α)
A critical input for the aethalometer model is the α values used for the traffic-and woodburning-sourced EBC (Harrison et al., 2013). The α values used here were previously derived from comparing 14 C and EC observations (Sect. 2.2; Zotter et al. (2017)), but to validate if the values were sensible for the sites analysed, α was calculated for each hourly absorption observation, and their distributions were investigated. Here, α was calculated by the curve fitting of all absorption wavelengths (370, 470, 520, 590, 660, 880, and 950 nm) with exponential regression models.
The distributions of the calculated α values were consistent with those reported by Zotter et al. (2017), which were 0.9 and 1.68 for α TR and α WB respectively. Fig. 4 shows the α distributions for three of the monitoring sites along a continuum of increasing woodburning activity (from left to right). The distributions for α at each of the monitoring sites peaked at 1.2, the value used for α TR in the aethalometer model. However, as the monitoring sites' were exposed to a progressively higher woodburning PM load, the distribution become positively skewed. Magadino-Cadenazzo for example, a site which is exposed to wood smoke, has a tail which extends beyond 1.68, while the urban-traffic Bern-Bollwerk site's α distribution contains only a small component at, and beyond, the α WB value of 1.68. Therefore, the features of the α distributions can be rather informative for diagnostic purposes.
San Vittore's calculated α values were unique because they displayed a bimodal distribution which was not present at the other monitoring sites (Fig. 4). The presence of the bimodal distribution indicates that San Vittore is exposed to an additional source that the other monitoring sites are not. Although this extra BC source is unknown, we speculate that it is very likely to be freshly emitted wood smoke which  had been emitted near the monitoring location and rapidly transported to the monitoring site before chemical processing had the opportunity to act on the aerosol. The importance of this feature is discussed further in Sect. 3.3. An additional check on the observations was performed in the form of a three-factor receptor model with the multilinear engine (ME-2) as implemented by the EPA PMF (positive-matrixfactorisation) tool (Norris et al., 2014;Brown et al., 2015). However, in this case the factor analysis was unable to resolve the three sources and did not offer any additional con-tribution to the data analysis. Time series of the calculated α values for the six EBC sites are also available in Fig. A4.

Diurnal cycles
Diurnal plots were used to validate if the two components' behaviour were consistent with what is expected of their emission source behaviour after the EBC fractions were calculated with the aethalometer model, and the model's coefficients were evaluated (Fig. 5). With the exception of Rigi-Seebodenalp, an elevated rural-mountain site, both the traffic and woodburning components' diurnal cycles could be readily understood in terms of source activity influencing the monitoring sites. This gives support to the plausibility of the aethalometer model's source apportionment ability.
Bern-Bollwerk, an urban-traffic site displayed a clear morning rush hour peak for EBC TR with a decline in concentrations in the late morning, followed by a smaller increase in the afternoon-evening rush hour (Fig. 5). This traffic-forced pattern was also present at Payerne and Zürich-Kaserne, rural and urban-background sites respectively, but at lower concentrations. Magadino-Cadenazzo also displayed a clear twin-peak EBC TR diurnal cycle. Here however, the evening peak had approximately the same concentration as the morning peak which was not observed at the other monitoring locations. The strong evening peak was only observed in winter at Magadino-Cadenazzo. This feature was most likely driven by very stable and stagnant atmospheric conditions in the morning and evening because this monitoring site (and surrounding population centres) is located in a deep valley where the day length is short at these times of the year due to confining terrain. This is confirmed by mean bivariate polar plots for Magadino-Cadenazzo where elevated EBC WB concentrations were found when wind speeds were very low (≤ 2 m s −1 ), while EBC TR sources were identified for a diverse range of wind speeds and directions (Fig. 6).
For most monitoring sites, the EBC WB components displayed elevated concentrations in the late evening with a decline in the early morning, consistent with domestic heating demands (Grange et al., 2013). In most cases, this evening peak was followed by a decrease and a lesser morning increase in EBC WB , likely due to the reignition of woodburning appliances in the morning. However, at the Bern-Bollwerk urban-traffic site, the morning EBC WB peak was reminiscent of the EBC TR peak, indicating a slight contamination of EBC WB by EBC TR which the aethalometer model was unable to separate. At Magadino-Cadenazzo, the rural monitoring location which is exposed to a high load of woodburning PM, the diurnal cycle was very clear and strong.
Rigi-Seebodenalp's EBC diurnal cycles were different than the other monitoring locations (Fig. 5). This monitoring site is located at an elevation of 1031 m, is isolated from significant local emissions, and is intermittently in the boundary layer and therefore, at times is not influenced by surface source activities . These site attributes resulted in EBC diurnal cycles being driven primarily by boundary layer evolution rather than local source strength and activities. When the convective boundary layer grew in vertical extent and exceeded the site's elevation, EBC was mixed to the site, but when the stratified nocturnal boundary layer decoupled the site from surface-based emissions, concentrations remained low (Fig. 5). This feature was clearer for EBC TR than EBC WB , suggesting that traffic, rather than woodburning emissions, are more influential at Rigi-Seebodenalp.

Aethalometer model failure
San Vittore is a small, rather rural settlement located in the Mesolcina Valley in the south of Switzerland (Fig. 2), an area which experiences high concentrations of PM and BC from residential woodburning. The diurnal plots for San Vittore demonstrate, on average, negative mass contributions for EBC TR when EBC WB concentrations were high in the evening (averages shown in Fig. 5). This is an implausible situation and indicates the failure of the aethalometer model to correctly apportion the traffic and woodburning EBC sources. When the time series of the EBC TR and EBC WB are plotted, the negative contributions are very clear during the colder months, November to December (Fig. 7).
These negative contributions can be managed to some extent by altering the woodburning Ångström exponent (α WB ), but only with justification. In the case of San Vittore, the α WB was increased to ≈ 2.0 (a very high value) in an attempt to eradicate the negative EBC TR mass contributions (Fig. 7). This solution only partially resolves the negative EBC TR concentrations but also results in the EBC TR being estimated at implausibly high concentrations for this location. This is most evident during the cooler periods where EBC WB concentrations are high due to a heavy burden of wood smoke, but EBC TR concentrations are also high. Based on San Vittore's characteristics, we concluded that the aethalometer model is unable to separate the two EBC components in this time series.
The San Vittore example represents a very clear example of when the aethalometer model fails and the results cannot be used further for data analysis. The aethalometer model fails in this example because the monitoring site is exposed to freshly emitted wood smoke from the nearby residential properties. Woodburning aerosol which has not been exposed to any or very little atmospheric ageing has different special properties than that exposed to atmospheric processing (Jimenez et al., 2009; which was indicated in Fig. 4) and is the explanation for the aethalometer model failure. The aethalometer model's two-source approach is most likely insufficient for the freshly emitted wood smoke which may need to be treated as a third distinct source. A receptor modelling approach (such as PMF) may be a technique which could resolve these three sources, but an investigation of this was outside of the scope of this study. Observations from the San Vittore site were excluded from any further analysis but have been presented to demonstrate that although the aethalometer model can be a useful technique, it is pragmatic and is not appropriate in all cases.

EBC dependence on air temperature
BC emissions from traffic and woodburning sources are not only different in respect to their timing as discussed in Sect. 3.2 but also regarding their dependence on ambient air temperature. Although vehicle emissions of PM are  known to be temperature dependent (Jamriska et al., 2008;Weilenmann et al., 2009), this effect can be expected to be less pronounced when compared to residential heating emissions which have a quite distinct "heating threshold temperature" where emissions will increase as the air temperature decreases further to meet an increasing space heating energy demand. Additionally, in warmer periods, there is very little or no emissions of BC from residential heating.
These different source patterns were clearly demonstrated for the Swiss EBC monitoring sites with Magadino-Cadenazzo EBC WB concentrations being negatively correlated with air temperature and with increases in mean EBC WB concentrations once air temperatures were below 12 • C (Fig. 8). Temperatures of 17 • C and above resulted in EBC WB being less than 0.15 µg m −3 at the same location. Payerne's EBC components showed patterns similar to Magadino-Cadenazzo, while the sites which were influenced more by traffic emissions, Bern-Bollwerk and Zürich-Kaserne, still demonstrated a greater ambient temperature dependence for EBC WB than EBC TR despite EBC TR concentrations being higher for most temperatures.
Due to Rigi-Seebodenalp's elevated location, the patterns observed are different than at the other monitoring sites; in particular as temperature increased so too did EBC TR con- Figure 7. EBC TR and EBC WB hourly time series between 2014 and 2018 at San Vittore using two values for α WB . The α WB value of 1.68 was used for the analysis but showed negative contributions of EBC TR when EBC WB concentrations were high, and using an α WB value of 2.0 showed implausible EBC TR concentrations at this site. Figure 8. EBC TR and EBC WB dependence on air temperature for five equivalent black carbon (EBC) monitoring sites in Switzerland. CI: confidence interval. centrations (Fig. 8). This was a result of the mixing of trafficsourced pollutants to the monitoring location and shows that Rigi-Seebodenalp is influenced more by traffic emissions than those originating from woodburning activities, which was also suggested by Fig. 5.

Trend analysis
After the aethalometer model split EBC into the two components, EBC TR was found to be significantly decreasing across all monitoring sites in Switzerland between 2008 and 2018, with the exception of the rural-mountain monitoring site, Rigi-Seebodenalp (Fig. 9). Bern-Bollwerk, an urban-traffic site was the most polluted EBC monitoring site (Table A2) and saw the greatest reduction in EBC TR with a decrease of −0.13 µg m −3 year −1 . Magadino-Cadenazzo, Payerne, and Zürich-Kaserne demonstrated modest significantly decreasing EBC TR trends ranging from −0.028 to −0.059 µg m −3 year −1 (Fig. 9).
EBC WB trends were different than those observed for EBC TR , even when comparing the same monitoring sites (Fig. 9). The urban monitoring locations, Bern-Bollwerk and Zürich-Kaserne, demonstrated significantly decreasing EBC WB trends but at a smaller magnitude when compared to the EBC TR components at the same sites. Payerne's EBC WB concentrations were also found to be decreasing significantly, but at a minute rate (−0.005 µg m −3 year −1 ). The monitoring site which experiences the greatest woodburning PM load, Magadino-Cadenazzo, south of the Alps (Fig. 2), showed no significant trend in the EBC WB component (Fig. 9). This is in contrast to the traffic-sourced BC. Figures 5 and 8 suggested that Rigi-Seebodenalp's isolation from localised BC sources made the site distinct from the other EBC monitoring sites in Switzerland. The trends for Rigi-Seebodenalp were also distinct and showed no significant trend in both the EBC TR and EBC WB components (Fig. 9).
EBC TR trends at Payerne, Magadino-Cadenazzo, and potentially Zürich-Kaserne suggest that breakpoints in the observations are present at the middle of 2013 (Fig. 9). Curiously, these changes were not observable in the absorption observations themselves but could be detected in the calculated α value with adaptive Kolmogorov-Zurbenko filters (KZA filters; Zurbenko et al., 1996; Fig. A4). The dates of these small breakpoints could not be robustly traced to operational activities, but they might have been caused by different batches of aethalometer filter tape which were continuously introduced across the monitoring network in 2013. Therefore, these data suggest that these breakpoints are instrument artefacts. The identification of these features reinforces that the use of the aethalometer model is very useful; it is a pragmatic technique which requires careful evaluation.

EBC/PM 2.5 ratios
BC is almost exclusively an anthropogenic pollutant with generally minor and intermittent contributions from wildfires. It is therefore useful to explore the contribution of BC emission sources on ambient PM concentrations.
The seasonal EBC/PM 2.5 ratios for the five Swiss EBC monitoring sites between 2014 and 2018 are shown in Fig. 10 and show that the ratios were variable among the monitoring sites and seasons. Bern-Bollwerk, the urban-traffic site, had the highest EBC/PM 2.5 ratio, while Payerne and Rigi-Seebodenalp had the lowest. Again, Fig. 10 gave plausibility to the aethalometer model's source apportionment abilities because features such as EBC WB had very low contributions in the summer for most sites, even for the biomass-burning-dominated Magadino-Cadenazzo location. However, as suggested by Fig. 5, Bern-Bollwerk's EBC WB contribution to PM 2.5 was most likely too high, indicating that the aethalometer model was unable to fully separate the two EBC TR and EBC WB fractions completely for this particular monitoring site.
On first inspection, the EBC TOT contribution to PM 2.5 mass seems modest, ranging from 6 % to 14 % depending on site and season. Compare these values to the largest European urban areas, Paris and London, where BC/PM 2.5 ratios (at single heavily trafficked locations) of 43 ± 20 % and ≈ 50 % respectively have been reported (Ruellan and Cachier, 2001;Grange et al., 2016).
The discussion above is focused on the contribution of EBC to PM. In addition, BC should be to reduced to the lowest possible level because of its classification as a Group 1 carcinogen by the International Agency for Research on Cancer (IARC; entered as soot; International Agency for Research on Cancer, 2019). The reduction of soot emissions also has co-benefits for reducing climate warming effects, and this can also be used for the motivation for the implementation of soot control at the source.

Conclusions
Using aethalometers for EBC monitoring allows for the application of the aethalometer model to split EBC into EBC TR and EBC WB components. The aethalometer model is a useful and pragmatic data processing technique, but it requires evaluation before using the model's outputs for further data analysis activities. Based on the results presented using Swiss aethalometer observations from six sites between 2008 and 2018, these recommendations can be offered to other data users: -Values used for the mass absorption cross-section (MAC) coefficients and the Ångström exponents (α) should be checked with the observational record, which is being analysed to ensure they are plausible for use with the particular data set. EC observations are how-  ever required for the calculation of MACs, and therefore, the expansion of monitoring networks to include regular EC samples is recommended for better exploitation of aethalometer absorption data.
-Evaluate the calculated EBC TR and EBC WB by investigating the presence of negative mass contributions, their diurnal cycles, and ambient temperature dependence to ensure these features are consistent with what is known about the sites' PM load and relevant physical and chemical atmospheric processes.
-If available, put the EBC TR and EBC WB in context by using PM mass concentrations.
These recommendations will help those who are interested in applying the aethalometer model, which has a rather low barrier to entry.
Despite the failure of the aethalometer model to produce value outputs for one of the six monitoring sites analysed, EBC TR and EBC WB were successfully analysed separately at five locations, and different trends were observed for the two EBC components. EBC TR concentrations significantly decreased across all monitoring sites between 2008 and 2018, with the exception of an isolated rural-mountain monitoring location at a maximum rate of −0.13 µg m −3 year −1 . The trends of EBC WB were more variable, and despite some significantly deceasing trends of some monitoring sites, a site known to be heavily burdened by wood smoke showed no significant trend over the 2008 and 2018 monitoring period. The EBC/PM 2.5 ratios for the five monitoring sites showed location and seasonal variability, and the EBC contribution to PM 2.5 was between 6 % and 14 %. This EBC/PM 2.5 ratio is low when compared to other heavily trafficked sites in Paris and London. The aethalometer model is a pragmatic and useful data processing technique, but full evaluation is needed to ensure the results are suitable for further data analysis.
Appendix A Table A1. Factory mass absorption cross-section (MAC) coefficients for the AE33 aethalometer (Magee Scientific, 2016    Code and data availability. The data sources used in this work are described and are publicly accessible in a persistent data repository (https://doi.org/10.5281/zenodo.3626658; . Additional information about these data is available from the authors on reasonable request.
Author contributions. SKG and CH developed the research questions and conducted the data analysis. HL and AF collected and supplied additional data. SKG and CH prepared the paper with input from LE.
Competing interests. The authors declare that they have no conflict of interest.