3.1. PM10 Levels at the Stations Graz Don Bosco (DB) and Graz Süd (GS)
Daily PM
10 levels recorded at the stations DB (urban-traffic) and GS (urban-background) exceed the 50 µg/m
3 on 242 and 196 days, respectively, within the six-year period from January 2013 to December 2018, thereby exceeding the 35 day/year limit in six years and one year, respectively. The time series are plotted in
Figure 2 showing that maximum PM
10 values, and hence, exceedance of limit values mostly occur during wintertime. This pattern of the seasonal variation is observed at both stations and is characteristic for stations in this region. The increase of PM concentrations during winter can on one hand be attributed to an enhanced influence of inorganic secondary aerosols due to lower air temperature and increased humidity and also to a stronger influence of wood combustion used for residential heating, as well as regional transport of PM [
22,
24]. Furthermore, meteorological conditions with wintertime inversions favor an enrichment of PM concentrations in the boundary layer. During wintertime inversion conditions are very likely in the area of Graz due to the orographic situation of the region, allowing air pollutants to concentrate over several days until the inversion is cleared out and concentrations drop [
17].
On top of these seasonal variations the influence of long-range transported desert dust may occur. As an example, an intensive desert dust episode in April 2016 is marked by a red star in
Figure 2. The identification, intensity and duration of this event, covering large parts of Austria, was already thoroughly described by Baumann-Stanzer et al. [
12]. Still a number of less pronounced events took place. Such events, identified via the use of WRF-Chem model forecasts [
12], FLEXTRA back-trajectories [
25] as well as the dust ensemble forecasts and forecast comparisons provided by the WMO sand and dust storm warning advisory and assessment system (
https://sds-was.aemet.es/forecast-products/dust-forecasts), are marked in orange in
Figure 2. Additionally optical aerosol properties measured at the high mountain global GAW station Hoher Sonnblick in the Austrian Alps are used for the identification of desert dust (DD) occurrence [
14]. The aim of this study was the quantification of DD contributions for days with PM
10 exceedances only. An additional influence of desert dust might have occurred during periods of lower PM
10 concentrations as well.
3.2. Determination of Background Loads According to the EC Methodology
The EC methodology [
10] for the evaluation of desert dust (DD) contributions is based on two steps. Firstly, dates affected by DD transport have to be identified and secondly, the net dust load (NDL) on the daily PM
10 record has to be quantified. The identification of DD events is based on several supporting information like trajectory analysis, model forecasts and/or satellite data, which have to be screened and interpreted. A number of models, back-trajectory calculations or data bases for satellite data are suggested by the EC methodology, but later evaluations included other sources as well [
7,
11] showing the technical and organizational change in that field. In case the various information portals give diverse results for the identification of a DD event further evaluations and comparisons are needed. The second step, the quantification phase, requires a continuous PM
10 monitoring at a single selected regional background site representative for the region under investigation, providing the reference or background load (BGL). The BGL represent a moving 30-day average value (mean, median or 40th percentile as suggested by the EC) of the 15 days before and 15 days after the investigated day, where days influenced by DD are excluded. The NDL at the regional background site, representing the contribution of DD to PM
10, is then estimated as difference between the measured daily PM
10 load and the calculated BGL:
Within a first approach to estimate the influence of DD for Austria we do not exclude days with DD influence for calculating the NDL. This is a difference to the EC methodology and can be regarded as a more conservative approach. Doing so, we want to consider two main facts. Even if DD is a driver for BGLs at the respective sites and respective days in Austria, other sources might contribute markedly as well and an exclusion of days with influence of DD would eliminate these effects as well. Secondly, modeling tools used for the characterization of DD reveal some uncertainty, which influences the determination of days which could be rejected. This uncertainty might be more pronounced for Austria than for regions closer to the source. The uncertainty is supposed to increase with longer transport time because atmospheric dispersion of air pollutants is a very complex process with very high uncertainties in chemical transport models. Thus, we did not put this elimination step to the beginning of our evaluation procedure. As a matter of fact we will get a potential overestimation of the BGL at the regional background station and hence an underestimation of the according NDL for the respective day. Thus, a potentially lower value of a DD contribution is subtracted compared to reality. Some discussion of the systematical error introduced by this approach will be given later.
Regarding the quantification phase, we compared the 30-day mean, median and 40th percentile for the computation of the BGL at three potential regional background stations (compare
Figure 3). The calculated BGL based on the mean is generally higher than the BGL based on the median or the 40th percentile, regarding all three stations. As may be expected from Equation (1), the estimated NDL decreases as the percentile for the computation of the BGL increases. Again, we decided to stick to a more conservative approach and use the 30-day mean for further computation of the BGL. Using the median or the 40th percentile would reduce the average BGL by 4–11% and 13–21%, respectively, depending on the station used (MB, BB or LB). The BGL calculated based on the 30-day mean varies between 10–50 µg/m
3 for BB and LB and between 4–21 µg/m
3 for MB and is reduced to 5–45 µg/m
3 (BB and LB) and 3–19 µg/m
3 (MB) when the 40th percentile criterion is applied.
3.3. Selection of a Suitable Regional Background Station
According to Equation (1) the NDLs of the respective days, calculated as the difference between PM10 and BGL, will vary around zero. Negative values represent days where the BGL is larger than the central PM10 value and vice versa.
In
Figure 3 and
Figure 4 it can be seen that the time series of MB feature different characteristics than the one of BB and LB. PM
10 values of BB and LB tend to cover a wider range compared to MB and show a seasonality similar to DB and GS with increased values especially during winter. Correspondingly, BGLs of BB and LB (
Figure 3) show this seasonality as well and NDLs (
Figure 4) yield a much larger scatter during the wintertime, than during summer. Negative NDL outliers at BB and LB were found to occur mostly during winter and could be traced back to days at the end of an inversion weather situation where the PM
10 concentration suddenly drops due to an air mass exchange. Days with inversion weather conditions, meaning an inversed vertical temperature profile with lower values at lower altitudes compared to higher ones, were identified using the daily mean temperature at the TAWES stations at the airport of Graz (337 m a.s.l.) and on top of Schöckel (1445 m a.s.l.). If the daily mean temperature at the airport of Graz was lower than the one at Schöckel the day was classified as inversion day. Inversion days occur all years whereas their number varies from year to year. Thus, for single years and during short time measurement campaigns BB and LB might serve as suitable background stations for Graz, but not over a longer period of several years. Positive NDL outliers at BB and LB are of course identic with PM
10 outliers and occur mostly during winter, again reflecting the inversion weather situation where PM
10 concentrations are enhanced. Due to this seasonal dependence, BB and LB are ineligible as regional background stations to compute the NDL. The calculated NDL at these stations cannot be interpreted as such due to the substantial influence of inversion weather conditions. It rather reflects a “pollution load” pointing to a variety of influencing factors from anthropogenic and natural sources.
NDLs of MB are similar between winter and summer and, as a matter of fact, for the whole data set. Negative NDL outliers were found to be associated to days with precipitation rather than inversion weather situations shortly before or after. During precipitation the PM10 concentration suddenly drops and becomes lower than during the previous and subsequent days representing the BGL. Positive outliers of the NDL are again identic with outliers of PM10, but were quite equally distributed between the different seasons. The single outlier for the BGL was observed on July 25, 2013 related to a Saharan dust episode a few days later. Due to averaging over 30 days, also other days are influenced by this Saharan dust episode. Their BGL values were found to be only slightly smaller than the upper whisker and are therefore not highlighted.
Results show that MB represents a suitable background station because it is largely independent of inversion weather situations and hence seasonality. Based on an urban climate analysis study of Graz [
17] it can be assumed that the representativeness of the background condition observed on MB is applicable for the stations in Graz. Thus, we investigate the contribution of desert dust on PM
10 measurements using MB to compute the BGL based on the mean of a 30-day period without extracting desert dust affected periods. This approach represents a modified methodology compared to the one proposed by the EC (40th percentile, exclusion of DD days) which was established for Spain, but the demand for its validation for other countries is clearly highlighted.
Due to its definition (see Equation (1)) the NDL is computed via a subtraction of the BGL from the daily PM
10 value and is likely to vary around zero. BGL values (representing a 30-day average PM
10 load) might be higher than the PM
10 load of a single day due to precipitation or increased vertical mixing after a period of stable atmospheric conditions on that special day, leading to negative NDL values. These negative values are just a consequence of the day to day variability and are, therefore, set to zero for further evaluation. Hence, the NDL of MB can be interpreted as a reduction potential in µg/m
3 due to the influence of desert dust lowering PM
10 levels observed at other stations in the investigated region.
Table 2 gives an overview of the main statistical parameters describing the NDL before and after negative values were set to zero.
3.4. Application of the Modified Methodology
Table 3 shows the number of days exceeding the 50 µg/m
3 daily PM
10 limit at the two investigated stations DB and GS as well as the number of days exceeding this limit after the subtraction of the NDL
MB. Throughout the six-year period 242 days exceed the daily limit value at the urban-traffic station DB, which comes up to five out of six years exceeding the 35 day/year limit. At the urban-background station GS 196 days exceed the daily limit value, whereas only in 2017 the 35 day/year limit was exceeded due to a longer lasting pollution episode in January and February.
Subtracting the NDL
MB without any further evaluation of the presence of DD, 40 days of the 242 days at DB fall below the 50 µg/m
3 limit, reducing the years which exceed the 35 day/year limit to three years (2013, 2017 and 2018) in contrast to the five years given before. For GS 35 days of 196 days fall below the 50 µg/m
3 limit, not influencing the exceedance of the 35 day/year limit. After an identification of DD days based on model result, back trajectories and optical aerosol measurements as described earlier, results reveal that of the 40 and 35 days at DB and GS, only 20 (DB) and 15 (GS) days show an influence of DD (compare
Table 3, columns DD days). Obviously only for those days a subtraction of the NDL
MB is allowed and the solely determination of NDL
MB is not sufficient.
The selection of DD days could be included before subtracting the NDL
MB, corresponding to the EC methodology. This would identify 51 out of 242 days at DB and 37 out of the 196 days at GS with a possible influence of DD. Subtractions of NDL
MB for these special days would again give the number of reductions listed in
Table 3 (5th column for DB and 9th column for GS).
As already discussed, the identification of DD days using models together with back trajectories is a very subjective process. Barnaba et al [
7] developed an automated and user-independent process for the identification of DD days, still using model calculations. As such, a user-independent process was not yet established for Austria. Therefore, we evaluate a method for the identification of DD days based on the time series of the PM
10 measurements at the background site only.
Considering only days with a possible contribution of DD similar features of NDLMB and BGLMB become visible, which define additional criteria to avoid undue reductions even when only NDLMB values are considered. The resulting thresholds (BGLMB: 12 µg/m3 and NDLMB: 10 µg/m3) are based on a descriptive approach, but give a better understanding of the several factors influencing NDLMB.
BGLMB concentrations, determined as a 30-day average, should not fall below 12 µg/m3, a value close to the annual average concentration of PM10 determined at Masenberg. We want to point out that monthly averaged PM10 values of the whole period of 6 years vary between 6–8 µg/m3 for November, December and January and between 10–14 µg/m3 for all other months. Setting a threshold for BGLMB at 12 µg/m3 will preclude the determination of a DD event valid for reduction, during the winter month and also during most of the year 2017. Still a lower threshold value would lead to the identification of events not associated with the influence of DD, but by other sources or general differences between the background site at 1180 m asl and the urban sites in Graz. Furthermore, the calculated NDLMB should be at least 10 µg/m3 and thus account for a substantial part of the PM10 concentration, at least at the background site. This can be regarded as a precondition as the method should identify DD events which influence the air quality in Austria markedly. Ohterwise, low concentration values more affected by random variations could account for undue reductions of PM10 levels. This is especially important during days when daily limit values at polluted sites are exceeded only slightly. The identified NDLMB threshold is close to the 95th percentile of 10.9 µg/m3.
These thresholds were applied to the 242 days exceeding the 50 µg/m
3 daily PM
10 limit within the six-year period at DB only, since this is the more critical station regarding the 35 day/year limit. Now the amount of days still exceeding the limit value is very well comparable as if DD days were identified using model and trajectory analysis. An influence on the exceedance of the 35 day/year limit remains scarce (compare column 3 and 5 in
Table 4) and occurs only in one year (2015) when the adapted methodology using the thresholds led to no reduction below the 35 day/year limit, whereas it would fall below if the identified DD days were used.
The general agreement and disagreement with the evaluation of identified DD days is given as contingency table (see columns 6 to 10 in
Table 4), while evaluations of the single days are given in the supplement (
Table S1).
In total 16 days (6.6%) were identified as true positive (TP) DD days, whereas 184 days (76.0%) were identified as true negatives (TN). This means, that 82% of the days exceeding the daily limit were correctly classified in DD and non-DD days using the thresholds. Regarding the remaining days, five days (2%) were identified as false positive (FP) DD days, while 34 days (13.6%) were found to be false negatives. For 20 of these days some uncertainty of the classification of DD days based on model and trajectory analysis remains, i.e., the visualizations showed slightly different results or the influence can be expected to be very small. The other 14 days show a clear influence of DD. Still all 34 days need further investigations, as this classification as false negatives could easily be an effect of the systematic error introduced by our calculation of the BGL. Adjusting the calculation procedure for every DD event could overcome this limitation, but this adjustment could also make the identification of DD days slightly more subjective.
For 15 out of these 34 days the NDLMB was below 2 µg/m3 while PM10 concentrations at DB ranged from 51 to 139 µg/m3 with a mean value of 64 µg/m3, indicating on one hand that DD is not the main influence on PM10 concentrations at this site, on the other hand that an increase of NDLMB via a reduction of the BGL would allow to reach a reduction below the daily limit value in single cases. A repeated evaluation of the model and trajectory results revealed that for these 15 days the DD identification could easily be biased by subjectivity since results from the models and trajectories do not show a clear picture. An additional four days show an NDLMB below 4 µg/m3 while PM10 concentrations at DB range from 56 to 92 µg/m3—again pointing to both statements. These days could be clearly identified as DD days and thus point to a limitation of the threshold method. Thus, the evaluation of single events remains of great importance, but should be supported by a chemical analysis of the crustal loads as given below exemplarily for two cases.
3.5. Validation of the NDLMB Based on Two Case Studies in 2016
In order to validate the NDL
MB chemical analyses of filter samples collected at three sites (Graz Don Bosco, DB; Graz Ost, GO; Gratwein, GW; compare
Figure 1) were used to determine the crustal loads. This analysis was performed for three days (23.02.2016, 05.04.2016 and 06.04.2016) featuring desert dust intrusion in was thoroughly described earlier.
For the strong transport event of desert dust in April, measured CLs match the calculated NDLMB very well, indicating that the NDLMB is mainly composed of mineral matter. CLs were found to account for more than 70% of the NDLMB with values ranging from 41.2 µg/m3 at GO to 39.0 µg/m3 at GW for the 5th of April 2016, compared to the NDLMB of 50.6 µg/m3 and from 31.5 µg/m3 at GO to 25.3 µg/m3 at GW for the 06th of April 2016, compared to the NDLMB of 35.3 µg/m3.
For the dust event in February, measured CLs are more than 50% lower (17.7 µg/m3 at DB, 16.5 µg/m3 at GO and 10.9 µg/m3 at GW) than the calculated NDLMB (37.1 µg/m3) pointing to additional contributions of other particulate matter sources. This result underlines the need of the threshold criteria. Using those the NDLMB is not allowed to be considered for a reduction, because the BGL is too low. Based on WRF-Chem model forecasts is can also be expected that the SD intrusion during this event was rather weak and other, maybe anthropogenic sources, dominate the NDL.