3.1. PM2.5 Concentrations at the AQMS in Western Macedonia
Figure 3 and
Figure 4 illustrate the box plots of the annual mean PM
2.5 and PM
10 concentrations as registered by the AQMS in the region of Western Macedonia Lignite Center for the years 2019, 2020, 2021, and 2022. The concentrations of PM
2.5 and PM
10 varied among the AQMS. Notably, AQMS S7 and S8 registered the highest PM
2.5 values. As for the years 2020 and 2021, when the lockdown restrictions were implemented, there was a general decrease in the mean, median, and interquartile range of PM
2.5 concentrations at all the AQMS, except for AQMS S6 and S10, where a slight increase was detected in 2020. Similarly, the mean, median, and interquartile range of PM
10 concentrations generally decreased in 2020 at all the AQMS except for AQMS S6.
In
Figure 3, the box plots of PM
2.5 concentrations for the years 2019, 2020, 2021, and 2022 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10 are shown. In 2019, the mean annual PM
2.5 concentrations were 13.7 μg/m
3, 12.8 μg/m
3, 17.0 μg/m
3, 12.7 μg/m
3, 11.9 μg/m
3, 22.9 μg/m
3, 30.5 μg/m
3, 11.9 μg/m
3, and 12.5 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In 2020, the mean PM
2.5 concentrations were 12.9 μg/m
3, 11.9 μg/m
3, 15.2 μg/m
3, 11.6 μg/m
3, 12.8 μg/m
3, 18.9 μg/m
3, 28.0 μg/m
3, 10.2 μg/m
3, and 12.9 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. A minor decreasing trend followed in the year 2021. In 2021, the mean PM
2.5 concentrations were 12.8 μg/m
3, 10.6 μg/m
3, 13.7 μg/m
3, 10.3 μg/m
3, 12.2 μg/m
3, 20.5 μg/m
3, 23.3 μg/m
3, 10.8 μg/m
3, and 11.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In 2022, a slight increase in PM
2.5 concentrations was detected, which is not apparent for all AQMS. Specifically, the mean PM
2.5 concentrations in 2022 were 12.0 μg/m
3, 11.7 μg/m
3, 14.9 μg/m
3, 11.6 μg/m
3, 12.0 μg/m
3, 18.2 μg/m
3, 26.6 μg/m
3, 11.4 μg/m
3, and 12.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively.
As we mentioned above, the AQMS S7 and S8 were constantly recording the highest values. The AQMS S8 was located in Meliti for the years under study after its relocation from Vevi in 2018, while the AQMS S7 was located in Florina. In both areas, district heating networks are currently not available. This means that traditional heating systems such as fireplaces and wood-burning stoves are prevalent, which are, of course, major sources of air pollution in these areas [
12]. It is well known that combustion appliances that burn fossil fuels, wood, or biomass emit a large number of air pollutants (CO, PM, NO
X, volatile organic compounds (VOCs), and polycyclic aromatic hydrocarbons (PAHs)) as well as CO
2 emissions [
32,
33].
Given that the amount of energy needed for heating applications depends on the prevailing weather conditions and that, due to its geographic location, the cold period in Western Macedonia is longer in comparison to the rest of the country, the region is the most energy-consuming region in Greece in terms of heating degree days [
34]. According to Zoras et al. (2007), the climatological monthly mean temperature in Kozani and Florina ranged from values below 3 °C in January to almost 24 °C in July, while the precipitation amounts were higher in November and December [
16].
In contrast to Florina, Kozani, Ptolemaida, Amyntaio, and Filotas are equipped with district heating networks that utilise the thermal load of the neighbouring lignite-fired power stations [
35]. The district energy systems have the advantage of producing heat and/or power with limited environmental impacts and reduced CO
2 and other GHG (greenhouse gas) emissions, as well as having economic benefits because they are highly efficient systems and minimise energy waste [
36]. Pitoska et al. (2021) conducted a questionnaire survey in the city of Ptolemaida and found that the participants recognised the high efficiency of district heating in reducing pollutants and protecting the environment [
35].
Figure 4 displays the box plots of PM
10 concentrations for the years 2019, 2020, 2021, and 2022 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10. In 2019, the mean annual PM
10 concentrations were 22.9 μg/m
3, 21.1 μg/m
3, 26.4 μg/m
3, 21.2 μg/m
3, 18.2 μg/m
3, 29.9 μg/m
3, 40.1 μg/m
3, 22.2 μg/m
3, and 27.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In 2020, the mean PM
10 concentrations were 19.5 μg/m
3, 19.8 μg/m
3, 22.4 μg/m
3, 20.0 μg/m
3, 19.3 μg/m
3, 25.9 μg/m
3, 36.5 μg/m
3, 17.4 μg/m
3, and 23.1 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In most of the AQMS, there was a decrease in annual mean and median PM
10 concentrations in 2021, followed by elevated values in 2022. Specifically, the mean PM
10 concentrations in 2021 were 19.5 μg/m
3, 18.3 μg/m
3, 20.0 μg/m
3, 18.6 μg/m
3, 20.0 μg/m
3, 28.9 μg/m
3, 33.8 μg/m
3, 19.0 μg/m
3, and 20.7 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. The mean PM
10 concentrations in 2022 were 19.3 μg/m
3, 20.9 μg/m
3, 22.4 μg/m
3, 19.5 μg/m
3, 18.4 μg/m
3, 26.9 μg/m
3, 32.6 μg/m
3, 22.4 μg/m
3, and 18.5 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively.
In general, the air pollutants in the region of Western Macedonia exhibit clear seasonality patterns, with their highest values during the winter months and their lowest values during the summer period. The great variations in the air pollutant values were largely affected by the emissions sources in the region and the meteorological conditions. As previously mentioned, Western Macedonia is largely dominated by mining activities and the generation of electricity in lignite-fired power plants. The activities in the open-pit lignite mines, such as mining, transportation of soil and coal, and movement of trucks on unpaved roads, are major sources of PM
10 in the region [
13,
14,
15,
17,
19,
20]. Prior studies have identified the role of these activities in the overall air quality in the region [
13,
14,
15,
16,
17,
18,
19,
20].
Figure 5a,b shows the monthly mean variations in PM
2.5 and PM
10 concentrations and the 95% confidence interval in the mean for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10 for the years 2019, 2020, 2021, and 2022. The monthly mean PM
2.5 and PM
10 concentrations for the months of lockdown (from March to May 2020) and the corresponding months of 2019, 2021, and 2022 are shown in
Table 2 and
Table 3.
The higher mean monthly PM
2.5 concentrations were detected during the cold period (October to April), while during the warm period (May to September), the lowest mean monthly PM
2.5 concentrations were observed (
Figure 5a and
Table 2). As a general trend for the period 2019–2022, the mean monthly PM
2.5 values among the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10 varied greatly. The highest mean monthly PM
2.5 concentrations were observed at the AQMS S7 and S8. In March, April, and May of 2020, the mean monthly PM
2.5 concentrations were 19.8 μg/m
3, 16.5 μg/m
3, and 9.5 μg/m
3 for the AQMS S7 and 34.5 μg/m
3, 28.0 μg/m
3, and 11.2 μg/m
3 for the AQMS S8. For the corresponding period of 2019, the mean monthly PM
2.5 concentrations were generally higher. Specifically, the corresponding values in March, April, and May 2019 for the AQMS S7 were 22.6 μg/m
3, 18.1 μg/m
3, and 9.2, and for the AQMS S8 were 42.7 μg/m
3, 28.4 μg/m
3, and 15.5 μg/m
3.
As expected, there is a significant decrease in PM
2.5 concentrations in 2020 compared to 2019 (
Figure 5a and
Table 2). As for the AQMS S1, S2, S4, S5, S6, S9, and S10, the mean monthly PM
2.5 concentrations also decreased in March and April of 2020 compared to March and April of 2019 (
Table 2). However, for the majority of the AQMS, in May 2020, there was an increase in PM
2.5 concentrations in terms of monthly values compared to the previous year. Specifically, the mean monthly PM
2.5 concentrations in May 2020 were 10.9 μg/m
3, 9.0 μg/m
3, 10.2 μg/m
3, 9.0 μg/m
3, 10.0 μg/m
3, 9.5 μg/m
3, 11.2 μg/m
3, 7.3 μg/m
3, and 9.1 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In contrast, in May 2019, the mean monthly values were 8.1 μg/m
3, 6.9 μg/m
3, 8.3 μg/m
3, 7.1 μg/m
3, 7.0 μg/m
3, 9.2 μg/m
3, 15.5 μg/m
3, 7.9 μg/m
3, and 7.5 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. The slight increase in concentrations could be attributed to the fact that from 4 May 2020, after the lockdown period of March and April 2020, the restrictions on movement and business activity began to gradually lift.
As for PM10 concentrations, the mean monthly values exhibit spatial heterogeneity, and no seasonality patterns can be detected at all AQMS. Solely, AQMS S7 and S8 exhibit clear seasonality patterns, with the highest PM concentrations in the winter period and the lowest concentrations in the summer period. As mentioned above, PM concentrations in Meliti and Florina are highly dependent on heating systems. Additionally, high PM10 concentrations during the winter period are detected at the AQMS S4.
In general, research on air quality in Western Macedonia prior to the COVID-19 pandemic has also found that the highest PM
10 concentrations were registered in the warm period of the year [
13,
14,
15,
19,
20]. For instance, Triantafyllou (2000) studied the PM
10 concentrations in the southern part of the Eordea Basin from January 1991 to December 1994 and found that PM
10 concentrations are higher during the summer and early autumn and lower during the spring [
13]. This pattern in PM
10 concentrations could be attributed to multiple causes, including seasonal changes in the atmospheric dispersion characteristics and the absence of scavenging by precipitation [
13,
14,
15,
19,
20]. Additionally, wind-induced resuspension is a secondary source of PM
10 in the open-pit coal mines during the warm period [
19,
20].
As we can see from
Figure 5b and
Table 3 at all the AQMS, the mean monthly PM
10 concentrations decreased in March and April of 2020 compared to March and April of 2019. In detail, the mean monthly PM
10 concentrations in March 2020 were 19.2 μg/m
3, 18.1 μg/m
3, 22.3 μg/m
3, 17.5 μg/m
3, 18.6 μg/m
3, 24.6 μg/m
3, 38.6 μg/m
3, 15.5 μg/m
3, and 17.3 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In April 2020, the mean monthly PM
10 concentrations were 18.7 μg/m
3, 16.2 μg/m
3, 19.3 μg/m
3, 14.9 μg/m
3, 17.7 μg/m
3, 21.3 μg/m
3, 33.0 μg/m
3, 14.3 μg/m
3, and 17.1 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. As for March and April of 2019, the mean monthly PM
10 concentrations were considerably higher. The mean monthly PM
10 in March 2019 was 26.2 μg/m
3, 22.5 μg/m
3, 28.2 μg/m
3, 20.7 μg/m
3, 21.4 μg/m
3, 34.1 μg/m
3, 55.0 μg/m
3, 25.7 μg/m
3, and 23.3 μg/m
3, and in April 2019 was 24.7 μg/m
3, 23.3 μg/m
3, 27.9 μg/m
3, 20.3 μg/m
3, 20.4 μg/m
3, 28.5 μg/m
3, 38.6 μg/m
3, 22.4 μg/m
3, and 24.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively.
On the other hand, at all the AQMS, the mean monthly PM
10 concentrations in May 2020 were higher compared to May 2019 (
Table 3). In May 2019, the mean monthly PM
10 concentrations were 15.3 μg/m
3, 11.7 μg/m
3, 15.0 μg/m
3, 12.0 μg/m
3, 10.9 μg/m
3, 13.7 μg/m
3, 23.4 μg/m
3, 18.1 μg/m
3, and 17.0 μg/m
3, while in May 2020, the mean monthly PM
10 concentrations were 24.2 μg/m
3, 22.7 μg/m
3, 23.1 μg/m
3, 24.4 μg/m
3, 23.5 μg/m
3, 22.7 μg/m
3, 19.9 μg/m
3, 20.4 μg/m
3, and 24.5 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively.
As for the years 2021 and 2022, it seems that the concentrations of PM remained at the levels of 2020 (
Figure 3,
Figure 4 and
Figure 5a,b and
Table 2 and
Table 3). This could not be attributed merely to the COVID-19 pandemic but also to other factors. In the period between the end of 2019 and mid-2020, the units I and II of the power station Kardia and the units I and II of the power station Amyntaio ceased to operate. The closure of these units had an obvious effect on the improvement of air quality in the region in the last decade, as reported in prior studies [
19,
20]. Sachanidis et al. (2022) reported that the improvement in ambient air quality in the Western Macedonia Lignite Center is correlated with the reduction in excavated rock volumes and the lignite amount produced [
25]. Importantly, they highlighted that the overall better air quality is measured higher in terms of the number of exceedances of PM limit values. The PM exceedances are mainly correlated with air pollution episodes attributed to lignite mining activities, prevailing meteorological conditions, and long-range dust transport seasonal phenomena [
25]. The emissions from the combustion of lignite in the power stations, the mining activities of the lignite coal, and the transport of fugitive dust sources and fly ash cause local air pollution phenomena [
14]. Matthaios et al. (2017) investigated the occurrence of extreme PM
10 pollution episodes in Greece, including the region of Western Macedonia, and found that the local sources contribution to PM
10 reached up to 64%–74%, which is the highest among the regions under study [
37].
As for the improvement of air quality during the lockdown in Greece, previous studies have also reported a decrease in air pollutant concentrations, but additional factors also influence air pollutant concentrations. Kotsiou et al. (2021) examined air quality during the pandemic in Volos, a coastal port city in Greece with almost 86,000 inhabitants, in accordance with the national census data of 2021 [
38]. They found that the lockdown resulted in a 37.4% reduction in mean daily PM
2.5 in 2020 compared to 2019 levels, but during the strictest lockdown (from the 23 March to the 4 May), the occurrence of high levels of PM
2.5 was not avoided, even though there were restrictions in human activity patterns [
38].
In our study, the period of the first national lockdown that began on the 23 March 2020 and ended in the first days of May 2020 is shown in
Figure 6 and
Figure 7. More specifically,
Figure 6a shows the daily mean PM
2.5 concentrations for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10 for the period from the 1 February 2019 to the 1 June 2019, and
Figure 6b shows the corresponding period in 2020. Similarly,
Figure 7a,b shows the daily mean PM
10 concentrations for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10 for the above-mentioned periods. In
Table 4, the average PM
2.5 and PM
10 concentrations (μg/m
3) over the lockdown period (23/3/2020–04/05/2020) and the corresponding period of the previous year (23/3/2019–04/05/2019) are shown.
During the first lockdown period (from the 23 March to the 4 May 2020) the daily PM
2.5 concentrations decreased. Reductions in daily mean PM
2.5 concentrations were observed at all AQMS in the weeks before the lockdown and during the 6-week lockdown (from the 23 March to the 4 May 2020) (
Figure 6b). During the 6-week lockdown, the mean PM
2.5 concentrations were 13.1 μg/m
3, 11.6 μg/m
3, 15.0 μg/m
3, 10.7 μg/m
3, 12.9 μg/m
3, 17.1 μg/m
3, 29.8 μg/m
3, 8.6 μg/m
3, and 10.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In this period, there was a decrease relative to 2019 of −16%, −27%, −21%, −25, −9%, −1%, −3%, −36%, and −19% in the average PM
2.5 concentrations for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. During the corresponding period in 2019, the mean PM
2.5 concentrations were 15.6 μg/m
3, 15.8 μg/m
3, 19.1 μg/m
3, 14.3 μg/m
3, 14.2 μg/m
3, 17.2 μg/m
3, 30.9 μg/m
3, 13.5 μg/m
3, and 12.4 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively.
Specifically, during the period from the 23 March 2019 to the 4 May 2019, the maximum daily PM2.5 concentrations were above 25 μg/m3, while the daily PM2.5 concentrations reached up to maximum daily values of 55 μg/m3 at the AQMS S8. In 2020, high levels of PM2.5 concentrations were also observed, with the daily maximum PM2.5 reaching up to 44.1 μg/m3 and 52.2 μg/m3 at the AQMS S7 and S8, respectively.
As we can see in
Figure 7a,b, the daily mean PM
10 concentrations during the first lockdown period (from the 23 March to the 4 May 2020) also decreased. During the 6-week lockdown, the average PM
10 concentrations were 18.7 μg/m
3, 16.9 μg/m
3, 20.0 μg/m
3, 15.6 μg/m
3, 18.2 μg/m
3, 21.9 μg/m
3, 34.4 μg/m
3, 14.5 μg/m
3, and 16.8 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. In this 6-week period, there was a decrease relative to 2019 of −25%, −28%, −30%, −26, −13%, −18%, −17%, −40%, and −30% in the average PM
10 concentrations for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. During the corresponding period in 2019, the mean PM
10 concentrations were 24.9 μg/m
3, 23.4 μg/m
3, 28.4 μg/m
3, 21.2 μg/m
3, 20.9 μg/m
3, 26.6 μg/m
3, 41.6 μg/m
3, 24.3 μg/m
3, and 24.0 μg/m
3 for the AQMS S1, S2, S4, S5, S6, S7, S8, S9, and S10, respectively. However, it is worth mentioning that in April 2019, there was a 5-day period of extremely high PM
10 values with spatial homogeneity (
Figure 7a). From the 24 April to the 27 April 2019, the PM
10 concentrations gradually increased at all the AQMS under study. This event, which occurred in the second half of April 2019, is attributed to a large-scale Saharan dust episode over Greece and Europe [
39]. In Greece, dust episodes are generally common phenomena during spring and autumn and mainly affect the southern part of the country. Similarly, a long-lasting Saharan dust episode occurred in Greece from the 14 May to the 20 May 2020, with high values of PM
10 concentrations at all the AQMS under study (
Figure 6b) [
40]. As has been mentioned above, long-range transport is an important mechanism for particle pollution episodes [
8].
3.2. PM Concentrations in Correlation with Meteorological Parameters and Air Pollutants
Figure 8a,b shows the mean monthly variations in daily PM
2.5, PM
10, SO
2, NO
2, NO, and NO
X concentrations (μg/m
3) along with the mean air temperature (°C) for the AQMS S6 and S7 for the years 2019 and 2020. AQMS S6 is located in Amyntaio, which is a town with 4348 inhabitants, and AQMS S7 is located in Florina, which is a town with a population of 17,188 inhabitants [
41].
As we previously discussed, the AQMS S7 registered high PM
2.5 and PM
10 concentrations because it is highly affected by the local sources of air pollution from traditional heating systems. The highest air pollutant concentrations were observed during the winter months, when the lowest air temperatures were observed. For example, in January of 2019, the mean air temperature was below 0 °C (T
mean = −1.8° C), while PM
10 and PM
2.5 reached up to 82.6 μg/m
3 and 80.1 μg/m
3, respectively. The PM
2.5/PM
10 ratio reached a value of 1.0, indicating the major contribution of fine particles attributable to anthropogenic air pollution sources. Based on previous studies, the high ratios that have been found in Florina indicate the great contribution of individual household heating systems (e.g., fireplaces and woodstoves) and biomass burning [
20].
In March and April of 2020, the air temperature was lower compared to the same months in 2019. The monthly mean air temperature was 8.4 °C and 11.4 °C in March and April of 2020, respectively, and 10.4 °C and 12.2 °C for the same months in 2019, respectively (
Figure 8b). However, PM concentrations were lower in 2020 compared to 2019, indicating that lockdown restrictions influenced air pollutant concentrations. As we have previously mentioned, PM
10 concentrations in March and April of 2020 were considerably lower compared to the corresponding months in 2019.
At the AQMS S7, the mean monthly SO
2, NO, NO
2, and NO
X concentrations were 19.6 μg/m
3, 4.8 μg/m
3, 13.9 μg/m
3, and 18.7 μg/m
3, respectively, in March of 2019, and 7.5 μg/m
3, 12.9 μg/m
3, 8.8 μg/m
3, and 21.7 μg/m
3, respectively, in April of 2019. In 2020, the air pollutant concentrations were lower. In detail, in March 2020, the mean monthly SO
2, NO, NO
2, and NO
X concentrations were 8.3 μg/m
3, 3.5 μg/m
3, 5.5 μg/m
3, and 9.0 μg/m
3, respectively, while in April 2020 the concentrations were 8.6 μg/m
3, 3.2 μg/m
3, 3.9 μg/m
3, and 7.1 μg/m
3, respectively (
Figure 8b).
At the AQMS S6, the mean monthly SO
2, NO, NO
2, and NO
X concentrations were 7.2 μg/m
3, 1.6 μg/m
3, 51.4 μg/m
3, and 53.0 μg/m
3, respectively, in March 2019, and 3.2 μg/m
3, 1.40 μg/m
3, 29.9 μg/m
3, and 31.3 μg/m
3, respectively, in April 2019 (
Figure 8a). Specifically, the mean monthly concentrations decreased sharply from May 2019 onwards to average monthly SO
2, NO, NO
2, and NO
X concentrations of 3.3 µg/m
3, 1.6 µg/m
3, 4.0 µg/m
3, and 5.6 µg/m
3, respectively. The decreased concentrations are attributed to zero emissions of PM, SO
2, and NO
X from May 2020 and onwards due to the closure of PS4 (
Table 5).
As for the meteorological conditions at the AQMS S6, the air temperature is relatively low but at higher levels compared to the AQMS S7. The monthly mean air temperature in January of 2019 was also below 0 °C, while in March and April of 2019 it was at 10.3 °C and 11.8 °C, respectively. During the lockdown period in 2020, the monthly mean air temperature in March and April was 8.2 °C and 11.2 °C, respectively.
Figures S1a–l and S2a–l in the
Supplementary Material show the seasonal correlation matrices with Pearson correlation coefficients between daily air pollutant concentrations and meteorological parameters for the AQMS S6 and AQMS S7 for the years 2019, 2020, and 2021. Although there are differences among the Pearson correlation coefficients through the years and seasons, very strong positive correlations were found between PM
2.5 and PM
10, as expected. In general, differences in correlation coefficients were found between the AQMS S6 and AQMS S7.
At the AQMS S6, the Pearson correlation coefficients between PM
2.5 and PM
10 are strong and very strong positive in all seasons and through the years (
Figure S1a–l, Supplementary Material). For example, these correlation values were r = 0.937 (
p-value < 0.001) for spring, r = 0.766 (
p-value < 0.001) for summer, r = 0.874 (
p-value < 0.001) for autumn, and r = 0.924 (
p-value < 0.001) for winter in 2019. In 2020, the correlation coefficients between PM
2.5 and PM
10 are r = 0.777 (
p-value < 0.001) for spring, r = 0.836 (
p-value < 0.001) for summer, r = 0.850 (
p-value < 0.001) for autumn, and r = 0.947 (
p-value < 0.001) for winter. As expected, the r values are higher in winter due to similar emission sources (e.g., domestic heating and the increased electricity demand from the lignite-fired power plants). In the spring of 2020, when the lockdown was implemented, the r value was lower compared to 2019 because of the reduction in emissions during the COVID-19 restrictions. In 2021, the correlation coefficients between PM2.5 and PM10 were r = 0.796 (
p-value < 0.001) for spring, r = 0.849 (
p-value < 0.001) for summer, r = 0.844 (
p-value < 0.001) for autumn, and r = 0.927 (
p-value < 0.001) for winter.
Additionally, at the AQMS S6, the correlation coefficients of PM
2.5 and PM
10 with NO
X, NO
2, NO, and SO
2 varied among the seasons. Generally, the r values between the above-mentioned variables are moderate, ranging between 0.4 and 0.6 in spring, autumn, and winter, but in summer, the r values are also negative in certain cases. For example, the correlation coefficients of PM
2.5 and PM
10 with NO
2 are r = 0.549 (
p-value < 0.001) and r = 0.642 (
p-value < 0.001) for the winter of 2019, r = 0.494 (
p-value < 0.001) and r = 0.486 (
p-value < 0.001) for the winter of 2020, while for the winter of 2021 the correlation coefficients are r = −0.272 (
p-value = 0.011) and r = −0.248 (
p-value = 0.021), respectively. From these correlation coefficients, we can conclude that during the winter period, NO
2 and both PM
2.5 and PM
10 came from the same source, which is mainly the power plant PS4 (
Table 5), while in 2021, when the PS4 was not operated, the correlation was negative but not significant. During the lockdown (spring 2020), the Pearson correlation coefficients between PM
2.5 and PM
10 and NO
X, NO
2, NO, and SO
2 were moderate or even negative. For example, the PM
2.5 and PM
10 correlation coefficients with the NO
2 correlation coefficient were r = 0.078 (
p-value = 0.462) and r = 0.334 (
p-value = 0.001) for the spring of 2020, respectively. In contrast, in the spring of 2019, the correlation coefficient of PM
2.5 and NO
2 was r = 0.569 (
p-value < 0.001), and the correlation coefficient of PM
10 and NO
2 was r = 0.476 (
p-value < 0.001). In the spring of 2021, PM
2.5 and NO
2 were correlated with r = −0.388 (
p-value < 0.001), and PM
10 and NO
2 were correlated with r = −0.041 (
p-value = 0.701).
In accordance with the Pearson correlation coefficient, the meteorological factors affecting the air pollutants varied in different seasons at the AQMS S6. In general, the air temperature was negatively correlated with PM
2.5 and PM
10 in winter, but in summer, the correlation coefficients were positive. The low temperatures in winter lead to increased emission rates from domestic heating and electricity generation from power plants. In the summer, the emissions from the power plants decreased, as we can also see from
Table 5. For example, PM
2.5 and air temperature exhibited a negative correlation in the winters of 2019, 2020, and 2021, with r = −0.126 (
p-value = 0.259), r = −0.299 (
p-value = 0.004), and r = −0.162 (
p-value = 0.134), respectively. On the other hand, in the summer, the correlations were reversed. The correlation coefficient was r = 0.310 (
p-value = 0.003) in the summer of 2019, r = 0.459 (
p-value < 0.001) in the summer of 2020, and r = 0.565 (
p-value < 0.001) in the summer of 2021. Similar values were found for temperature and PM
10, with correlation coefficients of r = 0.574 (
p-value < 0.001), r = 0.677 (
p-value < 0.001), and r = 0.716 (
p-value < 0.001) for the summers of 2019, 2020, and 2021. Additionally, the Pearson correlation coefficients between relative humidity, wind speed, and wind direction and PM
2.5 and PM
10 exhibit differences among the seasons. Wind speed and direction are negatively or weakly correlated with PM
2.5 and PM
10 in all seasons. In general, negative correlations between wind speed and air pollutant concentrations indicate the horizonal dispersion and dilution of air pollutant concentrations under the influence of strong winds.
At the AQMS S7, very strong positive Pearson correlation coefficients were generally found for PM
2.5 and PM
10 in all seasons and years (
Figure S2a–l, Supplementary Material). Very strong correlation coefficients were found for the winter period. The r values for PM
2.5 and PM
10 were r = 0.986 (
p-value < 0.001) for the winter of 2019, r = 0.981 (
p-value < 0.001) for the winter of 2020, and r = 0.959 (
p-value < 0.001) for the winter of 2021. Additionally, spring, summer, and autumn exhibited a high degree of correlation, as indicated by the r values ranging from 0.7 to 0.9. The correlation coefficient between PM
2.5 and PM
10 was r = 0.925 (
p-value < 0.001) in the spring of 2019, r = 0.668 (
p-value < 0.001) in the spring of 2020, and r = 0.782 (
p-value < 0.001) in the spring of 2021. A previous study also found a high coefficient (R
2) between PM
2.5 and PM
10 equal to 0.92 during the winter period in Florina (AQMS S7), which suggests that PM
2.5 and PM
10 came from similar emission sources [
20]. Similar to the AQMS S6, the correlation coefficients between PM
2.5 and PM
10 decreased during the months of lockdown in the spring of 2020.
As for the correlation coefficients of PM2.5 and PM10 with NOX, NO2, NO, and SO2 at the AQMS S7, the r values varied among the seasons. In general, in winter, the correlation coefficients for PM2.5 with NO2 are strong or even very strong, with values reaching up to 0.9, while relatively low or even negative values are found in summer. However, there are no major differences before, during, or after COVID-19 lockdown for the correlation coefficients of PM2.5 and PM10 with NOX, NO2, and NO. The Pearson correlation coefficients of PM2.5 with NO2 were r = 0.647 (p-value < 0.001) in the spring of 2019, r = 0.626 (p-value < 0.001) in the spring of 2020, and r = 0.848 (p-value < 0.001) in the spring of 2021. In general, it is expected that PM2.5 and NO2 exhibit strong correlations because both PM2.5 and NO2 are generated during the combustion process with household stoves, space heaters, furnaces, fireplaces, and boilers.
As for the meteorological parameters at the AQMS S7, strong but negative correlation coefficients were found for temperature with PM2.5 and PM10 during all seasons except summer. For example, PM2.5 and air temperature exhibited a negative correlation in the winters of 2019, 2020, and 2021, with r = −0.429 (p-value < 0.001), r = −0.396 (p-value < 0.001), and r = −0.401 (p-value < 0.001), respectively. On the other hand, in the summer, their correlations were reversed. The correlation coefficients were r = 0.195 (p-value = 0.062), r = 0.385 (p-value < 0.001), and r = 0.435 (p-value < 0.001) in the summer of 2019, 2020, and 2021, respectively. Similar values were found for temperature and PM10, with correlation coefficients r = 0.446 (p-value < 0.001), r = 0.564 (p-value < 0.001), and 0.612 (p-value < 0.001) for the summers of 2019, 2020, and 2021, while in the winters of 2019, 2020, and 2021, the correlation coefficients were r = −0.357 (p-value < 0.001), r = −0.356 (p-value < 0.001), and r = −0.264 (p-value = 0.013), respectively. In general, the Pearson correlation coefficients at the AQMS S7 also explain the high dependence of PM on meteorological conditions. The negative correlation between temperature and PM at the AQMS S7 indicates that lower temperatures led to higher PM concentrations. Similar to the AQMS S6, at the AQMS S7, the Pearson correlation coefficients between relative humidity, wind speed, and wind direction and PM2.5 and PM10 exhibit differences among the seasons. Wind speed and direction are negatively or weakly correlated with PM2.5 and PM10 in all seasons.
Previous studies have also used Pearson correlation coefficients to investigate the effects of meteorological parameters on air pollution during the COVID-19 pandemic [
43,
44]. In Italy, Cucciniello et al. (2022) evaluated the air quality during lockdown in the city of Avellino [
44]. The climatic conditions in the region favour the usage of domestic heating systems (boilers, fireplaces, and pellet stoves) during the winter period up to the end of April, which is comparable with Western Macedonia. The Pearson correlations of the atmospheric pollutants’ concentrations in Avellino also showed strong associations between PM
2.5 and PM
10, while a slight decrease was detected in correlations for NO
2 and both PM
10 and PM
2.5 during the lockdown. In accordance with Cucciniello et al. (2022), this could be attributed to the pollutants’ source (e.g., the decreased vehicular traffic during the lockdown). In Poland, Górka-Kostrubiec and Dudzisz (2023) analysed the effect of pandemic restrictions on air pollution in Warsaw and Krakow by examining the correlations of the average monthly concentrations of pollutants in particular months of 2019, 2020, and 2021 [
45]. The significant correlation coefficients for PM
2.5 and PM
10 suggest the same source (e.g., central heating system) of particle pollution, while in Warsaw, the concentrations of PM
2.5 and PM
10 with NO
X seemed to be less correlated during the pandemic restrictions. The difference in correlation coefficients between the cities may be attributed to the fact that in Kracow, the dominant stack emissions from heating systems that use mainly coal and wood during the winter months may have masked the reduction in NO
X during lockdowns [
45]. This is comparable with our study, where the local sources of air pollution from the traditional heating systems dominate in Florina. In Victoria, Mexico, Tello-Leal and Macías-Hernández (2021) used the Pearson correlation analysis to estimate the relationships between the concentrations of air pollutants and meteorological variables during the lockdown and found very strong positive correlations between PM factions (r =0.99,
p-value < 0.001) and consistently moderate to very strong negative correlations (−0.92 ≥ r ≥ −0.45) between temperature and all the air pollution variables [
43]. In Changchun, China, there are many industrial sources that are common sources for air pollutants [
46]. So, there are positive correlations between PM
2.5, PM
10, SO
2, NO
2, and CO before and after lockdown, while the negative correlation between temperature and the above-mentioned pollutants before lockdown and during the first phase of lockdown (25 January–25 February 2020) could be attributed to low temperatures that increase the emission sources (such as heating in coal-fired power plants).
In Volos, Greece, Kotsiou et al. (2021) found that PM
2.5 concentrations were negatively correlated with temperature (r = −0.47,
p = 0.001) and positively correlated with humidity (r = 0.37,
p = 0.011) during lockdown [
38]. Moreover, in Athens, Greece, Grivas et al. (2020) compared the observed changes in pollutant levels during the lockdown (23 March 2020–10 May 2020) to changes and the potential impact of prevailing weather conditions and concluded that the lockdown was the dominant factor for this drop in average NO
2 levels [
47]. This contribution is estimated to reach 65%, while the remaining 35% of the mean reduction in NO
2 levels can be attributed to the different meteorological conditions during the lockdown when compared to the pre-lockdown period [
47]. In Athens, during the lockdown period in March, vehicular traffic was significantly reduced (almost 50%) compared to the same period in 2019. Similarly, in London, UK, it was reported that there was a decrease of around 32% in total traffic and an overall 22.5% reduction in traffic during the 15 months of social restrictions [
48].
Generally, the traffic reduction correlated with decreased air pollutant concentrations and particularly strong reductions of NO
2 emissions and concentrations, while an upward trend during the lockdown has been observed for O
3. This inverse relationship between O
3 and NO
2 values has been detected in various studies [
8,
9,
48]. NO
X is a major precursor of O
3 as well as a quencher of O
3 through NOx titration. High concentrations of NO locally scavenge O
3 and lead to the formation of NO
2, and high concentrations of NO
2 block the oxidation step of VOCs by producing nitric acid, which prevents the net formation of O
3 [
9]. In parallel, the prevailing meteorological conditions (local winds, air temperature, and humidity) affect these chemical processes that result in an increase in O
3 because of a decrease in NO
X and low VOC/NO
X ratios. For instance, in the UK, the unseasonably warm start and end to 2020 and the reduction in NO
X caused increased O
3 production [
48]. In the highly polluted metropolitan cities in India, a remarkable drop in the mean concentration of the AQI (Air Quality Index) was detected during the COVID-19 lockdown, while O
3 concentrations significantly increased as the NO
X concentrations increased [
49]. Importantly, the inverse correlation between O
3 and NO
X is evident through the Pearson correlation coefficient, given that O
3 and AQI were negatively correlated [
49].
In our study, NO
X, NO
2, and NO were negatively correlated with air temperature during winter, spring, and autumn. For example, at the AQMS S7, the correlation coefficients for NO
2 and air temperature were r = −0.373 (
p-value < 0.001) for the spring of 2020, r = 0.092 (
p-value = 0.380) for the summer of 2020, r = −0.727 (
p-value < 0.001) for the autumn of 2020, and r = −0.021 (
p-value = 0.847) for the winter of 2020. The corresponding correlation coefficients for NO
2 and air temperature before the COVID-19 pandemic were r = −0.482 (
p-value < 0.001) for the spring of 2019, r = 0.119 (
p-value = 0.260) for the summer of 2019, r = −0.350 (
p-value < 0.001) for the autumn of 2019, and r = −0.383 (
p-value < 0.001) for the winter of 2019. Lower NO
2 concentrations usually occur at higher temperatures because the photochemical reaction of NO
2 and the vertical dispersion increase at higher temperatures [
50]. This leads to higher O
3 concentrations. High NO
2 concentrations generally occur in the winter when the temperature is lower and the emissions from heating systems increase, as in our study.
Furthermore, many studies have analysed the association of air pollution and meteorological variables with the incidence of the COVID-19 pandemic, while highly polluted regions seem to be correlated with COVID-19 cases [
51,
52]. This hypothesis arose in the early stages of COVID-19’s emergence [
52]. Western Macedonia was highly affected by COVID-19 in the early stages of the pandemic. So, this study provides an opportunity to further examine the relationships between air pollution and meteorology in the context of human health in the region.