The Patterns of Migration of Potentially Toxic Elements from Coal Mining Subsidence Areas and Associated Soils to Waterlogged Areas

It is crucial for effectively controlling potentially toxic element (PTE) pollution to understand the pollution situation, ecological risks, health risks, and migration patterns of PTEs. However, currently, no research has been conducted on the migration patterns of soil PTEs from coal mining subsidence areas to waterlogged areas under different restoration modes. In this study, a total of 15 sediment samples and 60 soil samples were collected from landscaped wetlands, aquaculture wetland, fish–photovoltaic complementary wetland, photovoltaic wetland, and waterlogged areas with untreated coal mining subsidence. The PTE pollution status, ecological risks, health risks, migration patterns, and the important factors influencing the migration were analyzed. The results indicated that the comprehensive pollution level of PTEs in waterlogged areas with coal mining subsidence can be reduced by developing them into landscaped wetlands, aquaculture wetlands, fish–photovoltaic complementary wetlands, and photovoltaic wetlands. Additionally, the closer to the waterlogged area, the higher the Cu content in the subsidence area soil is, reaching its peak in the waterlogged area. The Cd was influenced positively by SOC and pH. The research results were of great significance for formulating reclamation plans for waterlogged areas and controlling PTE pollution.


Introduction
Long-term underground coal mining has brought about a series of ecological problems, restraining the sustainable development of society [1].The most prominent among them is the issue of coal-mining-related subsidence and water accumulation [2,3].The scale of high-groundwater-level waterlogged areas with coal mining subsidence in China reaches up to 2000 km 2 .The water bodies in the subsided areas are generally closed or semi-closed systems with poor water circulation and limited self-purification capacity, resulting in the disruption of surface water systems and degradation of water quality [4].Contaminated subsidence water further enters the agricultural system through activities like irrigation [5].The subsidence water has a high content of potentially toxic elements (PTEs), including metals with a density greater than 4.5 g cm −3 , such as cadmium (Cd), lead (Pb), chromium (Cr), copper (Cu), zinc (Zn), and mercury (Hg), as well as arsenic (As) and other metalloids [6,7].PTEs in the polluted irrigation water accumulate in the soil or vegetation influenced by the soil's cultivated and physiochemical condition (pH, Eh et al.), severely affecting local food security, ecological safety, and social stability [8][9][10].Restoration is urgently needed for waterlogged areas with coal mining subsidence.
To address the environmental issues in waterlogged areas with coal mining subsidence and make the best use of natural resources, large-scale restoration efforts are currently underway [11,12].This includes activities such as aquaculture development and the creation of landscaped wetlands within the waterlogged areas [13].Furthermore, under the pressure of "carbon neutrality" and "carbon emission reduction", an increasing number of waterlogged areas are being developed as photovoltaic wetlands or fish-photovoltaic complementary wetlands [2,14].These development approaches can enhance the environmental aesthetics of the waterlogged areas and provide economic and energy value [15,16].However, it cannot be overlooked that these approaches can introduce additional environmental pressures on aquatic ecosystems [17].For instance, aquaculture feed is rich in PTEs such as Cu and Zn [18].Purified urban sewage contains PTEs such as Cd and Cr [19].Will the introduction of feed or purified sewage into the stagnant water area exacerbate the environmental pressure on the aquatic ecosystem?Do photovoltaic components disrupt the natural equilibrium of the aquatic ecosystem and potentially diminish its self-purification capacity for PTEs?It is crucial for us to investigate the extent of environmental stress that different restoration methods may impose on the aquatic ecosystem.
Furthermore, the main sources of PTEs in soil are coal mining operations, agricultural nonpoint source pollution, transportation, atmospheric deposition, and so on [20][21][22].As a result, soil PTEs enter the waterlogged areas with coal mining subsidence through surface runoff and agricultural drainage, accumulating in the sediment.Nonpoint source pollution and coal mining pollution exhibit characteristics such as randomness and fuzziness [23], making the study of them challenging.Currently, there is no research conducted on the migration patterns of soil PTEs from coal mining subsidence areas to waterlogged areas.
PTEs have a high toxicity, are difficult to degrade, and can disrupt the ecological balance, posing a threat to human health [24][25][26].Min Tan [27] and Lingli Zhou [28] found that the contents of Cd, Pb, Cr, Cu, Zn, Hg, and As were relatively high in the mining area environment.It is crucial to understand the pollution levels, ecological risks, health risks, and migration patterns of these seven PTEs in waterlogged areas with coal mining subsidence under different restoration modes.This understanding is vital for ensuring the health of residents living nearby.
For the aforementioned reasons, this study, for the first time, investigated the migration patterns and influencing factors of soil PTEs into subsidence waterlogged areas in a highgroundwater-level coal mining area under different remediation modes.The research aimed to (1) analyze the PTE pollution, ecological risks, and health risks in waterlogged water areas with coal mining subsidence under different remediation modes, (2) investigate the migration patterns of soil PTEs from the surrounding areas of coal mining subsidence to waterlogged areas, and (3) identify the key factors influencing the migration of PTEs.The research findings can provide data support for the management of PTEs and aid in the selection of future management decisions in subsidence areas.

Sample Collection and Analysis
The study area is located in the high-water-table region in eastern China.It has a temperate semi-humid monsoon climate with an average annual temperature of 14.2 • C, an average annual sunshine duration of 2307.9 h, and an average annual precipitation of 816.4 mm.The area is crisscrossed with rivers, and water resources are relatively abundant.The terrain is higher in the southwest and lower in the northeast, forming a typical alluvial plain.The soil is developed from yellow alluvial sediments.The research area has eight mines and is the only coal production base in Jiangsu Province.The proven coal reserves are 2.37 billion tons.Since 1977, a total of 230 million tons of raw coal has been extracted, resulting in significant land subsidence.

Assessment of Trace Metal Pollution in Sediment
The single-factor index is a measure of the pollution level of soil PTEs, which pares the measured data with standard data [33].The specific calculation formul follows: where CF is the single pollution index for PTE i, Ci is the measured value of PTE i, is the standard value of PTE i.The pollution level was classified based on the magn of the CF value: CF ≤ 1, no pollution; 1 < CF ≤ 2, mild pollution; 2 < CF ≤ 3, mo pollution; 3 < CF, severe pollution.The Pollution Load Index (PLI) is used to evaluate the environmental pollution sure of soil PTEs in the study area [34].When PLI is less than 0, it is taken as 0. Wh ≤ 1, it indicates the absence of environmental pollution.When PLI > 1, it indicat presence of environmental pollution.The specific calculation formula for PLI is as fo Three sample plots were set up in each of the five waterlogged areas, and three sediment samples at a depth of 0-20 cm were randomly collected in each plot.Additionally, surface soil samples (0-20 cm depth) were collected at distances of 5 m, 100 m, 200 m, and 300 m from each sediment plot in the farmland.All samples were cleaned of stones, gravel, and dead leaves, mixed thoroughly, and then divided into two sterile sealed bags.One bag was used for determining the physicochemical indicators, including PTEs, and the other bag was kept as a backup.A total of 15 sediment samples and 60 soil samples were collected in this study.All samples were transported to the laboratory on the day of collection and stored in a dark and cool environment.

Assessment of Trace Metal Pollution in Sediment
The single-factor index is a measure of the pollution level of soil PTEs, which compares the measured data with standard data [33].The specific calculation formula is as follows: where CF is the single pollution index for PTE i, C i is the measured value of PTE i, and Si is the standard value of PTE i.The pollution level was classified based on the magnitude of the CF value: CF ≤ 1, no pollution; 1 < CF ≤ 2, mild pollution; 2 < CF ≤ 3, moderate pollution; 3 < CF, severe pollution.The Pollution Load Index (PLI) is used to evaluate the environmental pollution pressure of soil PTEs in the study area [34].When PLI is less than 0, it is taken as 0. When PLI ≤ 1, it indicates the absence of environmental pollution.When PLI > 1, it indicates the presence of environmental pollution.The specific calculation formula for PLI is as follows: where CF i represents the single-factor pollution index of the PTE.

Sediment Quality Guidelines (SQGs)
The sediment quality guidelines (SQGs) establish threshold values for specific chemical substances in sediments that do not pose a risk of harm to benthic organisms or other relevant water body functions [35].Referring to the guidelines proposed by Macdonald et al. [35], the threshold effects level (TEL), probable effects level (PEL), and the biological effects induced by PTEs in this study were evaluated.According to the values, they were divided into three groups: rarely (below TEL), occasionally (between TEL and PEL), and frequently (above PEL) [36,37].

Ecological Risk Assessment
The potential ecological risk index (RI) is an ecological risk assessment method established by Swedish scientist Hakanson in 1980 [33].This method takes into account factors such as PTE content and level of toxic pollution.It has been widely applied in the assessment of environmental pollution and ecological risk in sediments, soils, and other areas [33,38].The specific calculation formula is as follows.
In the equation, E i r represents the potential ecological risk factor, and T i t represents the toxicity response coefficient of PTE i.The coefficients for Cu, Zn, Cr, Cd, Hg, As, and Pb were 5, 1, 2, 30, 40, 10, and 5, respectively [39].CF i represents the single-factor pollution index of the PTE, and RI represents the potential ecological risk index.The pollution levels corresponding to the E i r and RI are shown in Table 1.

Human Health Risk Assessment
The health risks of PTEs to the human body can be broadly categorized into two major types: carcinogenic health risks and noncarcinogenic health risks [40].The most widely used model for assessing health risks associated with PTEs is the health risk assessment model recommended by the U.S. Environmental Protection Agency (USEPA).The pathways of exposure to PTEs in the environment include ingestion through food and water, dermal contact, and inhalation through respiration.
The health hazards of PTEs to the human body depend on factors such as duration of exposure, route of exposure, and age of the exposed individual [41][42][43].Based on the methodology proposed by the Risk Assessment Information System (RAIS) in the United States, the health risks for children and adults were assessed.The average daily intake (ADI) (mg kg −1 d −1 ) was calculated using the following formula: where c is the PTE content, mg kg −1 , and EF represents the annual average exposure days, d a −1 .ED represents the exposure years, calculated separately for children and adults as 6a and 26a, respectively.IngR is the ingestion rate, calculated separately for children and adults as 200 mg d −1 and 100 mg d −1 , respectively.RBA represents the relative bioavailability, with an As value of 0.6 and the remaining PTEs calculated as 1.0 [44,45].AT represents the number of days, d; BW is the average body weight, calculated separately for children and adults as 15 kg and 70 kg, respectively.The hazard quotient (HQ) and the hazard index (HI) for PTE exposure in children and adults were calculated using the following formulas.If HQ and HI values were greater than 1, it indicated that PTEs may have adverse effects on health.If HQ and HI values were less than 1, it suggested that the risks of noncarcinogenic effects on humans are minimal.
where RfDi is the reference dose for PTE i, which can be determined based on the US EPA's Risk Information System.The reference doses for Cd, Cr, Cu, Zn, Pb, and As are 0.0003, 1.5, 0.04, 0.3, 0.00014, and 0.001, respectively [46].The sampling map in this study was created using ArcGIS (version 10.2).The data were tested for normality using the Shapiro-Wilk test, and the results indicated that the data followed a normal distribution.Data analysis and data visualization were conducted in Matlab (version R2023a), Sigmaplot (version 14.0), and Microsoft Excel (version 2021).

Differences in PTE Content among the Five Water Areas with Coal Mining Subsidence
The average content of Cu, Cr, Hg, Pb, Zn, Cd, and As in the sediment of the five stagnant water areas are shown in Table 2.The results indicated that among all the stagnant water areas, Zn had the highest average content, followed by Mn, Zn, Cr, Cu, Pb, As, Cd, and Hg.Comparing the five stagnant water areas, untreated coal mining subsidence water area E had higher contents of Hg, Pb, and Cd compared with areas A, B, C, and D. The highest contents of Cu and Zn were found in the aquaculture pond B. We evaluated the PTE pollution in the soil the study area with reference to the average values of chemical composition in the Upper Continental guidelines (UCG) and sediment quality standards, respectively.In addition, it is also useful to assess the level of PTE pollution in a specific area using the background values of PTEs as standards.The soil PTE background values represent the baseline values of the soil environment in that area that is either unaffected or minimally affected by human activities.The distribution of evaluation indices is shown in Figure 2. Based on the average values of the single-factor pollution index, when referring to UCG, the ranking of single-factor pollution indices for the seven studied PTEs was as follows: Cd In terms of the PLI, the values of A, B, and C were all below 1 when referring to the UCG, indicating the absence of environmental pollution.However, the PLI values of D and E were 1.02 and 1.61, respectively, indicating the presence of some degree of overall PTE pollution.When referring to the background values in Jiangsu Province, all sampling points had PLI values higher than 1, indicating the presence of PTE pollution.Among them, sampling point E had the highest PLI value, which was 2. 13.
With reference to the SQGs, it can be observed that for the sediment of Anguo Landscaped Wetland (A), the level of Cu was below the TEL, indicating that adverse biological effects were unlikely to occur.The Cu levels in the sediment of the other four waterlogged areas fell between the TEL and PEL, suggesting that occasional adverse biological effects may occur.
The Cr levels in the sediment of sampling points A, B, and C were below the TEL, which indicated that adverse biological effects were unlikely to occur.However, in the sediment of C and D, the Cr contents were between the TEL and PEL, suggesting that occasional adverse biological effects may occur.In the sediment of all five waterlogged areas, the levels of Hg, Pb, and Zn were below the TEL, indicating that adverse impacts were unlikely to occur.Additionally, the Cd level in all waterlogged areas (except E) was below the TEL, suggesting a low likelihood of adverse effects.As for As, the values in the other sampling zones except for B were above the TEL but below the PEL, indicating that occasional adverse biological effects may occur.

Ecological Risk Assessment
The ecological risk of PTEs was evaluated with reference to the UCG.The In terms of the PLI, the values of A, B, and C were all below 1 when referring to the UCG, indicating the absence of environmental pollution.However, the PLI values of D and E were 1.02 and 1.61, respectively, indicating the presence of some degree of overall PTE pollution.When referring to the background values in Jiangsu Province, all sampling points had PLI values higher than 1, indicating the presence of PTE pollution.Among them, sampling point E had the highest PLI value, which was 2. 13.
With reference to the SQGs, it can be observed that for the sediment of Anguo Landscaped Wetland (A), the level of Cu was below the TEL, indicating that adverse biological effects were unlikely to occur.The Cu levels in the sediment of the other four waterlogged areas fell between the TEL and PEL, suggesting that occasional adverse biological effects may occur.
The Cr levels in the sediment of sampling points A, B, and C were below the TEL, which indicated that adverse biological effects were unlikely to occur.However, in the sediment of C and D, the Cr contents were between the TEL and PEL, suggesting that occasional adverse biological effects may occur.In the sediment of all five waterlogged areas, the levels of Hg, Pb, and Zn were below the TEL, indicating that adverse impacts were unlikely to occur.Additionally, the Cd level in all waterlogged areas (except E) was below the TEL, suggesting a low likelihood of adverse effects.As for As, the values in the other sampling zones except for B were above the TEL but below the PEL, indicating that occasional adverse biological effects may occur.

Ecological Risk Assessment
The ecological risk of PTEs was evaluated with reference to the UCG.The E i r and RI are shown in Figure 3a.The E i r values of Cu, Cr, Hg, Pb, Zn, and As were below 40 in all five waterlogged areas, indicating that they were at a weak ecological risk level.The E i r value of Cd in sampling point A fell between 40 and 80, suggesting that Cd was at a moderate ecological risk level.However, the E i r values of Cd in B, C, and D were 99.51, 123.77, and 124.39, respectively, indicating that Cd was at a high ecological risk level.The E i r value of Cd in E was as high as 346.50, indicating that Cd was at an extremely high level of ecological risk in the waterlogged area E with untreated coal mining subsidence.Specifically, the RI values of sediments A and B were below 150, indicating a weak ecological risk.The RI values of sediments C and D fell between 150 and 300, indicating a moderate ecological risk.The sediment at point E had the highest RI value, reaching 426.44, which indicates a high ecological risk level.

Health Risk Assessment
The HQ and HI values for children and adults for each site are shown in Figure 4.It can be seen that the highest HQ and HI values for children and adults were observed at sampling point S1 in landscaped wetland A. This was because S1 was the outlet of treated wastewater from a sewage treatment plant, resulting in a high PTE content.After the influence of subsurface flow wetlands, lateral flow wetlands, and other factors, the PTE content at the remaining two sampling points, S2 and S3 in landscaped wetland A, was relatively lower.The three sampling points S13, S14, and S15 in the waterlogged area E with untreated coal mining subsidence all had higher HQ and HI values.

Health Risk Assessment
The HQ and HI values for children and adults for each site are shown in Figure 4.It can be seen that the highest HQ and HI values for children and adults were observed at sampling point S1 in landscaped wetland A. This was because S1 was the outlet of treated wastewater from a sewage treatment plant, resulting in a high PTE content.After the influence of subsurface flow wetlands, lateral flow wetlands, and other factors, the PTE content at the remaining two sampling points, S2 and S3 in landscaped wetland A, was relatively lower.The three sampling points S13, S14, and S15 in the waterlogged area E with untreated coal mining subsidence all had higher HQ and HI values.

Health Risk Assessment
The HQ and HI values for children and adults for each site are shown in Figure 4.It can be seen that the highest HQ and HI values for children and adults were observed at sampling point S1 in landscaped wetland A. This was because S1 was the outlet of treated wastewater from a sewage treatment plant, resulting in a high PTE content.After the influence of subsurface flow wetlands, lateral flow wetlands, and other factors, the PTE content at the remaining two sampling points, S2 and S3 in landscaped wetland A, was relatively lower.The three sampling points S13, S14, and S15 in the waterlogged area E with untreated coal mining subsidence all had higher HQ and HI values.Table 3 displays the noncarcinogenic risk assessment results for Cu, Cr, Hg, Pb, Zn, Cd, and As in the sediment of each waterlogged area.The noncarcinogenic risk assessment for children indicated that Pb had a HQ value higher than the threshold value of 1, indicating a certain noncarcinogenic health risk for children.For other PTEs, the HQ values for children were below the threshold value of 1, indicating no noncarcinogenic health risks for children associated with those metals.The HQ values of Cu, Cr, Hg, Pb, Zn, Cd, and As for adults were below the threshold value of 1 in all the sampling sites, indicating no noncarcinogenic health risks for adults associated with these PTEs.However, due to the HQ value of Pb for children being higher than 1, the HI values (excluding C) for all waterlogged areas were also higher than 1 for children.Among them, area E had the highest HI value, indicating that it posed the greatest noncarcinogenic risk for children.

Variation Trends of Soil PTEs at Different Distances
The main source of PTEs in the sediment from waterlogged areas with coal mining subsidence is human activities, including mining activities and the use of pesticides and fertilizers [49].Due to the lower topography of the waterlogged areas with coal mining subsidence, PTEs from surrounding farmland soils enter these areas through surface runoff.In theory, the closer the distance to the waterlogged areas with coal mining subsidence, the higher the soil's PTE content.Therefore, the sediment in waterlogged areas with coal mining subsidence is expected to have a higher content of PTEs compared with the region overall.In this study, the variation trends of the PTE content in soil surrounding the waterlogged areas with coal mining subsidence are shown in Figure 5.It can be found that the contents of Cu, Pb, Zn, and As increase as the distance from the waterlogged areas decreases, which was in line with theoretical expectations.However, there was no clear pattern observed in the variation in the Cr, Hg, and Cd content in relation to the distance in the soil.
coal mining subsidence is expected to have a higher content of PTEs compared with the region overall.In this study, the variation trends of the PTE content in soil surrounding the waterlogged areas with coal mining subsidence are shown in Figure 5.It can be found that the contents of Cu, Pb, Zn, and As increase as the distance from the waterlogged areas decreases, which was in line with theoretical expectations.However, there was no clear pattern observed in the variation in the Cr, Hg, and Cd content in relation to the distance in the soil.cantly negatively correlated with distance (p < 0.05).In Study Area C, Cu and Cr we significantly negatively correlated with distance (p < 0.05), and Cd was significantly ne atively correlated with distance (p < 0.1).In Study Area D, Cu, Zn, and Hg were signi cantly negatively correlated with distance (p < 0.05).In Study Area E, Cu, Pb, Zn, Cr, C and Hg were significantly negatively correlated with distance (p < 0.05).
As shown in Figure 6b, in Study Area A, Pb, As, Zn, and Hg were significantly po tively correlated with SOC (p < 0.05).In Study Area B, As and Cd were significantly po tively correlated with SOC at the 0.05 significance level.In Study Area C, Cu and Zn we significantly positively correlated with SOC (p < 0.05).In Study Area D, As and Zn we significantly positively correlated with SOC at the 0.05 significance level.In Study Ar E, Pb and Hg were significantly positively correlated with SOC (p < 0.05).As shown Figure 6c, in Study Area C, As was significantly positively correlated with pH (p < 0.0 In Study Area E, Pb and As were significantly positively correlated with pH (p < 0.05), an Cd was significantly positively correlated with pH at the 0.1 significance level.To further analyze the relationship between PTEs and the distance, soil SOC, and p in the surrounding soil of water areas with coal mining subsidence, we conducted a line correlation analysis of PTEs with the distance, SOC, and pH using significant correlatio tests.The results are shown in Figure 7.In Study Area A, the explanatory power of t distance on the Cu and Zn content was 0.2473 and 0.5934, respectively.The explanato power of the SOC on the Pb and As content was 0.3312 and 0.3807, respectively.In Stud Area B, the explanatory power of the distance on the Cu content was 0.5645, and the e planatory power of the SOC on the Cd content was 0.5162.In Study Area C, the explan tory power of the distance on the Cu content was 0.3854, and the explanatory power the SOC on the Zn content was 0.7405.In Study Area D, the explanatory power of t distance on the Zn content was 0.3548, and the explanatory power of the SOC on the Z content was 0.3969.In Study Area E, the explanatory power of the distance on the Cu, H Pb, and Zn content was 0.3352, 0.7039, 0.6165, and 0.4995, respectively.The explanato power of the SOC on the Pb and Zn content was 0.5101 for both.Overall, the distance ha a relatively high explanatory power for Cu and Zn, and a relatively low explanato power for Cr and Hg.To further analyze the relationship between PTEs and the distance, soil SOC, and pH in the surrounding soil of water areas with coal mining subsidence, we conducted a linear correlation analysis of PTEs with the distance, SOC, and pH using significant correlation tests.The results are shown in Figure 7.In Study Area A, the explanatory power of the distance on the Cu and Zn content was 0.2473 and 0.5934, respectively.The explanatory power of the SOC on the Pb and As content was 0.3312 and 0.3807, respectively.
were often planted with aquatic plants such as reeds and cattails, which had the ability capability to remove PTEs [50][51][52].This It can help mitigate the level of PTE pollution in aquatic ecosystems to a certain extent.
The main reason for the reduction in PTE content in sediments was the transfer of PTEs in the aquatic food web in fish farming and fish-photovoltaic complementary wetlands [53,54].It is worth mentioning that aquaculture feeds had relatively high contents of Cu and Zn, which can introduce certain amounts of Cu and Zn into the system.However, compared with other PTEs, Cu and Zn led to lower ecological and health risks.Lastly, photovoltaic wetlands integrated solar photovoltaic power generation with wetland ecosystems, providing the dual benefits of water purification and electricity generation.Through these development approaches, it is possible to comprehensively reduce the PTE pollution level in waterlogged areas with coal mining subsidence and achieve sustainable resource utilization and environmental protection.

Cu Showed Clear Diffusion Patterns
The process of PTE migration was complex and influenced by various factors, including the characteristics of the PTEs themselves and environmental factors [55].During the migration of PTEs from soil to waterlogged areas, physical, chemical, and biological interactions occurred between the PTEs and the soil's active components [8][9][10].For example, PTEs can be adsorbed by soil colloids, undergo ion exchange with soil solutions, or be absorbed by vegetation.Relevant studies have shown that the metal content generally decreased with distance [56,57].The findings of this study were similar to the results obtained by Zhengxi Gao [58], where Cu exhibited clear diffusion patterns, with higher soil PTE content observed closer to waterlogged areas.The low-lying topography of the waterlogged area leads to the entry of PTEs into the soil through factors such as rainwater runoff.The smaller the distance to the waterlogged area, the higher the accumulation of PTEs in the soil.
Additionally, the SOC also played an important role in the adsorption and migration of PTEs in the soil.In this study, there was a strong correlation between Cd and the SOC.The SOC can affect the adsorption of Cd, reducing the activity of Cd ions.It can also form complexes or chelates with Cd, altering the speciation of Cd and influencing its migration and transformation.In future studies, researchers can provide clearer explanations for their migration patterns by analyzing the state of PTE occurrence.

Conclusions
This paper investigated the migration patterns and influencing factors of soil from coal mining subsidence areas to waterlogged areas under different remediation modes for the first time.The PTE content, such as Cu, in the soil increased with proximity to the waterlogged area.The migration patterns were influenced by SOC and pH.Landscape development, aquaculture, development of fish-photovoltaic complementary wetland, and photovoltaic wetland development can all reduce the overall PTE pollution in waterlogged areas with coal mining subsidence.Currently, large-scale restoration is urgently needed in waterlogged areas with coal mining subsidence, and these research findings can provide data support for the selection of restoration modes and effective control of heavy metals in such areas.

Figure 1 .
Figure 1.Location of 15 sediment samples and 60 soil samples in the surrounding subsidence waterbodies.

Figure 2 .
Figure 2. Evaluation index of PTE pollution in sediments.(a) with reference to UCG; (b) with reference to the background values in Jiangsu Province.

Figure 2 .
Figure 2. Evaluation index of PTE pollution in sediments.(a) with reference to UCG; (b) with reference to the background values in Jiangsu Province.

EFigure 3 .
Figure 3. Ecological risk of PTEs in sediment.(a) with reference to UCG; (b) with reference to the background values in Jiangsu Province.

Figure 3 .
Figure 3. Ecological risk of PTEs in sediment.(a) with reference to UCG; (b) with reference to the background values in Jiangsu Province.The E i r and RI of the PTEs in the sediment are shown in Figure 3b, with reference to the soil background values in Jiangsu Province.By comparing the contribution of the E i r of individual PTEs to RI, it can be observed that Cd made the highest contribution.In areas A, B, C, D, and E, the contributions of Cd were 42.09%, 71.25%, 78.25%, 69.13%, and 65.18%, respectively.The RI values of the five waterlogged areas showed the following pattern: E (620.18) > D (209.93) > C (184.52) > B (162.94) > A (113.14).The results indicated that area E was at a very high ecological risk level, while areas B, C, and D were at a moderate ecological risk level, and area A had a weak ecological risk.

Figure 3 .
Figure 3. Ecological risk of PTEs in sediment.(a) with reference to UCG; (b) with reference to the background values in Jiangsu Province.

Figure 4 .
Figure 4. Radar plot for hazard index (HI) of each sampling location in the study area for adults and children.

Figure 5 .
Figure 5. Variation trends of PTE content at each sampling point (distance from the waterlogged areas: 5 m, 100 m, 200 m, 300 m).

3. 5 .
Analysis of the Correlation between PTE Content and Distance, SOC, and pH The enrichment and migration of PTEs are influenced by the distance, SOC, and pH.To analyze their correlations, we used the Pearson correlation coefficient method to analyze the correlations between PTEs and the distance, SOC, and pH.The results of the analysis are shown in Figure 6.From Figure 6a, it can be observed that in Study Area A, Cu and Zn were significantly negatively correlated with distance (p < 0.1), while Cr and Hg were significantly negatively correlated with distance (p < 0.05).In Study Area Toxics 2023, 11, 888 11 of 17 B, Cu was significantly negatively correlated with distance (p < 0.1), and Cr and Hg were significantly negatively correlated with distance (p < 0.05).In Study Area C, Cu and Cr were significantly negatively correlated with distance (p < 0.05), and Cd was significantly negatively correlated with distance (p < 0.1).In Study Area D, Cu, Zn, and Hg were significantly negatively correlated with distance (p < 0.05).In Study Area E, Cu, Pb, Zn, Cr, Cd, and Hg were significantly negatively correlated with distance (p < 0.05).

Figure 6 .
Figure 6.Pearson correlation analysis of PTE contents and distance, SOC, and pH (** represen significance at the 0.05 level, * represents significance at the 0.1 level).(a) Distance; (b) SOC; (c) p

Figure 6 .
Figure 6.Pearson correlation analysis of PTE contents and distance, SOC, and pH (** represents significance at the 0.05 level, * represents significance at the 0.1 level).(a) Distance; (b) SOC; (c) pH.As shown in Figure 6b, in Study Area A, Pb, As, Zn, and Hg were significantly positively correlated with SOC (p < 0.05).In Study Area B, As and Cd were significantly positively correlated with SOC at the 0.05 significance level.In Study Area C, Cu and Zn were significantly positively correlated with SOC (p < 0.05).In Study Area D, As and Zn were significantly positively correlated with SOC at the 0.05 significance level.In Study Area E, Pb and Hg were significantly positively correlated with SOC (p < 0.05).As shown in Figure 6c, in Study Area C, As was significantly positively correlated with pH (p < 0.05).In Study Area E, Pb and As were significantly positively correlated with pH (p < 0.05), and Cd was significantly positively correlated with pH at the 0.1 significance level.To further analyze the relationship between PTEs and the distance, soil SOC, and pH in the surrounding soil of water areas with coal mining subsidence, we conducted a linear correlation analysis of PTEs with the distance, SOC, and pH using significant correlation tests.The results are shown in Figure7.In Study Area A, the explanatory power of the distance on the Cu and Zn content was 0.2473 and 0.5934, respectively.The explanatory power of the SOC on the Pb and As content was 0.3312 and 0.3807, respectively.

Table 1 .
Classification of E i r and RI.

Table 2 .
Content of PTEs in sediments.
UCG: upper continental guidelines, TEL: threshold effect level, PEL: probable effect level.