Abstract

To study the remediation effect of hydroxyapatite with different particle sizes, a field in situ experiment was carried out by adding conventional hydroxyapatite (0.25 mm) and microhydroxyapatite (3 μm) and nanohydroxyapatite (40 nm) to the contaminated soil and planting Elsholtzia splendens. The distribution and migration of copper (Cu) and cadmium (Cd) in soil were investigated after 4 years. The results show that the application of three different particle sizes of hydroxyapatite significantly raise the soil pH, total phosphorus, and soil organic carbon. Moreover, the addition of hydroxyapatite can reduce the EXC fraction of Cu and Cd by 73.7%–80.1% and 20.8%–35.2%, respectively. In addition, the concentrations of Cu and Cd in >2 mm, 0.25–2 mm, 0.053–0.25 mm, and <0.053 mm aggregate are significantly increased. This improvement indicates that there are risks which may cause the increasing of total Cu and Cd in the soil where the pollution sources still exist. Furthermore, the content of soil colloid is significantly increased, and the colloidal Cu and Cd distribution percentage have been significantly increased by 49.9%–120% and 30.3%–181%. This result illustrates that the application of hydroxyapatite may greatly increase the possibility of colloid and dust migration of Cu and Cd.

1. Introduction

The contamination of heavy metals in soil has arisen much public concern both domestically and internationally [1, 2]. Previous studies have shown that activities such as mining and smelting can translocate heavy metals into the soil [3]. Meanwhile, excessive accumulation of heavy metals in soils not only reduces soil quality, microbial activity, and crop production but also threatens ecological security and human health [4]. During the past decades, more cost-effective and environmentally friendly techniques have been developed to reduce the harm of heavy metal contamination in soil on human and environment, including bioremediation and integrated remediation. Soil amendments have been widely used in in situ remediation of heavy metal contaminated soils due to their economic and in situ characteristics [5, 6].

A large number of studies have shown that phosphorus-containing materials (e.g., apatite, potassium dihydrogen phosphate, superphosphate, and hydroxyapatite) can effectively reduce the activities of heavy metals, such as Pb, Cd, and Cu in soil and wastewater [7]. Among them, hydroxyapatite has been widely used to remediate heavy metal contaminated soil and sediments owing to its strong adsorption and fixation ability for heavy metals [8, 9]. Although hydroxyapatite cannot reduce the total amount of heavy metal pollutants, it can be combined with heavy metals or transform them from active state to inactive state [10]. It has been demonstrated that the particle size of the material can significantly affect the bioavailability and geochemical stability of heavy metals in soil [11, 12]. However, the results about the influence of the particle size of hydroxyapatite on its passivation effect are controversial. Chen found that compared with phosphate rock powder with 35 μm, the phosphate rock powder with 133–266 μm can reduce the bioavailability of heavy metals in soil more effectively [13]. Dong’s research results suggested that micro-nano hydroxyapatite whose particle size between microhydroxyapatite and nanohydroxyapatite was more helpful to reduce the availability of Cu and Cd in contaminated soil than microhydroxyapatite and nanohydroxyapatite [14]. In addition, researches on the immobilization of heavy metals in soil by hydroxyapatite are mostly completed in the laboratory, but there is a lack of comparative study on field restoration experiments [15]. More importantly, most of the studies about hydroxyapatite only focus on its effect on the heavy metal activity and bioavailability, but there are few researches paying attention to its potential effect on the distribution and migration of heavy metals in soil.

In this study, three different particle sizes of hydroxyapatite (conventional particle size, micro- and nanohydroxyapatite) were selected as test materials to study the existing forms, distribution, and migration of heavy metals in soil after 4 years of remediation in field scale. The objectives of this study were to (1) compare the effects of three particle size hydroxyapatite on the existing forms and chemical fractions of Cu and Cd in soil and (2) investigate the changes of Cu and Cd in the distribution and migration in the smelter-impacted soil. The results would put insight into the stabilization effect of different particle size materials in field scale and remind us to pay attention to the risk management and control of soil after remediation.

2. Materials and Methods

2.1. Study Site

The study site is located in Guixi City, Jiangxi Province, China (116°55′ E, 28°12′ N). The area has a subtropical monsoon climate, with an average annual precipitation of 1808 mm. Farmers had used water containing heavy metals to irrigate for a long time, which led to heavy metal contamination in soil (mainly Cu and Cd) and Cd concentrations in rice exceeding the national food health standards (GB 15201-94). The soil texture is sandy loam, and the pH is 4.63. The contents of soil organic carbon, total nitrogen, total phosphorus, total potassium, available nitrogen, available phosphorus, and available potassium are 16.3 g·kg−1, 1.33 g·kg−1, 0.261 g·kg−1, 2.38 g·kg−1, 67.1 mg·kg−1, 186 mg·kg−1, and 54.8 mg·kg−1, respectively. The concentrations of soil total copper and total cadmium are 517 mg·kg−1 and 0.410 mg·kg−1.

2.2. Reagent and Plants

Hydroxyapatite (pH = 8.40, particle size 0.25 mm) was purchased from Nanzhang Lihua mineral powder factory, Hubei, China; microhydroxyapatite and nanohydroxyapatite (pH = 7.68, particle size = 3 μm; pH = 7. 72, particle size = 40 nm) were purchased from Emperor Nano Material Co. Ltd. (Nanjing, China). Cd and Cu concentration in hydroxyapatite, microhydroxyapatite, and nanohydroxyapatite were 1.18 mg·kg−1and 9.54 mg·kg−1, 3.83 × 10−2 mg·kg−1and 5.85 mg·kg−1, and 3.71 × 10−2 mg·kg−1 and 4.4 mg·kg−1, respectively. In our previous research, the transmission electron microscopy was used to analyze the morphology of these three different particle sizes of hydroxyapatite, and the transmission electron microscopy imaging of the different grain sizes of hydroxyapatite indicated that the shapes of HA, MHA, and NHA were bulk, spherical, and acicular, respectively (Figure 1) [16]. Elsholtzia splendens, which was a Cu-tolerant plant, was selected as the phytoextractor in this study.

2.3. Plot Design

The field experiments were conducted in triplicate and were designed with the land split into randomized plots. Each plot was 6 m2 (3 m × 2 m) and the plots were separated by plastic plates. The treatments applied were (1) 1% hydroxyapatite (HA) (w/w according to the mass of surface 17 cm soil, 22.3 t·ha−1, the same below), (2) 1% microhydroxyapatite (MHA), (3) 1% nanohydroxyapatite (NHA), and (4) the control (CK). The dosage chosen in the experiments was based on our previous experimental results [17]. The materials were applied in every plot on 26 December 2012. Then, they were mixed thoroughly with the soil by rotary tillage and then irrigated with tap water after the application of the materials (500 t·ha−1). Elsholtzia splendens were planted with a spacing of 30 cm × 30 cm (70 plants per plot) in every plot on April 26 in 2013, 2014, 2015, and 2016 with a compound fertilizer (N: P2O5: K2O = 15: 15: 15) before the plants were planted at a spreading of 0.833 t·ha−1. These hydroxyapatites were only applied at the start in 2012 during the study period 2012–2016.

2.4. Sample Collection

The test crop Elsholtzia splendens was harvested at the beginning of December every year. After the plant samples were harvested in 2016, five soil samples were collected from each plot in the 0–17 cm depth and fully mixed to form a composite soil sample (about 5 kilogram), the soil samples were divided into two portions: one portion was aired and ground in experiment for physicochemical and soil colloid analysis, and the other portion was aired for the aggregate analysis.

2.5. Sample Analysis

The soil pH was measured with a glass electrode at a water soil ratio of 2.5 : 1 (PHS-2CW-CN, Bante, Shanghai, China). Walkley–Black procedure was used to measure the soil organic carbon (SOC), total nitrogen (TN), and total phosphate (TP) [18]. Page’s method was used to measure the soil available phosphate (HA) and nitrogen (AN) contents [19]. Soil total potassium (TK) and available potassium (AK) were measured according to Olsen [20]. The soil cation exchange capacity (CEC) was measured according to the ammonium acetate method [21].

The soil aggregates were divided into dry aggregate and wet aggregate and measured by the methods of Panettieri and Elliott [22, 23]. The dry aggregates were divided into aggregates sized >5 mm, 2–5 mm, 1-2 mm, 0.5–1 mm, 0.25–0.5 mm, and <0.25 mm, and the wet aggregates were divided into aggregates sized >2 mm, 0.25–2 mm, 0.053–0.25 mm, and <0.053 mm [24].

The Cu and Cd in whole soil, aggregates, and colloid were measured by atomic absorption spectrophotometry (SpectrAA-220) after the samples were digested [25]. The fractions of Cu and Cd were determined by the sequential extraction procedure [26]. In order to ensure the reliability of the experimental data, a standard soil sample (GBW07405, National Research Center for Certified Reference Materials, China) was used during the experiment. The percentage recoveries of heavy metals ranged from 96% to 107%.

The soil colloid was extracted by siphon method, and the settling time of soil colloid (<2 μm) was calculated according to Stokes law [27]. After the soil colloid suspension was obtained, the content of soil colloid and heavy metal was determined by drying method. The settling time of soil colloid (<2 μm) was calculated based on Stokes law. The contents of soil colloid and heavy metal in soil colloid were determined by drying method of Hu [28]. Briefly, 30 g of air-dried soil, which had passed the 2 mm pore size sieve, was added into 250 mL pure water, dispersing it with ultrasonic wave for 30 min. The dispersed suspension was washed into a 1000 mL high type beaker with 750 mL pure water. The soil colloids of less than 2 μm were extracted by siphon method after stirring with glass rod for 1-2 min and settling for 6 h, 53 min, and 24 s. 20 mL soil colloid suspension was added into Teflon crucible by siphoning and dried at 50°C; then the content of soil colloid was determined by weighing method.

2.6. Statistical Analysis

SPSS20.0 (IBM SPSS, Somers, NY, USA) was used for one-way ANOVA and correlation analysis, and all the graphics were plotted by SigmaPlot 12.5.

3. Results and Discussion

3.1. Surface Soil Properties

The properties of surface soil in plots of the four treatments are summarized in Table 1, and there were no significant changes between different treatments in terms of the total N and K; however, the concentration of total P had been significantly increased by the application of HA, MHA, and NHA. This might be due to that as phosphorus-containing remediation materials, the application of three materials increased the content of total phosphorus in soil [25]. Meanwhile, the contents of available N, P, and K were significantly increased after the application of these three amendments. At the same time, the soil pH was greatly increased by 1.13–1.18 units after the application of HA, MHA, and NHA. Owing to hydroxyapatite being a kind of alkaline material which can consume a lot of H+ when dissolved [formula (1)], the presence of hydroxyapatite can decrease soil exchangeable acid and thus increase the soil pH [16]. The result was similar to a previous study, in which the addition of nanohydroxyapatite significantly increased soil pH by 1.8 units [29]. In addition, although the pH of HA is the highest, the results showed that the MHA had the best improvement on soil pH after three years remediation, followed by NHA and HA. This may be due to the fact that HA contains not only hydroxyapatite but also a certain amount of CaCO3, which can dissolve rapidly under acidic conditions, resulting in its less effective maintenance of soil pH than MHA and NHA [16]. The concentration of soil total nitrogen did not change significantly after remediation; however, the SOC contents in the 0–17 cm layer were significantly increased after the addition of HA, MHA, and NHA, and the highest content was found in the NHA (19.2 g·kg−1). This finding might be attributable to three reasons: (1) the vegetation restoration in this area increases vegetation coverage and reduces soil erosion and nutrient loss after the application of the three materials (Wang et al., 2018), (2) the presence of plant residues, roots, and root exudates also increased the input of organic matter to the soil [30], (3) the change in the microbial types and quantities is caused after cultivation that participated in soil organic matter degradation [31]. Thus, these changes could indirectly affect the content of the soil organic matter:

The total Cu content of soil did not significantly decrease after 3 years of remediation, while the concentration of total Cd in soil increased slightly. These results seemed to deviate from the goal of soil remediation, which was to reduce the concentration of heavy metals in soil, and similar results have been reported in our previous study in this region [25]. This might be due to the heavy metal pollution caused by dry and wet deposition in the air, although the area was no longer polluted by heavy metals brought by irrigation water. The research of Tao found that the contents of Cu and Cd in deposition in this area reached 1973 mg·m−2 and 15.2 mg·m−2, respectively [32]. In addition, by comparing heavy metal input (input of passivation material and atmospheric dry and wet deposition) and export (plant absorption, surface runoff and downward leaching), we found that the addition of passivation materials (conventional hydroxyapatite, lime, and charcoal) improved the adsorption and fixation capacity of Cu and Cd in soil, but the heavy metal output through surface runoff and leaching was reduced. As a result, the total Cd in the remediation surface soil was higher than that in the CK [33, 34]. Therefore, it is very important to completely cut off the input of foreign pollutants in the process of soil passivation remediation.

3.2. Chemical Fractions of Cd and Cu

Heavy metals in soil can be divided into solid and solution phases according to their existing forms, and they mainly exist in the solid phase [35]. Moreover, heavy metals in solid phase can be divided into different fractions according to their solubility, mobility, bioavailability, and potential environmental toxicity; therefore, a single extraction step cannot fully evaluate the toxicity and migration of heavy metals in soil [36]. In order to evaluate the migration and bioavailability of heavy metals in soil, a sequential extraction procedure was used in the present study. The different fractions of Cu and Cd in the soil are listed in Tables 2 and 3 , respectively. The distribution of five fractions of Cu was changed from RES > EXC > OM > Fe-Mn > CA to RES > Fe-Mn > OM > CA > EXC after the materials were applied. The application of the HA, MHA, and NHA could significantly decrease the EXC fraction of Cu (66.5%–80.1%) and increase the CA fraction and Fe-Mn fraction of Cu (94.1%–114% and 52.3%–71.5%). The application of the three materials had no significant effect on the OM fraction and RES fraction of Cu. This may be related to the mechanism of fixing heavy metals in soil by hydroxyapatite, which mainly reduces the activity of heavy metals through ion exchange, surface complexation, and coprecipitation [37, 38]. Similar to the chemical fraction of Cu, the application of materials significantly reduced the EXC fraction of Cd (20.8%–35.2%) and increased the Fe-Mn fraction of Cd (78.1%–127%). However, different from CA fraction of Cu, only MHA and NHA could significantly increase the CA fraction of Cd in soil (47.0% and 132%), and HA had no significant effect on the CA fraction of Cd.

The EXC fraction of Cu and Cd could be significantly reduced by MHA and NHA; the main reason might be that the three materials increased soil pH, reduced exchangeable acid, and exchangeable aluminum content and increased the fixation capacity of soil to heavy metals [16]. The results showed that there was a significant negative correlation between soil pH and the EXC fraction of Cu and Cd, indicating that the soil pH played an important role in affecting the activity of heavy metals in soil. In addition, hydroxyapatite could reduce the activity of heavy metals mainly through ion exchange, surface complexation, and coprecipitation [39]. More importantly, as a phosphorus-containing material, hydroxyapatite can reduce the availability of Cu and Cd in soil by forming phosphate precipitation [40]. Meanwhile, there are differences in the mechanism of adsorption and fixation of Cu and Cd by hydroxyapatite with different particle sizes; the surface and structural characteristics of MHA and NHA make it have more adsorption sites and improve the adsorption capacity of Cu and Cd [14]. Generally, the smaller the particle size, the larger the specific surface area and surface energy and the stronger adsorption and fixation ability to heavy metals [13]. However, there was no significant difference of activity reduce for Cu and Cd between MHA and NHA in our study.

The exchangeable fraction (EXC) is a suitable parameter to evaluate the bioavailability and environmental toxicity of heavy metals [41]. In this study, the exchangeable fraction (EXC) of Cu and Cd was significantly decreased, while the Fe-Mn oxides-bound and CA fraction of Cu and Cd were significantly increased after the application of these three amendments. The results indicated that all these three different particle sizes of hydroxyapatite could significantly decrease the bioavailability of Cu and Cd in this contaminated soil.

3.3. Distribution of Cu and Cd in Aggregate Fractions

Soil aggregate is made up of individual soil particles and organic matter whose size and stability can exert direct effects on soil’s physical, chemical, and biological characteristics and plant growth [42]. The composition of soil aggregates is listed in Table 4. For soil mechanical-stable aggregates, the application of the three kinds of amendments mainly increased the 1-2 mm aggregate (28.9%–57.2%) and reduced the <0.25 mm aggregate (2.62%–15.1%). For soil water-stable aggregates, the combined remediation mainly increased the >2 mm aggregate (31.5%–47.0%) and reduced the 0.053–0.25 mm aggregate (9.48%–27.2%). This study showed that the content of aggregates sized >2 mm was greatly increased after the remediation, and the increase of larger aggregate content was beneficial to improving the mechanical stability and water stability of soil. The stability of the soil aggregates was closely related to the soil organic carbon, which could promote the cementation of soil particles [43]. In this study, the application of these three amendments promoted the vegetation restoration which increased the amount of litter and root exudates, thus benefiting the accumulation of organic carbon (Table 1), and consequently improved the content of stability aggregates [44].

This study was mainly focused on examining the heavy metal in the wet soil aggregate fractions, which have been shown reliable to study changes in the distribution of heavy metals in long-term soil use [45]. Both maximum Cu and Cd concentration occurred in <0.053 mm aggregate, followed by the 0.25–2 mm, >2 mm, and 0.053–0.25 mm aggregate (Table 5). The results showed that the application of these three hydroxyapatites could increase the concentration of Cu and Cd in all particle size aggregates, and this change reached significant level in 0.25–2 mm (37.2%–63.5%) aggregate for Cu and in 0.25–2 mm (27.5%–104%) and <0.053 mm (26.2%–87.6%) aggregate for Cd.

The highest concentrations of heavy metals were found in the smallest aggregates (<0.053 mm), probably due to the larger surface area of the small aggregates [46]. The silt and clay were the main components in the <0.053 mm sized fraction, where metals might act as the binding agents for the clay-polyvalent metal-organic matter complexes [47]. Thus, heavy metals could easily accumulate on their large surfaces by adsorption, forming chelating complexes with the organic mineral colloidal particles in the finest fractions [48]. Furthermore, the Cu and Cd concentrations in all the wet soil aggregate fractions of the treatments were slightly increased, which might be attributed to the higher content of total Cu and Cd in the treated whole soil than CK (Table 1).

In addition, although the concentration of Cu and Cd in the <0.053 mm sized aggregates was the fundamental part, the mass loading levels were not predominant (Figure 2), and the 0.25–2 mm, 0.053–0.25 mm, and >2 mm aggregates played important roles as Cu and Cd reservoirs in aggregate. This was due partly to the higher aggregate content of the 0.25–2 mm, 0.053–0.25 mm, and >2 mm aggregates. Meanwhile, the content of heavy metals in the atmospheric deposition was high in this area, with some of the new input of heavy metals being first adsorbed on the surface of the microaggregates and then the fine soil particles, leading to a high concentration of heavy metals and the formation of macroaggregates under the action of SOC [49]. In our remediation, the content of SOC was increased because of the growth of plants, which subsequently promoted the transition of soil from the microaggregates to macroaggregates, eventually leading to the transfer of heavy metals from microaggregates into macroaggregates. This might be the possible reason for the increased mass loading levels of Cu and Cd in the >0.053 mm sized aggregates following the chemical-plant treatments. This process made the load capacity of heavy metals in larger aggregates increase gradually, while the mobility of heavy metals in larger aggregates was relatively weak [50], which reduced the migration ability of heavy metals in the environment. As a result, the accumulation of Cu and Cd in aggregates after application of amendments increased the total amount of heavy metals in soil due to the input of heavy metals in this area. Therefore, although the passivation remediation could reduce the mobility and biological activity of heavy metals in areas where pollution sources still existed, there was a risk of increasing the total amount of heavy metals in the soil.

3.4. Colloids and Cu and Cd in Colloids

The content of soil colloid was significantly increased by all the three particle size hydroxyapatite (Table 6), and the order of increase was MHA (130%) > NHA (124%) > HA (61.8%). The results illustrated that the addition of hydroxyapatite might affect the composition of soil particles and lead to the transformation between soil particles. Although the hydroxyapatite had similar chemical composition, there were big differences in their particle sizes. Due to the small particle size of microhydroxyapatite and nanohydroxyapatite, they might become a part of soil colloid after adding into soil, which might lead to the content of soil colloid in MHA and NHA treatment higher than that in HA treatment. However, because of the small size, large specific surface area, and the lack of adjacent coordination atoms on the surface, nanohydroxyapatite had high activity, which leads to the aggregation of nanohydroxyapatite [51]. Therefore, the colloid content of NHA treatment was slightly lower than that of MHA treatment.

Hydroxyapatite with three different particle sizes could significantly reduce the concentration of Cu in colloid, and the order of reduction is MHA (20.3%) > NHA (15.2%) > HA (14.2%) (Table 6). The concentration of Cd in soil colloid was significantly reduced by the addition of hydroxyapatite (22.5%), but MHA and NHA had no significant effect on the Cd concentration in soil colloid. The addition of amendments could reduce the concentration of Cu and Cd in the colloid, which might be due to the transformation of other particle size soils into smaller soil colloids. While the contents of Cu and Cd in larger particle size soils were lower, this transformation reduced the concentrations of Cu and Cd in soil colloids. More importantly, we found that the concentration of Cu and Cd in colloid was much higher than that in whole soil (Tables 1 and 6), and the concentration of Cu and Cd in colloid was 3.05–3.81 and 3.99–5.68 times higher than that in whole soil. This was similar to Ajmone-Marsan who found that Cr, Cu, Ni, Pb, and Zn in the soils of five European cities were mainly concentrated in soil particles less than 10 μm [52]. This might be due to the larger specific surface area and more active sites in soil colloids, which lead to the enhancement of adsorption capacity of soil colloids for heavy metals [43].

The product of soil colloid content and heavy metal concentration in colloid is regarded as the fixed component of soil colloid heavy metal, and its ratio to total heavy metal in soil is regarded as the colloidal heavy metal distribution percentage. All of three particle size hydroxyapatite could significantly increase the colloidal Cu and Cd distribution percentage, and the increase range was 49.9%–120% and 30.3%–181%, respectively (Figure 3). These results indicated that 3 kinds of hydroxyapatite lead more heavy metals to soil colloid, and the change of soil colloid content caused by hydroxyapatite might be the main reason for the different colloidal Cu and Cd distribution percentage. Particulate water and air transport are important ways of heavy metal migration in soil environment [53], which are closely related to water transport and workers’ physical contact, respectively [54]. The increase of colloidal Cu and Cd distribution percentage could greatly improve the colloid migration and dust migration possibility of heavy metal. Although many studies have shown that different particle sizes of hydroxyapatite have great passivation effect on heavy metals [55, 56], these results indicated that different particle sizes of hydroxyapatite can reduce the risk of heavy metals to organisms, but the possibility of heavy metals migration with soil colloids has not been considered. Therefore, it is necessary to pay attention to soil colloid content and colloidal heavy metal distribution percentage during soil remediation.

4. Conclusion

Heavy metal was distributed and migrated as an important component of the contaminated soil ecosystem, whose composition and fraction changed significantly in response to a one-time application of HA, MHA, and NHA combined with the planting of Elsholtzia splendens in the smelter-impacted soil. The application of three different particle sizes of hydroxyapatite significantly raised the soil pH, total phosphorus, and SOC. The addition of HA, MHA, and NHA could reduce the EXC fraction of Cu and Cd in soil, thus reducing the bioavailability of heavy metals in soil. In addition, the concentration of Cu and Cd in aggregates was significantly increased, which might cause the increasing of the total Cu and Cd in the soil where the pollution sources still existed. What is more, the content of soil colloid and the colloidal Cu and Cd distribution percentage had been significantly increased after 4 years of one-time addition of three different particle sizes of hydroxyapatite, which might greatly increase the colloid migration and dust migration possibility of Cu and Cd. Therefore, it is not only necessary to focus on the biological activity of heavy metals but also needed to pay attention to monitoring of soil colloid content and colloidal heavy metal distribution percentage, distribution of heavy metals in aggregates, and other indicators related to the distribution and migration of heavy metals, especially in areas where pollution sources still exist.

Data Availability

All the data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

L. X. (Associate Professor) conducted all the experiments and wrote the manuscript. XY X. (Associate Professor), JB P. (Associate Professor), D. L. (Associate Professor), and MF J. (Associate Professor) revised the manuscript. ZQ Z. (Engineer) and HB C. (Associate Professor) conducted part of the experiments. J. Z. (Professor) completed the conception of the article and the design of the experiment.

Acknowledgments

The research was funded by the PhD Special Project of Nanyang Normal University (2018ZX018), the Scientific and Technological Research Projects in Henan Province (212102310844), the National Natural Science Foundation of China (41701361) and the Key Scientific Research Project of Henan Institutions of Higher Learning (grant number 19B420002).