Front. Plant Sci. Frontiers in Plant Science Front. Plant Sci. 1664-462X Frontiers Media S.A. 10.3389/fpls.2023.1194083 Plant Science Original Research Secondary succession of shrub-herb communities in the hilly area of Taihang Mountain Liu Xiuping 1 † Zhou Wangming 2 † Li Xiaoxin 1 Zhang Yuming 1 Dong Wenxu 1 * 1 Key Laboratory of Agricultural Water Resources, Hebei Key Laboratory of Soil Ecology, Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang, China 2 School of life Sciences, Anqing Normal University, Anqing, China

Edited by: Qing Zhang, Inner Mongolia University, China

Reviewed by: Rentao Liu, Ningxia University, China; Yanjing Lou, Chinese Academy of Sciences (CAS), China

*Correspondence: Wenxu Dong, dongwx@sjziam.ac.cn

†These authors have contributed equally to this work

 08 09 2023 2023 14 1194083 26 03 2023 23 08 2023 Copyright © 2023 Liu, Zhou, Li, Zhang and Dong 2023 Liu, Zhou, Li, Zhang and Dong

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Introduction

To document the successional processes of shrub-herb communities after large-scale human disturbance, and understand how changing environmental conditions affect species replacement in semi-humid hilly areas.

 Methods

Utilizing the established permanent plots in the hilly area of Taihang Mountain, we evaluated temporal patterns of vegetation and soil following grass-to-shrub succession.

 Results and Discussion

Along secondary succession, Vitex negundo var. heterophylla gradually dominated in dry sunny slope and shared the dominance with Leptodermis oblonga in shaded slope. Herbaceous dominant species in shrub-herb communities switched from Themeda japonica, Bothriochloa ischaemum, Artemisia sacrorum, and Cleistogenes chinensis in 1986 census to B. ischaemum and A. sacrorum in 2008 census, but herb was no longer dominant in 2020 census. As succession progresses, species dominance increased while richness decreased generally, and herb cover and aboveground biomass decreased, whereas shrub height, cover, and aboveground biomass increased significantly. Soil organic matter (SOM), total nitrogen (TN), total phosphorus (TP), and total potassium (TK) in topsoil increased significantly while pH declined by 1.04 units over the past three decades. Plant communities transitioned from perennial herbs to shrub-herb and then shrub communities, and V. negundo var. heterophylla dominated in the succession of shrub-herb communities. Climate and soil properties, combined with plant attributes, together drive post-disturbance secondary succession. From a management perspective, the tight coupling between vegetation and soil under local climatic conditions should be considered to improve the fragile ecosystem in the hilly area of Taihang Mountain.

 Vitex negundo var. heterophylla secondary succession community structure soil property Taihang Mountain section-in-acceptance Functional Plant Ecology Introduction

Semi-humid ecosystems (precipitation 400-800 mm yr-1) cover about 6% of the global land area and provide crucial ecosystem services for human society (Calvo et al., 2002a; Wang et al., 2014). However, during several decades prior to the 1970s, large-scale human disturbances, such as deforestation, overgrazing, and excessive reclamation, have severely destroyed tree-dominated ecosystems and converted it into shrub-herb communities (Calvo et al., 2002b; Baudena et al., 2020). Even at deteriorated state, these shrub-herb communities provide at least partial ecosystem functions of primary forests and serve as biodiversity refuges, but their successional trajectories, especially in semi-humid hilly areas, depend on a variety of factors related to climate, topography, soils, and so on (Chua et al., 2013; Lai et al., 2020).

The successional trajectory of shrub-herb communities following disturbance is likely to vary widely across landscape, even among those with similar disturbance histories (Norden et al., 2015). Changing environmental conditions, such as temperature (Dodson and Root, 2013), precipitation (José Vidal-Macua et al., 2017; Liu et al., 2022), topography (Boukili and Chazdon, 2017), and soil properties (Engelbrecht et al., 2007; Liu et al., 2022), often influence community reassembly processes and alter its functions composition (Norden et al., 2015). For semi-humid areas, expected increase in drought may promotes leaf abscission (Nepstad et al., 2007), retards plant growth (José Vidal-Macua et al., 2017), and potentially shifts community composition over time (Baudena et al., 2020; Liu et al., 2022). In particular, severe drought combined with high temperature can potentially impede vegetation post-recovery via alterations in growth pattern (Dodson and Root, 2013; Liu et al., 2022), or even lead to the eventual replacement of established forests by drought-tolerant shrubs communities (José Vidal-Macua et al., 2017; Baudena et al., 2020). In addition, topographic heterogeneity can drive shifts in soil microclimate environments, which may potentially affect community structure, composition, and dynamics (Das et al., 2015; Sanaphre-Villanueva et al., 2016; José Vidal-Macua et al., 2017; Jucker et al., 2018). Plant response to these environmental fluctuations generally varies among climatic regions, functional groups, and ecosystems (Chang and Turner, 2019), so accurately understanding vegetation successional processes will assist in disentangling the mechanisms that drive variability in these ecosystems, and enhance community resilience to climatic change.

The Taihang Mountains run from north to south in northern China, and form a natural boundary between Loess Plateau and North China Plain (Yang et al., 2006). The original vegetation of the region was of broadleaved deciduous forest, but large-scale deforestation (e.g., cutting, tilling, logging, grazing, and so on) before the 1970s has severely destroyed and turned it into degraded shrub-herb community (Pérez-Devesa et al., 2008; Liu et al., 2012; Sun et al., 2014). Starting from 1980s, a series of vegetation restoration projects, such as afforestation, banning grazing, and returning slope cropland forest or grassland, were implemented to control soil erosion, reduce ecological degradation, and promote vegetation recovery (Liu et al., 2014; Guo et al., 2021). In addition, annual mean air temperature across the study area exhibited a gradually increasing trend (Liu et al., 2023). Therefore, studies on species composition and community structure of shrublands on time-scales can help explore plant defense responses to natural stresses, predict the potential direction of vegetation succession, and inform appropriate environmental management strategies (Xie and Tang, 2017).

To contribute to a better understanding the succession processes of shrub-herb communities under climate changes, we established permanent plots in the hilly area of Taihang Mountain, China, and monitored vegetation changes in relation to various environmental conditions. Vegetation and soil census was conducted in 1986 followed by a re-census in 2008 and again in 2020 to address the following issues: (1) quantify how vegetation and soil shifts during secondary succession; (2) examine which factors are closely related to the changes in community structure; and (3) elucidate the physio-morphological adaptations of Vitex negundo var. heterophylla to environmental changes during vegetation succession.

 Materials and methods Site description

This study was conducted in Niujiazhuang Catchment of middle Taihang Mountain, China (114°15′50″ E, 37°52′44″ N). Elevation ranges from 247 to 1040 m a.s.l. and slope varies from 20 to 45°. The climate is temperate continental monsoon, with warm summers and cool winters. Annual precipitation (1962-2020) averages 519 mm, ranging from 200 to 1129 mm, with 74% falling between June and September (Liu et al., 2023). The monthly mean air temperature (1962-2020) was -3.07°C in January and 26.8°C in July (Liu et al., 2023). The study site had a total area of 9.3 km2 and was located in an area of lowland hills with a relatively steep slope.

The soil parent materials in the area are mainly granite, granite porphyry, limestone, sandstone, and shale. Soils in the area are highly-weathered mountainous cinnamon soils (classified as Ustalf), mainly derived from gneissic granitoid, and partly come from limestone and shale (Liu et al., 2014). These soils are deep and well drained on north slopes but skeletal on south slopes, with limited organic matter and low water holding capacity.

Vegetation in the area is a mosaic of shrubs, herbs, plantation, deciduous and coniferous forests, and agricultural crops. Shrub-herb community originated by natural regeneration after the destruction of original forest vegetation, of which the dominant species Leptodermis oblonga is accustomed to growing in shaded slope, whereas V. negundo var. heterophylla and Ziziphus jujuba var. spinosa prefer dry sunny slope (Liu et al., 2014). As affected by climate related environmental change, V. negundo var. heterophylla gradually spread northward to occupy available ecological niches for L. oblonga. Natural grasslands are commonly found on gentle north-facing slopes, after 30 years of natural regeneration, grassland in the area has undergone species replacement, such as Bothriochloa ischaemum and Artemisia sacrorum are still dominant today, while Themeda japonica and Cleistogenes chinensis occur sporadically, but Pennisetum centrasiaticum and Achnatherum extremiorientale were no longer present in 2020 census. Robinia pseudoacacia used to be a common tree species in afforestation programs, and was widely planted on Taihang Mountain from 1980s, but now are sparsely distributed in this area and experienced extensive invasion of V. negundo var. heterophylla through natural succession (Liu et al., 2014).

 Field study

In 1986, 144 2 m × 2 m permanent plots were established along 42 30-m long transects (parallel to the slope) to monitor the post-disturbance succession processes of shrub-herb communities. These plots are stratified to capture changes in aspect, slope position, and soil type across the sites and to obtain a reasonable representation of the entire range (valley to ridge). Within each of these plots, we conducted vegetation and soil census in both 1986 and 2008, and in 2020, all plots recorded were grouped by vegetation type and geographic location and then were classified into V. negundo var. heterophylla shrubland (VS), L. oblonga shrubland (LS), R. pseudoacacia plantation (RP), and grassland (G). Finally, a total of 34 plots (8 VS, 8 LS, 8 RP, and 10 G) with expanded area up to 16 m2 were selected for re-census.

In each plot, the height, cover, and number of each shrub or herb species, and cover and aboveground biomass of all shrubs or herbs were recorded, and topographical factors such as elevation, slope gradient, slope aspect, and slope position, were also recorded for each sampling plot. The aboveground biomass of shrubs and herbs were collected within neighboring areas of each plot. The plant samples were oven-dried at 80°C for at least 72 hours to a constant weight, and then dry matter were weighed to determine total biomass. After biomass harvesting, three soil samples (0-20 cm) in each plot were collected along the transect of neighboring areas using a 5 cm diameter soil auger, all visible soil organisms, stones, and plant debris were removed, and the soils were sieved using a 2 mm sieve before laboratory analysis.

 Soil chemical properties

Soil pH was measured in a soil/water ratio of 1:2.5 using a pH meter (Bao, 2000), soil organic matter (SOM) was analyzed using the Walkley-Black method (Nelson and Sommers, 1996), total nitrogen (TN) was determined by the Kjeldahl method (Bremner, 1996), total phosphorus (TP) was digested by perchloric acid and determined by the molybdate colorimetric method (O’Halloran and Cade-Menun, 2006), and total potassium (TK) were analyzed by flame atomic absorption spectrophotometer.

 Data analysis

To account for the variation in community structure in terms of vegetation type and environmental conditions, all data collected were classified into four groups (VS, LS, RP, and G), and the differences in community compositions and soil properties were compared among plots and across the three census period.

In each plot, the importance value (IV) was calculated as the sum of relative cover, relative height, and relative frequency for herbs and shrubs (Do et al., 2019). Shannon-Wiener diversity index (H, H=-Σpi lnpi , pi is the proportion of individuals found in species i), species richness index (S, the number of species), and Simpson’s dominance index (D, D = 1–Σpi 2) were used to assess the changes in species diversity during succession.

Analysis of variance (ANOVA) and least significant difference (LSD) multiple range tested at an alpha-level of 0.05 were used to compare the difference in species diversity (S, D, and H), community structure (height, cover, and aboveground biomass), and soil properties (SOM, TN, TP, and TK) for the entire 0-20 cm soil layer across plots for census interval.

There is only one weather station in the study area and we could not obtain meteorological data for each plot, we used the equations of Guo (1994) and Yang et al. (2006) for the Taihang Mountains to estimate precipitation and temperature at each sampling plot.

P = 519.23 + 151.62 × ( E / 1000 ) − 43.26 × ( E / 1000 ) 2 T = 15.4 − 0.628 × ( L − 34.7 ) − 0.522 × ( E / 100 )

where P (mm) is precipitation, E (m) is elevation, T (°C) is annual average air temperature, L (°) is latitude. Since the study area is a small watershed, the latitudinal differences between sampling plots and weather station are small, we ignored the influence of latitude on temperature.

Multivariate analysis was performed through CANOCO 5 to identify the relationship between environmental variables (climate, topography, and soil properties) and vegetation growth, and examine which factors are closely related to the changes in community structure. At first, detrended correspondence analysis (DCA) was conducted to determine whether to use linear or unimodal numerical methods, and the length of the first DCA ordination axis we obtained was lower than 3 SD ( Table 1 ), indicating that linear model with redundancy analysis (RDA) was the most appropriate ordination method for direct gradient analysis. Accordingly, a final set of nine environmental (precipitation, temperature, slope aspect, slope position, pH, SOM, TN, TP, and TK) and nine vegetation variables (S, D, H, shrub height, shrub cover, shrub biomass, herb height, herb cover, and herb biomass) were selected to construct RDA analysis.

Redundancy analysis of shrubland from 1986 to 2008.

Axes 1 2 3 4
Lengths of gradient (Checked by DCA) 1.50 0.57 0.31 0.37
RDA
Eigenvalues 0.2238 0.0306 0.0226 0.0078
Explained variation (cumulative) 22.38 25.44 27.69 28.47
Pseudo-canonical correlation 0.6918 0.5184 0.6783 0.3784
Explained fitted variation (cumulative) 77.46 88.06 95.87 98.57
Results Species diversity

Along secondary succession, V. negundo var. heterophylla was the species with greatest important value in VS, RP, and G, and shared the dominance with L. oblonga in LS ( Figure 1 ). Herbaceous dominant species in shrub-herb communities switched from T. japonica, B. ischaemum, A. sacrorum, and C. chinensis in 1986 census to B. ischaemum and A. sacrorum in 2008 census, but herb was no longer dominant in 2020 census ( Figure 1 ).

Changes in importance values of major species with secondary succession. AE, Achnatherum extremiorientale; AHI, Arthraxon hispidus; AHE, Artemisia hedinii; AR, Allium ramosum; AS, Artemisia sacrorum; AT, Adenophora tetraphylla; BI, Bothriochloa ischaemum; CA, Chenopodium album; CC, Cleistogenes chinensis; CH, Carex humilis; CT, Cynanchum thesioides; DC, Dianthus chinensis; DI, Dendranthema indicum; EP, Euphorbia pekinensis; FO, Festuca ovina; IP, Ixeris polycephala; LB, Lespedeza bicolor; LO, Leptodermis oblonga; PC, Pennisetum centrasiaticum; PR, Patrinia rupestris; RC, Rubia cordifolia; SA, Scorzonera austriaca; SC, Salsola collina; SV, Setaria viridis; TJ, Themeda japonica; VN, Vitex negundo var. heterophylla; ZJ, Ziziphus jujuba var. spinosa.

There was no significant difference in diversity (H) of VS, RP, and G among three census periods (p > 0.05) ( Figure 2 ). Species richness (S) generally decreased while dominance (D) increased over the 34-year time period (p< 0.05) ( Figure 2 ).

Changes in species diversity of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA. (A) Species richness index (S); (B) Simpson’s dominance index (D); (C) Shannon-Wiener’s index (H); VS, Vitex negundo var. heterophylla shrubland; LS, Leptodermis oblonga shrubland; RP, Robinia pseudoacacia plantation, G, grassland.

 Community structure

Along secondary succession, shrub height in VS, LS, and RP significantly increased (p< 0.05), whereas herb height remained stable ( Figure 3 ). Shrub cover significantly increased after 34 years of vegetation recovery, whereas herb cover in VS and RP significantly decreased (p< 0.05) ( Figure 4 ). Shrub aboveground biomass significantly increased with succession, whereas herb aboveground biomass in VS, LS, and RP significantly decreased (p< 0.05) ( Figure 5 ). Plant communities transitioned from perennial herbs to shrub-herb and then shrub communities.

Changes in height of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA.

Changes in cover of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA.

Changes in aboveground biomass of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA.

In 2020, community mean height, cover, and aboveground biomass were 2.19m, 97.2%, and 15.7 Mg ha-1, whereas the values of V. negundo var. heterophylla were 2.37m, 78.0%, and 17.9 Mg ha-1, respectively ( Figures 3 – 5 ). V. negundo var. heterophylla was the most dominant species in the hilly area of Taihang Mountain.

 Soil property

Soil pH in topsoil declined significantly over the past three decades (p< 0.05), from 7.64 in 1986 to 6.60 in 2020, with a net decrease of about 1.04 units ( Figure 6 ). Long-term secondary succession increased SOM content in topsoil, on average, by 37.0% ( Figure 7 ). From 1986 to 2020, Surface TN and TP contents in VS, LS, and G increased significantly (p< 0.05), while in RP, the changes tended to be gradual ( Figure 7 ). Surface TK content increased quickly during the first 20 years of vegetation restoration and tended to be stable thereafter ( Figure 7 ).

Changes in soil pH in topsoil (0-20 cm) of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA. VS, Vitex negundo var. heterophylla shrubland; LS, Leptodermis oblonga shrubland; RP, Robinia pseudoacacia plantation; G, Grassland; T, Total.

Changes in soil organic matter (A), total N (B), total P (C), and total K (D) in topsoil (0-20 cm) of shrub-herb communities with succession. Different letters indicate significant differences between census years at p< 0.05 according to ANOVA. VS, Vitex negundo var. heterophylla shrubland; LS, Leptodermis oblonga shrubland; RP, Robinia pseudoacacia plantation; G, Grassland; T, Total.

 Effects of environmental variables

RDA analysis showed that there was a strong correlation between vegetation and environmental variables, with pseudo-canonical correlation of 0.692 on the first axis and 0.518 on the second axis ( Table 1 ). The cumulative percentage variance of vegetation-environment relations of the first and second axis were 77.5 (eigenvalue 0.224) and 88.1% (eigenvalue 0.031), respectively ( Table 1 ). The RDA results also indicated that temperature, TP, TK, and precipitation were the most important factors for community structure, contributing respectively 39.1 (p<0.01), 18.9 (p<0.01), 8.3 (p<0.05), and 7.7% (p<0.05) of the total variation ( Figure 8 , Table 2 ).

Redundancy analysis (RDA) of environmental and vegetation variables from 1986 to 2020. Prec, precipitation; Temp, temperature; Aspe, slope aspect; Posi, slope position; SOM, soil organic matter; TN, total N; TP, total P; TK, total K; S, species richness index; D, Simpson’s dominance index; H, Shannon-Wiener’s index; ShrH, shrub height; ShrC, shrub cover; ShrB, shrub biomass; HerH, herb height; HerC, herb cover; HerB, herb biomass.

Summary of the forward selection procedure in the redundancy analysis (RDA).

Name Explains (%) Contribution (%) F P
Temp 11.3 39.1 9.3 0.002
TP 5.5 18.9 4.7 0.002
TK 2.4 8.3 2.1 0.042
Prec 2.2 7.7 2.0 0.046
pH 2.1 7.2 1.9 0.050
Aspe 1.7 5.9 1.6 0.094
TN 1.5 5.3 1.4 0.094
Posi 1.4 5.0 1.3 0.110
SOM 0.7 2.5 0.7 0.436

Prec, precipitation; Temp, temperature; Aspe, slope aspect; Posi, slope position; SOM, soil organic matter; TN, total N; TP, total P; TK, total K.

 Discussion

V. negundo var. heterophylla is a highly water-consuming plant (Si et al., 2020), how did it become the most dominant species in semi-humid climate zone? Wang et al. (2017) and Zhu et al. (2016) reported that V. negundo var. heterophylla possesses a functionally dimorphic root system, which allows it to use both shallow and deep soil water, or even ground water (Moreira et al., 2003; Damascos et al., 2005). In addition, Du et al. (2010) and Zhang (2007) pointed out the sprouting time of V. negundo var. heterophylla could vary with environmental conditions, they sprouted new leaves in mid-April at mild-moderate drought while in early May at severe drought, but almost stagnant growth at extreme drought. Notably, Yan et al. (2000) found that the numerous epidermal hairs of V. negundo var. heterophylla could shielded most stomata and formed a relatively independent system for reducing water transpiration. Therefore, the high degree of morphological plasticity, conservative light utilization strategy, and strong recovery ability after disturbance make V. negundo var. heterophylla a widely predominant species in warm-temperate regions (Li et al., 2017).

Aboveground biomass accumulation increased during secondary succession ( Figure 5 ) (Castro and Freitas, 2009; Enes et al., 2020; Ma and Wang, 2020). Averaged over all plots, the aboveground biomass of shrubland increased over 5 times ( Figure 5 ), similar to the results of Yang et al. (2017) and Petrie et al. (2015). Despite differences in vegetation type, secondary shrublands following disturbance are substantial carbon sinks and that this capacity to store carbon increases with succession (Adhikari and White, 2016; Kong et al., 2022). In addition, previous studies reported species-rich stands had higher carbon stocks than stands with low richness (Chen et al., 2018; Liu et al., 2018), and Chen et al. (2018) proposed high plant diversity increased above- and belowground biomass as well as the resistance of productivity to climate extremes. While, in our study, the relationship between aboveground biomass and species diversity weakened with succession ( Figures 2 and 5 ), Potter and Woodall (2014) and Marin-Spiotta et al. (2007) also obtained similar results. This demonstrated that a limited number of species well-adapted to the local conditions promoted vegetation recovery and accumulated biomass in semi-humid areas (El-Sheikh, 2005; Potter and Woodall, 2014).

Soil pH in the hilly area of Taihang Mountain declined significantly over the past three decades, with an overall decrease of 1.04 units, similar to those reported by Yang et al. (2012) and Zhang et al. (2022) for China’s soils. Soil acidification is the result of atmospheric deposition, plant growth, and soil forming processes (Ritter et al., 2003; De Schrijver et al., 2006; Hong et al., 2018). Undoubtedly, atmospheric deposition reduced soil pH (De Schrijver et al., 2012; Yang et al., 2015). Whereas, plant metabolism could neutralize soil pH (Ovington, 1953; De Schrijver et al., 2012; Chen et al., 2019), and the minerals/ions released from granite and gneiss through weathering could buffer soil acidification (Rhoades and Binkley, 1996; Drohan and Sharpe, 1997). However, soil pH in our study area declined significantly. Limited precipitation may contribute to the absence of soil animals and microorganisms, which retard litter decomposition and mineral weathering, and consequently result in soil acidification (Shu et al., 2019).

Despite the relatively decrease of soil pH, soil nutrients in topsoil accumulated gradually as succession progresses. We found SOM content increased on average by 37% from 1986 to 2020 ( Figure 7 ), Piche and Kelting (2015) and Sokolowska et al. (2020) obtained similar results. This suggested that organic matter inputs exceeded decomposition outputs over time of secondary succession (Nadal-Romero et al., 2016; Sokolowska et al., 2020). As with SOM, TN content in topsoil significantly increased, similar to the observations of Deng et al. (2014) and Segura et al. (2020). The apparent increase in topsoil TN could be attributed to the release of N from plant residues (Johnson et al., 2016), the return of N from subsoil (Smal et al., 2019), and atmospheric N deposition (Falkengren-Grerup et al., 2006). Similar to SOM and TN, TP and TK contents in the top 20 cm of soil increased progressively with shrub age, this agreed with the results of Du et al. (2007). The commonly observed increase in topsoil P may be achieved by atmospheric P deposition (Cao et al., 2011) and P translocation from depth in the soil profile to the surface soil (Sullivan et al., 2019). Similarly, high levels of K in the 0-20 cm soil may come from the K uptake by roots from deeper soil layers to the topsoil (Segura et al., 2020). Consequently, vegetation succession in the hilly area of Taihang Mountain significantly increased most of the measured soil parameters, and consequently improved soil quality.

Climate varied with altitude results in different vegetation and soil, and soil, in turn, influences plant growth (Miki and Kondoh, 2002; Du et al., 2014). Among the environmental variables, temperature, TP, TK, and precipitation were important factors in explaining variations in community structure of Taihang Mountain. Several reports indicated that precipitation and temperature were the most important factors influencing plant growth, community structure and function (Yang et al., 2006; Nunes et al., 2019; Zhou et al., 2019). For example, Zhou and Zhang (1996) reported that widespread drought driven by high temperatures or low precipitations could lead to substantial forest decline, or replacement of drought-vulnerable with drought-tolerant species (Falk et al., 2019). As temperature rise and water availability decreases (Liu et al., 2023), the more drought-tolerant V. negundo var. heterophylla becomes the most dominant species in the hilly area of Taihang Mountain. As expected, vegetation recovery generally improved soil quality, in turn, soil properties also determine plant community composition (Zhang et al., 2006; Zhao et al., 2017). The increased soil nutrients, in particular P and K, could facilitate shrub vegetation restoration in this region. Thus, to improve and restore the fragile ecological ecosystem in the hilly area of Taihang Mountain, we should consider the tight coupling between vegetation and soil under local climatic conditions.

 Conclusion

Along secondary succession, V. negundo var. heterophylla gradually dominated in dry sunny slope and shared the dominance with L. oblonga in shaded slope. Herbaceous dominant species in shrub-herb communities switched from T. japonica, B. ischaemum, A. sacrorum, and C. chinensis in 1986 census to B. ischaemum and A. sacrorum in 2008 census, but herb was no longer dominant in 2020 census. As succession progresses, species dominance increased while richness decreased generally, and herb cover and aboveground biomass decreased, whereas shrub height, cover, and aboveground biomass increased significantly. Long-term secondary succession increased topsoil SOM, TN, TP, and TK significantly while declined pH by 1.04 units. Plant communities transitioned from perennial herbs to shrub-herb and then shrub communities, and V. negundo var. heterophylla dominated in succession of shrub-herb communities. Climate and soil properties, combined with plant attributes, together drive vegetation post-recovery.

 Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

 Author contributions

XPL and WZ: Investigation, Formal analysis, Writing - Original Draft; XXL and YZ: Methodology, Data Curation, Writing - Review & Editing; WD: Resources, Supervision, Project administration. All authors contributed to the article and approved the submitted version.

 Funding

This research was supported by the National Key Research and Development Program of China (2021YFD1901104, 2022YFD1901604), the “Strategic Priority Research Program” of the Chinese Academy of Sciences (XDA28020303, XDA26040103), the Key Research and Development Program of Hebei Province (22326412D), and the Natural Science Foundation of Hebei Province (D2021503009).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Adhikari A. White J. D. (2016). Climate change impacts on regenerating shrubland productivity. Ecol. Model. 337, 211–220. doi: 10.1016/j.ecolmodel.2016.07.003 Bao S. (2000). Soil and Agricultural Chemistry Analysis (Beijing: China Agriculture Press). Baudena M. Santana V. M. Baeza M. J. Bautista S. Eppinga M. B. Hemerik L. . (2020). Increased aridity drives post-fire recovery of Mediterranean forests towards open shrublands. New Phytol. 225, 1500–1515. doi: 10.1111/nph.16252 Boukili V. K. Chazdon R. L. (2017). Environmental filtering, local site factors and landscape context drive changes in functional trait composition during tropical forest succession. Perspect. Plant Ecol. 24, 37–47. doi: 10.1016/j.ppees.2016.11.003 Bremner J. M. (1996). “Nitrogen-total” in Methods of Soil Analysis. Eds. Sparks D. L. Page A. L. Helmke P. A. Loeppert R. H. Soltanpour P. N. Tabatabai M. A. Johnston C. T. Sumner M.E. (Madison, WI: Soil Science Society of America, American Society of Agronomy), 1085–1121. Calvo L. Tárrega R. de Luis E. (2002a). Secondary succession after perturbations in a shrubland community. Acta Oecol. 23, 393–404. doi: 10.1016/S1146-609X(02)01164-5 Calvo L. Tárrega R. de Luis E. (2002b). The dynamics of mediterranean shrubs species over 12 years following perturbations. Plant Ecol. 160, 25–42. doi: 10.1023/A:1015882812563 Cao Z. Yang Y. Lu J. Zhang C. (2011). Atmospheric particle characterization, distribution, and deposition in Xi’an, Shaanxi Province, Central China. Environ. pollut. 159, 577–584. doi: 10.1016/j.envpol.2010.10.006 Castro H. Freitas H. (2009). Above-ground biomass and productivity in the Montado: From herbaceous to shrub dominated communities. J. Arid. Environ. 73, 506–511. doi: 10.1016/j.jaridenv.2008.12.009 Chang C. C. Turner B. L. (2019). Ecological succession in a changing world. J. Ecol. 107, 503–509. doi: 10.1111/1365-2745.13132 Chen L. Liu L. Qin S. Yang G. Fang K. Zhu B. . (2019). Regulation of priming effect by soil organic matter stability over a broad geographic scale. Nat. Commun. 10, 5112. doi: 10.1038/s41467-41019-13119-z Chen S. Wang W. Xu W. Wang Y. Wan H. Chen D. . (2018). Plant diversity enhances productivity and soil carbon storage. Proc. Natl. Acad. Sci. U.S.A. 115, 4027–4032. doi: 10.1073/pnas.1700298114 Chua S. C. Ramage B. S. Ngo K. M. Potts M. D. Lum S. K. Y. (2013). Slow recovery of a secondary tropical forest in Southeast Asia. For. Ecol. Manage. 308, 153–160. doi: 10.1016/j.foreco.2013.07.053 Damascos M. A. Prado C. H. B. A. Ronquim C. C. (2005). Bud composition, branching patterns and leaf phenology in Cerrado woody species. Ann. Bot. 96, 1075–1084. doi: 10.1093/aob/mci258 Das A. Nagendra H. Anand M. Bunyan M. (2015). Topographic and bioclimatic determinants of the occurrence of forest and grassland in tropical montane forest-grassland mosaics of the Western Ghats, India. PloS One 10, e0130566. doi: 10.1371/journal.pone.0130566 Deng L. Wang K. Shangguan Z. (2014). Long-term natural succession improves nitrogen storage capacity of soil on the Loess Plateau, China. Soil Res. 52, 262–270. doi: 10.1071/sr12377 De Schrijver A. De Frenne P. Staelens J. Verstraeten G. Muys B. Vesterdal L. . (2012). Tree species traits cause divergence in soil acidification during four decades of postagricultural forest development. Global Change Biol. 18, 1127–1140. doi: 10.1111/j.1365-2486.2011.02572.x De Schrijver A. Mertens J. Geudens G. Staelens J. Campforts E. Luyssaert S. . (2006). Acidification of forested podzols in North Belgium during the period 1950–2000. Sci. Total Environ. 361, 189–195. doi: 10.1016/j.scitotenv.2005.06.015 Do H. T. T. Grant J. C. Trinh N. B. Zimmer H. C. Tran L. D. Nichols J. D. (2019). Recovery of tropical moist deciduous dipterocarp forest in Southern Vietnam. For. Ecol. Manage. 433, 184–204. doi: 10.1016/j.foreco.2018.10.062 Dodson E. K. Root H. T. (2013). Conifer regeneration following stand-replacing wildfire varies along an elevation gradient in a ponderosa pine forest, Oregon, USA. For. Ecol. Manage. 302, 163–170. doi: 10.1016/j.foreco.2013.03.050 Drohan J. R. Sharpe W. E. (1997). Long-term changes in forest soil acidity in Pennsylvania, U.S.A. Water Air Soil Poll. 95, 299–311. doi: 10.1023/A:1026469814427 Du N. Guo W. Zhang X. Wang R. (2010). Morphological and physiological responses of Vitex negundo L. var. heterophylla (Franch.) Rehd. to drought stress. Acta Physiol Plant 32, 839–848. doi: 10.1007/s11738-010-0468-z Du F. Liang Z. Xu X. Shan L. Zhang X. (2007). Community biomass of abandoned farmland and its effects on soil nutrition in the Loess hilly region of Northern Shaanxi, China. Acta Ecol. Sin. 27, 1673–1683. doi: 10.1016/S1872-2032(07)60038-9 Du B. Liu C. Kang H. Zhu P. Yin S. Shen G. . (2014). Climatic control on plant and soil δ13C along an altitudinal transect of Lushan mountain in subtropical China: characteristics and interpretation of soil carbon dynamics. PloS One 9, e86440. doi: 10.1371/journal.pone.0086440 El-Sheikh M. A. (2005). Plant succession on abandoned fields after 25 years of shifting cultivation in Assuit, Egypt. J . Arid. Environ. 61, 461–481. doi: 10.1016/j.jaridenv.2004.10.006 Enes T. Lousada J. Fonseca T. Viana H. Calvão A. Aranha J. (2020). Large scale shrub biomass estimates for multiple purposes. Life 10, 33. doi: 10.3390/life10040033 Engelbrecht B. M. J. Comita L. S. Condit R. Kursar T. A. Tyree M. T. Turner B. L. . (2007). Drought sensitivity shapes species distribution patterns in tropical forests. Nature 447, 80–82. doi: 10.1038/nature05747 Falk D. A. Watts A. C. Thode A. E. (2019). Scaling ecological resilience. Front. Ecol. Evol. 7, 275. doi: 10.3389/fevo.2019.00275 Falkengren-Grerup U. ten Brink D. J. Brunet J. (2006). Land use effects on soil N, P, C and pH persist over 40-80 years of forest growth on agricultural soils. For. Ecol. Manage. 225, 74–81. doi: 10.1016/j.foreco.2005.12.027 Guo Y. C. (1994). Calculation methods for precipitation in the Taihang Mountains and the Yanshan Mountains. Geogr. Terr. Res. 10, 35–39. Guo Y. Gheyret G. Liu T. Zhang Y. Kang M. Mohhamot A. . (2021). Distribution patterns and climate limitations of typical shrublands in northern China. Sci. Sin. Vitae 51, 346–361. doi: 10.1360/SSV-2020-0186 Hong S. Piao S. Chen A. Liu Y. Liu L. Peng S. . (2018). Afforestation neutralizes soil pH. Nat. Commun. 9, 520. doi: 10.1038/s41467-018-02970-1 Johnson D. W. Trettin C. C. Todd D. E. (2016). Changes in forest floor and soil nutrients in a mixed oak forest 33 years after stem only and whole-tree harvest. For. Ecol. Manage. 361, 56–68. doi: 10.1016/j.foreco.2015.11.012 José Vidal-Macua J. Ninyerola M. Zabala A. Domingo-Marimon C. Pons X. (2017). Factors affecting forest dynamics in the Iberian Peninsula from 1987 to 2012. The role of topography and drought. For. Ecol. Manage. 406, 290–306. doi: 10.1016/j.foreco.2017.10.011 Jucker T. Bongalov B. Burslem D. F. R. P. Nilus R. Dalponte M. Lewis S. L. . (2018). Topography shapes the structure, composition and function of tropical forest landscapes. Ecol. Lett. 21, 989–1000. doi: 10.1111/ele.12964 Kong L. Lu F. Rao E. Ouyang Z. (2022). Carbon sink under different carbon density levels of forest and shrub, a case in Dongting Lake Basin, China. Remote Sens. 14, 2672. doi: 10.3390/rs14112672 Lai H. R. Chong K. Y. Yee A. T. K. Tan H. T. W. van Breugel M. (2020). Functional traits that moderate tropical tree recruitment during post-windstorm secondary succession. J. Ecol. 108, 1322–1333. doi: 10.1111/1365-2745.13347 Li B. Zeng T. Ran J. Yue B. Zhang M. Shang T. . (2017). Characteristics of the early secondary succession after landslides in a broad-leaved deciduous forest in the south Minshan Mountains. For. Ecol. Manage. 405, 238–245. doi: 10.1016/j.foreco.2017.09.020 Liu X. Trogisch S. He J. Niklaus P. A. Bruelheide H. Tang Z. . (2018). Tree species richness increases ecosystem carbon storage in subtropical forests. Proc. R. Soc B. 285, 20181240. doi: 10.1098/rspb.2018.1240 Liu X. P. Zhang W. J. Hu C. S. Tang X. G. (2014). Soil greenhouse gas fluxes from different tree species on Taihang Mountain, North China. Biogeosciences 11, 1649–1666. doi: 10.5194/bg-11-1649-2014 Liu D. Zhang C. Ogaya R. Estiarte M. Zhang X. Pugh T. A. M. . (2022). Delayed and altered post-fire recovery pathways of Mediterranean shrubland under 20-year drought manipulation. For. Ecol. Manage. 506, 119970. doi: 10.1016/j.foreco.2021.119970 Liu X. Zhang W. Yang F. Zhou X. Liu Z. Qu F. . (2012). Changes in vegetation-environment relationships over long-term natural restoration process in Middle Taihang Mountain of North China. Ecol. Eng. 49, 193–200. doi: 10.1016/j.ecoleng.2012.06.040 Liu X. Zhou W. Wang X. Wu H. Dong W. (2023). Microbial gradual shifts during the process of species replacement in Taihang Mountain. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1158731 Ma X. Wang X. (2020). Biomass partitioning and allometric relations of the Reaumuria soongorica shrub in Alxa steppe desert in NW China. For. Ecol. Manage. 468, 118178. doi: 10.1016/j.foreco.2020.118178 Marin-Spiotta E. Ostertag R. Silver W. L. (2007). Long-term patterns in tropical reforestation: Plant community composition and aboveground biomass accumulation. Ecol. Appl. 17, 828–839. doi: 10.1890/06-1268 Miki T. Kondoh M. (2002). Feedbacks between nutrient cycling and vegetation predict plant species coexistence and invasion. Ecol. Lett. 5, 624–633. doi: 10.1046/j.1461-0248.2002.00347.x Moreira M. Z. Scholz F. G. Bucci S. J. Sternberg L. S. Goldstein G. Meinzer F. C. . (2003). Hydraulic lift in a neotropical savanna. Funct. Ecol. 17, 573–581. doi: 10.1046/j.1365-2435.2003.00770.x Nadal-Romero E. Cammeraat E. Pérez-Cardiel E. Lasanta T. (2016). Effects of secondary succession and afforestation practices on soil properties after cropland abandonment in humid Mediterranean mountain areas. Agr. Ecosyst. Environ. 228, 91–100. doi: 10.1016/j.agee.2016.05.003 Nelson D. W. Sommers L. E. (1996). “Total carbon, organic carbon, and organic matter,” in Methods of Soil Analysis. Eds. Sparks D. L. Page A. L. Helmke P. A. Loeppert R. H. Soltanpour P. N. Tabatabai M. A. Johnston C. T. Sumner M.E. (Madison, WI: Soil Science Society of America, American Society of Agronomy), 539–594. Nepstad D. C. Tohver I. M. Ray D. Moutinho P. Cardinot G. (2007). Mortality of large trees and lianas following experimental drought in an Amazon forest. Ecology 88, 2259–2269. doi: 10.1890/06-1046.1 Norden N. Angarita H. A. Bongers F. Martínez-Ramos M. Granzow-de la Cerda I. van Breugel M. . (2015). Successional dynamics in Neotropical forests are as uncertain as they are predictable. Proc. Natl. Acad. Sci. U. S. A. 112, 8013–8018. doi: 10.1073/pnas.1500403112 Nunes A. Köbel M. Pinho P. Matos P. Costantini E. A. C. Soares C. . (2019). Local topographic and edaphic factors largely predict shrub encroachment in Mediterranean drylands. Sci. Total Environ. 657, 310–318. doi: 10.1016/j.scitotenv.2018.11.475 O’Halloran I. P. Cade-Menun B. J. (2006). “Total and organic phosphorus,” in Soil Sampling and Methods of Analysis. Eds. Carter M. R. Gregorich E. G. (Boca Raton, FL, USA: CRC Press Taylon and Francis), 265–291. Ovington J. D. (1953). Studies of the development of woodland conditions under different trees: I. Soils pH. J. Ecol. 4, 13–34. doi: 10.2307/2257160 Pérez-Devesa M. Cortina J. Vilagrosa A. Vallejo R. (2008). Shrubland management to promote Quercus suber L. establishment. For. Ecol. Manage. 255, 374–382. doi: 10.1016/j.foreco.2007.09.074 Petrie M. D. Collins S. L. Swann A. M. Ford P. L. Litvak M. E. (2015). Grassland to shrubland state transitions enhance carbon sequestration in the northern Chihuahuan Desert. Global Change Biol. 21, 1226–1235. doi: 10.1111/gcb.12743 Piché N. Kelting D. L. (2015). Recovery of soil productivity with forest succession on abandoned agricultural land. Restor. Ecol. 23 (5), 645–654. doi: 10.1111/rec.12241 Potter K. M. Woodall C. W. (2014). Does biodiversity make a difference? Relationships between species richness, evolutionary diversity, and aboveground live tree biomass across US forests. For. Ecol. Manage. 321, 117–129. doi: 10.1016/j.foreco.2013.06.026 Rhoades C. Binkley D. (1996). Factors influencing decline in soil pH in Hawaiian Eucalyptus and Albizia plantations. For. Ecol. Manage. 80, 47–56. doi: 10.1016/0378-1127(95)03646-6 Ritter E. Vesterdal L. Gundersen P. (2003). Changes in soil properties after afforestation of former intensively managed soils with oak and Norway spruce. Plant Soil 249, 319–330. doi: 10.1023/A:1022808410732 Sanaphre-Villanueva L. Dupuy J. M. Andrade J. L. Reyes-García C. Jackson P. C. Paz H. (2016). Patterns of plant functional variation and specialization along secondary succession and topography in a tropical dry forest. Environ. Res. Lett. 12, 055004. doi: 10.1088/1748-9326/aa6baa Segura C. Navarro F. B. Noelia Jimenez M. Fernandez-Ondono E. (2020). Implications of afforestation vs. secondary succession for soil properties under a semiarid climate. Sci. Total Environ. 704, 135393. doi: 10.1016/j.scitotenv.2019.135393 Shu M. Zhao Q. Li Z. Zhang L. Wang P. Hu S. (2019). Effects of global change factors and living roots on root litter decomposition in a Qinghai-Tibet alpine meadow. Sci. Rep. 9, 16924. doi: 10.1038/s41598-019-53450-5 Si M. Cao J. Yang H. Zhu C. (2020). Soil water variation of different vegetation community in Taihang Mountain Area. Chin. J. Eco Agric. 28, 1766–1777. Smal H. Ligeza S. Pranagal J. Urban D. Pietruczyk-Poplawska D. (2019). Changes in the stocks of soil organic carbon, total nitrogen and phosphorus following afforestation of post-arable soils: A chronosequence study. For. Ecol. Manage. 451, 117536. doi: 10.1016/j.foreco.2019.117536 Sokolowska J. Józefowska A. Woźnica K. Zaleski T. (2020). Succession from meadow to mature forest: Impacts on soil biological, chemical and physical properties-Evidence from the Pieniny Mountains, Poland. Catena 189, 104503. doi: 10.1016/j.catena.2020.104503 Sullivan B. W. Nifong R. L. Nasto M. K. Alvarez-Clare S. Dencker C. M. Soper F. M. . (2019). Biogeochemical recuperation of lowland tropical forest during succession. Ecology 100, e02641. doi: 10.1002/ecy.2641 Sun J. Wang X. Cheng G. Wu J. Hong J. Niu S. (2014). Effects of grazing regimes on plant traits and soil nutrients in an alpine steppe, Northern Tibetan Plateau. PloS One 9, e108821. doi: 10.1371/journal.pone.0108821 Wang J. Fu B. Lu N. Zhang L. (2017). Seasonal variation in water uptake patterns of three plant species based on stable isotopes in the semi-arid Loess Plateau. Sci. Total Environ. 609, 27–37. doi: 10.1016/j.scitotenv.2017.07.133 Wang C. Wang X. Liu D. Wu H. Lü X. Fang Y. . (2014). Aridity threshold in controlling ecosystem nitrogen cycling in arid and semi-arid grasslands. Nat. Commun. 5, 4799. doi: 10.1038/ncomms5799 Xie Z. Tang Z. (2017). Studies on carbon storage of shrubland ecosystems in China. Chin. J. Plant Ecol. 41, 1–4. doi: 10.17521/cjpe.2017.0012 Yan C. Han X. Chen L. (2000). The relationship between the ecophysiological feature and leaf characteristics of some woody plants in Beijing mountain zone. Acta Ecol. Sin. 20, 53–60. Yang X. Guo Y. Mohhamot A. Liu H. Ma W. Yu S. . (2017). Distribution of biomass in relation to environments in shrublands of temperate China. Chin. J. Plant Ecol. 41, 22–30. doi: 10.17521/cjpe.2016.0199 Yang Y. Ji C. Ma W. Wang S. Wang S. Han W. . (2012). Significant soil acidification across northern China’s grasslands during 1980s-2000s. Global Change Biol. 18, 2292–2300. doi: 10.1111/j.1365-2486.2012.02694.x Yang Y. Li P. He H. Zhao X. Datta A. Ma W. . (2015). Long-term changes in soil pH across major forest ecosystems in China. Geophys. Res. Lett. 42, 933–940. doi: 10.1002/2014GL062575 Yang Y. Watanabe M. Li F. Zhang J. Zhang W. Zhai J. (2006). Factors affecting forest growth and possible effects of climate change in the Taihang Mountains, northern China. Forestry 79, 135–147. doi: 10.1093/forestry/cpi062 Zhang C. (2007). Study on bioecological characteristics of soil and water conservation shrub-Vitex Negundo (Beijing: Beijing Forestry University). Zhang Q. Jianxing Z. Wang Q.-F. Xu L. Li M. Dai H. . (2022). Soil acidification in China’s forests due to atmospheric acid deposition from 1980 to 2050. Sci. Bull. 67, 914–917. doi: 10.1016/j.scib.2022.01.004 Zhang J. T. Xi Y. Li J. (2006). The relationships between environment and plant communities in the middle part of Taihang Mountain Range, North China. Community Ecol. 7, 155–163. doi: 10.1556/ComEc.7.2006.2.3 Zhao H. Wang Q.-R. Fan W. Song G.-H. (2017). The relationship between secondary forest and environmental factors in the southern Taihang Mountains. Sci. Rep. 7, 16431. doi: 10.1038/s41598-017-16647-0 Zhou L. Shen H. Chen L. Li H. Zhang P. Zhao X. . (2019). Species richness and composition of shrub-encroached grasslands in relation to environmental factors in northern China. J. Plant Ecol. 12, 56–66. doi: 10.1093/jpe/rtx062 Zhou G. Zhang X. (1996). Study on climate-vegetation classification for global change in China. Acta Bot. Sin. 38, 8–17. Zhu Y. Wang G. Li R. (2016). Seasonal dynamics of water use strategy of two salix shrubs in alpine sandy land, Tibetan Plateau. PloS One 11, e0156586. doi: 10.1371/journal.pone.0156586