3.1. PM2.5 and PAHs concentrations
The concentrations of PM2.5 and 19 individual PAH compounds in the collected samples were analyzed. Mean values of PM2.5 concentrations and individual PAH compounds in the four seasons are shown in Table 1. PM2.5 concentrations ranged from 8 to 291 μg/m3 with a mean concentration of 60.21 ± 53.96 μg/m3 and the total concentration of 19 PAH compounds ranged from 0.32 to 61.38 ng/m3 with a mean concentration of 4.65 ± 8.54 ng/m3.
The mean PM2.5 mass concentrations in all seasons exceeded target value of 25 µg/m3 given by the WHO Directive (WHO 2005), whereas the national standard value is 35 µg/m3.
Seasonal variations of the PAH concentrations during the sampling period showed the highest and the lowest concentrations in fall and summer, respectively. One of the main reasons of increasing concentration in fall season could be due to temperature inversion which usually enhance the diffusion of PM2.5 and consequently affect PAH concentrations (Chen et al. 2020).
The dependence of PAHs concentration on atmospheric temperature and the increase in particulate PAHs concentration during the cold season have been reported by others (Chen et al. 2020; Vlachou et al. 2019; Tsapakis and Stephanou 2005; Li et al. 2006). Seasonal variations of emission sources, meteorological conditions, and the difference of gas-to-particle partitioning may result in a difference of PAHs concentration in cold and warm seasons (Tan et al. 2006). Reduced atmospheric dispersion as well as reduced photochemical reactions in low atmospheric temperatures can lead to higher pollutant concentrations during cold season (Chen et al. 2020; Vlachou et al. 2019). In contrast, the increase in ambient temperature during the warm season may result in an increased evaporation of the particle phase of PAHs to the gas phase. In addition, the PAHs degradation by photochemical or thermal reactions in the atmosphere especially in the warmer seasons is well known (Dörr et al. 1996). Furthermore, increasing the fossil fuel consumption is another reason for the higher concentrations of PAHs in the cold seasons (Wu et al. 2014).
Table 2 shows the mean concentration of PM2.5 and ƩPAHs in different cities worldwide. The mean concentration of PM2.5 in Isfahan were greater than those in Tehran, Iran (Ali-Taleshi et al. 2020; Kermani et al. 2017), Zaragoza, Spain (Callén et al. 2014), Kanazawa, Japan (Xing et al. 2020), São Paulo, Brazil (Bourotte et al. 2005), Venice, Italy (Masiol et al. 2012), Islamabad, Pakistan (Mehmood et al. 2020), Thessaloniki, Greece (Tolis et al. 2015) and Taiwan, China (Chen et al. 2016) but lower than those in Anshan, China (Wang et al. 2020b), Beijing, China (Feng et al. 2018) and Guangzhou, China (Liu et al. 2015). The concentration of ƩPAHs in Isfahan was higher than those in Bangi, Malaysia (Khan et al. 2015), Zaragoza, Spain (Callén et al. 2014), Kanazawa, Japan (Xing et al. 2020) and Taiwan, China (Chen et al. 2016) but lower than those in Tehran, Iran (Ali-Taleshi et al. 2020; Taghvaee et al. 2018; Kermani et al. 2017), Jinan, China (Zhang et al. 2019), Seoul, Korea (Kang et al. 2020), Anshan, China (Wang et al. 2020b), Beijing, China (Feng et al. 2018) Thessaloniki, Greece (Tolis et al. 2015), Islamabad, Pakistan (Mehmood et al. 2020), São Paulo, Brazil (Bourotte et al. 2005), Venice, Italy (Masiol et al. 2012) and Guangzhou, China (Liu et al. 2015).
3.2. Distribution of PAHs along seasons
According to the number of rings, PAHs are classified into five groups including 2-rings, 3-rings, 4-rings, 5-rings, and 6-rings. The analyzed PAHs in this study have been classified as follows: 2-rings (Nap); 3-rings (Acy, Ace, Flu, Dbt, Phe and Ant); 4-rings (Fl, Pyr, BaA and Chr); 5-rings (BbF, BkF, BeP, BaP, Per and DahA) and 6-rings (BghiP and Ind). They have been further divided into low molecular weight (LMW, 2- and 3-rings PAHs) and high molecular weight (HMW, 4-, 5- and 6-rings PAHs). The concentration of high molecular weight (HMW) PAHs were significantly higher than of the low molecular weight (LMW) PAHs in all seasons (Fig. 2).
A similar composition pattern over the four seasons were observed with the highest contribution from the 5–6-ring (79.5 %), followed by the 4-ring (18.5 %) and the 3-ring PAHs (2 %). The LMW PAHs can be formed in the pyrolysis of uncombusted fossil fuels, but the HMW PAHs mainly originate from high-temperature combustion processes such as vehicular exhaust (Dachs et al. 2002). Therefore, predomination of HMW PAHs in the study area indicated that PAHs mainly originate from combustion sources.
3.3. Source apportionment of PAHs
3.3.1. Diagnostic Ratios (DRs)
The ratios of LMW/HMW, An/(An+Phe), (anthracene to anthracene plus phenanthrene), Flu/(Flu+Pyr), (fluoranthene to fluoranthene plus pyrene), BaA/(BaA+Chr), (Benzo[a]anthracene to Benzo[a]anthracene plus Chrysene) and IP/(IP+Bghi), (Indeno[1,2,3-cd]pyrene to Indeno[1,2,3-cd]pyrene plus Benzo[ghi]perylene) were used for identification of PAHs sources. Fig. 3 shows cross-plots of PAHs ratios illustrating the different source types. Similar distribution patterns were observed along different seasons. For An/An+Ph, 0.10 is taken as a threshold to discriminate petrogenic from combustion sources. Values < 0.1 is usually considered as a petroleum source, while values > 0.1 indicates a combustion source (Yunker et al. 2002; Chen et al. 2013). The calculated ratio of An/(An+Ph) ranged from 0.09 to 0.34 with a mean value of 0.18. This means that the main source of PAHs in the city is from combustion process.
The Flu/Flu+Pyr ratio is another ratio that has been used to determine the emission source of PAHs. For this ratio, values less than 0.40 are distinctive of petroleum sources, values between 0.4 and 0.5 are distinctive of liquid fossil fuel combustion, and values above than 0.50 is distinctive of biomass or coal combustion (Yunker et al. 2002). These ratio have been also used to distinguish between gasoline and diesel emissions; where the values lower and higher than 0.5 refer to gasoline and diesel combustion, respectively. (Ravindra et al. 2006). In this study the Flu/(Flu+Pyr) ratio was from 0.32 to 0.58 with a mean of 0.43. This shows the importance of gasoline and diesel emissions to Isfahan’s atmosphere.
The ratio BaA/BaA+Chr is also declarative of the PAHs sources. For this ratio, values lower than 0.20 indicates a petroleum source, ranging from 0.2 to 0.35 suggest either petroleum or combustion source and values higher than 0.35 implies a combustion source (Yunker et al. 2002; Akyuz and Cabuk 2010). Values for this ratio reported in the literature are: 0.22-0.55 for gasoline, 0.38-0.64 for diesel (Simcik et al. 1997) and 0.5 for coal (Tang et al. 2005). The calculated ratio of BaA/(BaA+ Chr) was from 0.17 to 0.43 with a mean of 0.34. This information reflects the significance of gasoline and diesel emissions to Isfahan city. Likewise, PAHs sources may be identified by the ratio of IP/(IP+Bghi). For this ratio, values lower than 0.20 indicates a petrogenic source, ranging from 0.20 to 0.50 suggests liquid fossil fuel combustion, and values greater than 0.50 indicates grass, wood, or coal combustion (Yunker at al. 2002; Chen et al. 2005). For this ratio some researchers documented 0.21-0.22 values for gasoline emissions (Rogge et al. 1993), 0.35–0.70 for diesel emissions (Pio et al. 2001; Alves et al. 2001; Tang et al. 2005; Grimmer et al. 1983). In this work the IP/(IP+Bghi) ratio was from 0.22 to 0.41 with a mean of 0.31 which confirms the importance of gasoline and diesel emissions.
The ratio of LMW/HMW in all samples was less than 1 indicating overall combustion sources of PAHs. In conclusion of the PAHs diagnostic ratios analysis, combustion of fossil fuels was found to be the main source of PAHs in ambient PM2.5 in Isfahan.
Obviously, diagnostic ratios are a useful technique to identify the sources of PAHs. But, due to difficulty to distinguish between different emission sources, this technique should be used with a great caution. Furthermore, degradation (e.g. photolysis) and reaction with other species in the atmosphere can alter distribution pattern of PAHs during their transfer from emission sources to receptor sites (Ravindra et al. 2008; Pongpiachan 2014). Also should be aware that coals and fossil fuels from diverse origins and countries can generate different proportions of PAHs (Masclet et al. 1987). However, interpretation of PAHs sources must be done based on such numerical data and on the knowledge available about the historical background of the study area.
3.3.2. Source identification using PMF model
The datasets included concentrations and uncertainties of PAHs were loaded into the US EPA PMF 5.0 model. Six factors (sources) of PAHs including: diesel combustion, gasoline combustion, industrial activities, natural gas combustion, evaporative-uncombusted and other sources were revealed by employing the model. The source profiles and the mass contributions of PAHs are shown in Fig. 4. The predicted and measured PAHs concentration showed a significant correlation (r2 = 0.99, p < 0.01).
Factor 1 explained 10.52% of the sum of measured PAHs. This factor was dominated by high and low molecular weight PAHs. Compounds such as Flu, Chr, Pyr, BaA, BeP, BaP, Ind and BghiP were predominant in this factor. These types of PAHs are related with industrial activities (Lin et al. 2011). Steel and iron industry can produce PAHs through different processes such as sintering, casting and cooling and coke manufacturing (Yang et al. 2002). Ciaparra et al. (2009) demonstrated that low and moderate molecular weight PAHs (Flu, Phe, Ant, FluA and Pyr) explained by coke making and HMW PAHs (BaP, IcdP, DahP, DacP and BghiP) arising from the sintering process. As several types of industries, including the steel and iron industries, power plants, brick and cement factories are located around the city, therefore, we attributed this factor to “Industrial activities”.
Factor 2 contributed 12.95% of all the total concentration of PAHs. This factor was influenced mainly by Nap and moderately by Flu and Phe. LMW PAHs such as Nap, Flu and Phe were introduced as markers of uncombusted petroleum (Liu et al. 2015, 2009; Marr et al. 1999). The higher contribution of this factor in the warm seasons and the association with LMW PAHs leads us to the conclusion that this factor arises from “evaporative-uncombusted” sources, mainly from the gasoline emissions from storage tanks and pumps in gasoline stations.
Factor 3 was heavily dominated by HMW PAHs. PAHs compound such as BghiP, Ind, Pyr and BeP were reported as chemical tracers of gasoline combustion by some researchers (Khan et al. 2015; Guo et al. 2003; Schauer et al. 2002). Gasoline is widely used in Isfahan and gasoline-fueled cars and motorcycles have been found to be important sources of hydrocarbons in Isfahan. In recent decades, gasoline has been used primarily as a fuel in automobiles. Vehicles, especially those for private use have increased since past few decades in Isfahan city. We therefore defined this factor as “gasoline combustion” emissions. The factor contributed to 28 % of total PAHs on average.
Due to ambiguous characteristics, the forth factor was left as “other sources”. A mixture of minor sources such as wood, biomass and coal combustion, natural dust and etc. seems to explain this factor.
The dominant PAHs in Factor 5 were BkF, BaP, Ind and BghiP. These compounds have been shown to be tracers of “diesel combustion” (Fang et al. 2016; Harrison et al. 1996; Wang et al. 2009). BbF, BkF, BaP and Ind have been reported as chemical tracers of diesel combustion in numerous studies (Yang et al. 2013; Chen et al. 2011; Lin et al. 2011; Khalili et al. 1995; Ravindra et al. 2008; Ma et al. 2014; Harrison et al. 1996; Wang et al. 2014, 2015). Diesel fuel has been widely used in trucks and public transportation in Isfahan city. Moreover, the influence of power plants, steel and iron industries using diesel as fuel may be also reflected in this factor. This factor contributed 22.4 % of total PAHs.
Factor 6 was explained predominantly by Phe, Flu, BaA, Chr and Pyr. Jamhari et al. (2014), Khan et al. (2017) and Simcik et al. (1999) applied these markers to identify the emission sources of natural gas and coal combustion. As the natural gas is the main source of heating and cooking in Isfahan, this factor defined as “natural gas combustion”. This factor contributed to 16.75 % of total PAHs.
The contributions of each PMF factor to the level of total PAHs in different seasons are shown in Fig. 5. The sources fingerprints in different seasons are shown in Fig. S1. For the overall data, the contributions are dominated by gasoline combustion (28%) followed by diesel combustion (22%), natural gas combustion (17%), evaporative-uncombusted (13%), industrial activities (11%), and other sources or unidentified sources (9%). Since the natural gas are used in heating systems of houses and industries particularly in cold seasons (fall and winter), it showed the higher contribution in comparison to warm seasons (spring and summer). The natural gas is mainly used by industrial plants as fuel in spring and summer in the region. Since the main public transportation service in Isfahan city is provided by buses, they constitute together with trucks the major source of diesel combustion.
To the best of our knowledge, no source apportionment studies have so far been performed to identify sources and quantify their contributions to ambient PMs in Isfahan. PAHs source apportionment in respirable particles in Tehran, the capital of Iran, revealed five main sources including: diesel combustion (56.3 %), gasoline combustion (15.5 %), wood combustion and incineration (13.0 %), industry (9.2 %), and road soil particle (6.0 %) (Moeinaddini et al. 2014). In other study in Tehran, source apportionment of PM2.5-bound PAHs by PMF model revealed five main sources including: diesel exhaust (22.3 %), unburned petroleum (15.6 %), industrial (7.5 %), gasoline exhaust (30.9%) and coal/biomass and natural gas combustion (23.6%) (Ali-Taleshi et al. 2020).