Next Article in Journal
Prediction of the Impact of Land Use and Soil Type on Concentrations of Heavy Metals and Phthalates in Soil Based on Model Simulation
Next Article in Special Issue
Spatial and Temporal Variations of Heavy Metals’ Bioavailability in Soils Regulated by a Combined Material of Calcium Sulfate and Ferric Oxide
Previous Article in Journal
Arsenic Release from Various Sizes of Excavated Crushed Rock Mixed with Soil at Various Ratios under Three Mass Water Contents for Cropland
Previous Article in Special Issue
Isolation and Identification of Mercury-Tolerant Bacteria LBA119 from Molybdenum-Lead Mining Soils and Their Removal of Hg2+
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Source Apportionment of Potentially Toxic Elements in Baoshan Soils Employing Combined Receptor Models

1
Yunnan Agricultural University, Kunming 650201, China
2
Yunnan Laboratory of Improvement of Soil Fertility and Pollution Remediation, Kunming 650201, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxics 2023, 11(3), 268; https://doi.org/10.3390/toxics11030268
Submission received: 17 January 2023 / Revised: 9 March 2023 / Accepted: 10 March 2023 / Published: 14 March 2023

Abstract

:
Arable soils are crucial for national development and food security; therefore, contamination of agricultural soils from potentially toxic elements (PTEs) is a global concern. In this study, we collected 152 soil samples for evaluation. Considering the contamination factors and using the cumulative index and geostatistical methods, we investigated the contamination levels of PTEs in Baoshan City, China. Using principal component analysis, absolute principal component score-multivariate linear regression, positive matrix factorization, and UNMIX, we analyzed the sources and quantitatively estimated their contributions. The average Cd, As, Pb, Cu, and Zn concentrations were 0.28, 31.42, 47.59, 100.46, and 12.36 mg/kg, respectively. The Cd, Cu, and Zn concentrations exceeded the corresponding background values for Yunnan Province. The combined receptor models showed that natural and agricultural sources contributed primarily to Cd and Cu and As and Pb inputs, accounting for 35.23 and 7.67% pollution, respectively. Industrial and traffic sources contributed primarily to Pb and Zn inputs (47.12%). Anthropogenic activities and natural causes accounted for 64.76 and 35.23% of soil pollution, respectively. Industrial and traffic sources contributed 47.12% to pollution from anthropogenic activities. Accordingly, the control of industrial PTE pollution emissions should be strengthened, and awareness should be raised to protect arable land around roads.

1. Introduction

The expansion of cities and rapid development of modern agriculture have led to increased environmental pollution in cultivated soils. Some of the contaminants of greatest concern are potentially toxic elements (PTEs) [1,2], which not only destroy the quality of cultivated land, but also they indirectly or directly cause damage to human health [3,4]. Therefore, remediation of arable soil PTE contamination has become one of the most urgent problems in environmental science. However, PTEs in soils do not originate from human activities such as mining, metal smelting, pesticide/fertilizer usage, or automobilism [5,6,7]. Furthermore, several researchers have shown that the mineralization of soil parent material is the main cause of PTE enrichment [8,9].
Yunnan province has a complex geological structure and rich mineral resources. According to a 2020 report on the quality of arable land in the Baoshan region, medium-grade cultivation was performed on 263,500 ha of low-grade cultivated land, which accounted for 79.31% of the cultivated land area [10]. Currently, soil conditions are poor and heavy metal pollution is becoming an increasingly significant problem in the area [11,12,13,14,15]. Moreover, the PTE background value in Baoshan is significantly higher than the national average, making PTE pollution a serious hazard in the region [16,17].
Zhang et al. have found severe PTE enrichment in soils around the Baoshan mining area and have reported the soil to be at risk [18]. Other factors affecting Baoshan soil PTE enrichment include geogenesis and carbonate and basalt parent material differentiation. The Baoshan agricultural environmental protection monitoring station that evaluates the quality of farmland soil indicated that the arable soil in the study area is contaminated to varying degrees by different sources, the most prominent ones being Cd, Pb, As, Cu, and Zn PTEs [19,20,21]. Importantly, the quality and safety of arable soil determines the quality and safety of agricultural products. Cicchella et al. has recommended paying attention to the nature and sources of contaminants [22]. Still, previous studies have only evaluated the risk of soil contamination in Baoshan and have not distinguished the sources of contamination [23]. Thus, we must understand the sources of PTE pollution in arable land to protect it [24,25].
The methods of pollution source analysis are divided mainly into source identification and source apportionment. Scholars are increasingly using receptor models for source analyses, with the most widely used ones being principal component analysis (PCA), absolute principal component score–multivariate linear regression (APCS-MLR), UNMIX, and positive matrix factorization (PMF). Li et al. used PCA to determine potentially toxic elements in Gansu cropland soils; they found three sources of contamination and evaluated soil data [26]. The same group also used APCS-MLR to analyze the PTE sources in urban farmland. They combined the total hazard index (THI) and total carcinogenic risk (TCR) and found that APCS-MLR analyzed both the pollution sources and pollution status [27]. This model only required the composition of the receptor emission source and did not require the accurate source component spectrum data. Still, the estimated source component spectrum and contribution values were often observed. Chen et al. combined UNMIX with PMF for source analysis of heavy metal data in a suburban area of Kaifeng and found good agreement between the two models for source assignment [28]. The UNMIX model overcomes the drawback of negative source contribution, which does not require prior knowledge of the number of sources and identifies these sources via different source identifiers to infer them. However, it cannot determine the source components, contribution rates, or contribution values. Liu et al. used PMF to determine the pollution sources in a farmland in Wenzhou, Zhejiang Province, and combined it with multivariate statistical methods, which demonstrated the feasibility of PMF models for source analysis [29]. Lv used APCS/MLR and PMF, and both provided three identical sources and were consistent with the mapping of parent material distribution [30]. PMF avoids negative values in the results, thus providing interpretability and clear physical meaning to the obtained source component spectra and source contributions. In addition, PMF does not require the measurement of source component spectra and uses error estimates for each individual data point to deal with missing and imprecise data [31].
Therefore, the main objectives of this study were to determine the concentrations of five soil PTE elements (Cd, Pb, As, Cu, and Zn) and to investigate the sources of PTEs in soil samples from Baoshan cropland. This would allow us to better understand the specific sources of PTEs contamination in Baoshan cropland and to provide a reliable basis for improving pollution in the area. Specifically, our aims were to: (1) investigate the pollution status of five soil PTEs in Baoshan soil; (2) map the distribution of these soil PTEs using GIS; (3) identify the sources of contamination using PCA, APCS-MLR, PMF, and UNMIX models; and (4) compare the contributions of the four receptor models.

2. Materials and Methods

2.1. Study Area

Baoshan is located in the western part of the Yunnan Province of China (98°25′–100°02′ E, 24°08′–25°51′ N). The entire study area has a subtropical plateau climate, with the local annual temperature difference being slight but the daily temperature difference being significant. The average annual temperature is 15.5 °C, the average temperature of the coldest month is 8.2 °C, and the average temperature of the hottest month is 21 °C. The highest precipitation occurs over the western and southern areas. The terrain is high in the northwest and low in the southeast. The city is bordered by Dali, Lincang, the Salween River, and the Dehong Dai and Jingpo Autonomous Prefecture. The city spans three major river systems, namely, the Lancang, Salween, and Longchuan. The region is rich in terms of land, forest, and mineral resources and natural gas reserves; moreover, it is a hub for hydroelectric power generation and has several tourist attraction spots. The types of cultivated land are mainly dryland and a small number of paddy fields. The cultivated land in the study area is shown in Figure 1.
Many experts and scholars investigated the origin of the parent rock mass in Baoshan City and found severe PTE enrichment in soils around the Baoshan mining area. [17,32] The geology of Baoshan City mainly includes carbonate rocks and basalt parent materials; the geological composition of the full study area is shown in Figure 2. In addition, the mining of the rich mineral resources in Baoshan has caused contamination and parent material enrichment of the soil nearby. Moreover, the recent economic development of Baoshan City has driven the expansion of non-ferrous metal smelting, real estate development, animal husbandry, old and outdated production facility construction, and desulfurization, resulting in contaminated farmland.

2.2. Sample Collection and Preparation and Quality Control

We randomly collected 152 topsoil samples (at 0–20 cm depth) from agricultural land in Baoshan in August 2022, using a five-point sampling method with plum blossom. Moreover, 1 kg of soil samples was collected by mixing and dividing the samples into four parts. The sampling points were located by GPS. All the soil samples were stored in plastic bags and transported to the laboratory for air drying and removing plant roots, residues, and visible invaders. Samples were passed through a nylon sieve with a 2 mm pore size to remove soil samples greater than 2 mm. Then, the sieved samples were ground and passed through a nylon sieve with a 0.149 mm pore size (100 mesh), mixed, and prepared for analyses of the soil PTEs, Cd, Pb, Cu, Zn, and As. Table 1 shows the specific methods used for the determination of potentially toxic elements. A flame atomic absorption spectrometer and graphite furnace atomic absorption spectrometer (AA6880, Shimadzu, Kyoto, Japan) and atomic fluorescence spectrophotometer (AFS-230E, Haikou Instruments, Beijing, China) were used for element detection. National standard soil samples (GSS-25) were used for quality control. The test samples were analyzed three times each, and their relative standard deviation was ≤5%. Blank samples were tested for the determination of each heavy metal. Heavy metal standard solutions were used for each batch of sample reagents, and the element recovery rate was 90–110%.

2.3. Pollution Assessment

2.3.1. Pollution Factors

The pollution level of the study area was evaluated by pollution factor evaluation [33,34,35].
PF = C i B i
where Ci is the soil PTE concentration (mg/kg), and Bi indicates the corresponding background values of soil PTEs in Yunnan Province (mg/kg) (Cd 0.22, Pb 40.6, Cu 38.38, Zn 89.7, As 18.4). The PF values’ environmental indices could be divided into four classes to interpret the pollution levels of soil PTEs (Table 2) [36].

2.3.2. Geoaccumulation Index

The geoaccumulation index (Igeo) was developed by Muller (1969) to assess the level of heavy metal and metalloid elements in the sediment by comparing the status of the current concentration with the pre-industrial level. This method assesses heavy metal pollution related to anthropogenic activity, as well as the different rock geology regions, and the natural formation of heavy metal pollution [37].
I geo = log 2 C n 1.5 × B n
where Cn is the soil PTE concentration (mg/kg), and Bn is the corresponding background value of soil PTEs in Yunnan Province (mg/kg). The obtained Igeo values were classified into seven groups based on Categories (Table 3).

2.4. Receptor Models

Based on the receptor model, the contributions of various soil pollution sources were analyzed quantitatively by mathematical method according to the source and concentration of the receptor. Employing PMF, UNMIX, and APCS-MLR, several factors were extracted and identified according to source types, and their contributions were estimated.

2.4.1. APCS-MLR

Principal component analysis is a multivariate statistical analysis method that selects a small number of important variables by linear transformation. Using the linear combination of the original variables after standardization to form the principal component, PCA model can be transformed into several complementary and related comprehensive indexes. The method is used widely in, e.g., demographics, quantitative geography, molecular dynamics simulation, mathematical modeling, and mathematical analysis. The APCS-MLR model was proposed by Thurston [38]. After the data were standardized, the principal component factor was transformed into an absolute principal component score (APCS) by factor analysis. Subsequently, multiple linear regression analyses were performed to determine the content of each selected heavy metal, and the contribution rate of each factor to the pollution source was further calculated.
(1)
Standardizing raw data
Z i j = C i j C i ¯ σ i
where Zij is the standardized factor score, and Cij is the soil PTEi concentration (mg/kg). C i ¯ is the average soil PTEi concentration (mg/kg); σ i is the standard deviation of soil PTEi (mg/kg).
(2)
Introducing a factor with a concentration of 0
Z i 0 = 0 C i ¯ σ i
where APCS for each heavy metal element is obtained using ZijZi0. Using the obtained APCS, multiple linear regression analysis is conducted to obtain the regression coefficient:
C i = b i 0 + p = 1 p b p i × A P C S p
where Bi0 is the constant obtained from multivariate linear regression; bpi is the regression coefficient of source p to the soil PTE; and bpi × APCSp is the source contribution to Ci.

2.4.2. PMF

Positive matrix factorization (PMF) analysis, proposed by Paatero et al. [39], is a source analysis method used widely by the United States Environmental Protection Agency (USEPA). The method is often used in sediment, atmospheric, and soil contamination source analyses. The PMF model is least squares through multiple iterations to minimize the objective function Q for obtaining the optimal factor matrix and source profile.
Q = i = 1 n j = 1 m x i j k = 1 p G i k F k j u i j 2 Q = i = 1 n j = 1 m e i j u i j 2
where xij is the content of the j heavy metal element in the i sample; Gik is the contribution of source k to the i sample; Fkj is the content of the j heavy metal element in the k source; uij is the measured uncertainty (mg/kg); and eij is the model uncertainty. The sample chemical type uncertainty file is calculated as follows:
u i j =   5 6 × M L D , x i j < M L D x i j × E r r o r   F r a c t i o n 2 + M L D × 0.5 2   ,   x i j > M L D
where Cij is the concentration of the j sample chemical type of the i sample, and MDL is the species-specific method detection limit. The error fraction is a percentage of measurement uncertainty.

2.4.3. UNMIX Model

The UNMIX is a receptor model that is based on the pollutant concentration. The major sources and contributions of hand-held pollutants can be obtained directly through extremes. The UNMIX6.0 software be identified to Klstrong in data quality that signal-to-noise (S/n) was more than 2, and the fitting concentration was more than 0.8. Simple software operation, the results do not need their own analyses but cannot evaluate the source component, contribution rate, or contribution value. We used Equation (8), recommended by the USEPA, for performing source apportionment for potentially toxic elements.
C i j = l = 1 p k = 1 p U i k D k l V l j + ε i j
where U, D, and V are the np diagonal matrix, p × p diagonal matrix, and p × m matrix, respectively. εij is the error term that contains the variability of Cij and excludes the first major component p.

2.5. Data Treatment with Computer Software

The data was pre-processed using Excel 2010 software. The SPSS® Statistics 23.0 (Armonk, NY, USA) software has performed normal distribution test and statistical analysis of data of potentially toxic elements in soil. The figures were drawn by Origin 2021. The geostatistical analyst tools data exploration tool set analysis was performed by ArcGIS 10.6. The sampling point distribution map and heavy metal pollution spatial distribution feature map were drawn by ArcGIS 10.6. The analysis of PMF was conducted with PMF 5.0.

3. Results and Discussion

3.1. Pollution Characteristics of Soil PTEs

3.1.1. Soil PTE Concentration

Descriptive statistics for Cd, As, Pb, Cu, and Zn in soil PTEs are shown in Table 4. Compared with the background values for Yunnan Province [40], the average concentrations of PTEs Cd (0.28 mg/kg), Cu (47.59 mg/kg), and Zn (100.46 mg/kg) in soils were all higher, and the median values were 1.25, 1.24, and 1.12 times higher, respectively, than the corresponding regional background values. The average concentrations of soil PTEs did not exceed the national soil pollution risk screening values (GB15618-2018), indicating an overall low degree of pollution in the cultivated land of the study area. The box plot in Figure 3 shows the relationship between single soil sample and background values, as well as the risk screening value in Yunnan Province. The Cd content in 25–75% of samples was higher than the background and risk screening values. The contents of the potentially toxic elements As, Pb, Cu, and Zn in 25–75% of the samples were higher than the background values, and the contents of Cu and Zn in 54.6–48.6% of samples were higher than the background values.
In 35.53, 21.05, 54.61, 48.68, and 20.39% of the soil samples, the Cd, Pb, Cu, Zn, and As contents, respectively, exceeded the background values. The Cd, Pb, Cu, and As contents exceeded the risk screening values in 21.34, 2.63, 5.26, and 5.26% of the samples, respectively. These findings indicate the contamination of arable land in the study area and thus emphasize the need for further analysis of the contamination status.
The coefficient of variation (CV) reflects the degree of soil disturbance caused by anthropogenic activities. Table 4 shows the high variability of species elements, with variabilities of 0% < CV ≤ 12%, 12% < CV ≤ 31%, and >31% indicating low, moderate, and high variabilities, respectively [41]. The variations of Cd and As were 91.3 and 81.5%, respectively, indicating high degrees of variation and disturbance by anthropogenic activities. These activities include smelting and sugar production and plant fuel burning and sewage discharge causing Cd and As soil pollution. As Baoshan is located in the Cu–Pb–Zn belt, mining of these metals causes severe Cd, As, Cu, Zn, and Pb pollution of the surrounding farmland [42]. However, agricultural activities also lead to increasing Cd and As contents in farmland [43].
Relatively low variations were observed for Zn, Pb, and Cu. Both carbonate and basalt parent materials are developed in the Baoshan area, and soils with basalt parent materials such as granite are rich in Cd, Cu, and Zn. Soils derived from carbonate rocks are rich in Cd, Pb, Cu, and Zn, and their natural genesis results in high Zn and Cu contents and low coefficients of variation [44,45]. However, the Zn and Cu contents of the samples (48.6 and 54.6%, respectively) were higher than the background values. Duan et al. [46] found that improper agricultural practices could lead to an excess of Cu and Zn in soils, and mining could cause Zn, Cu, and Pb pollution. The pollution levels in the overall study area were low, but the soil PTE coefficient of variation was high, which could be owing to the fact that soil pollution is complex and can be caused by a combination of natural and anthropogenic factors. Accordingly, a policy for preventing further risk and a control strategy should be adopted, implying that the source of regional pollution must be identified to facilitate timely prevention and control.

3.1.2. Assessment of Soil PTE Pollution

To further advance our understanding of the extent of soil contamination, the PTEs in the study area were assessed using Igeo and PF, as shown in Figure 4. The results showed that the mean Igeo median and mean values of PTEs in five soils in the study area were all < 0, indicating non-pollution (Figure 4a); however, the soil PTE, Cd (10.53%), Pb (3.29%), Cu (23.68%), Zn (23.68%), and As (7.24%) samples showed weak to moderate pollution. A small proportion of soil was moderately polluted by PTEs (Cd, 11.8%; Cu, 2.6%). Similarly, as shown in Figure 4b, the average PF values of five PTEs in soils in the study area were all less than 3, indicating that the overall pollution in the study area was moderate. The average PF value of Cd (1.25) was the highest, followed by those of Cu (1.24), Zn (1.12), Pb (0.77), and As (0.67). The PF median values of <1 for Cd and Zn indicate moderate pollution, but the mean values of >1 for Cd and Zn indicate low pollution implying that most cultivated areas were slightly polluted, whereas other areas were weakly to moderately contaminated, with PTE contributions of 23.68, 21.05, 51.97, 48.68, and 20.39% for Cd, Pb, Cu, Zn, and As, respectively. A small part of the soil showed moderately contaminated, and the values for Cd (11.84%) and Cu (2.63%) indicated severe pollution from these elements in the farmland soil.
The relationship between the two pollution indices and soil PTEs is shown in Figure 4c. As the values of the five types of soil PTEs were less than 0, the smaller the PTE, the greater would be its proportion. Therefore, the largest contribution of soil PTEs to Igeo was of Zn and Cd, followed by that of Cu, As, and Pb. The greatest contribution of PTEs to PF was of Cd and Cu, followed by that of Zn, Pb, and As. As the results of the various pollution assessment methods differed, neither the degree of soil PTE pollution nor the pollution sources could not be judged directly. Therefore, assessing soil PTE pollution and pollution concentrations must be combined with spatial analyses to gradually explore the source of PTEs for related risk assessment.

3.2. Spatial Distribution of Soil PTEs

Analyses of the distribution of soil PTEs could provide information on the extent of soil pollution in the study area and improve our understanding of the sources of soil PTEs [47]. The data were interpolated using ArcGIS with the Kriging method, as shown in Figure 5. The distribution of Zn and Cu in soil was partly similar, with the high concentration area located near the northwest of China, where numerous coal mines are situated. Tu et al. [48] found that the development of the coal industry led to Zn and Cu pollution in the surrounding farmland. High concentrations of Cd, Pb, and Cu were found in the northeast of China, with the high Pb values being closer to the South, which are the main concentration areas of cities and towns. Lai et al. [49] indicated that automobile exhaust gas was the main source of Pb pollution. The high Cd and Cu value area is near the North. The investigation showed tungsten and copper deposits in the area, which are mostly located in the high mountain area. As and Zn were mainly concentrated in the southeast, with the spatial distribution characteristics of potentially toxic elements in soils being similar. Hu et al. [50] showed that the use of herbicides, chemical fertilizers, and animal manure could cause As and Zn pollution in soils. Pollution from Cu is distributed widely, and the other four soils (Cd, Zn, Pb, and As) exhibit concentrated distributions of PTEs. Therefore, step-by-step exploration of the sources of PTE contamination was necessary to analyze the source of soil PTE pollution based on a combination of receptor models and geostatistics.

3.3. Source Apportionment of Soil PTEs

3.3.1. Multivariate Statistical Analysis

Based on the correlation between PTEs, we used Pearson correlation analysis to determine whether the sources of soil PTEs were consistent [51]. A significant positive correlation between the elements indicates a similar source between the elements; a significant negative correlation between the elements indicates a difference in source between the elements [52]. The results of the correlation coefficient analysis are shown in Figure 6. Positive correlations were found between Pb–Cd, Cd–Cu, Cu–Pb, As–Pb, As–Cu, and As–Pb, and negative correlations were found between Pb–Zn and As–Zn. The Cu–Pb correlation coefficient was 0.589 (p < 0.01), indicating that Cu–Pb comes from the same pollution source, possibly industrial or traffic pollution. The following observations were made: (1) Pb is correlated with As, Cu, and Cd, which may be owing to the high proportion of Pb in different pollution sources; (2) the sources of Cd, Pb, Cu, and Zn could be similar; (3) Zn–Pb and As–Zn were correlated negatively with each other, indicating that they could have the same source but in opposite proportions. The source of each element and the correlation between them could be verified further with PCA.

3.3.2. APCS-MLR Model

We used SPSS Statistics 23.0 to analyze the soil heavy metal concentration. The extraction characteristics of the first two principal components were more than 1; the Kaiser–Meyer–Olkin measure of sampling adequacy and Bartlett values were 0.53, and the Bartlett spherical test p value was 0.00 (p < 0.05). These results showed a strong correlation among the potentially toxic elements, which could be analyzed using PCA. The cumulative contribution of variance was 58.25%, which could explain most of the information on the soil PTEs. Potentially toxic elements with a higher factor load under the same principal component have the same origin [53,54]. The first principal component (PC1) was Pb, Zn, and As, the variance contribution was 30.35%, and Zn was the main factor. Although the average values of Pb and As were lower than the background values in Yunnan Province, the samples of Pb (21%) and As (20%) sites exceeded the background values. As shown in Figure 7, the areas with high Pb concentrations were mainly located close to the Longyang District towns and cities, with the As concentrations deriving mainly from fertilizers, fossil fuels, and anthropogenic activities [55]. Investigating the mineral distribution in Baoshan showed the presence of numerous iron and coal mines in the region, the mining of which results in Pb, Zn, As, and Cu pollution [56]. Agricultural activities cause Zn and Cu pollution because of unreasonable fertilizer use, excessive application of conditioning agents, and agricultural machinery activities [57,58]. Therefore, Factor 1 could be considered anthropogenic pollution.
The second component (PC2), with a higher driving, included Cd and Cu, of which the variance contribution was 27.91%, with Cd being the main factor. Significant correlation was found between Cd and Cu, indicating that both pollutants were likely derived from the same source. Moreover, Baoshan is located on carbonate and basalt parent material, and carbonate parent material can weather Cd, Cu, Pb, Zn, and As. The basalt matrix could be divided into Cd, Cu, and Zn [59], and Cu and Cd could be differentiated from the mixed carbonate matrix [60]. Yunnan Province is characterized by a complex geological structure and rich metal deposits. The main distribution area of Cd and Cu is alpine cultivated land. Therefore, Factor 2 could be artificially induced as a natural cause.
The contribution rates of soil PTEs used to calculate the APCS score and the APCS-MLR receptor model of each soil PTE were obtained by using regression analysis (Figure 7). The main contributing factors were As, Cu, and Pb [61], which could be considered as agricultural pollution, with the main contributing factors being Cd, Zn, and Cu [62]. These results suggested that Factor 2 was a natural cause, and Factor 3 was the main contributing factor of Pb, Cd, and Zn, i.e., industrial and traffic pollution. Our results were consistent with those of [63].

3.3.3. PMF Model

Employing the PMF model, we determined that the signal-to-noise ratio (S/N) of all chemicals was greater than the combined Q value, the resulting source profile, and the scale residuals. The model was run 20 times, resulting in the identification of three sources, as shown in Figure 8. The main contributing factors to Factor 1 were Cd, As, and Cu [64], whereas those to Factor 2 were As and Pb [65]. The results showed that As and Pb pollution could be caused by the application of chemical fertilizers and pesticides in agriculture, i.e., As agricultural pollution. These results were consistent with the analysis results of Han et al. using the PMF model [66]. The main contributing factors of Factor 3 were Zn, Cu, and Pb, which could be considered industrial and traffic pollution. The results were consistent with those of Wang et al. [65].

3.3.4. UNMIX Model

We used UNMIX 6.0 software to analyze the sample concentration data. The results showed that Min-Rsq was 0.82, i.e., higher than the threshold of 0.8, and S/N was 2.21, greater than the system requirement maximum threshold of 2. Therefore, the approach was considered successful, and three sources could be detected, as shown in Figure 8. Factor 1 was mainly contributed to by Cu, As, and Pb [67], with the assumption being that Factor 1 is agricultural pollution. The main contributing factors to Factor 2 were Cd, Zn, and Cu [68], with Factor 2 considered natural causes. The main contributing factors to Factor 3 were Pb, Cd, and Zn, which could be considered industrial and traffic pollution. These results were consistent with those of Luo et al. [69].

3.4. Model Evaluation

The results of the three models were consistent. As shown in Figure 9a, the main sources of pollution derived from industrial activities and traffic, followed by natural causes and agricultural activity. Clearly, the pollution originated from human activity. Combining the results from the three receptor models enabled us to compile the information more comprehensively. As shown in Figure 9b, the main source of pollution in the cultivated land was disturbance by humans, which severely affected the quality of the land. Combined with Figure 7, Figure 8 and Figure 9, the four models identified three common source categories other than PCA, namely, natural sources (35.23%), agricultural sources (17.67%), and industrial and transportation sources (47.12%). The difference between the model results was that UNMIX yielded a low resolution of industrial and transportation sources and a high resolution of natural and agricultural sources, whereas APCS-MLR and PMF and resolution of the same contribution rate. In addition to UNMIX, the APCS-MLR and PMF models revealed that the agricultural source was dominated mainly by As, whereas the natural source was dominated primarily by Cd. We employed the UNMIX model to analyze the Cu with the highest contribution rate to As, with As being considered the main factor in agricultural pollution. The results of the other three models accounted for a substantial proportion of Zn and Pb in the analysis of industrial and traffic pollution. The potential variables could not be estimated by PCA, and the negative values limited the results of the analysis. Although APCS-MLR evolved from PCA, negative numbers remained in the calculation process. In contrast, UNMIX is a simple model that does not require setting the number of pollution sources and well explains the species concentrations [70]. Employing PMF, we could estimate the error at each site. The reasonable treatment of missing and imprecise data is used widely in soil and air pollution analyses; however, the results are affected by uncertainties of the data and model structure [29]. Three source apportionment models could be used to quantitatively identify pollution, with the contribution rate being the only difference [71,72]. Accordingly, an error that could be caused by employing a single model could be counteracted by the comparative analysis of several receptor models [73,74].

4. Conclusions

The pollution characteristics and source contributions of PTEs in cultivated soils in Baoshan were investigated, and three receptor models were used to analyze and compare the pollution sources in the study area. The concentrations of Cd, Cu, and Zn in soil were higher than the corresponding background values in Yunnan province. According to the Igeo and PF values, the cultivated soil was polluted to different degrees. In particular, pollution from Cd and Cu was severe. The analytical results of the APCS-MLR, PMF, and UNMIX models indicated that As and Pb pollution was caused by agricultural activities, with As pollution deriving mostly from the improper use of agricultural chemicals, together with Cu and Zn pollution. Moreover, the use of agricultural machinery in the region led to Pb pollution. The natural source is mainly caused by the elements Cd and Cu. As the study site is a carbonate-weathering area, Cd is the most significant controlling factor. Automobile exhaust emissions were the main source of Pb and Zn pollution. The emissions of Cu and Cd were the main source of Zn pollution. Overall, 64.76% of the pollution in Baoshan derived from agriculture, of which 47.12% of the pollution derived from industrial transport. Accordingly, this local study should pay more attention to the pollution by industrial production, and more attention is required to mitigate pollution caused by industrial production and traffic emissions.

Author Contributions

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

Funding

This research was funded by the Studies on the transformation, translocation regularity, and control mechanisms of cadmium for high geological background and anthropogenic pollution soils in Yunnan Province, grant number U200220167 and Key Technologies of green production of agricultural products quality safety control, grant number 202002AE320005.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by the Studies on the transformation, translocation regularity, and control mechanisms of cadmium for high geological background and anthropogenic pollution soils in Yunnan Province (U200220167), Key Technologies of green production of agricultural products quality safety control (202002AE320005), and the Transport and accumulation of cadmium and gene expression in low-cumulative cadmium in soil-corn system (201801YB00003).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barsby, A.; McKinley, J.M.; Ofterdinger, U.; Young, M.; Cave, M.R.; Wragg, J. Bioaccessibility of trace elements in soils in Northern Ireland. Sci. Total Environ. 2012, 433, 398–417. [Google Scholar] [CrossRef] [Green Version]
  2. Ahado, S.K.; Nwaogu, C.; Sarkodie, V.Y.O.; Borůvka, L. Modeling and Assessing the Spatial and Vertical Distributions of Potentially Toxic Elements in Soil and How the Concentrations Differ. Toxics 2021, 9, 181. [Google Scholar] [CrossRef]
  3. Hou, D.; O’Connor, D.; Nathanail, P.; Tian, L.; Ma, Y. Integrated GIS and multivariate statistical analysis for regional scale assessment of heavy metal soil contamination: A critical review. Environ. Pollut. 2017, 231, 1188–1200. [Google Scholar] [CrossRef] [PubMed]
  4. Jin, Y.; O’Connor, D.; Sik Ok, Y.; Tsang, D.C.W.; Liu, A.; Hou, D. Assessment of sources of heavy metals in soil and dust at children’s playgrounds in Beijing using GIS and multivariate statistical analysis. Environ. Int. 2019, 124, 320–328. [Google Scholar] [CrossRef]
  5. Jadoon, S.; Muhammad, S.; Hilal, Z.; Ali, M.; Khan, S.; Khattak, N. Spatial distribution of potentially toxic elements in urban soils of Abbottabad city, (N Pakistan): Evaluation for potential risk. Microchem. J. 2020, 153, 104489. [Google Scholar] [CrossRef]
  6. Moghtaderi, T.; Shakeri, A.; Rodríguez-Seijo, A. Potentially Toxic Element Content in Arid Agricultural Soils in South Iran. Agronomy 2020, 10, 564. [Google Scholar] [CrossRef]
  7. Wang, X.; Wang, L.; Zhang, Q.; Liang, T.; Li, J.; Bruun, H.H.C.; Sabry, M.S.; Vasileios, A.; Nanthi, B.; Jörg, R. Integrated assessment of the impact of land use types on soil pollution by potentially toxic elements and the associated ecological and human health risk. Environ. Pollut. 2022, 299, 118911. [Google Scholar] [CrossRef] [PubMed]
  8. Zuzolo, D.; Cicchella, D.; Lima, A.; Guagliardi, I.; Cerino, P.; Pizzolante, A.; Thiombane, M.; De Vivo, B.; Albanese, S. Potentially toxic elements in soils of Campania region (Southern Italy): Combining raw and compositional data. J. Geochem. Explor. 2020, 213, 106524. [Google Scholar] [CrossRef]
  9. Nogueira, T.A.R.; Abreu-Junior, C.H.; Alleoni, L.R.F.; He, Z.; Soares, M.R.; Vieira, C.d.S.; Lessa, L.G.F.; Capra, G.F. Background concentrations and quality reference values for some potentially toxic elements in soils of São Paulo State, Brazil. J. Environ. Manag. 2018, 221, 10–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Available online: http://www.baoshan.gov.cn/info/egovinfo/1001/zfxxgkpt/zfxxgkptzn-content/01525525-X-/2021-1229001.htm (accessed on 12 December 2022.).
  11. Li, Q.; Hu, Q.; Zhang, C.; Jin, Z. Effects of Pb, Cd, Zn, and Cu on Soil Enzyme Activity and Soil Properties Related to Agricultural Land-Use Practices in Karst Area Contaminated by Pb-Zn Tailings. Pol. J. Environ. Stud. 2018, 27, 2623–2632. [Google Scholar] [CrossRef]
  12. Jaradat, Q.M.; Massadeh, A.M.; Momani, K.A.; Al Saleem, M.A. The Spatial Distribution of Pb, Cd, Zn, and Cu in Agricultural Roadside Soils. SoilSediment Contam. 2010, 19, 58–71. [Google Scholar] [CrossRef]
  13. Zhou, R.; Liu, X.; Luo, L.; Zhou, Y.; Wei, J.; Chen, A.; Tang, L.; Wu, H.; Deng, Y.; Zhang, F.; et al. Remediation of Cu, Pb, Zn and Cd-contaminated agricultural soil using a combined red mud and compost amendment. Int. Bioeterioration Biodegrad. 2017, 118, 73–81. [Google Scholar] [CrossRef]
  14. Liu, Y.; Liu, D.; Zhang, W.; Chen, X.; Zhao, Q.; Chen, X.; Zou, C. Health risk assessment of heavy metals (Zn, Cu, Cd, Pb, As and Cr) in wheat grain receiving repeated Zn fertilizers. Environ. Pollut. 2020, 257, 113581. [Google Scholar] [CrossRef]
  15. Liao, Z.; Chen, Y.; Ma, J.; Islam, M.S.; Weng, L.; Li, Y. Cd, Cu, and Zn Accumulations Caused by Long-Term Fertilization in Greenhouse Soils and Their Potential Risk Assessment. Int. J. Environ. Res. Public Health 2019, 16, 2805. [Google Scholar] [CrossRef] [Green Version]
  16. Zhan, J.; Li, X.; Christie, P.; Wu, L. A review of soil potentially toxic element contamination in typical karst regions in southwest China. Curr. Opin. Environ. Sci. Health 2021, 23, 100284. [Google Scholar] [CrossRef]
  17. Zhang, L.; Yang, Z.; Peng, M.; Cheng, X. Contamination Levels and the Ecological and Human Health Risks of Potentially Toxic Elements (PTEs) in Soil of Baoshan Area, Southwest China. Appl. Sci. 2022, 12, 1693. [Google Scholar] [CrossRef]
  18. Zhang, C.; Wang, Z.; Liu, L.; Liu, Y. Source Analysis of Soil Heavy Metals in Agriccultural Land the Ming Area Based on APCS-MLR Receptor Model and Geostatistical Method. Environ. Sci. 2022, 1–14. [Google Scholar] [CrossRef]
  19. Khorasanipour, M.; Karimabadi, F.; Espahbodi, M.; Ebrahimnejad, M. The effect of flotation desulfurization on the trace element geochemistry of Sarcheshmeh mine tailings, SE of Iran: Recycling and the environmental opportunities. Environ. Earth Sci. 2021, 80, 420. [Google Scholar] [CrossRef]
  20. Liu, Y.; Zhang, W.; Yang, W.; Bai, Z.; Zhao, X. Chemical Compositions of PM2.5 Emitted from Diesel Trucks and Construction Equipment. Aerosol Sci. Eng. 2018, 2, 51–60. [Google Scholar] [CrossRef] [Green Version]
  21. Rozanski, S.; Castejon, J.; McGahan, D. Child risk assessment of selected metal(loid)s from urban soils using in vitro UBM procedure. Ecol. Indic. 2021, 127, 107726. [Google Scholar] [CrossRef]
  22. Cicchella, D.; Zuzolo, D.; Albanese, S.; Fedele, L.; Di Tota, I.; Guagliardi, I.; Thiombane, M.; De Vivo, B.; Lima, A. Urban soil contamination in Salerno (Italy): Concentrations and patterns of major, minor, trace and ultra-trace elements in soils. J. Geochem. Explor. 2020, 213, 106519. [Google Scholar] [CrossRef]
  23. Li, J.; Wu, J.; Jiang, J.; Teng, Y.; He, L.; Song, L. Review on Source Apportionment of Soil Pollutants in Recent Ten Years. Soil Sci. 2018, 49, 232–242. [Google Scholar] [CrossRef]
  24. Wu, J.; Long, J.; Liu, L.; Li, J.; Liao, H.; Zhang, M.; Zhao, C.; Wu, Q. Risk Assessment and Source Identification of Toxic Metals in the Agricultural Soil around a Pb/Zn Mining and Smelting Area in Southwest China. Int. Environ. Res. Public Health 2018, 15, 1838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Negahban, S.; Mokarram, M. Potential Ecological Risk Assessment of Ni, Cu, Zn, Cd, and Pb in Roadside Soils. Earth Space Sci. 2021, 8, e2020EA001120. [Google Scholar] [CrossRef]
  26. Li, Y.; Gou, X.; Wang, G.; Zhang, Q.; Su, Q.; Xiao, G. Heavy metal contamination and source in arid agricultural soil in central Gansu Province, China. Environ. Sci. 2008, 20, 607–612. [Google Scholar] [CrossRef]
  27. Li, N.; Li, Y.; Wang, G.; Zhang, H.; Zhang, X.; Wen, J.; Cheng, X. The sources risk assessment combined with APCS/MLR model for potentially toxic elements in farmland of a first-tier city, China. Environ. Sci. Pollut. Res. 2022, 29, 50717–50726. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, Z.; Ding, Y.; Jiang, X.; Duan, H.; Ruan, X.; Li, Z.; Li, Y. Combination of UNMIX, PMF model and Pb-Zn-Cu isotopic compositions for quantitative source apportionment of heavy metals in suburban agricultural soils. Ecotoxicol. Environ. Saf. 2022, 234, 113369. [Google Scholar] [CrossRef] [PubMed]
  29. Liu, H.; Anwar, S.; Fang, L.; Chen, L.; Xu, W.; Xiao, L.; Zhong, B.; Liu, D. Source Apportionment of Agricultural Soil Heavy Metals Based on PMF Model and Multivariate Statistical Analysis. Environ. Forensics 2022, 1–9. [Google Scholar] [CrossRef]
  30. Lv, J. Multivariate receptor models and robust geostatistics to estimate source apportionment of heavy metals in soils. Environ. Pollut. 2018, 244, 72–83. [Google Scholar] [CrossRef]
  31. Wang, J.; Yang, J.; Chen, T. Source appointment of potentially toxic elements (PTEs) at an abandoned realgar mine: Combination of multivariate statistical analysis and three common receptor models. Chemosphere 2022, 307, 135923. [Google Scholar] [CrossRef]
  32. Zhang, L.; McKinley, J.; Cooper, M.; Peng, M.; Wang, Q.; Song, Y.; Cheng, H. A regional soil and river sediment geochemical study in Baoshan area, Yunnan province, southwest China. Geochem. Explor. 2020, 217, 106557. [Google Scholar] [CrossRef]
  33. Lei, M.; Li, K.; Guo, G.; Ju, T. Source-specific health risks apportionment of soil potential toxicity elements combining multiple receptor models with Monte Carlo simulation. Sci. Total Environ. 2022, 817, 152899. [Google Scholar] [CrossRef] [PubMed]
  34. Candeias, C.; Avila, P.; da Silva, E.; Teixeira, J. Integrated approach to assess the environmental impact of mining activities: Estimation of the spatial distribution of soil contamination (Panasqueira mining area, Central Portugal). Environ. Monit. Assess. 2015, 187, 135. [Google Scholar] [CrossRef]
  35. Turan, D.; Kocahakimoglu, C.; Kavcar, P.; Gaygısız, H.; Atatanir, L.; Turgut, C.; Sofuoglu, S. The use of olive tree (Olea europaea L.) leaves as a bioindicator for environmental pollution in the Province of Aydın, Turkey. Env. Sci. Pollut. Res. 2011, 18, 355–364. [Google Scholar] [CrossRef] [Green Version]
  36. Rehman, I.; Ishaq, M.; Ali, L.; Khan, S.; Ahmad, I.; Din, I.; Ullah, H. Enrichment, spatial distribution of potential ecological and human health risk assessment via toxic metals in soil and surface water ingestion in the vicinity of Sewakht mines, district Chitral, Northern Pakistan. Ecotoxicol. Environ. Saf. 2018, 154, 127–136. [Google Scholar] [CrossRef] [PubMed]
  37. Mukhopadhyay, S.; Chakraborty, S.; Bhadoria, P.B.S.; Li, B.; Weindorf, D. Assessment of heavy metal and soil organic carbon by portable X-ray fluorescence spectrometry and NixProTM sensor in landfill soils of India. Geoderma Reg. 2020, 20, e00249. [Google Scholar] [CrossRef]
  38. Thurston, G.; Spengler, J. A quantitative assessment of source contributions to inhalable particulate matter pollution in metropolitan Boston. Atmos. Environ. 1985, 19, 9–25. [Google Scholar] [CrossRef]
  39. Paatero, P.; Tapper, U. Positive matrix factorization: A non-negative factor model with optimal utilization of error estimates of data values. Environmetrics 1994, 5, 111–126. [Google Scholar] [CrossRef]
  40. CNEMC. Chinese Soil Element Background Value; Environmental Science Press: Beijing, China, 1990. [Google Scholar]
  41. Gholizadeh, A.; Borůvka, L.; Saberioon, M.M.; Kozák, J.; Vašát, R.; Němeček, K. Comparing Different Data Preprocessing Methods for Monitoring Soil Heavy Metals Based on Soil Spectral Features. Soil and Water Research. 2015, 10, 218–227. [Google Scholar] [CrossRef] [Green Version]
  42. Xing, R.; Wu, Z.; Du, G. Risk assessment and source analysis soil heavy metal pollution in Xuanzhou District, Anhui Province. East China Geol. 2022, 43, 336–344. [Google Scholar] [CrossRef]
  43. Gao, Y.; Li, F.; Mao, L.; Gu, B.; Peng, C.; Yang, Q.; Lu, L.; Chen, X.; Zhang, D.; Tao, H. Potential Loss of Toxic Elements from Slope Arable Soil Erosion into Watershed in Southwest China: Effect of Spatial Distribution and Land-Uses. Minerals 2021, 11, 1422. [Google Scholar] [CrossRef]
  44. Xiao, J.; Chen, W.; Wang, L.; Zhang, X.; Wen, Y.; Bostick, B.; Wen, Y.; He, X.; Zhang, L.; Zhuo, X.; et al. New strategy for exploring the accumulation of heavy metals in soils derived from different parent materials in the karst region of southwestern China. Geoderma 2022, 417, 115806. [Google Scholar] [CrossRef]
  45. Cabral Pinto, M.; Silva, M.; Ferreira da Silva, E.; Dinis, P.; Rocha, F. Transfer processes of potentially toxic elements (PTE) from rocks to soils and the origin of PTE in soils: A case study on the island of Santiago (Cape Verde). Geochem. Explor. 2017, 183, 140–151. [Google Scholar] [CrossRef]
  46. Duan, B.; Qiang, F. Comparison of the Potential Ecological and Human Health Risks of Heavy Metals from Sewage Sludge and Livestock Manure for Agricultural Use. Toxics 2021, 9, 145. [Google Scholar] [CrossRef] [PubMed]
  47. Francos, N.; Gholizadeh, A.; Ben Dora, E. Spatial distribution of lead (Pb) in soil: A case study in a contaminated area of the Czech Republic. Geomat. Nat. Hazards Risk 2022, 12, 610–620. [Google Scholar] [CrossRef]
  48. Tu, C.; Yang, K.; He, C.; Zhang, L.; Li, B.; Wei, Z.; Jian, X.; Yang, M. Sources and Risk Assessment of Heavy Metals in Small Watersheds in Typical Coal Mining Areas of Eastern Yunnan. East China Geol. 2022, pp. 1–16. Available online: http://kns.cnki.net/kcms/detail/11.1167.P.20220330.1029.002.html (accessed on 16 January 2023).
  49. Lai, S.; Cao, R.; Tan, G. Study on accumulation characteristics of heavy metal elements in topsoil of southern Longhai City, Fujian Province. East China Geol. 2021, 42, 29–36. [Google Scholar] [CrossRef]
  50. Hu, Y.; He, K.; Sun, Z.; Chen, G.; Cheng, H. Quantitative source apportionment of heavy metal(loid)s in the agricultural soils of an industrializing region and associated model uncertainty. Hazard. Mater. 2020, 391, 122244. [Google Scholar] [CrossRef]
  51. Gong, C.; Wang, S.; Wang, D.; Lu, H.; Dong, H.; Liu, J.; Yan, B.; Wang, L. Ecological and human health risk assessment of heavy metal(loid)s in agricultural soil in hotbed chives hometown of Tangchang, Southwest China. Sci. Rep. 2022, 12, 8563. [Google Scholar] [CrossRef]
  52. Singh, H.; Pandey, R.; Singh, S.K.; Shukla, D.N. Assessment of heavy metal contamination in the sediment of the River Ghaghara, a major tributary of the River Ganga in Northern India. Appl. Water Sci. 2017, 7, 4133–4149. [Google Scholar] [CrossRef]
  53. Jin, G.; Fang, W.; Shafi, M.; Wu, D.; Li, Y.; Zhong, B.; Ma, J.; Liu, D. Source apportionment of heavy metals in farmland soil with application of APCS-MLR model: A pilot study for restoration of farmland in Shaoxing City Zhejiang, China. Ecotoxicol. Environ. Saf. 2019, 184, 109495. [Google Scholar] [CrossRef]
  54. Samsudin, M.; Azid, A.; Khalit, S.; Saudi, A.; Zaudi, M. River water quality assessment using APCS-MLR and statistical process control in Johor River Basin, Malaysia. Adv. Appl. Sci. 2017, 4, 84–97. [Google Scholar] [CrossRef]
  55. Jiang, H.; Cai, L.; Wen, H.; Hu, G.; Chen, L.; Luo, J. An integrated approach to quantifying ecological and human health risks from different sources of soil heavy metals. Total Environ. 2020, 701, 134466. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, X.; Lei, M.; Zhang, S.; Zhang, D.; Guo, G.; Zhao, X. Apportionment and Spatial Pattern Analysis of Soil Heavy Metal Pollution Sources Related to Industries of Concern in a County in Southwestern China. Environ. Res. Public Health 2022, 19, 7421. [Google Scholar] [CrossRef]
  57. Li, C.; Zhang, C.; Yu, T.; Liu, X.; Xia, X.; Hou, Q.; Yang, Y.; Yang, Z.; Wang, L. Annual net input fluxes of cadmium in paddy soils in karst and non-karst areas of Guangxi, China. Geochem. Explor. 2022, 241, 107072. [Google Scholar] [CrossRef]
  58. Wu, J.; Chen, Y.; Ma, J.; Cao, J.; Jiang, Y. Sustainable Strategies for the Agricultural Development of Shaanxi Province Based on the Risk Assessment of Heavy Metal Pollution. Foods 2022, 11, 1409. [Google Scholar] [CrossRef]
  59. Wu, H.; Wang, J.; Guo, J.; Hu, X.; Bao, H.; Chen, J. Record of heavy metals in Huguangyan Maar Lake sediments: Response to anthropogenic atmospheric pollution in Southern China. Sci. Total Environ. 2022, 831, 154829. [Google Scholar] [CrossRef] [PubMed]
  60. Xiao, G.; Chen, J.; Bai, B.; Li, Y.; Zhu, N. Content Characteristics and Risk Assessment of Heavy Metals in Soil of Typical High Geological Background Areas, Yunnan Province. Geol. Explor. 2021, 57, 1077–1086. [Google Scholar]
  61. Zhang, C.; Zou, X.; Yang, H.; Liang, J.; Zhu, T. Bioaccumulation and Risk Assessment of Potentially Toxic Elements in Soil-Rice System in Karst Area, Southwest China. Front. Environ. Sci. 2022, 10, 866427. [Google Scholar] [CrossRef]
  62. Li, Y.; Xu, Z.; Ren, H.; Wang, D.; Wang, J.; Wu, Z.; Cai, P. Spatial Distribution and Source Apportionment of Heavy Metals in the Topsoil of Weifang City, East China. Front. Environ. Sci. 2022, 10, 893938. [Google Scholar] [CrossRef]
  63. Wang, J.; Wu, H.; Wei, W.; Xu, C.; Tan, X.; Wen, Y.; Lin, A. Health risk assessment of heavy metal(loid)s in the farmland of megalopolis in China by using APCS-MLR and PMF receptor models: Taking Huairou District of Beijing as an example. Total Environ. 2022, 835, 155313. [Google Scholar] [CrossRef]
  64. Chen, Z.; Xu, J.; Duan, R.; Lu, S.; Hou, Z.; Yang, F.; Peng, M.; Zong, Q.; Shi, Z.; Yu, L. Ecological Health Risk Assessment and Source Identification of Heavy Metals in Surface Soil Based on a High Geochemical Background: A Case Study in Southwest China. Toxics 2022, 10, 282. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, S.; Cai, L.; Wen, H.; Luo, J.; Wang, Q.; 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]
  66. Han, L.; Xu, X. Quantitative Evaluation of Human Health Risk of Heavy Meyals in Soil Based on Postitive Matrix Factorization Model and Geo-statistics. Environ. Sci. 2020, 41, 5114–5124. [Google Scholar] [CrossRef]
  67. Atafar, Z.; Mesdaghinia, A.; Nouri, J.; Homaee, M.; Yunesian, M.; Ahmadimoghaddam, M.; Mahvi, A. Effect of fertilizer application on soil heavy metal concentration. Environ Monit Assess. Environ. Monit. Assess. 2010, 160, 83. [Google Scholar] [CrossRef]
  68. Zhang, F.; Peng, M.; Wang, H.; Ma, H.; Xu, R.; Cheng, X.; Hou, Z.; Chen, Z.; Li, K.; Cheng, H. Ecological risk assessment of heavy metals at township scale in thehigh background of heavy metals, southwestern China. Environ. Sci. 2020, 41, 4197–4209. [Google Scholar] [CrossRef]
  69. Luo, X.; Wu, C.; Lin, Y.; Li, W.; Deng, M.; Tan, J.; Xue, S. Soil heavy metal pollution from Pb/Zn smelting regions in China and the remediation potential of biomineralization. Environ. Sci. 2022, 125, 662–677. [Google Scholar] [CrossRef] [PubMed]
  70. Li, Y.; Kuang, H.; Hu, C.; Ge, G. Source Apportionment of Heavy Metal Pollution in Agricultural Soils around the Poyang Lake Region Using UNMIX Model. Sustainability 2021, 13, 5272. [Google Scholar] [CrossRef]
  71. Zhang, S.; Wang, L.; Zhang, W.; Wang, L.; Shi, X.; Lu, X. Pollution assessment and source apportionment of trace metals in urban topsoil of Xi’ a City in Northwest China. Arch. Environ. Contam. Toxicol. 2019, 77, 575–586. [Google Scholar] [CrossRef]
  72. 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. Clean. Prod. 2019, 237, 117792. [Google Scholar] [CrossRef]
  73. Yang, B.; Zhou, L.; Xue, N.; Li, F.; Li, Y.; Vogt, R.; Cong, X.; Yan, Y.; Liu, B. Source apportionment of polycyclic aromatic hydrocarbons in soils of Huanghuai Plain, China: Comparison of three receptor models. Sci. Total Environ. 2013, 443, 31–39. [Google Scholar] [CrossRef]
  74. Jin, Z.; Lv, J. Integrated receptor models and multivariate geostatistical simulation for source apportionment of potentially toxic elements in soils. CATENA 2020, 194, 104638. [Google Scholar] [CrossRef]
Figure 1. Sampling locations of the study area.
Figure 1. Sampling locations of the study area.
Toxics 11 00268 g001
Figure 2. Geological maps of the study area.
Figure 2. Geological maps of the study area.
Toxics 11 00268 g002
Figure 3. Box plot of soil PTEs in the study area. The red line of the box plot of soil PTE content represents the background value, and the purple line represents the risk screening value.
Figure 3. Box plot of soil PTEs in the study area. The red line of the box plot of soil PTE content represents the background value, and the purple line represents the risk screening value.
Toxics 11 00268 g003
Figure 4. Violin with box plot of soil PTEs in the study area (a) Igeo, (b) PF, and (c) the relationship between soil PTEs and environmental indices.
Figure 4. Violin with box plot of soil PTEs in the study area (a) Igeo, (b) PF, and (c) the relationship between soil PTEs and environmental indices.
Toxics 11 00268 g004
Figure 5. Spatial distribution of soil PTEs in the study area.
Figure 5. Spatial distribution of soil PTEs in the study area.
Toxics 11 00268 g005aToxics 11 00268 g005b
Figure 6. Correlation analysis for soil PTEs. (* indicates significant correlation.)
Figure 6. Correlation analysis for soil PTEs. (* indicates significant correlation.)
Toxics 11 00268 g006
Figure 7. Source profiles and source concentrations of soil potentially toxic elements from PCA. ((a) is the source profiles and source concentrations of soil potentially toxic elements for factor 1; (b) is the source profiles and source concentrations of soil potentially toxic elements for factor 2).
Figure 7. Source profiles and source concentrations of soil potentially toxic elements from PCA. ((a) is the source profiles and source concentrations of soil potentially toxic elements for factor 1; (b) is the source profiles and source concentrations of soil potentially toxic elements for factor 2).
Toxics 11 00268 g007
Figure 8. Source contribution of each factor derived from PMF, UNMIX, and APCS/MLR.
Figure 8. Source contribution of each factor derived from PMF, UNMIX, and APCS/MLR.
Toxics 11 00268 g008
Figure 9. Estimated average source contribution (%).
Figure 9. Estimated average source contribution (%).
Toxics 11 00268 g009
Table 1. Detection methods.
Table 1. Detection methods.
Test ItemsDetection MethodGB NumberDetection Limits
Copper (Cu)
Zinc (Zn)
Furnace atomic absorption spectrophotometryGB/T17138-19970.01
0.5
Cadmium (Cd)
lead (Pb)
Graphite furnace atomic absorption spectrophotometryGB/T17141-19970.01
0.2
Arsenic (As)Atomic fluorescence SpectrophotometerGB/T221050.02
Table 2. PF with different classifications.
Table 2. PF with different classifications.
PFLevel of Pollution
PF ≤ 1Low pollution
1 < PF ≤ 3Moderate pollution
3 < PF ≤ 6Considerable pollution
PF > 6Very high pollution
Table 3. Geoaccumulation index with different classifications.
Table 3. Geoaccumulation index with different classifications.
LevelIgeoLevel of Pollution
Igeo ≤ 0Not to weakly contaminated
0 < Igeo ≤ 1Weakly to moderately contaminated
1 < Igeo ≤ 2Moderately contaminated
2 < Igeo ≤ 3Moderately to strongly contaminated
3 < Igeo ≤ 4Strongly contaminated
4 < Igeo ≤ 5Strongly to extremely contaminated
Igeo ≥ 5Extremely contaminated
Table 4. Descriptive statistics of soil PTEs in the study area (mg/kg). Min = minimum value; Max = maximum value; Med = median value; AM = arithmetic mean; SD = standard deviation; CV = coefficient of variation; Sk = skewness; Ku = kurtosis; k-s test = Kolmogorov–Smirnov test; BG = background values of Yunnan Province (CNEMC, 1990); RSV = risk screening values of the soil PTEs (GB15618-2018).
Table 4. Descriptive statistics of soil PTEs in the study area (mg/kg). Min = minimum value; Max = maximum value; Med = median value; AM = arithmetic mean; SD = standard deviation; CV = coefficient of variation; Sk = skewness; Ku = kurtosis; k-s test = Kolmogorov–Smirnov test; BG = background values of Yunnan Province (CNEMC, 1990); RSV = risk screening values of the soil PTEs (GB15618-2018).
PTEsMinMaxMedAMSDCV (%)SkKuk-s TestBGRSV
Cd0.021.060.190.280.25191.31.9212.80600.220.3
Pb15.1511228.14 31.42 15.28 49.02.1687.119040.6120
Cu5181.344.647.5929.12961.21.5754.018038.38100
Zn15.4423987.74100.4639.66439.50.493−0.177089.7250
As1.2354.469.3912.3610.07681.52.0164.63018.430
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dong, C.; Zhang, H.; Yang, H.; Wei, Z.; Zhang, N.; Bao, L. Quantitative Source Apportionment of Potentially Toxic Elements in Baoshan Soils Employing Combined Receptor Models. Toxics 2023, 11, 268. https://doi.org/10.3390/toxics11030268

AMA Style

Dong C, Zhang H, Yang H, Wei Z, Zhang N, Bao L. Quantitative Source Apportionment of Potentially Toxic Elements in Baoshan Soils Employing Combined Receptor Models. Toxics. 2023; 11(3):268. https://doi.org/10.3390/toxics11030268

Chicago/Turabian Style

Dong, Chunyu, Hao Zhang, Haichan Yang, Zhaoxia Wei, Naiming Zhang, and Li Bao. 2023. "Quantitative Source Apportionment of Potentially Toxic Elements in Baoshan Soils Employing Combined Receptor Models" Toxics 11, no. 3: 268. https://doi.org/10.3390/toxics11030268

APA Style

Dong, C., Zhang, H., Yang, H., Wei, Z., Zhang, N., & Bao, L. (2023). Quantitative Source Apportionment of Potentially Toxic Elements in Baoshan Soils Employing Combined Receptor Models. Toxics, 11(3), 268. https://doi.org/10.3390/toxics11030268

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