Functional Diversity of Plant Communities in Relationship to Leaf and Soil Stoichiometry in Karst Areas of Southwest China

Assessment of relationships between functional diversity and ecological stoichiometry in plant communities can aid in determining the relative variability and ecological complementarity of functional attributes among species, which is a better approach to understanding ecosystem processes and functions than studying species taxonomic diversity. Here, we analyzed the relationships among community weighted means of functional traits, functional diversity, and leaf and soil chemical properties of plant communities during various stages of vegetation restoration in Mao Lan National Karst Forest Nature Reserve, located in humid subtropical Guizhou of China. Our results showed significant changes in four weighted functional traits of plant communities at different restoration stages, namely, plant height, leaf width to leaf length ratio, and leaf area. Additionally, with the progression of the recovery of plant communities, functional richness, functional separation, and quadratic entropy, the coefficient tended to increase. Functional divergence tended to gradually decrease. The association of functional diversity with soil chemical properties was stronger than that with leaf ecological stoichiometry. Regarding leaf and soil chemical properties, soil phosphorus content and leaf C:P were particularly important in influencing functional diversity. Our overall findings indicate that functional traits shift from “acquisitive” to “conservative” as the community is restored. Karst plant communities reduce interspecific resource competition as restoration proceeds, thereby increasing functional overlap effects.


Introduction
The relationship between biodiversity and ecosystem function is currently one of the focal issues in the field of ecology [1]. The ecosystem function needs to consider the different functional properties of each species within the community [2]. Functional diversity of plants represents the range of values and variations in functional traits of organisms within communities and ecosystems [3,4] as well as compositions and variations of morphological, physiological, or other traits that influence plant function [5] and highlights differences in functional traits of species within communities [6][7][8]. Species are likely to exhibit significant differences in functional diversity due to variations in their functional traits [3,9]. Therefore, an increasing number of scholars have advocated the study of the diversity of functional traits instead of the taxonomical diversity of plant species [10][11][12]. Thus, the establishment of the relative variability and ecological complementarity of functional attributes among species that have a more direct role in determining ecosystem function is of considerable research interest [9,[13][14][15]. It was found, for instance, that the relationship between functional diversity and succession was not significant in wetland ecosystems [16]. Other studies found that functional diversity declined as plant species richness and community mean specific leaf area declined over the course of the succession study, and all of these trends  Table 1. Sample site information overview. "-" represents the same as the previous column.

Sample Plot, Sample Square Setting and Plant Community Survey
The degradation of karst forests in the study area is considered to be an important cause of disturbance, which is mainly in the form of fire, raising livestock, grazing, and woodcutting. In 1984, there was the establishment of a nature reserve in the area to stop disturbing the plants in the reserve, thus creating a natural recovery [37]. All of our sample

Sample Plot, Sample Square Setting and Plant Community Survey
The degradation of karst forests in the study area is considered to be an important cause of disturbance, which is mainly in the form of fire, raising livestock, grazing, and woodcutting. In 1984, there was the establishment of a nature reserve in the area to stop disturbing the plants in the reserve, thus creating a natural recovery [37]. All of our sample sites are in the Maolan Karst Nature Reserve. Due to the different degrees of human disturbance, the succession sequence of plant community is not consistent in the restoration process. The interference is larger in the test area, followed by that in the buffer area and least in the core area. Therefore, the early restoration grassland stage and grass irrigation stage were mainly distributed in the experimental area with large disturbance. In the buffer zone with moderate disturbance, shrub stage, tree stage, and tree stage were distributed. The ultimate phase is mainly distributed in the core region with less disturbance. Previous studies have shown that the method of "space instead of time" [38] can be used to study the plant community at each stage, and the vegetation restoration process in this area can be divided into six stages: herbaceous stage (HE), grass and shrub transition stage (HS), shrub stage (SH), tree and shrub transition stage (TS), tree stage (TR) and ultimate stage (CL). In this study, three sample plots were selected as replicates of each stage, leading to a total of 18 standard plots for analysis. To facilitate the plant survey, small sample squares were set up at each site.
The sample squares in TS, TR, and CL stages were set up in the same way; among them, 9 tree layer small sample squares with an area of 10 m × 10 m were set up. One 4 m × 4 m shrub layer small sample was set in each tree small square. One 1 m × 1 m herbaceous layer small sample square was set in each shrub layer small sample square. Therefore, 9 tree layer samples, 9 shrub layer samples, and 9 herb layer samples can be set up in each stage of TS, TR, and CL. Thus, there are 27 small samples in each stage of TS, TR, and CL. Since there are three replicate sample plots, there are 243 small samples in these three stages of TS, TR, and CL.
The sample plots in the HS and SH stages were set up in the same way; among them, 9 shrub layer mini samples with an area of 4 m × 4 m were set up. One 1 m × 1 m herbaceous sample was set in each shrub layer sample. Thus, 9 shrub subsamples and 9 herb subsamples can be set up in each stage of HS and SH. Thus, there are 18 small sample squares in each stage of HS and SH, and since there are 3 replicate sample plots, there are 108 small sample squares in these two stages of HS and SH.
In the HE stage, 10 small sample squares of 1 m × 1 m were set, and since there were 3 replicate sample plots, there were 30 small sample squares in the HE stage.
Therefore, we surveyed a total of 381 small sample squares in 6 restoration stages. Finally, the names, plant height, diameter at breast height, and crown width of trees and shrubs were recorded, along with the names, plant numbers, average height, and cover of herbaceous plants, to facilitate the determination of weighted functional traits of the community (Table 2).

Sample Collection and Determination
The method of measurement sampling was divided into two phases. The first phase involved selecting all the plants in the sample square, cutting off the branches of the sampling plant in the four directions of the south, east north, and west of the canopy with high pruning shears, picking about 20 healthy disease-free leaves on each branch, and mixing them into one sample. The second phase was to select plants with the top two importance values (dominant species) in the sample square as the sampling plants and mix each dominant species from each stage into one sample. The importance value of the tree layer was calculated as: (relative abundance + relative frequency + relative dominance based on diameter at breast height)/3, and shrub and grass layers as (relative abundance + relative frequency + relative cover)/3. The "S" method [39] was used in sample plots of each stage to sample the soil from 0 to 20 cm within a small square (less than 20 cm to the actual depth). Soil was sampled from 0 to 20 cm (if <20 cm, the actual depth was determined) and equal volumes mixed into a single sample. We took 5 soil samples from each sample plot in each stage, with 3 replicate plots per stage, so there were 15 soil samples per stage; for the 6 stages, the total consisted of 90 soil samples. Several leaf samples from the first phase, 36 leaf samples from the second phase, and 18 soil samples were collected in total. Plant samples were desiccated at 105 • C for 2 h, followed by incubation at 75 • C until a consistent weight was obtained, while soil samples were air-dried at room temperature. All samples were finely ground, sieved through 60-mesh screen, and stored in a ventilated area for nutrient analysis.

Functional Trait Selection and Sample Analysis
Six quantitative traits representing functional diversity, including plant height (PLH), leaf thickness (LT), chlorophyll content (CHL), leaf dry matter content (LDMC), leaf area (LA), ratio of leaf length to leaf width (RLW) and specific leaf area (SLA), were selected for analysis. The new global manual for standardized measurement of plant functional traits was utilized [40]. Leaf thickness was measured using electronic vernier calipers (Deli, DL91150, Ningbo Deli Group Limited, Ningbo, China). Leaf length and area were calculated using a scanner in combination with Photoshop software (HP, HPScanJetN92120, Qingdao Lekai Office Culture Supplies Co., Qingdao, China), and chlorophyll content measured using a specific chlorophyll meter (Linde, LD-YD, Shandong Lainde Only Technology Co., Weifang, China). PLH was recorded based on actual measurements in the field and LDMC obtained from the ratio of dried leaf weight to fresh leaf weight. Organic carbon content (SCC) was determined by oxidation with the potassium dichromate external heating method. Plant samples were decocted using the H 2 SO 4 -H 2 O 2 method and total nitrogen (LNC) (NY/T2017-2011) and total phosphorus (LPC) (NY/T2017-2011) contents assessed using indophenol blue colorimetric and molybdenum antimony anti-colorimetric procedures. Soils were analyzed via Kjeldahl nitrogen (SNC) (LY/T1228-2015) and NaOH fusion-molybdenum antimony anti-colorimetric method for the determination of total phosphorus (SPC) (LY/T1232-2015) methods. These measurement standards are derived from the Chinese forestry industry standard, which has been cited by some researchers. (http://www.gov.cn/fuwu/bzxxcx/bzh.htm, accessed on 10 July 2021).

Calculation of Community-Weighted Functional Traits
Karst plant community weighted functional traits (CWM) were obtained by weighting the functional trait values and relative abundance of species [41] to obtain average values for each functional trait at the community level [42]. The parameter was calculated using the formula: where S represents the number of community species, P i the relative richness of species i, and V i the characteristic value of a functional trait of species i.

Functional Diversity Calculation
Five functional diversity indices with three dimensions (functional richness, functional evenness, and functional dispersion) were selected to characterize the functional diversity of communities [43]. Functional richness quantifies the size of community functional niches and reflects the degree of spatial resource utilization by communities, which is characterized by the functional richness index (FRic). Functional evenness indicates the uniformity of spatial distribution of community functional traits, which is expressed by the multidimensional functional evenness index (FEve). Functional dispersion represents the dispersion of community functional trait values and characterizes the degree of interspecific ecological niche complementarity. Functional dispersion consists of functional divergence index (FDiv), functional separation index (FDis) and quadratic entropy coefficient (RaoQ). FDiv functional divergence reflects the overall dispersion of community traits, and the degree of dispersion is inversely related to the strength of interspecific resource competition. FDis functional separation concentrates on the degree of complementarity of plant community niches, and the degree of functional separation in an ecosystem is inversely related to the degree of niches differentiation, and positively related to the effect of overlapping niches and the degree of resource competition [7,[43][44][45].
where SF ic represents the ecological niche space occupied by the community, R c represents the ecological niche space occupied by trait C in the community; S is the species richness, PEW i is the local weighted evenness of species I; a j is the multiplicity of species j, zj is the distance from species j to the weighted center of mass; d ij is the Euclidean distance, 0 ≤ d ij ≤ 1, indicating the dissimilarity of species i and j in a set of trait spaces; C i represents the value of the ith functional trait, A i represents the relative richness of the ith functional trait, and lnx represents the weighted mean of the natural logarithm of the species trait values. Microsoft Excel 2016 was used to conduct preliminary collation of the data. Before analysis, normality and homogeneity of variance of the data were tested. CWM and functional diversity index were calculated by FD software package in R4.1.2, and SPSS25.0 statistical software was used for data analysis. One-way ANOVA and Tukey HSD were used to compare the differences in weighted functional traits, functional diversity index, leaf and soil nutrient contents, and ecological stoichiometry in different restoration stages. Pearson correlation analysis was used to reveal the interaction between these indicators, and the data expression form was mean ± standard value. In order to further verify the changes in functional diversity and leaf and soil stoichiometry, "WGCNA", "Corrplot" and "GGplot2" program packages in R4.1.2 were used for network analysis and correlation analysis of functional diversity and ecological stoichiometry. Finally, we scatterplotted two variables that are correlated and then performed a linear correlation analysis of these two variables using Spearman analysis. The essence of this linear correlation analysis is that when one variable changes, the other variable changes correspondingly in approximately the same way, and the two variables are said to be linearly correlated. Linear correlation analysis was used to further reveal which stoichiometry has a greater effect on functional diversity.

Weighted Functional Traits of Communities at Different Restoration Stages
As shown in Figure 2, the weighted functional traits of the community changed at different stages with the progress of the restoration. PLH showed an increasing trend and reached the maximum at the ultimate stage, and at the late recovery stage (TS, TR, CL) was significantly higher than that at the SH stage, HS stage, and HE stage. LDMC increased first and then decreased, and the LDMC in the SH stage was significantly higher than that in other stages. There was no significant difference in CHL in HE, HS, SH, and CL stages, while CHL in TS and TR was significantly higher than that in other stages. LT showed a decreasing trend, and LT in HE and HS stages was significantly higher than that in other stages. RLW increased first and then decreased, and reached the maximum at the shrub stage, which was significantly different from other stages. LA was the largest in the grassland stage and significantly higher than that in other stages. SLA changes in a fluctuating pattern. HE stage is significantly lower than other stages, and the CL stage is significantly higher than other stages.

Weighted Functional Traits of Communities at Different Restoration Stages
As shown in Figure 2, the weighted functional traits of the community changed at different stages with the progress of the restoration. PLH showed an increasing trend and reached the maximum at the ultimate stage, and at the late recovery stage (TS, TR, CL) was significantly higher than that at the SH stage, HS stage, and HE stage. LDMC increased first and then decreased, and the LDMC in the SH stage was significantly higher than that in other stages. There was no significant difference in CHL in HE, HS, SH, and CL stages, while CHL in TS and TR was significantly higher than that in other stages. LT showed a decreasing trend, and LT in HE and HS stages was significantly higher than that in other stages. RLW increased first and then decreased, and reached the maximum at the shrub stage, which was significantly different from other stages. LA was the largest in the grassland stage and significantly higher than that in other stages. SLA changes in a fluctuating pattern. HE stage is significantly lower than other stages, and the CL stage is significantly higher than other stages.

Functional Diversity of Forest Plant Communities at Different Restoration Stages
As shown in Figure 3, a comparison of functional diversity at different recovery stages showed that FRic and FDis in the ultimate stage tended to increase gradually as recovery proceeded and both FRic (8.20) and FDis (2.35) values at this stage were significantly higher relative to the other stages. FDiv tended to decrease gradually and the FDiv value (0.86) at the herbaceous stage was significantly higher than that recorded during the rest of the stages. We observed no significant differences in FEve among all the stages. RaoQ showed a gradually increasing trend, with a markedly higher value (6.65) at the ultimate stage compared to the remaining stages.
As shown in Figure 3, a comparison of functional diversity at different recovery stages showed that FRic and FDis in the ultimate stage tended to increase gradually as recovery proceeded and both FRic (8.20) and FDis (2.35) values at this stage were significantly higher relative to the other stages. FDiv tended to decrease gradually and the FDiv value (0.86) at the herbaceous stage was significantly higher than that recorded during the rest of the stages. We observed no significant differences in FEve among all the stages. RaoQ showed a gradually increasing trend, with a markedly higher value (6.65) at the ultimate stage compared to the remaining stages.

Changes in Leaf and Soil Nutrient Contents at Different Recovery Stages
As shown in Table 3, soil organic carbon and total nitrogen contents were significantly higher in the TS and CL stages while soil total phosphorus was markedly higher in the CL stage. The highest coefficients of variation of soil C, N, and P were obtained at the HS stage, TS stages. The total nitrogen content in leaves at the HS stage was significantly higher than that at the HE stage, TS stage, and TR stage. There was no significant difference between the SH stage and CL stage. The variation coefficient of leaf carbon content was the highest in the TS stage, while the variation coefficient of leaf nitrogen and phosphorus content was the highest in the SH stage.

Changes in Leaf and Soil Nutrient Contents at Different Recovery Stages
As shown in Table 3, soil organic carbon and total nitrogen contents were significantly higher in the TS and CL stages while soil total phosphorus was markedly higher in the CL stage. The highest coefficients of variation of soil C, N, and P were obtained at the HS stage, TS stages. The total nitrogen content in leaves at the HS stage was significantly higher than that at the HE stage, TS stage, and TR stage. There was no significant difference between the SH stage and CL stage. The variation coefficient of leaf carbon content was the highest in the TS stage, while the variation coefficient of leaf nitrogen and phosphorus content was the highest in the SH stage. Table 3. Characteristics of leaf and soil C, N, and P contents at different recovery stages. Different lowercase letters represent significant differences. p < 0.05; CC-organic carbon content, NC-total nitrogen content, PC-total phosphorus content; CV(cc)-coefficient of variation of organic carbon, CV(nc)-coefficient of variation of total nitrogen, CV(pc)-coefficient of variation of total phosphorus; CL-ultimate stage, TR-tree stage, TS-tree and shrub transition stage, SH-shrub stage, HS-grass and shrub transition stage, HE-herbaceous stage.

Leaf and Soil Eco-Ecological Stoichiometry Characteristics at Different Restoration Stages
Soil C:N ratio was significantly higher in the shrub stage relative to all the other stages ( Figure 4). Maximal soil C:P was obtained in the tree and shrub transition stage, which was markedly higher than that in the shrub, tree, and ultimate stages. Soil N:P was maximum in the Tree and shrub transition stage and the difference relative to its value during the other stages was significant. Leaf C:N reached a maximum value in the tree and shrub transition stage and was significantly higher relative to the grass and shrub transition stage. Leaf C:P was highest in the herbaceous stage to a significant extent compared to grass and shrub transition and tree stages. Maximal leaf N:P value was reached in the ultimate stage, which was significantly higher compared to tree and shrub transition and tree stages.

Leaf and Soil Eco-Ecological Stoichiometry Characteristics at Different Restoration Stages
Soil C:N ratio was significantly higher in the shrub stage relative to all the other stages ( Figure 4). Maximal soil C:P was obtained in the tree and shrub transition stage, which was markedly higher than that in the shrub, tree, and ultimate stages. Soil N:P was maximum in the Tree and shrub transition stage and the difference relative to its value during the other stages was significant. Leaf C:N reached a maximum value in the tree and shrub transition stage and was significantly higher relative to the grass and shrub transition stage. Leaf C:P was highest in the herbaceous stage to a significant extent compared to grass and shrub transition and tree stages. Maximal leaf N:P value was reached in the ultimate stage, which was significantly higher compared to tree and shrub transition and tree stages.

Network Relationships between the Functional Diversity of Plant Communities and Leaf and Soil Chemical Properties
The correlation network analysis of functional diversity and leaf-soil stoichiometry is presented in Figure 5A. Measurements with p < 0.05 were connected by correlation line segments, with the red lines representing positive correlations and blue lines negative correlations. At the same time, measurements were divided into three major categories, specifically, functional diversity, leaf ecological stoichiometry, and soil chemical properties, and the correlations among these parameters were determined ( Figure 5B). The size of the circular bubbles indicates the number of correlation lines between the three major categories and the specific numerical values are labeled. The overall number of correlation line segments was 43. The correlation segments between functional diversity and soil stoichiometry were the most numerous (10 lines, accounting for 23.26%). Eight correlation lines between functional diversity and leaf ecological stoichiometry were observed (accounting for 18.60%). The correlation lines between related measurements were 21 (7 + 9 + 5), accounting for 48.84%. (Related measurement data are data that represent the same type of index composition, for example: SCC, SNC, SPC, SCN, SCP, SNP is a class of data; FEve, FRic, FDiv, FDis, RaoQ is a class of data; LCC, LNC, LPC, LCN, LCP, LNP is a class

Network Relationships between the Functional Diversity of Plant Communities and Leaf and Soil Chemical Properties
The correlation network analysis of functional diversity and leaf-soil stoichiometry is presented in Figure 5A. Measurements with p < 0.05 were connected by correlation line segments, with the red lines representing positive correlations and blue lines negative correlations. At the same time, measurements were divided into three major categories, specifically, functional diversity, leaf ecological stoichiometry, and soil chemical properties, and the correlations among these parameters were determined ( Figure 5B). The size of the circular bubbles indicates the number of correlation lines between the three major categories and the specific numerical values are labeled. The overall number of correlation line segments was 43. The correlation segments between functional diversity and soil stoichiometry were the most numerous (10 lines, accounting for 23.26%). Eight correla-tion lines between functional diversity and leaf ecological stoichiometry were observed (accounting for 18.60%). The correlation lines between related measurements were 21 (7 + 9 + 5), accounting for 48.84%. (Related measurement data are data that represent the same type of index composition, for example: SCC, SNC, SPC, SCN, SCP, SNP is a class of data; FEve, FRic, FDiv, FDis, RaoQ is a class of data; LCC, LNC, LPC, LCN, LCP, LNP is a class of data. The data that make up these three broad categories become measurement data because they are obtained in actual measurements or calculations). specifically, functional diversity, leaf ecological stoichiometry, and soil chemical properties, and the correlations among these parameters were determined ( Figure 5B). The size of the circular bubbles indicates the number of correlation lines between the three major categories and the specific numerical values are labeled. The overall number of correlation line segments was 43. The correlation segments between functional diversity and soil stoichiometry were the most numerous (10 lines, accounting for 23.26%). Eight correlation lines between functional diversity and leaf ecological stoichiometry were observed (accounting for 18.60%). The correlation lines between related measurements were 21 (7 + 9 + 5), accounting for 48.84%. (Related measurement data are data that represent the same type of index composition, for example: SCC, SNC, SPC, SCN, SCP, SNP is a class of data; FEve, FRic, FDiv, FDis, RaoQ is a class of data; LCC, LNC, LPC, LCN, LCP, LNP is a class of data. The data that make up these three broad categories become measurement data because they are obtained in actual measurements or calculations).

Correlation Analysis of Functional Diversity of Plant Communities with Leaf and Soil Chemical Properties
The indices composed of functional diversity and leaf and soil chemical properties were further analyzed to determine the correlation between functional diversity and leafsoil chemical properties, as shown in Figure 6. Linear plots of the correlation between functional diversity and leaf and soil properties in Figure 7. Soil total phosphorus (SPC) content was significantly correlated with quadratic entropy index, functional divergence, functional separation, and functional richness. The functional dispersion index showed the highest correlation with leaf and soil stoichiometry. Six indices RaoQ was significantly negatively correlated with SCP and LCP, significantly positively correlated with SPC, LPC, and SCC, and had the largest positive correlation with SPC and the largest negative correlation with LCP. FDiv was significantly positively correlated with LCP, significantly negatively correlated with LNC, SPC, and LPC, and showed the largest positive correlation with LCP and largest negative correlation with SPC. FDis was significantly negatively correlated with LCP and SCP, significantly positively correlated with SCC, SNC, SPC, and LNC, and showed the largest positive correlation with SPC and largest negative correlation with LCP. FRic was positively correlated with SCC, SNC, and most significantly, SPC. correlation with LCP. FDiv was significantly positively correlated with LCP, significantly negatively correlated with LNC, SPC, and LPC, and showed the largest positive correlation with LCP and largest negative correlation with SPC. FDis was significantly negatively correlated with LCP and SCP, significantly positively correlated with SCC, SNC, SPC, and LNC, and showed the largest positive correlation with SPC and largest negative correlation with LCP. FRic was positively correlated with SCC, SNC, and most significantly, SPC. Figure 6. Correlations of functional diversity with leaf and soil chemical properties. SCC-soil organic carbon content, SNC-soil total nitrogen content, SPC-soil total phosphorus content, LCCleaf organic carbon content, LNC-leaf total nitrogen content, LPC-leaf total phosphorus content, SCN-soil C:N, SCP-soil C:P, SNP-soil N:P, LCN-leaf C:N, LCP-leaf C:P, LNP-leaf N:P. FRic-Functional richness index, FEve-Functional evenness index, FDiv-Functional divergence index, FDis-Functional separation index, RaoQ-Quadratic entropy coefficient. *, p < 0.05, significant correlation; **, p < 0.01, highly significant correlation; ***, p < 0.001, highly significant correlation. The amount of colored area filled within the small circle represents the correlation size.

Response of Community-Weighted Functional Traits to Different Recovery Stages
The recovery of plant communities governs the process of species replacement over time [46]. As restoration proceeds, abiotic conditions gradually improve [47] and the environment favors species or traits suitable for survival and maintenance in the community [48], resulting in changes in plant community composition and configuration. Our experiments showed that karst plant communities respond differently to weighted functional traits at each restoration stage (Figure 2). Community restoration was not achieved by actions on each trait but through combined effects of a series of traits. With progressive restoration, PLH gradually increased, and community-weighted functional traits showed a trend of decreasing LA and increasing LDMC, consistent with the trend of functional traits at the species level [49][50][51]. Both LDMC and SLA reached maximal values at the ultimate stage, indicating that with recovery, the environment becomes less stressful to the plant community with increasing soil organic matter and soil moisture during succession [52].

Changes in Functional Diversity Patterns at Different Restoration Stages
The main objective of studying the changing characteristics of functional diversity at different restoration stages in karst areas is to explore the occupation of ecological niche space by plant communities and the distribution pattern of functional traits in ecological niche space in different restoration stages [44,53]. The functional richness index is generally positively correlated with species richness [54,55], reflecting the occupancy of niche space by existing species. Higher richness is associated with a more fully occupied niche space, productive community, and stable ecosystem function [7]. In our study, the functional richness (FRic) increased gradually during the restoration process ( Figure 3A), indicating a progressive increase in the occupation of resources by karst plant communities. Our data suggest that plant communities can make fuller use of ecological space in the later stages of recovery. The functional richness of the ultimate stage was significantly higher than that of the other stages, which mostly comprised shade-tolerant species with low light compensation points, high saturation points, and high drought tolerance. The habitat gradually changed from a drastic heterogeneous environment to a moderate mesophytic environment [56], implying that functional richness has a certain response to habitat changes. Functional evenness measures the distribution of species traits within the occupied trait space and indicates the degree of resource use by the community [57,58] and functional uniformity. When the distribution of community traits is high, the distribution of community traits is more regular and uniform, and the community makes full use of resources. Conversely, low FEve values indicate that species and abundance in the community are in scattered clusters within the trait space [43]. The functional evenness (FEve) estimates in this study ( Figure 3B) were not significantly different during the course of restoration, indicating that the efficiency of resource use by plant communities in karst areas did not vary extensively over this time. Functional dispersion reflects the overall dispersion degree of community traits, and the degree of dispersion is inversely related to the strength of interspecific resource competition. This metric is measured by a combination of the degree of community resource variation, degree of competition, and dominance of extreme species in the community, and only has an effect on the weighted average of functional traits, independent of species richness. Higher functional dispersion increases the ecosystem function of communities and leads to more efficient use of resources while lower functional dispersion implies lower differentiation of community niches, underutilization of resources, and intense competition for resources [7,43]. Our data showed significantly higher FDis and RaoQ values in the ultimate stage (CL) compared with the earlier stages ( Figure 3D,E), clearly suggesting that the degree of ecological niche differentiation gradually decreased and the intensity of resource competition weakened as restoration proceeded. The degree of functional separation in an ecosystem is inversely proportional to the degree of niche differentiation and positively proportional to the effect of niche overlap and the degree of resource competition [59]. In the current study, FDiv showed a gradual decrease, with a markedly lower value in the ultimate stage relative to earlier stages ( Figure 3C). Accordingly, we propose that during the restoration of karst plant communities, habitat resources are gradually enriched, and species types increase while the degree of ecological niche differentiation decreases over time. Simultaneously, the interspecific resource competition of plant communities gradually decreases, the ecological niche overlap increases, the functional representation of species in the community decreases, and the functional redundancy increases. Although our metrics do not define functional redundancy, FDiv can represent functional overlap. When FDiv decreases, the overlap effect between functions increases. The more functions overlap, the more redundancy exists. As a result, functional redundancy increases.

Response of Community Functional Diversity to Leaf and Soil Chemical Properties
Exploration of the relationship between functional diversity of plant communities and ecological stoichiometry of leaf and soil can provide information on the correlation between ecological niche differences of functional attributes in ecosystems and plant and soil nutrients. In this study, soil C, N, and P contents gradually increased with recovery progress while leaf C, N, and P showed a trend of increase followed by a decrease (Table 3), indicating that karst areas maintain the rise of plant communities by improving soil nutrient content in the process of positive recovery. In addition, both leaf C:N and C:P were greater than global levels [60] (C:N = 22, C:P = 232) (Figure 4), suggesting that the karst plant communities are less efficient in utilizing N and P [61]. The soil C:N value showed an increase and subsequent decrease as restoration proceeded, indicating a higher rate of organic matter mineralization [62]. Soil C:P was significantly lower than the national level [63] (C:P = 136), signifying a lack of P content available for uptake in the study area.
Our experiments showed that functional diversity was correlated to a greater extent with soil chemical properties than leaf ecological stoichiometry ( Figure 5). RaoQ, FDiv, FDis, and FRic levels were highly correlated with soil phosphorus content (SPC) and leaf C: P (LCP) (Figures 6 and 7). Soil P content is mainly controlled by the weathering stage of soil-forming parent material [63], which, in turn, is influenced by a combination of different climates, topography, vegetation, microorganisms, and human activities. The study area is developed from a typical karst landscape with a shallow soil surface under a humid climate, where precipitation reduces the P content in the soil through surface runoff and leaching. The significance of plant C:P can be summarized into two aspects: (1) to characterize the efficiency of plant P uptake and C assimilation [64] and (2) to reflect the availability of soil P. In addition, the correlation between soil chemical properties and functional diversity was stronger than that of leaf ecological stoichiometry. The possible reasons underlying this finding are as follows: first, ecological stoichiometry characteristics of soils are influenced by species, restoration stage, and habitat [65] and functional diversity represents the magnitude and range of functional traits while soil moisture, total N content, bulk weight, and organic matter play major roles in functional trait alterations during restoration [52]. Second, plant communities were affected by habitat heterogeneity during restoration and soil C, N, and P contents were low in the early recovery period. Plants occupied ecological space through intense resource competition and ecological niches were severely differentiated, showing efficient use of P elements. The C:P content in this study (689) was significantly higher than the global average (232) [60], suggesting higher utilization of plant P within the study area [66]. The relative scarcity of nutrient elements in soil samples from the karst region has been documented, indicating that nutrient utilization can be improved when soil nutrient elements are scarce as a survival strategy for plants to adapt to nutrient-poor conditions [67]. Therefore, the P required for plant growth and development is more scarce, which may be one of the reasons for the higher P utilization efficiency of plants in this region. Overall, as restoration proceeds, the habitat gradually turns homogeneous and the plant community predominantly shows a decrease in ecological niche differentiation, stronger functional overlap effect, and an increase in functional redundancy.

Conclusions
Our collective findings demonstrate that leaf and soil chemical properties and functional diversity interact through a complex network. As the functional traits of plant communities shift from "acquisitive" to "conservative" [52], functional diversity adapts through changes in the level of ecological niche differentiation. Karst plant communities evolve to reduce interspecific resource competition with progressive restoration, thereby increasing the functional overlap effect.