Combination of Different Approaches to Infer Local or Regional Contributions to PM 2.5 Burdens in Graz, Austria

: In early 2017 high particulate matter (PM) levels were observed across mid ‐ Europe, including Austria. Here we characterize PM pollution in the city of Graz during January to March 2017, a period with substantial exceedances (34 days) of the European Union (EU) PM 10 short time limit value. This study evaluates whether the observed exceedances can be attributed to the accumulation of pollutants emitted by local sources or to a larger scale pollution episode including transport. The analyses are based on the ratios of PM 10 concentrations determined at an urban and background site, and the analyses of chemical composition of PM 2.5 samples (i.e., water soluble ions, organic and elemental carbon, anhydro ‐ sugars, humic ‐ like substances, aluminum, and polycyclic aromatic hydrocarbons). Source apportionment was realized using a macro ‐ tracer model. Overall, the combination of different approaches (PM 10 ratios, chemical composition, and macro ‐ tracer derived source apportionment) enabled a conclusive identification of time periods characterized by the accumulation of emissions from local sources or regional pollution episodes.


Introduction
Today it is well understood that exposure to air pollution, especially particulate matter (PM), has negative impacts on human health. Epidemiological and toxicological studies report a relation between exposure to elevated PM concentrations and respiratory and cardiovascular diseases [1]. To protect human health limit values for PM2. 5 and PM10 have been included in the European Directive 2008/50/CE. Even lower threshold values are recommended by the World Health Organization [2]. In addition to direct limit values for PM2.5 and PM10, the European Directive 2004/107/CE sets target values for several toxic substances related to PM, including polycyclic aromatic hydrocarbons (PAHs). Moreover the Sustainable Development Goals (SDGs) defined by the United Nations to achieve a better and more sustainable future for the whole globe by 2030, list air quality within the targets and explicitly mention PM concentrations in urban environments as an indicator to evaluate the current situation [3].
Despite reductions in the PM burden during recent decades, many regions within the European Union are still struggling to attain the PM10 limit value and the violations of air quality regulations pose a challenge on both the local and regional scale. To reduce ambient air pollution, it is necessary to understand the emitter-receptor relationship and to have detailed knowledge about individual source contributions. PM is a complex mixture, and the ambient PM concentration at a specific location depends on both local emissions and the background pollution load [4,5]. In January and February 2017 high PM levels were observed throughout mid-Europe [6]. Enhanced PM burdens were also observed in Austria. Within Austria the city of Graz is known as a 'PM hotspot', prone to enhanced PM levels due to the combination of orography, in-situ emissions, and long-range transport of air masses [7].
During the first months of 2017, the European Union (EU) short time limit value for the PM10 concentrations (50 μg/m 3 as daily mean) was exceeded several times in Graz. Peak abundances even exceeded 100 μg/m 3 . This study evaluates whether the observed exceedances can be attributed mainly to local effects (i.e., an urban impact) or regional pollution transport leading to a high initial pollution load and how the chemical composition of PM and the contribution of main aerosol source to total PM mass are affected by this interplay, aiming at health related issues as well as possible reduction measures. All analyses are based on input data readily available from a state-of-the-art air quality monitoring network, covering both polluted and background sites, thus allowing a transfer of the analytical approach presented to other regions.

Sampling Site, Particulate Matter Sampling, and Gravimetric Mass Determination
The city of Graz is located southeast of the Alps, within the valley of the river Mur, surrounded by mountains. The PM filter samples analyzed in this study were collected within the Styrian air quality network, maintained by the provincial government of Styria (Department 15). The prime monitoring sites of interest are Graz Don Bosco (DB, LON 15.41643°, LAT 47.05702°, 358 m asl) and Bockberg (BB, LON 15.49583°, LAT 46.87139°, 449 m asl). While DB represents an urban environment, BB is not directly influenced by anthropogenic emissions and thus represents a regional background site. At DB PM10 and PM2.5 were sampled with Digitel high volume samplers (air flow at operating condition: 30 m 3 /h) and PM mass was determined gravimetrically according to EN 1234:2014. PM2.5 collected on quartz fiber filters (Pall Life Sciences) was provided for further chemical analyses. At BB PM10 mass was determined with a beta attenuation mass monitor (MetOne BAM 1020; EN 16450:2017) with a flow rate of 1 m 3 /h. This study focuses on PM abundances during the time period of 2 January-31 March 2017. For chemical analyses 24 h samples were grouped to sample pools comprising 2 to 10 days with comparable DB PM10 concentrations. The pools were formed on the basis of PM10 mass concentrations determined at DB, and the ratio of PM10 concentrations between DB and BB. PM10 concentrations were only used to select sample pools or to calculate PM10 mass ratios to compare the urban (DB) and background (BB) site. Chemical analyses were solely performed on PM2.5 filters collected at DB.

Statistical Analyses
All statistical analyses described here were performed with the statistical software R (Version 3.6.3). Averages refer to time weighted means of the pooled samples; median values and standard deviation of all pooled samples are presented. Correlation coefficients between analytes were determined using the Spearman's rank correlation method and considered as statistically significant with a p-value less than 0.05.

Source Apportionment
The definition of PM sources follows the review of Belis et al. [11] who lists six major source categories for PM in Europe (atmospheric formation of secondary inorganic aerosols, traffic, resuspension of crustal/mineral dust, biomass burning, (industrial) point sources, and sea/road salt). As the limited number of samples analyzed for Graz (14 sample pools covering 64 days) does not permit source apportionment based on a multivariate technique like positive matrix factorization, an exploratory method, a macro-tracer approach, was used to investigate the contributions of the main PM sources in Austria. These were identified during previous work [4,12,13]. This macro-tracer approach is based on concentrations of marker compounds and respective conversion factors optimized for the region of interest [14][15][16].
Secondary inorganic aerosols (SIA) include nitrate, sulfate, and ammonium generated from precursor gases NOx, SO2 and NH3. A conversion factor of 1.1 is applied to the sum of respective ions (∑(NO3 − , SO4 2− , NH4 + )) to account for water uptake due to humidity. SIA cannot be attributed to a specific source of PM but indicate an influence of regional-scale or long-range transboundary anthropogenic pollution [11,17]. The contribution of wood burning aerosol (WS) to PM was determined based on the concentrations of the anhydro-sugar levoglucosan (Lev) as proposed by Simoneit et al. [18] and since then widely deployed. As emissions depend on the type of wood and stove used, as well as combustion conditions, the conversion factor needs to be adjusted to local conditions: yielding 10.7 [15]. The related emissions of elemental carbon (i.e., ECWS) were assumed to be 5% of WS, which is within the lower range of data reported in the literature [19], but matched to observed EC values. A proxy for PM emissions attributable to traffic was deduced from ECtr (ECtr = EC − ECWS) concentrations. 'Traffic primary' was calculated as the sum of traffic exhaust (ECtr multiplied with a conversion factor of 1.33, accounting mainly for organics) and the contribution of abrasion (which was estimated to be 30% of traffic exhaust). The respective conversion factors used are based on tunnel measurements in Austria [20,21]. We highlight that this proxy, called 'traffic primary', does not include secondary inorganic aerosols related to traffic (e.g., particulate nitrate formed of precursor NOx), as those contributions are already included in SIA, which denote several types of regional-scale anthropogenic pollution. In Austria traffic emissions accounted for 51% of the total NOx emissions in 2017 [22], resulting in a pronounced share of traffic-related SIA. We note in passing that this approach neglects possible contributions of coal and oil to overall EC concentrations, which is justifiable according to the total emissions reported for Austria [23]. Furthermore, mineral dust (MD) is separately identified. MD is emitted by both anthropogenic and natural sources. MD is related to resuspension of split, erosion, construction work or agricultural activities, and is preferably found in PM10 samples [10]; thus, only a minor contribution is expected for our data set. MD was determined based on the analysis of calcium and aluminum and the respective conversion factors and geogenic ratios to account for silicates and carbonates [24]. We note in passing that a fraction of MD might overlap with traffic abrasion, but due to the overall minor contribution of MD to PM2.5 no further considerations were made.
Regarding the carbonaceous fraction, WS has been identified as a major contributor to organic carbon during the cold season [13]. OC contributed by WS is also scaled via levoglucosan concentrations [13,15]. Further contributions are formed by HULIS, which indicate the presence of secondary organic aerosols [9]. A partial overlap of HULIS and WS cannot be excluded, although based upon emission profiles [15,25], the contribution of HULIS in wood smoke is below 5.8% and decreases due to newer combustion technologies. While the OC contribution by HULIS was measured directly, OC concentrations related to traffic were scaled by ECtr concentrations [20,21]. After subtracting these contributions, the remaining OCnd (OC not identified) was converted to OMnd (organic material not defined) using a conversion factor accounting for heteroatoms occurring in the organic compounds. Based on studies from Europe and the US [26,27], a conversion factor of 1.5 was applied.
The relationships and factors detailed above are also provided in Table S1. Figure 1 shows the daily PM concentrations at the urban sampling station DB over the study period. The PM10 concentrations observed between 2 January-31 March 2017 ranged from 9.9 to 140 μg/m 3 . On 1 January 2017 (following fireworks on New Year's Eve) a maximum concentration of 157 μg/m 3 of PM10 mass was measured, but due to the influence of fireworks this day is excluded from further evaluations in this study. During the sampling period 34 days occurred with a PM10 concentration exceeding 50 μg/m 3 and thus violating the EU short time PM10 limit value. Of these exceedance days, 21 already occurred in January 2017. Considering the whole year of 2017, 54 days violating the EU short time PM10 limit value occurred, indicating a substantial PM burden in contrast to the general decreasing trend of annually occurring days of exceedances [28]. PM2.5 concentrations at the urban sampling station DB followed the trend of PM10, with concentrations between 5.6 and 120 μg/m 3 , giving an average PM2.5 to PM10 ratio of 0.7.

PM2.5 and PM10 Concentrations
As chemical analysis could not be performed for the whole data set, selected filters were pooled to represent characteristic time periods. The pools were formed on the basis of PM10 concentrations, and time periods selected for chemical analysis are marked with grey shading in Figure 1. In total, 14 sampling pools were analyzed, comprising 64 (out of 89) days of the sampling period. Overall, the highest PM concentrations occurred during January and are represented by six sampling pools (comprising 19 days). February and March 2017 are represented by five and three sampling pools, covering 25 and 20 days, respectively.   Based on the DB/BB PM10 ratios, which indicate differences or conformity between the PM concentrations determined within the city and the background site, we postulate a subdivision of these six periods into two scenarios: (1) DB/BB PM10 ratios larger than 2 represent a pollution period with local character (i.e., a marked impact of urban emissions); (2) DB/BB PM10 ratios smaller or equal to 1.5 represent comparable PM concentrations on both sampling sites and thus indicate larger scale pollution. In the following we will refer to these scenarios as 'scenario 1' and 'scenario 2', respectively. Table 1 provides detailed information about the six time periods of interest (i.e., average temperature, average DB PM10 and DB PM2.5 concentrations, and DB/BB PM10 ratios). PM ratios smaller than 1.5 or larger than 2.0 were found in March as well but were not evaluated in detail because of PM10 concentrations below the daily limit value. Concentrations of both PM10 and PM2.5 at DB were higher during periods with temperatures below 0 °C, independent of the scenario. The PM2.5/PM10 ratio at DB was found to be lower during time periods characterized as scenario 1 (average PM2.5/PM10 ratio 0.68) than for scenario 2 (average PM2.5/PM10 ratio 0.82).

Chemical Composition
The mean values, standard deviations, and ranges for quantified PM2.5 constituents are summarized in Table 2. A list showing ambient concentrations for all analytes and sample pools is given in Table S2. Mean concentrations of ionic species, as well as standard deviations, medians, and ranges are summarized in Table 2. Quantified anions contributed an average of 27% to the total PM2.5 mass, with maximum contributions of 45% (27 January-30 January 2017) and minimum contributions of 18% (27 March-31 March 2017). Total cation mass contributions were 12% on average and ranged between 8% and 14% of the total PM2.5 concentration. Overall equivalent concentrations of cations and anions are well balanced with an average value of 0.97 (Figure 3a), but show a distinct temporal evolution (Figure 3b). At the beginning of the sampling period a cation deficit emerges, which vanishes in early February with changing ambient conditions. This shift is primarily caused by decreasing concentrations of anions and ammonium, while the concentrations of the other cations remain rather constant over time, or in the case of magnesium and calcium, even show an increasing tendency. We link these changes to an increasing influence of resuspension of mineral dust, which will occur predominantly in the coarse fraction, but slightly affects PM2.5 as well [4]. Regarding the anion deficit we assume the occurrence of different anions like carbonate, but also organic acids and phosphate, which were not quantified in this study.
Elevated chloride concentrations were only found during January, concurrent with above average sodium concentrations. During January the molar ratio of chloride to sodium showed excess chloride with an average molar ratio of 1.4, while the opposite became true in February and March, where average molar ratios of 0.52 and 0.15, respectively, were observed. While the concurrent increase of sodium and chloride might indicate the use of de-icing salt applied to the urban environment, the influence of combustion sources can also contribute to the excess chloride [29,30]. Nitrate was identified as the dominant anion over the whole sampling period, with concentrations ranging between 1.6 and 20.0 μg/m 3 , while concentrations of sulfate ranged between 0.73 and 12.9 μg/m 3 and concentrations of ammonium, the most abundant cation, ranged between 0.12 and 8.0 μg/m 3 . Nitrate concentrations emerge as highly correlated with sulfate (r = 0.93) and ammonium (r = 0.97), indicating the formation of (NH4)2SO4 and semi-volatile NH4NO3 [31]. Neutralization ratios based on equivalent concentrations of ammonium and sulfate showed excess ammonium for all but two sample pools with rather low PM mass in March (1 March-10 March 2017; 27 March-31 March 2017). If nitrate is included, the neutralization ratios (NH4 + /∑(SO4 2− + NO3 − )) show a maximum value of 0.83, indicating that a fraction of the nitrate present is associated with other ions such as potassium (r = 0.75). These conditions allow to estimate the concentrations of (NH4)2SO4 and NH4NO3 based on the molar concentrations of sulfate and excess ammonium (each multiplied with the respective molar mass of the salt). During time periods with average ambient temperature below 0 °C the contribution of NH4NO3 to major ions (ammonium, nitrate, and sulfate) was on average 33% (up to a maximum of 47%). Conversely, during warmer periods (average temperature above 0 °C) the contribution of NH4NO3 never exceeded 26% and appeared even as negligible during some periods. The contribution of (NH4)2SO4 to major ions showed no variation with temperature (average contribution of 37% during time periods with temperatures below 0 °C and 33% during time periods above freezing, respectively). Figure 4 shows the temperature dependency of major ions presented as mass concentrations ( Figure 4a) and relative contributions to PM2.5 mass (Figure 4b). The sampling periods representing the two scenarios defined before are indicated with black squares (scenario 1) and black triangles (scenario 2). Summed ammonium, nitrate, and sulfate contributed between 17% and 51% of PM2.5. A weak tendency is found towards higher contributions during periods with lower air temperatures (Figure 4b). Independently from a temperature dependence, pronounced differences in the contribution of major ions to the overall mass of PM2.5 were observed, when classification according to the DB/BB PM10 ratios is considered. During time periods with DB/BB PM10 ratios ≤ 1.5 (scenario 2) contributions were rather high (36-40%), compared to periods with DB/BB PM10 ratios > 2 (scenario 1; 17-19%). Even during the coldest time span of the study period (11 January-12 January 2017, average T= −7.8 °C), the contribution of nitrate, sulfate, and ammonium yielded only 17%. The classifications according to the DB/BB PM10 ratio indicate a more pronounced influence of secondary inorganic ions during time periods characterized by transport of air pollutants. The differentiation according to the DB/BB PM10 ratios is also reflected by the neutralization ratios (NH4 + /∑(SO4 2− + NO3 − ). These are higher during periods characterized by scenario 2 (average value of 0.72) than scenario 1 (average neutralization ratio of 0.49). Likewise, the contribution of NH4NO3, calculated as above, to major ions differs. It was found to be higher during time periods of scenario 2 (30%) than during time periods of scenario 1 (16%), although the time periods representing scenario 1 were characterized by lower temperatures. Average concentrations of sodium, calcium, and chloride are slightly higher during time periods with DB/BB PM10 ratios > 2 (i.e., when urban sources dominate). Overall ion concentrations observed during winter 2017 showed slightly lower contributions of nitrate and sulfate to PM2.5 (15% and 8%, respectively) than previously observed [12]. During time periods with elevated PM2.5 concentrations (>50 μg/m 3 ), the average contributions of nitrate to PM2.5 mass were in a comparable range for the present and previous study (accounting to 18% in 2017 and <20% of PM2.5 mass in 2000/2001). For sulfate, lower contributions to PM2.5 (11%) were observed during time periods in 2017, while higher contributions (>20%) were reported in 2000/2001 [12]. This difference can be attributed to generally decreasing SO2 emissions [32].

Carbonaceous Compounds
Average OC and EC concentrations as well as the standard deviations, median values, and ranges are summarized in Table 2. Figure 5 shows the correlation between OC and EC concentrations, as the OC/EC ratio is commonly used to distinguish between different sources of carbonaceous material [33]. Previous studies reported lower OC/EC ratios (up to 2.3) to be associated with fossil fuel burning, while OC/EC ratios associated to emissions from biomass combustion are higher and exhibit a wider range (up to 16) [29,33,34]. The OC/EC ratios found in this study ranged from 2.4 to 11. The classification according to DB/BB PM10 ratios characterizes differences in chemical composition and possible aerosol sources ( Figure 5). Scenario 2 (DB/BB PM10 ratios ≤ 1.5) shows higher OC/EC ratios (5.2 to 11), characteristic for aerosols related to biomass burning, while scenario 1 (DB/BB PM10 ratios > 2) yields samples with lower OC/EC ratios (3.2 to 4.6) indicating the influence of biomass burning as well as fossil fuel combustion and thereby a plausible, stronger urban impact. The highest OC/EC ratios (>5) occurred when ambient temperatures were rather low (until mid-February), pointing to a pronounced influence of biomass combustion. Yet, this was not accompanied with lower relative contributions of inorganic ions, as could be expected because of the increased source for carbonaceous aerosols. Figure 6 shows the correlation between OC/EC ratios and the contribution of ammonium, nitrate, and sulfate (summarized as SIA, see Section 2.4) to PM2.5. Scenario 1 shows generally lower OC/EC ratios, and also lower contributions of SIA. The opposite can be observed during time periods of scenario 2 where the more regional character of the pollution event characterized by similar PM concentrations at the urban and the background site leads to higher OC/EC ratios and an above average contribution of SIA. The influence of biomass combustion indicated by the OC/EC ratios was further investigated using the molecular tracer levoglucosan. Emission studies linked levoglucosan (Lev), mannosan (Man), and galactosan (Gal) to biomass combustion processes [18,35], whereas the ratios of Lev/Man and Lev/Gal provide information about the combustion material [13,30]. Lev was found to be the most abundant saccharide during the sampling period with an average ambient concentration of 0.86 μg/m 3 (0.22 to 2.6 μg/m 3 ), contributing 1% to 4% of total PM2.5 mass. Atmospheric concentrations of Lev were higher during cold periods and decreased with time and increasing temperatures. Man concentrations were on average eight times lower than those of Lev, with concentrations ranging from 0.03 to 0.35 μg/m 3 . Average Lev/Man ratios yielded 7.9 indicative for the use of softwood and hardwood as fuel for biomass combustion [15]. Gal concentrations were mainly below the limit of detection. Figure 7a shows a high correlation between Lev and OC concentrations (r = 0.93) for the whole data set, suggesting a major fraction of OC originating from biomass burning across all samples. The correlation between Lev concentrations and OC/EC ratios (Figure 7b) shows larger scatter (r = 0.74) and some difference between samples representing different DB/BB PM10 ratios. Samples representing scenario 1 (black squares) show variable Lev concentrations but lower OC/EC ratios, while time periods of scenario 2 (black triangles) show higher OC/EC ratios despite Lev concentrations similar to scenario 1. This is mainly driven by the relative contribution of EC and the occurrence of other EC sources. During scenario 2 the average contribution of EC to PM2.5 mass is 4%, while it averages at 7% during scenario 1. In addition, the correlation of OC with potassium (r = 0.76, Figure 7c) suggests biomass burning as a possible source of particulate matter [35,36]. Two data points associated with scenario 2 do not follow the general trend. Period F shows rather high potassium concentrations and during that time period the influence of long-range transport of desert dust appears likely [28]. No explanation can be given for the rather low potassium concentrations determined for data point E.

PM2.5 Associated PAHs
Ten PAHs including fluoranthene, benzo(a)anthracene, chrysene, benzo(b,k)fluoranthene, benzo(e)pyrene, benzo(a)pyrene, perylene, indeno(1,2,3-cd)pyrene, benzo(ghi)perylene, and dibenz(ah)anthracene were quantified for all time periods. The sum of all 10 quantified PAHs is hereinafter abbreviated as PAHtotal. Table 3 presents average PAH concentrations as well as the standard deviations, median values, and ranges for each congener associated to PM2.5 samples at Graz DB. Concentrations of benzo(a)pyrene (BaP), the primary representative of PAHs and a known carcinogen (class 1 by the International Agency for Research on Cancer (IARC) [37]), ranged between 0.51 and 8.4 ng/m 3 (2.2 ng/m 3 ). The highest concentrations of BaP was measured during 2 January-3 January 2017, a time period still showing influence of local activities at New Year's Eve, where PAHtotal concentrations reached 102 ng/m 3 . PAHtotal concentrations varied from 4.6 to 56 ng/m 3 (17 ng/m 3 ) over the study period and are comparable to recent results reported for other European cities [38,39]. With proceeding time, a decreasing trend of PAHtotal concentrations was observed and highest concentrations could be observed for benzo(b,k)fluoranthene (B(b,k)F) followed by benzo(a)anthracene (BaA) and BaP. The toxicity of PAHs is known to differ dependent on molecular weight, and many of the heaviest PAHs are categorized by IARC as already known (class 1), probable (class 2A) or possible (group 2B) human carcinogens [37]. Therefore, we grouped the quantified PAHs in middle molecular weight PAHs (MMW, containing 4-ring congeners) and high molecular weight PAHs (HMW, containing 5-and 6-ring congeners). PAH profiles are given in Figure 8. Due to their high volatility, low molecular weight PAHs containing 2-and 3-ring congeners are not included in this work. During all time periods HMW PAHs were identified as the major components with PAH concentrations ranging from 2.9 to 42 ng/m 3 corresponding to 60-79% of the PAHtotal concentrations. The ratio of PAHtotal to PM2.5 ranged between 0.19 and 0.88 ng/μg and the corresponding ratios of BaP/PM2.5 varied between 0.02 and 0.13 ng/μg during the whole sampling period. With increasing time not only the observed PAH patterns changed, also the PAHtotal/PM2.5 ratios decreased (Figure 8). Generally increasing temperatures result in lower PAH emissions, due to reduced residential heating activities and higher degradation of PAHs. This is already visible towards the end of the study period, which agrees with previously reported differences between summer and winter [38,39]. A clear difference in PAHtotal/PM2.5 ratio was observed when a filter pool belonging to scenario 1 was analyzed. During that time period (21 February-23 February 2017) PAH concentrations within the city were found to be much higher than at the background site, leading to an increase in PAH/PM2.5 ratio.

Source Apportionment with a Macro-Tracer Approach
The concentrations of the PM source groups defined in Section 2.4 are shown in Figure 9. For all sample pools WS and SIA emerge as the dominant sources, although with alternating dominance. WS contributions to PM2.5 mass ranged between 10% and 43% (2.4-27.7 μg/m 3 ) and contributions from SIA ranged between 19% and 57% (2.8-44.9 μg/m 3 ). Contributions of 'traffic primary' ranged between 1.6% and 16% (1.4-6.0 μg/m 3 ). Contributions of HULIS and MD showed rather low average contributions to PM2.5 mass of 2.3% (1.0 μg/m 3 ) and 4.1% (2.2 μg/m 3 ), respectively. Uncertainties in conversion factors lead to under-and over-estimations of PM sources in certain cases. However, during the study period the carbonaceous fraction was slightly overestimated (1.8 μg/m 3 ) only between 2 January-3 January 2017 via this macro-tracer approach. As the discussions of the chemical composition revealed marked differences between the two scenarios, source apportionment via the macro-tracer approach is shown in more detail for these samples ( Figure 10). Time periods of scenario 1 (DB/BB PM10 ratios > 2; Figure 10a) showed higher average contributions of 'traffic primary', contributing 7-11% to PM2.5 mass (average value of 10%), while contributions of SIA ranged from 19-21% (average value of 20%). Contributions of WS ranged between 28% and 43% (average value of 35%). Furthermore, the mass of not identified PM was found quite high. During time periods of scenario 2 (DB/BB PM10 ratios  1.5; Figure 10b), SIA increased and ranged from 40-44%, even though these days were not the coldest during the study period. During this scenario the contribution of 'traffic primary' decreased to an average value of 4%. Contributions of WS ranged from 15-30% (average value of 20%), and thus are lower than during samples of scenario 1. The time period of 9 February-17 February 2017 shows elevated contributions of MD and also OMnd, which can be linked to an input of desert dust identified during this time period [28]. Generally speaking, MD did not show pronounced variations during either scenario considered (4% during scenario 1 and 5% during scenario 2). The same was found for HULIS, with average contributions of 2% and 3%, respectively. Contributions of OMnd were found to be higher during time periods of scenario 2 (9%) than during scenario 1 (2%).
Thus, the simple classification based on the difference in DB/BB PM10 ratios (local vs. regional) is reflected by chemical characterization and findings of source apportionment.

Supporting Information from Meteorological Analyses
The chemical analyses described above detail the PM source apportionment to individual subgroups/components of interest. To supplement the findings of the chemical analyses regarding contributions from local and regional emissions to PM burdens meteorological data is analyzed. Today it is well understood that ambient meteorological conditions strongly modulate PM burdens and determine particulate accumulation and decay [40][41][42][43]. High PM burdens coincide normally with one or multiple of the following meteorological phenomena: stagnant weather conditions, calm winds and limited ventilation, cold ambient temperatures leading to a shallow boundary layer, and temperature inversions inhibiting vertical mixing.
To characterize the meteorological conditions on a regional basis we used the WLKC733 weather classification [44], which provides information on dominant wind direction, cyclonicity at 925 hPa and 500 hPa, and dominance of wet or dry conditions. To characterize local meteorological conditions, we considered measurements of ambient air temperature and wind direction and wind speed. Of these covariates, temperature is most important as it is relatively variable over the study period. Wind speeds are generally modest in Graz during winter times. Over the study period daily mean wind speed was found for most days to be below 1.5 m/s (and daily maxima below 3 m/s). As the city of Graz is prone to temperature inversions, due to its valley location, we diagnosed the presence of a local temperature inversion by determining the air temperature gradient from measurements at the sites DB and Kalkleiten (at ~12 km distance located on a hill encompassing the basin).
From a meteorological perspective stagnant conditions were dominating during January and February 2017. Furthermore, extended periods with cold temperatures (<0 °C) occurred, particularly during January, favoring PM accumulation within the shallow boundary layer (Figure 11b). Chemical source apportionment indicates predominantly widespread regional pollution events besides the sample pools spanning 2 January-3 January 2017, 11 January-12 January 2017, and 21 February-23 February 2017. The regional uniform evolution of ambient air temperatures indicates a relatively homogeneous boundary layer as its depth is determined by ambient temperatures. This is mirrored in the parallel evolution of PM burdens across Styrian air quality monitoring sites (Figure 11a). The highest PM burden is found for the PM pool sampled during the time period of 22 January-23 January 2017. Meteorological data indicate a particularly shallow boundary layer in combination with stagnation. In addition, meteorological data for the city of Graz indicates the presence of a temperature inversion (red highlighted symbols in Figure 11a). Contrasting PM data for DB with those of other Styrian and Austrian sites indicate however, that this local inversion is not the dominant driver of PM pollution. In fact, widespread increases in PM abundance have been reported on a larger scale due to enhanced European pollution backgrounds [6]. In contrast, the PM pools sampled during 2 January-3 January 2017, 11 January-12 January 2017, and 21 February-23 February 2017 carry the signature of an enhanced contribution of local urban emissions to the Graz PM burden. For both periods a local temperature inversion is diagnosed for Graz. Contrasting PM burdens for DB with data from the site BB and Hartberg indicate that this temperature inversion has been a localized phenomenon trapping urban emissions within the boundary layer. For the data pool sampled during 9 February-17 February 2017 stagnant conditions were diagnosed. The regionally homogeneous cold temperatures indicate a shallow boundary layer and combined consideration of urban, rural, and background monitoring sites indicates a regionally enhanced PM background.

Summary and Conclusions
The study presents results of the chemical characterization and source apportionment of atmospheric PM2.5 at an urban monitoring station in Graz, Austria during January to March 2017. The chosen sampling station is a known hotspot for elevated PM10 concentrations in Austria. Given the limited number of samples, source apportionment was performed via an exploratory method using specific marker compounds and conversion factors optimized for the region of interest. Building on this study, future work is suggested, which should include a detailed analysis of longer time periods to allow the application of more comprehensive receptor models such as positive matrix factorization.
Based on the comparison of the PM ratios between an urban (DB) and background (BB) site (DB/BB PM10 ratios), two scenarios were identified: (1) time periods with DB/BB PM10 ratios > 2, indicating dominant urban impact; (2) time periods with DB/BB PM10 ratios  1.5, which characterize pollution episodes observed across a wider area. Pronounced differences in the chemical characterization of the respective samples and in source apportionment show the feasibility of this simple 'dual-site' approach.
Time periods with DB/BB PM10 ratios > 2: Several PM2.5 associated components (e.g., EC, PAHs) were found at elevated fractions during these time periods indicating a dominant influence of urban sources. Higher EC contributions as well as ratios of individual PAHs/PM2.5 related to traffic emissions and combustion pointed towards a pronounced local influence. According to source apportionment the influence of 'traffic primary' is markedly higher, but also WS is a dominant source. This is also indicated by OC/K + correlations and levoglucosan concentrations.
Time periods with DB/BB PM10 ratios  1.5: Time periods have been identified as regional pollution episodes, indicated, for example, by a strong increase of the contribution of SIA. The largest ammonium and nitrate inputs were observed within these time periods, with elevated levels not only found during cold spells. Samples representing scenario 2 also showed differences in OC/K + correlation and the macro-tracer derived source apportionment showed a clear dominance of SIA. WS emerged as the second most important aerosol source.
The comparison of the results of the chemical analyses with supporting meteorological information showed the presence of inversions during time periods with DB/BB PM10 ratios > 2. These local inversions favored the accumulation of locally produced PM. On the other hand, the comparison of meteorological data across a wider area indicated that during other time periods high PM background concentrations occurred. During these episodes, local inversion, even if present, does not appear as the primary driver of elevated PM concentrations.
Supplementary Materials: The following are available online at www.mdpi.com/2076-3417/10/12/4222/s1, Table  S1: List of conversion factors to calculate the source contribution; Table S2: Summary of ambient concentrations of quantified analytes related to PM2.5 at Graz DB.
Author Contributions: Conceptualization of the study was done by B.K. and A.K.-G.; chemical analyses were done by B.K., M.K., T.S., C.H., and A.L.; meteorological data was interpreted and prepared by C.S. and H.R.; writing and original draft preparation was done by B.K.; writing-review and editing was done by A.K.-G. and H.R.; All authors have commented, read, and agreed to the published version of the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research was partly funded by the provincial government of Styria. Open Access Funding by TU Wien.

Conflicts of Interest:
The authors declare no conflicts of interest.