Health risk assessment of heavy metal pollution and its sources in agricultural soils near Hongfeng Lake in the mining area of Guizhou Province, China

Background Accelerated modern industrial processes, extensive use of pesticides and fertilizers and remaining issues of wastewater irrigation have led to an increasingly severe composite pollution of heavy metals in arable land. Soil contamination can cause significant damage to ecological environments and human health. Mineral resource mining can result in varying degrees of heavy metal pollution in surrounding water systems and soil. As a plateau lake, Hongfeng Lake has a fragile watershed ecosystem. Coupled with the rapid development of the current socio-economy and the ongoing activities of mining, urbanization and agricultural development, the water and soil environment of the lake and arable land are facing serious heavy metal pollution. Therefore, the situation warrants attention. Methods This study focused on characterizing soil types and conducted sampling and laboratory testing on the farmland soil in Hongfeng Lake. The integrated Nemero comprehensive pollution assessment and potential ecological pollution assessment methods were used to evaluate the heavy metal pollution status. The APCS-MLR model was employed to explore the sources of heavy metal pollution. In addition, the human health risk model was used to analyze the association between heavy metal content in cultivated land and human health risks. Results The single-factor pollution of each element was ranked in descending order: Hg > As > Pb > Cr > Cd, with Hg being the main pollutant factor. The entire area was subjected to mild pollution according to the pollution index. Pollution source analysis indicated two main pollution sources. Hg, As, Pb and Cr pollution mainly resulted from Source 1 (industrial and natural activities), accounting for 71.99%, 51.57%, 67.39% and 68.36%, respectively. Cd pollution was mainly attributed to Source 2 (agricultural pollution source), contributing 84.12%. The health risk assessment model shows that heavy metals posed acceptable carcinogenic risks to humans rather than non-carcinogenic risks. As was the main non-carcinogenic risk factor, while Cr was the main carcinogenic risk factor, with higher risks in children than adults. Conclusion Our study identified the heavy metal pollution in farmland soil in Hongfeng Lake, evaluated and analyzed the pollution sources and identified the heavy metal elements in cultivated lands that have the greatest impact on human health risks. The aim of this study is to provide a scientific basis for soil heavy metal pollution control.


Introduction
Accelerated modern industrialization, extensive pesticide and fertilizer use and longstanding issues of sewage irrigation have led to increasingly severe heavy metal complex pollution in agricultural soils (1)(2)(3)(4).Mining activities can also cause varying levels of heavy metal pollution in water systems and agricultural soils.For example, research on a copper mine in the Jiulong section of the Yalong River basin found that the Cu element was more likely to accumulate and cause pollution in some upstream water systems.These systems were closer to the mining point, at higher elevations, in steeper terrain and with smaller water convergence (5).In addition, it has also been found that heavy metal contents in farmland soils surrounding the Nhue River mining area in Vietnam far exceeded background values, and Cd content exceeded the heavy metal standard for Vietnamese agricultural soils (6).Fry found that soils at a copper smelter in Namibia were contaminated with As, Cu, and Pb, and surface dust on buildings posed a health risk to humans (7).Therefore, it is necessary to assess heavy metal pollution in soils surrounding mining areas in watersheds in order to address pollution-induced risks.
This study took the Hongfeng Lake basin as the study area.Hongfeng Lake was formed in 1958 as an artificial lake to construct the Maotiao River cascade hydropower station with fragile watershed ecosystems.It is the largest source of drinking water in Guiyang City.With the rapid development of the current socio-economy, mining, urbanization construction and agricultural development activities, the water environment of the lake is in a critical situation.Hongfeng Lake is mainly fed by the Maiweng River, the Yangchang River, the Maxian River, and the Houshui River.There are many factories and enterprises in the upstream watershed, such as chemical, metallurgical and mining enterprises.The watershed has been receiving industrial and mining wastewater for a long time, and the watershed soils have been polluted to a certain extent (8,9).Soil heavy metals can migrate to aquatic ecosystems, threatening the drinking water of watersheds.Since soil environmental conditions interact with water quality in watersheds, soils in watersheds serve as a region of aquatic ecosystems.Excessive soil heavy metals can inhibit soil functions, alter soil physicochemical properties, disrupt nutrient supply and balance and impact the ecological environment and crop growth (10)(11)(12).This can threaten crop growth and development, agricultural production security and human health through various exposure pathways during human labor.For the sake of human society, life and ecology, it is vital to evaluate and explore the soil pollution status and sources in the Hongfeng Lake basin.We should pay attention to the agricultural soil safety in the Hongfeng Lake basin.Studies should focus on refining soil types and selecting cultivated land soils that have the most significant impact on ecosystems and human health in the watershed.The pollution levels and sources of the five heavy metals (Hg, As, Pb, Cr, and Cd) in the cropland soil of the watershed should also be clarified.
Heavy metal pollution sources have been widely explored using receptor and diffusion models.The diffusion models require information on discharge data and related information of each pollution source to predict the spatiotemporal changes of individual pollutants.The receptor models require the chemical content and related microscopic analysis of each collected sample.Establishing the corresponding relationship between pollutant sources and pollution status can help to identify the main pollution sources and their contributions to the receptors (13).The receptor models are more widely used with fewer restrictions, mainly including the chemical mass balance model (CMB), principal component analysis (PCA) and positive matrix factorization (PMF) models (14,15).In particular, the absolute principal component score combined with multiple linear regression analysis (APCS-MLR) focuses on receptor pollution sources.This technique does not require knowledge of the number of pollution sources and can easily and objectively determine the specific contribution of pollution sources to the receptors with simple operating conditions (16,17).The method first proposed by Thurston and Spengler in 1985 focuses on obtaining APCS of factors through factor analysis and then calculating the contribution rate of common factors to pollutants using APCS-MLR (18).This model has been widely used to analyze pollution sources in air, surface water and sediment (19,20).
Therefore, this study refined the soil types and selected cropland soils affecting ecology and human health the most as the research object.The pollution level of five soil heavy metals in the cropland of the watershed was evaluated.The APCS-MLR receptor model was used to analyze soil heavy metal sources.Combined with the spatial distribution and correlation analysis of heavy metal pollution, pollution levels, sources and distribution characteristics of heavy metals in the cropland in the mining area of the Hongfeng Lake basin were clarified.The human health risk impact assessment method defined by the US Environmental Protection Agency (USEPA) was used to evaluate the health risks of heavy metals to humans.This study can provide a reference and scientific basis for preventing and controllinh farmland heavy metal pollution and health risks in Hongfeng Lake.

Study area
The Hongfeng Lake basin is located in the suburbs of Qingzhen City with more than 30 types of mineral deposits.These deposites include bauxite, hematite, pyrite, coal and marble.Bauxite is one of the primary mineral sources of Guizhou Aluminum Corporation.Hongfeng Lake is mainly fed by the Maiweng River, the Yangchang River, the Maxian River and the Houshui River.The reservoir has an area of 57.2 km 2 , a water level of about 1,240 m, a dead water level elevation of 1,127.5 m and a capacity of 752.9 million m 3 .There are 67 water supply points.Hongfeng Lake supplies 250,000 tons of water daily to Guiyang City.

Sampling, processing, and analysis
The sampling point is shown in Figure 1.From August to December 2020, after preliminary investigations, sampling points were arranged radially from Shangjian Street in the upstream station of Hongfeng Lake to Hongfeng Town, Gaofeng Town and Machang Town downstream.The sampling points were set up based on the geographical location and climatic conditions of the mining area.These points were in the farmland in the vicinity of the mine, tailings and beneficiation plants using the plum blossom pattern.Individual soil points consisting of 5 sub-samples were taken at 20 cm using a small shovel.A total of 190 soil samples were collected.Next, the soil samples were placed in a polyethylene bag, crushed and sealed into bottles for subsequent analysis.Heavy metals were determined primarily using reagent blanks and national standard soil samples (GSS-3).One sample was randomly selected from every 30 samples for triple parallel analysis to ensure quality control.Cd, As, Pb and Cr contents were measured using an ICAP RQ inductively coupled plasma mass spectrometer (ICP-MS).The microwave digestion was performed using aqua regiahydrofluoric acid with high chloric acid.The Hg content was determined using an AFS-8220 atomic fluorescence spectrophotometer.The digestion was conducted in a boiling water bath for 2 h after the addition of aqua regia.Mercury vapor was carried by the carrier gas into the atomizer for measurement.Standard quality assurance and quality control were strictly adhered to throughout the study.The recovery rate of heavy metals was between 70 and 125%, and the experimental results complied with the quality control standards.A single pollutant can be quantitatively evaluated to determine the degree of exceedance of individual pollutants relative to background values and clarify pollution levels and main pollutants:
An evaluation method based on a single factor pollution index can consider both the maximum and average values of individual pollutants.This method can comprehensively evaluate the overall pollution situation in the study area: where P i represents the single factor pollution index of pollutant i in soils; C i is the measured contamination value of pollutant i; S i denotes the assessment standard of pollutant i (21, 22); P n is the comprehensive pollution index for heavy metals in Inner Mongolia; P imax and P ave are the maximum and average value among single-factor pollution indices of each pollutant, respectively.In this study, P i was calculated for the five elements (Cd, Hg, Pb, As, and Cr).

Ecological risk assessment
The potential ecological risk index method proposed by Lars (23) can comprehensively consider the ecotoxicity of pollutants and ecological environmental factors.This method can also reflect the combined effect of various pollutants on ecology.This index can be used to quantitatively analyze and predict potential ecological risks: where RI is the comprehensive potential ecological risk index, E r i is the individual potential ecological risk index of heavy metal i, T i is the toxicity coefficient of heavy metal i, P i is the single factor index value for heavy metal i.The toxicity response coefficients of Cd, Hg, As, Pb, and Cr were 30, 40, 10, 5, and 2, respectively.The classification based on the heavy metal pollution assessment method is shown in Table 1.

APCS-MLR receptor model
The APCS-MLR receptor model can be obtained by subtracting the principal component score of a zero concentration sample from the principal component scores from principal component analysis (24,25).The calculation process includes three steps: Step (1): to standardize the raw data, introduce the zero concentration factor and calculate the APCS: where Z ij is the standardized value (dimensionless); C ij is the content of element j in the i sample (mg/kg); C J and σ j are the mean and standard deviation of C ij (mg/kg).
Step (2): multiple linear regression (MLR) analysis is performed using the APCS and heavy metal content as the independent and dependent variables, respectively.The contribution rate of individual pollutant sources is then calculated using regression coefficients.

C b APCS a
where C i is the content of the measured heavy metal element (mg/ kg); a im represents the regression coefficient of source m on heavy metal element i; APCS im represents the absolute principal component score of heavy metal element i for all samples; b i is a constant term in multiple linear regression; APCS im × a im represents the contribution of source m to C i in the soil sample.
Step (3): to avoid negative values affecting the results in the APCS-MLR modeling, the contribution rate of pollutant source m to heavy metal element i is calculated using the absolute value: The contribution rate of a known pollutant source is expressed as The contribution rate of an unknown pollutant source is expressed as where PC im is the contribution rate for heavy metal element i and pollutant source m and APCS im represents the average of the absolute principal component scores of the entire sample for heavy metal element i.

Heavy metal exposure models and parameters
Hazardous substances can enter human bodies through different transmission pathways.Health risk assessment is an effective method of determining the health risk level of various pollutants.It can evaluate the nature and possibility of the adverse effects on the human body caused by the chemical content of certain media in the current or future environment (26).Referring to the health risk model of the USEPA, the health risks of residents around the mining area were evaluated.The exposure pathways included oral ingestion, respiratory intake and skin contact (27)(28)(29)(30).Human health risks were classified into noncarcinogenic and carcinogenic risks.This study evaluated the chronic non-carcinogenic risk of Cd, Hg, As, Pb and Cr and the carcinogenic risk of Cd, As, and Cr.The average daily intake of heavy metals by humans under three exposure pathways can be calculated using Oral ingestion: Inhalation: Dermal contact: where C i is the content of heavy metal i in soils (mg/kg); other parameters refer to the USEPA health risk assessment standards and relevant research worldwide and are listed in Table 2.The chronic non-carcinogenic risk for human bodies exposed to soil heavy metals is calculated using: where i is the heavy metal species, j represents different exposure pathways, ADD (average daily dose) refers to the longterm daily exposure dose of heavy metals, RfD (reference dose) refers to the reference dose of non-carcinogenic heavy metals for daily intake, HQ represents the non-carcinogenic risk of potential toxic heavy metal elements in soils to human bodies and is obtained from the ratio of ADD of each heavy metal element to the corresponding RfD and HI is the total non-carcinogenic risk value for multiple heavy metals or one heavy metal under different exposure pathways.

Carcinogenic risk
The human carcinogenic risk of heavy metals in soils is calculated as follows: where ADI (carcinogenic heavy metals) is the long-term daily average exposure dose of heavy metals [mg/(kg•d)]; SF (slope factor) refers to the slope factor of carcinogenic heavy metals; CR is the carcinogenic risk caused by the potential toxic heavy metal elements in soils to human bodies and can be obtained by multiplying ADD of each heavy metal element by its corresponding SF; TCR represents the total carcinogenic risk of multiple heavy metals or a single heavy metal under multiple exposure pathways (dimensionless).The values of RfD and SF are shown in Table 3 (30-34).

Statistical analysis and geospatial distribution
The SPSS 21.0 software was used for descriptive statistics, Pearson correlation analysis, principal component analysis and operation of the APCS-MLR receptor model on heavy metal data

Soil characterization
The sample size at the sampling points was n = 190.The pH value of agricultural soils in the study area ranged from 4.5 to 8.24, with an average of 6.20, indicating an overall acidic soil.According to Table 5, except for Cd, all other elements exceed the soil background values of Guizhou Province, with Hg having the highest exceedance (89.47%).The exceedance of Hg, As, and Pb was all higher than 50%.In contrast, the exceedance of Cd and Cr was relatively low.The coefficients of variation for individual elements were as follows: Hg (0.55) > As (0.65) > Cd (0.61) > Pb (0.31) > Cr (0.29).Except for the element Cr, the coefficients of variation for the other elements were all greater than 0.3, indicating moderate to high variability.The coefficients of variation for Cd, Hg and As fell between 0.5 and 1, indicating a higher degree of variability.This suggests that the variations in the Cd, Hg and As content in the study area were most significantly influenced by anthropogenic activities.

Pollution assessment results
The assessment results of heavy metal pollution in the cropland of the study area are shown in Table 6.The single factor pollution index method indicates that Hg was the main pollutant in the study area.The averages of the one-way pollution index of various elements were Hg (2.15) > As (1.31) > Pb (1.11) > Cr (0.99) > Cd (0.67).Except for the Cd element, the heavy metal contents at the sampling points lay within the clean range.However, the contents of other elements exceeded the background value to some extent, with As, Pb and Cr exhibiting slight exceedance.The mean values of the potential ecological pollution indices for individual heavy metals were Hg (86.01) > Cd (20.03) > As (13.08) > Pb (5.56) > Cr (1.97).The mean value of the comprehensive potential ecological risk index for the entire area was 126.65.This indicates that the Hg content in the study area will pose potential moderate ecological risks, and the other four types of elements will pose slight  ecological risks.The comprehensive potential ecological risk of the entire area was classified as slight.

Pearson correlation matrix
In order to understand the correlation between various heavy metals in the soil samples, Pearson's correlation analysis was conducted on the raw data from the sampling points.Table 7 shows that the correlation coefficients between Hg-As-Cr and As-Pb-Cr, two groups of combined elements, were above 0.5.They all exhibited a significant positive correlation (p < 0.01) (35).This indicates that Hg, As, Pb, and Cr had similar pollution sources.

Source apportionment using the APCS-MLR model
Principal component analysis was conducted using the SPSS software.The results of the Bartlett's spherical test and the KMO measure were 0.000 and 0.65, respectively.This indicates that it was suitable for principal component analysis.Combining the correlation analysis results, the standardized factors of the two extracted principal components were subjected to orthogonal rotation.The cumulative contribution rate reached 72.528%, which can represent most of the information on the research object.

APCS-MLR model
Through principal component analysis, the two principal component scores were converted into APCS.Then, in order to obtain the contribution rates of different pollution sources, the APCS and heavy metal content were used as independent and dependent variables, respectively, to conduct multiple linear regression analysis.The APCS-MLR receptor model results for individual heavy metal elements are shown in Table 8.

Health risk assessment 3.4.1. Non-carcinogenic risk
Table 9 shows that the HI and HQ values for both adults and children were less than 1.Therefore, the exposure to heavy metals in the study area posed a negligible non-carcinogenic risk to the human body within the acceptable range.

Carcinogenic risk
Table 10 shows that the CR and TCR values for both adults and children were less than 1.00E-04.Therefore, the exposure to heavy metals in the study area posed a negligible carcinogenic risk to the human body, within the acceptable range.

. Analysis of heavy metal pollution
In the preliminary analysis, this study performed a descriptive analysis on the heavy metal contents in the soil samples (Table 5).The results show that the variation coefficients of each element were in the order of Hg > As > Cd > Pb > Cr.Except for the element Cr, the variation coefficients of all other elements exceeded 0.3, indicating medium or high variations.The variation coefficients of Cd, Hg and As were between 0.5 and 1, indicating a high degree of variation (36).This suggests that the variations in the Cd, Hg and As contents in the study area were most significantly affected by human activities.
The single factor pollution index of individual elements was interpolated by inverse distance weighting to depict the pollution status of each element (Figure 2). Figure 2 shows that there were high pollution areas with point-like distributions for Cd, Hg, As and Pb in the pollution status map of individual elements.This indicates that these four elements were enriched due to human activities.In the study area, 84.74% of the soil samples were clean in Cd, and 15.26% exceeded the Cd standard limit and were in a polluted state with point-like distributions.Most of the soil samples did not exceed the standard limit of Cd.Regarding Hg, 89.47% of the samples were in a polluted state, among which 5.26% were severely polluted.High pollution points exhibited patchy distributions, indicating the largest influence of human activities.In terms of Pb, 58.95% of the sampling points were categorized as "slightly polluted, " most of which were distributed in the upper reaches of the river.In terms of As, 45.26% of the sampling points were in a clean state, most of which were slightly polluted.The clean sampling points accounted for 56.84% with respect to Cr, with moderate coefficients of variation.This indicated fewer influences from anthropogenic pollution.The average pollution index for the entire area was 2.50, indicating mild pollution.Hg was the biggest contributing factor to the comprehensive pollution level in this area, with a contribution rate of 36.90%.
According to the contribution rates of each element to the slight ecological risk level in the region (Hg, 67.91%; Cd, 15.81%; As, 10.33%; Pb, 4.39%; Cr, 1.56%) in Table 6, Hg was the main factor causing the potential ecological risk in the region.Therefore, relevant departments should pay attention to the potential slight ecological risk that the Hg content may pose to the area.
The spatial distribution of potential ecological risks in the study area is shown in Figure 3.According to the RI value, there were patch-shaped areas with high potential ecological risks, most of which were unevenly distributed.This indicates that the soil was contaminated and enriched by human activities.The RI level analysis shows that the proportion of areas with slight pollution in the entire area was the largest.In contrast, the proportion of areas with very serious potential ecological risks was the smallest.The area with the largest mining area was the most severely contaminated in the entire area, and the spatial distribution showed a trend of gradually increasing potential ecological risks from southwest to northeast.The investigation and analysis also show that the areas with severe pollution were the mining areas and regions with developed transportation.Thus, it can be concluded that industrial and mining activities are one of the causes of the spatial distribution of potential ecological risks in this region.

Analysis of heavy metal pollution sources 4.2.1. Pearson correlation matrix
The correlation analysis in Table 7 shows that stronger correlations indicated higher similarity of pollution sources.The correlation between Cd and other heavy metals was insignificant, suggesting that Cd may result from other pollution sources.The correlation results of this study were similar to those of previous studies on the sources of heavy metals in typical farmland soils in Huibu Town, Changshan County, West Zhejiang.There were two types of pollution sources, both due to frequent industrial and mining activities in the previous and present study areas.However, the present study area is distinctive because it has many roads and railways, indicating that the pollution sources were more likely to come from industrial and mining activities.

PCA-APCS-MLR
The Hongfeng Lake watershed is located in Guizhou Province, which is rich in mineral resources.The surrounding area of the Hongfeng Lake watershed is primarily the alumina ore mining development zone in the Qiannan region (37).Many industrial and mining enterprises release heavy metals into the environment through exhaust gas, wastewater and waste residues.Mineral exploitation can lead to the accumulation and pollution of heavy metal elements in the soil (38,39).Heavy metal pollution has been identified in the sediments of the Hongfeng Lake watershed, with the pollution level of each element as follows: Hg > Pb > As, and the potential ecological risk pollution level was moderate.The main pollution source for the primary pollutant, Hg, is the wastewater pollution generated by chemical plants and industrial enterprises around the lake (40).The farmland soil in the watershed is part of the aquatic ecosystem, and the farmland soil conditions interact with water quality.Through agricultural activities such as ecological migration, planting and irrigation, the deposition and pollution of heavy metals can occur in the farmland soil.Therefore, based on previous research results, the exploration results of the sources of heavy metal pollution in the surrounding farmland of the watershed are as follows.The principal component analysis (Table 11) shows that the cumulative variance contribution rate of PC1 was 50.402%.The variables with the highest factor loadings on PC1 were Hg, As, Pb, and Cr, with loadings of 0.802, 0.842, 0.730, and 0.790, respectively.The evaluation and analysis of heavy metal pollution show that the Hg, As, and Pb contents were higher than the background values of soils.Hg had a slight cumulative pollution, and they all had medium-high variability.This indicates that these heavy metals were greatly affected by human activities and had similar pollution sources.Moreover, the pollution distribution map (Figure 2) of individual elements shows that these three heavy metals had similar local high-pollution distributions, mostly located in areas with intensive mining activities.Previous studies have shown that mining activities in mining areas may aggravate the contamination of heavy metals such as Hg, As and Pb in the region.During mineral mining, rock blasting, grinding and ore beneficiation can cause the migration of As and Pb fine powders to surrounding soils, significantly increasing the As and Pb contents in the soil (35,41,42).Artisanal smelting and tailings are the main Hg sources.These metal elements pollute the soil through atmospheric deposition, weathering and surface runoff, and exhaust emissions from trucks transporting minerals contribute to the Pb pollution in farmland soils (3,43,44).Cr did not exceed the background value, with a low coefficient of variation, and no cumulative pollution occurred.Thus, Cr was less affected by human activities.Combined with previous studies, Cr may be influenced by natural factors such as parent materials and bedrocks.Therefore, the pollution sources of the first principal component include industrial and mining activities and natural sources.
The cumulative contribution rate of PC2 was 72.528%.The main loading factor was Cd, with a coefficient of 0.973.Based on the analysis of heavy metal contents and pollution status, the Cd content in the study area did not exceed the standard limit.This indicated a clean status.The coefficient of variation was 0.61, indicating mediumhigh variations and more significant influences of anthropogenic factors.Studies have also shown that fertilizers contain a large amount of Cd, and pesticides and fertilizers can cause Cd pollution (45,46).This indicates that the pollution sources of the second principal component are agricultural sources.
The contribution rates of heavy metals based on the APCS-MLR model combined with principal component analysis are shown in Figure 4. Factor 1 was attributed to industrial and mining activities and parent material pollution, including Hg, As, Pb, and Cr, with contribution rates of 71.99, 51.57, 67.39, and 68.36%, respectively.Factor 2 was attributed to agricultural pollution, mainly due to Cd, with a contribution of 84.12%.It should be noted that As and Cr accounted for a relatively large proportion in Factor 2, with contribution rates of 43.11 and 27.31%, respectively.Factor 3 was attributed to unknown sources of heavy metal pollution, including Cd (4.95%), Hg (22.98%),As (5.32%), Pb (23.73%), and Cr (4.32%).

Health risk of heavy metals in locals
Tables 9, 10 show that the HI and HQ values for adults and children were both less than 1, and the TCR and CR values were both below 1.00E-04.This indicates that the health risks posed by heavy metals to the human body were within an acceptable range.Among the non-carcinogenic risks, the mean HQ values of As for adults (5.06E-02) and children (9.02E-02) were the highest.In contrast, the mean HQ values of Cd for adults (5.14E-04) and children (4.95E-04) were the lowest.In terms of carcinogenic risk, the mean CR values of Cr for adults (3.87E-05) and children (5.07E-05) were the highest.In contrast, the mean CR values of Cd for adults (1.57E-06) and children (2.78E-06) were the lowest.Therefore, As was the main contributor to non-carcinogenic risks in the human body, and Cr was the main contributor to carcinogenic risks.Cd posed the lowest health risk to the human body.It is particularly noteworthy that the maximum values of As (1.23E-04) and Cr (1.07E-04) in carcinogenic risks exceeded 1.00E-04, which may pose significant carcinogenic risks to the human body.This indicates that the human health risks in the study area are generally acceptable.A few high-value points that pose significant carcinogenic risks cannot be ignored.The average amount of heavy metals ingested by the human body through different exposure pathways is shown in Figure 5.
Figure 5 shows that the amount of heavy metals in each ingestion pathway is in the following order: oral ingestion > skin absorption > inhalation.Oral ingestion is the main pathway for heavy metals to pose health risks to the human body, consistent with previous research.Children have certain physiological and behavioral differences compared to adults, such as weaker detoxification ability, finger-sucking behavior, greater skin exposure area and higher vulnerability (30)(31)(32)(33)(34)47).Therefore, we should pay attention to children's exposure pathways to heavy metal.Health education and publicity should be strengthened to enhance their hygiene awareness.The total non-carcinogenic risk value (HI) in the study area was less than 1.The total carcinogenic risk value (TCR) was between 1.00E-06 and 1.00E-04.This indicates that the health risks of soil heavy metals to the human body in this area are within an acceptable range.However, potential health hazards posed by the primary elements causing health risks should be prevented.

Conclusion
The development of the mining industry and agricultural activities in the upstream of Hongfeng Lake aggravated the heavy metal pollution in agricultural soil.The development of the mining industry resulted in high Hg accumulation in the soil.Multiple high pollution points were distributed in areas with the mining industry and road traffic.The overall pollution level of the study area was classified as mild, with Hg being the primary pollutant.As and Pb were secondary pollutants and were classified as slight pollution.Cd and Cr were in a clean state, without pollution.The overall potential ecological risk level in the entire area was considered slight.Hg was the main factor causing potential ecological risks in this region.
Descriptive data shows that the coefficients of variation for Cd, Hg and As ranged between 0.5 and 1, indicating a high degree of variability.They were more greatly influenced by human activities.The pollution source analysis using PCA combined with the APCS-MLR model reveals that the main pollution sources for Hg, As, Pb, and Cr were related to industrial and mining activities and natural sources.The main source of Cd pollution was agricultural activities.Agricultural sources also contributed significantly to As and Cr.Therefore, in future development, attention should be paid to the combined water and soil pollution.
Heavy metals did not pose a non-carcinogenic risk to human health but had an acceptable carcinogenic risk.As and Cr were the main non-carcinogenic and carcinogenic factors, respectively.Children had a higher hazard index (HI) and target cancer risk (TCR) values than adults due to physiological and behavioral differences.For example, children are more susceptible to the health risks of heavy metals.This is because they tend to suck their fingers, come into contact with substances and have lower detoxification abilities.

Limitations and future perspectives
This study highlights the impact of heavy metal contents in agricultural soils on human health, with an emphasis on cropland soils.It does not cover all soil types.For a more in-depth exploration, other soil types could be included in future research.There are multiple variables involved in heavy metal pollution sources and human health risk assessment.For example, there is a lack of detailed pollution inventories for comparison of pollution sources.Human health risk studies are based on exposure levels to heavy metals in soils.The subjectivity of individuals is overlooked.These can limit the variability of the study.Exploring pollution sources and research methods for human health risk assessment will also involve studying more advanced and timely approaches.In future studies, other research methods will be incorporated for comparison to better support the conclusions and provide a more comprehensive understanding.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.Health risk values under different pathways.

FIGURE 1
FIGURE 1Schematic diagram of sampling points.

FIGURE 2
FIGURE 2Single factor pollution distribution maps of heavy metals.

FIGURE 3 Maps
FIGURE 3Maps of potential ecological risks.

FIGURE 4
FIGURE 4Contributions of heavy metal pollution sources.

TABLE 1
Heavy metal evaluation index.

TABLE 3
RfD and SF values of heavy metals under different exposure pathways.
Considering the actual situation of the area and relevant USEPA regulations, the health risk model evaluation classification values are shown in Table4.

TABLE 2
Parameters for soil heavy metal exposure risk.The ArcMap 10.8 software was used to draw a schematic diagram of the sampling points and a distribution map of heavy metal pollution.Origin 2019 was used to illustrate the relevant statistical results.

TABLE 4
Health risk model evaluation criteria.

TABLE 5
Descriptive statistics of soil heavy metals.

TABLE 6
Calculated pollution indices due to heavy metal concentrations.

TABLE 7
Correlation analysis of heavy metals.

TABLE 8
Linear regression models.

TABLE 9
Non-carcinogenic risks to humans from soil heavy metals.

TABLE 11
Rotational component matrix of heavy metals.