Vegetation Restoration in Karst Southwest China: Effects of Plant Community Diversity and Soil Physicochemical Properties on Soil Cadmium

TL;DR

Restoring vegetation in rocky desert areas of China helps reduce harmful cadmium in soil, with one plant mix being most effective.

Contribution

This study identifies optimal plant community structures for cadmium stabilization in karst soils through empirical analysis.

Findings

PMC restoration reduced total and available soil cadmium by nearly 50% and 60%, respectively.

Soil nitrogen availability and plant diversity were key drivers of cadmium regulation.

Vegetation coverage and soil properties like compaction influenced cadmium stabilization.

Abstract

In southwest China, vegetation restoration is widely used in karst rocky desertification control projects. However, mechanistic evidence explaining how plant community composition and species diversity regulate cadmium (Cd) bioavailability remains limited. Here, the plant community’s species diversity, soil properties, Cd, and available Cd contents were evaluated. Four plant community types, NR (natural recovery), PMC (Pistacia weinmannifolia + Medicago sativa + Chrysopogon zizanioides), and PME (Pistacia weinmannifolia + Medicago sativa + Eragrostis curvula), were selected as the research objects. The species composition was recorded, and dominant plant species and soil samples were collected to analyze Cd accumulation characteristics. Relative to NR, composite restorations increased plant diversity and soil nutrient availability and reduced soil compaction, with PMC showing the strongest remediation, decreasing total Cd by 49.4% and available Cd by 59.5%. Model-averaged regression and hierarchical partitioning analyses further identified nitrogen availability and community structure as the dominant drivers. Specifically, available nitrogen (AN), vegetation coverage, Margalef species richness (DMG), ammonium nitrogen (NH4+–N), and total N (TN) were the main factors of soil total Cd, and BD, TN, nitrate nitrogen (NO3−–N), mean crown diameter (MCD), and Shannon–Wiener index (H′) were the main factors of soil available Cd. The results indicate that PMC provides a plant community structure configuration decisions of a scalable, site-adaptable strategy for durable Cd stabilization and soil conservation in thin, carbonate-rich karst soils.

Figures

Full text

1. Introduction

Cadmium (Cd) stands out as the most hazardous heavy metal, displaying a non-degradable nature and an extended half-life in soil [1]. Excessive accumulation of heavy metals in soils can deteriorate soil fertility and reduce its productive capacity [2]. It can harm soil microbes and plant growth, and may pose a threat to human health through the food chain [3]. Cd poses long-term risks to ecosystem functioning and human health due to its bioaccumulative nature [4]. These risks can be categorized into two types: anthropogenic contamination and naturally high background regions [5].

Karst regions characterized by high geochemical backgrounds exemplify the latter, where rock types, diagenetic environments, and tectonic activity have resulted in markedly enriched levels of heavy metal content in soils and rocks, particularly in carbonate rock formations [6]. Compared with non-karst areas, karst ecosystems are highly fragile, owing to their weak rock–soil system and coupled surface–subsurface hydrological framework [7]. Recent evidence further indicates that karst systems are exceptionally sensitive to climate extremes, including intensified diurnal temperature ranges and altered wind regimes, which can amplify soil moisture fluctuations and erosion processes [8]. Under such climate stress, vegetation cover becomes a critical regulator of surface–subsurface hydrological connectivity, thereby influencing Cd mobilization and lateral transport in karst soils [9]. Research in high-geochemical-background regions highlights that the high Cd content is mainly derived from the parent rock [10]. Vegetation restoration has been widely applied in southwest China to mitigate rocky desertification. This approach improves soil aeration and porosity, enhances soil nutrient levels, and reduces heavy metal concentrations, thereby promoting the recovery of ecosystem structure and function. [3,11]. Multiple studies have reported that vegetation recovery can substantially reduce both total and extractable Cd in karst soils [3,12,13]. Other studies have reported different degrees of soil heavy metal pollution across different land utilization and vegetation categories within karst areas [14,15]. Vegetation restoration not only conserves soil and water resources and reduces soil erosion, but also improves soil quality via soil–plant interaction, which is an essential measure for rocky desertification management and ecological reconstruction [16]. Plant community diversity can influence soil Cd content and bioavailability through vegetation-driven changes in the cover, biomass, and litter inputs that affect nutrient cycling [17,18]. It can also influence soil bulk density [19], porosity, organic matter, and nitrogen status, thereby shifting Cd between solid-bound and extractable forms [20,21]. A previous study showed that shrub ecosystems in karst areas exhibit a greater capacity for improving soil quality, thereby promoting the restoration of ecosystem structure and function [22]. Current studies on soil heavy metal contamination in karst regions of southwest China mainly focus on the effects of heavy metal pollution on soil properties [23] and the distribution and remediation [24] of heavy metals under different land-use types [25]. However, most research focuses on single-species responses, and mechanistic evidence on how plant community composition and species diversity influence Cd mobility and extractability in karst soils is limited [26]. Therefore, it is essential to investigate how plant communities regulate soil total and bioavailable Cd by modulating physical structure, soil chemical properties, plant uptake, and functional–complementarity pathways in karst areas.

The karst area of southwest China is characterized by a subtropical monsoon climate with ample rainfall and thermal resources, though these are irregularly distributed spatially and temporally [27]. The area is also characterized by high rock exposure rates, thin and discontinuous soil layers, and sticky and heavy soil textures of low fertility [28]. The Yilong Lake Basin, situated in Shiping County, south-central Yunnan (23°36′–23°48′ N, 102°20′–102°41′ E; Figure 1), spans 360.4 km^2^, with an average elevation of 1414 m above sea level [29]. The basin experiences a subtropical monsoon climate, characterized by a distinct dry season (November–June) and wet season (July–October), with a mean annual precipitation of 928.3 mm, evapotranspiration of 1908.6 mm, and an average temperature of 18 °C. The Yilong Lake watershed primarily comprises Yilong, Baoxiu, and Bazin towns of Shiping County, extending in an east–west orientation [30]. The principal soil types are Dystric Cambi soils, Cumulic Anthro soils, Haplic Acri soils, and Humic Acri soils [30]. Accordingly, the northern lakeshore area of Baxin Town, a representative karst region surrounding Yilong Lake, was selected as the study site to (1) analyze the effects of different restoration measures on community characteristics, soil physicochemical properties, and the concentrations of total and available cadmium in soils, and (2) quantify relationships among plant community structures, soil physical properties, soil chemical factors, and Cd bioavailability. This study provides practical recommendations for ecological restoration in karst areas. It offers perspectives to improve restoration strategies and enhance the efficiency of eco-engineering efforts, enhancing ecological engineering efficiency while promoting the control, stabilization, and sustainable use of Cd-contaminated soils in karst regions, with potential applicability to similar environments worldwide.

2. Materials and Methods

2.1. Study Area

The experimental site of this study is located on the northern shore of Yilong Lake, Baxin Town, Dalongjing Mountain, Shiping County, Yunnan Province, China (23°40′12″ N, 102°36′36″ E) (Figure 1). The area represents a karst mountainous landscape characterized by sparse vegetation and severe soil erosion, with an elevation of 1452 m. The region experiences a subtropical monsoon climate, with a mean annual temperature of 18.3 °C, average annual sunshine duration of 2145.2 h, and annual precipitation ranging between 786 mm and 1116 mm. Rainfall varies greatly among years, with abundant sunlight throughout the four seasons, concentrated precipitation in summer and autumn, and relatively little rainfall in winter. The mountain soil is primarily sandy. The vegetation in the mining area is sparse, mainly including shrubs and herbs. The dominant plants include Imperata cylindrica (L.) Raeusch., Miscanthus sinensis Anderss., Bidens pilosa L., Heteropogon contortus (L.) Beauv. ex Roem. & Schult., Pistacia weinmanniifolia J. Poiss., Leucaena leucocephala (Lam.) de Wit, Photinia serratifolia (Desf.) Kalkman, and Salix myrtillacea Andersson.

According to the soil conditions of the northern bank of Lake Yilong and the physiological characteristics of the selected plants, Pistacia weinmannifolia J. Poiss., Medicago sativa cv. WL525HQ, Chrysopogon zizanioides (L.) Roberty, and Eragrostis curvula (Schrad.) Nees were chosen to establish a two-layer “shrub + herb” restoration model. Seeds of Medicago sativa and Eragrostis curvula were purchased from Beijing Zhengdao Seed Co., Ltd. (Beijing, China), while seedlings of Pistacia weinmannifolia and Chrysopogon zizanioides were obtained from the Kunming Zhengxin Seedling Base (Kunming, China). The experimental site consisted of nine plots (each 4 m × 6 m) with slopes of 22–26°. Plots were separated by asbestos boards buried underground to prevent hydrological interference between adjacent plots (Figure 2). To minimize potential edge effects and hydrological interference, asbestos boards were installed between adjacent plots and buried belowground.

Three planting models were established: PMC (Pistacia weinmannifolia + Medicago sativa + Chrysopogon zizanioides), PME (Pistacia weinmannifolia + Medicago sativa + Eragrostis curvula), and a control treatment (NR) representing the original natural shrub–grass community under natural growth without intervention. Pistacia weinmannifolia individuals with similar growth vigor were selected and planted at a spacing of 1.5 × 1 m. Planting and seeding were completed in June 2022. Chrysopogon zizanioides was interplanted between shrub rows in nine planting holes per plot, with five seedlings per hole. Medicago sativa and Eragrostis curvula were sown as a mixture (1:1 ratio) with a seeding rate of 2 kg per 666.7 m^2^ and a row spacing of 30 cm, ensuring good seed–soil contact. Soils were tilled and leveled before planting to ensure uniform conditions. Each model was replicated three times. Plot arrangements are shown in Table 1.

P. weinmanniifolia J. Poiss, a shrub species belonging to the family Anacardiaceae, is native to Yunnan Province, China. It typically grows in limestone mountain forests or shrublands and is a heliophilous species [31]. M. sativa L. is a perennial leguminous herb with strong ecological adaptability and well-developed root system possesses nitrogen-fixing capability, thereby enhancing soil fertility [32,33]. C. zizanioides (L.) Roberty, a gramineous xerophytic herb, is characterized by rapid growth, an extensive root system, and strong slope-stabilization capacity. This species exhibits broad environmental adaptability, high stress tolerance, and the ability to thrive in nutrient-poor soils [34]. E. curvula (Schrad.) Nees is a perennial grass species in the Poaceae family, known for its high drought tolerance and resistance to soil infertility. Owing to its extensive root system and strong regenerative capacity, it commonly colonizes disturbed or degraded sites and plays an important role in vegetation restoration [35].

2.2. Analysis of Community Traits

Vegetation investigations were conducted in October 2024, recording shrubs and herbs in each plot. Three 1 × 1 m quadrats were randomly selected per plot [36]. For shrub and grass layers, the number of individuals, average height, and canopy width were measured in each quadrat. Height and canopy width were measured with tape lines. Vegetation cover is the percentage of the vertical projection of the plant parts; slope was measured with a handheld GPS device placed on the ground. Species identification was based on regional floras and the Flora República Popularis Senica website [37]. Taxonomic identification was conducted using field classification, literature manuals and monographs, and professional expert identification by the Kunming Institute of Botany. Plant coverage of each plot was visually estimated by three investigators. Basic plot information is shown in Table 2.

2.3. Sample Collection and Preparation

Soils at 0–20 cm were collected by the five-point method in each plot and mixed into one composite. In each plot, two dominant species were sampled, with two mature plants per species representing spontaneous growth. All plant and soil samples were sealed in clean PE bags and delivered to the laboratory [38]. After removing litter, a cutting ring 100 cm^3^ collected soils at 0–20 cm horizontally. Soils were naturally dried, ground, and crushed, passed through 20- and 100-mesh sieves, and sealed at 0–4 °C. Plant samples were divided into shoots and roots. The soil attached to the plant samples was rinsed off with tap water, and subsequently the plants were washed with deionized water 2–3 times. Shoots and roots were deactivated at 70 °C for 30 min, dried at 55 °C to constant weight, pulverized, passed through 100-mesh sieves, and stored sealed for subsequent analysis [39].

2.4. Soil and Plant Chemical and Metal Analysis

BD, non-capillary porosity (NCP), capillary porosity (CP), and total porosity (TPO) were determined by the cutting ring method [40]. Soil organic matter (SOM) was determined by the sulfuric acid–potassium dichromate method (NY/T 1121.6-2006). Total nitrogen (TN) was measured by the Kjeldahl method. Total phosphorus (TP) was determined by the molybdenum–antimony colorimetric method. Available nitrogen was determined by the alkaline hydrolysis diffusion method, and available phosphorus was measured by atomic absorption spectrophotometry method (PerkinElmer Inc., Waltham, MA, USA) [40]. Available Cd was measured by CaCl_2_ extraction and atomic absorption spectrometry [41]. Soil samples were digested in an acid mixture of HNO_3_–HClO_4_–HF (5:4:1) with a graphite digestion oven (CEM Corporation, Matthews, NC, USA) [42]. Plant Cd was determined through the HNO_3_–HClO_4_ mixed-acid method. Cd in digests was quantified by graphite furnace atomic absorption; method setup followed HJ 832-2017 and GB 5009.15-2014. The LODs for soil and plant Cd were 0.01 mg·kg^−1^ and 0.02 mg·kg^−1^, respectively. All samples were analyzed in triplicate, and the relative standard deviation (RSD) of replicate measurements was less than 5%. No significant contamination was detected in blank samples. Two calibration standards were tested for every 10 samples, and the calibration curves showed coefficients of determination (R^2^) greater than 0.990.

2.5. Calculation Methods for Species Diversity and Accumulation Characteristics

Species diversity parameters were calculated as follows [43]: Ecological dominance (Do), Margalef species richness (DMG), Shannon–Wiener index (H′), Simpson’s diversity index (Dsi), and Pielou’s evenness index (Epi).

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$ $$[eqn]$$

In these equations, S is the total number of species, n_i_ is the number of individuals of the i-th species, N represents the total number of individuals, and P_i_ = n_i_/N denotes the relative abundance of the i-th species. The natural logarithm (ln) is used in all logarithmic calculations.

The importance value (IV) of each plant species is calculated as follows [44]:

$$[eqn]$$

In this equation, RD is the relative density of species and RC is the relative coverage of the species. Species with the highest IVs within each treatment were identified as the dominant species and selected for subsequent biomass and cadmium accumulation analyses. Detailed IV values of all recorded species under different treatments are provided in Appendix A (Table A1).

The bioaccumulation factor (BCF) is determined by the ratio of heavy metal concentration in the shoot parts to that in the soil [45]. The translocation coefficient (TC) is calculated as the heavy metal concentration in shoots divided the concentration in roots [46].

2.6. Statistics Analysis

Differences in species richness, evenness index, soil physicochemical properties, total and available cadmium contents, and plant cadmium concentrations among different treatments were examined using one-way ANOVA. Pearson correlation coefficients (significance p < 0.05, very significant p < 0.01) were applied to examine the relationships between species diversity and soil properties, verified using the R packages spaa and corrgram (version 4.4.3). Variables that showed significant correlations with total and available cadmium in the Pearson correlation analysis were selected as predictor variables in the models. Multiple regression models were applied to evaluate the effects of species diversity and soil properties on total and available soil cadmium contents, respectively. All predictors and response variables were standardized using Z-scores to interpret parameters on a common scale. Using the “MuMIn” package [47], all possible combinations of the initial predictors were used to generate a set of candidate models. These models were ranked according to the Akaike information criterion (AIC) based on maximum-likelihood estimation in R. Models with ΔAIC < 2 were selected, and model averaging was applied to estimate parameters and associated p-values using the MuMIn function model.avg. The relative contribution of each independent variable was then calculated based on its parameter estimate relative to the total model effect. Given the limited number of replicates (n = 9), the regression and model averaging analyses were applied in an exploratory manner to evaluate the relative importance of candidate predictors, rather than to build predictive models.

3. Results

3.1. Plant Community Characteristics

A total of 17 plant species belonging to 13 families and 15 genera are recorded across the experimental plots (Table 2). The dominant herbaceous species included Imperata cylindrica (L.) Raeusch., Miscanthus sinensis Anderss., Bidens pilosa L., Medicago sativa cv. WL525HQ, Chrysopogon zizanioides (L.) Roberty, and Eragrostis curvula (Schrad.). Plant abundance ranged from 23 to 93 individuals·m^−2^, showing higher values in the PMC and PME plots than in NR. Diversity parameters varied significantly among the three community types (Figure 3). Species abundance, richness (DMG), and Shannon–Wiener index (H′) were the highest in PMC, intermediate in PME, and the lowest in NR (p < 0.05). Conversely, ecological dominance (Do) was the greatest in NR, indicating that the community was mainly composed of a few dominant species. Simpson’s diversity index (Dsi), Pielou evenness index (Epi), coverage, mean height, and mean crown diameter followed the same pattern as richness and diversity, all showing significantly higher values in PMC than in NR (p < 0.05). Overall, the results suggest that vegetation restoration significantly increased plant abundance and diversity while reducing dominance. The community structure in PMC and PME plots shifted toward a more diverse and evenly distributed pattern compared with NR plots.

3.2. Soil Properties and Cd Concentrations

The physicochemical properties of soils in the three treatments (NR, PMC, and PME) are summarized in Table 3. BD differed significantly among treatments (p < 0.05), being the highest in NR (1.64 ± 0.12 g·cm^−3^) and the lowest in PMC (1.27 ± 0.06 g·cm^−3^). Available nitrogen (AN) and nitrate nitrogen (NO_3_^−^–N) followed the same pattern, while available phosphorus (AP) remained relatively stable across plots. Total and available Cd concentrations varied significantly among plots (p < 0.05), both being the lowest in PMC and the highest in NR. Specifically, total Cd decreased with a reduction of 49.4%, while PME showed a 35.3% decrease relative to NR. Similarly, available Cd declined from 0.042 ± 0.007 mg·kg^−1^ in NR to 0.017 ± 0.003 mg·kg^−1^ in PMC, representing a reduction of 59.5%, while a reduction of 35.7% was observed in PME compared with NR. These results indicate that vegetation restoration significantly reduced both total and bioavailable Cd in the soil, suggesting improved soil environmental quality under PMC conditions. Overall, PMC exhibited lower BD and Cd contents but higher nutrient availability, suggesting a potential improvement in soil structure and environmental quality.

3.3. Biomass Cd Content and Accumulation of Dominant Plants

The highest biomass was observed in PMC C. zizanioides for both roots and shoots, whereas B. pilosa of NR and E. curvula of PME exhibited the lowest values. There was no significant difference in root and shoot biomass of M. sativa under PMC and PME treatments. The roots of PMC C. zizanioides contained the highest Cd concentration (0.167 ± 0.02 mg kg^−1^), while NR I. cylindrica had the lowest (0.032 ± 0.002 mg kg^−1^) (Figure 4B). Cd concentrations in roots were higher than those in shoots across all species. Cd accumulation (Figure 4C) exhibited a similar trend to content, with the greatest values of C. zizanioides in PMC of both root and shoot (p < 0.05). Cd accumulation (Figure 4C) exhibited a similar pattern to Cd content, with the greatest values observed in PMC C. zizanioides for both roots and shoots (p < 0.05) (2.603 ± 0.627 mg plot^−1^ and 1.255 ± 0.177 mg plot^−1^, respectively).

The bioconcentration factor (BCF) and translocation factor (TF) values of Cd for dominant species under different treatments are presented in Figure 5. The BCF of Cd ranged from 0.21 to 2.10 among all species, indicating large differences in Cd uptake ability. BCF and TF of Cd showed significant differences among the dominant species (p < 0.05). C. zizanioides in the PMC treatment exhibited the highest BCF (2.10 ± 0.23), suggesting a strong Cd accumulation potential in roots, while the lowest in NR I. cylindrica. There were no significant differences between M. sativa cv. WL525HQ of PMC, E. curvula, and M. sativa cv. WL525HQ of PME. The TF of Cd varied between 0.26 and 0.65 across the tested species. B. pilosa and C. zizanioides had the highest TF values with 0.65 ± 0.04 and 0.62 ± 0.05, respectively, both significantly greater than in others (p < 0.05), indicating efficient translocation of Cd from roots to shoots. In contrast, I. cylindrica (0.33 ± 0.02), E. curvula (0.29 ± 0.03), and M. sativa cv. WL525HQ (0.31 ± 0.04) had lower TF values, suggesting that most Cd was retained in the roots.

3.4. Correlations Between Environmental Variables and Cd Contents

Species richness, DMG, H′, and coverage were positively correlated with AN, NH_4_^+^–N, NO_3_^−^–N (p < 0.05), indicating that diverse plant communities improved soil nutrient status (Figure 6). By contrast, these diversity variables were negatively correlated with BD (p < 0.05), suggesting that greater vegetation diversity enhanced soil porosity. Moreover, richness, DMG, coverage, mean crown diameter, and H′ showed negative correlations with Total Cd and Available Cd (p < 0.05), implying that higher vegetation diversity was associated with lower Cd enrichment and mobility in soils. Except for BD, both total Cd and available Cd were highly negatively correlated with most soil physicochemical properties, including TPO, CP, AN, NH_4_^+^–N, NO_3_^−^–N. BD showed a significant positive correlation with both total Cd and available Cd, indicating that soil compaction enhanced Cd accumulation and increased its bioavailability.

Multiple linear regression (MLR) and hierarchical partitioning were used to explore potential predictors of total and available Cd concentrations in soil (Figure 7). For total Cd, AN, vegetation coverage, DMG, NH_4_^+^–N, and TN repeatedly appeared in the best-supported MLR models (ΔAIC < 2), indicating their relatively higher importance among the considered variables. The regression equation is expressed as follows:

$$[eqn]$$

Among these variables, AN and DMG exhibited the strongest associations, showing opposite directions: AN negatively associated with soil total Cd, while NH_4_^+^–N and TN positively associated with soil total Cd. The hierarchical partitioning model yielded an adjusted R^2^ of 0.836, suggesting that the selected predictors jointly explained 83.6% of the total variance in total Cd (Figure 7A). Soil variables (AN, NH_4_^+^–N, and TN) contributed 49.21% of the explained variance, while plant community properties (DMG and coverage) accounted for 34.35%. These results indicate that soil N fractions explained a relatively large proportion of the variation in total Cd, while plant diversity exerted a moderate regulatory influence.

For available Cd, the best model incorporated BD, total N, nitrate nitrogen (NO_3_^−^–N), mean crown diameter (MCD), and Shannon–Wiener index (H′) (p < 0.01). The regression equation is as follows:

$$[eqn]$$

In this model, BD, H′, and MCD showed positive effects on Cd bioavailability, while NO_3_^−^–N and TN exhibited negative effects (Figure 7D). The hierarchical partitioning analysis explained 89.3% of the total variation in available Cd (Adj. R^2^ = 0.893; Figure 7C). Soil properties contributed approximately 58.82% of the explained variance, suggesting that soil physicochemical factors played a major role in explaining Cd availability, whereas plant structural and diversity attributes contributed 30.44% of the variation. Together, these results suggest that soil N fractions and plant community composition jointly regulate both Cd accumulation and bioavailability in the restoration system.

4. Discussion

4.1. The Transport and Absorption of Cd by the Dominant Species

The dominant species varied among treatments, with B. pilosa and Imperata cylindrica in NR, C. zizanioides and Medicago sativa in PMC, and M. sativa and E. curvula in PME. The dominance of B. pilosa and I. cylindrica in NR suggests that under nutrient-poor conditions without human intervention, the plant community is dominated by drought- and stress-tolerant weeds. By contrast, both PMC and PME introduced a legume–grass mixture that may increase nutrient acquisition and biomass production, thereby potentially strengthening plant-mediated Cd capture and stabilization. As a leguminous plant, Medicago sativa can improve soil nitrogen availability through symbiotic N fixation with rhizobia, converting atmospheric N_2_ into ammonium [48], and potentially facilitating neighboring species by improving soil fertility and rhizosphere activity. Moreover, the rhizodeposition of amino acids and organic acids by M. sativa has been suggested to provide nutrients to neighboring and may contribute to the alleviation of Cd toxicity through root–root interactions, according to previous studies [49]. Medicago sativis [50], a taproot species that interacts, respectively, with the fibrous-rooted C. zizanioides [51] and E. curvula [52] may improve resource utilization and increase plant biomass through niche complementarity, and ultimately improve Cd pollution remediation [53]. This process may enhance its own biomass accumulation and resistance to heavy metal stress through the synthesis of chelating peptides and activation of antioxidant systems, as reported in previous studies [54]. For instance, Xiong et al. [48] found that the intercropping of M. sativa with Ricinus communis significantly improved the latter’s Cd and Zn extraction efficiency due to enhanced soil fertility, consistent with the facilitative role of M. sativa observed here. In this study, Cd concentrations were generally higher in roots than in shoots for all species, suggesting that root retention was a major pathway in the restored karst soils. Among the dominant species, C. zizanioides in PMC exhibited the highest Cd accumulation, with root and shoot Cd concentrations of 0.167 ± 0.020 mg/kg and 0.103 ± 0.011 mg/kg, respectively, and BCF and TF values of 2.10 ± 0.23 and 0.62 ± 0.05, respectively. This result suggests strong Cd enrichment in roots, consistent with Phusantisampan et al. [55], who emphasized that vetiver is an effective phytostabilizer capable of immobilizing Cd in roots. The large and fibrous root system of C. zizanioides provides vast surface area for metal adsorption and sequestration, which is consistent with its recognized phytostabilization potential [56]. In addition to physical adsorption, dense fibrous roots may enhance rhizosphere water retention and dampen short-term hydraulic fluctuations, thereby potentially limiting Cd leaching during rainfall pulses. This process is particularly relevant in karst regions, which are highly sensitive to climate extremes such as intensified diurnal temperature variability and altered wind regimes [8]. Recent evidence further suggests that Cd exposure may reshape the vetiver rhizosphere microbiome and stress-responsive transcription, supporting a rhizosphere-assisted stabilization mechanism rather than simple passive adsorption alone [57]. Therefore, the dense fibrous root system and enhanced rhizosphere recruitment of C. zizanioides facilitated stronger Cd stabilization and a more pronounced reduction in available Cd, which is likely associated with an advantage over PME. The BCF value of C. zizanioides in this study was higher than that reported by Ng and Abas [58], which may be explained by differences in soil pollution levels or co-cultivation with M. sativa, which may enhance Cd uptake via rhizosphere facilitation.

In conclusion, revegetation may promote community complexity and functional root traits that are important for Cd stabilization in the thin carbonate-rich soils typical of karst landscapes, where metal mobility is largely governed by rhizosphere processes [59] and soil structure rather than bulk chemistry alone [60]. The legume M. sativis can supply nutrients and rhizosphere carbon that may stimulate microbial activity and root growth, while C. zizanioides provides a dense root matrix that is favorable for Cd adsorption and physical immobilization.

4.2. Response of Community Functional Diversity to Soil Chemical Properties and Cd Concentration

Vegetation restoration markedly altered the community structure and diversity, which in turn was associated with improvements in soil physicochemical properties and reduced Cd contamination [61,62]. Compared with the NR, the composite restoration treatments (PMC and PME) supported higher plant abundance, richness, and evenness, indicating the establishment of a more complex and stable shrub–grass community structure. In particular, PMC exhibited the greatest species richness and H′, accompanied by the lowest ecological dominance, suggesting strong interspecific complementarity among shrub, leguminous, and grass species.

The results of the Pearson correlation heatmap showed that species richness, H′, and vegetation coverage were positively correlated with SOM, TN, and AP, whereas all diversity parameters were negatively correlated with BD, total cadmium, and available Cd. These relationships suggest that higher diversity may enhances soil nutrient and structure (lower BD), thereby potentially reducing Cd availability—an effect that is especially important in karst ecosystems characterized by thin soils, severe nutrient deficiency, and high vulnerability to degradation [63]. The negative correlations between all diversity parameters and BD (p < 0.05) indicate that the underground intricate root network of shrubs and herbs interpenetrate the soil of vegetation community structure, increasing soil pores and providing preferential flow paths, thereby reducing compaction and potentially improving water content [64]. This mechanism is consistent with the findings of Su et al. [65], who demonstrated that soil water retention characteristics and structural heterogeneity strongly regulate water and solute movement. The roots of shrubs are well developed [66], whereas those of herbaceous plants are dense [67], and both can continuously release various chemical substances (such as sugars and organic acids) throughout the year, significantly influencing soil aggregate structure and microbial activity [42]. Root exudates, litter inputs, and microbial activation are widely reported to promote aggregate formation and organic matter accumulation, which may contribute to reducing BD and improving porosity and aeration [68]. Root exudates can regulate the composition and functional performance of the soil microbial community, promoting the production of specific functional molecules (such as glomalin-related proteins and extracellular polymers), thereby enhancing the binding among soil particles [69]. Improvements in soil physical structure further facilitate nutrient retention and cycling [70]. Moreover, plant diversity increasing plant diversity was associated with enhanced biomass accumulation and distribution, as well as litter and root inputs and decomposition [71], thereby influencing soil nutrient status and fertility [72]. In the PMC and PME treatments, leguminous species enhanced NH_4_^+^–N and NO_3_^−^–N via N fixation, while grasses increased soil carbon inputs and were associated with Cd immobilization, which may have enhanced extracellular enzyme activity and nitrogen mineralization rate, thereby helping to explaining the simultaneous increase in NH_4_^+^–N and NO_3_^−^–N concentrations [73]. According to Wang et al. [74], vegetation restoration improves soil physicochemical properties and biotic functions; their study of spontaneous restoration of Pb–Zn mining wastelands showed that revegetation increased SOM and TN and decreased heavy metals, indicating that plant establishment and microbial interactions are associated with enhanced nutrient cycling and reduced metal mobility. Moreover, both total and available Cd concentrations decreased by nearly half under PMC compared with NR, indicating an apparent improvement in soil environmental quality following vegetation restoration and community diversification. However, this pronounced and rapid decline cannot be attributed solely to plant uptake and rhizosphere-mediated stabilization. In karst landscapes, strong spatial heterogeneity in parent material and soil geochemical background may generate substantial variability in baseline Cd concentrations, which could partially contribute to the observed treatment differences [75]. In addition, enhanced vegetation cover under composite restoration may reduce surface Cd concentrations through erosion control and lateral redistribution, particularly during intense rainfall events characteristic of subtropical karst regions [3]. The accumulation of soil organic matter following revegetation may also induce a dilution effect, whereby increases in soil mass and carbon inputs lower Cd concentrations on a mass basis without proportional changes in total Cd stocks [61]. Taken together, the pronounced decline in Cd under PMC likely reflects the combined effects of biological stabilization, altered soil physical structure, and landscape scale redistribution processes, rather than a single dominant mechanism.

Cd concentrations are closely linked to rhizosphere processes and microbial regulation [61]. Wu et al. [3] reported that vegetation restoration improved soil structure by increasing soil aggregates and was associated with reduced heavy metal concentrations through plant uptake and microbial biosorption. As suggested by previous studies, in karst soils, vegetation restoration may modified rhizosphere chemistry, shifted Cd toward more stable residual and carbonate-bound fractions [76]. Previous studies have shown that diverse plant communities also increased microbial metabolism and enzyme activities (urease, phosphatase, catalase), thereby potentially enhancing Cd adsorption and deposition [74]. Similarly, Wu et al. [3] confirmed that root exudate and aggregate formation were jointly associated with reduced Ni, Pb, and Cd activity and improved soil stability. Total and available Cd decreased by 50% in PMC, highlighting the potential phytostabilization and ecological risk mitigation potential of diverse shrub–grass systems in karst environments. Our findings also align with other karst restoration studies showing that vegetation recovery may reduce soil bulk density and improve soil nutrient status, microbial functioning, and Cd immobilization more effectively than natural succession in rocky desertification landscapes [11,40,77]. However, microbial-mediated mechanisms remain underexplored. Insights from disturbance–recovery studies indicate that microbial communities may maintain functional resilience under coupled climatic and edaphic stress. For instance, Liu et al. [78] demonstrated that ecosystem recovery following environmental perturbation is strongly mediated by microbial functional stability, offering a useful parallel for understanding Cd immobilization in stressed karst ecosystems.

Beyond edaphic constraints, karst ecosystems are climatically vulnerable, and vegetation–soil feedbacks may modulate drought–flood alternation by improving infiltration and water holding capacity, thereby influencing nutrient cycling and, potentially, metal mobility [63]. At broader spatial scales, vegetation restoration may contributes to climate and hydrology feedbacks through altered evapotranspiration and terrestrial water cycling, as demonstrated by global greening studies [9,79], providing an important climatic context for interpreting local Cd stabilization.

Model-averaged regression and hierarchical partitioning identified the dominant drivers, indicating that soil properties were associated with a stronger control of total and available Cd concentrations than plant community attributes. Nutrient availability also appeared to modulated Cd accumulation. MLR identified AN, coverage, DMG, NH_4_^+^–N, and TN as significant predictors for total Cd (p < 0.01). The negative effect of AN on total Cd may be related to active N availability stimulates microbial processing and extracellular enzyme activity that can foster cadmium complexation with organic ligands and mineral surfaces, thereby potentially promoting Cd into less mobile soil pools [61,80]. The positive effect of NH_4_^+^–N on total Cd may be attributed to the nitrification process induced by NH_4_^+^–N, which releases H^+^ and decreases rhizosphere pH, thereby potentially causing Cd to desorb from carbonate or oxide binding sites [81]. For available Cd, the best MLR model included BD, TN, NO_3_^−^–N, MCD, and H′ (p < 0.01). This observation may be explained mechanistically: NO_3_^−^–N uptake by plants is commonly associated with rhizosphere alkalinization (relative to NH_4_^+^–N uptake), and a higher rhizosphere pH can reduce Cd solubility and favors adsorption to mineral and organic binding sites, thereby reducing the content of available Cd [81]. For example, Hu et al. [82] and Bai et al. [83] reported nitrate-facilitated Cd uptake or increased plant Cd accumulation for specific species, which is consistent with nitrate increasing root uptake fluxes, even while reducing extractable soil pools via plant removal. The role of nitrogen forms in Cd mobility has been widely discussed, and comparisons within coupled nutrient–metal systems provide additional mechanistic insight. Wang et al. [84] showed that coordinated regulation of nitrogen transformations can simultaneously suppress arsenic mobilization and N_2_O emissions, offering a valuable conceptual framework for interpreting coupled N–Cd dynamics in restored soils. Higher BD reduces macro porosity and gas diffusion, creating a low-oxygen environment that may promote microbially mediated reduction and dissolution of Fe/Mn oxides, thereby potentially releasing adsorbed Cd into the soil and increasing the content of available Cd [70,85]. Thus, the coupling of plant diversity, soil nutrient dynamics likely jointly governs Cd mobility and stabilization in restored karst soils, reinforcing the ecological benefits of multi-species revegetation.

From an application perspective, this provides clear practical meaning for karst restoration: PMC-type shrub–legume–vetiver assemblies may represent implementable nature-based solutions that can help stabilize fragile soils and mitigate Cd ecological risk by lowering phytoavailable Cd rather than relying solely on harvest-based removal.

5. Conclusions

This study shows that vegetation restoration was associated with higher plant species diversity, vegetation cover, and the productivity of newly established plant communities, as well as higher soil nutrient levels and lower concentrations of total and available Cd. These patterns suggest a close relationship between plant species diversity and ecological restoration in karst regions. Compared with natural recovery, both shrub–grass mixtures were associated with improved soil conditions and reduced Cd, with PMC (Pistacia weinmannifolia + Medicago sativa + Chrysopogon zizanioides) performing the best, reducing total Cd by 49.4% and available Cd by 59.5%. Soil AN, NH_4_^+^–N, TN, vegetation coverage, and Margalef species richness (DMG) were the main factors affecting soil total Cd content. Soil BD, TN, NO_3_^−^–N, MCD, and H′ were the main factors affecting soil available Cd content. Based on these results, PMC may represent a promising plant community configuration for the long-term phytostabilization of Cd-contaminated karst soils. Future studies should improve research quality by integrating advanced analytical approaches. Synchrotron-based spectroscopic techniques (X-ray absorption fine structure), high-throughput sequencing of rhizosphere microbial communities, and isotope tracing (^15^N) can be applied to identify Cd binding forms in soils, clarify microbial contributions to Cd stabilization, and quantify vegetation-mediated nitrogen transformation, respectively. In addition, long-term monitoring and field-scale validation should be incorporated to enhance early-warning capacity and ensure the sustainable application of phytostabilization in karst ecosystems.