Next Article in Journal
Pilot Ecological Momentary Assessment Study of Subjective and Contextual Factors Surrounding E-Cigarette and Combustible Tobacco Product Use among Young Adults
Next Article in Special Issue
A Comprehensive 2018-Based Vehicle Emission Inventory and Its Spatial–Temporal Characteristics in the Central Liaoning Urban Agglomeration, China
Previous Article in Journal
Smoking a Dangerous Addiction: A Systematic Review on an Underrated Risk Factor for Oral Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Concentrations, Source Characteristics, and Health Risk Assessment of Toxic Heavy Metals in PM2.5 in a Plateau City (Kunming) in Southwest China

1
Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650500, China
2
Faculty of Environmental Science and Engineering, Kunming University of Science and Technology, Kunming 650500, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(21), 11004; https://doi.org/10.3390/ijerph182111004
Submission received: 2 July 2021 / Revised: 14 October 2021 / Accepted: 18 October 2021 / Published: 20 October 2021

Abstract

:
To explore the mass concentration levels and health risks of heavy metals in the air in dense traffic environments, PM2.5 samples were collected at three sites in the city of Kunming in April and October 2013, and January and May 2014. Ten heavy metals––V, Cr, Mn, Co, Ni, Cu, Zn, As, Cd and Pb––were analyzed by ICP–MS, and the results showed PM2.5 concentrations significantly higher in spring and winter than in summer and autumn, especially for Zn and Pb. The concentration of heavy metals on working days is significantly higher, indicating that vehicle emissions are significant contributors. An enrichment factor analysis showed that Cr, Mn, Ni, Cu, Zn, As, Cd and Pb come mainly from anthropogenic sources, while V and Co may be both anthropogenic and natural. The correlation and principal component analysis (PCA) showed that Ni, Cu, Zn, Cd and Pb mainly come from vehicles emissions and metallurgical industries; Cr and Mn, from vehicles emissions and road dust; and As, mainly from coal combustion. The health risk assessment shows that the non-carcinogenic risk thresholds of the heavy metals in PM2.5 to children and adult men and women are all less than 1. The carcinogenic risk of Cr for men and women in traffic-intensive areas exceeds 10−4, reaching 1.64 × 10−4 and 1.4 × 10−4, respectively.

1. Introduction

In recent decades, the number of motor vehicles and total energy consumption have increased in China, as a result of rapid economic development and urbanization. Today, atmospheric particulate matter (PM) is one of the most significant air contaminants [1,2,3]. PM less than 2.5 µm (PM2.5) comprises heavy metals, bacteria, toxins and carcinogens due to its small particle size and large specific surface area [4,5]; consequently, it has received extensive attention from researchers. Studies have shown the adverse health effects from exposure to PM and heavy metals [6,7,8,9], and epidemiological studies over the past few decades identified a close relationship between increased atmospheric heavy metals and increased mortality and morbidity [10,11,12].
The heavy metals come from natural and anthropogenic sources, including coal combustion, traffic emissions, and smelting. Because of their high toxicity, persistence and bioaccumulation, research into the health effects of long or short-term exposure to traffic-related air pollutants has become the focus of attention [13,14,15,16]. Fang et al. [17] proposed that heavy metals in road dust are not easy to remove, and may diffuse into the atmosphere and persist in the environment for a long time. He et al. [18] believe that large and medium-sized cities show obvious effects of vehicle pollution, and a large number of studies found high concentrations of Cr, Cu, Zn, Cd and Pb. These metals can cause respiratory and cardiovascular disease, and increase the risk of cancer [19,20,21,22,23,24]. Recent studies also showed that exposure to As, Cd, Hg and Pb in early or late pregnancy is associated with childhood asthma [25].
According to statistics from 2014 to 2019, deaths from respiratory diseases in Yunnan Province were 115.26 per 100,000 people, much higher than the national rate of 69.47/100,000 [26,27,28,29]. We conducted this study to find out whether heavy metals of PM2.5 had an impact on the mortality rate of people in Yunnan. However, most atmospheric pollution studies in China focus on the Beijing–Tianjin–Hebei region and coastal areas [30,31]; only a few have investigated PM2.5 pollution in Yunnan, and data on heavy-metal pollution is lacking. Therefore, we collected PM2.5 samples at three sites in Kunming from 2013 and 2014 to explore the contribution of heavy metals to environmental pollution and human health.
Kunming is in the middle of the Yunnan–Guizhou Plateau in southwestern China (24°23′–26°22′ N, 102°10′–103°40′ E) and the capital of Yunnan Province. Its population is 6.95 million, one of the largest cities on the plateau. Kunming is to the northeast of the Dianchi Basin and surrounded by mountains on three sides, with an average altitude of 1891 m. The overall terrain gradually decreases from north to south in a step-like manner. Kunming has a subtropical–plateau mountain monsoon climate and the annual average urban temperature is about 15 °C. The year-round temperature difference is small, but the daily temperature varies greatly, and the ultraviolet intensity is high, which means inversion layers form easily, hindering the diffusion of pollutants. The mountainous regions of the Yunnan–Guizhou Plateau are rich in minerals, and the densely distributed non-ferrous metallurgical enterprises around Kunming City are a source of atmospheric heavy metal pollution. In addition, Kunming is the most important transportation hub in Yunnan Province. As the economy opened up to Southeast Asia car ownership rose to 1.8 million in 2013, and in 2019 it was 58.29% higher [32,33]. Traffic density in Kunming in 2013 was 1570 vehicles/km2, comparable to that of Beijing (2030) [34], Chengdu (1320) [35], and Zhengzhou (2070) [36].
Because industries have moved to industrial parks away from urban areas, motor vehicle PM emissions seem to be the critical source of pollutants. Therefore, we conducted field measurements of the seasonal variations of PM2.5 heavy metals in Kunming to quantify the contribution of various sources, particularly motor vehicles, and assess the associated health risks. This study provides essential information for developing control strategies and protecting the health of people on the Yunnan–Guizhou Plateau.

2. Materials and Methods

2.1. Sampling Sites Description

The sampling of PM2.5 was performed at three sites in Kunming (Figure 1). Dongfeng East Road (DR, 25°24′ N, 102°42′ E; 1900 m AMSL), near a main road with a traffic volume of 3436 vehicles·h−1; Jinding Mountain (JDM, 25°0′ N, 102°41′ E; 1920 m AMSL), which is near steel and glass industries and where traffic volume is about 562 vehicles·h−1; and West Mountain (WM, 24°57′ N, 102°37′ E; 2300 m AMSL), a clean reference.

2.2. Sampling Schedule and Method

PM2.5 samples for 28 days were collected during 19–25 April and 24–30 October 2013, and 9–15 January and 21–28 May 2014. Daily samples were collected simultaneously at the sites for 23 ± 1 h from 9:00 a.m. to 9:00 a.m. During the sampling period, three standard medium-volume PM2.5 samplers (Wuhan Tianhong Intelligent Instrument, Model TH150) were employed at a flow rate of 100 L·min−1. Before sampling, we calibrated the flow rate with a soap film flowmeter (Qingdao Laoshan application technology research institute, Model 7030). Airborne particles were collected on 90 mm-diameter polypropylene-fiber filters, which had been annealed at 75 °C for 0.5 h to remove any impurities. The filters were stored in a desiccator at 25 °C and 45% relative humidity for 48 h before and after sampling and then weighted, packed in sample cells and stored at 15–30 °C till extraction and analysis. All filters were analyzed within 2 weeks. Table 1 lists the detailed sampling time and weather conditions at the three sites. Wind speed and temperature data are derived from https://www.wunderground.com/ (accessed on 23 August 2021).

2.3. Heavy Metal Analysis

Inductively Coupled Plasma Mass Spectrometry (ICP–MS, Agilent 7500a, Agilent, Palo Alto, CA, USA) analysis was employed to determine concentrations of 10 kinds of heavy metal elements: V, Cr, Mn, Co, Ni, Cu, Zn, As, Cd and Pb. Aluminum (Al) was analyzed to calculate the enrichment factors as the reference element. The scrap sampled filters were immersed in mixture of HCl and HNO3 (HCl: 16.7 mL, HNO3: 5.5 mL, ultrapure water: 150 mL) in a Teflon beaker. Then the beaker was covered by a lid, and the solution was placed on a 220 °C temperature-controlled electric hot plate for 2.5 h of circumfluence. As the liquid cooled, we washed out the inner wall of the beaker with ultrapure water, stewing for 0.5 h for extraction, filtration, then at constant volume to 50.0 mL under test.
Blanks (including filters) and duplicate sample analyses were performed for approximately 10% of all the samples. Certified reference materials (CRM) were used to ensure accuracy and precision (National Research Center of CRM, China). The field blank samples were analyzed using these same procedures. The results of the blank analyses were used to correct the corresponding samples.

2.4. Method

2.4.1. Enrichment Factor (EF)

EF is an important indicator of the extent of disturbance to the natural environment caused by human activity. By comparing measured values with soil background values, we can see the influence of human activity on particulates [37]. EF is calculated by the following equation:
EF = (Ci/Cn)sample/(Ci/Cn)soil background
where Ci represents the concentration of heavy metals; Cn represents the concentration of the reference element (we choose Al as the reference element in this study); (Ci/Cn)sample is the concentration ratio of the samples; (Ci/Cn)soil background is the concentration ratio of the corresponding element in the crust. If EF is close to 1, the element can be considered to originate mainly from soil particles; if EF > 10, the element mainly originates from human activity [38].

2.4.2. Principal Component Analysis (PCA)

PCA is a resource allocation method approved and recommended by the U.S. Environmental Protection Agency (EPA). It is an important multivariate statistical tool that can reduce the dimensionality of large datasets and extract the number of principal components needed to explain all the variances of such datasets, which is much less than the original number of variables [39,40]. This method establishes the orthogonal distribution between the components. No matter how many variables were included in this study [41], the regression adjustment results for each factor were simple and stable. During the analysis, the factors were determined by selecting principal components with eigenvalues greater than 1 according to previous studies [42,43]. They effectively explained and represented the source characteristics of each heavy metal in PM2.5, guided the understanding and control of heavy metal pollution emissions in the study area, and helped improve the health of local populations.

2.4.3. Health Risk Assessment (HRA)

This research used the health risk assessment model provided by the U.S. EPA as the basic framework, which was used for assessing the health risk of heavy metals in PM2.5 [44]. The target population of this study was children, and adult men and women. Then, based on the concentration levels of heavy metals in the study area, we calculated and assessed the health effects according to the health risk assessment model [45,46]. Because the heavy metals in PM2.5 enter the body through the respiratory system, this study focused on assessing the health risks of inhalation exposure, and the dose of this route was calculated by the following formula:
ADD(LADD)Inhale = (C × IRi × EF × ED × CF)/(BW × AT)
where C (ng·m−3) represents the concentration of heavy metal elements in PM2.5; the variable ADD (mg·(kg·d)−1) is normal average daily dose for non-carcinogenic elements; and LADD (mg·(kg·d)−1) is lifetime average daily dose for carcinogenic elements. The exposure parameters are mainly based on the exposure parameter manual of the population in China [47] as shown in Table 2.
According to the EPA’s Integrated Risk Information System (IRIS) and the International Agency for Research on Cancer (IARC), pollutants can be divided into carcinogens and non-carcinogens [48]. In this study, V, Mn, Cu, Zn and Pb are non-carcinogens; Cr, Co, Ni, As and Cd are carcinogens. The non-carcinogenic risk quotient (HQ) and carcinogenic risk (CR) of heavy metals were calculated as follows:
HQ = ADD/Rfd
CR = LADD × SF
In the formula, Rfd (mg·(kg·d)−1) is the reference dose of each heavy metal; SF ((kg·d)·mg−1) is the carcinogenic slope factor. The specific parameters were shown in Table 3 [49,50,51]. If HQ were less than 1, it meant that the non-carcinogenic risk was small or negligible; if it were greater than or equal to 1, it meant that there was a non-carcinogenic risk, which increased as the HQ value increased; CR represents the carcinogenic risk of heavy metals: when CR was between 10−6 and 10−4, it was within the acceptable range. If CR ≥ 10−4, there was a significant risk [52].

2.4.4. HYSPLIT4 Model

The HYSPLIT4 model is a professional model jointly developed by the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory (ARL) and the Australian Bureau of Meteorology (BOM) over the past 20 years to calculate and analyze the transport and diffusion of atmospheric pollutants. The model has a relatively complete transport, diffusion and sedimentation model that can handle multiple meteorological element input fields and physical processes, and different types of pollutant emission sources. It is widely used in the study of the transmission and diffusion of multiple pollutants in various regions. In this study, the independent version of the backward trajectory model and the auxiliary software package were used (GUI; Ghostscript; ImageMagick). Meteorological data were obtained through NCEP (National Centers for Environmental Prediction) and GDAS (Global Data Assimilation System).

3. Results and Discussion

3.1. Concentrations of Heavy Metals in PM2.5 in Kunming

3.1.1. Concentration Variation for Heavy Metals in PM2.5 on Working Days and Rest Days

During the four sampling campaigns, the average mass concentrations of PM2.5 at the three Kunming sites (DR, JDM, WM) were 125.16 ± 71.27 μg·m−3, 170.10 ± 104.10 μg·m−3 and 114.00 ± 73.99 μg·m−3, respectively. The overall concentration of the 10 heavy metals accounted for about 1.26% of PM2.5. Pb and Zn accounted for the highest proportions, 0.32 and 0.38%, respectively. Because most large industries in Kunming had moved to industrial parks far from urban areas, only a few industries were located around the JDM site. These industries were in a continuous operation, which shows that the concentrations of discharged heavy metals may not have fluctuated much between working days and rest days. Heavy-duty diesel vehicles are common around the site and emit a large amount of PM2.5 heavy metals on workdays. Figure 2a shows the concentration variation of 10 heavy metals on working days and rest days. The total concentrations ranged from 481.70 to 2052.65 ng·m−3 on working days, and 199.97 to 1105.21 ng·m−3 on rest days. The concentration levels of heavy metal markers (Cr, Mn, Cu, Zn, and Pb) from traffic sources are significantly higher on working days, especially in spring and winter. Figure 2b shows the concentration variation of heavy metals at the three sampling sites. It was found that almost all of the concentrations fluctuated wildly on working and rest days, especially at the DR and JDM sites, which suggest that heavy metals emitted by the traffic sources make an important contribution to ambient air.

3.1.2. Seasonal and Spatial Variations of Heavy Metals in Kunming

Seasonal and spatial variations in PM2.5 heavy metals in Kunming are shown in Figure 3, which also shows the concentration levels of all 10 heavy metals and their standard variance. The concentrations followed the order Zn > Pb > Mn > As > Cu > V > Cr > Ni > Cd > Co. The concentration of Cr, As and Cd exceeded the secondary level of National Ambient Air Quality Standard of China (GB 3095-2012, the limit concentrations of As, Cr and Cd are 6.00, 0.025 and 5.00 ng·m−3). The highest concentration was found in autumn at JDM with 996.6 ± 791.56 ng·m−3 (Zn), and the lowest concentration was found in autumn at WM with 0.52 ± 0.75 ng·m−3 (Co). Overall, the concentrations of Zn and Pb were much higher than other heavy metals, and the concentrations of most were consistent with the seasonal distribution of PM2.5, showing a pattern of high concentrations in winter and spring, and low concentrations in summer and autumn. Due to the intensive traffic flow and frequent industrial activity at the two sampling sites DR and JDM, concentrations of Ni, Cu, Zn, Cd and Pb were much higher than that at WM. Concentrations of Cr reached the highest value at DR because of heavy summer traffic (147.80 ± 247.85 ng·m−3), but summer concentrations of As and Pb were significantly higher at WM than at DR and JDM, possibly because WM is downwind from the petroleum, steel, phosphate, and salt industries at the Anning Industrial Park. However, there are many high mountains around the WM site that block the transmission of pollutants to Kunming. Therefore, the pollutants from Anning Industrial Park may not have had a significant impact on the DR and JDM sampling sites.

3.1.3. The Transport of Heavy Metals

To understand the transport of heavy metals in PM2.5 from distant sources, 72 h backward trajectories starting at 1000 m at the center of Kunming during the four sampling periods were calculated every 6 h using the Hybrid Single Particle Lagrangian Integrated Trajectory 4.0 (HYSPLIT4) model. Meteorological data came from the Global Data Assimilation System (GDAS). Figure 4 shows the clustering of air mass trajectories at 1000 m in Kunming. It was found that the air masses mainly originated in Myanmar in the southwest during the spring, summer and winter, and 47.02% of autumn air masses originated in Guizhou.
Figure 5 shows the concentration variations of total heavy metals and PM2.5 in Kunming (Figure 5a) and at WM (Figure 5b) during the four sampling periods. Combined with Figure 4, we can see that high concentrations occur in spring, possibly caused by the long-distance transmission of dense biomass burning from Myanmar. In summer, due to rainfall, concentrations were relatively low, and in autumn, clean air from Guizhou did not increase concentrations. In winter, however, concentrations of total heavy metals were the highest, which may be the result of coal burning by local and regional residents or of the long-distance transmission of pollution from India.
Figure 5b shows that the concentration variation of total heavy metals and PM2.5 at WM during all sampling periods followed the order spring > summer > winter > autumn. In summer, the high concentrations of PM2.5 may have been caused by forest fires in the region of Anning during early sampling, or the valley wind that brought pollution up to the WM site on the top of the mountain.

3.1.4. Comparisons with the Concentration of Heavy Metals in Other Cities

As shown in Table 4, the concentrations of most heavy metals in PM2.5 were generally higher than for most other cities, but they were lower than in Xi’an [53], which is an important industrial base and comprehensive transportation hub in Northwest China. The concentrations of some carcinogenic heavy metals (Cr, Ni and As) were comparable to those reported in other cities. Kunming’s As and Pb levels were 2.68 and 1.88 times higher, respectively, than those of Beijing [54]. Nanjing, Chengdu, Chongqing and Guangzhou [55,56,57] are all provincial capitals or municipalities with a dense transportation network, but the concentrations of heavy metals are much lower than for Kunming. Compared with Barcelona [58], all heavy metal concentrations for Kunming were higher: Cu, Zn and Pb were more than twice as high and Cr was 4.85 times higher. In Istanbul [59], the level of Cr was 4.18 times higher, but other heavy metals were lower. Similarly, as a plateau mountain city, the city of Lanzhou [60] also has very high concentrations of heavy metals in PM2.5. These results may explain to a certain extent the reason why the death rate due to problems with the respiratory system in Yunnan is higher than the national average. However, some cities have higher concentrations of heavy metals and lower death rates, which implies that there are other atmospheric pollutants in the environment that affect the respiratory system, such as persistent organic matter, etc.

3.2. Sources of Heavy Metals

3.2.1. Enrichment Factor (EF) and Correlation Analysis

In this study we used an enrichment factor (EF) and PCA to analyze possible sources of the 10 heavy metals in PM2.5 in Kunming. Figure 6 shows that the EFs of Cr, Mn, Ni, Cu, Zn, As, Cd and Pb were higher than 10, which suggested that these heavy metals were significantly affected by anthropogenic sources. The EFs of V and Co were lower than 10 but higher than 1, which suggests that the two elements may have been affected by natural and anthropogenic sources.
The correlation among the elements can indicate the sources of the heavy metals. In this paper, we used SPSS17 software to calculate the correlation coefficient between elements. The coefficient matrix is shown in Table 5.
Table 5 shows that Cr and Mn had the highest correlation coefficient of 0.915. Previous studies showed that Cr, Mn and Ni were the main pollutants in road dust [61,62]. Manoli et al. [63] also reported that road dust had high levels of Cr and Mn, indicating that this might be the source.
V, Ni, Cu and Pb were strongly correlated, and some studies showed that they are usually related to the metallurgical industries and vehicle exhaust emissions [64,65,66]. Cu and Pb may also come from brake and tire wear [67]. V may be influenced by the type of pavement or industrial dust [64,68]. Therefore, V, Ni, Cu and Pb may mainly come from vehicle emissions and metallurgical industries.
There was no significant correlation between As and the other heavy metals since it was commonly used as a coal combustion indicator [69]. The table shows that As in Kunming PM2.5 came mainly from coal combustion and did not highly correlate with vehicle emissions.

3.2.2. PCA of Heavy Metals

We use PCA to analyze heavy metal concentration data to identify possible sources of heavy metals in Kunming. The heavy metals data of all samples were input to the PCA model and analysis was performed. As shown in Table 6, three factors loading for heavy metals from PCA were listed, with varimax rotation of the data matrix. Factors 1–3 explained 49.43%, 18.73% and 12.61% of the total variance in the data set, respectively. In Factor 1, the heavy metals Ni, Cu, Zn, Cd and Pb have higher scores, which are 0.892, 0.862, 0.724, 0.912 and 0.880, respectively. The high load values of Ni, Cu, Zn, Cd and Pb indicate that Factor 1 is mainly from vehicle emissions and metallurgical industry sources [70,71]. The heavy metals Cr, Mn and Zn in Factor 2 are markers of traffic and road dust emission sources, so it is inferred that Factor 2 represents vehicle emissions and road dust [72]. In the Factor 3, heavy metal As has the highest score of 0.995. As is a marker component of coal, so it is inferred that Factor 3 is from coal combustion sources [48,73]. According to these results, we can know that vehicle emissions and metallurgical industry sources make the biggest important contribution to heavy metals in Kunming.

3.3. Health Risk Assessment of Heavy Metal Elements Exposure

In Kunming, the non-carcinogenic risk of heavy metals in PM2.5 at three sampling sites was analyzed through respiration. The results during the entire sampling period are shown in Table 7. Among the three sampling sites, JDM had the highest non-carcinogenic risk of heavy metals, followed by DR and WM, but the trends for all three were the same. The risk levels of non-carcinogenic heavy metals for children, and men and women were Mn > Pb > Cu > V > Zn; and the order is children > men > women. The non-carcinogenic risk for each heavy metal was less than 1, indicating that it was within the acceptable range of respiratory exposure. However, most were affected by Mn and Pb.
The carcinogenic risk of each heavy metal to each groups through respiration is shown in Table 8. At the three sampling sites, the carcinogenic risk of heavy metals for children was DR > WM > JDM; the carcinogenic risk of heavy metals for adults was DR > JDM > WM. The risk levels of carcinogenic heavy metals were Cr > As > Co > Cd > Ni. The carcinogenic risk of Cr to men and women in the DR sampling site exceeded 10−4. The carcinogenic risk of other heavy metals in all populations at the three sampling sites was between 10−6 and 10−4, which was within the acceptable range. Cr posed the highest carcinogenic risk to all people through respiration, and it should receive more attention.

3.4. Discussion

We conducted this study to find out whether heavy metals from PM2.5 had an impact on the death rate of the respiratory system in Yunnan. In the past few years, the number of deaths due to respiratory diseases in Yunnan Province is much higher than for the whole country, which implies that the environmental pollution has a significant impact on people in Yunnan. Heavy metals emitted from transportation sources may cause harm to the health of urban populations. Car ownership in Kunming has risen sharply in recent years, which is comparable to that of other medium and large cities in China. Therefore, we collected PM2.5 samples in Kunming to explore the contribution of heavy metals to environmental pollution and human health. The results are in line with the original assumptions. First, we analyzed the spatial and temporal distribution characteristics of heavy metal concentrations, and found that at sampling site where traffic source emissions are more obvious, the changes in heavy metal concentrations between working days and rest days fluctuate significantly. At the same time, we used PCA to analyze the source of heavy metals. The results show that Cr, Mn, Cu, Zn and Pb are all related to traffic source emissions. In the health risk assessment of heavy metal exposure, Cr has a significant carcinogenic risk for people in Kunming, which may explain to a certain extent the reason why the death rate of the respiratory system in Yunnan is higher than the national average.
There are some limitations to this study. Firstly, there is inadequate explanation regarding the death rate due to problems with the respiratory system in Yunnan, indicating that more variables need considering, for example, the level of persistent organics in the atmosphere. Secondly, the data is a little old, so health risks may not reflect current reality well. Thirdly, due to the complexity of the terrain of the plateau city, we only studied the characteristics of heavy metal pollution in Kunming, which is far from being able to clarify the entire mechanism of atmospheric heavy metal pollution in the plateau area.

4. Conclusions

The concentrations of PM2.5 at three sites (DR, JDM, WM) were 125.16 ± 71.27, 170.10 ± 104.10 and 114.00 ± 73.99 μg·m−3, respectively, and as for the concentrations of the 10 heavy metals (Zn > Pb > Mn > As > Cu > V > Cr > Ni > Cd > Co), Cr, As and Cd exceeded the second level of China’s national standard for ambient air quality (GB 3095-2012).
An EF analysis showed that Cr, Mn, Ni, Cu, Zn, As, Cd and Pb came mainly from anthropogenic sources, and that V and Co may have come from a mixture of natural and anthropogenic sources. For external sources, the clustering of air masses in Kunming showed that they mainly originated in Myanmar during the sampling periods, and the long-distance transmission of dense biomass burning in Myanmar in the spring affected the concentration levels of heavy metals.
We analyzed the concentration of heavy metals on working and rest days. Almost all concentrations fluctuated intensely on working and rest days at all three sampling sites, especially DR and JDM, which suggested that heavy metals from traffic made an important contribution to the ambient air. At the same time, the PCA results show that Ni, Cu, Zn, Cd and Pb came mainly from vehicle and metallurgical industry emissions (49.43% of the total variance), and Cr, Mn and Zn came from vehicle emissions and road dust (18.73% of the total variance). These results all showed that motor vehicles were already a critical air pollution source of heavy metals.
The total non-carcinogenic risk of heavy metals was less than 1, which was in the acceptable range of respiratory exposure. The risk levels of carcinogenic heavy metals were in the order: Cr > As > Co > Cd > Ni. The carcinogenic risk of Cr to men and women at the DR sampling site was greater than 10−4.
To reduce urban traffic emissions and protect people’s health, it is recommended that preferential policies be enacted to encourage consumers to use clean energy vehicles on the Yunnan–Guizhou Plateau.

Author Contributions

Conceptualization, X.H. and J.S.; investigation, X.P. and Y.B.; resources, J.S.; data curation, S.L.; writing—original draft preparation, Z.L.; writing—review and editing, J.S. and X.H.; supervision, P.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Projects of China (grant number 2019YFC0214405) and the National Natural Science Foundation of China (grant number 21966016 and 21667014).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this paper can be provided by Jianwu Shi ([email protected]).

Conflicts of Interest

The authors declare that there are no competing financial interests that could inappropriately influence the contents of this manuscript.

Nomenclature

PM2.5Particulate matter with aerodynamic equivalent diameter less than or equal to 2.5 microns in ambient air
VVanadium
CrChromium
MnManganese
CoCobalt
NiNickel
CuCopper
ZnZinc
AsArsenic
CdCadmium
PbLead
DRDongfeng East Road sampling site
JDMJinding mountain sampling site
WMWest mountain sampling site
SDStandard deviation
ICP-MSInductively coupled plasma mass spectrometry
AMSLAbove mean sea level
CRMCertified reference materials
EFEnrichment factor
PCAPrincipal component analysis
EPAU.S Environmental Protection Agency
HRAHealth risk assessment
IRISIntegrated Risk Information System
IARCInternational Agency for Research on Cancer
NOAANational Oceanic and Atmospheric Administration
ARLAir Resources Laboratory
BOMBureau of Meteorology
NCEPNational Centers for Environmental Prediction
GDASGlobal Data Assimilation System

References

  1. Zong, Z.; Wang, X.P.; Tian, C.G.; Chen, Y.J.; Qu, L.; Ji, L.; Zhi, G.R.; Li, J.; Zhang, G. Source apportionment of PM2.5 at a regional background site in North China using PMF linked with radiocarbon analysis: Insight into the contribution of biomass burning. Atmos. Chem. Phys. 2016, 16, 11249–11265. [Google Scholar] [CrossRef] [Green Version]
  2. Li, L.; Wang, W.; Feng, J.L.; Zhang, D.P.; Li, H.J.; Gu, Z.P.; Wang, B.J.; Sheng, G.Y.; Fu, J.M. Composition, source, mass closure of PM2.5 aerosols for four forests in eastern China. J. Environ. Sci. 2010, 22, 405–412. [Google Scholar] [CrossRef]
  3. Lu, D.B.; Xu, J.H.; Yang, D.Y.; Zhao, J.N. Spatio-temporal variation and influence factors of PM2.5 concentrations in China from 1998 to 2014. Atmos. Pollut. Res. 2017, 8, 1151–1159. [Google Scholar] [CrossRef]
  4. Coronas, M.V.; Rocha, J.A.V.; Salvadori, D.M.F.; Vargas, V.M.F. Evaluation of area contaminated by wood treatment activities: Genetic markers in the environment and in the child population. Chemosphere 2016, 144, 1207–1215. [Google Scholar] [CrossRef]
  5. Shi, L.N.; Xu, X.; Dou, X.Y.; Zhao, X.D. The research progress of heavy metals in PM2.5. Adv. Mater. Res. 2014, 955, 1397–1404. [Google Scholar] [CrossRef]
  6. Zhang, N.; Han, B.; He, F.; Xu, J.; Niu, C.; Zhuo, J.; Kong, S.F.; Bai, Z.P.; Xu, H. Characterization, health risk of heavy metals, and source apportionment of atmospheric PM2.5 to children in summer and winter: An exposure panel study in Tianjin, China. Air Qual. Atmos. Health 2015, 8, 347–357. [Google Scholar] [CrossRef]
  7. Lai, C.H.; Lin, C.H.; Liao, C.C. Respiratory deposition and health risk of inhalation of particle-bound heavy metals in the carbon black feeding area of a tire manufacturer. Air Qual. Atmos. Health 2017, 10, 1281–1289. [Google Scholar] [CrossRef]
  8. Heys, K.A.; Shore, R.F.; Pereira, M.G.; Jones, K.C.; Martin, F.L. Risk assessment of environmental mixture effects. RSC Adv. 2016, 6, 47844–47857. [Google Scholar] [CrossRef] [Green Version]
  9. Rodriguez-Espinosa, P.F.; Flores-Rangel, R.M.; Mugica-Alvarez, V.; Morales-Garcia, S.S. Sources of trace metals in PM10 from a petrochemical industrial complex in Northern Mexico. Air Qual. Atmos. Health 2017, 10, 69–84. [Google Scholar] [CrossRef]
  10. Pope, C.A.; Burnett, R.T.; Thun, M.J.; Calle, E.E.; Krewski, D.; Ito, K.; Thurston, G.D. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. Jama 2002, 287, 1132–1141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Dockery, D.W.; Pope, C.A. Acute respiratory effects of particulate air pollution. Ann. Rev. Public Health 1994, 15, 107–132. [Google Scholar] [CrossRef]
  12. Pope, C.A.; Dockery, D.W. Health effects of fine particulate air pollution: Lines that connect. J. Air Waste Manag. Assoc. 2006, 56, 709–742. [Google Scholar] [CrossRef]
  13. Ahmad, H.R.; Mehmood, K.; Sardar, M.F.; Maqsood, M.A.; Rehman, M.Z.U.; Zhu, C.X.; Li, H.N. Integrated risk assessment of potentially toxic elements and particle pollution in urban road dust of megacity of Pakistan. Human Ecol. Risk Assess. An. Int. J. 2020, 26, 1810–1831. [Google Scholar] [CrossRef]
  14. Khreis, H.; Warsow, K.M.; Verlinghieri, E.; Guzman, A.; Pellecuer, L.; Jones, A.F.I.; Heinen, E.; Rojas-Rueda, D.; Mueller, N.; Schepers, P.; et al. The health impacts of traffic-related exposures in urban areas: Understanding real effects, underlying driving forces and co-producing future directions. J. Transport. Health 2016, 3, 249–267. [Google Scholar] [CrossRef]
  15. Zamora, M.L.; Pulczinski, J.C.; Johnson, N.; Garcia-Hernandez, R.; Rule, A.; Carrillo, G.; Zietsman, J.; Sandragorsian, B.; Vallamsundar, S.; Askariyeh, M.H.; et al. Maternal exposure to PM2.5 in south Texas, a pilot study. Sci. Total Environ. 2018, 628–629, 1497–1507. [Google Scholar] [CrossRef] [PubMed]
  16. Sohrabi, S.; Zietsman, J.; Khreis, H. Burden of disease assessment of ambient air pollution and premature mortality in urban areas: The role of socioeconomic status and transportation. Int. J. Environ. Res. Public Health 2020, 17, 1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Fang, L.P.; Xu, C.H.; Li, J.; Borggaard, O.K.; Wang, D.S. The importance of environmental factors and matrices in the adsorption, desorption, and toxicity of butyltins: A review. Environ. Sci. Pollut. Res. 2017, 24, 9159–9173. [Google Scholar] [CrossRef] [PubMed]
  18. He, K.B.; Huo, H.; Wang, Q.D.; Yao, Z.L. Method and Application of Road Vehicle Emission Model.; Science Press: Beijing, China, 2014; pp. 1–377. (In Chinese) [Google Scholar]
  19. Lu, X.W.; Wang, L.J.; Li, L.Y.; Lei, K.; Huang, L.; Kang, D. Multivariate statistical analysis of heavy metals in street dust of Baoji, NW China. J. Hazard. Mater. 2010, 173, 744–749. [Google Scholar] [CrossRef]
  20. Apeagyei, E.; Bank, M.S.; Spengler, J.D. Distribution of heavy metals in road dust along an urban-rural gradient in Massachusetts. Atmos. Environ. 2011, 45, 2310–2323. [Google Scholar] [CrossRef]
  21. Cobbina, S.J.; Chen, Y.; Zhou, Z.X.; Wu, X.S.; Feng, W.W.; Wang, W.; Mao, G.G.; Xu, H.; Zhang, Z.; Wu, X.Y.; et al. Low concentration toxic metal mixture interactions: Effects on essential and non-essential metals in brain, liver, and kidneys of mice on sub-chronic exposure. Chemosphere 2015, 132, 79–86. [Google Scholar] [CrossRef]
  22. Wang, C.M.; Ouyang, H.; Wang, J.D.; Liu, J.S.; Zhang, X.L.; Wang, Y. Impact of lead pollution in environment on children’s health in Shenyang city. Environ. Sci. 2003, 24, 17–22. (In Chinese) [Google Scholar]
  23. Du, Y.R.; Gao, B.; Zhou, H.D.; Ju, X.X.; Hao, H.; Yin, S.H. Health risk assessment of heavy metals in road dusts in urban parks of Beijing, China. Procedia Environ. Sci. 2013, 18, 299–309. [Google Scholar] [CrossRef] [Green Version]
  24. Qi, J.Y.; Zhang, H.L.; Fang, J.D.; Ren, M.Z.; Lin, B.K.; Wang, X.; Zhang, F.F. Characterization of airborne particulate metals in the surroundings of a municipal solid waste incinerator (MSWI) in the Nanshan Shenzhen and health risk assessment via inhalation exposure. China Environ. Sci. 2013, 33, 113–118. (In Chinese) [Google Scholar]
  25. Hsieh, C.Y.; Jung, C.R.; Lin, C.Y.; Hwang, B.F. Combined exposure to heavy metals in PM2.5 and pediatric asthma. J. Allergy Clin. Immunol. 2020, 147, 2171–2180. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, Y.F.; Shao, Y.; Chen, Y.; Xiao, Y.Z.; Qing, M.F. Analysis on causes of death among residents in Yunnan surveillance-spots in 2014. Chin. J. Health Educ. 2016, 32, 490–494. (In Chinese) [Google Scholar]
  27. Yang, Y.F.; Shao, Y.; Chen, Y.; Qing, M.F.; Xiao, Y.Z. Monitoring results of causes of death among residents in disease surveillance spots of Yunnan Province in 2015. Occup. Health 2017, 33, 1793–1796. (In Chinese) [Google Scholar]
  28. Yang, Y.F.; Qing, M.F.; Chen, Y.; Shi, Q.P.; Tang, X.; Yu, Q.L.; Shao, Y. Analysis on epidemiological characteristics of deaths caused by respiratory disease among residents, Yunnan province, 2016–2019. Prev. Med. Trib. 2021, 27, 276–280. [Google Scholar]
  29. SCP. China Cause of Death Surveillance Data Set; Science and Technology of China Press: Beijing, China, 2019. (In Chinese) [Google Scholar]
  30. Wang, H.L.; Qiao, L.P.; Lou, S.R.; Zhou, M.; Ding, A.J.; Huang, H.Y.; Chen, J.M.; Wang, Q.; Tao, S.K.; Chen, C.H.; et al. Chemical composition of PM2.5 and meteorological impact among three years in urban Shanghai, China. J. Clean. Prod. 2016, 112, 1302–1311. [Google Scholar] [CrossRef]
  31. Zhu, W.H.; Xu, X.D.; Zheng, J.; Yan, P.; Wang, Y.J.; Cai, W.Y. The characteristics of abnormal wintertime pollution events in the Jing-Jin-Ji region and its relationships with meteorological factors. Sci. Total Environ. 2018, 626, 887–898. [Google Scholar] [CrossRef]
  32. CSP. Statistical Bureau of Yunnan Province. Yunnan Statistical Yearbook; China Statistics Press: Beijing, China, 2020. (In Chinese) [Google Scholar]
  33. CSP. Statistical Bureau of Yunnan Province. Yunnan Statistical Yearbook; China Statistics Press: Beijing, China, 2014. (In Chinese) [Google Scholar]
  34. CSP. Beijing Statistical Yearbook; China Statistics Press: Beijing, China, 2014. (In Chinese) [Google Scholar]
  35. XPH. Chengdu Yearbook; Xinhua Publishing House: Beijing, China, 2014. (In Chinese) [Google Scholar]
  36. CSP. Henan Statistical Yearbook; China Statistics Press: Beijing, China, 2014. (In Chinese) [Google Scholar]
  37. Loska, K.; Wiechula, D.; Korus, I. Metal contamination of farming soils affected by industry. Environ. Int. 2004, 30, 159–165. [Google Scholar] [CrossRef]
  38. Ji, Y.; Zhu, T.; Feng, Y.C.; Bai, Z.P.; Liu, H.X. Evaluate the elements pollution of urban soil wind dust using enrichment factor method. Acta Sci. Nat. Univ. Nankaiensis 2006, 39, 94–99. (In Chinese) [Google Scholar]
  39. Jain, S.; Sharma, S.; Mandal, T.; Saxena, M. Source apportionment of PM10in Delhi, India using PCA/APCS, UNMIX and PMF. Particuology 2018, 37, 107–118. [Google Scholar] [CrossRef]
  40. Dallarosa, J.; Teixeira, E.C.; Meira, L.; Wiegand, F. Study of the chemical elements and polycyclic aromatic hydrocarbons in atmospheric particles of PM10and PM2.5in the urban and rural areas of South Brazil. Atmos. Res. 2008, 89, 76–92. [Google Scholar] [CrossRef]
  41. Gupta, S.; Gadi, R.; Sharma, S.; Mandal, T. Characterization and source apportionment of organic compounds in PM10using PCA and PMF at a traffic hotspot of Delhi. Sustain. Cities Soc. 2018, 39, 52–67. [Google Scholar] [CrossRef]
  42. Henry, R.C.; Lewis, C.W.; Hopke, P.K.; Williamson, H.J. Review of receptor model fundamentals. Atmos. Environ. 1984, 18, 1507–1515. [Google Scholar] [CrossRef]
  43. Henry, R.C. Current factor analysis receptor models are ill-posed. Atmos. Environ. 1987, 21, 1815–1820. [Google Scholar] [CrossRef]
  44. EPA. Risk Assessment Guidance for Superfund Volume I: Human Health Evaluation Manual; EPA: Washington, DC, USA, 1989. [Google Scholar]
  45. Holgate, S.T.; Smaet, Z.M.; Koren, H.; Maynard, R.L. Air Pollution and Health; Academic Press: London, UK, 1999. [Google Scholar]
  46. Yu, Y.J.; Zhang, Y.; Che, F.; Sun, P.; Cai, Y.; Xiang, M.D.; Wang, Q. Health risk assessment of environmental pollution and the application. J. Environ. Occup. Med. 2011, 28, 309–313. (In Chinese) [Google Scholar]
  47. MEP. Chinese Population Exposure Parameter Manual; Chinese Environmental Science Press: Beijing, China, 2013. [Google Scholar]
  48. Tian, H.Z.; Wang, Y.; Xue, Z.G.; Cheng, K.; Qu, Y.P.; Chai, F.H.; Hao, J.M. Trend and characteristics of atmospheric emissions of Hg, As and Se from coal combustion in China, 1980–2007. Atmos. Chem. Phys. 2010, 10, 11905–11919. [Google Scholar] [CrossRef] [Green Version]
  49. Gao, Y.; Guo, X.Y.; Li, C.; Ding, H.J.; Tang, L.; Ji, H.B. Characteristics of PM2.5 in Miyun, the northeastern suburb of Beijing: Chemical composition and evaluation of health risk. Environ. Sci. Pollut. Res. 2015, 22, 16688–16699. [Google Scholar] [CrossRef] [PubMed]
  50. Betha, R.; Pradani, M.; Lestari, P.; Joshi, U.M.; Reid, J.S.; Balasubramanian, R. Chemical speciation of trace metals emitted from Indonesian peat fires for health risk assessment. Atmos. Res. 2013, 122, 571–578. [Google Scholar] [CrossRef]
  51. Chen, P.F.; Bi, X.H.; Zhang, J.Q.; Wu, J.H.; Feng, Y.C. Assessment of heavy metal pollution characteristics and human health risk of exposure to ambient PM2.5 in Tianjin, China. Particuology 2015, 20, 104–109. [Google Scholar] [CrossRef]
  52. USEPA. Child Specific Exposure Factors Handbook; National Center for Environmental Assessment: Washington, DC, USA, 2002. [Google Scholar]
  53. Liu, P.P.; Lei, Y.; Ren, H.R.; Gao, J.J.; Xu, H.M.; Shen, Z.X.; Zhang, Q.; Zheng, C.L.; Liu, H.X.; Zhang, R.J.; et al. Seasonal variation and health risk assessment of heavy metals in PM2.5 during winter and summer over Xi’an, China. Atmosphere 2017, 8, 91. [Google Scholar] [CrossRef] [Green Version]
  54. Gao, J.J.; Wang, K.; Wang, Y.; Liu, S.H.; Zhu, C.Y.; Hao, J.M.; Liu, H.J.; Hua, S.B.; Tian, H.Z. 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]
  55. Wang, H.B.; Qiao, B.Q.; Zhang, L.M.; Yang, F.M.; Jiang, X. Characteristics and sources of trace elements in PM2.5 in two megacities in Sichuan Basin of southwest China. Environ. Pollut. 2018, 242, 1577–1586. [Google Scholar] [CrossRef]
  56. Ming, L.L.; Jin, L.; Li, J.; Fu, P.Q.; Yang, W.Y.; Liu, D.; Zhang, G.; Wang, Z.F.; Li, X.D. PM2.5 in the Yangtze River delta, China: Chemical compositions, seasonal variations, and regional pollution events. Environ. Pollut. 2017, 223, 200–212. [Google Scholar] [CrossRef]
  57. Tao, J.; Zhang, L.M.; Cao, J.J.; Zhong, L.J.; Chen, D.S.; Yang, Y.H.; Chen, D.H.; Chen, L.G.; Zhang, Z.S.; Wu, Y.F.; et al. Source apportionment of PM2.5 at urban and suburban areas of the Pearl River Delta region, south China—With emphasis on ship emissions. Sci. Total Environ. 2017, 574, 1559–1570. [Google Scholar] [CrossRef]
  58. Moreno, T.; Querol, X.; Alastuey, A.; Viana, M.; Salvador, P.; de la Campa, A.S.; Artiñano, B.; de la Rosa, J.; Gibbons, W. Variations in atmospheric PM trace metal content in Spanish towns: Illustrating the chemical complexity of the inorganic urban aerosol cocktail. Atmos. Environ. 2006, 40, 6791–6803. [Google Scholar] [CrossRef]
  59. Onat, B.; Sahin, U.A.; Akyuz, T. Elemental characterization of PM2.5 and PM1 in dense traffic area in Istanbul, Turkey. Atmos. Pollut. Res. 2013, 4, 101–105. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, Y.N.; Jia, C.H.; Tao, J.; Zhang, L.M.; Liang, X.X.; Ma, J.M.; Gao, H.; Huang, T.; Zhang, K. Chemical characterization and source apportionment of PM2.5 in a semi-arid and petrochemical-industrialized city, Northwest China. Sci. Total Environ. 2016, 573, 1031–1040. [Google Scholar] [CrossRef] [PubMed]
  61. Lee, B.K.; Hieu, N.T. Seasonal variation and sources of heavy metals in atmospheric aerosols in a residential area of Ulsan, Korea. Aerosol Air Qual. Res. 2011, 11, 679–688. [Google Scholar] [CrossRef]
  62. Chang, S.H.; Wang, K.S.; Chang, H.F.; Ni, W.W.; Wu, B.J.; Wong, R.H.; Lee, H.S. Comparison of source identification of metals in road-dust and soil. Soil Sediment. Contam. 2009, 18, 669–683. [Google Scholar] [CrossRef]
  63. Manoli, E.; Voutsa, D.; Samara, C. Chemical characterization and source identification/apportionment of fine and coarse air particles in Thessaloniki, Greece. Atmos. Environ. 2002, 36, 949–961. [Google Scholar] [CrossRef]
  64. Cheng, X.; Huang, Y.; Long, Z.J.; Ni, S.J.; Shi, Z.M.; Zhang, C.J. Characteristics, sources and health risk assessment of trace metals in PM10 in Panzhihua, China. Bullet. Environ. Contam. Toxicol. 2017, 98, 76–83. [Google Scholar] [CrossRef] [PubMed]
  65. Xia, L.; Gao, Y. Characterization of trace elements in PM2.5 aerosols in the vicinity of highways in northeast New Jersey in the US east coast. Atmos. Pollut. Res. 2011, 2, 34–44. [Google Scholar] [CrossRef] [Green Version]
  66. Xu, H.M.; Cao, J.J.; Ho, K.F.; Ding, H.; Han, Y.M.; Wang, G.H.; Chow, J.C.; Watson, J.G.; Khol, S.D.; Qiang, J.; et al. Lead concentrations in fine particulate matter after the phasing out of leaded gasoline in Xi’an, China. Atmos. Environ. 2012, 46, 217–224. [Google Scholar] [CrossRef]
  67. Li, Y.P.; Zhang, Z.S.; Liu, H.F.; Zhou, H.F.; Fan, Z.G.; Lin, M.; Wu, D.L.; Xia, B.C. Characteristics, sources and health risk assessment of toxic heavy metals in PM2.5 at a megacity of southwest China. Environ. Geochem. Health 2016, 38, 353–362. [Google Scholar] [CrossRef]
  68. Fiala, M.; Hwang, H.M. Influence of highway pavement on metals in road dust: A case study in Houston, Texas. Water Air Soil Pollut. 2021, 232, 1–12. [Google Scholar] [CrossRef]
  69. Tian, H.Z.; Wang, Y.; Xue, Z.G.; Qu, Y.P.; Chai, F.H.; Hao, J.M. Atmospheric emissions estimation of Hg, As, and Se from coal-fired power plants in China, 2007. Sci. Total Environ. 2011, 409, 3078–3081. [Google Scholar] [CrossRef]
  70. Taiwo, A.M.; Harrison, R.M.; Shi, Z. A review of receptor modeling of industrially emitted particulate matter. Atmos. Environ. 2014, 97, 109–120. [Google Scholar] [CrossRef] [Green Version]
  71. Guo, G.H.; Lei, M.; Chen, T.B.; Song, B.; Li, X.Y. Effect of road traffic on heavy metals in road dusts and roadside soils. Acta Sci. Circumstantiate 2008, 28, 1937–1945. (In Chinese) [Google Scholar]
  72. Yatkin, S.; Bayram, A. Elemental composition and sources of particulate matter in the ambient air of a Metropolitan City. Atmos. Res. 2007, 85, 126–139. [Google Scholar] [CrossRef]
  73. Miguel, E.D.; Iribarren, I.; Chacon, E.; Ordonez, A.; Charlesworth, S. Risk-based evaluation of the exposure of children to trace elements in playgrounds in Madrid (Spain). Chemosphere 2007, 66, 505–513. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location map of the sampling sites in Kunming (WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; km: kilometer; N in the upper left corner of the picture: North arrow; N on the ordinate in the figure: North latitude; E on the abscissa in the figure: East longitude).
Figure 1. Location map of the sampling sites in Kunming (WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; km: kilometer; N in the upper left corner of the picture: North arrow; N on the ordinate in the figure: North latitude; E on the abscissa in the figure: East longitude).
Ijerph 18 11004 g001
Figure 2. (a) Variation in heavy metal concentrations on working and rest days with the seasons; (b) Variation of heavy metals concentration on working and rest days at different sampling sites. (V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; ng·m−3: Concentration unit, indicating the number of nanograms per cubic meter; SD: Standard deviation; n: Number of samples).
Figure 2. (a) Variation in heavy metal concentrations on working and rest days with the seasons; (b) Variation of heavy metals concentration on working and rest days at different sampling sites. (V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; ng·m−3: Concentration unit, indicating the number of nanograms per cubic meter; SD: Standard deviation; n: Number of samples).
Ijerph 18 11004 g002
Figure 3. Variation of heavy metals concentrations in three sites of Kunming (V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; ng·m−3: Concentration unit, indicating the number of nanograms per cubic meter; SD: Standard deviation; n: Number of samples).
Figure 3. Variation of heavy metals concentrations in three sites of Kunming (V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; ng·m−3: Concentration unit, indicating the number of nanograms per cubic meter; SD: Standard deviation; n: Number of samples).
Ijerph 18 11004 g003
Figure 4. The clustering of air mass in Kunming during the four sampling periods (hPa: Hectopascal; Age Hour: Indicates 72 h before the selected time point of the backward trajectory; N on the upper left corner of the figure: North arrow; N on the left axis: North latitude; E on the upper axis: East longitude; 1 and 2 represent the clustering lines of the air mass trajectories).
Figure 4. The clustering of air mass in Kunming during the four sampling periods (hPa: Hectopascal; Age Hour: Indicates 72 h before the selected time point of the backward trajectory; N on the upper left corner of the figure: North arrow; N on the left axis: North latitude; E on the upper axis: East longitude; 1 and 2 represent the clustering lines of the air mass trajectories).
Ijerph 18 11004 g004
Figure 5. Concentration variation of total heavy metals and PM2.5 (a) in Kunming and (b) at the WM site during all sampling periods (WM: West Mountain, Clean reference sampling site; μg·m3: Concentration unit, indicating the number of micrograms per cubic meter; SD: Standard deviation; n: Number of samples).
Figure 5. Concentration variation of total heavy metals and PM2.5 (a) in Kunming and (b) at the WM site during all sampling periods (WM: West Mountain, Clean reference sampling site; μg·m3: Concentration unit, indicating the number of micrograms per cubic meter; SD: Standard deviation; n: Number of samples).
Ijerph 18 11004 g005
Figure 6. Enrichment factors of heavy metals in PM2.5 in Kunming (EFs: Enrichment factors; SD: Standard deviation; n: Number of samples; V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead).
Figure 6. Enrichment factors of heavy metals in PM2.5 in Kunming (EFs: Enrichment factors; SD: Standard deviation; n: Number of samples; V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead).
Ijerph 18 11004 g006
Table 1. Sampling time and weather conditions at three sampling sites in Kunming.
Table 1. Sampling time and weather conditions at three sampling sites in Kunming.
SeasonSampling TimeSampling SitesWeatherWind Speed (m·s−1)Concentration of PM2.5 (μg·m3, Mean ± SD, n = 3)Air Temperature (°C)
WMJDMDR
Spring4/18/2013-4/19/2013, (23 ± 1 h)YYYCloudy7.25162.32 ± 46.3618.00
4/19/2013-4/20/2013, (23 ± 1 h)YYYCloudy7.77267.42 ± 127.1518.00
4/20/2013-4/21/2013, (23 ± 1 h)YYYClear4.78242.5 6± 64.6720.50
4/21/2013-4/22/2013, (23 ± 1 h)YYYClear6.67244.74 ± 41.3520.50
4/22/2013-4/23/2013, (23 ± 1 h)YYYClear8.01293.83 ± 32.7820.50
4/23/2013-4/24/2013, (23 ± 1 h)NNYClear6.71232.4421.00
4/24/2013-4/25/2013, (23 ± 1 h)NNYCloudy6.82187.8021.00
Summer5/21/2014-5/22/2014, (23 ± 1 h)YYYShower8.72101.02 ± 4.0022.50
5/22/2014-5/23/2014, (23 ± 1 h)YYYShower7.38117.67 ± 46.9021.50
5/23/2014-5/24/2014, (23 ± 1 h)YYYShower7.29158.44 ± 84.0022.00
5/24/2014-5/25/2014, (23 ± 1 h)YYYCloudy6.80109.13 ± 51.2821.50
5/25/2014-5/26/2014, (23 ± 1 h)YYYCloudy7.0679.69 ± 8.1322.00
5/26/2014-5/27/2014, (23 ± 1 h)NNYShower2.9553.6222.00
5/27/2014-5/28/2014, (23 ± 1 h)NNYShower4.2550.6521.50
Autumn10/23/2013-10/24/2013, (23 ± 1 h)YYYModerate rain5.3647.49 ± 17.4113.50
10/24/2013-10/25/2013, (23 ± 1 h)YYYLight rain2.9185.37 ± 35.0113.50
10/25/2013-10/26/2013, (23 ± 1 h)YYYCloudy2.28122.62 ± 45.6213.00
10/26/2013-10/27/2013, (23 ± 1 h)YYYCloudy2.9589.31 ± 21.7912.50
10/27/2013-10/28/2013, (23 ± 1 h)YYYShower3.6751.79 ± 15.9512.00
10/28/2013-10/29/2013, (23 ± 1 h)NNYShower4.9226.6713.00
10/29/2013-10/30/2013, (23 ± 1 h)NNYShower4.7854.0814.00
Winter1/9/2014-1/10/2014, (23 ± 1 h)YYYClear5.50101.42 ± 46.568.50
1/10/2014-1/11/2014, (23 ± 1 h)YYYClear8.09131.08 ± 74.1311.50
1/11/2014-1/12/2014, (23 ± 1 h)YYYClear6.88127.56 ± 69.1210.50
1/12/2014-1/13/2014, (23 ± 1 h)YYYClear4.8367.16 ± 1.729.50
1/13/2014-1/14/2014, (23 ± 1 h)YYYSleet4.25152.41 ± 40.894.00
1/14/2014-1/15/2014, (23 ± 1 h)NNYCloudy4.47132.297.00
1/15/2014-1/16/2014, (23 ± 1 h)NNYCloudy6.93113.707.50
WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; m·s−1: Speed unit, expressed in meters per second; μg·m−3: Concentration unit, indicating the number of micrograms per cubic meter; SD: Standard deviation; n: Number of samples; °C: Degree Celsius; h: Hour; Y: Sampling; N: No sampling.
Table 2. The exposure parameter values.
Table 2. The exposure parameter values.
ParameterPhysical SignificanceValueUnit
ChildrenMenWomen
EFExposure relative frequency365365365d·a−1
EDExposure duration63030a
CFConversion coefficient10−610−610−6kg·mg−1
IRiInhalation rate515.211.3m3·d−1
BWAverage body rate1562.754.4kg
ATAveraging time (non-carcinogens)219010,95010,950d
Averaging time (carcinogens)25,50025,50025,500d
d: Day; a: Year: kg: Kilogram; mg: Milligram; m: Meter.
Table 3. Response parameters of elements entering through the respiratory system [49,50,51].
Table 3. Response parameters of elements entering through the respiratory system [49,50,51].
ElementRiskRfd (mg·(kg·d)−1)SF ((kg·d)·mg−1)
VNon-carcinogenic7 × 10−3-
MnNon-carcinogenic3 × 10−4-
CuNon-carcinogenic1.43 × 10−2-
ZnNon-carcinogenic0.3-
PbNon-carcinogenic3.5 × 10−3-
CrCarcinogenic-42
CoCarcinogenic-32
NiCarcinogenic-0.84
AsCarcinogenic-20.07
CdCarcinogenic-8.4
V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; Rfd: The reference dose of each heavy metal; SF: The carcinogenic slope factor; kg: Kilogram; mg: Milligram; d: Day; “-” represents no data.
Table 4. Concentrations of heavy metals in PM2.5 in Kunming compared with other cities (Unit: ng·m−3).
Table 4. Concentrations of heavy metals in PM2.5 in Kunming compared with other cities (Unit: ng·m−3).
CitiesAltitude (m)Sampling SitesYearSeasonPM2.5 (μg·m−3)VCrMnCoNiCuZnAsCdPbRef.
Kunming (Mean±SD, n = 68)1891urban2013~2014Four seasons130.45 ± 86.3540.30 ± 92.0829.10 ± 32.21155.79 ± 68.712.20 ± 2.0620.50 ± 14.7678.80 ± 127.37327.00 ± 233.7826.76 ± 21.926.00 ± 4.42281.50 ± 225.87This study
Beijing43.5urban2014Four seasons126.00-30.0070.00-40.00200.00310.0010.00-150.00[54]
Nanjing8.9urban2013~2014Summer97.80 ± 40.509.8813.20--9.3024.70247.00--90.90[55]
Chengdu505.9urban2014~2015Four seasons-1.905.6033.80-2.1018.70238.0010.80-55.40[56]
Chongqing259.1urban2014~2015Four seasons--11.1037.70-4.2011.30113.00--50.40[56]
Guangzhou6.6urban2014Four seasons48.00 ± 22.009.009.0034.00-4.0037.00225.00--77.00[57]
Lanzhou1520urban2014Four seasons88.90 ± 52.00--120.00-33.80106.00237.0010.50-491.00[60]
Xi’an396.9urban2008~2009Summer, winter--152.0095.00-19.0041.001775.00127.008.00408.00[53]
Barcelona20urban2001Four seasons--6.0014.00--32.00160.00--120.00[58]
Istanbul7urban2010Summer40.50 ± 13.70 2.54121.7042.100.52-19.60384.70---[59]
SD: Standard deviation; n: Number of samples; μg·m−3: Concentration unit, indicating the number of micrograms per cubic meter; ng·m−3: Concentration unit, indicating the number of nanograms per cubic meter; V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; Ref.: References; “-” represents no data.
Table 5. Correlation coefficient of heavy metals in PM2.5 of Kunming (n = 68).
Table 5. Correlation coefficient of heavy metals in PM2.5 of Kunming (n = 68).
VCrMnCoNiCuZnAsCdPb
V1
Cr0.2671
Mn0.3050.915 **1
Co−0.0580.4520.1981
Ni0.881 **0.2310.309−0.2161
Cu0.821 **0.5610.598 *−0.1520.870 **1
Zn0.042−0.0190.113−0.1830.2330.2181
As−0.090−0.144−0.065−0.191−0.008−0.118−0.0121
Cd0.366−0.194−0.162−0.2060.596 *0.3160.2950.2131
Pb0.791 **0.5650.5650.2130.777 **0.828 **0.1530.1510.3941
n: Number of samples; V: Vanadium; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; * Showed significant correlation at 0.05 level (double side); ** Showed significant correlation at 0.01 level (double side).
Table 6. PCA analysis of heavy metal elements in PM2.5 of Kunming (n = 68).
Table 6. PCA analysis of heavy metal elements in PM2.5 of Kunming (n = 68).
CrMnNiCuZnAsCdPbVariance (%)Cumulative (%)
Factor 10.1100.0180.8920.8620.7240.2250.9120.88049.4349.43
Factor 20.7350.864−0.1290.3130.529−0.032−0.0890.25318.7368.16
Factor 3−0.019−0.019−0.097−0.038−0.0070.9950.1570.04712.6180.77
n: Number of samples; Cr: Chromium; Mn: Manganese; Co: Cobalt; Ni: Nickel; Cu: Copper; Zn: Zinc; As: Arsenic; Cd: Cadmium; Pb: Lead; %: Percent sign.
Table 7. Non-carcinogenic risk of heavy metals in PM2.5 in different groups in Kunming.
Table 7. Non-carcinogenic risk of heavy metals in PM2.5 in different groups in Kunming.
ElementsDRJDMWM
ChildrenMenWomenChildrenMenWomenChildrenMenWomen
V5.46 × 10−43.97 × 10−43.4 × 10−49.77 × 10−47.11 × 10−46.09 × 10−43.89 × 10−42.83 × 10−42.43 × 10−4
Mn1.72 × 10−11.25 × 10−11.07 × 10−12.11 × 10−11.53 × 10−11.31 × 10−11.37 × 10−19.93 × 10−28.51 × 10−2
Cu8.2 × 10−45.96 × 10−45.11 × 10−44.78 × 10−33.47 × 10−32.98 × 10−35.58 × 10−44.06 × 10−43.48 × 10−4
Zn3.12 × 10−42.27 × 10−41.94 × 10−46.38 × 10−44.64 × 10−43.97 × 10−41.79 × 10−41.3 × 10−41.11 × 10−4
Pb2.06 × 10−21.5 × 10−21.29 × 10−24 × 10−22.91 × 10−22.49 × 10−22.31 × 10−21.68 × 10−21.44 × 10−2
Sum1.94 × 10−11.41 × 10−11.21 × 10−12.57 × 10−11.87 × 10−11.6 × 10−11.61 × 10−11.17 × 10−11 × 10−1
WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; V: Vanadium; Mn: Manganese; Cu: Copper; Zn: Zinc; Cd: Cadmium; Pb: Lead; Sum: Summation.
Table 8. Carcinogenic risk of heavy metals in PM2.5 posed to different groups in Kunming.
Table 8. Carcinogenic risk of heavy metals in PM2.5 posed to different groups in Kunming.
ElementsDRJDMWM
ChildrenMenWomenChildrenMenWomenChildrenMenWomen
Cr4.5 × 10−51.64 × 10−41.4 × 10−42.61 × 10−59.47 × 10−58.12 × 10−52.69 × 10−59.78 × 10−58.38 × 10−5
Co1.49 × 10−65.42 × 10−64.65 × 10−63.61 × 10−51.31 × 10−51.23 × 10−51.14 × 10−64.16 × 10−63.57 × 10−6
Ni4.09 × 10−71.49 × 10−61.27 × 10−67.84 × 10−72.85 × 10−62.44 × 10−62.82 × 10−71.03 × 10−68.78 × 10−7
As9.06 × 10−63.29 × 10−52.82 × 10−52.21 × 10−58.03 × 10−56.88 × 10−51.49 × 10−55.41 × 10−54.64 × 10−5
Cd9.75 × 10−73.55 × 10−63.04 × 10−61.97 × 10−67.18 × 10−66.15 × 10−67.7 × 10−72.8 × 10−62.4 × 10−6
Sum5.69 × 10−52.07 × 10−41.77 × 10−45.45 × 10−51.98 × 10−41.7 × 10−44.4 × 10−51.6 × 10−41.37 × 10−4
WM: West Mountain, Clean reference sampling site; JDM: Jinding Mountain, Industrial sampling site; DR: Dongfeng East Road, Traffic-intensive sampling site; Cr: Chromium; Co: Cobalt; Ni: Nickel; As: Arsenic; Cd: Cadmium; Sum: Summation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Han, X.; Li, S.; Li, Z.; Pang, X.; Bao, Y.; Shi, J.; Ning, P. Concentrations, Source Characteristics, and Health Risk Assessment of Toxic Heavy Metals in PM2.5 in a Plateau City (Kunming) in Southwest China. Int. J. Environ. Res. Public Health 2021, 18, 11004. https://doi.org/10.3390/ijerph182111004

AMA Style

Han X, Li S, Li Z, Pang X, Bao Y, Shi J, Ning P. Concentrations, Source Characteristics, and Health Risk Assessment of Toxic Heavy Metals in PM2.5 in a Plateau City (Kunming) in Southwest China. International Journal of Environmental Research and Public Health. 2021; 18(21):11004. https://doi.org/10.3390/ijerph182111004

Chicago/Turabian Style

Han, Xinyu, Shuai Li, Zezheng Li, Xiaochen Pang, Yuzhai Bao, Jianwu Shi, and Ping Ning. 2021. "Concentrations, Source Characteristics, and Health Risk Assessment of Toxic Heavy Metals in PM2.5 in a Plateau City (Kunming) in Southwest China" International Journal of Environmental Research and Public Health 18, no. 21: 11004. https://doi.org/10.3390/ijerph182111004

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