Abstract
Nitrogen and phosphorus are commonly recognized as causing eutrophication in aquatic systems, and their transport in subsurface environments has also aroused great public attention. This research presented four natural clay minerals (NCMs) evaluated for their effectiveness of NH4 + and PO4 3- adsorption from wastewater. All the NCMs were fully characterized by X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), scanning electron microscopy (SEM), BET analysis, and adsorption kinetics and isotherms to better understand the adsorption mechanism-property relationship. The results show that the adsorption efficiency of the four NCMs for phosphate was better than that for ammonia nitrogen. The removal rate of phosphate was higher than 65%, generally in the range of 80%-90%, while the removal rate of ammonia nitrogen was less than 50%. The adsorption kinetic behavior followed the pseudo-second-order kinetic model. The ammonia nitrogen adsorption isotherm was in good agreement with the Freundlich isotherm equilibrium model, and the phosphate adsorption isotherm matched the Langmuir model. Among all the NCMs studied, bentonite (7.13 mg/g) and kaolinite (5.37 mg/g) showed higher adsorption capacities for ammonia nitrogen, while zeolite (0.21 mg/g) and attapulgite (0.17 mg/g) showed higher adsorption capacities for phosphate. This study provides crucial baseline knowledge for the adsorption of nitrogen and phosphate by different kinds of NCMs.
1. Introduction
With the development of human society, the impacts on water quality from industrial and human activity continue to increase [1, 2]. Industrial wastewater and domestic sewage containing nitrogen and phosphorus are still often directly discharged into natural water bodies without effective treatment [3, 4]. Natural reservoirs, lakes, and oceans are facing severe water eutrophication problems [5, 6]. In particular, the high concentration of nitrogen and phosphorus in natural water bodies is the main cause of eutrophication of water bodies [7–10]. How to deal with nitrogen- and phosphorus-polluted water has become a matter of widespread concern.
For the removal of excess nitrogen and phosphorus in eutrophic water bodies, many methods have been developed over several decades of investigations, including physical treatment and chemical treatment, biodegradation, and selective ion exchange [11–14]. For practical applications, adsorption is the most advantageous remediation method due to its low cost, small energy consumption, and relative ease in implementation [15, 16]. Among various adsorbents, clay minerals are highly favorable due to their advantages of easy availability, low cost, and lack of secondary contamination to the ecological environment, and they are key components of soil [6, 17–19]. Four representative minerals including kaolinite, zeolite, attapulgite, and bentonite could be commonly found in soil and sediments. These minerals have good ion exchange capacity and exhibit a certain adsorption effect on pollutants in water or soil, especially on nutrients such as nitrate, ammonia, and phosphate [20, 21]. Clay minerals are the most important adsorbents which have made a significant contribution to the natural self-purification capacity of water and soil environments.
Many scholars focus on the clay mineral’s adsorption property for treating nitrogen and phosphorus compounds, such as nitrate in groundwater and ammonia in surface water, wastewater, or soil [22–26]. For example, natural zeolite has been used to remediate the sediments in eutrophic lakes because it has a high ion exchange capacity and high selectivity for NH4 + [27]. It is proved to be feasible to use bentonite, illite, and zeolite as capping materials to stabilize nutrients and interrupt their release from contaminated lake sediments. Bentonite was especially used to capture phosphate released from sediments [28].
Several factors such as concentration, pH, and temperature could affect the adsorption behavior, especially the concentration of raw water. We could divide the nitrogen or phosphorus concentration of raw water into three types: high level (hundreds of mg/L), medium level (dozens of mg/L), and low level (a fraction of mg/L). For instance, the PO4 3- adsorption capacity of TiO2/Yemeni natural zeolite may reach 37.6 mg/g from the raw water of 80 mg/L phosphate, and the Freundlich model provided the best fit for the equilibrium data [29]. For eutrophic water with 0.1 mg/L phosphate and 2 mg/L ammonia nitrogen, one novel bentonite-humic acid composite Bephos™ could show phosphate adsorption capacities of 26.5 mg/g and ammonia nitrogen adsorption capacities of 202.1 mg/g, respectively. The Langmuir and Freundlich models both fitted well with the experimental results, and the pseudo-second-order kinetic model correlated well with the experimental results [30]. The ammonia nitrogen adsorption of halloysite was saturated with 1.66 mg/g at 303 K, pH of 5.6 in 600 mg/L solution ammonia nitrogen, and the ammonia nitrogen isothermal adsorption of halloysite matched the Langmuir and Freundlich isotherms, and the adsorption process conformed closely to the pseudo-second-order kinetic equation [31].
A previous study mainly adopted modified clay minerals to maximize their adsorption capacity for nutrients. Although a relatively high treatment effect could be achieved, modified clay minerals or materials were commonly applied in the lab or small-scale application, but they were expensive and infeasible to mass production and large-scale use to treat nutrients in natural eutrophic lakes. By comparison, the adsorption behavior of natural clay minerals would have more practical applications, especially for eutrophic water bodies with slightly higher nitrogen of 0.2–0.3 mg/L and phosphorus concentrations of 0.01–0.02 mg/L. Moreover, most existing research studies were independent [32]. There was a lack of systematic and comprehensive comparative studies under the same conditions. Therefore, we focused on the natural clay minerals and chose eutrophic waters with slightly higher nitrogen and phosphorus concentrations, which helped us to understand the structure and properties of NCMs and their adsorption effects on NH4 + and PO4 3-.
This study evaluated four common natural clay minerals, namely, kaolinite, zeolite, attapulgite, and bentonite, as adsorbents under the same experimental conditions, and studied the adsorption behavior of ammonium ions and phosphate on the surface of NCMs. We made a comprehensive comparison of adsorption capacities, adsorption kinetics, adsorption isotherms, and different influencing factors such as pH, adsorbent dose, and temperature. The NCM samples before and after adsorption were characterized by XRD, FTIR, SEM, and BET analysis to determine the mineralogical properties. This research could provide meaningful references for the exploration of the self-purification function of NCMs.
2. Materials and Methods
2.1. Characterization of Natural Clay Minerals
Kaolinite, zeolite, attapulgite, and bentonite were the tested adsorbents in this paper. They were purchased from Gongyi City YuanXiang Water Purification Materials Factory (China). All natural clay minerals are in powder form and in grey. They are aluminosilicate minerals, with cavities and pores with a certain pore size, and have good adsorption and ion exchange properties.
The surface elemental composition of NCMs was determined by XRD (SmartLab SE, Science, Japan). The SEM (FlexSEM 1000, Hitachi, Japan) was used to examine the morphology of NCMs. The functional groups of NCMs were analyzed by FTIR (Nicolet iS50, Thermo Fisher Scientific, USA). The specific surface area was measured by BET analysis (ASAP 2460 3.01, Micromeritics, USA).
The XRD crystallinity of four NCMs was calculated based on the percentage of the contribution of the crystal region in the total area of the XRD intensity peak of the clay mineral. The four NCMs in the samples were identified using the MDI Jade software package (version 6.5).
2.2. Batch Adsorption Experiments
To simulate a water body with a high degree of eutrophication, the concentrations of ammonia nitrogen and phosphate in the experiment were selected to be 20 mg/L and 10 mg/L, respectively.
The adsorption isotherm was determined by using 25 milliliters of phosphate aqueous solutions with different initial concentrations (from 1 to 50 mg/L) and 0.3 g of the sample. The experiment was performed at room temperature (
The adsorption isotherm was determined by batching 25 milliliters of ammonium chloride solutions with different initial concentrations (from 5 to 100 mg/L) and 0.3 g of the sample. The experiment was carried out at room temperature (
The adsorption kinetic data of phosphate and ammonia nitrogen on natural clay minerals at contact times of 5 to 240 min were studied. The experimental conditions were as follows: pH 7, phosphate concentration 10 mg/L, ammonia nitrogen concentration 20 mg/L, adsorbent dosage 0.3 g, and temperature
To investigate the effect of pH on adsorption, adsorption experiments were carried out at pH 2, 4, 6, 8, and 10. The other conditions were as follows: phosphate concentration 10 mg/L, ammonia nitrogen concentration 20 mg/L, adsorbent dosage 0.3 g, and temperature
Experiments were also set up with distinct adsorptive doses, weighing 0.1 g, 0.2 g, 0.3 g, 0.4 g, and 0.5 g of natural clay minerals for ammonia nitrogen and phosphate adsorption experiments, with the following experimental conditions: pH 7, phosphate concentration 10 mg/L, ammonia nitrogen concentration 20 mg/L, and temperature
To evaluate the effects of different temperatures, batch experiments were performed at 25, 35, 45, and 55°C. The experiment was performed in a temperature-controlled water bath shaker to determine the absorption of phosphate and ammonia nitrogen at different temperatures at
2.3. Data Analyses
The removal efficiency (%) and adsorption capacity (
To evaluate the adsorption mechanism of ammonium and phosphate, four linear models of the pseudo-first-order kinetic model, the pseudo-second-order kinetic model, the Elovich model, and the intraparticle diffusion model were adopted [35–38].
To further analyze the adsorption mechanism, the adsorption isotherm model was introduced, including Langmuir II, Freundlich, and Temkin isotherm models [39].
3. Results and Discussion
3.1. Characterization
3.1.1. XRD Characterization Analysis
The XRD patterns show that the four NCMs were almost entirely composed of natural minerals and a small number of impurities, such as quartz, dolomite, and kaolinite (Figure 1(a)). Analysis showed that kaolinite contains a high content of quartz and a small amount of kaolinite; zeolite contains certain amounts of plagioclase and albite; attapulgite contains quartz and a small amount of dolomite; and bentonite contains a large amount of quartz. Therefore, the high-purity natural clay mineral sample reaches more than 90–95%, which can make an important contribution to the adsorption of nitrogen and phosphorus [40].

XRD patterns of four NCMs (a) and bentonite after adsorption of ammonia nitrogen and phosphate (b). Q = quartz; K = kaolinite; D = dolomite; Alb = albite; Mic = microcline; M = montmorillonite.
Among the four NCMs, kaolinite, attapulgite, and bentonite were similar in composition. They all contain quartz, and the characteristic peaks of quartz appear at 20.8°, 26.5°, and 68.8°, indicating that the three NCMs have similar properties. However, zeolite had certain differences due to its composition.
Compared with the XRD patterns before and after the adsorption of ammonia nitrogen, the characteristic peaks of the main reflection of natural clay minerals did not change greatly, and their intensity was reduced due to the adsorption of NH4 +. This result shows that the structure of NCMs did not change much in the adsorption experiment, and most of the adsorbent was adsorbed or exchanged on the surface of NCM particles [40]. In particular, the characteristic peak intensity of the main reflection of zeolite and attapulgite becomes weaker, while that of kaolinite is at 26.6°, and the peak intensity of bentonite at 50.2° also increases. These observations indicate that cation exchange was involved in NH4 + adsorption on NCM.
Compared with the XRD patterns before and after the adsorption of phosphate, the most significant change on NCM after adsorption could be observed markedly on bentonite (Figure 1(b)). The positions of main reflections of bentonite vary remarkably, which increase from 50.1 to 50.2. The significant change in interlayer distance of bentonite further confirms that more PO4 3- is adsorbed [26]. The crystallinity calculation results of clay minerals after adsorbing ammonia nitrogen and phosphate changed, indicating that nitrogen and phosphorus entered the clay mineral interlayer [41]. The crystallinity of attapulgite before adsorption was 61%, and the change of crystallinity after adsorption of ammonia nitrogen and phosphate was 58% and 81%.
3.1.2. FTIR Characterization Analysis
The FTIR spectra before and after the adsorption of ammonia nitrogen and phosphate from the four natural clay minerals are shown in Figure 2. In the FTIR spectrum of kaolinite, three strong absorption peaks appeared at 3436, 3621, and 3695 cm-1, which were the telescopic vibrations of hydroxyl (–OH) in kaolinite, as well as a strong absorption band at 1644 cm-1, which was the absorption peak of H–OH vibrations in H–O–H in kaolinite. Peaks from 400–600 cm-1 belong to Si–O curved vibration absorptions, and those from 1000–1200 cm-1 were from Si–O telescopic vibration absorptions. After adsorption of ammonia nitrogen, the kaolinite spectrum showed signs of weakening absorptions, indicating that the kaolinite has lost part of the crystalline water, and after adsorption of phosphate, the kaolinite spectrum presented obvious signs of absorption enhancement, indicating that the kaolinite contained partially crystalline water. However, the other shapes are basically similar, indicating that the adsorption process did not change the internal structure of kaolinite [42].

FTIR spectra of NCMs before and after adsorption of ammonia nitrogen (N) and phosphate (P).
In the zeolite FTIR spectrum, in the region of 400–1000 cm−1, the bands correspond to the Si2O stretching bands. The general structure of zeolite contains exchangeable cations filled between thin layers of silicate, and the composition of its basic structure indicates that zeolite has a large ion exchange performance. The band at 1672.2 cm−1 corresponds to single –OH stretching and bending, and the slight change suggests the involvement of hydroxyl groups in the adsorption process [43].
In the attapulgite FTIR spectrum, the typical high-frequency peaks at 3413 cm-1, 3556 cm-1, and 3612 cm-1 are caused by different types of hydroxyl telescopic vibrations in attapulgite crystals. More characteristic structural peaks of silicate minerals were observed in the low and medium frequencies. After adsorption of ammonia nitrogen and phosphate, there were obvious signs of weakening of the absorption peak of the attapulgite map features, indicating that the attapulgite participated in the adsorption reaction [44, 45].
In the bentonite FTIR spectrum, the characteristic absorption peaks of bentonite showed signs of weakening after the adsorption of ammonia nitrogen, indicating that the bentonite lost part of its crystalline water. In addition, two characteristic absorption peaks at 529 cm-1 and 464 cm-1 disappeared, which may be due to the adsorption of ammonia nitrogen. Moreover, after adsorption of phosphates, the peak intensity at 529 cm-1 and 464 cm-1 decreased, possibly as a result of interactions between phosphate root ions and hydrogen bonds and ligand exchange on bentonite [24, 46].
3.1.3. SEM Characterization Analysis
Figure 3 shows the morphological analysis comparison diagram of kaolinite, zeolite, attapulgite, and bentonite before and after the adsorption of ammonia nitrogen and phosphate. The typical characteristic morphology of each clay mineral can be observed, namely, layered or tubular slices. The stratification and abundant folds of clay minerals are conducive to the penetration of water, enabling facile adsorption of nitrogen and phosphorus. When the aqueous solution diffuses into the clay minerals, nitrogen and phosphorus are readily captured by the active sites of the clay minerals so that the adsorption system reaches equilibrium in a short time. It can be seen from the figure that the structure of natural zeolite is tightly arranged, and the surface of the zeolite after adsorption of ammonia nitrogen is hollow, loose, and porous, with more pores extending into the interior. This result shows that the zeolite has adsorbed ammonia nitrogen. Attapulgite has a large specific surface area and a rich porous structure. Bentonite has a smooth surface and relatively small pores. Zeolite appears loose and porous, its block structure decreases, the surface is uneven and rough, and the number of pores per unit area is large [47].

SEM patterns of four NCMs before and after adsorption of ammonia nitrogen (a) and phosphate (b).
The surface shape of the four clay minerals before adsorption is irregular and rough. After the adsorption of ammonia nitrogen, these materials are closely arranged and have many fine particles. This effect may be caused by the adsorption of ammonia nitrogen. Recent studies have found that the adsorption capacity of clay minerals is affected not only by their chemical location but also by their surface morphology [48]. These studies show that surface morphology plays an important role in the process of ammonia nitrogen adsorption.
After the adsorption of phosphate, due to the agglomeration of phosphate on the surface of the adsorbent, the gap between the particles is reduced, and the surface of the adsorbent becomes smooth and dense. This observation indicates that adsorption of phosphate by clay minerals has occurred.
3.1.4. BET Characterization Analysis
The N2 adsorption/desorption isotherms and the BJH pore size distributions of attapulgite are presented in Figure 4. Its characteristic feature is a hysteresis loop, mostly attributed to mesoporous materials. From these results, the BET specific surface area, total pore volume, and average pore diameter of the attapulgite were calculated to be 92.91 m2/g, 0.13 cm3/g, and 5.71 nm, respectively. In addition, the BET analysis results of the other natural clay minerals are shown in Table 1. The specific surface area of attapulgite is the largest, followed by bentonite and kaolinite, and that of zeolite is the smallest [29].

The adsorption/desorption isotherms on samples of attapulgite.
The BET surface area, total pore volume, and average pore diameter of four NCMs.
3.2. Factors Affecting Adsorption
3.2.1. pH
The adsorption of ammonia nitrogen and phosphate on NCMs at different pH values was studied (Figure 5(a)). The results showed that the removal rate of attapulgite and bentonite for ammonia nitrogen was greater than that of kaolinite and zeolite, and the removal efficiency of NH4
+ increased with pH from 2.0 to 6.0 and then decreased to 10.0. At lower pH (

The removal rate of ammonia nitrogen (a–c) and phosphate (d–f) under different pH values, dosages, and temperatures.
3.2.2. Adsorbent Dose
The adsorption of ammonia nitrogen and phosphate on NCMs under different adsorbent dosages was evaluated. As the adsorbent dosage increases from 0.1 g to 0.5 g, the removal efficiency of the four kinds of NCMs for NH4 + increases (Figure 5(b)). When the amount of adsorbents increases, the abundance of active sites available for adsorption increases accordingly, so the capacity of adsorbed ammonia nitrogen increases, thereby increasing the removal rate of ammonia nitrogen. However, the addition to the adsorption capacity of the four NCMs for ammonia nitrogen is different. As shown in Figure 5(e), the adsorption of phosphate by the other three NCMs except zeolite increases with increasing adsorbent dosage. When the adsorbent dosage increases from 0.1 to 0.3 g, the adsorption is relatively rapid. This effect can be attributed to the increase in the specific surface area and available binding sites of the adsorbent [29]. With the increase in the amount of adsorbents, the available specific surface area increases, and the adsorption interaction between the two ions and the number of sites also increase, thereby increasing the overall removal efficiency, but the static unit volume of adsorbent dosage and adsorption capacity decrease. Generally, for a constant concentration of the target pollutant, increasing the adsorbent dose increases the number of contact points, thereby enhancing the adsorption of clay minerals. However, increasing the adsorbent dose to the peak results in decreased activity of the adsorbent on the specific surface area per unit volume, which in turn decreases the number of binding sites and reduces the removal rate [50].
3.2.3. Temperature
The adsorption of ammonia nitrogen and phosphate by four clay minerals at different temperatures was evaluated. As shown in Figure 5(c), the adsorption capacity of attapulgite and bentonite for ammonia nitrogen is greater than that of kaolinite and zeolite. As the temperature increases, the adsorption of attapulgite for ammonia nitrogen decreases. As shown in Figure 5(f), with increasing temperature, the adsorption of phosphate by the four NCMs does not change significantly, indicating that the adsorption process does not have much influence in this temperature range. The adsorption of salt increases, but its removal capacity is lower than that of the other three NCMs. The removal efficiency of NH4 + decreases with increasing temperature, indicating that the adsorption of NH4 + is not affected by the increase in temperature and is exothermic during the adsorption process. The comparison temperature has almost no effect on the adsorption of PO4 3-. It is possible that the different valence states of NH4 + and PO4 3- lead to different changes in the surface charge and adsorption mechanism [38].
3.3. Adsorption Kinetic Analysis
The figures and correlation coefficients of the kinetic model are shown in Figures 6 and 7. The fitting effect of the pseudo-second-order kinetic model is good, and the correlation coefficient (

Four linear kinetic models for the adsorption of ammonia nitrogen on NCMs. (a) Pseudo-first-order kinetic model. (b) Pseudo-second-order kinetic model. (c) Elovich model. (d) Intraparticle diffusion model.

Four linear kinetic models for the adsorption of phosphate on NCMs. (a) Pseudo-first-order kinetic model. (b) Pseudo-second-order kinetic model. (c) Elovich model. (d) Intraparticle diffusion model.
Linear form kinetic model parameters of ammonia nitrogen adsorption on NCMs.
Linear form kinetic model parameters of phosphate adsorption on NCMs.
The values obtained by the Elovich model indicate that diffusion may be part of the rate-determining step, confirming that chemical adsorption may be part of the rate-limiting step [34].
Neither the pseudo-first-order equation nor the pseudo-second-order equation can explain the behavior of ammonia nitrogen and phosphate adsorption on natural clay minerals. Therefore, the intraparticle diffusion model was used to evaluate the diffusion mechanism and rate control steps that affect the adsorption kinetics of ammonia nitrogen [25]. The drawn intercept line of the intraparticle diffusion model cannot pass through the origin, indicating that although the intraparticle process of ammonia nitrogen and phosphate on natural clay minerals is a control step, it is not the only rate control step [52]. For simulation, the adsorption rate of natural clay minerals includes film diffusion, surface adsorption, and intraparticle diffusion for simulation [53].
3.4. Adsorption Isotherm Analysis
The adsorption isotherm of ammonia nitrogen was in good agreement with the Freundlich isotherm equilibrium model, and the adsorption isotherm of phosphate was in good agreement with the Langmuir II isotherm equilibrium model (Figure 8). Among all the NCMs studied, bentonite (7.13 mg/g) and kaolinite (5.37 mg/g) showed higher adsorption capacities for ammonia nitrogen, while zeolite (0.21 mg/g) and attapulgite (0.17 mg/g) showed higher adsorption capacities for phosphate.

Three models for the adsorption of isotherms of NCMs on ammonia nitrogen ((a) Langmuir II model, (b) Freundlich model, and (c) Temkin model) and phosphate ((d) Langmuir II model, (e) Freundlich model, and (f) Temkin model).
The Freundlich model assumes that the surface of the adsorbent is not uniform and that there may be multiple layers of adsorption. The Freundlich model is more suitable for the adsorption isotherms of kaolinite and bentonite for ammonia nitrogen and fits the adsorption of phosphate by zeolite well (
Fitting results of the ammonia nitrogen adsorption isotherm.
Fitting results of the phosphate adsorption isotherm equation.
Attapulgite is a hydrated aluminum-magnesium silicate mineral with a structure consisting of parallel ribbons of 2 : 1 layers. Natural zeolites are porous crystalline aluminosilicate minerals whose framework structure is negatively charged and can adsorb exchangeable cations such as Na+, K+, Ca2+, and Mg2+ [55]. Some studies have found that Ca2+ in attapulgite can form a precipitate with PO4 3- to achieve the effect of removing phosphate [56]. Kaolinite is a layered silicate clay mineral, and bentonite generally consists mostly of calcium montmorillonite, with permanent negative charges on its interlamellar sites. To balance this negative charge, cations in water such as Ca2+, Na+, and Mg2+ are readily used as balance electrons, and it also gives ion exchangeability to bentonite [57]. The removal of NH4 + is the ion exchange of cations such as K+, Mg2+, and Na+ for NH4 + in wastewaters [58].
4. Conclusion
Four NCM-based NH4 + and PO4 3- adsorbents were studied and compared, and their adsorption behaviors were elucidated. Intermittent adsorption experiments show that the contact time, initial concentration, adsorbent dosage, and solution pH have significant effects on adsorption. Kinetic studies show that the adsorption process follows the pseudo-second-order kinetic model, and the adsorption isotherm is in good agreement with the Langmuir isotherm equilibrium model. In all NCMs studied, the adsorption capacity for ammonia nitrogen is bentonite>kaolinite>zeolite>attapulgite, and the adsorption capacity for phosphate is zeolite>attapulgite>kaolinite>bentonite. In particular, these four natural clay minerals have better adsorption effects on phosphate than on ammonia nitrogen and are more suitable for the treatment of phosphorus-containing wastewater. The use of NCMs reduces the treatment cost of nitrogen- and phosphorus-polluted water, with a certain reference value for NCMs in natural environment utilization. These research results provide a data reference for the prediction and evaluation of soil and lake eutrophication.
Footnotes
Data Availability
All data generated or analyzed during this study are included in this published article.
Additional Points
Highlights. (1) The adsorption of nitrogen and phosphorus by natural clay minerals (NCMs) was studied comprehensively. (2) The four NCMs had better adsorption efficiency for phosphate. (3) NCMs had application potential in removing nitrogen and phosphate in the natural water bodies or wastewater treatment.
Conflicts of Interest
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Acknowledgments
The authors extend their appreciation to the support by the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research (Grant No. IWHR-SKL-201907), China Energy Investment Corporation 2030 Pilot Project (Grant No. GJNY2030XDXM-19-03.2), National Key Research and Development Program of China (2020YFC1908601), National Natural Science Foundation of China (Grant No. 51878004), Anhui Provincial Outstanding Young Backbone Talents Visiting and Studying Program at Home and Abroad (Grant No. gxgwfx2019011), State Key Laboratory of Safety and Health for Metal Mines (Grant No. 2020-JSKSSYS-02), and Support Program for Excellent Talents in Universities of Anhui Province (Grant No. gxyqZD2021129). This study was also supported by the University Synergy Innovation Program of Anhui Province (Grant No. GXXT-2020-075). Thanks are due for the program from the Collaborative Innovation Center of Recovery and Reconstruction of Degraded Ecosystem in Wanjiang Basin Co-founded by Anhui Province and Ministry of Education, Anhui Normal University, and State Key Laboratory of Safety and Health for Metal Mines, Sinosteel Maanshan General Institute of Mining Research Company Limited.
