Next Article in Journal
YoutubeTM Content Analysis as a Means of Information in Oral Medicine: A Systematic Review of the Literature
Next Article in Special Issue
Characterization of Atmospheric Fine Particles and Secondary Aerosol Estimated under the Different Photochemical Activities in Summertime Tianjin, China
Previous Article in Journal
Sedentary Behaviour and Its Relationship with Early Vascular Ageing in the General Spanish Population: A Cross-Sectional Study
Previous Article in Special Issue
Characteristics, Source Contributions, and Source-Specific Health Risks of PM2.5-Bound Polycyclic Aromatic Hydrocarbons for Senior Citizens during the Heating Season in Tianjin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of PM2.5 in an Industrial City of Northern China: Mass Concentrations, Chemical Composition, Source Apportionment, and Health Risk Assessment

1
State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China
2
College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China
3
National Institute of Metrology, Beijing 100029, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(9), 5443; https://doi.org/10.3390/ijerph19095443
Submission received: 26 February 2022 / Revised: 25 April 2022 / Accepted: 27 April 2022 / Published: 29 April 2022

Abstract

:
Urban and suburban PM2.5 samples were collected simultaneously during selected periods representing each season in 2019 in Zibo, China. Samples were analysed for water-soluble inorganic ions, carbon components, and elements. A chemical mass balance model and health risk assessment model were used to investigate the source contributions to PM2.5 and the human health risks posed by various pollution sources via the inhalation pathway. Almost 50% of the PM2.5 samples exceeded the secondary standard of China’s air quality concentration limit (75 µg/m3, 24 h). Water-soluble inorganic ions were the main component of PM2.5 in Zibo, accounting for 50 ± 8% and 56 ± 11% of PM2.5 at the urban and suburban sites, respectively. OC and OC/EC decreased significantly in the past few years due to enhanced energy restructuring. Pearson correlation analysis showed that traffic emissions were the main source of heavy metals. The Cr(VI) concentrations were 1.53 and 1.92 ng/m3 for urban and suburban sites, respectively, exceeding the national ambient air quality standards limit of 0.025 ng/m3. Secondary inorganic aerosols, traffic emissions, and secondary organic aerosols were the dominant contributors to PM2.5 in Zibo, with the total contributions from these three sources accounting for approximately 80% of PM2.5 and the remaining 20% attributed to traffic emissions. The non-carcinogenic risks from crustal dust for children were 2.23 and 1.15 in urban and suburban areas, respectively, exceeding the safe limit of 1.0 in both locations, as was the case for adults in urban areas. Meanwhile, the carcinogenic risks were all below the safe limit, with the non-carcinogenic and carcinogenic risks from traffic emissions being just below the limits. Strict control of precursor emissions, such as SO2, NOx, and VOCs, is a good way to reduce PM2.5 pollution resulting from secondary aerosols. Traffic control, limiting or preventing outdoor activities, and wearing masks during haze episodes may be also helpful in reducing PM2.5 pollution and its non-carcinogenic and carcinogenic health impacts in Zibo.

1. Introduction

Particulate matter with an aerodynamic diameter ≤ 2.5 μm, known as PM2.5, is suspended in the air and is composed of various chemicals and particle sizes. PM2.5 resides in the atmosphere for long periods and can penetrate deep into the lungs [1,2]. It has attracted significant attention due to its adverse effects on climate, visibility, and especially human health [3,4,5,6]. Long-term exposure to PM2.5 has been proven to increase the risk of lung cancer, dementia, cardiovascular disease, chronic respiratory diseases, and depressive-like responses [7,8,9,10,11,12].
Because PM2.5 poses human health risks, it is necessary to explore its toxicity. Carbonaceous species, water-soluble ionic species (WSI), and mineral dust are the major components of PM2.5 [13,14]. Several studies suggest that to reduce the health risks associated with PM2.5, chemical components and their sources should be considered, alongside mass concentrations [15,16,17]. Given that a substantial body of research highlights the toxicological importance of metals in PM2.5, an investigation into metal levels is necessary to understand their risks, and identification of their emission sources is needed to develop effective risk management plans [1].
Zibo is a highly industrialized cluster city of about 6000 km2 with a population of 4.7 million in the central area of Shandong Province. Its many heavy and light industries, such as ceramic and petrochemical industries, make it a highly polluted area [18,19,20,21,22]. Enrichment factor analysis, principal component analysis, and the PMF model were used to analyse the potential sources of PM2.5 and related metal elements in Zibo in recent years [19,23,24,25]. However, receptor models have not often been used in previous studies, and sampling stations have been concentrated in urban areas. Therefore, the main objectives of this work were to investigate (a) the characteristics of PM2.5 and its chemical composition (water-soluble ionic species, carbonaceous species, and element species) in urban and suburban areas, (b) the main source apportionment of PM2.5 using the CMB model, and (c) the source-specific health risks for adults and children induced by PM2.5-bound elements via inhalation.

2. Materials and Methods

2.1. Site Description and Sample Collection

The surrounding environments of the sampling sites, namely Nanding (urban area) and Lishan (suburban area), are shown in Figure 1. The urban site was located approximately 6 m above the ground on the second-storey rooftop of a local monitoring station (36°48′10.67″ N, 118°01′28.19″ E) and was surrounded by residential and administrative buildings. The suburban sampling site was approximately 12 m above the ground on the fourth-storey rooftop of a building (36°11′13.76″ N, 118°10′26.08″ E) and was surrounded by residential and commercial buildings.
A total of 45 PM2.5 samples were simultaneously collected at the two sites each day during the periods of 8–22 January (winter), 16–25 April (spring), 30 July–8 August (summer), and 16–25 October (autumn) in 2019. Each day’s sampling time was defined as running from 10:00 a.m. until 9:00 a.m. the next day (23 h). The samplers (ZR-3930; Zhongrui Inc., Qingdao, China) had two parallel channels and a flow rate of 16.7 L min−1 and used Teflon filters (diameter = 47 mm; Whatman Inc., Buckinghamshire, UK) and quartz filters (diameter = 47 mm; Pall, San Diego, CA, USA). The PM2.5 manual sampling method was in accordance with the latest national standard, HJ 656-2013 [26].

2.2. Gravimetric and Sample Analysis

Before sampling, quartz filters were baked in a muffle furnace at 550 °C for 3 h to remove carbonaceous pollutants. Before and after sampling, the filters were balanced under constant temperature (20 ± 1 °C) and humidity (50 ± 5%) for more than 24 h and weighed using an automatic filter-weighing system (AWS-1; Comde-Derenda GmbH, Germany, readability = 1 μg). The difference between two consecutive measurements was not more than 40 µg. After sampling, the filters were stored at −20 °C until analysis to prevent volatilization of volatile components [27]. The concentrations of PM2.5 were determined based on the values analysed from the Teflon filters.
Quartz filters were used to analyse the carbon contents and water-soluble inorganic ions (WSIIs). Organic carbon (OC) and elemental carbon (EC) were measured according to interagency monitoring of protected visual environments (IMPROVE) thermal optical reflectance (TOR; Model 2001A, Desert Research Institute, Reno, NV, USA). Because the oxidation temperatures of OC and EC are different, the IMPROVE temperature-programmed method was used to analyse the samples in this study. The heating program was divided into two stages: in the first stage, in the oxygen-free and pure-helium environment, the quartz filter was heated to 140 °C, 280 °C, 480 °C, and 580 °C to obtain the four components of OC (OC1, OC2, OC3, and OC4, respectively); in the second stage, in a helium atmosphere containing 2% oxygen, the temperature was gradually increased to 580 °C, 740 °C, and 840 °C to obtain the three components of EC (EC1, EC2, and EC3, respectively). Because OC forms optical pyrolyzed carbon (OPC) during the carbonization process, according to the IMPROVE analysis protocol, OC and EC were defined according to Equations (1) and (2), respectively [28].
OC = OC1 + OC2 + OC3 + OC4 + OPC
EC = EC1 + EC2 + EC3 − OPC
Water-soluble inorganic ions (Na+, NH4+, K+, Mg2+, Ca2+, Cl, NO3, and SO42−) were measured using ion chromatography analysers (DIONEX ICS-1100 and ICS-2100, Thermo Company, Waltham, MA, USA). The samples were soaked in deionized water (10 mL), shaken well, leached under ultrasonic vibration for 15 min, and let stand for 5 min; then, the supernatant was injected for ion chromatographic analysis [21].
To analyse the concentrations of elements, one 1/2 of a Teflon filter was cut into fragments and solubilized in a 5 mL HF acid solution in a PTFE crucible. The solution was then progressively heated to 160 °C on an electric stove for 2 h and leached in 10 mL of 2% HCl. After cooling, the solution was moved into a plastic colorimetric tube before analysis of eleven trace elements (As, Cd, Co, Cr, Cu, Mn, Ni, Pb, Sr, V, and Zn) using inductively coupled plasma-mass spectrometry (ICP-MS, Agilent 7500a, Santa Clara, CA, USA). The other 1/2 of Teflon filter was put in a Nickel crucible and progressively heated to 300 °C to ash in a muffle furnace and kept in the furnace under constant temperature for about 20 min. After cooling, pure ethanol and 0.2 g NaOH were added to crucible and melted for 10 min. Following the extraction from boiling water, the solution was moved into a plastic tube with 1 mL HCl (12 mol/L). After dilution, the other trace elements (Al, Mg, Ca, Fe, and Ba) were analysed by an inductively coupled plasma-optical spectrometer (ICP-OES, Agilent, Santa Clara, CA, USA) [29].

2.3. Source Apportionment Model

Source apportionment was conducted with a CMB model. Considering that a large number of samples (usually >100) is required for statistical analysis and there some subjectivity and uncertainty are involved in the determination of the number of source types and the discrimination in PMF model [30,31], we chose a CMB model rather than a PMF model. The theory behind the CMB model has is introduced in detail in [31]. The fundamental relationship between the concentration at a receptor site and source information in the CMB model can be expressed as Equation (3):
C i = j = 1 p F i j × S j  
where Ci is the concentration (µg m−3) of elemental component i at a receptor site, p represents the number of sources, Fji is the fraction of the ith element from the jth source, and Sj is the concentration of the contribution from the jth source to the receptor site.
In our study, 90 ambient PM2.5 samples and source profiles were input to CRAES CMB 1.0 to obtain the source apportionment results and contributions of various PM2.5-bound metals. CRAES CMB 1.0 is a Chinese-language version of the model derived from EPA CMB 8.2 by the Chinese Research Academy of Environmental Sciences, which includes an enumeration method [32]. Based on the emission inventory and industrial structure of Zibo, we chose the potential sources and conducted source sampling and analysis work. The squared regression coefficient (R2), sum of residual squares value (Χ2), degrees of freedom (DF), and percentage of explained mass to sample total mass were used to evaluate the validity of source contribution measures [33]. Table 1 summarizes the performance measures of the CMB model used in the present study.

2.4. Health Risk Assessment Model

The investigated model was introduced by the United States Environmental Protection Agency to quantify health risks (carcinogenic and non-carcinogenic) and human exposure to elements in dust [34]. The three main exposure pathways in humans are direct ingestion, inhalation, and dermal contact. PM2.5-bound metals cause human health hazards mainly via the respiratory system [31], so we only considered the inhalation exposure pathway in this study. Exposure is expressed in terms of a daily dose (mg kg−1 day−1), and the doses received through inhalation were calculated with Equations (4) and (5) [35,36]. Ba, Co, Cr(VI), and Mn exhibited non-carcinogenic risks, whereas As, Cd, Co, Cr(VI), and Ni exhibited carcinogenic risks in this study [37,38]. Because the toxicity of Cr is attributed to its hexavalent state, its concentrations were scaled by a factor of 1/7, given that atmospheric Cr(VI) is reported to have a ratio of 1/6 relative to Cr(III) [39,40].
ADD inh = C × InhR × EF × ED PEF × BW × AT
LADD = C × EF AT × PEF × ( InhR child × ED child BW child + InhR adult × ED adult BW adult )
where ADDinh (mg kg−1 day−1) is the dose obtained via inhalation, and LADD (mg kg−1 day−1) is the lifetime average daily dose of carcinogens. The exposure parameters are listed in Table 2.
To assess the health risks due to trace element exposure in PM2.5, the hazard quotient (HQ), hazard index (HI), and carcinogenic risk (RI) techniques were employed and calculated with Equations (6)–(8) [35,36].
HQ = ADD × RfD
HI =   HQ i
RI = LADD × SF
where RfD (mg kg−1 day−1) is the corresponding reference dose, and SF is the corresponding slope factor. The values of RfD and SF used in this study are listed in Table 3.

3. Results and Discussion

3.1. PM2.5 Mass Concentrations

Figure 2 shows the PM2.5 mass concentrations in winter, spring, summer, and autumn during the sampling period at the urban and suburban sites in Zibo. The minimum, median, and maximum concentrations of urban PM2.5 were 20, 56, and 313 µg/m3, respectively, with 48.89% of the values exceeding the secondary standard of China’s air quality concentration limit (75 µg/m3, 24 h) [44]. The suburban minimum, median, and maximum concentrations of PM2.5 were 20, 65, and 212 µg/m3, respectively, with 42.22% of data representing concentrations > 75 µg/m3. Furthermore, 80% of the values at both sites exceeded the primary standard limit (35 µg/m3, 24 h), and 100% exceeded the WHO air quality guidelines for PM2.5 (15 µg/m3, 24 h) [44,45]. By comparing the statistics of the two sites, we can see that the median and maximum concentrations at the urban site were greater than those at the suburban site, especially the maximum values. Compared with sites in the North China Plain (NCP; northwest of Zibo), the proportion of polluted weather (PM2.5 > 75 µg/m3) at the urban (67%) and suburban (67%) sites during the winter sampling period were greater than that at the NCP sites (50%). However, the overall concentrations in summer were <75 µg/m3 at both the urban and suburban sites, whereas the NCP sites experienced haze pollution during their intensive summer observation period [46]. In other words, Zibo experienced greater winter pollution than the NCP but less summer pollution.

3.2. Chemical Composition Levels

3.2.1. WSII Levels

The annual mean WSII concentrations of PM2.5 at the urban and suburban sites (Table 4) were 42.2 ± 35.0 and 44.0 ± 31.3 µg/m3, accounting for 50 ± 8% and 56 ± 11% of the total PM2.5, respectively, which is close to the results for 2018 [27]. This result shows that WSIIs were the main components of PM2.5 in both urban and suburban areas. The secondary water-soluble inorganic ions composed of sulphate (SO42−), nitrate (NO3), and ammonium (NH4+) aerosol (SNAs) were the major contributors to WSIIs, and their annual average concentrations were 37.8 ± 32.7 µg/m3 (approximately 44% of the PM2.5 in urban site) and 39.7 ± 29.8 µg/m3 (approximately 50% of the PM2.5 in suburban site). SNAs are formed in the atmosphere from gaseous precursors, for example, sulphur dioxide (SO2), nitrogen oxides (NOx), and ammonia (NH3), through complex chemical reactions [47]. The ratio of SNA/PM2.5 was highest in autumn (56% for both urban and suburban sites), which is indicative of intense chemical reactions in this season.
The ratio of NO3 to SO42− ([NO3/SO42−]) has been used as an indicator of the relative importance of mobile vs. stationary sources of nitrogen and sulphur in the atmosphere [48,49,50]. Ratios > 1.0 indicate that the particle sources at the sampling site are likely dominated by mobile sources, whereas stationary sources predominate at ratios < 1.0. In this study, the ratios of NO3/SO42− were 2.5, 2.2, 0.5, and 3.9 (2.5, 2.7, 0.8 and 4.2) in winter, spring, summer, and autumn, respectively, with an annual average ratio of 2.3 (2.5) at the urban site (suburban site, Table 4). Note that, in summer, a lower NO3/SO42− ratio may be ascribed to high temperatures, which lead to the evaporation of NH4NO3, whereas high ratios in other seasons are more likely related to traffic emissions from roads near the sampling sites [51].
There was a typical haze episode from 10 to 14 January 2019 in Zibo according to the definition of haze episode [27]; Figure 3 gives the daily concentrations of PM2.5, SNA, and NO3/SO42− during this episode. The average concentrations of PM2.5 were 197 µg/m3 and 180 µg/m3 at the urban and suburban sites, respectively, and the peak values were 312 and 212 µg/m3, respectively. Obviously, SNAs were the main components in PM2.5 during this episode, and the peak values of SNA/PM2.5 were over 60%. SNA components ranked in the order of NO3 > SO42− > NH4+, which was different from Beijing, where the SO42− ranked first [50]. In addition, NO3/SO42− was above 1.0 and increased during the accumulation period. The above information demonstrates that NO3 played a predominant role in the formation of heavy haze, and mobile sources had a relatively important contribution.

3.2.2. OC and EC Levels

OC and EC are also major components of PM2.5, and their daily concentrations are shown in Figure 4. The annual concentrations of OC in PM2.5 in urban and suburban areas were 10.9 and 11.2 µg/m3, respectively, which are lower than those in other Chinese cities, e.g., 13.8 µg/m3 in Changzhou [51] and 14.4 µg/m3 in Wuhan [52], and are comparable to those measured in Guangzhou (9.4 µg/m3) [5]. The annual EC concentrations were 4.7 µg/m3 in urban and 5.1 µg/m3 in suburban areas, which are close to those of Changzhou (5.4 µg/m3) [51], Wuhan (5.2 µg/m3), and Guangzhou (6.2 µg/m3) [5]. Obvious similar seasonal variations in OC and EC were observed in both urban and suburban areas, with the highest average concentrations occurring in winter: 18.6 and 6.8 µg/m3 in urban and 19.1 and 8.2 µg/m3 in suburban areas, respectively. Moreover, the ratios of OC/PM2.5 were highest in winter, accounting for 19% and 20% of the PM2.5 in urban and suburban areas, respectively. The lowest concentrations of OC and EC were observed in summer.
The OC/EC ratio was used to indicate the existence of secondary organic carbon (SOC). More than 60% of OC/EC ratios in this study exceeded 2.0, which reflects the fact that SOC is a significant component of PM2.5 [53]. However, the annual average OC/EC ratio in this study was only 2.4, which is much lower than that reported in a previous study (OC/EC = 6.68, OC = 31.97 ± 20.43 μg/m3, EC = 5.01 ± 2.68 μg/m3; sampling time = 2006–2007) [25]. Clearly, the distinction between annual average OC concentrations caused this result. This might be attributed to the partial replacement of coal-fired energy with clean energy over the past decade, especially after the Chinese State Council implemented the Air Pollution Prevention and Control Action Plan in September 2013. According to a study by Li et al., study [54], OC emissions decreased by 31.7% in Beijing and surrounding areas (including Zibo) during the period of 2014–2017 as a result of strict governance of residential fuel and enhanced energy restructuring by using cleaner natural gas and electricity replacing coal.

3.2.3. Element Levels

The seasonal average concentrations of various elements in PM2.5 are shown in Figure 5. The total elemental concentrations accounted for 1.96% and 1.33% of the mean PM2.5 concentrations at the urban and suburban site, respectively. In urban PM2.5 samples, Ca exhibited the highest levels (453.9), followed by Al (384.7), Zn (234.6), Fe (192.0), Mg (127.5), Mn (34.5), Pb (34.0), Ba (12.4), Cu (11.4), Cr (10.7), As (6.0), Ni (4.1), Sr (3.6), V (1.7), Cd (1.6), and Co (0.3; median values, ng/m3). Similarly, Ca exhibited the highest levels in the suburban PM2.5 samples (261.7), followed by Al (152.0), Fe (116.4), Zn (106.7), Mg (100.1), Mn (24.8), Pb (18.4), Cr (13.5), Cu (5.3), As (3.8), Ni (3.0), Ba (2.6), V (1.2), Sr (1.1), Cd (0.7), and Co (0.2; medians, ng/m3). The annual median concentrations of urban and suburban heavy metals Pb, Cd, and As were lower than the national ambient air quality standards limits of 500 ng/m3, 5 ng/m3, and 6 ng/m3, respectively. However, the Cr(VI) concentrations exceeded the limit of 0.025 ng/m3 [45]. Shandong province, where the Zibo is located, has the highest Cr emissions among all provinces in China due to its high capacity for steel production and coal consumption of industrial boilers [55].
Figure 6 shows the annual trends of 16 metal elements in urban and suburban areas, and their Pearson correlations are shown in Figure 7. The urban and suburban sites have several common characteristics. Al, Ca, and Fe were significantly correlated with one another (p < 0.01, r > 0.48), and their concentrations increased significantly in spring, which is the season of frequent sand/dust storm weather in northern China [14], increasing these crustal elements. Additional, Mn, Co, Cu, Zn, Cd, and Pb were also significantly correlated with one another (p < 0.01, r > 0.86), all of which are typically elements of motor vehicle emissions sources or classic traffic-source-related elements [56,57]; most showed a significant increase in winter and autumn. The probable reasons are that vehicle emissions have strong time sensitivity, being highest in winter and lowest in summer [58], with clean travel modes, such as electric bicycles, being more commonly used in summer. Arsenic, a representative coal-combustion element, showed higher concentrations in winter than in other seasons due to increased demands for space heating [13].
However, there were some distinctions between the urban and suburban sites; e.g., V and Ni were not correlated with any of the other 15 elements at the suburban site, whereas the concentration of Ni was significantly lower in winter. This is probably due to the fact that agricultural and construction activities are more common in suburban areas than in urban areas and are less prevalent in winter, resulting in a significant decrease in the use of related off-road machinery and corresponding exhaust emissions [59,60]. Sr and Ba, which are typical alkaline-earth metals emitted by industrial processes such as alloy smelting, ceramic production, etc., were significantly related to each other in urban samples (p < 0.01, r = 0.81).

3.3. Source Apportionment Using CMB

3.3.1. CMB Source Profiles

Chemical mass balance models for source apportionment require information about the chemical characteristics of the sources that are likely to affect pollutant concentrations at receptor sites [33]. The weight percentages of the PM2.5-bound chemical species in all sources (crustal dust, coal boilers, graphite production, aluminium smelting, glass manufacturing, iron smelting, traffic emissions, secondary inorganic aerosols, and secondary organic aerosols) are presented in Figure 8.
It is evident from the source profiles that the major components of crustal dust are Ca, OC, SO42−, Si, and Al (Ca, Si, and Al are crustal elements). The high percentage of OC is an indication of the mixing of dust and organic matter during aging or by entrainment of organic materials from soils [61]. The coal boiler profile is dominated by SO42−, OC, NH4+, Na, Cl, Si, and Al, with abundances of 23.9%, 18.0%, 9.9%, 3.4%, 2.5%, 2.4%, and 2.0%, respectively. This is consistent with the results of a study by Han et al. [62] in which OC, EC, SO42−, Ca, Si, Fe, and Al were found to be the abundant species in coal combustion fly ash from Xining City, China. The source profiles for graphite production and glass manufacturing are similar; therefore, SO42−, NH4+, and OC were the major components, and Zn [63] and Na [64] were also enriched in graphite production and glass manufacturing, respectively. Aluminium smelting and iron smelting are characterized by large amounts of Al (19.6%) and Fe (9.6%) [1], respectively. Carbon components are the most abundant species in traffic emissions. It is worth noting that the abundance of EC in traffic emissions is 14.3%, which is 5–27 times that in other profiles. EC was significantly correlated with traffic-source-related elements (Mn, Co, Cu, Zn, Cd, and Pb) (p < 0.01) in ambient PM2.5 samples. Hence, EC can be used as a key marker to distinguish traffic emissions from other sources, as it primarily arises from engines [13].

3.3.2. CMB Source Apportionment

The CMB results identified nine sources of PM2.5 during the winter, spring, summer, and autumn sampling periods, including crustal dust (CD), coal boilers (CB), graphite production (GP), aluminium smelting (AS), glass manufacturing (GM), iron smelting (IS), traffic emissions (TE), secondary inorganic aerosols (SIA, including NO3, SO42−, and NH4+), and secondary organic aerosols (SOA; Table 5). The top annual contributors were ranked as SIA > TE > SOA, with percentage contributions of 43.34%, 24.03%, and 9.78% in urban and 49.44%, 20.49%, and 10.57% in suburban areas, respectively. The percentage contributions of SIA were nearly 50%, which is consistent with the high concentrations of SNAs in PM2.5. The high contribution of traffic emissions might be attributable to the rapid increase in motor vehicle numbers over the past 10 years (from 443,179 to 1,136,377; Zibo Statistical Yearbook, 2010, 2019; http://tj.zibo.gov.cn/col/col886/index.html; accessed on 31 October 2021). Apart from TE and SIA, SOAs also played an important role as a PM2.5 source in Zibo. Zibo is an important petrochemical industry base of China, and there were amounts of volatile organic compounds (VOCs) emissions (100 kt/year) [18] that could result in the formation of SOAs [25,65]. In addition, the percentage contribution of crustal dust was lower than that in in Tianjin (26.4%), Langfang (21.4%), Baoding (16.9%) [66], and Zibo in 2006 (29.2%) [25]. Combined with the low percentage of elements, the main reason should be that the Zibo government implemented many strict dust management policies in the last few years [67]. Although the percentage contributions of coal boilers, graphite production, aluminium smelting, glass manufacturing, and iron smelting were all lower than 3%, their contributions to gaseous precursors, such as SO2 and NOx, could not be ignored because of the high contribution of SIA. The contribution rates were similar between the urban and suburban sites.

3.4. Risk Characterization

Because using the total concentration of heavy metals not only overestimated the human health risk but also inaccurately identified the high-risk pollution sources [68], we chose a source-specific human health risk assessment model based on the inhalation mode of exposure. The results are shown in Table 6.
Crustal dust had the highest non-carcinogenic risk (1.26 for adults and 2.23 for children at the urban site vs. 0.65 for adults and 1.15 for children at the suburban site), and except for 0.65, the other three values exceeded the safe limit 1.00. This was caused by the high HQ values of Mn (Table 6). Meanwhile, Mn had the highest HQ values of all sources, which is a similar result to that reported in a previous study [34]. Furthermore, the high HI values for children from crustal dust indicate that children are more vulnerable to non-carcinogenic health risks [34]. For the other four sources, both the HQ and accumulative HI values of individual emission sources identified in this study were lower than the safe limit of 1.00 for both adults and children, which indicates that there were no adverse non-carcinogenic health risks. However, the HI values of traffic emissions were 0.78 and 0.86 for children in urban and suburban areas, respectively, which are only a little lower than 1. For cancer risk (RI), the model estimated a non-significant cancer risk for both adults and children, with RI values lower than 1.00 × 10−6, suggesting that the carcinogenic risk from these elements is negligible. In addition, the RI values were highest for traffic emissions at both the urban and LS sit (1.60 × 107 and 1.33 × 10−7, respectively).

4. Conclusions

In this study, we investigated the characteristics of PM2.5 chemical composition, source apportionment, and health risk in Zibo. The results show that (1) almost 50% of PM2.5 samples exceeded the secondary standard of China’s air quality concentration limit (75 µg/m3, 24 h). Water-soluble inorganic ions were the main component of PM2.5 in Zibo, accounting for 50 ± 8% and 56 ± 11% at the urban and suburban site, respectively. Except in summer, the ratios of NO3 to SO42− exceeded 1, implying important contributions from mobile sources. The annual average concentrations of OC accounted for 15% of PM2.5, and enhanced energy restructuring resulted in a decrease in OC and OC/EC in the past few years. Ni exhibits significant seasonal variation due to agricultural and construction activities in suburban areas. The total elemental concentrations accounted for 1.96% and 1.33% of the mean PM2.5 concentrations at the urban and suburban sites, respectively, and the Cr(VI) concentrations were 1.53 and 1.92 ng/m3, respectively, which exceed the national ambient air quality Standards limit of 0.025 ng/m3 due to its large Zibo’s high capacity for coal consumption of industrial boilers. (2) The source apportionment results show that secondary inorganic aerosols, traffic emissions, and secondary organic aerosols were the dominant contributors to PM2.5 in Zibo. The total contributions of the three sources accounted for approximately 80% of PM2.5. Strict control of precursor emissions, such as SO2, NOx, and VOCs, is a good way to reduce the PM2.5 pollution resulting from secondary aerosols. Traffic emissions contributed more than 20% to PM2.5 and were the main source of heavy metals, according to combined Pearson correlation analysis of elements, so proper control measures should be taken. (3) The non-carcinogenic risks from crustal dust were 2.23 and 1.15 in urban and suburban areas, respectively, exceeding the safe limit for children and the limit for adults at the urban site. Meanwhile, the carcinogenic risks were all below the safe limit. However, non-carcinogenic and carcinogenic risks from traffic emissions were just below the safe limit. These results may help local authorities to design strategies for reducing PM2.5 pollution and limiting health impacts in Zibo during haze episodes; for example, by traffic control, by limiting or preventing outdoor activities, and by wearing masks.

Author Contributions

Conceptualization, W.B.; methodology, X.Z. and W.B.; software, X.Z.; validation, L.G. and W.B.; formal analysis, W.B.; investigation, B.Y.; resources, W.B.; data curation, W.B.; writing—original draft preparation, W.B.; writing—review and editing, X.W.; visualization, W.B.; supervision, W.Z. and X.W.; project administration, W.B.; funding acquisition, W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Program for Key Issues in Air Pollution Control (Granted No. DQGG2021202 and DQGG0107-20).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors gratefully acknowledge the support of the National Research Program for Key Issues in Air Pollution Control (DQGG2021202, DQGG0107-20).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, I.; Lee, K.; Lee, S.; Kim, S.D. Characteristics and health effects of PM2.5 emissions from various sources in Gwangju, South Korea. Sci. Total Environ. 2019, 696, 133890. [Google Scholar] [CrossRef] [PubMed]
  2. Khan, M.F.; Latif, M.T.; Saw, W.H.; Amil, N.; Nadzir, M.S.M.; Sahani, M.; Tahir, N.M.; Chung, J.X. Fine particulate matter in the tropical environment: Monsoonal effects, source apportionment, and health risk assessment. Atmos. Chem. Phys. 2016, 16, 597–617. [Google Scholar] [CrossRef] [Green Version]
  3. Wang, Y.; Wang, M.; Zhang, R.; Ghan, S.J.; Lin, Y.; Hu, J.; Pan, B.; Levy, M.; Jiang, J.H.; Molina, M.J. Assessing the effects of anthropogenic aerosols on Pacific storm track using a multiscale global climate model. Proc. Natl. Acad. Sci. USA 2014, 111, 6894–6899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Xu, H.-H.; Pu, J.-J.; He, J.; Liu, J.; Qi, B.; Du, R.-G. Characteristics of Atmospheric Compositions in the Background Area of Yangtze River Delta during Heavy Air Pollution Episode. Adv. Meteorol. 2016, 2016, 7240486. [Google Scholar] [CrossRef]
  5. Tao, J.; Zhang, L.; Ho, K.; Zhang, R.; Lin, Z.; Zhang, Z.; Lin, M.; Cao, J.; Liu, S.; Wang, G. Impact of PM2.5 chemical compositions on aerosol light scattering in Guangzhou—the largest megacity in South China. Atmos. Res. 2014, 135–136, 48–58. [Google Scholar] [CrossRef]
  6. Lelieveld, J.; Evans, J.S.; Fnais, M.; Giannadaki, D.; Pozzer, A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 2015, 525, 367–371. [Google Scholar] [CrossRef]
  7. Carey, I.M.; Anderson, H.R.; Atkinson, R.W.; Beevers, S.D.; Cook, D.G.; Strachan, D.P.; Dajnak, D.; Gulliver, J.; Kelly, F.J. Are noise and air pollution related to the incidence of dementia? A cohort study in London, England. BMJ Open 2018, 8, e022404. [Google Scholar] [CrossRef] [Green Version]
  8. Pope, C.A.; Dockery, D.W. Health Effects of Fine Particulate Air Pollution: Lines that Connect. J. Air Waste Manag. Assoc. 2012, 56, 709–742. [Google Scholar] [CrossRef]
  9. Chu, C.; Zhang, H.; Cui, S.; Han, B.; Zhou, L.; Zhang, N.; Su, X.; Niu, Y.; Chen, W.; Chen, R.; et al. Ambient PM2.5 caused depressive-like responses through Nrf2/NLRP3 signaling pathway modulating inflammation. J. Hazard Mater. 2019, 369, 180–190. [Google Scholar] [CrossRef]
  10. Manisalidis, I.; Stavropoulou, E.; Stavropoulos, A.; Bezirtzoglou, E. Environmental and Health Impacts of Air Pollution: A Review. Front. Public Health 2020, 8, 14. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, G.; Li, Y.; Zhou, J.; Xu, J.; Yang, B. PM2.5 deregulated microRNA and inflammatory microenvironment in lung injury. Env. Toxicol Pharm. 2022, 91, 103832. [Google Scholar] [CrossRef] [PubMed]
  12. Hao, G.; Zuo, L.; Xiong, P.; Chen, L.; Liang, X.; Jing, C. Associations of PM2.5 and road traffic noise with mental health: Evidence from UK Biobank. Environ. Res. 2022, 207, 112221. [Google Scholar] [CrossRef] [PubMed]
  13. Tian, S.L.; Pan, Y.P.; Wang, Y.S. Size-resolved source apportionment of particulate matter in urban Beijing during haze and non-haze episodes. Atmos. Chem. Phys. 2016, 16, 1–19. [Google Scholar] [CrossRef] [Green Version]
  14. Yu, Q.; Chen, J.; Qin, W.; Cheng, S.; Zhang, Y.; Ahmad, M.; Ouyang, W. Characteristics and secondary formation of water-soluble organic acids in PM1, PM2.5 and PM10 in Beijing during haze episodes. Sci. Total Environ. 2019, 669, 175–184. [Google Scholar] [CrossRef] [PubMed]
  15. Park, M.; Joo, H.S.; Lee, K.; Jang, M.; Kim, S.D.; Kim, I.; Borlaza, L.J.S.; Lim, H.; Shin, H.; Chung, K.H.; et al. Differential toxicities of fine particulate matters from various sources. Sci. Rep. 2018, 8, 17007. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, L.; Liu, G.; Lin, Z.; Wang, Y.; He, H.; Liu, T.; Kamp, D.W. Pro-inflammatory response and oxidative stress induced by specific components in ambient particulate matter in human bronchial epithelial cells. Environ. Toxicol. 2016, 31, 923–936. [Google Scholar] [CrossRef]
  17. Bell, M.L.; Ebisu, K.; Leaderer, B.P.; Gent, J.F.; Lee, H.J.; Koutrakis, P.; Wang, Y.; Dominici, F.; Peng, R.D. Associations of PM2.5 constituents and sources with hospital admissions: Analysis of four counties in Connecticut and Massachusetts (USA) for persons >/= 65 years of age. Environ. Health Perspect 2014, 122, 138–144. [Google Scholar] [CrossRef] [Green Version]
  18. Li, K.; Wang, X.; Li, L.; Wang, J.; Liu, Y.; Cheng, X.; Xu, B.; Wang, X.; Yan, P.; Li, S.; et al. Large variability of O3-precursor relationship during severe ozone polluted period in an industry-driven cluster city (Zibo) of North China Plain. J. Clean. Prod. 2021, 316, 128252. [Google Scholar] [CrossRef]
  19. Li, H.; Ma, Y.; Duan, F.; He, K.; Zhu, L.; Huang, T.; Kimoto, T.; Ma, X.; Ma, T.; Xu, L. Typical winter haze pollution in Zibo, an industrial city in China: Characteristics, secondary formation, and regional contribution. Environ. Pollut. 2017, 229, 339–349. [Google Scholar] [CrossRef]
  20. Wu, L.; Li, M.; Zhang, X.; Wang, X.; Geng, C.; Zhao, X.; Yang, W.; Yin, B. Spatial-Temporal Characteristics of Ground Ozone in Zibo City from 2016 to 2019. Res. Environ. Sci. 2021, 34, 9. [Google Scholar] [CrossRef]
  21. Wu, L.; Li, L.; Zhang, X.; Wang, X.; Ma, Y.; Geng, C.; Xu, B. Characteristics of Atmospheric Particulate Matters during Winter Sandstorm Period in 2018 in Zibo City. Res. Environ. Sci. 2021, 34, 11. [Google Scholar] [CrossRef]
  22. Gan, H.; Xu, B.; Zhang, X.; Geng, C.; Yang, W.; Wang, X.; Yin, B. Characteristics and Health Risk Assessment of Atmospheric VOCs in a Chemical Industrial Park of Zibo in Summer. Res. Environ. Sci. 2022, 35, 10. [Google Scholar] [CrossRef]
  23. Dai, J.; Yu, C.; Zhang, M. Distribution Characteristics and Sources of Heavy Metals in Urban Atmospheric Particulate Matter in Zibo City. J. Jilin Univ. 2018, 48, 1201–1211. [Google Scholar] [CrossRef]
  24. Liu, X.; Wang, D.; Meng, C. Pollution Status and Source Analysis of PM2.5 in Zibo City, 2019. Mod. Prev. Med. 2021, 48, 23–26+43. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=XDYF202101006&DbName=DKFX2021 (accessed on 10 January 2021).
  25. Luo, Y.; Zhou, X.; Zhang, J.; Xiao, Y.; Wang, Z.; Zhou, Y.; Wang, W. PM2.5 pollution in a petrochemical industry city of northern China: Seasonal variation and source apportionment. Atmos. Res. 2018, 212, 285–295. [Google Scholar] [CrossRef]
  26. HJ 656-2013; Technical Specifications for Gravimetric Measurement Methods for PM2.5 in Ambient Air; Ministry of Environmental Protection (MEP) of the People’s Republic of China: Beijing, China, 2013.
  27. Li, M.; Wu, L.; Zhang, X.; Wang, X.; Bai, W.; Ming, J.; Geng, C.; Yang, W. Comparison of PM2.5 Chemical Compositions during Haze and Non-haze Days in a Heavy Industrial City in North China. Aerosol Air Qual. Res. 2020, 20, 1950–1960. [Google Scholar] [CrossRef]
  28. Yao, L.; Yang, L.; Yuan, Q.; Yan, C.; Dong, C.; Meng, C.; Sui, X.; Yang, F.; Lu, Y.; Wang, W. Sources apportionment of PM2.5 in a background site in the North China Plain. Sci. Total Environ. 2016, 541, 590–598. [Google Scholar] [CrossRef]
  29. Kong, S.; Han, B.; Bai, Z.; Chen, L.; Shi, J.; Xu, Z. Receptor modeling of PM2.5, PM10 and TSP in different seasons and long-range transport analysis at a coastal site of Tianjin, China. Sci. Total Environ. 2010, 408, 4681–4694. [Google Scholar] [CrossRef]
  30. Norris, G.; Duvall, R.; Brown, S.; Bai, S. EPA Positive Matrix Factorization (PMF) 5.0 Fundamentals and User Guide; U.S. Environmental Protection Agency: Washington, DC, USA, 2014.
  31. Chen, C.R.; Lai, H.C.; Liao, M.I.; Hsiao, M.C.; Ma, H.W. Health risk assessment of trace elements of ambient PM2.5 under monsoon patterns. Chemosphere 2021, 264, 128462. [Google Scholar] [CrossRef]
  32. Zhao, X.; Wang, X.; Zhu, Y.; Yang, W.; Ren, L.; Bai, Z. Characterization of Chemical Composition Mass Balance and Source Appointment of Ambient PM2.5 in Xinzhou City. Chin. J. Environ. Eng. 2017, 11, 4660–4668. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=HJJZ201708030&DbName=CJFQ2017 (accessed on 5 August 2017).
  33. Gupta, A.K.; Karar, K.; Srivastava, A. Chemical mass balance source apportionment of PM10 and TSP in residential and industrial sites of an urban region of Kolkata, India. J. Hazard. Mater. 2007, 142, 279–287. [Google Scholar] [CrossRef] [PubMed]
  34. Ali, M.U.; Liu, G.; Yousaf, B.; Abbas, Q.; Ullah, H.; Munir, M.A.M.; Fu, B. Pollution characteristics and human health risks of potentially (eco)toxic elements (PTEs) in road dust from metropolitan area of Hefei, China. Chemosphere 2017, 181, 111–121. [Google Scholar] [CrossRef] [PubMed]
  35. EPA. Risk Assessment Guidance for Superfund (RAGS): Part F; U.S. Environmental Protection Agency: Washington, DC, USA, 2011.
  36. Zhou, Q.; Zheng, N.; Liu, J.; Wang, Y.; Sun, C.; Liu, Q.; Wang, H.; Zhang, J. Residents health risk of Pb, Cd and Cu exposure to street dust based on different particle sizes around zinc smelting plant, Northeast of China. Environ. Geochem. Health 2015, 37, 207–220. [Google Scholar] [CrossRef] [PubMed]
  37. Ferreira-Baptista, L.; De Miguel, E. Geochemistry and Risk Assessment of Street Dust in Luanda, Angola: A Tropical Urban Environment. Atmos. Environ. 2005, 39, 4501–4512. [Google Scholar] [CrossRef] [Green Version]
  38. EPA. IRIS Assessments. Available online: https://iris.epa.gov/AtoZ/?list_type=alpha (accessed on 1 January 2022).
  39. Dahmardeh Behrooz, R.; Kaskaoutis, D.G.; Grivas, G.; Mihalopoulos, N. Human health risk assessment for toxic elements in the extreme ambient dust conditions observed in Sistan, Iran. Chemosphere 2021, 262, 127835. [Google Scholar] [CrossRef]
  40. Hsu, C.Y.; Chiang, H.C.; Lin, S.L.; Chen, M.J.; Lin, T.Y.; Chen, Y.C. Elemental characterization and source apportionment of PM10 and PM2.5 in the western coastal area of central Taiwan. Sci. Total Environ. 2016, 541, 1139–1150. [Google Scholar] [CrossRef] [PubMed]
  41. Zhao, X.; Duan, X.; Wang, B. Environmental Exposure Related Activity Patterns Survey of Chinese Population (Children); Ministry of Environmental Protection (MEP) of the People’s Republic of China: Beijing, China, 2016.
  42. Duan, X. Report of Environmental Exposure Related Activity Patterns Research of Chinese Population; Ministry of Environmental Protection (MEP) of the People’s Republic of China: Beijing, China, 2013.
  43. EPA, U.S. Exposure Factors Handbook 2011 Edition (Final Report); U.S. Environmental Protection Agency: Washington, DC, USA, 2011.
  44. GB 3095-2012; Ambient Air Quality Standards. Ministry of Environmental Protection (MEP) of the People’s Republic of China: Beijing, China, 2012.
  45. WHO. Ambient (Outdoor) Air Pollution. Available online: https://www.who.int/news-room/fact-sheets/detail/ambient-(outdoor)-air-quality-and-health (accessed on 1 January 2022).
  46. Li, Q.; Zhang, H.; Jin, X.; Cai, X.; Song, Y. Mechanism of haze pollution in summer and its difference with winter in the North China Plain. Sci. Total Environ. 2022, 806, 150625. [Google Scholar] [CrossRef]
  47. Peng, J.; Hu, M.; Shang, D.; Wu, Z.; Du, Z.; Tan, T.; Wang, Y.; Zhang, F.; Zhang, R. Explosive Secondary Aerosol Formation during Severe Haze in the North China Plain. Environ. Sci. Technol. 2021, 55, 2189–2207. [Google Scholar] [CrossRef]
  48. Arimoto, R.; Duce, R.A.; Savoie, D.L.; Prospero, J.M.; Talbot, R.; Cullen, J.D.; Tomza, U.; Lewis, N.F.; Ray, B.J. Relationships among aerosol constituents from Asia and the North Pacific during PEM-West A. J. Geophys. Res. Atmos. 1996, 101, 2011–2023. [Google Scholar] [CrossRef]
  49. Huebert, B.J.; Wang, M.X.; Wei-Xiu, L. Atmospheric nitrate, sulfate, ammonium and calcium concentrations in China. Tellus B 1988, 40, 260–269. [Google Scholar] [CrossRef] [Green Version]
  50. Wang, X.; Wei, W.; Cheng, S.; Yao, S.; Zhang, H.; Zhang, C. Characteristics of PM2.5 and SNA components and meteorological factors impact on air pollution through 2013–2017 in Beijing, China. Atmos. Pollut. Res. 2019, 10, 1976–1984. [Google Scholar] [CrossRef]
  51. Ye, Z.; Liu, J.; Gu, A.; Feng, F.; Liu, Y.; Bi, C.; Xu, J.; Li, L.; Chen, H.; Chen, Y.; et al. Chemical characterization of fine particulate matter in Changzhou, China, and source apportionment with offline aerosol mass spectrometry. Atmos. Chem. Phys. 2017, 17, 2573–2592. [Google Scholar] [CrossRef] [Green Version]
  52. Xiong, Y.; Zhou, J.; Schauer, J.J.; Yu, W.; Hu, Y. Seasonal and spatial differences in source contributions to PM2.5 in Wuhan, China. Sci. Total Environ. 2016, 577, 155–165. [Google Scholar] [CrossRef] [PubMed]
  53. Jiang, N.; Yin, S.; Guo, Y.; Li, J.; Kang, P.; Zhang, R.; Tang, X. Characteristics of mass concentration, chemical composition, source apportionment of PM 2.5 and PM 10 and health risk assessment in the emerging megacity in China. Atmos. Pollut. Res. 2018, 9, 309–321. [Google Scholar] [CrossRef]
  54. Li, H.; Cheng, J.; Zhang, Q.; Zheng, B.; Zhang, Y.; Zheng, G.; He, K. Rapid transition in winter aerosol composition in Beijing from 2014 to 2017: Response to clean air actions. Atmos. Chem. Phys. 2019, 19, 11485–11499. [Google Scholar] [CrossRef] [Green Version]
  55. Liu, Y.; Xing, J.; Wang, S.; Fu, X.; Zheng, H. Source-specific speciation profiles of PM2.5 for heavy metals and their anthropogenic emissions in China. Environ. Pollut. 2018, 239, 544–553. [Google Scholar] [CrossRef] [Green Version]
  56. Hulskotte, J.H.J.; Roskam, G.D.; Denier van der Gon, H.A.C. Elemental composition of current automotive braking materials and derived air emission factors. Atmos. Environ. 2014, 99, 436–445. [Google Scholar] [CrossRef]
  57. De Silva, S.; Ball, A.S.; Huynh, T.; Reichman, S.M. Metal accumulation in roadside soil in Melbourne, Australia: Effect of road age, traffic density and vehicular speed. Environ. Pollut. 2016, 208, 102–109. [Google Scholar] [CrossRef]
  58. Jiang, P.; Zhong, X.; Li, L. On-road vehicle emission inventory and its spatio-temporal variations in North China Plain. Environ. Pollut. 2020, 267, 115639. [Google Scholar] [CrossRef]
  59. Cong, Z.; Kang, S.; Zhang, Y.; Li, X. Atmospheric wet deposition of trace elements to central Tibetan Plateau. Appl. Geochem. 2010, 25, 1415–1421. [Google Scholar] [CrossRef]
  60. Zhao, S.; Tian, H.; Luo, L.; Liu, H.; Wu, B.; Liu, S.; Bai, X.; Liu, W.; Liu, X.; Wu, Y.; et al. Temporal variation characteristics and source apportionment of metal elements in PM2.5 in urban Beijing during 2018–2019. Environ. Pollut. 2021, 268, 115856. [Google Scholar] [CrossRef] [PubMed]
  61. Kuhn, N.J. Erodibility of soil and organic matter: Independence of organic matter resistance to interrill erosion. Earth Surf. Processes Landf. 2007, 32, 794–802. [Google Scholar] [CrossRef]
  62. Han, J.; Han, B.; Li, P.; Kong, S.; Bai, Z.; Han, D.; Dou, X.; Zhao, X. Chemical Characterizations of PM10 Profiles for Major Emission Sources in Xining, Northwestern China. Aerosol Air Qual. Res. 2014, 14, 1017–1027. [Google Scholar] [CrossRef] [Green Version]
  63. Dhinasekaran, D.; Soundharraj, P.; Jagannathan, M.; Rajendran, A.R.; Rajendran, S. Hybrid ZnO nanostructures modified graphite electrode as an efficient urea sensor for environmental pollution monitoring. Chemosphere 2022, 296, 133918. [Google Scholar] [CrossRef] [PubMed]
  64. Papučová, I.; Plško, A.; Pagáčová, J. Corrosion Layers on the Glass Surface of Na2O-K2O-CaO-ZrO2-SiO2 System. Procedia Eng. 2016, 136, 287–293. [Google Scholar] [CrossRef] [Green Version]
  65. Chen, T.; Chu, B.; Ma, Q.; Zhang, P.; Liu, J.; He, H. Effect of relative humidity on SOA formation from aromatic hydrocarbons: Implications from the evolution of gas- and particle-phase species. Sci. Total Environ. 2021, 773, 145015. [Google Scholar] [CrossRef]
  66. Gao, J.; Wang, K.; Wang, Y.; Liu, S.; Zhu, C.; Hao, J.; Liu, H.; Hua, S.; Tian, H. Temporal-spatial characteristics and source apportionment of PM2.5 as well as its associated chemical species in the Beijing-Tianjin-Hebei region of China. Environ. Pollut. 2018, 233, 714–724. [Google Scholar] [CrossRef]
  67. Housing and Urban-Rural Development Bureau. Implementation Plan for the Special Treatment of Dust Pollution in the Construction Field of Zibo City. Available online: http://js.zibo.gov.cn/gongkai/channel_c_5f9fa491ab327f36e4c1305d_n_1605682494.016/doc_5fb4744399d584448a0dcf0a.html (accessed on 1 January 2022).
  68. Liu, X.; Ouyang, W.; Shu, Y.; Tian, Y.; Feng, Y.; Zhang, T.; Chen, W. Incorporating bioaccessibility into health risk assessment of heavy metals in particulate matter originated from different sources of atmospheric pollution. Environ. Pollut. 2019, 254, 113113. [Google Scholar] [CrossRef]
Figure 1. Geographical information of: (a) North China Plain, and (b) sampling sites of Nanding (urban site) and Lishan (suburban site) in Zibo.
Figure 1. Geographical information of: (a) North China Plain, and (b) sampling sites of Nanding (urban site) and Lishan (suburban site) in Zibo.
Ijerph 19 05443 g001
Figure 2. Sampled concentrations of PM2.5 in urban (red) and suburban areas (blue) during each season. Each subfigure shows the daily mean values (dots), median (horizontal line), central 50% of data (25–75th percentile; box), and the central 90% of data (5–95th percentile; whiskers).
Figure 2. Sampled concentrations of PM2.5 in urban (red) and suburban areas (blue) during each season. Each subfigure shows the daily mean values (dots), median (horizontal line), central 50% of data (25–75th percentile; box), and the central 90% of data (5–95th percentile; whiskers).
Ijerph 19 05443 g002
Figure 3. Variation of SNA during haze episode (left: urban; right: suburban).
Figure 3. Variation of SNA during haze episode (left: urban; right: suburban).
Ijerph 19 05443 g003
Figure 4. Concentrations of OC and EC and OC/EC ratio in urban (n = 45) and suburban (n = 45) areas during the sampling period.
Figure 4. Concentrations of OC and EC and OC/EC ratio in urban (n = 45) and suburban (n = 45) areas during the sampling period.
Ijerph 19 05443 g004
Figure 5. Boxplots showing the descriptive statistics of the sampled concentrations of PM2.5-metals. Each subfigure shows the mean (dot), median (horizontal line), central 50% of data (25–75th percentile; box), and central 90% of data (5–95th percentile; whiskers).
Figure 5. Boxplots showing the descriptive statistics of the sampled concentrations of PM2.5-metals. Each subfigure shows the mean (dot), median (horizontal line), central 50% of data (25–75th percentile; box), and central 90% of data (5–95th percentile; whiskers).
Ijerph 19 05443 g005
Figure 6. Daily variation series of PM2.5 metals (urban (ac); suburban (df)).
Figure 6. Daily variation series of PM2.5 metals (urban (ac); suburban (df)).
Ijerph 19 05443 g006
Figure 7. Pearson correlation of elements in urban (n = 45) and suburban (n = 45) samples.
Figure 7. Pearson correlation of elements in urban (n = 45) and suburban (n = 45) samples.
Ijerph 19 05443 g007
Figure 8. PM2.5 source profiles (weight percent by mass) in this study.
Figure 8. PM2.5 source profiles (weight percent by mass) in this study.
Ijerph 19 05443 g008
Table 1. Chemical mass balance model performance test results during the study period.
Table 1. Chemical mass balance model performance test results during the study period.
Sampling
Duration
Sampling
Site
R2CHI2DFPercentage of Explained Mass to Sample Total Mass
Winterurban (n = 15)0.940.39583.9
suburban (n = 15)0.890.82583.0
Springurban (n = 10)0.881.06882.7
suburban (n = 10)0.930.69788.4
Summerurban (n = 10)0.941.57584.1
suburban (n = 10)0.831.48880.7
Autumnurban (n = 10)0.980.24582.3
suburban (n = 10)0.940.65584.4
Annualurban (n = 10)0.930.34582.8
suburban (n = 10)0.930.39587.3
Note: The criteria for acceptable CMB results included the square regression coefficient, R2 ≥ 0.8; the sum of residual square value, CHI2 ≤ 4; the degree of the freedom, DF ≥ 5; and the percentage of explained mass to sample total mass, ranging from 80 to 120%.
Table 2. Values of exposure parameters for children and adults.
Table 2. Values of exposure parameters for children and adults.
ParameterDefinitionUnitValueReference
ChildrenAdult
CExposure-point concentrationmg kg−1Present study a
InhRInhalation ratem3 day−19.016.1Refs. [41,42]
EFExposure frequencyday year−1350350Ref. [34]
EDExposure durationyear624Ref. [34]
PEFParticle emission factorm3 kg−11.36 × 1091.36 × 109Ref. [34]
BWAverage body weightkg20.565.0Refs. [41,42]
ATAveraging timeday365 × ED b
365 × 78 c
365 × ED b
365 × 78 c
Ref. [43]
Note: a By using the concentration of the individual metal (profile) of each source for C. b For non-carcinogenic risk assessment. c For carcinogenic risk assessment.
Table 3. Values of RfD and SF used in the present study [37,43].
Table 3. Values of RfD and SF used in the present study [37,43].
ElementRfDSF
As 1.51 × 101
Ba1.43 × 10−4
Cd 6.30
Co5.71 × 10−69.80
Cr2.86 × 10−54.20 × 101
Mn1.43 × 10−5
Ni 8.40 × 10−1
Table 4. Seasonal average WSII concentrations and PM2.5 ratios in urban and suburban areas.
Table 4. Seasonal average WSII concentrations and PM2.5 ratios in urban and suburban areas.
Winter n = 15Spring n = 10Summer n = 10Autumn n = 10Annual n = 45
UrbanMeanSDMeanSDMeanSDMeanSDMeanSD
PM2.5 (μg/m3)112±7561±2332±897±3979±58
WSIIs (μg/m3)59.1±46.831.6±14.114.8±4.954.6±23.342.2±35.0
Cl (μg/m3)4.5±2.41.4±1.20.1±0.11.8±0.92.2±2.3
Ratios
OC/EC2.7±0.32.6±1.22.0±0.31.9±0.22.3±0.7
NO3/SO42−2.5±0.92.2±1.20.5±0.33.9±1.62.3±1.6
WSIIs/PM2.5 (%)49±1051±746±556±550±8
SNA/PM2.5 (%)42±1146±639±651±844±9
OC/PM2.5 (%)19±413±516±410±215±5
EC/PM2.5 (%)7±26±38±25±17±2
Winter n = 15Spring n = 10Summer n = 10Autumn n = 10Annual n = 45
SuburbanMeanSDMeanSDMeanSDMeanSDMeanSD
PM2.5 (μg/m3)109±6058±2232±789±3076±49
WSIIs (μg/m3)59.2±39.337.2±19.317.3±4.054.7±22.144.0±31.3
Cl (μg/m3)3.3±1.21.8±2.10.4±0.21.3±0.71.9±1.7
Ratios
OC/EC2.5±0.53.4±1.31.9±0.42.0±0.72.4±1.0
NO3/SO42−2.5±0.82.7±1.30.8±0.44.2±1.72.5±1.6
WSIIs/PM2.5 (%)52±962±1654±760±656±11
SNA/PM2.5 (%)46±1153±1048±856±750±10
OC/PM2.5 (%)20±615±313±411±315±6
EC/PM2.5 (%)8±35±17±26±27±3
Table 5. Percentage contributions of each source type to PM2.5 concentrations (%) in each season.
Table 5. Percentage contributions of each source type to PM2.5 concentrations (%) in each season.
Sampling
Duration
Sampling
Site
Percentage Contribution of Each Source Type
CDCBGPASGMISTESIASOAOthers
8–22 January
(winter)
urban6.42.92.71.71.51.6126.043.812.60.8
suburban6.12.61.31.11.21.422.545.413.45.0
16–25 April (spring)urban13.62.12.20.81.21.918.842.28.29.0
suburban7.21.52.30.30.60.914.552.510.99.4
30 July–8 August
(summer)
urban10.11.81.81.61.82.227.834.89.98.2
suburban6.61.21.11.11.01.325.735.86.819.4
16–25 October
(autumn)
urban5.81.91.91.31.51.524.446.35.59.9
suburban4.61.92.10.41.01.423.151.47.36.7
Annualurban7.82.32.21.41.80.924.043.39.86.4
suburban5.61.61.90.71.01.220.549.410.67.5
Table 6. Source-specific non-carcinogenic and carcinogenic risks of heavy metals in ambient PM2.5 collected from urban and suburban areas in Zibo.
Table 6. Source-specific non-carcinogenic and carcinogenic risks of heavy metals in ambient PM2.5 collected from urban and suburban areas in Zibo.
Non-Carcinogenic Risk: Hazard Quotient (HQ) of Each Heavy Metal and Hazard Index (HI)
AdultsChildren
CDCBGPGMTECDCBGPGMTE
Urban site
HQCr1.72 × 10−32.74 × 10−34.46 × 10−30.004.38 × 10−23.05 × 10−34.85 × 10−37.90 × 10−30.007.77 × 10−2
HQMn1.202.18 × 10−12.83 × 10−10.003.70 × 10−12.123.86 × 10−15.02 × 10−10.006.57 × 10−1
HQCo7.36 × 10−36.55 × 10−37.97 × 10−30.005.12 × 10−31.30 × 10−21.16 × 10−21.41 × 10−20.009.08 × 10−3
HQBa4.94 × 10−21.31 × 10−33.92 × 10−30.001.92 × 10−28.76 × 10−22.32 × 10−36.95 × 10−30.003.40 × 10−2
HI1.262.28 × 10−13.00 × 10−10.004.39 × 10−12.234.05 × 10−15.31 × 10−10.007.77 × 10−1
Suburban site
HQCr8.28 × 10−42.99 × 10−33.06 × 10−33.33 × 10−34.54 × 10−21.47 × 10−35.31 × 10−35.42 × 10−35.91 × 10−38.05 × 10−2
HQMn6.36 × 10−12.70 × 10−12.22 × 10−19.53 × 10−24.29 × 10−11.134.79 × 10−13.94 × 10−11.69 × 10−17.60 × 10−1
HQCo2.17 × 10−34.37 × 10−33.35 × 10−38.69 × 10−43.26 × 10−33.85 × 10−37.75 × 10−35.95 × 10−31.54 × 10−35.78 × 10−3
HQBa8.04 × 10−34.42 × 10−48.97 × 10−41.47 × 10−46.72 × 10−31.43 × 10−27.83 × 10−41.59 × 10−32.61 × 10−41.19 × 10−2
HI6.47 × 10−12.78 × 10−12.30 × 10−19.97 × 10−24.84 × 10−11.154.93 × 10−14.07 × 10−11.77 × 10−18.59 × 10−1
Carcinogenic Risk: Carcinogenic Risk (RIi) of Each Heavy Metal and Total Carcinogenic Risk (RI)
Urban SiteSuburban Site
CDCBGPGMTECDCBGPGMTE
RICr2.62 × 10−94.17 × 10−96.79 × 10−90.006.68 × 10−81.26 × 10−94.56 × 10−94.66 × 10−95.08 × 10−96.92 × 10−8
RICo5.22 × 10−104.65 × 10−105.66 × 10−100.003.64 × 10−101.54 × 10−103.10 × 10−102.38 × 10−106.17 × 10−112.31 × 10−10
RINi3.18 × 10−101.29 × 10−91.35 × 10−90.004.98 × 10−101.50 × 10−101.37 × 10−99.04 × 10−103.14 × 10−105.05 × 10−10
RIAs2.66 × 10−105.46 × 10−107.23 × 10−100.009.15 × 10−88.50 × 10−113.91 × 10−103.29 × 10−106.06 × 10−106.30 × 10−8
RICd3.45 × 10−102.30 × 10−101.08 × 10−80.009.21 × 10−101.11 × 10−101.34 × 10−104.08 × 10−93.12 × 10−104.90 × 10−10
RI4.07 × 10−96.70 × 10−92.03 × 10−80.001.60 × 10−71.76 × 10−96.77 × 10−91.02 × 10−86.37 × 10−91.33 × 10−7
Note: The values larger than the threshold value (HQ, HI, 1.00; RI, 1.00 × 10−6) were marked in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bai, W.; Zhao, X.; Yin, B.; Guo, L.; Zhang, W.; Wang, X.; Yang, W. Characteristics of PM2.5 in an Industrial City of Northern China: Mass Concentrations, Chemical Composition, Source Apportionment, and Health Risk Assessment. Int. J. Environ. Res. Public Health 2022, 19, 5443. https://doi.org/10.3390/ijerph19095443

AMA Style

Bai W, Zhao X, Yin B, Guo L, Zhang W, Wang X, Yang W. Characteristics of PM2.5 in an Industrial City of Northern China: Mass Concentrations, Chemical Composition, Source Apportionment, and Health Risk Assessment. International Journal of Environmental Research and Public Health. 2022; 19(9):5443. https://doi.org/10.3390/ijerph19095443

Chicago/Turabian Style

Bai, Wenyu, Xueyan Zhao, Baohui Yin, Liyao Guo, Wenge Zhang, Xinhua Wang, and Wen Yang. 2022. "Characteristics of PM2.5 in an Industrial City of Northern China: Mass Concentrations, Chemical Composition, Source Apportionment, and Health Risk Assessment" International Journal of Environmental Research and Public Health 19, no. 9: 5443. https://doi.org/10.3390/ijerph19095443

APA Style

Bai, W., Zhao, X., Yin, B., Guo, L., Zhang, W., Wang, X., & Yang, W. (2022). Characteristics of PM2.5 in an Industrial City of Northern China: Mass Concentrations, Chemical Composition, Source Apportionment, and Health Risk Assessment. International Journal of Environmental Research and Public Health, 19(9), 5443. https://doi.org/10.3390/ijerph19095443

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