Changes in Soil C, N, and P Concentrations and Stoichiometry in Karst Trough Valley Area under Ecological Restoration: The Role of Slope Aspect, Land Use, and Soil Depth

: This study aims to investigate the roles of slope aspect, land use and soil depth in altering the soil organic carbon (C), total nitrogen (N), and total phosphorus (P) traits in the karst trough valley area experiencing extensive ecological restoration. A total of 54 soil samples were collected at 0–10, 10–20, and 20–30 cm soil depths from secondary forest, plantation forest, and grassland on the relatively more shaded east-facing slope and the contrasting west-facing slope, respectively. The independent and interactive effects of slope aspect, land use, and soil depth on soil C, N, and P concentrations and stoichiometry were determined. The results show that soil C and N concentrations were markedly higher on the east-facing slope than on the west-facing slope, and soil P concentrations showed an opposite trend, leading to signiﬁcant differences in soil C:P and N:P but not in C:N ratios between the two aspects. Soil C and N concentrations were not affected by land use, and soil P concentration was signiﬁcantly higher in plantation forest than in secondary forest and grassland. Soil C and N concentrations signiﬁcantly decreased with increasing soil depth, but soil P concentration presented no signiﬁcant changes with soil depth. Both the land use and soil depth did not differ in terms of their elemental stoichiometry. There were no signiﬁcant interactive effects of slope aspect, land use and soil depth on soil C, N, and P traits. Our results indicate that soil C, N, and P changes are more sensitive to slope aspect rather than land use and soil depth in the karst trough valley area under ecological restoration.


Introduction
The biogeochemical cycles of carbon (C), nitrogen (N), and phosphorus (P) are usually interconnected through their coupled responses to primary production, respiration, and decomposition [1,2]. The C:N:P stoichiometric ratio is the significant indicator for characterizing the nutrient cycling, and its shifts affect the capacity to regenerate the natural vegetation after the anthropogenic abandonment [3][4][5]. Recently, it has been suggested that C, N, and P cycles vary greatly in soils because of the effects by various factors such as topography, vegetation types, and human activities [6,7]. Due to the effects of those influencing factors, soil C, N, and P concentrations are markedly changed, and thus their coupled stoichiometric ratios are probably disrupted, indicating the potential differences occurring in the ecological interactions and processes in different soil ecosystems [8][9][10].
Slope aspect is widely regarded as the crucial topography factor altering soil C, N, and P concentrations and stoichiometry [11][12][13] as a result of the distinct microclimates in different aspects [14]. Nevertheless, the effects of slope aspect on C, N, and P traits are still unclear. For example, Sidari et al. [15] found lower soil C and N contents on the south-facing slope rather than the relatively more shaded northern aspect, and the C:N ratio declined dramatically along the soil profile on the north-facing aspect, but did not significantly differ on the south-facing slope in a Pinus laricio forest ecosystem

Study Area
The study area is situated in the Longtan Trough Valley (LTV) (108.96 • E, 29.01 • E), Youyang Tujia, and Miao Autonomous County, Chongqing Municipality, Southwest China ( Figure 1). The regional climate is characterized by a subtropical monsoon humid climate with a mean annual air temperature of 14.6 • C. The average annual precipitation is 1200 mm, with over 80% occurring in the wet season from April to October [29].

Study Area
The study area is situated in the Longtan Trough Valley (LTV) (108.96° E, 29.01° E), Youyang Tujia, and Miao Autonomous County, Chongqing Municipality, Southwest China ( Figure 1). The regional climate is characterized by a subtropical monsoon humid climate with a mean annual air temperature of 14.6 °C. The average annual precipitation is 1200 mm, with over 80% occurring in the wet season from April to October [29]. LTV is located on the southeast wing of the Tongmaling anticline, where its axial area is dominated by Cambrian dolomite and Ordovician limestone, and the areas in the two wings of the anticline are Silurian sand shale and Permian and Triassic limestone, respectively [30]. LTV runs from northeast to southwest, with a total length of approximately 35 km and an area of 1677 km 2 , which is dominated by karst landform accounting for 76% of the total area [30]. The bottom area of LTV has a width ranging from 1 to 4 km with an elevation of 300-500 m above sea level, and the elevation of east to west wings of LTV varies from 900 to 1100 m above sea level with the slope gradient of 15-30° [30]. The bedrock stratum tendency on either side of the LTV is the same, while the slope aspect is the opposite. The relatively more shaded east-facing slope of LTV has the same dip direction with the bedrock stratum tendency, while the west-facing slope of LTV has inverse dip direction with the bedrock stratum tendency [31]. Since the severe water and soil losses usually occur on the slope presenting the same dip direction with bedrock stratum tendency [31], the eastfacing slope has a higher rock bareness rate than that in the west-facing slope, which together results in severe rocky desertification accounting for 11% of the total area.
The overlying soils in the study area are mainly classified as yellow calcareous soil corresponding to mollic Inceptisols according to the United States Department of Agriculture (USDA) Taxonomy [32]. Soils in both east-facing and west-facing hillslopes are discontinuously distributed within the cavities and grooves between outcrop rocks, and the thickness of soil layers is heterogeneous and not in excess of 60 cm [30]. Although the study area is attributed to subtropical evergreen broadleaf forest zone, most native plants on both east-facing and west-facing hillslopes were removed by deforestation to develop LTV is located on the southeast wing of the Tongmaling anticline, where its axial area is dominated by Cambrian dolomite and Ordovician limestone, and the areas in the two wings of the anticline are Silurian sand shale and Permian and Triassic limestone, respectively [30]. LTV runs from northeast to southwest, with a total length of approximately 35 km and an area of 1677 km 2 , which is dominated by karst landform accounting for 76% of the total area [30]. The bottom area of LTV has a width ranging from 1 to 4 km with an elevation of 300-500 m above sea level, and the elevation of east to west wings of LTV varies from 900 to 1100 m above sea level with the slope gradient of 15-30 • [30]. The bedrock stratum tendency on either side of the LTV is the same, while the slope aspect is the opposite. The relatively more shaded east-facing slope of LTV has the same dip direction with the bedrock stratum tendency, while the west-facing slope of LTV has inverse dip direction with the bedrock stratum tendency [31]. Since the severe water and soil losses usually occur on the slope presenting the same dip direction with bedrock stratum tendency [31], the east-facing slope has a higher rock bareness rate than that in the west-facing slope, which together results in severe rocky desertification accounting for 11% of the total area.
The overlying soils in the study area are mainly classified as yellow calcareous soil corresponding to mollic Inceptisols according to the United States Department of Agriculture (USDA) Taxonomy [32]. Soils in both east-facing and west-facing hillslopes are discontinuously distributed within the cavities and grooves between outcrop rocks, and the thickness of soil layers is heterogeneous and not in excess of 60 cm [30]. Although the study area is attributed to subtropical evergreen broadleaf forest zone, most native plants on both east-facing and west-facing hillslopes were removed by deforestation to develop the local economy over the past decades. Recently, in response to the local policy of ecological restoration, the arable lands have been converted to the plantation forest (PF) and grassland (GL) through revegetation since 2014. The PF is dominated by Zanthoxylum bungeanum Maxim. with a stand density of about 1100 tree ha −1 , and the GL is dominated by herbs such as Oxalis corniculata (L.) and Eleusine indica (L.) Gaertn. Additionally, the secondary forest (SF) has also been found on the hillslopes since 1990s and is dominated by macrophanerophytes, such as Pinus massoniana (L.), Cunninghamia Lanceolata (L.) Hook (Table 1). Approximately, 500-1650 kg ha −1 water-soluble compound fertilizers with high potassium sulfate concentration are applied to the PF in February and July in each year, but human management activities are rarely conducted in SF and GL, according to the local farmers.

Soil Sampling and Measurements
The soil sampling was conducted in the middle of June 2018, at a time in which there was no rain for a week before the sampling days. In situ, three land uses including SF, PF, and GL were all identified on the east-facing and west-facing slopes, where the slope gradients varied between 23-28 degrees, respectively ( Figure 1). For each land use, three sampling sites (10 m × 10 m for each site) in the top, middle and foot slopes were randomly chosen, respectively. Five soil sampling points in each site were assigned by drawing an 'S' shape and soil cores were collected from 0-10, 10-20, and 20-30 cm soil depths, respectively, using a 5 cm diameter soil auger in each point. These soil cores sampled from the five sampling points in each soil depth of each site corresponded to each land use, and were first mixed thoroughly and then subsampled for a composite sample weighing approximately 1 kg. In total, 54 soil samples (2 slope aspects × 3 land uses × 3 sampling sites × 3 soil depths) were thus collected. These soil samples were air-dried and passed through a 2 mm sieve to remove the roots, rocks and other debris. Additionally, undisturbed soil cores corresponding to a 0-20 cm soil depth were collected to determine soil bulk density (BD) and porosity.
Soil organic C concentration was determined by the Walkley-Black dichromate oxidation method [33]. Soil N concentration was measured by the Kjeldahl digestion procedure. Soil P concentration was determined using sodium hydroxide fusion combined with molybdenum antimony colorimetry. Soil pH (1:5 soil/water ratio) was measured with a pH meter. Fraction size distributions of soil micro-aggregates and mass percentages of soil sand, silt, and clay were analyzed using the pipette method. Soil BD was determined using the weighing method. Soil porosity was determined after complete saturation by capillarity. The measurements of the above variables (i.e., soil N, P, pH, micro-aggregates, sand, silt, clay, BD, and porosity) followed the standard Chinese protocols [34]. Moreover, the elevation gradient usually plays an important role in altering soil C, N, and P concentrations and stoichiometry as well as other soil physicochemical properties [35]. Nevertheless, in this study, the elevation gradients for sampling sites were the same between land uses for each slope aspect, and were not in excess of 60 m between slope aspects for each land use; the elevation effects on analyzed soil properties are negligible.

Statistical Analysis
A Kolmogorov-Smirnov (K-S) test was conducted to examine the normality of variables and the log-transformation (base 10) was performed when the data were not normally distributed. A one-way analysis of variance (ANOVA) with least significant difference (LSD) test was used to determine the differences in C, N, and P concentrations and stoichiometry, and selected soil physicochemical variables among land uses and soil depths, respectively. The Independent Samples T test was performed to examine those differences between slope aspects. Multiple-way ANOVA was adopted to examine the effects of slope aspect, land use, soil depth and their possible interactions on soil variables. The mean values calculated from measured data (27 samples for each slope aspect, 18 samples both for each land use and soil depth) for the study soil variables were displayed visually. Pearson correlation analysis was employed to disclose the correlations between study soil variables. Simple regression analysis was applied to determine the relationships among C, N, and P concentrations and their stoichiometric ratios, and the mean values calculated from three soil depths of each land use were presented in terms of different slope aspects. A significance level of p < 0.05 was specified in this study. All statistical analyses abovementioned were implemented using SPSS 22.0 (IBM, Armonk, New York, NY, USA).

Soil C, N, and P Concentrations
Soil C, N, and P concentrations were all significantly affected by slope aspect (Figure 2a,d,g). Specifically, soil C and N concentrations were 57.6% and 61.3% higher on the east-facing slope than those on the west-facing slope, while soil P concentration was 40.7% lower on the east-facing slope relative to the west-facing slope.
Soil C and N concentrations were not significantly affected with land use (Figure 2b,e). However, soil P concentration was significantly higher in PF than in SF and GL (Figure 2h).
Soil C and N concentrations varied among soil depths, and generally decreased with increasing soil depths (Figure 2c,f). Soil P concentrations did not vary significantly with soil depth (Figure 2i). There were no significant interactive effects of slope aspect, land use, and soil depth on C, N, and P concentrations (Supplementary Materials Table S1).

Figure 2.
Soil C (a-c), N (d-f), and P (g-i) concentrations affected by slope aspect, land use and soil depth, respectively. Error bars represent standard deviation of means. "**" marked in the middle indicates significant difference between slope aspects. Different lowercase letters denote significant differences among land uses and soil depths, respectively. SF, secondary forest; PF, plantation forest; GL, grassland.

Soil C, N, and P Stoichiometry
Soil C:N ratio was not significantly affected by slope aspect, while soil C:P and N:P ratios were both 4-fold greater on the east-facing slope than those on the west-facing slope (Figure 3a,d,g). Soil C:N, C:P, and N:P ratios were not affected by land use nor by soil depth (Figure 3b,c,e,f,h,i). There were no significant interactive effects of slope aspect, land use, and soil depth on C:N, C:P, and N:P ratios (Supplementary Materials Table S1). ), and N:P (g-i) ratios affected by slope aspect, land use and soil depth, respectively. Error bars represent standard deviation of means. "**" marked in the middle indicates significant difference between slope aspects. "ns" represents no significance. Same lowercase letters denote no significant difference among land uses and soil depths, respectively. SF, secondary forest; PF, plantation forest; GL, grassland.

Selected Soil Physical and Chemical Variables
Soil pH, 1000-250 µm, 250-50 µm, and <50 µm fractions differed significantly with slope aspect (Figure 4a,d,g,j). While soil pH, 1000-250 µm, and <50 µm fractions were affected with land use, and 250-50 µm fractions were not (Figure 4b,e,h,k). Soil pH and three micro-aggregate fractions were not affected by soil depth (Figure 4c,f,i,l). There were no significant interactive effects of slope aspect, land use, and soil depth on soil pH and three micro-aggregate fractions (Supplementary Materials Table S1).
Higher soil BD and lower porosity were found on the west-facing slope compared to those of the eastern aspect, while sand, silt and clay percentages did not differ with slope aspect (Table 2, Figure 5). These variables were also not significantly affected with land use, though soil BD and porosity did differ significantly between PF and GL, and they were not significantly different from SF on the west-facing slope (Table 2, Figure 5). There were no significant interactive effects of slope aspect and land use on soil BD, porosity, sand, silt, and clay percentages ( Table 2).  (g-i), and <50 µm (j-l) affected by slope aspect, land use, and soil depth, respectively. Error bars represent standard deviation of means. "**" marked in the middle indicates significant difference between slope aspects. Different lowercase letters denote significant differences among land uses and soil depths, respectively. SF, secondary forest; PF, plantation forest; GL, grassland.

Relationships between Selected Soil Physicochemical Variables and C, N, and P Traits
Soil pH did not vary significantly with C, N, and P traits ( Table 3). The 1000-250 µm fraction was significantly positively correlated with C, N concentrations and C:P, N:P ratios, and the 250-50 µm fraction was significantly negatively correlated with these variables. The <50 µm fraction exhibited significant negative correlations with C concentration and C:P ratio. Soil BD showed a negative relationship with C concentration, while a positive relationship between porosity and C concentration was also found. Soil clay percentage presented a significant negative relationship with C:N ratio. Further, soil C concentration was positively correlated with soil N concentration on the west-facing slope, but was not on the east-facing slope ( Figure 6). Soil C:P ratio was positively correlated with N:P ratio on both slope aspects.

Effects of Slope Aspect, Land Use, and Soil Depth on C, N, and P Concentrations
Our results indicated that slope aspect significantly affected soil C, N, and P concentrations, showing higher C and N concentrations on the east-facing slope in comparison to the west-facing slope (Figure 2). This was probably because the east-facing slopes received less solar radiation and the darker, moister, and cooler microclimate depressed microbial activities, and thus suppressed the organic matter decomposition, leading to increase in C and N accumulation in soils [13,14,36]. Higher P was found in the relatively more sunlit western slope in this study, which was in contrast to some observations of higher P concentrations in soils of relatively more shaded aspects [12,37], but was similar to other previous findings [38,39]. This suggested that the vegetation on the east-facing slope took up relatively more P for acquiring more abundant above-and/or below-ground biomass [11,40]. Additionally, soil erosion may vary with differences in environmental conditions, such as vegetation density and runoff discharge, in distinct slope aspects [13,41], which could induce losses in C and nutrient elements, and thus cause significant differences in C, N, and P concentrations between slope aspects [42].
Land use did not significantly alter the soil C and N concentrations, supporting the previous findings [10,18,20]. This was probably due to the overlapping of plant species among SF, PF, and GL, resulting in the marginal differences in soil C and N concentrations [43]. Land use differed significantly soil P concentration, showing that PF exhibited the highest soil P concentration in comparison to SF and GL. The reason for this was that the P inputs through fertilization in PF markedly increased the soil P concentration [10], though the P accumulation by plant communities in SF and GL might also contribute to the increase in soil P concentration [43]. Plant communities regulated the soil C, N, and P concentrations through the plant-soil feedbacks [44], which were potentially affected by the slope aspects [13]. However, the responses of plant communities in different land uses to the slope aspects and their influences on soil properties were not detected in this study, and further exploration on this issue was required. Since the changes in soil P concentration are the comprehensive responses to the dynamics of P inputs and outputs [12], the results indicated the potential in reduction of fertilization in the PF in the course of human management activities. Conversely, Liu et al. [24] found that soil P concentration did not differ significantly between land uses in some specified soil layers in an agro-pastoral region of Northwest China. This resulted from the existence of soluble crystal lattices that could conserve most soil P [45].
Meanwhile, soil C and N concentrations varied significantly with soil depth, showing a decreased trend with increasing soil depth, while soil P concentrations were not significantly differed among soil depths. Luo et al. [46] also found soil C and N concentrations were significant higher in the top soil layer and decreased with increasing soil depth, while soil P concentration did not differ significantly among most soil layers among primary broadleaf forest, plantation, and secondary forest in subtropical China. This was likely due to the greater soil nutrient inputs in the top layer increased the C and N concentrations, and the weathering of parent material caused the distinct variations of soil C, N, and P concentrations in response to the downward profile [18]. Further, there were no significant interactive effects of slope aspect, land use, and soil depth on soil C, N, and P concentrations, possibly suggesting the interaction among these three influencing factors counteracted their independent effects on soil C and nutrient concentrations [13,24].

Effects of Slope Aspect, Land Use, and Soil Depth on C, N, and P Stoichiometry
Soil C:N:P stoichiometry was not significantly affected by independent and interactive slope aspect, land use, and soil depth, except that soil C:P and N:P varied with slope aspect (Supplementary Materials Table S1, Figure 3). Gebrelibanos and Assen [13] reported that soil C:N was not significantly different between slope aspects in the northern highlands of Ethiopia, and the same result was also found in our case. This was probably caused by the proportional shifts in C and N concentrations as these two elements were significantly affected by slope aspect in this study. The proportional shifts in C and N concentrations were also responsible for the slight variation in C:N ratio affected by land use and soil depth [21]. Further, soil P is mainly sourced from weathering of parent material and fertilization, and consumed through plant uptake [12,47], showing the different cycling pathways from soil C and N. This might lead to the decoupling of linkages between soil P and soil C as well as N in the undisturbed sites [9,28], which explained the distinct responses of soil C:P and N:P ratios to slope aspect, land use and soil depth in this study.
Our results showed higher C:N ratio (average 20.4), but lower C:P ratio (average 86) and N:P ratio (average 8.5) in PF soils relative to the averaged ratios (13.4 for C:N, 202.2 for C:P, and 14.1 for N:P) given by Li et al. [28]. This was likely due to the increase in soil P through the fertilization in the PF. Moreover, the current study also found higher C:N ratio (average 21.5), C:P (average 160.5) and lower N:P (average 7.9) in comparison to the Chinese averaged ratios (14.4 for C:N, 136 for C:P, and 9.3 for N:P) in the 0-10 cm soil depth [23]. This suggested that accelerated organic matter accumulation rather than decomposition existed in the study area [48]. Since the N:P ratio was not affected by soil depth, the lower average N:P ratio from three soil depths in comparison to the Chinese average likely meant that the N limitation potentially occurred in the study area, which had been widely reported in other soil ecosystems [5,18,24].

Effects of Slope Aspect, Land Use, and Soil Depth on Selected Physicochemical Properties and Their Linkages with C, N, and P Traits
Soil pH was significantly less on the east-facing slope than on the relatively more sunlit west-facing slope (Figure 4), implying accelerated leaching occurred on the eastern slope [11]. Markedly lower soil BD and higher porosity were recorded on the east-facing slope (Table 2, Figure 5), which was in support of previous findings [13]. This might be a result of higher soil organic carbon on the relatively more shaded slope [49]. Mohammadi and Motaghian [50] found that aggregate organic carbon was the highest in the mesoaggregate fraction (1000-2000 µm), followed by the micro-aggregate fraction (<1000 µm), and was the lowest in the macro-aggregate fraction (>2000 µm). This result supported our findings that the relative bigger micro-aggregate fraction changed consistently with soil organic carbon, showing profoundly higher 1000-250 µm, and lower 250-50 µm and <50 µm fractions on the relatively more shaded east-facing aspect. This was because the degree of stable micro-aggregate commonly played an important and direct role in reflecting the sequestration of soil organic carbon [51]. Weak variations on sand, silt, and clay percentages between slope aspects were possibly attributed to the slight difference in the parent material of these two slopes [52]. Land use remarkably affected soil pH, 1000-250 µm, and <50 µm fractions, but did not affect other selected soil variables, due largely to the different plant species in each land use that shaped distinct soil properties [10,24,43]. Soil depth solely, and interactively with slope aspect and land use, did not significantly affect selected soil physicochemical properties. This was likely due to the soil formation and later development, such as vertical exchanges of soil water, nutrients, and soil organic matter between soil layers, jointly homogenized these soil traits [22], which overrode the effects of slope aspect and land use on soil variables in this study.
Soil C, N, and P concentrations and their ratios showed similar relationships with selected soil physicochemical properties in each slope aspect (e.g., P concentration versus pH, C:P versus pH, and N:P versus pH), while these relationships shifted when two aspects were synthetically assessed (Table 3). Soil C varied significantly with soil N, but not with soil P (Figure 6). This was in good agreement with previous findings given by Cleveland and Liptzin [9], who suggested that soil C and N increasingly decoupled from soil P as soil organic matter accumulated. Meanwhile, Li et al. [28] also demonstrated that accumulation of soil C and N associated with the increasingly conservative P cycling as a consequence of plant uptake during succession or forest development decouple the linkages between C, N, and P, and simultaneously lead to the increases in C:P and N:P ratios. This might explain the significant relationships between C:P and N:P ratios, but not between other stoichiometric ratios in this study ( Figure 6). Additionally, the results were also partially supported by their uniformly close linkages with soil pH, micro-aggregate, BD, and porosity, because soil physical variables are generally associated with chemical indicators due to the processes of soil formation, and thus close links among those properties could be detected in various soil ecosystems [5,36]. These results emphasized that the slope aspect, not land use changes, show a conclusive role in shaping the soil physicochemical features in the study area, as indicated by stronger variations in above relationships among soil physicochemical properties affected by slope aspect. While the land use changes have been demonstrated to strongly alter the soil physical and chemical linkages in other studies [18,24], more data are required to examine the effects of land use changes on these interrelationships among soil C and nutrients and other physicochemical properties for understanding their roles in affecting soil ecosystem functions [10,46].

Implications
Our results suggest that soil C, N, and P traits, overall, presented marginal differences among three land uses in this study. This was in contrast to the previous studies showing that land use significantly changed the soil C, N, and P concentrations, and thus affected their stoichiometric ratios [17,24,28]. Wang et al. [53] found that soil C and N concentrations were significantly higher in SF than PF and GL, and soil P concentration was not significantly different among these land uses, while C:N:P ratios were significantly different among those three land uses. This was probably attributed to the differences in restoration ages, which could be the major factor affecting the soil C:N:P values [54]. Since soil C, N, and P concentrations and stoichiometry potentially characterized the soil ecosystem structures and functions [55,56], the subtle differences in the C, N, and P changes and interactions and the other soil properties (BD, porosity, sand, silt, clay, and 250-50 µm micro-aggregate) among SF, PF, and GL seemed to indicate the equal effectiveness of these ecological measures for restoring the soil ecosystem in such karst trough valley area. Overall, slope aspect played a more important role in changing soil C, N, and P traits in comparison to land use and soil depth, suggesting that the influences of slope aspect should be strongly taken into account when optimizing the ecological restoration strategies in the karst trough valley area.

Conclusions
Soil C and N concentrations were higher on the east-facing slope than on the relatively more sunlit west-facing slope, while soil P concentration showed the opposite trend, leading to significant differences in soil C:P and N:P but not in C:N between two aspects. Soil C and N concentrations were not significantly affected with land use, while soil P concentration was significantly higher in PF than in SF and GL, resulting from the application of fertilizers. Soil C and N concentrations significantly decreased as soil depth increased, but soil P concentration was not significantly different from soil depth. However, both land use and soil depth did not differ in terms of their elemental stoichiometry. There were no significant interactive effects of slope aspect, land use and soil depth on soil C, N, and P concentrations and their stoichiometric ratios. Soil micro-aggregate fractions were, overall, closely linked to C and nutrient traits in this study. The results indicated that SF, PF and GL displayed similar benefits for restoring the soil ecosystem, and slope aspect had more important effects on soil C, N, and P changes in comparison to land use and soil depth. Our results will potentially provide useful information for optimizing sustainable ecological restoration strategies in karst trough valley area.

Informed Consent Statement: Not applicable.
Data Availability Statement: This study did not report any data.