Next Article in Journal
Zircon Chemistry and Oxidation State of Magmas for the Duobaoshan-Tongshan Ore-Bearing Intrusions in the Northeastern Central Asian Orogenic Belt, NE China
Next Article in Special Issue
Asbestiform Amphiboles and Cleavage Fragments Analogues: Overview of Critical Dimensions, Aspect Ratios, Exposure and Health Effects
Previous Article in Journal
DFT Simulation of the Water Molecule Interaction with the (00l) Surface of Montmorillonite
Previous Article in Special Issue
Characterization of Particulates from Australian Underground Coal Mines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Health Risks and Source Analysis of Heavy Metal Pollution from Dust in Tianshui, China

1
Key Laboratory of Oasis Ecology, Xinjiang University, Urumqi 830046, China
2
Key Laboratory of Smart City and Environment Modelling of Higher Education Institute, College of Resources and Environment Sciences, Xinjiang University, Urumqi 830046, China
3
College of Resources and Environmental Engineering, Tianshui Normal University, Tianshui 741000, China
*
Author to whom correspondence should be addressed.
Submission received: 8 April 2021 / Revised: 3 May 2021 / Accepted: 7 May 2021 / Published: 9 May 2021

Abstract

:
The purpose of this study is to explore the degree and spatial distribution of dust heavy metal pollution in Tianshui City, the health risks, and the sources of heavy metals. The geoaccumulation index and health risk assessment are used to study pollution levels and human health risks, and Cu, Zn, and Pb pollution are found to be serious. The total exposure of children to dust and heavy metals is 8.329 × 10−3 mg·kg−1·d−1, which is 4.66 times that of adults. The effect of carcinogenic heavy metal exposure is more significant for adults than for children. The total non-carcinogenic risk quotient to children via multiple pathways is 2.1690, which is higher than that of adults. Children’s Pb non-carcinogenic risk quotient is 4.79 times that of adults, and children are more sensitive than adults to the health risks of Pb pollution. The GeoDetector and Unmix 6.0 models are used for source analysis, revealing that Zn, Pb, and As pollution originate primarily from urban transportation systems, V is sourced from soil-forming parent materials, and Mn, Ni, Cu, and Co arise from mixed sources. Therefore, the treatment of heavy metal pollution in cities needs to focus more on the urban transportation system.

Graphical Abstract

1. Introduction

Heavy metal pollution is an important component of urban environmental pollution research. Because the urban environment is greatly affected by population concentration, developed industry and commerce, traffic congestion, and human activities, urban heavy metals are not only high in content but also present in many types. Moreover, as toxic and harmful pollutants that are difficult to degrade, heavy metals entering the soil will not only cause soil quality degradation but will also negatively affect human health through food chain accumulation, inhalation of dust, and skin contact and induce deadly diseases such as cancer [1], especially for sensitive groups such as the elderly and children [2]; thus, research on urban heavy metal pollution is necessary.
The main body of urban heavy metal pollution research includes urban surface dust [3,4], urban soil [5], green spaces [6,7], and other urban components. The research content covers content analysis, spatial variability, risk assessment [8,9,10,11], and systematic and in-depth research on the temporal and spatial properties of urban heavy metal pollution. Among this research, a study on the hazards of heavy metals in urban dust found that heavy metals in dust pose a greater risk to human health than do heavy metals in urban soil [12,13,14,15]. Among types of urban dust, road dust is an important source and sink of urban heavy metals [16], which accumulate from various sources, such as soil, air, industrial production, transportation, and coal burning [16,17,18,19,20]. The distribution and source of dust on urban roads vary from city to city, showing different characteristics [14,17,21]. Previous studies on heavy metal pollution in different types of cities mainly considered industrial cities [22,23] and mixed-use cities [24,25]. The characteristics of urban heavy metal pollution and possible pollution sources have been studied, but the traditional diffusion model and receptor model were mainly used in research on developed large-scale cities. Insufficient research exists on small and medium-sized cities. In addition, the analysis of pollution sources is mainly based on principal component analysis, correlation analysis, and cluster analysis. It is difficult to analyze the categories of all pollution sources, and the results are not thorough enough. It is difficult to judge the cumulative effect of natural sources and environmental particles on the content of heavy metal elements. Therefore, further in-depth research on the sources of heavy metal pollution in urban dust needs to be conducted.
Among current source analysis methods, the Unmix model, a new source analysis method, has a visual graphical interface and diagnostic tools. In contrast to the chemical mass balance (CMB) model, the Unmix model does not need to determine the source component spectrum data in advance, thus avoiding component spectrum collection and the disadvantages of high difficulty and a heavy workload [6,7]. In contrast to the PMF model, the Unmix model does not need to set the number of pollution sources and does not need to know the uncertainty of the data, which reduces the impact caused by human factors. The model obtains the results by itself according to the selected components, and the uncertainty is shown in the analysis results [26]. The Unmix model is widely used in air pollutant source analysis [27], and its accuracy and comprehensiveness have allowed it to be gradually applied to soil and sediment [28,29]. Due to the extensive sources of heavy metals in urban dust [30], traceability analysis requires the identification and analysis of multifactor sources, and geographic detectors are highly practical for such problems. Geographic detectors are based on the spatial heterogeneity of geographic phenomena, and it is assumed that if geographic factor A is controlled by geographic factor B, then B will show a spatial distribution similar to A [31,32]. Therefore, the explanatory power of influencing factors on the differentiation of heavy metals can be judged through the correlation and similarity of the influencing factors and the airborne distribution of heavy metals [33,34], and the main influencing factors can be determined. This paper uses geographic detectors in traceability analysis to reveal the influence of multiple factors.
Tianshui City is a livable construction city in Gansu Province. However, due to the development of urbanization and industrialization in recent years, the discharge of dust pollutants in the environment has increased, which directly affects the environmental quality and the health of urban residents. Compared with the sporadic distribution of super-large and large cities, such as Beijing, Shanghai, Guangzhou, and provincial capital cities, the urban system constitutes the majority of small and medium cities. Therefore, this article chose Tianshui City, an underdeveloped area in western China, as the research area. Nine heavy metal elements, Cr, Mn, Ni, Cu, Zn, As, Pb, V, and Co, were selected as the research object. The pollution degree and spatial distribution characteristics of dust heavy metals were investigated, the health risk level was assessed, and GeoDetector and the Unmix model were used to analyze the dust heavy metal pollution sources in detail. The results reveal the pollution levels and sources of heavy metals in a tourist city and provide a scientific reference for dust environmental management and heavy metal pollution control in cities.

2. Materials and Methods

2.1. Sample Collection in the Study Area

This study took the main urban areas (Maiji District and Qinzhou District) of Tianshui City as the research object. The built-up area of Qinzhou District and Maiji District is 16.4 square kilometers. The sampling points were mainly divided by the grid method, supplemented by the serpentine method, in an attempt to make the sampling points evenly distributed in the grid. Taking into account the flow of people, traffic, and other factors, a total of 51 sampling sites were set up, covering residential areas, commercial squares, transportation hubs, hospitals, hotels, schools, and other places. In autumn 2018, sufficient surface dust samples were collected with plastic brushes and dustpans at various points near each site (approximately 10 m in radius from the center of the site), mixed, put into sample bags, and labelled for storage. A total of 51 dust samples were collected. The samples were sieved through a 180-mesh nylon sieve and then packaged in a sealed bag. The location of Tianshui City and the distribution of sampling points are shown in Figure 1. The sampling points in Qinzhou District are represented by A1, A2, A3, ... (25 in total), and those in Maiji District are represented by B1, B2, B3, ... (26 in total).

2.2. Laboratory Determination

The dust samples were transported to the laboratory, air-dried, passed through a 0.15 mm copper sieve, and mixed evenly. Fifteen grams of sample was ground with a tungsten carbide mill (ZHM-1A) to a particle size of less than 200 mesh. The ground samples were dried at 105 °C. A total of 4 g was placed in a sample preparation mold, boric acid edging backing was added, and a 30 t semiautomatic sample press (ZHY-401A) was used to press a sample with an inner diameter of 32 mm, after which the sample was placed in a dryer until testing. A sequential wavelength dispersion-type X-ray fluorescence spectrometer (model: Axios, place of production: Almelo, Netherlands) was used to determine the element content, and the analysis software was SuperQ Version 5.0. This study selected 9 heavy metal elements (Cr, Mn, Ni, Cu, Zn, As, Pb, V, and Co) for analysis.

2.3. Method

2.3.1. Geoaccumulation Index

The geoaccumulation index (Igeo) was proposed by German scientists in 1969 [35]. The geoaccumulation index method takes into account the impact of human activities, environmental geochemical background values, and natural diagenesis disturbances in the evaluated system. Therefore, the geoaccumulation index method is now often used in the evaluation of soil heavy metal pollution.
The calculation formula is as follows:
I g e o = l o g 2 ( C n / K B n )
In the formula, Igeo is the geoaccumulation index, Cn is the heavy metal concentration in the sample, and Bn is the background value of the local soil environment (due to the differences in the background values of heavy metal elements in different soil types, this article uses the heavy metal contents of limestone soil in the Longzhong area where Tianshui City is located). These element contents are used as the background values in this study [36]. K is a coefficient representing the variation in the background value caused by the difference in rocks in various places. The value of K is 1.5 in this paper [37]. The classification of the geoaccumulation index and the corresponding pollution degrees are shown in Table 1 [23,37].

2.3.2. Health Risk Assessment Method

Based on the soil health risk model proposed by the United States (US) Environmental Protection Agency (EPA) [38], the non-carcinogenic risk of 9 heavy metals at sampling points in the main urban area of Tianshui City was quantitatively estimated. Cr, Ni, and As also have carcinogenic risks. An assessment was conducted to reveal the health risks of heavy metals in dust.
The exposure calculation formula is as follows:
A D D L n g = C × L n g R × E F × E D B W × A T × 10 6 ,
A D D L n g = C × L n g R × E F × E D B W × A T × 10 6 ,
A D D d e r m a l = C × S A × S L × A B S × EF × ED BW × A T × 10 6 ,
L A D D L nh = C × C × EF PEF × A T × ( L nh R c h i l d × E D c h i l d B W c h i l d + L nh R a d u l t × E D a d u l t B W a d u l t ) .
where   A D D L n g , A D D L nh , and   A D D d e r m a l are the average daily dust exposure via the human digestive system, respiratory system, and skin contact, respectively.   L A D D L nh is the average daily dust exposure to carcinogenic heavy metals through the respiratory system. The unit of average daily dust exposure is mg/(kg·d); see Table 2 for the meaning and value of each parameter.
The non-carcinogenic and carcinogenic risks of heavy metals in dust are calculated based on the reference dose of chronic poisoning and the average lifetime daily exposure, respectively. The formulas are as follows:
H Q i = A D D / R f D i ,
H I = H Q i ,
R i s k i = L A D D i × S F i ,
R i s k T = R i s k i ,
where i is a specific heavy metal; H Q i is the non-carcinogenic risk quotient of heavy metal i through a single route; and A D D is the average daily dust. R f D i is the reference dose of a single route [38,44], which represents the maximum amount of pollutant that will not affect human health per unit time and unit body mass, mg/(kg·d).
H I is the sum of the non-carcinogenic risk quotients of multiple heavy metals and multiple exposure pathways; S F i is the slope coefficient of heavy metal I [38], which refers to the maximum probability of developing cancer due to human exposure to a certain dose of the heavy metal, mg/(kg·d); and R i s k is the carcinogenic risk quotient, indicating the possibility of cancer. If the value of R i s k is between 10−6 and 10−4 (that is, one cancer patient per 104 to 106 people), the substance is considered to have no obvious risk of cancer [20]. See Table 3 for the meaning and value of each parameter. It is generally believed that when H Q < 1 or H I < 1, the non-carcinogenic risk is small or negligible; when H Q ≥ 1 or H Q ≥ 1, it is considered that there is a non-carcinogenic risk.

2.3.3. GeoDetector

GeoDetector is a set of statistical analysis methods that detect the spatial heterogeneity of variables, reveal the correlation between variables, and reflect the driving force behind them [31,45]. In recent years, it has also been used to study the driving factors of the spatial variability of soil heavy metals [46,47]. The GeoDetector model is composed of four submodels: a factor detector, interaction detector, risk detector, and ecological detector. This article mainly uses the core module of GeoDetector, the factor detector, to quantitatively study the effect of natural and socioeconomic factors on the heavy metal loading of dust in the study area. The influence of differentiation and the calculation method of the factor detector q value are detailed in the literature [31].

2.3.4. Unmix

Unmix 6.0 is a multireceptor model developed by the US EPA. It aims to solve the general mixing problem; that is, the data are a linear combination of contributions from a number of unknown sources, and the contributions of these unknown components to each sample are unknown. In this model, the singular value decomposition (SVD) method is used to reduce the dimensionality of the data space to determine the number of sources, the source composition, and the contribution rate of each source to each sample [48]. The model formulas are detailed in the literature [48].

2.4. Other Data Sources

Considering that dust heavy metals mainly come from the soil-forming parent material and human activities, the source factors of heavy metals in this study mainly include the following categories: anthropogenic activity factors (daily average bus traffic (X1), population density (X2), distance from industrial and mining enterprises (X3), road density (X4), machine repair facility density (X5) and natural environmental factors (distance from the river (X6)), altitude (X7), slope (X8), Fe2O3 (X9), MgO (X10), pH (X11), and soil texture (X12)). Among all factors, X1 and X4 are factors related to the transportation system; X2, X3, and X5 are socioeconomic factors; X6, X7, and X8 are geographic factors, and X9, X10, X11, and X12 are soil parent material factors. X1 is based on a Baidu map query and actual statistical results. X2 comes from the Strategic Advanced Technology Products and Linkages Development Agency (https://www.worldpop.org/project/categories?id=18, accessed on 22 June 2020). X3 and X5 are extracted based on the place of interest (POI) data obtained from Baidu mapping and the results of buffer and nuclear density analysis. X4 is based on OpenStreetMap (www.openstreetmap.org, accessed on 8 April 2020) road data obtained by nuclear density analysis. X6 is obtained by buffer analysis based on OpenStreetMap river data. X7 and X8 are based on 30-m-resolution digital elevation model (DEM) data extracted from the geospatial data cloud hosted on the computer network of the Chinese Academy of Sciences Information Center (http://www.gscloud.cn/, accessed on 1 March 2021). Since X9 and X10 are indicator compounds of the soil parent material [49], they are used as influencing factors, and the data are based on the monitoring results for the samples in this study. X11 and X12 come from the China Ecosystem Assessment and Ecological Security Database (http://www.ecosystem.csdb.cn/index.jsp, accessed on 20 October 2020).

3. Results

3.1. Pollution Status Evaluation

3.1.1. Analysis of Pollution Statistics

The results of the geoaccumulation index of heavy metal elements at the sampling points of the study area are shown in Figure 2. From the location of the box-shaped distribution in the figure, there is a large difference in the heavy metal Igeo in the main urban area of Tianshui City. From the calculation results, the order of Igeo for each element in the dust is Zn > Pb > Cu > As > Cr > Mn > V > Ni > Co. Among these elements, Mn, V, and Co are in a non-polluting state, Ni is in a light–medium pollution state, and the rest are in a safe state. The pollution statuses of Cr and Ni are similar. Cu and As predominantly exhibit mild and mild–strong pollution, accounting for 80.7% and 53.8%, respectively, and the remaining samples are in a pollution-free state. The Igeo values of Zn and Pb range from 1.0224–5.6786 and 0.6899–4.7120, respectively, with mild–strong pollution as the main pollution degree. Zn pollution is more serious than Pb pollution. The average Igeo values for Cr, Mn, Ni, Cu, Zn, As, Pb, V, and Co are −0.0625, −0.7191, −0.8102, 0.5390, 2.2386, 0.1449, 2.1256, −0.8126, and −1.1417, among which Cu and As exhibit light and moderate pollution, and Zn and Pb exhibit light and strong pollution.

3.1.2. Spatial Distribution of Heavy Metal Pollution

The calculated Igeo value of each sample point was loaded into ArcGIS 10.2 and the Igeo index was spatially interpolated and mapped through Kriging analysis and histogram equalization. The results are shown in Figure 3.
It can be seen from Figure 3 that there are significant differences in the Igeo values of Cr, Mn, Ni, Cu, Zn, As, Pb, V, and Co in Qinzhou District (A) and Maiji District (B), and heavy metal pollution also exists in the same region. The Cr pollution in Qinzhou District is scattered, mainly concentrated in the eastern region in the areas near A21, A23, A26, A27, and A34. Cr pollution is concentrated in the Maiji District and distributed in patches. The pollution level in the central region is high, and the eastern and western regions are more polluted, with B3, B11, B12, B14, B19, and the surrounding sample sites as the main pollution points forming high pollution areas. Mn in Qinzhou District and Maiji District is characterized by strong pollution in the east and is concentrated in patches, and the western and central regions exhibit weak and strong local distribution characteristics, respectively. For Ni in Qinzhou district, A4, A10, A19, and A23 are the main high value points, and the accumulation at other sample sites is weak. The cumulative distribution of Ni in the Maiji area is similar to the distribution of Cr. The central part has a strong distribution, while the eastern and western parts have weak distributions. Cu pollution is scattered in Qinzhou District, mainly with A15, A19, A21, and A23 as strong pollution points, while Maiji District has a scattered and heavily polluted area with the western sites B9 and B5 as the main pollution locations and some polluted areas in the eastern part. Two strong pollution points are located at B24 and B27, and a weakly polluted area is located in the central part. The spatial distributions of Zn, As, and Pb pollution in Qinzhou District are similar, indicating that their sources are similar from west to east; A1, A33, A29, A12, A15, A19, A23, A26, and A27 constitute the high pollution area, while A2, A3, A4, A5, A6, and A7 form a low pollution group, and there are small differences in local samples. As a whole, the accumulation area presents a characteristic distribution of alternating strong and weak points from west to east. The pollution of As is considerably stronger than that of Zn and Pb in the central area. The Igeo values of V and Co reach the pollution level, but the spatial distributions of these elements in the Qinzhou and Maiji districts differ, indicating that there are large differences in their sources. The V content at sampling points A18, A19, A21, A23, A33, A34, B4, B9, B10, B17, B19, and B20 is high, and the Co content is high in samples A4, A24, A26, A27, A30, B3, B5, B11, B12, and B14, indicating that the possibility of contamination is high.

3.2. Health Risk Assessment

Based on the calculation of the non-carcinogenic exposure to heavy metals in dust via the digestive system, respiratory system, and skin contact and the carcinogenic exposure to some heavy metals (Table 4), the non-carcinogenic and carcinogenic risk parameters of dust and heavy metals (Table 5) were obtained. The risk value is calculated based on the average content. It should be noted that this article assumes that deposited dust is the source of inhaled dust to obtain the greatest possible health risk from inhaled deposited dust. However, the dust that naturally exists in the air is also an important source of dust inhaled by humans. In addition, this article does not consider the difference in absorbability due to the size of dust particles and uses the dust content as the inhalation amount to make a rough estimate of the risk assessment. Therefore, the results in this paper are overestimated, which means that they represent the greatest possible occurrence of health risks. The actual exposure through the respiratory system may be different from the calculated results, and the uncertainty of this assessment method needs further discussion. Although the results are exaggerated using the maximum possibility as the evaluation result, they have more value for risk warning for pollution prevention and control decision makers.
Table 4 shows that the total human exposure to heavy metals in dust follows the order As < Ni < Co < Cr < Pb < Cu < Mn < Zn < V and the exposure to non-carcinogenic heavy metals is greater than that to carcinogenic heavy metals. Among carcinogenic heavy metals, Cr has the largest exposure (65.59% of total carcinogenic heavy metal exposure), followed by Ni, and As has the least exposure (accounting for 13.99%). From the perspective of the main exposed population, the total exposure of children to heavy metals in dust is much greater than that of adults. The total exposure of children is 8.329 × 10−3 mg·kg−1·d−1, which is 4.66 times that of adults (1.786 × 10−3 mg·kg−1·d−1). The exposure of children to heavy metals in dust is more significant than that of adults. Children’s carcinogenic heavy metal exposure accounts for 0.903% of their total exposure, while adult carcinogenic heavy metal exposure accounts for 3.516%. The effect of carcinogenic heavy metal exposure is more significant for adults than for children; the proportions of Mn, Cu, Zn, Pb, V, and Co exposure relative to the total non-carcinogenic exposure are 38.86%, 3.82%, 39.25%, 12.20%, 4.42%, and 0.59%, respectively, for children and 38.38%, 3.77%, 38.77%, 12.04%, 4.37%, and 0.580% for adults. The effect of non-carcinogenic heavy metals on children and adults is relatively small. The exposure through the three different pathways for the different groups of people shows that the digestive system is responsible for much more exposure than skin contact and the respiratory system.
According to the risk quotients in Table 5, among children and adults, the non-carcinogenic risk quotient for exposure to all nine heavy metals via the digestive system is the highest (2.0885 for children, accounting for 96.29% of the total multipathway non-carcinogenic risk quotient, and 1.4702 for adults, accounting for 97.89%). The risk generated by the respiratory system is relatively small. The total multipathway non-carcinogenic risk quotient is 2.1690 for children and 1.5020 for adults, both of which are greater than 1, indicating that heavy metals from dust in the main urban area of Tianshui have strong non-carcinogenic effects on children and adults. The carcinogenic risk for adults is more significant than the non-carcinogenic risk, and the opposite trend is observed for children. The total non-carcinogenic risks of the nine heavy metals in descending order are Ni < Cu < Zn < Co < V<As < Mn < Pb < Cr. The worst situation is that all Cr in the dust is hexavalent. In this case, the non-carcinogenic risks of Cr to children and adults are 1.629 and 1.362, respectively, which are greater than the warning value of 1, and there is an obvious non-carcinogenic risk. The impact on human health is the strongest. However, in real conditions, the Cr is more likely to be trivalent chromium. Although excessive trivalent chromium in the human body can also cause poisoning, its toxicity is negligible compared with hexavalent chromium [50,51]. Therefore, the non-carcinogenic risk of Cr in this article is the maximum value in the theoretical sense rather than the actual risk value of dust.
This is followed by Pb, with a non-carcinogenic risk much greater than those of other elements. The non-carcinogenic risk quotient of Pb in children is 4.79 times that in adults; thus, children are more sensitive than adults to Pb pollution. Therefore, although the non-carcinogenic risk of Pb in children is lower than the warning value, there is the possibility of affecting health. Other heavy metals have a low non-carcinogenic risk quotient and have little impact on human health.
In this study, the total carcinogenic risk of the carcinogenic heavy metals Cr, Ni, and As is 5.515 × 10−9, while the total carcinogenic risk of single heavy metals and polymetallic metals is less than 106, a relatively small value.

4. Discussion

Determining the source of heavy metal pollution is very important for the prevention and control of heavy metal pollution. Therefore, this article first uses geographic detectors to detect the relevance of each influencing factor. Factor detection can analyze whether influencing factors have an impact on the distribution of heavy metals in dust and the size of the impact and by determining their relative importance, the main influencing factors can be screened. In this paper, factor detection was carried out on the content of heavy metals and the selected 12 factors, and the results are shown in Table 6.
The factor detection results show that most of the influencing factors passed the significance test, and the average explanatory power q was 0.1177. The distribution of heavy metals in the study area and the spatial distribution of influencing factors have macroscopic trends. The main influencing factors were extracted, and it was found that X10, X1, X12, X9, X2, and X7 have great explanatory power for the spatial distribution of nine heavy metals, with cumulative explanatory power values of 2.2747, 2.2598, 1.9255, 1.8406, 1.0291, and 0.9931. The content of heavy metals in the urban environment is the result of both natural background values and human activities. Therefore, from the results of factor detection analysis, the factors related to the parent material (X10, X12, X9, and X7) and socioeconomic factors (X1 and X2) have a significant impact on the content of heavy metals.
Ranking of the cumulative explanatory power of heavy metal elements by various influencing factors reveals that, in addition to the soil-forming parent material factor, the transportation system factor X1 is the factor with the most frequent occurrences and the strongest explanatory power (cumulative explanatory power of 2.2598); furthermore, the traffic system factor X4 also has strong explanatory power for the spatial differentiation of Cr, Ni, and Co content. The explanatory power of most heavy metals in terms of X3 and X5 is significant, but the values are all less than 0.1, at a weak level, not enough for the heavy metal content to exhibit strong space differentiation. The above results indicate that in addition to natural sources, the sources of heavy metals in the main urban area of Tianshui City are mainly from the urban transportation system and the cumulative heavy metals produced by this system are mainly due to daily operating activities such as automobile exhaust emissions and vehicle component wear.
The geographic detector effectively identified the main influencing factors of heavy metal content accumulation, but they were not enough to reveal the quantitative contribution rate of specific pollution sources to heavy metal content. Therefore, Unmix 6.0 was used to quantitatively analyze the source composition and, combined with the geographic detector factor detection results, determine the source of heavy metals.
Unmix 6.0 software was used for source analysis. When all nine heavy metals were analyzed, the software prompted that there was no solution. After removing Cr according to the software recommendations, the model ran normally, and the results show that the system has four sources and that the model Min R2 = 0.97. A total of 97% of the species variance can be explained by the model, which is greater than the minimum value required by the system (Min R2 > 0.8). Min Sig/Noise = 2.09, which is greater than the minimum value required by the system (Min Sig/Noise > 2). The analytical results from these four sources are credible. The results of geographic detection in the previous section indicate that the content of heavy metals is mainly affected by natural sources (soil-forming parent material) and transportation system factors. Therefore, the soil-forming parent material and transportation system factors were used as the main analytical factors in the source analysis process, and other factors were used as references in conjunction with the source contribution data. The similarity between the spatial distributions of sources and high-load pollutants was used to determine the sources of heavy metal pollution. The spatial distribution of source contributions is shown in Figure 4.
Source 1 has the highest Zn loading in the composition spectrum, reaching 83.7%. Source 1 causes Zn accumulation. Factor detection analysis shows that apart from natural sources of Zn, the daily average bus traffic has the strongest explanatory power, with the q value reaching 0.2024, indicating that operation of the transportation system has the greatest influence on the accumulation of Zn. Zn-containing castings are important materials for automobile transmission parts, engine parts, and body components in the automobile industry. Zn-containing roadside soil produced by vehicle wear and tear accumulates in surface dust, causing Zn pollution [45]. Moreover, the contribution rate of source 1 to sample points 29, 39, and 44 in Figure 5 is significantly greater than that for other sample points. These sample sites are adjacent to the main roads in Maiji District. These roads are the first-class roads in the region and bear the main passenger and freight transportation in the city. The roads have a large traffic capacity and are a section that is prone to congestion during rush hours. Passenger and freight vehicles emit large amounts of exhaust gas, mechanical friction and tire wear are serious, and Zn accumulation is serious. Studies have shown that Zn may come from roofs and gutters [52,53], which may be related to the Zn-containing materials used in roof construction, especially galvanized materials and coatings [54]. Zn accumulates in drains along with roof runoff. However, the roofing materials used in Tianshui City are mainly expanded polystyrene boards, polymer cement waterproof coatings, and asphalt. There are fewer Zn-containing materials, and the possibility of Zn coming from roofs and drainage ditches is less likely. Therefore, source 1 represents the accumulation of heavy metals generated by the wear of automobile parts in the transportation system.
The loading ratios of Mn, Zn, and Pb in the Source 2 composition spectrum are 28.1%, 32.0%, and 28.0%, respectively. Source 2 is the primary source of Pb pollution and has a dominant contribution to Pb content. In the transportation system, automobile exhaust is a direct source of Pb [55]. At present, the fuel used by Chinese cars is #92 or #95 gasoline and #0 diesel. Pb is the most important heavy metal in the exhaust of four vehicle types [49]. The heavy metals in exhaust gas also include Mn, Zn, Ni, As, Cr, Co, and others. Existing studies have shown that the content of Mn in exhaust gas from #92 and #95 gasoline is significantly higher than that of other heavy metals [55]. Exhaust gas is an important source of Mn in dust. Therefore, it is believed that source 2 is automobile exhaust in the transportation system.
The load of Mn in the composition spectrum of source 3 is 61.8%, the Cu load is 12.2%, and the V load is 6.8%; these values are far greater than the contribution rates of other sources to these three heavy metals, so source 3 is the dominant source of the above three elements. The Zn and Pb loadings in source 3 are both less than 2%, indicating that traffic system pollution sources account for a relatively small proportion in source 3; thus, traffic system pollution sources are excluded. The factor detection analysis results (Table 6) show that the primary factors leading to the accumulation of Mn, Cu, and V are the soil parent material factors (soil texture or Fe2O3) and that the main influencing factors of V are the soil parent material factors and other factors. Natural factors are characterized by strong natural sources. These three elements mainly come from natural environmental conditions, and thus source 3 can be identified as a natural source.
Source 4 has the largest loading of Mn and Zn in the composition spectrum, but its contribution rate to the accumulation of Ni is the highest among the four sources, which is 3.3%. In a study of the sources of Ni in cities, it was found that Ni was mainly derived from natural sources such as soil-forming parent material [56] and industrial activities [57]. The factor detection results (Table 6) show that the content of the soil-forming parent material factor MgO and socioeconomic factors such as daily average bus traffic, road density, distance to industrial and mining enterprises, and population density have significant explanatory power for the accumulation of Ni. The cumulative explanatory power of the soil-forming parent material factor and the socioeconomic factor are 0.6485 and 0.5198, respectively. The difference in influence intensity is small. The soil parent material and urban socioeconomic activities, including industrial activities, are contributing factors.
Regarding the high contribution rate of Mn and Zn, it can be seen from the factor detection results that Mn and Zn cover natural and anthropogenic sources. When the dominant source is excluded, the cumulative effect of heavy metals produced by other mixed sources also appears. In the source 4 component spectrum, the Cu load is 2.6%, second only to that of source 3.
However, existing urban heavy metal traceability studies have shown that the sources of Cu are diverse, including factors such as man-made emissions [58], brake materials [59], construction activities [55], and atmospheric dust fall [60]. The factor detection results (Table 6) indicate that X1 and X2 have strong explanatory power for Cu accumulation. It is believed that Cu is not only from natural sources but also from urban construction and development and urban daily economic activities. Therefore, source 4 is considered to be a mixed source composed of natural and manmade sources.
The loading of Co in each source component spectrum shows that the loadings of Co in source 3 and source 4 are both 0.009, so source 3 and source 4 are the pollution sources of Co. It is generally believed that in addition to natural accumulation during soil formation, Co mainly comes from industrial production, such as alloy production, electroplating, glass manufacturing, dyeing and other industries, with a wide range of sources, so Co has certain loadings in natural sources and mixed sources.
Since Cr was excluded in the Unmix model analysis, principal component analysis was used to assess the traceability of Cr, and the results were used to verify the above traceability analysis results. Three principal components were extracted by SPSS principal component analysis, and the cumulative contribution rate of the three principal components was 78.088%. The three principal components basically reflect the variability of the nine elements studied.
It can be seen from Table 7 that the elements with higher loads in Factor 1 are Zn, Pb, and As; the elements with higher loads in Factor 2 are Cr, Ni, and Co; and the elements with higher loads in Factor 3 are Mn, V, and Co. Since heavy metals with higher factor loadings under the same principal component have the same source, according to the source types of heavy metal elements in the analysis results of the Unmix model, Factor 1 can be identified as the impact of the transportation system, Factor 2 is the impact of natural and manmade mixed factors, and Factor 3 is the influence of natural factors. The distribution of each element among Factors and sources is consistent, indicating that the division of source types is credible.
In the principal component analysis results for Cu, Factor 1 and Factor 2 have significant effects on Cu, and the eigenvalues of Factor 1 and Factor 2 differ by 0.047, which is at a low level. Therefore, it is believed that the accumulation of Cu comes from transportation system factors and mixed factors. The Unmix model analysis results show that Cu comes from natural sources and mixed sources. In view of the diversity of Cu sources, the analysis results are considered credible, and the principal component analysis results can be used as supplements. The sources of Cu include transportation systems, natural sources, and mixed sources.
Regarding our final research goal, if we want to effectively prevent heavy metal pollution, we need to clearly identify the sources of heavy metals. At present, with regard to the composition of mixed sources, except for soil parent materials and traffic factors, the impact of other human activities is not clear. Land use type, as an important form of human activity, has a significant impact on the accumulation of heavy metals. The study area in this paper is the main urban area of Tianshui City. The land use in the study area is mainly for transportation, residential land, cultural and educational land, park green space, commercial land, and industrial land [14,17,22,24]. To meet the needs of people’s daily activities, the types and functions of land use are divided, the functions are diverse, and the boundaries are unknown. At the same time, using planar land-use patterns to reveal the sources of heavy metal pollution at individual sampling points, results in a certain degree of uncertainty.
Therefore, this paper collects POI data in the study area, replaces land use types with POI functional attributes, and uses elements (Cr, Ni, Cu, Co) mainly from mixed sources as research cases. A buffer with a radius of 160 m (the minimum distance between sampling points) is created around the sampling points, POI samples are extracted in the buffer, and we analyze the correlation between the number of POIs of different functions and the pollution index, exploring the impact of land use. The correlation analysis results in SPSS show that the correlation between the Ni and Co contents and the number of industrial POIs passed was significant at the 0.01 and 0.05 levels, although the correlation coefficients were 0.38 and 0.34, respectively, i.e., less than 0.6. Moreover, Cr and Cu failed the significance test, and there was no correlation. Therefore, the type of land use has little effect on the existence of heavy metals from mixed sources, and it appears as a single weak effect on industrial land. This may be another manifestation of the small explanatory power (0.083 and 0.0396) of the distance from the factory (X3) to Ni and Co in the geographic detection results.

5. Conclusions

First, the accumulation characteristics and health risks of heavy metals in the main urban areas of Tianshui City were evaluated by the geoaccumulation index and health risk index, and heavy metals were determined. Then, GeoDetector, the Unmix model, and principal component analysis were used to comprehensively identify the sources of heavy metals. The research results show that there are significant differences in the spatial distribution of different heavy metal pollution levels in Qinzhou District and Maiji District of Tianshui City and that there are also large differences in the pollution levels of different heavy metals in the same area, with serious Cu, Zn, and Pb pollution. The exposure to non-carcinogenic heavy metals in the study area is greater than that to carcinogenic heavy metals. The exposure to Cr is the largest. The total exposure of children to heavy metals in dust is much greater than that of adults, which is consistent with the conclusions of other researchers [61,62]. The exposure of adults to carcinogenic heavy metals is more significant than that of children. The exposure effects of non-carcinogenic heavy metals to different groups are different. The digestive system is the main pathway of exposure and has the highest non-carcinogenic risk quotient. Dust and heavy metals have a strong non-carcinogenic risk to children and adults; children have a more significant non-carcinogenic risk than carcinogenic risk, and they are more sensitive to the health risk of Pb pollution. Under the same pollution conditions, children’s hand-to-mouth frequency (playing and unconscious behaviors) and respiration rate per unit weight are higher than adults [63,64]; thus, children face greater health risks, and special attention should be given to their health. Furthermore, in view of children’s sensitivity to Pb pollution, when severe Pb accumulation occurs in cities, it is recommended to test Pb levels in children’s blood.
The traceability analysis results show that the soil parent material factor has the strongest explanatory power for the accumulation of heavy metal elements, followed by the transportation system factor. Industrial and other social factors have explanatory power but are at a low level. Zn, Pb, and As mainly come from the transportation system; V comes from the soil-forming parent material; and Mn, Ni, Cu, and Co come from a mixed source composed of anthropogenic and natural sources. The anthropogenic source is mainly composed of population activity intensity, urban construction, urban economic activities, and industrial production, reflecting the impact of a single industrial system. Compared with cities where heavy metal pollution originates primarily from industrial and agricultural production [17,22,47,61,62,65], the heavy metal pollution in the main urban area of Tianshui City arises mainly from the transportation system. For this type of city, in which urban traffic is the main source of heavy metal pollution, more efforts may be needed to optimize the layout and power of the urban transportation system, promote clean energy vehicles, reduce vehicle emissions and brake wear, and control the accumulation of heavy metal elements.
In addition, some key sources of heavy metals may vary in time. Therefore, it is necessary to continue to explore the time-varying characteristics of heavy metal pollution and sources to adopt more targeted prevention and control strategies to reduce the risk of heavy metals in dust [66].

Author Contributions

Conceptualization, B.T. and H.W.; methodology, B.T. and X.W.; software, C.M.; investigation, B.T. and X.W.; visualization, B.T. and J.Z.; writing—review and editing, B.T.; supervision, X.D.; funding acquisition, H.W. and X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (grant numbers 41861037).

Data Availability Statement

Third party data, not applicable.

Acknowledgments

We thank the editor and the reviewers for their insightful and valuable suggestions, which greatly improved the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hu, B.; Jia, X.; Hu, J.; Xu, D.; Xia, F.; Li, Y. Assessment of heavy metal pollution and health risks in the soil-plant-human system in the Yangtze River Delta, China. Int. J. Environ. Res. Public Health 2017, 14, 1042. [Google Scholar] [CrossRef]
  2. Yang, Q.; Li, Z.; Lu, X.; Duan, Q.; Huang, L.; Bi, J. A review of soil heavy metal pollution from industrial and agricultural regions in China: Pollution and risk assessment. Sci. Total Environ. 2018, 642, 690–700. [Google Scholar] [CrossRef] [PubMed]
  3. Qadeer, A.; Saqib, Z.A.; Ajmal, Z.; Xing, C.; Khalil, S.K.; Usman, M.; Huang, Y.; Bashir, S.; Ahmad, Z.; Ahmed, S.; et al. Concentrations, pollution indices and health risk assessment of heavy metals in road dust from two urbanized cities of Pakistan: Comparing two sampling methods for heavy metals concentration. Sustain. Cities Soc. 2020, 53, 101959. [Google Scholar] [CrossRef]
  4. Yang, S.; Liu, J.; Bi, X.; Ning, Y.; Qiao, S.; Yu, Q.; Zhang, J. Risks related to heavy metal pollution in urban construction dust fall of fast-developing Chinese cities. Ecotoxicol. Environ. Saf. 2020, 197, 110628. [Google Scholar] [CrossRef] [PubMed]
  5. Morton-Bermea, O.; Hernández-Álvarez, E.; González-Hernández, G.; Romero, F.; Lozano, R.; Beramendi-Orosco, L.E. Assessment of heavy metal pollution in urban topsoils from the metropolitan area of Mexico City. J. Geochem. Explor. 2009, 101, 218–224. [Google Scholar] [CrossRef]
  6. Liu, L.; An, Y.; Ma, J.; Chen, Y.; Wu, Y.; Li, J.; Huang, J. Analysis of the source of heavy metals in Beijing urban parks based on UNMIX model. Environ. Sci. Res. 2020. [Google Scholar] [CrossRef]
  7. Liu, L.; Liu, Q.; Ma, J.; Wu, H.; Qu, Y.; Gong, Y.; Yang, S.; An, Y.; Zhou, Y. Heavy metal(loid)s in the topsoil of urban parks in Beijing, China: Concentrations, potential sources, and risk assessment. Environ. Pollut. 2020, 260, 114083. [Google Scholar] [CrossRef] [PubMed]
  8. Al-Momani, I.F. Trace elements in street and household dusts in Amman, Jordan. Soil Sediment Contam. 2007, 16, 485–496. [Google Scholar] [CrossRef]
  9. Hogervorst, J.; Plusquin, M.; Vangronsveld, J.; Nawrot, T.; Cuypers, A.; Van Hecke, E.; Roels, H.A.; Carleer, R.; Staessen, J.A. House dust as possible route of environmental exposure to cadmium and lead in the adult general population. Environ. Res. 2007, 103, 30–37. [Google Scholar] [CrossRef]
  10. Maharachpong, N.; Geater, A.; Chongsuvivatwong, V. Environmental and childhood lead contamination in the proximity of boat-repair yards in Southern Thailand—I: Pattern and factors related to soil and household dust lead levels. Environ. Res. 2006, 101, 294–303. [Google Scholar] [CrossRef]
  11. Naimabadi, A.; Gholami, A.; Ramezani, A.M. Determination of heavy metals and health risk assessment in indoor dust from different functional areas in Neyshabur, Iran. Indoor Built Environ. 2020. [Google Scholar] [CrossRef]
  12. Christoforidis, A.; Stamatis, N. Heavy metal contamination in street dust and roadside soil along the major national road in Kavala’s Region, Greece. Geoderma 2009, 151, 257–263. [Google Scholar] [CrossRef]
  13. Dehghani, S.; Moore, F.; Keshavarzi, B.; Hale, B.A. Health risk implications of potentially toxic metals in street dust and surface soil of Tehran, Iran. Ecotoxicol. Environ. Saf. 2017, 136, 92–103. [Google Scholar] [CrossRef] [PubMed]
  14. Khademi, H.; Gabarrón, M.; Abbaspour, A.; Martínez-Martínez, S.; Faz, A.; Acosta, J.A. Environmental impact assessment of industrial activities on heavy metals distribution in street dust and soil. Chemosphere 2019, 217, 695–705. [Google Scholar] [CrossRef] [PubMed]
  15. Yadav, I.C.; Devi, N.L.; Singh, V.K.; Li, J.; Zhang, G. Spatial distribution, source analysis, and health risk assessment of heavy metals contamination in house dust and surface soil from four major cities of Nepal. Chemosphere 2019, 218, 1100–1113. [Google Scholar] [CrossRef] [PubMed]
  16. Duzgoren-Aydin, N.S.; Wong, C.S.C.; Aydin, A.; Song, Z.; You, M.; Li, X.D. Heavy metal contamination and distribution in the urban environment of Guangzhou, SE China. Environ. Geochem. Health 2006, 28, 375–391. [Google Scholar] [CrossRef]
  17. Li, Z.; Feng, X.; Li, G.; Bi, X.; Zhu, J.; Qin, H.; Dai, Z.; Liu, J.; Li, Q.; Sun, G. Distributions, sources and pollution status of 17 trace metal/metalloids in the street dust of a heavily industrialized city of central China. Environ. Pollut. 2013, 182, 408–416. [Google Scholar] [CrossRef]
  18. Liu, C.H.; Cen, K. Chemical composition and possible sources of elements in street dusts in Beijing. Acta Sci. Circumstantiae 2007, 27, 1181–1188. [Google Scholar] [CrossRef]
  19. Zhao, X.; Li, Z.; Wang, D.; Tao, Y.; Qiao, F.; Lei, L.; Huang, J.; Ting, Z. Characteristics, source apportionment and health risk assessment of heavy metals exposure via household dust from six cities in China. Sci. Total Environ. 2020, 143126. [Google Scholar] [CrossRef]
  20. Zheng, N.; Liu, J.; Wang, Q.; Liang, Z. Health risk assessment of heavy metal exposure to street dust in the zinc smelting district, Northeast of China. Sci. Total Environ. 2010, 408, 726–733. [Google Scholar] [CrossRef] [PubMed]
  21. De Miguel, E.; Llamas, J.F.; Chacón, E.; Berg, T.; Larssen, S.; Røyset, O.; Vadset, M. Origin and patterns of distribution of trace elements in street dust: Unleaded petrol and urban lead. Atmos. Environ. 1997, 31, 2733–2740. [Google Scholar] [CrossRef]
  22. Sabouhi, M.; Ali-Taleshi, M.S.; Bourliva, A.; Nejadkoorki, F.; Squizzato, S. Insights into the anthropogenic load and occupational health risk of heavy metals in floor dust of selected workplaces in an industrial city of Iran. Sci. Total Environ. 2020, 744, 140762. [Google Scholar] [CrossRef]
  23. Sun, L.; Guo, D.; Liu, K.; Meng, H.; Zheng, Y.; Yuan, F.; Zhu, G. Levels, sources, and spatial distribution of heavy metals in soils from a typical coal industrial city of Tangshan, China. Catena 2019, 175, 101–109. [Google Scholar] [CrossRef]
  24. Acosta, J.A.; Faz, A.; Martinez-Martinez, S. Identification of heavy metal sources by multivariable analysis in a typical Mediterranean City (SE Spain). Environ. Monit. Assess. 2009, 169, 519–530. [Google Scholar] [CrossRef]
  25. Weerasundara, L.; Magana-Arachchi, D.N.; Ziyath, A.M.; Goonetilleke, A.; Vithanage, M. Health risk assessment of heavy metals in atmospheric deposition in a congested city environment in a developing country: Kandy City, Sri Lanka. J. Environ. Manag. 2018, 220, 198–206. [Google Scholar] [CrossRef]
  26. Zhang, C. Research on Source Resolution Technology Based on UNMIX Model. Master’s Thesis, Nankai University, Tianjin, China, 2014. [Google Scholar]
  27. Peter, A.E.; Nagendra, S.M.S.; Nambi, I.M. Comprehensive analysis of inhalable toxic particulate emissions from an old municipal solid waste dumpsite and neighborhood health risks. Atmos. Pollut. Res. 2018, 9, 1021–1031. [Google Scholar] [CrossRef]
  28. Guan, Q.; Zhao, R.; Pan, N.; Wang, F.; Yang, Y.; Luo, H. Source apportionment of heavy metals in farmland soil of Wuwei, China: Comparison of three receptor models. J. Clean. Prod. 2019, 237, 117792. [Google Scholar] [CrossRef]
  29. Li, F.; Liu, S.Y.; Li, Y.; Shi, Z. Spatiotemporal variability and source apportionment of soil heavy metals in a industrially developed city. Huan Jing Ke Xue 2019, 40, 934–944. [Google Scholar] [CrossRef]
  30. Zhang, C. Using multivariate analyses and GIS to identify pollutants and their spatial patterns in urban soils in Galway, Ireland. Environ. Pollut. 2006, 142, 501–511. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, J.F.; Li, X.H.; Christakos, G.; Liao, Y.L.; Zhang, T.; Gu, X.; Zheng, X.Y. Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun Region, China. Int. J. Geogr. Inf. Sci. 2010, 24, 107–127. [Google Scholar] [CrossRef]
  32. Wang, J.F.; Zhang, T.L.; Fu, B.J. A measure of spatial stratified heterogeneity. Ecol. Indic. 2016, 67, 250–256. [Google Scholar] [CrossRef]
  33. Gao, J.; Wang, H. Temporal analysis on quantitative attribution of Karst soil erosion: A case study of a peak-cluster depression basin in Southwest China. Catena 2019, 172, 369–377. [Google Scholar] [CrossRef]
  34. Zhou, C.; Chen, J.; Wang, S. Examining the effects of socioeconomic development on fine particulate matter (PM2.5) in China’s cities using spatial regression and the geographical detector technique. Sci. Total Environ. 2018, 619–620, 436–445. [Google Scholar] [CrossRef]
  35. Muller, G. Index of geoaccumulation in sediments of the Rhine River. J. Geol. 1969, 2, 108–118. [Google Scholar]
  36. Wang, S.; Lian, B. Characteristics and distribution of soil environmental background values in Gansu Province. Gansu Environ. Res. Monit. 1993, 6, 1–7. [Google Scholar]
  37. Guo, X.; Liu, C.; Zhu, Z.; Wang, Z.; Li, J. Evaluation methods for soil heavy metals contamination: A review. Chin. J. Ecol. 2011, 30, 889–896. [Google Scholar]
  38. US EPA. Risk Assessment Guidance for Superfund (Rags), Volume I: Human Health Evaluation Manual (Part E, Supplemental Guidance for Dermal Risk Assessment) Interim; Office of Solid Waste Emergency Response, Waste and Clean-Up Risk Assessment: Washington, DC, USA, 2009.
  39. US EPA. Risk Assessment Guidance for Superfund Volume I: Human Health Evaluation Manual Part. F, Supplemental Guidance for Inhalation Risk Assessment; Office of Superfund Remediation and Technology Innovation: Washington, DC, USA, 2009.
  40. Wang, Z.S.; Duan, X.L.; Liu, P.; Jing, N.; Huang, N.; Zhang, J.L. Discussion on Chinese residents’ exposure parameters in environmental health risk assessment. Environ. Sci. Res. 2009, 1164–1170. [Google Scholar] [CrossRef]
  41. 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]
  42. Li, R.Z.; Zhang, A.J.; Tong, F.; Wu, Y.; Zhang, P.; Yu, J. Distribution characteristics and environmental health risk assessment of heavy metals in surface dust in Hefei City. Environ. Sci. 2011, 2661–2668. [Google Scholar] [CrossRef]
  43. Li, Z.S.; Wu, T.; Duan, X.L.; Wang, S.; Zhang, W.; Wu, X.; Yu, Y. Research on the exposure parameters of Chinese residents’ respiratory rate in environmental health risk assessment. Environ. Sci. Res. 2009, 1171–1175. [Google Scholar] [CrossRef]
  44. Beijing Research Institute of Environmental Protection. DB 11 T656—2009 Site Environmental Assessment Guidelines; Beijing Municipal Bureau of Quality and Technical Supervision: Beijing, China, 2009. [Google Scholar]
  45. Wang, J.F.; Hu, Y. Environmental health risk detection with GeogDetector. Environ. Model. Softw. 2012, 33, 114–115. [Google Scholar] [CrossRef]
  46. Jiang, X.T.; Wang, Q.; Li, R. Investigating factors affecting carbon emission in China and the USA: A perspective of stratified heterogeneity. J. Clean. Prod. 2018, 199, 85–92. [Google Scholar] [CrossRef]
  47. Shi, T.; Hu, Z.; Shi, Z.; Guo, L.; Chen, Y.; Li, Q.; Wu, G. Geo-detection of factors controlling spatial patterns of heavy metals in urban topsoil using multi-source data. Sci. Total Environ. 2018, 643, 451–459. [Google Scholar] [CrossRef]
  48. Henry, R.C. Multivariate receptor modeling by N-dimensional edge detection. Chemom. Intell. Lab. Syst. 2003, 65, 179–189. [Google Scholar] [CrossRef]
  49. Wilson, M.J. Weathering of the primary rock-forming minerals: Processes, products and rates. Clay Miner. 2004, 39, 233–266. [Google Scholar] [CrossRef]
  50. Sharma, N.; Sodhi, K.K.; Kumar, M.; Singh, D.K. Heavy metals eco-toxicity with major concern to chromium and recent advancement in remediation technologies. Environ. Nanotechnol. Monit. Manag. 2020, 15, 100388. [Google Scholar] [CrossRef]
  51. Ccaa, B.; Mebg, C.; Ahz, B. Hexavalent chromium: Regulation and health effects. J. Trace Elem. Med. Biol. 2021. [Google Scholar] [CrossRef]
  52. Chang, M.; Mcbroom, M.W.; Beasley, R.S. Roofing as a source of nonpoint water pollution. J. Environ. Manag. 2005, 73, 307–315. [Google Scholar] [CrossRef]
  53. Gromaire, M.C.; Robert-Sainte, P.; Bressy, A.; Saad, M.; Gouvello, B.D.; Chebbo, G. Zn and Pb emissions from roofing materials--modelling and mass balance attempt at the scale of a small urban catchment. Water Sci. Technol. 2011, 63, 2590–2597. [Google Scholar] [CrossRef] [PubMed]
  54. Hackbarth, K.; Lorenz, W.; Wilke, E.; Roth, M.; Wark, R.; Müller, S. Composition and Process for Coating Metal Surfaces; Henkel AG & Co. KGaA: Dusseldorf, Germany, 2009. [Google Scholar]
  55. Wang, S.; Cai, L.M.; Wen, H.H.; Luo, J.; Wang, Q.S.; Liu, X. Spatial distribution and source apportionment of heavy metals in soil from a typical county-level city of Guangdong Province, China. Sci. Total Environ. 2019, 655, 92–101. [Google Scholar] [CrossRef] [PubMed]
  56. Meng, Q. Research on the Removal Method of Heavy Metal Pollutants in Automobile Exhaust. Master’s Thesis, Shenyang Normal University, Shenyang, China, 2017. [Google Scholar]
  57. Liu, R.; Bao, K.; Yao, S.; Yang, F.; Wang, X. Ecological risk assessment and distribution of potentially harmful trace elements in lake sediments of Songnen Plain, NE China. Ecotoxicol. Environ. Saf. 2018, 163, 117–124. [Google Scholar] [CrossRef] [PubMed]
  58. Gong, M.; Wu, L.; Bi, X.Y.; Ren, L.M.; Wang, L.; Ma, Z.D.; Bao, Z.Y.; Li, Z.G. Assessing heavy-metal contamination and sources by GIS-based approach and multivariate analysis of urban–rural topsoils in Wuhan, Central China. Environ. Geochem. Health 2010, 32, 59–72. [Google Scholar] [CrossRef]
  59. Gultekin, D.; Uysal, M.; Aslan, S.; Alaf, M.; Guler, M.O.; Akbulut, H. The effects of applied load on the coefficient of friction in Cu-MMC brake pad/Al-SiCp MMC brake disc system. Wear 2010, 270, 73–82. [Google Scholar] [CrossRef]
  60. Lindström, M. Urban land use influences on heavy metal fluxes and surface sediment concentrations of small lakes. Water Air Soil Pollut. 2001, 126, 363–383. [Google Scholar] [CrossRef]
  61. Keshavarzi, B.; Tazarvi, Z.; Rajabzadeh, M.A.; Najmeddin, A. Chemical speciation, human health risk assessment and pollution level of selected heavy metals in urban street dust of Shiraz, Iran. Atmos. Environ. 2015, 119, 1e10. [Google Scholar] [CrossRef]
  62. Shahab, A.; Zhang, H.; Ullah, H.; Rashid, A.; Xiao, H. Pollution characteristics and toxicity of potentially toxic elements in road dust of a tourist city, Guilin, china: Ecological and health risk assessment. Environ. Pollut. 2020, 115419. [Google Scholar] [CrossRef] [PubMed]
  63. Jiang, Y.; Chao, S.; Liu, J.; Yang, Y.; Chen, Y.; Zhang, A.; Cao, H. Source apportionment and health risk assessment of heavy metals in soil for a township in Jiangsu Province, China. Chemosphere 2017, 168, 1658–1668. [Google Scholar] [CrossRef]
  64. Mohmand, J.; Eqani, S.A.M.A.S.; Fasola, M.; Alamdar, A.; Mustafa, I.; Ali, N.; Liu, L.; Peng, S.; Shen, H. Human exposure to toxic metals via contaminated dust:Bio-accumulation trends and their potential risk estimation. Chemosphere 2015, 132, 142–151. [Google Scholar] [CrossRef] [PubMed]
  65. Jun, Z.; Jie, D.; Qingfang, L.; Ningning, Y.; Yani, G. Heavy Metal Pollution Characteristics and Influencing Factors in Baoji Arban Soils. Environ. Sci. 2019, 8, 3774–3784. [Google Scholar] [CrossRef]
  66. Men, C.; Liu, R.; Xu, L.; Wang, Q.; Shen, Z. Source-specific ecological risk analysis and critical source identification of heavy metals in road dust in Beijing, china. J. Hazard. Mater. 2019, 388, 121763. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Dust sampling sites at sampling stations in Tianshui City. (a) Qinzhou District; (b) Maiji District.
Figure 1. Dust sampling sites at sampling stations in Tianshui City. (a) Qinzhou District; (b) Maiji District.
Minerals 11 00502 g001
Figure 2. Distribution of the geoaccumulation indexes of heavy metals.
Figure 2. Distribution of the geoaccumulation indexes of heavy metals.
Minerals 11 00502 g002
Figure 3. Spatial distribution of urban dust heavy metal geoaccumulation index values and pollution types in Tianshui City. (A) Qinzhou District; (B) Maiji District.
Figure 3. Spatial distribution of urban dust heavy metal geoaccumulation index values and pollution types in Tianshui City. (A) Qinzhou District; (B) Maiji District.
Minerals 11 00502 g003
Figure 4. Contribution rates of different sources to heavy metal pollution in dust.
Figure 4. Contribution rates of different sources to heavy metal pollution in dust.
Minerals 11 00502 g004
Figure 5. Source contribution spatial distribution. (A) Qinzhou district; (B) Maiji district.
Figure 5. Source contribution spatial distribution. (A) Qinzhou district; (B) Maiji district.
Minerals 11 00502 g005
Table 1. Geoaccumulation index classification and pollution degree.
Table 1. Geoaccumulation index classification and pollution degree.
Igeo Value≤00–11–22–33–44–5≥5
Contamination degreeSafeMild–moderate pollutionMild–strong pollutionModerate pollutionStrong pollutionStrong–heavy pollutionHeavy pollution
Table 2. Calculation parameters of daily average dust exposure.
Table 2. Calculation parameters of daily average dust exposure.
ItemParameterMeaningUnitValueData Sources
ChildAdult
Basic parametersCMass fraction of heavy metalsmg/kgAverage value of study areaAverage value of study areaThis study
Exposure behavioral parametersEDExposure perioda624[39]
BWAverage body masskg20.358.5[40]
EFExposure frequencyd/a180180[41]
AT (non-carcinogenic risk)Action timed365 × ED365 × ED[20]
AT (carcinogenic risk)Action timed365 × 70365 × 70[20]
Digestive tractIngRDust intakemg/d250150[42]
InhalationInhRAir intakem3/d8.612.9[43]
Skin contactPEFParticulate emission factorm3/kg1.36 × 1091.36 × 109[20]
SLSkin attachment factormg/cm211[38]
SAExposed skin areacm2/d11502145[42]
ABSInhalation factor0.0010.001[41]
Table 3. Reference dose and slope factor of different exposure pathways of heavy metals (mg·kg−1·d−1).
Table 3. Reference dose and slope factor of different exposure pathways of heavy metals (mg·kg−1·d−1).
ItemCrMnNiCuZnAsPbVCo
RfDIng2.86 × 10−54.60 × 10−22.06 × 10−24.02 × 10−23.00 × 10−13.01 × 10−43.52 × 10−37 × 10−32 × 10−2
RfDInh3.00 × 10−31.40 × 10−52.00 × 10−24.00 × 10−23.00 × 10−13.00 × 10−43.50 × 10−37 × 10−31.6 × 10−2
RfDdermal6.00 × 10−51.80 × 10−35.40 × 10−31.20 × 10−26.00 × 10−21.23 × 10−43.25 × 10−47 × 10−35.06 × 10−6
SFInh6.30.8415.1
Table 4. Dust heavy metal exposure doses at bus stations in Tianshui City (mg·kg−1·d−1).
Table 4. Dust heavy metal exposure doses at bus stations in Tianshui City (mg·kg−1·d−1).
GroupItemCrMnNiCuZnAsPbVCoTotal
ChildrenADDlng4.65 × 10−53.22 × 10−31.45 × 10−53.16 × 10−43.25 × 10−39.83 × 10−61.01 × 10−33.67 × 10−44.87 × 10−58.29 × 10−3
ADDlnh1.18 × 10−98.154 × 10−83.66 × 10−108.011 × 10−98.23 × 10−82.49 × 10−102.56 × 10−89.28 × 10−91.23 × 10−92.09 × 10−7
ADDdermal2.14 × 10−71.48 × 10−56.65 × 10−81.457 × 10−61.50 × 10−54.52 × 10−74.65 × 10−61.69 × 10−62.24 × 10−73.85 × 10−5
AdultsADDlng3.87 × 10−56.79 × 10−41.20 × 10−56.594 × 10−56.78 × 10−48.18 × 10−62.11 × 10−47.63 × 10−51.01 × 10−51.77 × 10−3
ADDlnh2.45 × 10−94.24 × 10−87.61 × 10−104.170 × 10−94.28 × 10−85.18 × 10−101.33 × 10−84.83 × 10−96.41 × 10−101.12 × 10−7
ADDdermal5.53 × 10−79.60 × 10−61.72 × 10−79.43 × 10−73.32 × 10−61.17 × 10−71.03 × 10−63.74 × 10−74.97 × 10−81.61 × 10−5
LADDlnh3.62 × 10−91.12 × 10−97.66 × 10−106.04 × 10−8
Table 5. Health risks of heavy metals at bus stations in Tianshui City.
Table 5. Health risks of heavy metals at bus stations in Tianshui City.
Research ObjectItemCrMnNiCuZnAsPbVCoHI
ChildrenHQlng1.627.00 × 10−27.02 × 10−47.88 × 10−31.08 × 10−23.26 × 10−22.87 × 10−15.24 × 10−22.44 × 10−32.09
HQlnh3.92 × 10−75.82 × 10−31.83 × 10−82.00 × 10−72.74 × 10−78.29 × 10−77.31 × 10−61.33 × 10−67.70 × 10−85.83 × 10−3
HQderm3.56 × 10−38.23 × 10−31.23 × 10−51.21 × 10−42.49 × 10−43.68 × 10−31.43 × 10−22.41 × 10−44.43 × 10−27.47 × 10−3
HI1.638.41 × 10−27.14 × 10−48.00 × 10−31.11 × 10−23.63 × 10−23.02 × 10−15.27 × 10−24.67 × 10−22.17
AdultsHQlng1.351.46 × 10−25.85 × 10−41.64 × 10−32.26 × 10−32.72 × 10−25.98 × 10−21.09 × 10−25.07 × 10−41.470
HQlnh8.16 × 10−73.03 × 10−33.81 × 10−81.04 × 10−71.43 × 10−71.73 × 10−63.80 × 10−66.90 × 10−74.01 × 10−83.04 × 10−3
HQderm9.22 × 10−35.33 × 10−33.19 × 10−57.86 × 10−55.54 × 10−59.51 × 10−43.18 × 10−35.35 × 10−59.83 × 10−42.87 × 10−2
HI1.362.29 × 10−26.16 × 10−41.72 × 10−32.31 × 10−52.81 × 10−26.30 × 10−21.10 × 10−21.03 × 10−21.50
Risk2.28 × 10−89.47 × 10−101.16 × 10−83.53 × 10−8
Table 6. Explanatory power (q value) of the studied factors for eight heavy metals in soil.
Table 6. Explanatory power (q value) of the studied factors for eight heavy metals in soil.
FactorsX1X2X3X4X5X6X7X8X9X10X11X12
Cr0.3882 *0.1037 *0.0407 *0.097 *0.0603 *0.0505 *0.00050.0118 *0.0616 *0.4129 *0.0179 *0.0947 *
Mn0.2465 *0.1579 *0.01970.065 *0.0830.02250.2366 *0.07110.4515 *0.0949 *0.0654 *0.2919 *
Ni0.2746 *0.0786 *0.083 *0.0836 *0.0465 *0.0257 *0.00630.00270.0281 *0.6485 *0.0045 **0.0027
Cu0.1327 *0.165 *0.0152 *0.0286 *0.01290.01310.0728 *0.0473 *0.0546 *0.0433 *0.0200 *0.2490
Zn0.2024 *0.1053 *0.0146 *0.0251 *0.0467 *0.0414 *0.155 *0.0588 *0.16660.0188 *0.0841 *0.2850 *
As0.3751 *0.1663 *0.0178 *0.014 *0.0425 *0.0196 *0.2901 *0.0982 *0.42270.1864 *0.124 *0.5099 *
Pb0.1909 *0.1632 *0.0423 *0.0142 *0.0426 *0.0567 *0.1175 *0.0292 *0.2092 *0.1227 *0.0811 *0.2204 *
V0.0554 *0.03980.00850.00370.02670.02910.1211 *0.01030.4034 *0.1151 *0.1397 *0.2746 *
Co0.3941 *0.089 *0.0396 *0.0761 *0.04050.1038 *0.00060.00360.0429 *0.6323 *0.00340.0163
Note: the explanatory power (q value) is significant at p < 0.05 (**) or p < 0.01 (*).
Table 7. Principal component analysis of heavy metal contents in dust.
Table 7. Principal component analysis of heavy metal contents in dust.
ElementCrMnNiCuZnAsPbVCo
Factor 1−0.5130.412−0.5780.2450.7840.6860.778−0.012−0.673
Factor 20.6330.3910.7340.1980.450.5350.3660.2040.688
Factor 3−0.3510.706−0.067−0.516−0.185−0.095−0.1040.910.121
Cumulative contribution rate%32.99558.20578.08887.4991.80495.71197.88699.349100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, B.; Wang, H.; Wang, X.; Ma, C.; Zhou, J.; Dai, X. Health Risks and Source Analysis of Heavy Metal Pollution from Dust in Tianshui, China. Minerals 2021, 11, 502. https://0-doi-org.brum.beds.ac.uk/10.3390/min11050502

AMA Style

Tan B, Wang H, Wang X, Ma C, Zhou J, Dai X. Health Risks and Source Analysis of Heavy Metal Pollution from Dust in Tianshui, China. Minerals. 2021; 11(5):502. https://0-doi-org.brum.beds.ac.uk/10.3390/min11050502

Chicago/Turabian Style

Tan, Bo, Hongwei Wang, Xinmin Wang, Chen Ma, Jing Zhou, and Xinyan Dai. 2021. "Health Risks and Source Analysis of Heavy Metal Pollution from Dust in Tianshui, China" Minerals 11, no. 5: 502. https://0-doi-org.brum.beds.ac.uk/10.3390/min11050502

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