Next Article in Journal
E-Cigarette Use among Male Smokers in Al-Ahsa, Kingdom of Saudi Arabia: A Cross-Sectional Study
Previous Article in Journal
Towards an Evidence-Based Model of Workplace Postvention
Previous Article in Special Issue
The Comprehensive Health Risk Assessment of Polish Smelters with Ecotoxicological Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Volatility of PM2.5, PM10 and PM10-Bound B[a]P Concentrations and Assessment of the Exposure of the Population of Silesia in 2018–2021

Department of Air Protection, Silesian University of Technology, 22B Konarskiego St., 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(1), 138; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20010138
Submission received: 23 November 2022 / Revised: 15 December 2022 / Accepted: 19 December 2022 / Published: 22 December 2022

Abstract

:
Air pollution both indoors and outdoors is a major cause of various diseases and premature deaths. Negative health effects are more frequently observed in a number of European countries characterized by significant pollution. In Poland, especially in Upper Silesia, the most serious problem is the high concentration of particulate matter (PM) and PM10-bound benzo[a]pyrene (B[a]P). The main source of these two pollutants is so-called “low emissions” associated with the burning of solid fuels mainly in domestic boilers and liquid fuels in road traffic. This study examined the variability in the PM and PM10-bound B[a]P concentrations and their relationships with meteorological parameters, i.e., atmospheric pressure, air temperature and wind speed, in 2018–2021 at 11 monitoring stations. In many Silesian cities, the average annual concentrations of PM10, PM2.5 and B[a]P were much higher than those recorded in other European countries. At each station, the average daily PM10 concentrations were exceeded on 12 to 126 days a year. Taking into account the WHO recommendation for PM2.5, the highest recorded average daily concentration exceeded the permissible level by almost 40 times. The same relationships were observed in all measurement years: PM10 concentrations were negatively correlated with air temperature (R = −0.386) and wind speed (R = −0.614). The highest concentrations were observed in the temperature range from −15 °C to −5 °C, when the wind speed did not exceed 0.5 m·s−1. The calculated lifetime cancer risk (LCR) associated with the exposure to B[a]P in the Silesian Voivodeship suggested 30–429 cases per 1 million people in the heating season depending on the scenario used for the calculations (IRIS, EPA or WHO).

1. Introduction

Air pollution occurring indoors and outdoors is the most serious environmental health threat in the modern world [1]. The pollutants with the most serious impacts on human health are particulate matter, nitrogen dioxide and ground-level ozone. Their synergistic effect is particularly dangerous [2,3]. They are the main causes of illness and premature death (in 2016, it was estimated that air pollution was the cause of 4.2 million premature deaths [1]) caused by cardiovascular diseases, i.e., heart diseases and strokes, as well as respiratory diseases, predominantly lung diseases and lung cancer [4,5]. In the countries of the European Union (EU), approx. 300,000 premature deaths caused by fine particulate matter were recorded for the year of 2019 [6,7]. It is worth emphasizing that the onset of symptoms of the diseases caused by air pollution is delayed after exposure.
Research shows that the environmental burden of disease varies across Europe. In areas with significant environmental pollution, negative health effects are seen more often. Exceedances of the daily limit value recommended by the EU for PM10 (50 µg·m−3) occur throughout the continent; however, the highest concentrations are observed in Northern Italy, Croatia, Bulgaria, Serbia, Kosovo, Turkey, Bosnia and Herzegovina, North Macedonia and Poland [8,9,10,11]. Unfortunately, Poland is considered to be a country with some of the severest problems related to air pollution in the EU, mainly related to PM2.5, PM10 and benzo[a]pyrene pollutants. Both particulate matter and benzo[a]pyrene are classified by the International Agency for Research on Cancer (IARC) as carcinogenic to humans (Group 1) [12]. It should be noted that, for many years, B[a]P has been recognized as a good and sufficient indicator of human exposure to priority polycyclic aromatic hydrocarbons (PAHs) present in the air [13]. Studies conducted on animals confirmed that B[a]P can cause various types of neoplasms, such as gastrointestinal tract, liver, kidney, respiratory tract, pharynx and skin cancers [14]. In many countries, including Poland, the PM10-bound B[a]P concentration in the air is routinely monitored in order to ensure that its average annual concentration does not exceed the limit value (the established limit that the average annual concentration of B[a]P should not exceed is 1 ng·m−3 in most areas [15], and, according to the WHO recommendations, it is 0.12 ng·m−3 [16]). Over the last few decades, the interest in exposure to PM10-bound B[a]P and PM2.5-bound B[a]P has increased because they are associated with a broad range of health effects that have a major impact on public health.
In Poland, especially in Silesia, the highest share of the emissions of PM10, PM2.5 and B[a]P (77%, 87.9% and 97.8%, respectively) is attributed to communal and household sources [17]. Consequently, municipal emissions attributed to the combustion of fuels in domestic boilers (mainly coal and wood and, despite the ban [18], municipal waste [19,20,21]) are considered the main cause of exceeding the concentration limits of PM and B[a]P in the air in Poland. The evidence for this is that the PM and B[a]P concentrations are higher in the heating season than in the non-heating season [22,23,24]. According to the most recent available data, in 2018, almost half of the households in Poland, i.e., 45.4%, still used solid-fuel heating devices [25]. Other important sources of emissions include industrial activities [26], as well as agricultural and road transport, especially in urban areas [27,28]. According to the ranking prepared by the World Health Organization (WHO), 36 out of 50 cities with the worst air quality in the EU are located in Poland, most of them in Upper Silesia [29].
During the heating season in Silesia, so-called smog incidents [30] are also more frequent, during which the concentration of pollutants increases above the permissible standards. The formation of smog depends on the coexistence of a number of meteorological factors, i.e., low temperature, a low wind speed, thermal inversion and the height of the mixing layer. The occurrence of smog in Poland is primarily associated with excessive air pollution, especially the particulate matter of large particles (PM10 and PM5), the particulate matter of medium and small particles (PM2.5 and PM1) and polycyclic aromatic hydrocarbons (PAHs), and, to a lesser extent, sulfur oxides, nitrogen oxides and ozone or carbon monoxide. According to a recent study [31], Poland can be divided into three smog regions (Figure 1): Region I comprises the south and a part of Silesia, where smog incidents occur frequently and last the longest. This region includes cities and poorly ventilated areas, where a stable atmosphere with temperature inversion occurs frequently. Region II is a strip of land stretching from the Lublin Upland to the Wielkopolska Lowland, with a lower frequency of smog incidents. Finally, Region III comprises northern Poland, where smog incidents are the least frequent, which is attributed to the unstable atmosphere dominant in this area. Analyses of smog incidents in Silesia have shown significant increases in the number of registered patients for ambulatory care due to bronchitis and asthma exacerbation in response to increases in PM2.5 concentrations [32,33]. According to [32], the described dependencies are visible during the first few days of smog episodes, and they do not cease in the following days. In the case of hospitalization due to acute respiratory diseases, a greater number of registered events appeared with a two-day delay from the observation of a smog incident. During the periods of smog incidents, a significant increase in mortality was also observed [32].
Although various legal regulations and programs have been introduced in Poland [34,35,36], in the era of the energy crisis and increasing energy poverty [37,38], we can expect a deterioration in air quality in the near future. It should also be remembered that the quality of outdoor air also has a significant impact on the quality of indoor air [39], especially in places where natural ventilation is used. It has been observed that higher concentrations of pollutants characteristic of outdoor air occur inside buildings where windows and doors face directly onto streets, main communication routes or car parks [40,41]. The aim of this study was to examine the correlation between the concentrations of PM10, PM2.5, B[a]P and meteorological parameters and to analyze the risk of developing cancer resulting from inhalation exposure to B[a]P among the inhabitants of Upper Silesia (Poland).

2. Materials and Methods

2.1. Study Area

Silesian Voivodeship is a province located in the southern part of Poland, and it occupies 12,333 km2, with a population of about 4.47 million (363 inhabitants per one square kilometer [42]). It is one of the most urbanized and industrialized regions of Central Europe. There are 30 monitoring stations located in this region [17]. In the frame of this research, the daily concentrations of particulate matter PM2.5, PM10 and PM10-bound B[a]P were analyzed in all monitoring stations; however, only 11 localizations (Figure 2), in which PM and B[a]P were tested simultaneously in a given time period, were chosen for further analyses (Table S1). The locations of all stations fulfilled the conditions required for the so-called urban background site (UB). PM2.5 was initially only measured at 3 measuring stations, and from 2021, it was measured at 5 stations. Only those stations for which a measurement coverage of 80% was ensured each year were included in further analyses. The heating season and the non-heating season were analyzed separately. It was assumed that the heating season lasts from the beginning of October to the end of March each year. As a criterion for a particular day qualifying as a smog incident day, the exceedance of the average daily PM10 concentration of 50 μg·m−3 was assumed, in accordance with the literature [43] (the 24 h threshold value).

2.2. Sample Collection

The available data concerning 24 h PM2.5, PM10 and PM10-bound B[a]P ambient concentrations in the Silesian Voivodeship from the period of 2018–2021 were taken from the Chief Inspectorate of Environmental Protection (CIEP) website [44].
At the selected monitoring stations, suspended particulate matter was manually sampled for 24 h a day, using high-volume samplers with an air flow of 30 m3·h−1. Samplers were equipped with a separate PM2.5 or PM10 head and a quartz fiber filter container, introducing subsequent filters into the measurement duct automatically. The sampled dust mass was determined with the use of the reference gravimetric method according to the PN-EN 12341:2014-07 standard [45]. The detection limit of this method is 1 μg·m−3.
PM10-bound B[a]P was extracted with an organic solvent and then analyzed using high-performance liquid chromatography (HPLC) with a fluorescence detector or gas chromatography with a mass spectrometry detector (GC-MS) in accordance with the PN-EN 15549:2011 standard [46]. These methods are applicable to the measurement of B[a]P in the concentration range of 0.04 to 20 ng·m−3, and their detection limits are below 0.04 ng·m−3 [46].

2.3. Statistical Analysis

In order to present a concise characterization of the analyzed concentrations, the basic parameters of descriptive statistics were calculated for each station, i.e., mean, median, standard deviation, maximum and minimum values, and the 1st and 3rd quartiles. The Spearman’s rank correlation coefficient was used to analyze the data on PM10 and B[a]P concentrations due to the fact that the criterion of similarity between the distribution of the examined variables and the normal distribution was not met. The similarity of the distribution of the analyzed variables to the normal distribution was examined using the Shapiro–Wilk test. All analyses were performed using Statistica 13 software, with a significance level of 0.05.

2.4. Lifetime Cancer Risk Assessment

Benzo[a]pyrene is recognized as a marker of the carcinogenic potential of a mixture of polycyclic aromatic hydrocarbons [47]; for this reason, it was taken as the indicator of PAH carcinogenic risk in this study similarly to other research [48,49].
Carcinogenic risk is often associated with the toxicity equivalence factor (TEF) of each of the 16 PAHs compared to B[a]P [50,51,52,53,54]. The B[a]P international toxic equivalent (TEQ) concentration is calculated by multiplying the concentration of each PAH with its corresponding TEF. As this study focused on B[a]P, it only took the TEQ of B[a]P with its TEF as a unity into consideration. Furthermore, the inhalation unit risk for benzo[a]pyrene was derived with the intention that it would be paired with EPA’s relative potency factors for the assessment of the carcinogenicity of the PAH mixtures [55].
The lifetime cancer risk (LCR) from the inhalation exposure to B[a]P can be expressed as follows:
LCR = TEQB[a]P × IURB[a]P
where:
TEQB[a]P—B[a]P equivalent equal to the B[a]P concentration in this study, µg·m−3;
IURB[a]P—inhalation cancer unit risk factor of exposure to B[a]P, IUR = 6 × 10−4 (µg·m−3)−1.
This exposure duration scenario included full lifetime exposure (assuming a 70-year lifespan). The inhalation unit risk (IUR) is defined as the upper-bound excess lifetime cancer risk estimated to result from continuous exposure to an agent at a concentration of 1 µg·m−3 in the air [56]. The IURB[a]P values for the carcinogenic benzo[a]pyrene used in this study were extracted from the database provided by the Integrated Risk Information System (IRIS) [57,58]. Previously, according to the Office of Environmental Health Hazards Assessment of the California Environmental Protection Agency, IURB[a]P was taken as 1.1 × 10−6 (ng·m−3)−1 [59,60]; however, the WHO suggests an IURB[a]P value of 8.7 × 10−6 (ng·m−3)−1 [61].

3. Results and Discussion

Excessive air pollution in Poland has been present for decades, especially in Upper Silesia. Studies of total particulate matter, conducted as early as the 1970s, showed that, in most cities in this area, the average concentrations of TSP exceeded 300 μg·m−3 annually, and the B[a]P concentrations were higher than 200 ng·m−3 [62,63]. The currently monitored concentrations of PM2.5, PM10 and B[a]P have decreased compared to the last century, but they are still high [44]. The highest emissions of dust pollutants from particularly burdensome industrial plants in Poland have been recorded in the area of the Silesian Voivodeship for many years. In 2020, these amounted to 4.4 thousand tone in total [64].

3.1. Concentrations of PM2.5 and PM10

The PM2.5 and PM10 concentrations measured in the period of 2018–2021 at 11 monitoring stations located in the Silesian Voivodeship were characterized by some variability (Table S2). The average annual concentrations, depending on the location and year, ranged as follows: from 19.62 to 38.55 μg·m−3 for PM2.5 and from 22.87 to 54.88 µg·m−3 for particulate matter PM10. These concentrations are much higher than the concentrations recorded in other European countries [29,65]. This means that, in many Silesian cities, the average annual concentrations exceeded the Polish standards of permissible concentrations (20 μg·m−3 for PM2.5 and 40 μg·m−3 for PM10). In comparison with the WHO guidelines (5 μg·m−3 for PM2.5 and 15 μg·m−3 for PM10), the average annual concentrations exceeded the safe level at every station. The highest concentrations were observed in the cities of Rybnik, Żywiec and Godów, where the share of the so-called low emissions in the overall balance is high [66].
In relation to the permissible limits regulated by Polish law and recommended by the WHO, the average daily concentrations of PM2.5 and PM10 were exceeded at almost every station (Figure 3 and Figure 4). According to the standards, the average daily concentration of PM10 may be exceeded for only 35 days a year; however, at most of the measuring stations, the daily average concentrations were exceeded from 12 to 126 days a year. A great number of these exceedances (approx. 73% of the total number of days) took place in the heating season (Tables S3 and S4). The highest concentration of PM2.5 was noted in the city of Godów in 2019, where it amounted to 228.60 μg·m−3, and the recommended level was exceeded almost 40 times according to the WHO guidelines (15 μg·m−3). Comparable daily concentration levels have only been recorded in some European cities (Sarajevo [67] and Belgrade [68]). These high concentrations and numerous exceedances of permissible limits are the cause of the frequent smog incidents in Silesia.
Due to the negative health effects of PM2.5, the ratio between the mass concentrations of PM2.5 and PM10, i.e., PM2.5/PM10, was additionally calculated (Table 1). For the entire calendar year, the PM2.5/PM10 ratio ranged from 0.6533 to 0.8804. In the heating and the non-heating seasons, the PM2.5/PM10 ratios were very similar, ranging from 0.6573 to 0.9877 and from 0.6451 to 0.7557, respectively, and they are similar to those in the literature [69,70]. For example, in Melpitz, located in eastern Germany, the average PM2.5/PM10 ratios were 0.72 in the summer and 0.82 in the winter [71]. The maximum difference between the heating and non-heating seasons was around 0.15. The calculations show that, in Silesia, the fine-grained PM fraction clearly predominates over the course-grained PM fraction. The increase in the emissions of fine and sub-micron particles is driven by the intensification of various activities, such as home heating based on coal fuel, industrial fuel combustion processes and the low mixing height of the air, resulting in the accumulation of secondary organic aerosol precursors and re-emissions [70,71].

3.2. Relationship between Some Meteorological Parameters and Particulate Matter Concentrations

The PM10 concentration depends on many environmental and anthropogenic factors. Meteorological conditions, i.e., wind direction and speed, air temperature, atmospheric pressure, precipitation, humidity, solar radiation, atmospheric stability and the height of the mixing layer, are considered to be factors of key importance [72,73,74]. They affect the diffusion, deposition and dilution of PM and, thus, the spatial distribution of pollutants in the bottom part of the atmosphere. In this paper, the authors analyzed selected meteorological parameters (i.e., wind speed, air temperature and atmospheric pressure) that were measured at the station in Rybnik. They were juxtaposed with the high concentrations that occurred regularly during the heating season at this station.
For this purpose, Spearman’s correlation coefficients were calculated (Table 2). The same dependencies were observed in all analyzed years. The PM10 concentrations were negatively correlated with air temperature (R = −0.386). Higher air temperatures lead to an effective vertical dispersion of pollutants, which results in low concentrations of particulate matter [72,75]. At low air temperatures, the activity of heating systems increases, and the height of the planetary boundary layer (PBL) [76] decreases, which causes a sharp increase in PM10 concentrations [77,78]. The speed of the wind was shown to be a very important factor, as it was strongly and negatively correlated with PM10 (R = −0.614). The highest concentrations of PM10 occurred in windless or almost windless conditions, which favors the formation of temperature inversion in the PBL [76], and this may cause an additional accumulation of pollutants in the ground level of the troposphere [79].
The highest concentration of PM10 was observed in a temperature range from −15 °C to −5 °C, when the wind speed did not exceed 0.5 m·s−1 (Figure 5). Similar results were observed in the literature [72,73,78,80]. The lowest average temperature (4.23 °C) in the heating season was recorded in Rybnik in 2018. Moreover, in the period from March 1 to March 6, the daily concentrations of PM10 and B[a]P remained at very high levels, exceeding the alarm level (150 µg·m−3 for PM10 according to [81]). Cyclic high concentrations occurred during every period characterized by very low temperatures and very low wind speeds (Figure 6).

3.3. Concentration of PM10-Bound B[a]P

In the Silesian Voivodeship in 2018–2021, the average concentration of PM10-bound B[a]P ranged from 6.2 to 7.1 ng·m−3 (Table S2). The highest B[a]P concentrations were recorded in Rybnik; in the 2018, 2019, 2020 and 2021 heating seasons, the concentrations were 23.5, 27.7, 15.9 and 15.7 ng·m−3, and in the non-heating seasons, the concentrations were 0.8, 2.9, 2.2 and 1.3 ng·m−3, respectively. In Częstochowa, the average annual concentrations of B[a]P in the studied period were the lowest, but they were still three times higher than the permissible values (EU). Despite this, the recorded concentration of B[a]P in Częstochowa in 2018 was the lowest in the entire voivodship, but the permissible concentrations were still exceeded as many as 20.6 times according to the WHO recommendation level of 0.12 ng·m−3. The observed PM10-bound B[a]P levels showed a clear seasonal pattern (Figure 7), with the highest values recorded during the cold season (heating season), which was due to fuel burning in local households, as well as low wind speed and air temperature, a low planetary boundary layer height and a stable atmosphere conducive to the accumulation of pollutants. In the Silesian Voivodship, the PM10-bound B[a]P concentrations of 4.0 to 27.7 ng were recorded from October to the end of March. Earlier studies conducted by Kozielska et al. [22] also indicated high concentrations of B[a]P in the 2009–2010 heating season in Silesia, and they are comparable to those presented in this work (Table 3). Significantly higher concentrations of B[a]P have only been found in India (Rajpur—43–76 ng·m−3), Pakistan (Lahore—81 ng·m−3) and China (Tangshan—61.6 ng·m−3) (Table 3). It should also be noted that, in the heating season, there may be as much as 380 μg B[a]P per 1 g of PM10, and, for example, in a coking plant, these values are in the range of 6.7 to 3356 μg·g−1 [82].
The possible emission sources around individual measuring stations can be indicated by analyzing the Spearman’s linear rank correlation coefficient between PM10 and B[a]P. The Spearman’s rank correlation coefficients and their corresponding levels of significance are presented in Table 4. Very strong positive correlations (R = 0.843 and R = 0.873) were found for Rybnik and Żywiec, thus suggesting that B[a]P comes from the same source of pollution, such as the combustion of solid fuels. Significant but slightly lower correlations (R = 0.723 and R = 0.721) between PM10 and B[a]P were observed in urban areas with high anthropogenic activity, especially in urban areas with elevated traffic levels (Katowice, Częstochowa). The obtained results suggest that, at the measuring stations exposed to the strong impact of the so-called low emissions from the combustion of solid fuels, the correlation between PM10 and B[a]P is higher than that at the sites exposed to both traffic emissions and fuel combustion in local furnaces.

3.4. Health Risk Assessment

The lifetime cancer risk (LCR) of inhalation exposure was estimated on the basis of the concentration of PM10-bound B[a]P in the atmospheric air in the Silesian Voivodeship. The lifetime lung cancer risk values calculated using IRIS IURB[a]P are shown in (Table 5). The WHO estimated an acceptable LCR of 1 × 10−6 (1 in 1,000,000) to 1 × 10−4 (1 in 10,000) for carcinogens. Under most guidelines used for LCR assessment, an LCR between 1 × 10−6 and 1 × 10−4 suggests a potential cancer risk, while a potential risk is high at LCR > 10−4 [15]. The LCR values were calculated for three scenarios: (a) annual, (b) the heating season and (c) the non-heating season. In the area of the Silesian Voivodeship, the LCR of the inhalation exposure to B[a]P was at an acceptable level. The lowest risk was found in the non-heating season, with values in the range of 0.32 × 10−6 (Bielsko-Biała) to 1.08 × 10−6 (Rybnik). For the heating season, these values were in the range of 2.8 × 10−6 to 12 × 10−6. The highest LCR occurred in Rybnik, where the average concentration of B[a]P in the heating seasons was at a level of 20 ng·m−3 (15.6–24.7 ng·m−3). Widziewicz et al. [95] analyzed B[a]P concentrations in Poland in the period between 2010 and 2015 using inhalation cancer unit risk factors recommended by EPA and the WHO. According to the authors, LCR was above the acceptable limits, and it was 7.33 × 10−4 in most conservative scenarios (Opolskie Voivodeship). If we take into account the same conservative scenario as Widziewicz et al. [95] (WHO), then during the heating season in Godów, Rybnik and Żywiec, the LCRs were also above the acceptable limits at 1.16 × 10−4, 1.73 × 10−4 and 1.25 × 10−4, respectively [95].
Applying the suggested IRIS IURB[a]P for lifetime (70 years) B[a]P exposure, the corresponding lifetime cancer risk is 3.88 × 10−6 on average while by WHO IURB[a]P is 56.33 × 10−6 for the measurements in The Silesian Voivodeship (2018–2021) (Table 6). Thus if 1,000,000 people are exposed to 6.5 ng·m−3 of ambient B[a]P for 70 years, then lung cancer may develop in 17–248 people on average, depending on the adopted IURB[a]P values (Table 6). If 1,000,000 people are exposed to 11.2 ng·m−3 (heating seasons) of ambient B[a]P, then lung cancer may develop in 30,429 people.
The presented results for the cancer risk of the inhalation exposure to B[a]P are understated and probably underestimate the carcinogenic potential of airborne PAH mixtures. For example, in the research conducted by Kozielska et al. [22,23], it was found that, in urban areas in the Silesian Voivodeship, the share of B[a]P in the total of 16 PAHs was about 9% (6–12%) on average, and for the B[a]P equivalent (TEQ) calculated by Nisbet and LaGoy [50], the value was between 56 and 68%. In [88], the share of B[a]P in TEQ was 58–62% in Indian semi-urban and urban sites. The estimation of cancer risk also has a significant impact on the toxicity equivalence factor scheme [96]. Moreover, there are additional factors contributing to carcinogenicity. One of the complicating factors is that PAHs in the air are bound to particles that may cause adverse health effects themselves. The carcinogenic potential of PAHs may even be enhanced when combined with those particles [91]. Additionally, in particles of different particle fractions (PM10, PM2.5 and PM1), the share of individual PAHs can vary significantly [97,98,99]; therefore, it can be expected that carcinogenic potency also differs substantially for those particle fractions [28,91].

4. Conclusions

Currently, the quality of atmospheric air is one of the most important environmental problems. This especially concerns developed and developing countries, including Poland. Unfortunately, for years, Poland has occupied one of the first places in the classification of countries with the most polluted air in Europe, and air quality standards are exceeded in its dominant area, especially when it comes to particulate matter and B[a]P. As many as 36 cities with the worst air quality in Europe are located in Poland, and most of them are located in Upper Silesia. Considering the fact that particulate matter is characterized by different particle sizes and chemical compositions, depending on the place of occurrence and the season, it is the most serious health hazard. In the Silesian Voivodeship in 2018–2021, the concentrations of PM and PM10-bound B[a]P very often exceeded the permissible concentrations regulated by national and European regulations, and they were higher than the concentrations recorded in other European cities. Considering the WHO guidelines for PM, the average daily concentrations of PM2.5 and PM10 were exceeded for almost the entire heating season. The average number of exceedances per year was 210 for PM2.5 and 75 for PM10, respectively. In Silesia, the fine-grained PM fraction clearly predominated over the coarse-grained PM fraction (PM2.5/PM10 ratio was in the range of 0.6533 to 0.8804). High concentrations of PM were strongly negatively correlated with wind speed and air temperature. The highest temporary PM concentrations reaching over 800 μg·m−3 were observed in the temperature range of −15 °C to −5 °C when the wind speed did not exceed 0.5 m·s−1. Moreover, the daily concentrations remained at a very high level at that time, exceeding the alarm level (150 µg·m−3 for PM10).
The concentrations of B[a]P recognized as a marker of the carcinogenic potential of a PAH mixture ranged from 4 ng·m−3 to around 28 ng·m−3 from October to the end of March each year, while EU regulations recommend a value of 1 ng·m−3, and the WHO recommends a value of 0.12 ng·m−3. Importantly, even 1 g of PM10 can contain 380 μg B[a]P.
According to IRIS, the estimated lifetime cancer risk (LCR) due to the inhalation exposure to PM10-related B[a]P concentrations in Silesia is acceptable. If we take into account the WHO scenario, then in Godów, Rybnik and Żywiec, the LCRs are above the acceptable limits during the heating season. The calculated lifetime lung cancer risk associated with exposure to B[a]P in the measured heating seasons in the Silesian Voivodeship suggested 30–429 cases per 1 million people depending on the IRIS, EPA or WHO scenario used.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph20010138/s1, Table S1. Detailed information on the location of selected urban background air quality monitoring stations in the Silesian Voivodeship. Table S2. Descriptive statistics of concentrations of PM2.5, PM10 and PM-bound B[a]P in the period of 2018–2021. Table S3. Number of days exceeding the values of the daily mean concentrations of PM10 in the heating season and during the year. Table S4. Number of days exceeding the values of the daily mean concentrations of PM2.5 in the heating season and during the year acc. to the WHO regulations.

Author Contributions

Conceptualization, D.K. and B.K.; methodology, D.K. and B.K.; formal analysis, D.K. and B.K.; investigation, D.K.; resources, B.K.; data curation, D.K.; writing—original draft preparation, D.K. and B.K.; writing—review and editing, B.K. and D.K.; visualization, D.K. and B.K.; supervision, D.K.; project administration, B.K.; funding acquisition, D.K. and B.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Faculty of Energy and Environmental Engineering, Silesian University of Technology (statutory research).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used in the paper are publicly available. We may supply the data that we gathered from public sources upon reasonable request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Ambient (Outdoor) Air Pollution. Available online: https://www.who.int/en/news-room/fact-sheets/detail/ambient-(outdoor)-air-quality-and-health (accessed on 25 August 2022).
  2. Mauderly, J.L.; Samet, J.M. Is there evidence for synergy among air pollutants in causing health effects? Environ. Health Perspect. 2009, 117, 1–6. [Google Scholar] [CrossRef] [PubMed]
  3. Billionnet, C.; Sherrill, D.; Annesi-Maesano, I. Estimating the health effects of exposure to multi-pollutant mixture. Ann. Epidemiol. 2012, 22, 126–141. [Google Scholar] [CrossRef] [PubMed]
  4. WHO. Compendium of WHO and other UN Guidance on Health and Environment, 2022 Update; WHO: Geneva, Switzerland, 2022. Available online: https://apps.who.int/iris/handle/10665/352844 (accessed on 25 August 2022).
  5. Badyda, A.J.; Grellier, J.; Dąbrowiecki, P. Ambient PM2.5 exposure and mortality due to lung cancer and cardiopulmonary diseases in Polish cities. Adv. Exp. Med. Biol. 2017, 944, 9–17. [Google Scholar] [CrossRef] [PubMed]
  6. European Commission. Air. Available online: https://environment.ec.europa.eu/topics/air_pl?etrans=pl (accessed on 22 July 2022).
  7. European Environment Agency. 2021 Country Fact Sheets. Available online: https://www.eea.europa.eu/themes/air/country-fact-sheets/2021-country-fact-sheets (accessed on 22 July 2022).
  8. European Environment Agency (EEA). Air Quality in Europe—2018 Report. 2018. Available online: http://www.eea.europa.eu/publications/air-quality-in-europe-2018 (accessed on 25 August 2022).
  9. European Environment Agency (EEA). Air Quality in Europe—2020 Report. 2020. Available online: https://www.eea.europa.eu/publications/air-quality-in-europe-2020-report (accessed on 25 August 2022).
  10. European Environment Agency (EEA). Air Quality in Europe—2021 Report. 2021. Available online: https://www.eea.europa.eu/publications/air-quality-in-europe-2021/ (accessed on 20 October 2022).
  11. Pehnec, G.; Jakovljević, I.; Godec, R.; Štrukil, Z.S.; Žero, S.; Huremović, J.; Džepina, K. Carcinogenic organic content of particulate matter at urban locations with different pollution sources. Sci. Total Environ. 2020, 734, 139414. [Google Scholar] [CrossRef] [PubMed]
  12. IARC International Agency for Research on Cancer: Monographs on the Evaluation of Carcinogenic Risks to Humans. 2015. Available online: https://monographs.iarc.who.int/list-of-classifications (accessed on 21 July 2022).
  13. IARC. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans; IARC International Agency for Research on Cancer: Lyon, France, 2005; Available online: http://www-cie.iarc.fr/ (accessed on 6 September 2022).
  14. Supplemental Information Toxicological Review of Benzo[a]pyrene. 2017. Available online: https://iris.epa.gov/static/pdfs/0136tr.pdf (accessed on 6 September 2022).
  15. EC. Council Directive 2004/107/EC Relating to Arsenic, Cadmium, Mercury, Nickel and Polycyclic Aromatic Hydrocarbons in Ambient Air. Available online: http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2005:023:0003:0016:EN:PDF (accessed on 6 September 2022).
  16. World Health Organization. Air Quality Guidelines for Europe, 2nd ed; WHO Regional Office for Europe: Copenhagen, Denmark, 2000; Available online: https://apps.who.int/iris/handle/10665/107335 (accessed on 12 September 2022).
  17. Annual Assessment of Air Quality in the Śląskie Voivodeship; Voivodship Report for 2021; CIEP: Katowice, Poland, 2022.
  18. Ustawa z Dnia 14 Grudnia 2012 r. o Odpadach (Dz.U. z 2022 r. poz. 699). Available online: https://isap.sejm.gov.pl/isap.nsf/DocDetails.xsp?id=WDU20220000699 (accessed on 22 September 2022).
  19. Muzyka, R.; Chrubasik, M.; Pogoda, M.; Tarnowska, J.; Sajdak, M. Py–GC–MS and PCA analysis approach for the detection of illegal waste combustion processes in central heating furnaces. Chromatographia 2019, 82, 1101–1109. [Google Scholar] [CrossRef] [Green Version]
  20. Cieślik, E.; Fabiańska, M.J. Preservation of geochemical markers during co-combustion of hard coal and various domestic waste materials. Sci. Total Environ. 2021, 768, 144638. [Google Scholar] [CrossRef]
  21. Kozielska, B.; Żeliński, J.; Cieślar, M. The occurrence of polycyclic aromatic hydrocarbons in bottom ash from individual heating devices. Zeszyty Naukowe SGSP 2022, 83, 7–18. [Google Scholar] [CrossRef]
  22. Kozielska, B.; Rogula-Kozłowska, W.; Klejnowski, K. Seasonal variations in health hazards from polycyclic aromatic hydrocarbons bound to submicrometer particles at three characteristic sites in the heavily polluted polish region. Atmosphere 2015, 6, 1–20. [Google Scholar] [CrossRef] [Green Version]
  23. Kozielska, B.; Rogula-Kozłowska, W.; Klejnowski, K. Selected organic compounds in fine particulate matter at the regional background, urban background and urban traffic points in Silesia (Poland). Int. J. Environ. Res. 2015, 9, 575–584. [Google Scholar]
  24. Křůmal, K.; Mikuška, P.; Večeřa, Z. Characterization of organic compounds in winter PM1 aerosols in a small industrial town. Atmos. Pollut. Res. 2017, 8, 930–939. [Google Scholar] [CrossRef]
  25. CSO. Central Statistical Office. Household Energy Consumption in 2018. 2020. Available online: https://stat.gov.pl/obszary-tematyczne/srodowisko-energia/energia/zuzycie-energii-w-gospodarstwach-domowych-w-2018-roku,2,4.html (accessed on 21 July 2022).
  26. Křůmal, K.; Mikuška, P. Mass concentrations and lung cancer risk assessment of PAHs bound to PM 1 aerosol in six industrial, urban and rural areas in the Czech Republic, Central Europe. Atmos. Pollut. Res. 2020, 11, 401–408. [Google Scholar] [CrossRef]
  27. Rogula-Kozłowska, W.; Kozielska, B.; Rogula-Kopiec, P. Road traffic effects in size-segregated ambient particle-bound PAHs. Int. J. Environ. Res. 2016, 10, 531–542. [Google Scholar]
  28. Rogula-Kozłowska, W. Traffic-generated changes in the chemical characteristics of size-segregated urban aerosols. Bull. Environ. Contam. Toxicol. 2014, 93, 493–502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Myllyvirta, L.; Howard, E. Five Things We Learned from the World’s Biggest Air Pollution Database. Available online: https://unearthed.greenpeace.org/2018/05/02/air-pollution-cities-worst-global-data-world-health-organisation/ (accessed on 21 July 2022).
  30. Widziewicz, K.; Rogula-Kozłowska, W.; Loska, K.; Kociszewska, K.; Majewski, G. Health risk impacts of exposure to airborne metals and benzo(a)pyrene during episodes of high PM10 concentrations in Poland. Biomed. Environ. Sci. 2018, 31, 23–36. [Google Scholar] [CrossRef] [PubMed]
  31. Rawicki, K.; Czarnecka, M.; Nidzgorska-Lencewicz, J. Regions of pollution with particulate matter in Poland. E3S Web Conf. 2018, 28, 01025. [Google Scholar] [CrossRef] [Green Version]
  32. Kowalska, M.; Zejda, J.E. Relationship between PM2.5 concentration in the ambient air and daily exacerbation of respiratory diseases in the population of Silesian Voivodeship during winter smog. Med. Pr. 2018, 69, 523–530. [Google Scholar] [CrossRef] [PubMed]
  33. Kowalska, M.; Skrzypek, M.; Kowalski, M.; Cyrys, J.; Niewiadomska, E.; Czech, E. The relationship between daily concentration of fine particulate matter in ambient air and exacerbation of respiratory diseases in Silesian Agglomeration, Poland. Int. J. Environ. Res. Public Health 2019, 16, 1131. [Google Scholar] [CrossRef] [Green Version]
  34. Regulation of the Minister of Entrepreneurship and Technology of 21 February 2019 Amending the Regulation on Requirements for Solid Fuel Boilers; (Dz.U. 2019, pos. 363). Available online: https://isap.sejm.gov.pl/isap.nsf/DocDetails.xsp?id=WDU20190000363 (accessed on 22 September 2022).
  35. Act of 5 July 2018 Amending the Act on the Fuel Quality Monitoring and Scrutinizing System and the Act on the National Revenue Administration; (Dz.U. 2018, pos. 1654). Available online: https://isap.sejm.gov.pl/isap.nsf/DocDetails.xsp?id=WDU20180001654 (accessed on 22 September 2022).
  36. Resolution No. VI/21/12/2020 of the Sejmik of the Śląskie Voivodeship of 22 June 2020 on the Adoption of the “Air Protection Program for the Śląskie Voivodeship”; Katowice City Council: Katowice, Poland, 2020.
  37. Belaïd, F. Exposure and risk to fuel poverty in France: Examining the extent of the fuel precariousness and its salient determinants. Energy Policy 2018, 114, 189–200. [Google Scholar] [CrossRef]
  38. Primc, K.; Slabe-Erker, R.; Majcen, B. Constructing energy poverty profiles for an effective energy policy. Energy Policy 2019, 128, 727–734. [Google Scholar] [CrossRef]
  39. Kozielska, B.; Mainka, A.; Żak, M.; Kaleta, D.; Mucha, W. Indoor air quality in residential buildings in Upper Silesia, Poland. Build. Environ. 2020, 177, 106914. [Google Scholar] [CrossRef]
  40. Oliveira, M.; Slezakova, K.; Delerue-Matos, C.; Pereira, M.C.; Morais, S. Children environmental exposure to particulate matter and polycyclic aromatic hydrocarbons and biomonitoring in school environments: A review on indoor and outdoor exposure levels, major sources and health impacts. Environ. Int. 2019, 124, 180–204. [Google Scholar] [CrossRef] [PubMed]
  41. Kozielska, B.; Kaleta, D. Assessment of indoor benzene and its alkyl derivatives concentrations in offices belonging to University of Technology (Poland). Atmosphere 2021, 12, 51. [Google Scholar] [CrossRef]
  42. CSO. Central Statistical Office. Population. Population Status and Structure as Well as Vital Statistics in the Territorial Profile (as of June 30, 2021). 2021. Available online: https://stat.gov.pl/obszary-tematyczne/ludnosc/ludnosc/ludnosc-stan-i-struktura-ludnosci-oraz-ruch-naturalny-w-przekroju-terytorialnym-stan-w-dniu-30-06-2021,6,30.html (accessed on 21 July 2022).
  43. Aarnio, P.; Martikainen, J.; Hussein, T.; Valkama, I.; Vehkamäki, H.; Sogacheva, L.; Härkönen, J.; Karppinen, A.; Koskentalo, T.; Kukkonen, J.; et al. Analysis and evaluation of selected PM10 pollution episodes in the Helsinki Metropolitan Area in 2002. Atmos. Environ. 2008, 42, 3992–4005. [Google Scholar] [CrossRef]
  44. CIEP. Chief Inspectorate of Environmental Protection. Available online: https://powietrze.gios.gov.pl/pjp/archives (accessed on 21 June 2022).
  45. PN-EN 12341:2014-07; Ambient Air—Standard Gravimetric Measurement Method to Determine the Concentration of Mass Fractions PM10 or PM2.5 Particulate Matter. CEN: Brussels, Belgium, 2014.
  46. PN-EN 15549:2011; Air Quality—Standard Method for the Measurement of the Concentration of Benzo[a]pyrene in Ambient Air. CEN: Brussels, Belgium, 2011.
  47. Delgado-Saborit, J.M.; Stark, C.; Harrison, R.M. Carcinogenic potential, levels and sources of polycyclic aromatic hydrocarbon mixtures in indoor and outdoor environments and their implications for air quality standards. Environ. Int. 2011, 37, 383–392. [Google Scholar] [CrossRef] [PubMed]
  48. Zheng, H.; Xing, X.; Hu, T.; Zhang, Y.; Zhang, J.; Zhu, G.; Li, Y.; Qi, S. Biomass burning contributed most to the human cancer risk exposed to the soil-bound PAHs from Chengdu Economic Region, western China. Ecotoxicol. Environ. Saf. 2018, 159, 63–70. [Google Scholar] [CrossRef] [PubMed]
  49. Wu, M.; Luo, J.; Huang, T.; Lian, L.; Chen, T.; Song, S.; Wang, Z.; Ma, S.; Xie, C.; Zhao, Y.; et al. Effects of African BaP emission from wildfire biomass burning on regional and global environment and human health. Environ. Int. 2022, 162, 107162. [Google Scholar] [CrossRef]
  50. Nisbet, I.C.T.; LaGoy, P.K. Toxic Equivalency Factors (TEFs) for polycyclic aromatic hydrocarbons (PAHs). Regul. Toxicol. Pharmacol. 1992, 16, 290–300. [Google Scholar] [CrossRef]
  51. USEPA EPA. Health Effects Assessement Summary Tables; U.S. Environmental Protection Agency: Washington, DC, USA, 1993.
  52. Ramírez, N.; Cuadras, A.; Rovira, E.; Marcé, R.M.; Borrull, F. Risk assessment related to atmospheric polycyclic aromatic hydrocarbons in gas and particle phasesnear industrial sites. Environ. Health Perspect. 2011, 119, 1110–1116. [Google Scholar] [CrossRef] [Green Version]
  53. Masiol, M.; Hofer, A.; Squizzato, S.; Piazza, R.; Rampazzo, G.; Pavoni, B. Carcinogenic and mutagenic risk associated to airborne particle-phase polycyclic aromatic hydrocarbons: A source apportionment. Atmos. Environ. 2012, 60, 375–382. [Google Scholar] [CrossRef]
  54. Bootdee, S.; Chantara, S.; Prapamontol, T. Determination of PM2.5 and polycyclicaromatic hydrocarbons from incense burning emission at shrine for health risk assessment. Atmos. Pollut. Res. 2016, 7, 680–689. [Google Scholar] [CrossRef] [Green Version]
  55. U.S. Environmental Protection Agency (USEPA). Guidelines for Carcinogen Risk Assessment; (EPA/630/P-03/001F); USEPA: Washington, DC, USA, 2005. Available online: http://www2.epa.gov/osa/guidelines-carcinogen-risk-assessment (accessed on 21 July 2022).
  56. RAIS. The Risk Assessment Information System. Available online: https://rais.ornl.gov (accessed on 6 September 2022).
  57. IRIS. Integrated Risk Information System. 2017. Available online: https://iris.epa.gov (accessed on 6 September 2022).
  58. EPA. Toxicological Review of Benzo[a]pyrene—Exectuvie Summary. Available online: https://iris.epa.gov/static/pdfs/0136_summary.pdf (accessed on 6 September 2022).
  59. U.S. EPA. Health Effects Assessment for Benzo[a]pyrene; EPA/540/1-86/022 (NTIS PB86134335); U.S. Environmental Protection Agency: Washington, DC, USA, 1984.
  60. CEPA. Air Toxics Hot Spots Program Risk Assessment Guidelines; Office of Environmental Health Hazard Assessment: Oakland, CA, USA, 2005.
  61. WHO. Air Quality Guidelines for Europe, 2nd ed.; European Series, No. 91; World Health Organisation Regional Publications: Geneva, Switzerland, 2000. [Google Scholar]
  62. Rogula-Kozłowska, W.; Kozielska, B.; Błaszczak, B.; Klejnowski, K. The mass distribution of particle-bound PAH among aerosol fractions: A case-study of an urban area in Poland. In Organic Pollutants Ten Years after the Stockholm Convention-Environmental and Analytical Update; Puzyn, T., Mostrag-Szlichtyng, A., Eds.; InTech: Rijeka, Croatia, 2012; pp. 163–190. [Google Scholar]
  63. Kozielska, B.; Rogula-Kozłowska, W. Polycyclic aromatic hydrocarbons in particulate matter in the cities of Upper Silesia. Arch. Waste Manage. Environ. Prot. 2014, 16, 75–84. [Google Scholar]
  64. Annual Assessment of Air Quality in the Śląskie Voivodeship; Voivodship Report for 2020; CIEP: Katowice, Poland, 2021.
  65. Samek, L.; Stegowski, Z.; Styszko, K.; Furman, L.; Fiedor, J. Seasonal contribution of assessed sources to submicron and fine particulate matter in a Central European urban area. Environ. Pollut. 2018, 241, 406–411. [Google Scholar] [CrossRef] [PubMed]
  66. Annual Assessment of Air Quality in the Śląskie Voivodeship; Voivodship Report for 2018; CIEP: Katowice, Poland, 2019.
  67. Huremović, J.; Žero, S.; Bubalo, E.; Dacić, M.; Čeliković, A.; Musić, I.; Bašić, M.; Huseinbašić, N.; Džepina, K.; Cepić, M.; et al. Analysis of PM10, Pb, Cd, and Ni atmospheric concentrations during domestic heating season in Sarajevo, Bosnia and Herzegovina, from 2010 to 2019. Air Qual. Atmos. Health 2020, 13, 965–976. [Google Scholar] [CrossRef]
  68. Rajšić, S.; Mijić, Z.; Tasić, M.; Radenković, M.; Joksić, J. Evaluation of the levels and sources of trace elements in urban particulate matter. Environ. Chem. Lett. 2008, 6, 95–100. [Google Scholar] [CrossRef]
  69. Klejnowski, K.; Pastuszka, J.S.; Rogula-Kozłowska, W.; Talik, E.; Krasa, A. Mass size distribution and chemical composition of the surface layer of summer and winter airborne particles in Zabrze, Poland. Bull. Environ. Contam. Toxicol. 2011, 88, 255–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Han, Y.M.; Cao, J.J.; Lee, S.C.; Ho, K.F.; An, Z.S. Different characteristics of char and soot in the atmosphere and their ratio as an indicator for source identification in Xi‘an, China. Atmos. Chem. Phys. 2010, 10, 595–607. [Google Scholar] [CrossRef] [Green Version]
  71. Spindler, G.; Bruggemann, E.; Gnauk, T.; Gruner, A.; Muller, K.; Herrmann, H. A four-year size-segregated characterization study of particles PM10, PM2.5 and PM1 depending on air mass origin at Melpitz. Atmos. Environ. 2010, 44, 164–173. [Google Scholar] [CrossRef]
  72. Czernecki, B.; Półrolniczak, M.; Kolendowicz, L.; Marosz, M.; Kendzierski, S.; Pilguj, N. Influence of the atmospheric conditions on PM10 concentrations in Poznań, Poland. J. Atmos. Chem. 2017, 74, 115–139. [Google Scholar] [CrossRef] [Green Version]
  73. Wang, J.; Xie, X.; Fang, C. Temporal and spatial distribution characteristics of atmospheric particulate matter (PM10 and PM2.5) in Changchun and analysis of its influencing factors. Atmosphere 2019, 10, 651. [Google Scholar] [CrossRef] [Green Version]
  74. Pohjola, M.A.; Kousa, A.; Kukkonen, J.; Härkönen, J.; Karppinen, A.; Aarnio, P.; Koskentalo, T. The spatial and temporal variation of measured urban PM10 and PM2.5 in the Helsinki Metropolitan Area. Water Air Soil Pollut. Focus 2002, 2, 189–201. [Google Scholar] [CrossRef]
  75. Li, X.; Ma, Y.; Wang, Y. Temporal and spatial analyses of particulate matter (PM10 and PM2.5) and its relationship with meteorological parameters over an urban city in northeast China. Atmos. Res. 2017, 198, 185–193. [Google Scholar] [CrossRef]
  76. Oke, T.R. Boundary Layer Climates; Methuen: London, UK, 1995. [Google Scholar]
  77. Fiore, A.M.; Naik, V.; Spracklen, D.V.; Steiner, A.; Unger, N.; Prather, M.; Bergmann, D.; Cameron-Smith, P.J.; Cionni, I.; Collins, W.J.; et al. Global air quality and climate. Chem. Soc. Rev. 2012, 41, 6663–6683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Wang, J.; Ogawa, S. Effects of meteorological conditions on PM2.5 concentrations in Nagasaki, Japan. Int. J. Environ. Res. Public Health. 2015, 12, 9089–9101. [Google Scholar] [CrossRef] [PubMed]
  79. Anquetin, S.; Guilbaud, C.; Chollet, J.-P. Thermal valley inversion impact on the dispersion of a passive pollutant in a complex mountainous area. Atmos. Environ. 1999, 33, 3953–3959. [Google Scholar] [CrossRef]
  80. Wielgosiński, G.; Czerwińska, J. Smog Episodes in Poland. Atmosphere 2020, 11, 277. [Google Scholar] [CrossRef] [Green Version]
  81. Regulation of the Minister for the Environment of 24 August 2012 on Levels of Certain Substances in the Air; (Dz. U. 2021, pos. 845). Available online: https://isap.sejm.gov.pl/isap.nsf/DocDetails.xsp?id=WDU20210000845 (accessed on 22 September 2022).
  82. Kozielska, B.; Konieczyński, J. Polycyclic aromatic hydrocarbons in particulate matter emitted from coke oven battery. Fuel 2015, 144, 327–334. [Google Scholar] [CrossRef]
  83. Mehmood, T.; Zhu, T.; Ahmad, I.; Li, X. Ambient PM2.5 and PM10 bound PAHs in Islamabad, Pakistan: Concentration, source and health risk assessment. Chemosphere 2020, 257, 127187. [Google Scholar] [CrossRef]
  84. Kalim, I.; Zahra, N.; Saeed, M.K.; Gilani, R.; Munawar, A.; Jalees, M.I. Polycyclic aromatic hydrocarbons in the atmosphere of data Darbar Chowk of Lahore, Pakistan. Bangladesh J. Sci. Ind. Res. 2020, 55, 147–152. [Google Scholar] [CrossRef]
  85. Kaur, S.; Kumar, B.; Chakraborty, P.; Kumar, V.; Kothiyal, N.C. Polycyclic aromatic hydrocarbons in PM10 of a north-western city, India: Distribution, sources, toxicity and health risk assessment. Int. J. Environ. Sci. Technol. 2022, 19, 1041–1056. [Google Scholar] [CrossRef]
  86. Ramteke, S.; Patel, K.S.; Sahu, B.L.; Deb, M.K.; Giri, B.; Aggarwal, S.G.; Ren, H.; Fu, P. Spatial variation, distribution and source impacts in urban organic aerosols. Asian J. Chem. 2018, 30, 2582–2590. [Google Scholar] [CrossRef]
  87. Ray, D.; Chatterjee, A.; Majumdar, D.; Ghosh, S.K.; Raha, S. Polycyclic aromatic hydrocarbons over a tropical urban and a high altitude Himalayan Station in India: Temporal variation and source apportionment. Atmos. Res. 2017, 197, 331–341. [Google Scholar] [CrossRef]
  88. Kumar, A.; Ambade, B.; Sankar, T.K.; Sethi, S.S.; Kurwadkar, S. Source identification and health risk assessment of atmospheric PM2.5-bound polycyclic aromatic hydrocarbons in Jamshedpur, India. Sustain. Cities Soc. 2020, 52, 101801. [Google Scholar] [CrossRef]
  89. Islam, N.; Dihingia, A.; Khare, P.; Saikia, B.K. Atmospheric particulate matters in an Indian urban area: Health implications from potentially hazardous elements, cytotoxicity, and genotoxicity studies. J. Hazard. Mater. 2020, 384, 21472. [Google Scholar] [CrossRef]
  90. Fang, B.; Zhang, L.; Zeng, H.; Liu, J.; Yang, Z.; Wang, H.; Wang, Q.; Wang, M. PM2.5-bound polycyclic aromatic hydrocarbons: Sources and health risk during non-heating and heating periods (Tangshan, China). Int. J. Environ. Res. Public Health 2020, 17, 483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Pehnec, G.; Jakovljevic, I. Carcinogenic potency of airborne polycyclic aromatic hydrocarbons in relation to the particle fraction size. Int. J. Environ. Res. Public Health 2018, 15, 2485. [Google Scholar] [CrossRef] [Green Version]
  92. Tsiodra, I.; Grivas, G.; Tavernaraki, K.; Bougiatioti, A.; Apostolaki, M.; Paraskevopoulou, D.; Gogou, A.; Parinos, C.; Oikonomou, K.; Tsagkaraki, M.; et al. Annual exposure to polycyclic aromatic hydrocarbons in urban environments linked to wintertime wood-burning episodes. Atmos. Chem. Phys. 2021, 21, 17865–17883. [Google Scholar] [CrossRef]
  93. Skiba, A.; Styszko, K.; Furman, P.; Dobrowolska, N.; Kistler, M.; Kasper-Giebl, A.; Zięba, D. Polycyclic aromatic hydrocarbons (PAHs) associated with PM10 collected in Wadowice, South Poland. E3S Web Conf. 2019, 108, 02007. [Google Scholar] [CrossRef]
  94. Furman, P.; Styszko, K.; Skiba, A.; Zięba, D.; Zimnoch, M.; Kistler, M.; Kasper-Giebl, A.; Gilardoni, S. Seasonal variability of PM10 chemical composition including 1,3,5-triphenylbenzene, marker of plastic combustion and toxicity in Wadowice, South Poland. Aerosol Air Qual. Res. 2021, 21, 200223. [Google Scholar] [CrossRef]
  95. Widziewicz, K.; Rogula-Kozłowska, W.; Majewski, G. Lung cancer risk associated with exposure to benzo(a)pyrene in Polish Agglomerations, cities, and other areas. Int. J. Environ. Res. 2017, 11, 685–693. [Google Scholar] [CrossRef] [Green Version]
  96. Ayoko, G.; Lim, M.; Morawska, L. Assessing Health Risk Associated with Airborne Polycyclic Aromatic Hydrocarbons by Chemometrics and Toxic Equivalency Factors. In Proceedings of the 17th International Clean Air and Environment Conference, Hobart, Tasmania, 3–6 May 2005. [Google Scholar]
  97. Teixeira, E.C.; Agudelo-Castañeda, D.M.; Fachel, J.M.G.; Leal, K.A.; Garcia, K.O.; Wiegand, F. Source identification and seasonal variation of polycyclic aromatic hydrocarbons associated with atmospheric fine and coarse particles in the Metropolitan Area of Porto Alegre. RS. Brazil. Atmos. Res. 2012, 118, 390–403. [Google Scholar] [CrossRef]
  98. Rogula-Kozłowska, W.; Kozielska, B.; Klejnowski, K.; Szopa, S. Hazardous compounds in urban PM in the central part of Upper Silesia (Poland) in winter. Arch. Environ. Prot. 2013, 39, 53–65. [Google Scholar] [CrossRef] [Green Version]
  99. Agudelo-Castañeda, D.M.; Teixeira, E.C.; Schneider, I.L.; Lara, S.R.; Silva, L.F.O. Exposure to polycyclic aromatic hydrocarbons in atmospheric PM1.0 of urban environments: Carcinogenic and mutagenic respiratory health risk by age groups. Environ. Pollut. 2017, 224, 158–170. [Google Scholar] [CrossRef] [PubMed]
Figure 1. “Smoggy” areas with higher PM concentrations in air in Poland [31]. Region I—smog incidents occur frequently and last the longest; Region II—smog incidents occur with lower frequency; Region III—smog incidents occur rarely.
Figure 1. “Smoggy” areas with higher PM concentrations in air in Poland [31]. Region I—smog incidents occur frequently and last the longest; Region II—smog incidents occur with lower frequency; Region III—smog incidents occur rarely.
Ijerph 20 00138 g001
Figure 2. Location of selected air monitoring stations in Silesian Voivodeship.
Figure 2. Location of selected air monitoring stations in Silesian Voivodeship.
Ijerph 20 00138 g002
Figure 3. Number of days exceeding the values of the WHO-recommended daily mean concentrations of PM2.5.
Figure 3. Number of days exceeding the values of the WHO-recommended daily mean concentrations of PM2.5.
Ijerph 20 00138 g003
Figure 4. Number of days exceeding the values of the recommended daily average concentrations of PM10 (a) according to national standards and (b) according to the WHO in Poland.
Figure 4. Number of days exceeding the values of the recommended daily average concentrations of PM10 (a) according to national standards and (b) according to the WHO in Poland.
Ijerph 20 00138 g004
Figure 5. Changes in PM10 concentration depending on the (a) air temperature and (b) wind speed that occurred in the 2018 heating season in Rybnik.
Figure 5. Changes in PM10 concentration depending on the (a) air temperature and (b) wind speed that occurred in the 2018 heating season in Rybnik.
Ijerph 20 00138 g005
Figure 6. Distribution of 1 h PM10 concentrations and recorded meteorological parameters from 1–6 March 2018 in Rybnik.
Figure 6. Distribution of 1 h PM10 concentrations and recorded meteorological parameters from 1–6 March 2018 in Rybnik.
Ijerph 20 00138 g006
Figure 7. Seasonal variability of PM-bound B[a]P in individual weeks in 2018–2021.
Figure 7. Seasonal variability of PM-bound B[a]P in individual weeks in 2018–2021.
Ijerph 20 00138 g007
Table 1. The PM2.5/PM10 ratio in the period of 2018–2021 by heating and non-heating seasons.
Table 1. The PM2.5/PM10 ratio in the period of 2018–2021 by heating and non-heating seasons.
GodówKatowiceTarnowskie GóryBielsko BiałaCzęstochowa
2018
year0.74530.72100.7235--
heating season0.79990.80240.7629--
non-heating season0.68550.64510.6800--
2019
year0.78400.70670.6533--
heating season0.78820.75550.6573--
non-heating season0.77990.66240.497--
2020
year0.81040.69000.6833--
heating season0.83270.72250.6938--
non-heating season0.77610.65230.6705--
2021
year0.80630.88040.71780.78250.8066
heating season0.85590.98770.76170.85770.8229
non-heating season0.76340. 75570.67550.71580.7916
Table 2. Spearman rank correlation coefficients between PM10 and selected meteorological parameters in 2018–2021.
Table 2. Spearman rank correlation coefficients between PM10 and selected meteorological parameters in 2018–2021.
TauaPa
PM10−0.386 *−0.614 *0.164 *
Ta—air temperature, ua—wind speed, Pa—air pressure, * significance p < 0.001.
Table 3. Comparisons of concentrations of PM and PM-bound B[a]P at various locations.
Table 3. Comparisons of concentrations of PM and PM-bound B[a]P at various locations.
LocationSampling PeriodSampling PointFraction
PM
PM,
µg·m−3
B[a]P, ng·m−3Reference
Islamabad, Pakistan2017clean urban PM2.551.591.55[83]
PM10102.792.3
Lahore, Pakistan trafficPM10188.781.4[84]
Amritsar, IndiaXI 2013–I 2014urban PM10196–4621.6[85]
Raipur, IndiaXI 2013–I 2014urban PM10468
(291–783)
43–76[86]
Kolkata, IndiaX 2015–V 2016urbanPM10nd5.3[87]
Darjeeling, India 2.5
Jamshedpur, IndiaXII 2016–V 2017semi-urbanwinterPM2.5nd1.24[88]
summer2.45
urbanwinter4.74
summer3.42
ruralwinter1.95
summer4.39
Jorhat, IndiaII–III 2018urban PM2.5126.8-[89]
PM10260.6-
X–XI 2018PM2.5159.04.77
PM10274.25.48
Tangshan, China2014urbanPM2.523–3673.64 *
61.6 **
[90]
Mlada Boleslav, Czech Republic15–28 II 2013residential districtPM1260.78[24]
Zagrzeb, Croatia2014urbanwinterPM1nd2.23[91]
spring0.21
summer0.03
autumn0.59
winterPM2.52.39
spring0.27
summer0.03
autumn0.73
winterPM107.66
spring0.75
summer0.07
autumn1.35
Sarajevo, Bosnia and HerzegovinaXII 2017–II 2018urban backgroundPM10507.28[11]
Zagreb, CroatiaXII 2017–II 2018residential403.24
Athenus, GreekXII 2016–I 2018urban backgroundwinterPM2.5nd0.81[92]
spring0.02
summer0.04
autumn0.10
Silesia, PolandVIII 2009–XII 2010regional backgroundnon-heatingPM2.519.832.43[23]
heating31.127.02
urban backgroundnon-heating23.802.89
heating75.3925.40
traffic non-heating27.975.65
heating50.5614.05
Wadowice, PolandIII 2017urban backgroundPM1034.4
(10.8–116.7)
11.1[93]
VIII 20171.4
Wadowice, PolandII–X 2017urban backgroundnon-heatingPM1027.11.1[94]
heating 43.34.98
Silesia, Poland2018urban backgroundPM1043.27.1This study
201935.16.3
202030.56.2
202132.36.3
* non-heating period; ** heating period; nd—no data.
Table 4. Spearman’s rank correlation coefficients between PM10 and B[a]P in 2018–2021.
Table 4. Spearman’s rank correlation coefficients between PM10 and B[a]P in 2018–2021.
R Spearmant (N−2)p
Częstochowa0.72114.8550.000
Katowice0.72315.0420.000
Rybnik0.84322.5310.000
Żywiec0.87325.5450.000
R Spearman—correlation coefficient, p—significance level, t—significance test result for correlation coefficient.
Table 5. Descriptive data on lifetime cancer risk (LCR) from the inhalation exposure to PM10-bound B[a]P in ambient air according to IRIS [57].
Table 5. Descriptive data on lifetime cancer risk (LCR) from the inhalation exposure to PM10-bound B[a]P in ambient air according to IRIS [57].
Air Monitoring StationAnnualHeating SeasonNon-Heating Season
Bielsko-Biała2.36 × 10−64.05 × 10−60.32 × 10−6
Częstochowa1.66 × 10−62.79 × 10−60.33 × 10−6
Godów5.19 × 10−68.66 × 10−60.82 × 10−6
Katowice2.38 × 10−64.15 × 10−60.42 × 10−6
Rybnik6.89 × 10−611.95 × 10−61.08 × 10−6
Żywiec4.53 × 10−68.00 × 10−60.62 × 10−6
Table 6. Lifetime cancer risk and number of additional lung cancer cases due to exposure to long-term average (2018–2021) concentrations of B[a]P (ng·m−3) in the Silesian Voivodeship according to IRIS, EPA and WHO.
Table 6. Lifetime cancer risk and number of additional lung cancer cases due to exposure to long-term average (2018–2021) concentrations of B[a]P (ng·m−3) in the Silesian Voivodeship according to IRIS, EPA and WHO.
PeriodIRIS [57]EPA [60]WHO [61]
LCRannual3.88 × 10−67.12 × 10−656.33 × 10−6
heating season6.72 × 10−612.32 × 10−697.44 × 10−6
non-heating season0.60 × 10−61.10 × 10−68.70 × 10−6
Number of additional lung cancer casesannual1731248
heating season3054429
non-heating season3538
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kaleta, D.; Kozielska, B. Spatial and Temporal Volatility of PM2.5, PM10 and PM10-Bound B[a]P Concentrations and Assessment of the Exposure of the Population of Silesia in 2018–2021. Int. J. Environ. Res. Public Health 2023, 20, 138. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20010138

AMA Style

Kaleta D, Kozielska B. Spatial and Temporal Volatility of PM2.5, PM10 and PM10-Bound B[a]P Concentrations and Assessment of the Exposure of the Population of Silesia in 2018–2021. International Journal of Environmental Research and Public Health. 2023; 20(1):138. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20010138

Chicago/Turabian Style

Kaleta, Dorota, and Barbara Kozielska. 2023. "Spatial and Temporal Volatility of PM2.5, PM10 and PM10-Bound B[a]P Concentrations and Assessment of the Exposure of the Population of Silesia in 2018–2021" International Journal of Environmental Research and Public Health 20, no. 1: 138. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph20010138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop