Pollution Characteristics, Spatial Patterns, and Sources of Toxic Elements in Soils from a Typical Industrial City of Eastern China

: Soil pollution due to toxic elements (TEs) has been a core environmental concern globally, particularly in areas with developed industries. In this study, we sampled 300 surface (0–0.2 m) soil samples from Yuyao City in eastern China. Initially, the geo-accumulation index, potential ecological risk index, single pollution index, and Nemerow composite pollution index were used to evaluate the soil contamination status in Yuyao City. Ordinary kriging was then deployed to map the distribution of the soil TEs. Subsequently, indicator kriging was utilized to identify regions with high risk of TE pollution. Finally, the positive matrix factorization model was used to apportion the sources of the different TEs. Our results indicated that the mean content of different TEs kept the order: Zn > Cr > Pb > Cu > Ni > As > Hg ≈ Cd. Soil pollution was mainly caused by Cd and Hg in the soil of Yuyao City, while the content of other TEs was maintained at a safe level. Regions with high TE content and high pollution risk of TEs are mainly located in the central part of Yuyao City. Four sources of soil TEs were apportioned in Yuyao City. The Pb, Hg, and Zn contents in soil were mainly derived from trafﬁc activities, coal combustion, and smelting. Meanwhile, Cu was mainly sourced from industrial emissions and atmospheric deposition, Cr and Ni mainly originated from soil parental materials, and Cd and As were produced by industrial and agricultural activities. Our study provides important implications for improving the soil environment and contributes to the development of efﬁcient strategies for TE pollution control and remediation.


Introduction
Soil pollution by toxic metal (TEs) has attracted wide attention both domestically and internationally, with public concern mainly regarding the safety of agricultural products and human health [1][2][3][4][5][6]. Among the TEs, chromium (Cr), lead (Pb), cadmium (Cd), mercury (Hg), arsenic (As), copper (Cu), zinc (Zn), and nickel (Ni) were listed as priority monitoring pollutants regulated by the Chinese government. Owing to rapid advances in industrialization and urbanization during the last several decades, TE pollution in soil has become increasingly serious and widespread in China [7][8][9][10][11]. with an average annual temperature of 18.1°C and annual precipitation of 1835.2 mm in 2020. Yuyao City is one of the most developed cities or counties in China, featuring developed hardware products, home appliance manufacturing, and plastic product industries.

Study Area
Yuyao City is situated in the northeast of Zhejiang Province, eastern China (F with area of around 1500 km 2 . The Yuyao City climate has a subtropical monsoon with an average annual temperature of 18.1 ℃ and annual precipitation of 1835.2 2020. Yuyao City is one of the most developed cities or counties in China, featuring oped hardware products, home appliance manufacturing, and plastic product ind

Soil Sampling and Analysis
In this study, the process of soil sampling and analysis strictly refer to ational standard issued by the Ministry and Ecology and Environment of P Republic of China (The Technical Specification for soil Environmental moni HJ/T 166-2004, https://www.mee.gov.cn/ywgz/fgbz/bz/bzwb/jcffbz/200412/t2004 3367.shtml, accessed on 21 October 2021). We collected 300 surface soil sam -0.2 m) from Yuyao City in 2012 ( Figure 1) based on spatial coverage, land d soil groups. In some areas which have been known as heavily polluted re r an important production base for rice and others agricultural crops, we co d a more extensive survey, and thus more samples were collected from thes ns. The coordinates of the sampling sites were recorded with a portable GP soil samples were then transported to the laboratory and air-dried for severa

Soil Sampling and Analysis
In this study, the process of soil sampling and analysis strictly refer to the national standard issued by the Ministry and Ecology and Environment of People's Republic of China (The Technical Specification for soil Environmental monitoring, HJ/T 166-2004, https:// www.mee.gov.cn/ywgz/fgbz/bz/bzwb/jcffbz/200412/t20041209_63367.shtml, accessed on 21 October 2021). We collected 300 surface soil samples (0-0.2 m) from Yuyao City in 2012 (Figure 1) based on spatial coverage, land use, and soil groups. In some areas which have been known as heavily polluted regions or an important production base for rice and others agricultural crops, we conducted a more extensive survey, and thus more samples were collected from these regions. The coordinates of the sampling sites were recorded with a portable GPS. The soil samples were then transported to the laboratory and air-dried for several days at ambient temperature. After that, all soil samples were sieved with a 2 mm nylon sieve, meanwhile, plant residuals and gravels were removed. The concentrations of soil Cr, Pb, Cd, Cu, Zn, and Ni were measured by inductively coupled plasma mass spectrometry (ICP-MS). The contents of Hg and As in the soil were determined by atomic fluorescence spectrometry (AFS). Reagent blanks and standard reference materials (GBW07405, supplied by IGGE at the Chinese Academy of Geological Sciences) were used for quality assurance and quality control with recovery rates varying between 90% and 110%.

Geo-Accumulation Index (I geo )
The I geo is calculated to evaluate the enrichment status of soil TEs by comparing soil background values. I geo is determined by the following formula [45]: here, C i means the content of soil TE i, C 0 represents the corresponding soil background value [46]. The I geo is classified into seven classes: I geo < 0 (unpolluted), 0 < I geo ≤ 1 (unpolluted to moderately polluted), 1 < I geo ≤ 2 (moderately polluted), 2 < I geo ≤ 3 (moderately to heavily polluted), 3 < I geo ≤ 4 (heavily polluted), 4 < I geo ≤ 5 (heavily to extremely polluted), and I geo > 5 (extremely polluted).

Potential Ecological Risk Assessment (ER and RI)
The ER and RI combine the concentration of each TE with its toxicity efficiency [47]. The ER is calculated using the following formula below [48]: where C i f indicates the contamination factor for TE i, C i 0 denotes the content of soil TE i, C i n on behalf of the background value of TE i issued by the China National Environmental Monitoring Center [46], ER i means the single potential ecological risk of TE i l, and TR i is the toxic coefficient of TE i [47]. The toxic coefficients of Cr, Pb, Cd, Hg, As, Cu, Zn, and Ni were 2, 5, 30, 40, 10, 5, 1, and 3, respectively. RI is the composite potential ecological risk index, which is an integrated value of several TEs. ER is classified into five grades [48]: ER < 40 (low risk), 40 < ER < 80 (moderate risk), 80 < ER < 160 (considerable risk), 160 < ER < 320 (high risk), and ER > 320 (very high risk). The RI is classified into four categories [48]: RI < 150 (low risk), 150 < RI < 300 (moderate risk), 300 < RI < 600 (considerable risk), and RI > 600 (high risk).

Single Pollution Index (PI) and Nemerow Composite Pollution Index (NPI)
The PI and NPI were also computed to evaluate the pollution status of TEs by comparing the content of soil TEs with national screening values regulated by the Chinese government [49]. The PI was calculated using the formula: here, C i on behalf of the concentration of soil TE i, and S i means the national screening value of TE i [49].

Ordinary Kriging
The spatial variability of soil TEs is a core concern for soil TE pollution control and remediation. As one of the most classic interpolation algorithms, the OK method has been widely used to estimate the spatial variation of soil properties around the world [29]. In this study, the OK was deployed to map the content of TEs in the soil around the survey region. In OK, the experimental semi-variogram was fitted to reveal the spatial dependence of soil TEs using formula as follows [51]: where Z(x i ) is the observed TE content at sampling site i, and n is the number of observed soil samples. ϕ i denotes the optimal weight value that could contribute to the unbiased estimation of the TE content with minimum variance. Detailed information on OK is available from Webster and Oliver (2007) [52].

Indicator Kriging
Indicator kriging (IK) is an extension of OK that can estimate the probability that the target variable exceeds a critical threshold at unvisited points [32,53,54]. The IK was deployed on a binary-transformed sample dataset. First, the indicator is defined using the following formula [52]: Subsequently, using Equation (6), the original variable is transformed into a new set of binary indicator variable. Then, the experimental semi-variogram is as follows: where γ * 1 (h) denotes the experimental semi-variogram of the newly transformed binary indicator. N(h) means the number of pairs of indicator transformations, I(χ i ; Z k ) and I(χ i + h; Z k ), with a distance lag of h.
The conditional cumulative distribution function is produced by the IK at unvisited sites: here, I * (χ 0 ; Z k ) on behalf of the predicted indicator at the unvisited sites χ 0 , and λ 1 represents the weight assigned to the indicator in sites x i .

US-EPA Positive Matrix Factorization (PMF)
The PMF (EPA PMF version 5.0) model is a multivariate factor analysis method [55], which was used in current research to identify the potential sources of soil TEs in Yuyao City. In PMF, the soil sample data are decomposed into two matrices: factor contributions and factor profiles. After that, the users could appoint the sources of TEs that may contribute to the samples based on source profile information and emissions inventories. A detailed description of the PMF method is available in the study organized by

Data Analysis
In this study, descriptive statistical analysis and correlation analysis were finished by R software (R Development Core Team, 2017). OK and IK were carried out with ArcGIS

Summary Statistics of Soil TEs
The descriptive statistics of the TEs in Yuyao City are listed in Table 1. The averaged concentration of different TEs kept the order: Zn (97.61 mg/kg) > Cr (61.80 mg/kg) > Pb (34.43 mg/kg) > Cu (32.75 mg/kg) > Ni (26.52 mg/kg) > As (6.13 mg/kg) > Hg (0.20 mg/kg) ≈ Cd (0.20 mg/kg). The averaged concentration of all TEs were greater than corresponding soil background values in Zhejiang Province, in addition to As. In comparison with the national arithmetic averaged concentrations of different TEs in China reported by Teng et al. (2014) [59], the levels of Pb, Cd, Hg, Cu, and Zn in Yuyao City were clearly higher than their national values. Cr content was slightly higher than the corresponding national value (Table 1). Moreover, the contents of As and Ni were lower than their national values (Table 1). The coefficient of variation (CV, %) was used to measure variation in the target variable [60,61]. The CV of different TEs followed a decreasing order: Hg (70.00%) > Cd (55.00%) > Cu (49.47%) > Ni (41.78%) > Cr (35.79%) > As (29.36%) > Zn (28.57%) > Pb (4.59%), indicating that Hg and Cd exhibited strong variation and suggested significant spatial heterogeneity of Hg and Cd in Yuyao City [61]. Our results are consistent with the results obtained in many other similar studies, such as those by Xia

Correlation Analysis
In this study, we used the Spearman's correlation method to analyze the associations among the different TEs. As presented in Figure 2, most TEs under study showed a significant positive correlation with each other in Yuyao City. In particular, Cr had a strong relationship with Ni (0.94), Pb had a strong relationship with Hg (0.70), and Cu had a strong relationship with Zn (0.63). TEs are highly related to each other and may share similar sources which could contribute to the identification of sources of soil TEs and more efficient control of soil TE contamination [47]. similar sources which could contribute to the identification of sources of soil TEs and more efficient control of soil TE contamination [47].

Figure 2.
Relationships among TEs in soil (n = 300) (the symbol × means that the relationships among corresponding TEs was not significant at the 99% confidence level).

The Geo-Accumulation Index
The Igeo was calculated to quantitatively evaluate the accumulation of soils in the study area. A boxplot of Igeo is shown in Figure 3. The average Igeo values for different TEs followed the order: Cd (0.84) > Hg (0.27) > Cu (0.19) > Pb (-0.13) > Zn (-0.17) > Cr (-0.46) > Ni (-0.59) > As (-1.24). The results revealed that evaluated by the geo-accumulation index, the pollution status for Pb, Zn, Cr, Ni, and As was at an unpolluted level, while the pollution status for Cd, Hg, and Cu was unpolluted to moderately polluted.  Relationships among TEs in soil (n = 300) (the symbol × means that the relationships among corresponding TEs was not significant at the 99% confidence level).

The Geo-Accumulation Index
The I geo was calculated to quantitatively evaluate the accumulation of soils in the study area. A boxplot of I geo is shown in The results revealed that evaluated by the geo-accumulation index, the pollution status for Pb, Zn, Cr, Ni, and As was at an unpolluted level, while the pollution status for Cd, Hg, and Cu was unpolluted to moderately polluted.
Land 2021, 10, 1126 7 of 20 similar sources which could contribute to the identification of sources of soil TEs and more efficient control of soil TE contamination [47].

Figure 2.
Relationships among TEs in soil (n = 300) (the symbol × means that the relationships among corresponding TEs was not significant at the 99% confidence level).

The Geo-Accumulation Index
The Igeo was calculated to quantitatively evaluate the accumulation of soils in the study area. A boxplot of Igeo is shown in Figure 3 The results revealed that evaluated by the geo-accumulation index, the pollution status for Pb, Zn, Cr, Ni, and As was at an unpolluted level, while the pollution status for Cd, Hg, and Cu was unpolluted to moderately polluted.   Table 2 shows the proportions of different geological accumulation grades of different TEs in the whole survey region, referring to the I geo classification standard [35]. Obviously, the soil pollution caused by Cd and Hg is much more serious than that caused by Cr, Pb, As, Cu, Zn, and Ni. Regarding As, all soil samples were unpolluted, while most of the soil samples were at unpolluted levels for Cr, Cu, Zn, and Ni (Table 2). In terms of Pb, 47.33% of the soil samples belonged to the unpolluted to moderately polluted grade, while the remaining soil samples were unpolluted. As for Cd, 28.33%, 1.33%, and 0.33% of the soil samples were moderately polluted, moderately to heavily polluted, and heavily to extremely polluted, respectively. With regard to Hg, 28.00% and 2.33% of soil samples were moderately polluted and moderately to heavily polluted, respectively.  The ER values of different TEs and RI were used to quantitatively evaluate the potential ecological risk due to TE accumulation. The mean value of ER for different TEs followed the order: Hg (89.57) > Cd (86.32) > Cu (9.30) > Pb (7.26) > As (6.67) > Ni (3.23) > Cr (2.34) > Zn (1.38) (Figure 4). The results showed that soil in Yuyao City had considerable potential ecological risk results from Hg and Cd accumulation, while the potential ecological risk caused by other TEs was low (Figure 4). The mean RI value was 206.08, which is at the moderate risk grade.     The proportions of soil samples at different potential ecological risk levels for different TEs are listed in Table 3. The potential ecological risk of all soil samples was low for Cr, Pb, As, Cu, Zn, and Ni. With regard of Cd, the proportions of ER value of soil samples at low, moderate, considerable, high, and very high risk were 3%, 49%, 45.67%, 2%, and 0.33%, respectively. For Hg, the proportions of ER value of soil samples at low, moderate, considerable, high, and very high risk were 24.33%, 29.67%, 34.33%, 11.33%, and 0.33%, respectively. For the RI value, the proportions of soil samples at low, moderate, considerable, and high risk were 35.67%, 49.33%, 14.67%, and 0.33%, respectively.

Pollution Index
Moreover, we calculated the PI and NPI to evaluate the contamination status of TEs in the soil of Yuyao City. As presented in Figure 5, the averaged PI decreased by the following order: Cd (0.48) > Zn (0.42) > As (0.34) ≈ Hg (0.33) ≈ Pb (0.33) > Ni (0.29) > Cr (0.23) ≈ Cu (0.23). The averaged PI of all TEs were at a safe level, but it is worth noting that the PI of some soil samples were at polluted levels for Cd, Hg, Cu, Zn, and Ni. Additionally, the average NPI value was 0.46, which was also at a safe level.    Table 4. Similar to the I geo and ER indices, the PI of Cr, Pb, and As in all soil samples were at a safe level, while only a small part of the soil samples had PI values at polluted levels (slight pollution or more serious) for Cd (5.33%), Hg (3%), Cu (0.67%), Zn (1.33%), and Ni (0.33%). With regard to the NPI value, the proportion of soil samples with NPI values at safe, precautionary, slightly polluted, moderately polluted, and seriously polluted levels were 86.33%, 11%, 2.33%, 0%, and 0.33%, respectively.  Overall, Cd and Hg clearly accumulated in the soil in Yuyao City. Generally, the pollution was mainly caused by Cd and Hg, while the content of other TEs was kept at a safe level. This result is consistent with previous surveys conducted in Yuyao City [65][66][67]. Yuyao City is ranked 12th among the top 100 counties in China according to the level of economic development. Extensive industrial activities may lead to the accumulation of soil Cd and Hg [68,69].

Identify Regions with High Risk of TEs Contamination
Traditional spatial interpolation algorithms, such as OK and IDW, can only map the specific value of the target variable in unobserved locations. However, the predictions made by all interpolation methods are biased, which means that uncertainty exists in the predicted results. To fill this gap, IK was developed to map the probability that the target variable exceeds a specific threshold. In this study, we used IK to predict the probability of the occurrence of soil TE pollution in Yuyao City. For the geo-accumulation status, we focused on Cd and Hg because Cd and Hg are the most seriously polluted TEs in Yuyao City. Meanwhile, for the potential ecological risk and pollution index, we focused on NPI and RI because it can represent the composite pollution status and ecological risk caused by different TEs in the soil.
As shown in Figure 6a, the high pollution risk caused by Cd accumulation (probability of I ego > 1) was mainly expected in the central part of Yuyao City. In some regions, the probability of Cd pollution could exceed 90%. In terms of Hg, regions with a high pollution risk were situated in the center of Yuyao City, but most of the regions had a pollution risk lower than 50%. For the NPI value, areas with high composite pollution risk were mainly situated in the central part of Yuyao City, while the high value of probability of potential ecological risk (RI > 150) occurred mainly in the central and middle parts of Yuyao City.
In summary, the central part of Yuyao City had a relatively higher risk of soil TE contamination, especially for Cd and Hg. This indicates that these regions should be listed as priority regions for soil TE contamination control and remediation. There is an urgent need to take some measures to protect residents from the health risks caused by soil TE contamination. For example, the crop planting system could be adjusted in regions with high pollution risk, and economic plants could be planted take the place of food crops to reduce the accumulation of TEs in the human body via the food chain. Planting crops with low accumulation of TEs is another efficient way to reduce the risk of TE accumulation.  In summary, the central part of Yuyao City had a relatively higher risk of soil TE contamination, especially for Cd and Hg. This indicates that these regions should be listed as priority regions for soil TE contamination control and remediation. There is an urgent need to take some measures to protect residents from the health risks caused by soil TE contamination. For example, the crop planting system could be adjusted in regions with high pollution risk, and economic plants could be planted take the place of food crops to reduce the accumulation of TEs in the human body via the food chain. Planting crops with low accumulation of TEs is another efficient way to reduce the risk of TE accumulation.

Spatial Pattern of TEs
Obtaining accurate information on the spatial variability of soil TEs plays a critical role in controlling TE pollution. In current research, OK interpolation was conducted to estimate the spatial variation of TE concentration in the soil of Yuyao City (Figure 7a-h). As presented in Figure 7, the concentrations of Cr, Cd, Hg, As, Zn, and Ni showed similar spatial patterns, with high values detected in the central and southern parts of Yuyao. Pb and Cu also had similar spatial patterns, with high values majorly detected in the central part of Yuyao City. According to the map of land use of Yuyao City presented in Figure 7i, the land use types in the central part of Yuyao are mainly urban, rural, industrial, mining, and residential land. Anthropogenic activities, such as industrial activities (emission of solid waste and wastewater), traffic activities (vehicle emissions), and agricultural activities (application of chemical fertilizer and pesticide) discharge large amounts of waste in which large amounts of TEs exist and then accumulate in soil [8,[70][71][72][73][74].

Spatial Pattern of TEs
Obtaining accurate information on the spatial variability of soil TEs plays a critical role in controlling TE pollution. In current research, OK interpolation was conducted to estimate the spatial variation of TE concentration in the soil of Yuyao City (Figure 7a-h). As presented in Figure 7, the concentrations of Cr, Cd, Hg, As, Zn, and Ni showed similar spatial patterns, with high values detected in the central and southern parts of Yuyao. Pb and Cu also had similar spatial patterns, with high values majorly detected in the central part of Yuyao City. According to the map of land use of Yuyao City presented in Figure  7i, the land use types in the central part of Yuyao are mainly urban, rural, industrial, mining, and residential land. Anthropogenic activities, such as industrial activities (emission of solid waste and wastewater), traffic activities (vehicle emissions), and agricultural activities (application of chemical fertilizer and pesticide) discharge large amounts of waste in which large amounts of TEs exist and then accumulate in soil [8,[70][71][72][73][74].

Source Apportionment by PMF
PMF has been commonly used to apportion the source of soil TEs [57,[75][76][77][78]. In current research, the PMF method was deployed to compute the contribution of different sources to the accumulation of eight soil TEs in Yuyao City. As displayed in Figures 8 and

Source Apportionment by PMF
PMF has been commonly used to apportion the source of soil TEs [57,[75][76][77][78]. In current research, the PMF method was deployed to compute the contribution of different sources to the accumulation of eight soil TEs in Yuyao City. As displayed in Figures 8 and 9, four sources were identified by the PMF model.

Source Apportionment by PMF
PMF has been commonly used to apportion the source of soil TEs [57,[75][76][77][78]. In c rent research, the PMF method was deployed to compute the contribution of differ sources to the accumulation of eight soil TEs in Yuyao City. As displayed in Figures 8 a  9, four sources were identified by the PMF model.   Factor 1 was heavily loaded with Pb (59.16%), Hg (75.96%), and Zn (35.05%). As p sented in Figure 2, Pb was significantly and highly correlated with Hg and Zn, which al verified that Pb, Hg, and Zn were derived from similar sources in the soil of Yuyao Ci Numerous studies have verified that vehicle exhaust emissions are a critical source of accumulation in the soil [79,80]. As listed in Table 1, Hg had the highest CV value, whi means that anthropogenic activities are the major origin of soil Hg [21]. Industrial acti ties such as smelting, coal mining, and combustion are regarded as the main sources Hg pollution in the soil [21] In addition, coal mining and combustion emissions are i portant sources of Pb in the soil. Moreover, previous studies have also reported that Zn closely associated with anthropogenic actions, such as burning fossil fuels, smelting, a coal combustion emissions [81,82]. The high values of Pb, Hg, and Zn were mainly sit ated in the central part of Yuyao City, which overlays with the urban land. Thus, Facto could represent the source of traffic activities, coal combustion, and smelting ( Figure 6 Factor 2 was dominated by Cu, with a very high loading of 71.69%. Generally, Cu closely related to the emission of industrial activities like metal smelting, petrochemic cement production, and so on [83,84]. In addition, some researchers have also found th atmospheric deposition is another significant contributor to Cu accumulation in the s [57,85]. Therefore, Factor 2 was mainly derived from industrial emissions and atmo pheric deposition. Factor 3 was dominated by Cr and Ni with high factor loading values of 74.28% a 76.06%, respectively. As presented in Figure 2, the content of Cr was extremely high, lated to the content of Ni in the soil of Yuyao City, with a correlation coefficient of 0. This demonstrates that Cr and Ni in the soil originated from the same source. The me Factor 1 was heavily loaded with Pb (59.16%), Hg (75.96%), and Zn (35.05%). As presented in Figure 2, Pb was significantly and highly correlated with Hg and Zn, which also verified that Pb, Hg, and Zn were derived from similar sources in the soil of Yuyao City. Numerous studies have verified that vehicle exhaust emissions are a critical source of Pb accumulation in the soil [79,80]. As listed in Table 1, Hg had the highest CV value, which means that anthropogenic activities are the major origin of soil Hg [21]. Industrial activities such as smelting, coal mining, and combustion are regarded as the main sources of Hg pollution in the soil [21] In addition, coal mining and combustion emissions are important sources of Pb in the soil. Moreover, previous studies have also reported that Zn is closely associated with anthropogenic actions, such as burning fossil fuels, smelting, and coal combustion emissions [81,82]. The high values of Pb, Hg, and Zn were mainly situated in the central part of Yuyao City, which overlays with the urban land. Thus, Factor 1 could represent the source of traffic activities, coal combustion, and smelting ( Figure 6). Factor 2 was dominated by Cu, with a very high loading of 71.69%. Generally, Cu is closely related to the emission of industrial activities like metal smelting, petrochemical, cement production, and so on [83,84]. In addition, some researchers have also found that atmospheric deposition is another significant contributor to Cu accumulation in the soil [57,85]. Therefore, Factor 2 was mainly derived from industrial emissions and atmospheric deposition.
Factor 3 was dominated by Cr and Ni with high factor loading values of 74.28% and 76.06%, respectively. As presented in Figure 2, the content of Cr was extremely high, related to the content of Ni in the soil of Yuyao City, with a correlation coefficient of 0.94. This demonstrates that Cr and Ni in the soil originated from the same source. The mean content of Cr and Ni in Yuyao City was slightly higher than the corresponding background value, revealing that the I geo and ER values of Cr and Ni were also at the safe or low risk level, which means that the Cr and Ni in the soil of Yuyao City may result from the weathering of soil parental materials. Previously, many researchers, such as Salonen and Korkka-Niemi  [57] also reported that the contents of Cr and Ni were mainly sourced from soil parental materials. Therefore, Factor 3 represented natural sources, such as soil parent materials.
Factor 4 was mainly loaded with As and Cd with loadings of 50.84% and 40.56%, respectively. The enrichment of As and Cd has been shown to be associated with agricultural activities, like applying chemical fertilizers, herbicides, and pesticides [89]. In addition, the utilization of wastewater for irrigation could also give rise to the accumulation of As in the soil [89,90]. Some studies also reported that As could probably be sourced from poultry farming because organic As is widely used as a feed additive for livestock to promote animal growth, which could then result in a high content of As in livestock manure [91,92]. In terms of Cd, numerous studies have proved that Cd originates from industrial activities such as smelter plant emission, discharge of industrial waste, and so on. Cd is also usually used as an input for industrial production of various products such as plastics [93,94]. The plastic industry is a pillar industry in Yuyao. Thus, the plastic industry may also contribute to Cd accumulation in the soil of Yuyao City. Therefore, Factor 4 represents the source related to industrial and agricultural activities.
In summary, the main sources of TEs in the soil of Yuyao City were industrial activities, industrial emissions and atmospheric deposition, natural sources, and agricultural activities. The results of source apportionment of soil TEs were also verified by the results of summary statistics listed in Table 2, correlation analysis (Figure 2), and the map of the content of soil TEs in Yuyao City produced by OK ( Figure 6) interpolation and the map of pollution risk outputted by IK (Figure 7). This confirmed the feasibility of using PMF for quantitatively analyzing the sources of soil TEs. However, as reported by previous researchers, uncertainty still exists in the results produced by the PMF model, which highlights the importance of introducing and combined with other source apportionment methods such as cluster analysis [12], principal component analysis [90], machine learning methods [17], the Unmix model [95], the stable isotope ratio technique [96], and the finite mixture distribution model [97] to reduce the uncertainty of the results produced by PMF.

Implication and Policy Recommendation
The results we obtained in from presented research could provide important implication for decision-makers to make effective management measures to control soil TEs contamination and thus improve environmental quality. Firstly, regions with high pollution risk should be listed as priority areas for soil TEs contamination control and remediation. Some management practices, such as lime application and bioleaching, could be employed in acidic soils to regulate soil pH and reduce the availability of HMs [98][99][100][101]. Secondly, we could consider establishing a spatial buffer around industrial park or regions featured by high density of industrial companies. We should also keep a safe distance between polluted land or industrial park and residents. Thirdly, breeding crop cultivars with low accumulation is another promising way to mitigate soil TEs contamination in food crops. Crops with different cultivars showed strong variation in ability of absorbing TEs from soil [7]. Finally, according to the results presented in Section 3.6, burning fossil fuels, traffic emissions, and coal combustion emissions contributed to Pb, Hg and As accumulation in the study area. Therefore, the government could encourage adopting green travel methods like bike or public traffic as well as promote use of clean energy such as wind energy, hydropower, tidal energy, and solar energy. These measures are expected to reduce accumulation of TEs in soil through atmospheric deposition, emission of wastewater, and solid waste. In our further work, comparison with previous similar studies which made surveys in the same regions as well as conducting another soil sampling in the same sites is a promising way to detect the temporal change of TEs in soil.

Conclusions
Obvious accumulation of Cr, Pb, Cd, Hg, Cu, Zn, and Ni was detected in the soil of Yuyao City. The concentrations of Cr, Pb, As, Cu, Zn, and Ni in most of the soil samples were at safe levels, whereas clear pollution caused by Cd and Hg accumulation in soil was detected. The high contents of the different TEs were mainly situated in the central and middle parts of Yuyao City. In addition, there was a low pollution risk caused by Cr, Pb, As, Cu, Zn, and Ni in Yuyao City, while a high risk of Cd and Hg contamination was detected. The regions with a high risk of TE pollution were mainly distributed in central Yuyao City.
In terms of the sources of different TEs in the soil, Pb, Hg, and Zn were mainly derived from traffic activities, coal combustion, and smelting, Cu mainly sourced from industrial emissions and atmospheric deposition, Cr and Ni mainly came from soil parental materials, and Cd and As were attributed to industrial activities and agricultural activities. This study could provide critical information for improving soil environmental quality and contribute to the development of more efficient policies for reducing TE accumulation in soil.