Spatiotemporal distribution patterns and risk characteristics of heavy metal pollutants in the soil of lead–zinc mines

The current soil environmental assessment system is inadequate in terms of the spatiotemporal distribution of heavy metal pollutants. This study employed the numerical simulation technique to predict spatiotemporal distribution patterns of heavy metals within 50 days and to assess the soil risk characteristics of heavy metal pollution near a lead–zinc mine in Hunan Province, China. The spatiotemporal distribution results indicate that the soil in the sewage plant and mining areas served as the pollution center, exhibiting a ladder-shaped pollution diffusion trend outward. When the pollution migration time reached 20 days, pollutant migration and changes tended to remain stable, high-pollution areas exhibited no drastic changes within 10 m, and low-pollution and medium-pollution areas revealed obvious changes. Moreover, the low-pollution area width approached 2 m, the depth reached 2 m, the medium-pollution area width was close to 2.5 m, and the depth approached 4 m. The percentage of areas containing lead–zinc mine soil with high to extremely high risks reached 82.88%, and extremely high-risk farmland, mining and residential areas accounted for up to 100%, 95% and 90%, respectively, of the total area. Among the pollution sources, high-risk and extremely high-risk areas in regard to heavy metal Cd accounted for 13.51 and 49.55%, respectively, of the total area. This study provides new insights into the migration patterns and risk characteristics of pollutants to address soil environmental assessment system problems.


Background
After heavy metal pollutants originating from mining activities enter soil, these heavy metals are continuously accumulated in the soil through the processes of physical and mechanical absorption, retention, colloidal physical and chemical adsorption, chemical precipitation and biological absorption [1]. When the soil capacity is exceeded, the soil composition, structure, nature and function could deteriorate [2], but heavy metals cannot be decomposed by soil microorganisms and are easily enriched in organisms [3]. Mine-contaminated land mainly involves heavy metal pollution, and heavy metal pollution evaluation encompasses a complex process and represents a worldwide problem [4]. In order to reliably evaluate heavy metal pollution various aspects need to be considered. Contaminated area analysis, pollution source determination, and spatiotemporal distribution evaluation of pollutants are essential for soil environmental assessment [5]. Soil environmental assessment can determine the risk level in the contaminated area of a given mine and thus identify pollution sources, which facilitates soil ecological remediation [6]. Heavy metal pollution can directly threaten the lives and safety of residents Open Access *Correspondence: xiechengyu42@xtu.edu.cn 1 College of Environment and Resources, Xiangtan University, Xiangtan 411105, China Full list of author information is available at the end of the article in the functional area. Hence, mine soil environmental assessment requires more attention [7].
Environmental assessment of the soil environment reflects the extent to which heavy metal pollutants are harmful to the human body, and long-term exposure to heavy metal enrichment could lead to neurological disorders and kidney failure [8]. Mass heavy metal poisoning among children in Africa has been attributed to exposure to contaminated soil originating from mines, inhalation of contaminated soil particles, and ingestion of crops grown on contaminated mine soil [9]. An environmental assessment study of a gold mine in Ghana, Africa, confirmed the above claim [10], identifying potential health risks to the local population under the influence of carcinogenic mining practices and contamination of groundwater, agricultural soil and plants, thus indicating the importance of soil environmental assessment for human health. Environmental assessment studies of soils on a global scale can be classified in three categories. The first category involves index equation-based evaluation [11][12][13], in which an equation is employed to project current soil heavy metal pollution conditions, thereby considering the heavy metal content, pollution differences among multiple elements and corresponding risk factors in result calculation, but this approach cannot assess the uncertainty in the prediction results. A typical soil evaluation index is the environmental index [14,15], with objective evaluation results, which can accurately distinguish and quantify the soil risk and identify pollution sources and areas. However, this approach exhibits the disadvantage of not reflecting the spatiotemporal changes in heavy metal pollution. The second category encompasses matrix function evaluation [16][17][18], based on establishing an affiliation function and constructing a fuzzy relationship matrix, thereby determining the classical domain, nodal domain and coefficients to distinguish and quantify the degree of soil heavy metal pollution. However, there exist uncertainties in the factors, coefficients and matrices, the evaluation framework is thus difficult to establish, and the evaluation results are unreasonable. The third category entails geographic information evaluation [19][20][21], which considers spatial database information and employs remote sensing technology and geostatistics for grid establishment, resulting in training sets and spatial graphics. However, the assessment results are inaccurate, and the authenticity of the results cannot be verified. In addition, the numerical simulation technique, as a new type of mathematical model, can spatially characterize soil heavy metal pollution and incorporate temporal parameters to visualize the obtained results [22]. In the numerical simulation process of contaminated sites, Zhou et al. [23] found that pollutants were transported in soil along the direction of groundwater flow and concluded that numerical simulation could be applied to soil pollutants. Thereafter, Cao et al. [24] explored spatiotemporal distribution patterns of heavy metal pollutants in lead-zinc mine soil through numerical simulations and verified the applicability of the numerical simulation technique to heavy metal pollutants in soil.
Soil environmental assessment can reflect the degree of heavy metal damage to the environment, but cannot reflect the spatiotemporal variation characteristics of heavy metal pollution. Numerical simulations can capture soil heavy metal pollution in both two-and threedimensional space, accurately visualize the obtained results and suitably reveal the spatiotemporal distribution of heavy metals. The application of numerical simulation technology could overcome the shortcomings of soil environmental assessment and yield more comprehensive and reliable assessment results. Therefore, this study applied the numerical simulation technique combined with environmental assessment analysis, thereby providing new insights into mine soil evaluation and providing notable support for the prevention and control of heavy metal pollutants in mine soil, in addition to risk warning.

Sampling site
The study area is located in a lead-zinc mine within Hunan Province, China. The mine site exhibits a rugged terrain and variable geological structure, and many landslide disasters have occurred. The mine area is mainly dominated by red loam, followed by sandy loam (Fig. 1). The total mine area consists of six functional areas: mining area, tailings area, residential area, reclamation area, farmland area, and sewage plant area. The mining area is the functional area for ore mining and transportation [25], the tailings area is the functional area for waste ore or tailings stockpiling [26], the residential area is the functional area where local villagers live and work [27], the reclamation area is the functional area for soil ecological restoration and phytoremediation measures by the local government [28], the farmland area is the functional area for local villagers to conduct agricultural cultivation [29], and the sewage plant area is the functional area for the treatment of local mining and metallurgical wastewater and domestic sewage [30]. The present study focuses on profiling the typical heavy metals of Pb, Zn, and Cd within the soil of this lead-zinc mine.

Sampling and experimental methods
To comprehensively assess soil heavy metal pollution in the different functional areas of the mine, 20 soil sampling points were established in the reclamation area (113°04 ′ At each soil sampling point, 2.1 kg of soil was collected at a soil depth of 40 cm, and three replicates were obtained at each sampling point, which were mixed into one sample, placed in a sealed bag, marked, and transported to the laboratory for processing [31]. Considering the ruggedness of the mine roads and actual conditions, the sampling points were arranged via a combination of the S-based and random sampling methods for soil collection purposes [32]. The soil collected at the sampling sites was left outdoors to dry naturally for 7 days, and large sand grains and gravel particles were then separated with a 10-mesh sieve [33]. Individual samples weighed 140 g, which were dried at 104 °C for 12 h, ground, passed through a 100-mesh sieve, and thoroughly mixed in a blender [34]. The total heavy metal content was determined after the quadratic discard method was applied.

Heavy metal content assays
First, 0.5 g of mixed soil was removed via quadratic individual sampling, and the mixed soil sample was placed in a digestion tube. Dissolution of the soil samples was performed with a graphite digestion instrument. A blank group was set up during soil sample digestion, with the aim of reducing the experimental error. Certified materials (HJ/T166) were employed. The digestion process was as follows: 0.5 g of the soil sample was placed in the digestion tube, which was inserted into the digestion cavity of the graphite digestion instrument, after which a few drops of water were added for moistening. Then, 10 mL of concentrated HCl was added, and the controller was set to a low temperature (80-100 °C) for heating and evaporation until approximately 5 mL remained. Subsequently, 15 mL of HNO 3 was added, heat (100-120 °C) was applied until the sample became viscous, 10 mL of HF was added, heat (120 °C) was applied, and the ablation tube was shaken. Then, 5 mL of HClO 4 was added, the sample was heated (130 °C) until white smoke was emitted, and the final decomposition product was white or light yellow and viscous. After cooling, the inner wall of the tube and cap were rinsed with a dilute HNO 3 solution, the residue was dissolved under warming, and the volume was fixed to 50 mL after cooling [35,36]. The five-point calibration method involving standard samples was implemented before heavy metals Pb, Zn and Cd were determined via flame atomic absorption spectrophotometry [37]. The slopes of the calibration curves for Pb, Zn and Cd were 0.9998, 0.9999 and 0.9997, respectively. This sample measurement approach was highly accurate and conformed to the basic standards of the JJD 694-2009 National Metrology Calibration Protocol of the People's Republic of China.

Numerical simulation technology
COMSOL Multiphysics and MATLAB are highly advanced numerical simulation software packages. This study relied on COMSOL Multiphysics and MATLAB software to obtain numerical simulation results. Migration of heavy metal pollutants considers that the process should initially start as point pollution, spread through infiltration into surface pollution and finally reach the scope of pollution [38]. The pollution source in the numerical simulations is set to (15, − 2), and the boundary conditions of the contaminant migration behavior in soil are defined as boundaries (1), (2), and (3), of which boundary (1) is a fixed boundary, boundary (2) is the entrance, and boundary (3) is the exit.
The Richards equation [39] and governing equation for solute transport [40] were adopted to simulate the migration process of heavy metal pollutants via numerical  1) and (2), the migration changes in heavy metal pollutants in soil could be predicted via simulations: where B is the specific moisture capacity, Aa is the effective saturation, A is the storage coefficient, C p is the pressure head, t is the time, D is the hydraulic conductivity, and E is the direction: where A is the concentration, p a is the bulk density, C a is the weight, θ is the porosity, D a is the hydrodynamic dispersion tensor, C a is the reaction in water, C b is the reaction solute, and E a is the solute per unit time.

Environmental risk assessment methodology
Rapant et al. [41] proposed the environmental risk index method to characterize the environmental risk in contaminated environments, which specifies corresponding environmental risk classification criteria and can quantitatively measure the environmental risk level of soil contaminated with heavy metals. Based on Eqs. (3) and (4), the environmental risk of heavy metal pollutants in soil can be calculated.
where I ERi is the environmental risk index of the ith element at the supercritical limit, AC i is the analytical content of the ith element, RC i is the critical limit of the ith element, and I ER is the environmental risk of the test sample.

Environmental risk assessment
The mean value of the environmental risk index for the three heavy metals in the sampled mine soil exhibited the order of Pb (7. 57) > Cd (6. 60) > Zn (5. 51) > 5, and the I ER value reached 19.68 > 5. Hence, each of the three heavy metals posed a very high risk, indicating that these three heavy metals generated severe compound pollution, and during soil remediation, joint treatment of the above three heavy metals should be considered (Fig. 2). The tested mine soil as a whole occurred at a very high environmental risk level, indicating that the enrichment impact on the surrounding soil after years of high-intensity mineral resource remining and smelting activities in the mine was significant and required soil pollution prevention and control as soon as possible [42].
The mean soil environmental risk index for each functional area of the mine was as follows: mining area (49.24) > sewage plant area (18.65) > tailings area (17.94) > farmland area (15.46) > residential area (10.23) > 5 > reclamation area (4.69) > 3 (Table 1). These results demonstrate that the mining, sewage plant, tailings, farmland, and residential areas all occurred at the extremely high environmental risk level. Among the various functional areas of the mine, the environmental risk index for the mining area was several times that for the other areas, indirectly reflecting that the heavy metal pollution level in the soil of the mining area was extremely high, which is directly related to the mining activities in the mining area, including unregulated ore collection, smelting, rock weathering and reckless accumulation of waste rock tailings. This could lead to direct infiltration of heavy metal pollutants into the soil under migration [43]. The reclamation area exhibited a high environmental risk, which indicates that the local government should achieve a better ecological restoration effect in land reclamation and restoration.
In terms of the soil environmental risk level, the proportions of the risk-free, low-risk, medium-risk, highrisk, and extremely high-risk mine areas reached 8.11%, 4.50%, 4.50%, 6.31%, and 76.58%, respectively, of the total area. Areas with a high risk and above accounted for 82.88% of the total area, which constitutes a notable environmental warning signal to surrounding residents, and soil remediation of the mine site should be carried out as soon as possible (Table 2). Extremely highrisk farmland areas occupied 100% of the soil in the (5) Mean = X 1 + X 2 + · · · + X n n .  different functional areas. This phenomenon is closely related to the nature of soil in these farmland areas. The soil in these farmland areas belongs to clay soil with a high water content, in which heavy metal elements can easily accumulate [44]. In addition, under rainwater runoff, many pollutants are transported from the top of hills and subsequently accumulate in farmland areas [45]. Extremely high-risk mining areas accounted for 95%, extremely high-risk residential areas occupied 90%, extremely high-risk tailings areas accounted for 85%, extremely high-risk sewage plant areas occupied 65%, and extremely high-risk reclamation areas accounted for 35% of the total area, which indicates that the impact of the different production activities on heavy metal accumulation in soil was differentiated and priority should be given to the treatment and maintenance of farmland areas [46]. Prevention and control of soil environmental pollution should first determine the risk source areas and heavy metal risk sources for priority control determination. According to physical examination, most residents living near lead-zinc mines suffer certain diseases, which are directly related to heavy metals Pb, Zn and Cd in lead-zinc mine soil. Heavy metals Pb, Zn and Cd in the soil of the mining area indicate extremely high pollution levels, with risks reaching as high as 95, 80 and 80%, respectively, thus confirming that mining activities in the mining area can cause serious compound pollution in the soil environment within the region, which is highly difficult to manage (Table 3). Therefore, the mining area  should be prioritized as a risk source area [47]. The soil environmental risk levels in the sewage plant area, tailings area, farmland area, residential area, and reclamation area were highly variable, indicating that the soil in the mining area suffered great damage due to human factors and was contaminated to notably different degrees based on the different areas. Consequently, an uneven soil pollution distribution was observed [48]. The highrisk and extremely high-risk heavy metal pollution source areas in regard to the occurrence of Pb, Zn, and Cd in the mine soil accounted for 12.61 and 36.04%, 27.93 and 35.14%, and 13.51 and 49.55%, respectively, of the total area. Via comparison of the sources of these three heavy metal pollutants, it could be concluded that Cd contamination of the soil in the mining area was the most serious and that the study area was contaminated by Cd with the highest soil environmental risk coefficient. Moreover, the contaminated area was very extensive, the degree of contamination was high, and the toxic risk to humans occurred at a high level [49]. Therefore, according to the ranking of the soil pollution risk sources, i.e., Cd > Zn > Pb, heavy metal Cd should be considered the highest-priority soil pollution risk source for control in the study area.

Environmental risk spatial distribution prediction model
The spatial distribution of pollutants is an important method to reveal the risk of soil environmental pollution. There were mainly three high-pollution risk areas distributed in the north and east. The mine pollution distribution was concentrated, and the pollution trend basically exhibited an 8-type pattern (Fig. 3). The soil in the northern sewage treatment plant area and eastern mining area was highly polluted, and treatment and maintenance measures should be implemented in advance to prevent pollution migration from impacting the surrounding soil after the emergence of high soil environmental risks. The overall pollution situation demonstrated that the mining area exhibited a peak, there occurred abnormally high values, and the surrounding distribution revealed a gradually decreasing trend, indicating that production activities affected the accumulation of heavy metals in the mining area. Hence, it is necessary to rectify and plan human activities and the development and utilization of mineral resources in the mining area [50]. The  Fig. 3 Three-dimensional map of the environmental risk levels of the heavy metal pollutants in the different areas of the mine sewage plant area, mining area and tailings area pollution distributions are closely linked, the heavy metal pollution distribution within the region is relatively concentrated, and there exist differences in the pollution distribution. It could be found that these three areas served as the pollution center, while the surrounding pollution level exhibited a gradual weakening trend, which could lead to pollution penetration and diffusion in the different functional areas within the mine. This indicates that the pollution sources between these three areas are consistent and significantly correlated, which is related to early unreasonable mineral resource development, tailings piling and other highly correlated anthropogenic activities [51]. In addition, there occurs mutual contamination between mining, tailing and sewage disposal activities within the mine, thus altering the status of soil contamination between the different regions.
The soil in the sewage plant area exhibited two highly polluted blocks, while the center of pollution coincided with the center of the area, indicating wave-like outward extension. Most of the surrounding soil exhibited a step-type pollution distribution. On the one hand, this could be attributed to the unregulated discharge of smelting wastewater by sewage plants [52], while on the other hand, this could also be attributed to the improperly managed sewage plant drainage system [53], which indicates that the sewage plant severely polluted the surrounding soil, and sewage treatment measures should be implemented to prevent heavy metal pollution upon migration through water after diffusion. There were similar contamination distributions between the reclaimed and residential areas in the southern part of the mine, with a moderate risk of contamination, and these two areas are adjacent to the mining and tailings areas, thus confirming that the mining and tailings areas of the mine exerted a fundamental impact on the soil in adjacent areas. Therefore, effective soil contamination prevention and control of pollution sources in the mining and tailings areas could eliminate the associated pollution attributed to soil pollution migration and reduce the risk of pollution in the surrounding functional areas [54]. In addition, the large differences in soil pollution status between the different functional areas suggest that the variability in human activities in the different areas considerably influenced the degree of soil heavy metal accumulation.

Spatiotemporal distribution prediction model
The migration mechanism of heavy metal pollutants within mine soil entails a dynamic process [55], and the spatiotemporal distribution of heavy metal pollutants within soil can be addressed by applying simulation technology with the aim of exploring migration changes in heavy metal pollutants. During the migration process, heavy metal pollutants initially infiltrate into the soil from boundary (2) and continue to diffuse toward boundary (3) due to gravity and soil pore water transport, and the pollution diffusion direction near boundary (1) demonstrates a vertical diffusion trend, which is directly related to the angle of the mine soil slope (Fig. 4). In the vertical state, heavy metal pollutants in soil could appear to accumulate unilaterally, but the pollution infiltration direction remains oriented from the high-pollution area to the low-pollution area during migration [56]. In addition, based on the pollution migration direction between boundaries (2) and (3), it can be concluded that pollutant migration also exhibits the characteristic of diffusion along multiple lateral directions, which can be summarized as diffusion from the polluted side toward the right hemisphere, which is consistent with the conclusion of Jacques et al. [57] that the diffusion direction of heavy metal pollutants could exhibit irregular dispersal movement along all directions with increasing soil depth, but the change trend of pollutant infiltration always exhibits a semielliptical shape during migration. The spatiotemporal distribution of heavy metal pollutants in the mine soil was investigated (Fig. 5), which revealed the following: on day zero, the heavy metal pollutants gathered at the pollution source, the surface width reached 1 m, and no pollution infiltration changes occurred. On the 1st day, the heavy metal pollutants started to migrate, the pollutants gradually infiltrated, the high-pollution area width reached 2 m, the depth reached 3 m, and the pollution sources exhibited the low-pollution phenomenon, which confirms the theory that pollutants can continuously infiltrate into soil [58]. On the 2nd day, heavy metal pollutant migration drastically changed, with a high-pollution area width of 4 m and a depth of 4 m. On the 5th day, the migration rate of the heavy metal pollutants continuously increased, and the width of the highly contaminated area reached 2.5 m. However, if the width of the boundary was not set in the model, the width of this contaminated area could exceed the limit, and the depth could approach 6 m. On the 10th day, the migration capacity of the heavy metal pollutants reached its peak, the depth of the highly polluted area already exceeded 10 m, and the pollution source area gradually became a low-risk area. On the 15th day, with continuous downward movement of pollutants, the low-pollution area also expanded, with a width close to 1.5 m and a depth close to 1.5 m, indicating that if the subsequent environmental pollution could be controlled, the surface layer would not further expand. On the 20th day, the pollutant migration trend reached the steady state, no drastic change in the high-pollution area occurred within 10 m, obvious changes were observed in the low-and medium-pollution areas, the width of the lowpollution area approached 2 m, the depth was close to 2 m, the width of the medium-pollution area almost reached 2.5 m, and the depth approached 4 m. The 25th day was similar to the 20th day, which suggests that the spatiotemporal distribution of the pollutants within 10 m already approached the stable state and could not further produce large fluctuations without subsequent pollution [59]. To verify the validity of this statement, a 50 day numerical simulation experiment was carried out, and the final verification results were consistent. The migration trends on the 50th day were the same as those on the 25th day, confirming that the migration process of heavy metal pollutants within the 10 m depth approached stability, and the spatiotemporal distribution of the three pollution levels of high, medium and low risks did not significantly differ from that on the 20th day.
At the temporal level, the transport process of heavy metal pollutants in soil was enhanced with increasing number of days, from day zero to the 5th day. The migration trend of pollutants was very obvious, spreading outward in a semicircular wave, and a high soil pollution level always occurred within the depth range of 4 m, which indicates that the transport capacity of pollutants greatly fluctuated at the initial stage. Moreover, if the pollution sources could be controlled during this period, the pollution infiltration area could be greatly reduced. From the 10th day to the 15th day, the change trend of the pollutants gradually weakened, and the pollution depth exceeded 10 m, but the low-pollution area also continued to expand during pollutant infiltration. A pollution accumulation surface was often generated in the deep layers, and only a small fraction accumulated in surface soil [60]. This phenomenon indicates that groundwater could also suffer from serious heavy metal pollution. From the 15th day to the 20th day, the pollutant migration change trend approached stability, there occurred no drastic response fluctuations in the high-pollution area, and the distribution of low-and medium-pollution areas expanded to the 3.5 m depth. Therefore, 20 days could be adopted as a time node for pollution migration, and pollution sources should be identified and then mitigated within 20 days during pollution source control to avoid greater environmental pollution. From the 20th day to the 25th day, the migration change in pollutants was similar to that on day 20. At this time, two speculations could be made. First, the migration behavior of the pollutants reached the stable stage, and 20 days could be regarded as a turning point for migration capacity stabilization. Second, when pollutant migration exceeds a soil depth of 10 m, pollutants could continue to migrate downward, and pollution changes could continue [61]. There exists a specific period ranging from 20 to 25 days. By studying the period from 25-50 days, validation revealed that when the simulation time was extended by 25 days, the migration behavior of pollutants always exhibited a high similarity to that on day 25. Therefore, day 20 was the turning point in regard to stabilization of the migration capacity, but this did not indicate that pollution infiltration had ceased, and the heavy metal pollutants in soil could still affect a larger area of soil for a considerable period. At the spatial level, there occurred three stages, namely, before, during and after, and the first stage involved a high-risk area with a depth ranging from 0-3 m. The contamination level of heavy metal pollutants in the surface layer was serious, the soil ecological risk was extremely high, and soil control is needed targeting surface soil to prevent local residents from being harmed through the food chain [62]. In the medium term, the depth range from 3 to 10 m represented the high-risk area, and the pollution level in the middle soil layer was higher than that in the surface soil. To more effectively remediate this risk area, plants with a root system larger than 3 m should be selected as the first choice for mine soil remediation [63]. Later, when the high-risk area exceeds the 10 m depth, heavy metal pollution further infiltrates into the deep soil. When pollution reaches the deep soil, groundwater is very easily contaminated, which could lead to serious contamination of surrounding water sources [64]. Therefore, to effectively control pollution, it is necessary to control pollution infiltration in the medium term to prevent pollution expansion. Heavy metal pollutants diffuse from the source to surrounding areas in soil, and the percolation path can continuously change in the soil space and erratically diminish with increasing depth [65], but the soil at the different depth levels could exhibit varying pollutant migration behaviors due to the difference in soil texture and soil pore space [66]. The migration ability of heavy metal pollutants in the range from 4 to 6 m was obviously weakened, and the migration ability and pollution level were also much lower than those in the depth range from 0 to 4 m, which indicates that the infiltration rate of soil pollutants could decrease with increasing depth in the migration process of heavy metal pollutants, thus causing pollutants to gather in the middle soil layer. It could be inferred that when the depth exceeds 10 m, the migration ability of heavy metal pollutants relative to the 0-4 m depth could demonstrate a significant weakening trend, and the pollution level could continue to decrease. This is consistent with the results of Černík et al. [67], who found that the migration capacity of heavy metal pollutants is mainly influenced by convection, diffusion and dispersion, with the ranking of the pollution magnitude in the order of middle soil layer > surface soil layer > deep soil layer.
Migration of heavy metal pollutants involves a continuous seepage process, and heavy metal pollutants in soil can be regarded as solutes, which can continuously seep downward under certain external conditions [68]. The spatiotemporal distribution of soil heavy metals revealed that most pollutants accumulated in the middle soil layer at depths from 3 to 10 m. The migration change ability from 0 to 3 m was notable, but the pollution distribution was narrow, and at depths beyond 10 m, the heavy metal pollutant migration capacity decreased, and the pollution level also gradually decreased with increasing soil depth. It was concluded that most heavy metal pollutants in the soil accumulated at depths from 3 to 10 m, with less contamination at depths from 0 to 3 m, which indicates that the ecological risk level in the middle soil layer could be higher than that in the surface soil layer, and long-root plants should be cultivated to mitigate the environmental pollution encountered in the area [69] and reduce the soil environmental risk.

Conclusion
Improving the soil environment evaluation system is crucial for early warnings and prevention and control of heavy metal pollution problems. This study applied the environmental risk index to analyze pollution sources and pollution areas, proposed the integration of the numerical simulation technique to further explore the spatiotemporal distribution patterns of heavy metal pollutants and revealed the hazard level and associated risk control measures for the soil in each functional area within the lead-zinc mine. The research route differed from previous single-factor environmental assessment studies of soil, but numerical simulation and predictive analysis were employed to further improve the assessment system and comprehensively reveal the pollution risk and spatiotemporal distribution characteristics of lead-zinc mines.
The soil I ER value in the total mine area reached 19.68 > 5, which indicates a very high environmental risk level and requires soil pollution prevention and control as soon as possible. The soil I ER value in the mining area reached 49.24, the heavy metal pollution sources of Pb, Zn and Cd in the functional area indicated extremely high pollution levels, and the mining area should be given priority as a risk source area. The soil pollution risk sources were ranked as Cd > Zn > Pb, and heavy metal Cd should be considered a priority source of soil contamination risks for control in the study area.
There mainly occurred three high-pollution risk areas in the mine, distributed in the north and east. The mine pollution distribution was concentrated, and the pollution trend basically exhibited an 8-type pattern. The overall pollution conditions revealed the mining area as the peak area. There were abnormally high values, and the surrounding distribution exhibited a gradually decreasing trend. Treatment and maintenance measures should be implemented in advance to avoid high soil environmental risks after the migration of contaminants from surrounding soil.
The 20th day was identified as a turning point in regard to stabilization of the migration capacity. In pollution source control, pollution sources should be identified and then mitigated within 20 days to avoid greater environmental pollution, but this does not suggest that pollution infiltration has ceased, and heavy metal pollutants in soil could still affect a large area of soil for a considerable period.
Most heavy metal pollutants in the soil accumulated and were distributed at depths from 3 to 10 m, but the contamination distribution was narrow at depths from 0 to 3 m, which indicates that the ecological risk level in the middle soil layer could be higher than that in the surface soil, and long-root plants should be cultivated to mitigate environmental pollution in this area and reduce the soil environmental risk.
Numerical simulation technology methods are suitable for the prediction of the spatiotemporal distribution of soil heavy metal pollutants. However, soil samples are difficult to collect at depths from 3 to 10 m for quantitative verification of simulation results, so it is necessary to collect soil samples at depths from 3 to 10 m for subsequent scientific research on prediction verification. impact and human health risk assessment of heavy metals in surface soils from selected major mining areas in Ghana. Environ Geochem Health 41 (6)