Next Article in Journal
Aspen and Spruce Densities Affect Tree Size, Future Stand Volume, and Aboveground Carbon Following Precommercial Thinning
Previous Article in Journal
Characterization of Fomes fomentarius s.s. and F. inzengae in Belgian Beech Forests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau

1
School of Life Sciences, Ludong University, Yantai 264025, China
2
School of Life Sciences, Northwest University, Xi’an 710069, China
3
China ENFI Engineering Co., Ltd., Beijing 100038, China
4
Institute of Earth Environment, Chinese Academy of Sciences, Xi’an 710058, China
*
Authors to whom correspondence should be addressed.
Forests 2024, 15(2), 222; https://doi.org/10.3390/f15020222
Submission received: 5 December 2023 / Revised: 6 January 2024 / Accepted: 15 January 2024 / Published: 23 January 2024

Abstract

:
Shrubland is an important type of vegetation in the semi-arid region of the western Loess Plateau, and it is of great value to the maintenance of biodiversity and soil and water conservation. Equally, there may be significant differences in the level of diversity and the community assembly patterns across different types of shrublands (desert, alpine, and secondary shrublands). This study aimed to reveal the diversity and the community assembly patterns of different types of shrublands and the environmental explanations for these by using a taxonomic and phylogenetic analysis approach, as well as by considering soil and climate factors. The diversity level of the desert shrublands was low, and the habitats of the three types of shrubland were significantly different. Precipitation may be the main environmental factor driving the variation in and diversity of these types of shrubland. All three types of shrubland were strongly affected by environmental filtering and competitive exclusion, and their community assembly patterns were similar.

1. Introduction

The Loess Plateau experiences the most severe soil erosion of any region in the world, contributing to 97% of the Yellow River’s total sediment content every year [1]. Within it, the western region of the Loess Plateau experiences this more seriously due to its harsh habitat and the serious human interference, low vegetation coverage, and soil erosion present here [2]. Shrubland is the main vegetation type in the western region of the Loess Plateau, which is of great significance to ecological security and soil and water conservation in this region [3]. Therefore, a comparative study of the existing shrubland formation in the western region of the Loess Plateau can help us understand the ecological differences between different shrubland types, clarify the reasons for the formation of the corresponding shrubland communities and provide theoretical support for vegetation restoration and biodiversity protection in this ecologically fragile area [4].
Here, “shrubs” specifically refer to the desert shrubland that is distributed in arid (including extremely arid, arid, and semi-arid) areas with an average annual precipitation of less than 400 mm; the alpine shrubland that is distributed above the upper limit of the vertical distribution of forest; and the secondary shrubland that has formed as part of the community succession process due to high-intensity and large-scale human activities (such as the secondary shrubland communities formed after forest harvesting and the return of farmland to forest). In the western region of the Loess Plateau, these three types of typical desert, alpine, and secondary shrubland communities vary in their formation. Therefore, due to the different reasons for the formation of different types of shrublands, their diversity levels and ecological processes are obviously different. Biodiversity is a good indicator that characterizes the relationship between plants and environmental conditions [5]. It is necessary to explore the formation and maintenance mechanisms of the diversity of the three types of shrubland, clarify the environmental causes of the formation of these communities, and then reveal the development rules of the community in the secondary succession sequence of the Loess Plateau [6]. These studies may judge the community succession outcomes as well. However, there are few relevant studies, and the mechanism of assembling different shrub communities has not yet been discussed.
The phylogeny of a community can be regarded as a “proxy” for the functional traits of species in the absence of data on their functional traits; that is, it is assumed that closely related species in a community have similar functional traits and can be filtered by their environment [7]. Distantly related species can also occupy different ecological niches to avoid the limitations posed by similarity [8]. Therefore, studying diversity and community assembly mechanisms using phylogenetic methods has some advantages, such as high sensitivity and strong ecological significance. Meanwhile, species diversity can directly reflect the taxonomic differences between coexisting species in a community, thus revealing the spatiotemporal variation in species composition [9,10]. Thus, studies on community diversity and its formation and maintenance mechanisms should take into account taxonomic and phylogenetic diversity to clarify the pattern of diversity within a community from wider perspectives [11].
Based on vegetation maps and literature records, this study took the shrublands of the Shaanxi, Gansu, and Ningxia provinces of China in the western Loess Plateau as the research object. The taxonomic diversity, phylogenetic diversity, and community assembly patterns of desert shrublands, alpine shrublands, and secondary shrublands were explored. The diversity of the three shrublands and the mechanisms of maintaining this diversity were combined with geographical, climatic, and soil factors to explore the environmental causes of the formation of different types of communities. Based on the above research, we tried to answer the following questions: (1) Do environmental factors affect the distribution of shrubland communities? Which environmental factors dominate in the formation of shrubland community diversity? (2) Are there differences in the diversity patterns and assembly mechanisms of different types of shrublands?

2. Materials and Methods

2.1. Study Area

The research area of this study is located in the western Loess Plateau, covering three provinces, Gansu, Shaanxi, and Ningxia, in China. There are many different geographical units in the study area, such as the Tengger Desert, the Mu Us Desert, the Helan Mountains, the Liupan Mountains, the Qinling Mountains, the Ziwuling Mountains, the Huanglong Mountains, the loess hilly and gully region, etc., with a complex topography and diverse landforms. The altitude range is 550~3700 m, the mean annual temperature is −5~17 °C, and the mean annual precipitation is 130~1000 mm (Table 1 and Table A1). A total of 52 shrub types were selected for this study, including 24 desert shrub types, 12 alpine shrub types, and 18 secondary shrub types. The sites sampled in this study were selected to investigate the most typical areas of each type. The specific shrub types and study sites are shown in Table 1. In our study, we selected shrubland communities with typical single dominant species for analysis to determine whether the community type was alpine, desert, or secondary shrubland. If the dominant species belonged to typical alpine species (flora recorded as only distributed in alpine environments), the shrubland was defined as alpine shrubland. Where the dominant species were typical desert species (e.g., those belonging to the Amaranthaceae or Zygophyllaceae or those whose flora was clearly recorded as growing in deserts), the shrubland was defined as desert shrubland. The remaining types of shrubland, especially non-zonal vegetation, were defined as secondary shrubland.

2.2. Investigation Methods

From 2018 to 2020, a field survey was conducted during the plant growth period from July to August, a time when the leaves of the plants are already mature and capable of reaching the maximum photosynthetic rate. In order to minimize differences between different shrubland groups being caused by different sampling times, the sampling time was unified, and the sampling was concentrated in July–August of the same year. In 52 plots, three quadrats were arranged at sites with typical vegetation. The distribution area of the shrublands where the quadrats were located was no less than 25 × 25 m. The quadrats were 5 × 5 m, and the radius of the distribution of the three quadrats did not exceed 250 m. The geographical factors, such as latitude and longitude, elevation, slope, and slope direction, of each quadrat were recorded. The identified name, abundance, and coverage of the shrub species in the quadrat were recorded. For further analysis, species from the three quadrats were aggregated by plot. The actual survey area of each plot was 75 m2. A total of 191 shrub species were counted in the survey.

2.3. Environmental Factor Data

Geographical factors were recorded in the field survey. In addition, at least 3 points were selected along the diagonal of each 5 × 5 m quadrat. Then, the undecomposed and semi-decomposed litter layers were removed, and 0–20 cm soil samples were collected for each quadrat, mixed, and put into self-sealing bags. A total of 9 mixed soil samples were taken from each quadrat. These samples were brought back to the laboratory for analysis of the physicochemical properties of the soil. The soil’s pH and organic carbon (TOC), total nitrogen (TN), and total phosphorus (TP) contents were determined using standard methods [12]. The average value was taken as representative data on the soil properties of the whole quadrat.
According to the longitude and latitude information, a 1 × 1 km grid map depicting the climatic factors in the region was taken from the WorldClim database (http://www.worldclim.org/, accessed on 15 December 2022). As supported by previous studies, this study selected mean annual temperature (MAT) and mean annual precipitation (MAP) as the environmental factors for analysis [13]. The above data were obtained using the raster package in R 4.1.2 software (R core team, San Diego, CA, USA).

2.4. Data Analysis

In order to investigate whether similar types of shrubland have consistent preferences in terms of the habitats selected, principal component analysis (PCA) was used to treat the environmental factors of each community to identify the differences in habitat between different types of shrublands. The PCA and mapping were performed in Canoco 4.5 software (Wageningen UR, Wageningen, The Netherlands).
To analyze the diversity patterns of different shrubland communities, species richness and the Shannon–Wiener index were selected as the species diversity indexes in this study [14], and Faith’s PD (namely the evolutionary branch length of all species in the community [15]) was selected as an indicator of the phylogenetic diversity.
Before phylogenetic analysis, it is necessary to construct a suitable phylogenetic supertree to reach accurate conclusions. In this study, to avoid false positive results being caused by an excessively large pool of species, species lists of desert shrubs, alpine shrubs, and secondary shrubs were selected to construct species pools for desert, alpine, and secondary shrubs. Then, phylogenetic trees of the three different types of communities were constructed based on these data. The phylogenetic database of species reported by Zanne [16] was used for the phylogenetic tree construction to assign evolutionary branch lengths. Construction of the phylogenetic trees was carried out in the “phylocomr” package in R 4.1.2 software. The species richness and Shannon–Wiener index were calculated using the “vegan” package in R 4.1.2 software. Faith’s PD was calculated using the “phylo” package in R 4.1.2.
In order to clarify the formation and maintenance mechanisms of diversity, the net relatedness index (NRI) was the factor used to describe the phylogenetic structure. The specific calculation method of the NRI is as follows:
N R I S = 1 × M P D S M P D r n d s s d ( M P D r n d s )
where MPD is the mean lineage distance between species in a quadrat [17]. The corner label S represents the actual observed value, and “rnds” represents the random simulation value of the software. When the NRI > 0, the phylogeny of the community shows aggregation, while at an NRI < 0, the phylogeny shows divergence; if the NRI/NTI = 0, there is random phylogenetic relation in the community [18]. The aggregation and divergence in the phylogeny can explain the key processes that influence community assembly and thus elucidate the mechanisms that maintain the diversity of the community. The NRI calculations were performed in the “COMSTRUCT” module in Phylocom 3.0 software (Webb et al., Cambridge, MA, USA). The diversity and phylogenetic structure of different types of shrubs were compared using ANOVA. The ANOVA was conducted in Statistica 9.0 software (StatSoft Inc., Tulsa, OK, USA). Multiple comparison box diagrams were produced in Origin 8.5 Pro (OriginLab Co., Studio City, CA, USA).
To clarify whether environmental factors influence the diversity and community assembly between different communities, Pearson analysis was used to explore the relationship between different environmental factors and the diversity, phylogenetic structure, and diversity index. Then, MAP was selected as the main influencing factor based on this method. The diversity and phylogenetic structure data were linearly fitted to the MAP gradient. The Pearson correlation analysis was performed in SPSS 19.0 software (IBM, Armonk, NY, USA). The linear fitting and mapping were carried out in Origin 8.5 Pro.

3. Results

3.1. Differences between Environmental Factors among Different Types of Shrublands

In this study, the climate factors and soil factors of different types of shrublands were all involved in the principal component analysis. The results showed that principal component axes 1 and 2 can explain 92.6% of the total variation in the environmental factors (of which the PC1 axis can explain 82.4% and the PC2 axis can explain 10.2%; see Figure 1). In the PCA, the distribution range of the desert, alpine, and secondary shrublands was clear, indicating that the different types of shrublands had obvious environmental differences. The different types of shrublands were separated by the PC2 axis, which was mainly related to MAP. Meanwhile, the alpine shrublands and secondary shrublands were separated by the PC1 axis as well, which was related to the MAT, SOC, TN, etc., indicating that these environmental factors may be the reasons behind the colonization and growth of different types of shrublands.

3.2. Community Assembly Patterns of Different Types of Shrublands

Pearson analysis showed that the species richness was highly correlated with the Shannon–Wiener index (Pearson = 0.859, p ≈ 0). Meanwhile, the species richness and phylogenetic diversity (Pearson = 0.964, p ≈ 0) and the Shannon–Wiener index and phylogenetic diversity (Pearson = 0.838, p ≈ 0) showed significant correlations as well. The results showed that the variation patterns of several diversity indexes were consistent in different types of shrublands. In terms of taxonomic diversity, the species richness and the Shannon–Wiener index of the desert shrubland were the lowest, while there was no significant difference in taxonomic diversity between the alpine shrubland and secondary shrubland species (Figure 2 and Figure 3). The pattern in the phylogenetic diversity was similar to that in the taxonomic diversity. The phylogenetic diversity (PD) of the desert shrublands was significantly lower than that of the alpine shrublands and secondary shrublands (Figure 4). However, a comparison of the phylogenetic structure of the three types of shrub showed that there was no significant difference between the NRIs of the different types of shrublands, and the NRI of most communities was > 0, with no significant variation from 0 (Figure 5).

3.3. Diversity Patterns among Different Types of Shrublands along the Precipitation Gradient

The Pearson analysis showed that there was no significant correlation between the soil factors and the taxonomic diversity, phylogenetic diversity, and phylogenetic structure. Among the climatic factors, the mean annual precipitation (MAP) was related to the taxonomic and phylogenetic diversity but not to the phylogenetic structure. Therefore, this study further discussed the rules of the variation in diversity among different types along the MAP gradient. The results showed that the different types of shrublands had an obvious distribution range along the MAP gradient. The desert shrubland was only distributed where the MAP was < 450 mm, the alpine shrubland experienced a MAP between 470 and 630 mm, and the secondary shrubland relied on a MAP between 500 and 1000 mm. In addition, the taxonomic diversity and phylogenetic diversity increased with an increase in precipitation (Figure 6, Figure 7 and Figure 8).

4. Discussion

4.1. Variation in Diversity among Different Types of Shrublands

In this study, the taxonomic diversity and phylogenetic diversity among different types of shrublands showed a consistent pattern, and the species richness, Shannon–Weiner index, and Faith’s PD were highly correlated. These results are similar to those of previous studies [19], indicating that the selected indexes of taxonomic and phylogenetic diversity are appropriate.
Previous studies have found that with the progression of community succession, the community diversity may first increase and then decrease, and the inflection point appears at the shrub stage [20]. In one particular study, the secondary shrubland emerged due to the secondary succession sequence after being disturbed by humans, and there was not a climax community at the site, as it had a tendency toward forward succession [21]. Therefore, the diversity level of secondary shrublands may be high. In addition, due to the influence of human activities, more species closely related to human activities and even exotic species may be introduced into the secondary shrub, thus further increasing the diversity level of the secondary shrublands. In this study, it was also found that the diversity of the secondary shrublands at the taxonomic and phylogenetic levels was higher than that of the other shrublands, which can also confirm the previous hypothesis.
In this study, the taxonomic and phylogenetic diversity of the desert shrublands was significantly lower than that of the alpine shrublands and secondary shrublands. The habitats of desert shrublands are harsh, with many stress conditions (drought, high temperatures in summer, low temperatures in winter, high soil salinity, and low soil nutrients). Therefore, desert shrubland areas provide an ecological niche that is significantly less readily available than other shrubland areas of the same size. The high level of environmental filtering and interspecific niche competition in desert shrublands results in a significant decrease in the diversity level [22].

4.2. Precipitation Leads to Colonization and Diversity Differences in Different Types of Shrublands

The core argument of community ecology is to explore the relationship between the environment and community assembly [23]. Previous studies have found that the diversity level of a community is related to its local habitat. Factors such as the terrain, moisture and nutrients in the soil, and heat may affect the diversity level of a community [24,25,26,27]. However, in this study, the topographic factors (slope, aspect, elevation, etc.; see Table A2) were not related to the diversity level of the community, which was similar to the results of another study in this region [22], indicating that the topographic factors in the Loess Plateau region were not the main reasons for changes in the diversity level within the community. Although there were significant differences between the mean annual temperature of different places in this study (Table A1), there was no significant change trend in the community diversity pattern along the MAT gradient, and MAT did not seem to be related to the community assembly pattern. This result is different from previous studies [28,29], but other studies have shown that the levels and formation mechanisms of community diversity in alpine regions are not related to environmental factors such as heat [30]. In this study, the habitats of the three types of communities were all stressed by various conditions (such as drought, low temperature, nutrient deficiency, disturbance, etc.), and temperature may no longer be the main limiting factor leading to community changes. Similarly, although the nutrient gradient in the soil may also be the cause of community differences in this study (in the PCA analysis, three different types of communities could be separated by the PC1 and PC2 axes, and PC1 was closely related to soil indexes such as TOC and TN), a single soil factor could not explain the differences in the diversity levels.
In this study, MAP was the only environmental factor associated with changes in diversity levels. Previous studies have shown that differences in water availability can drive differences in diversity [31]. An increase in precipitation can lead to an increase in the woody species diversity in this study area and is a major factor leading to changes in community diversity levels [11]. The results of this study confirm the above conclusions. This indicates that water is still the main factor that limits the formation of shrub communities in the western part of the Loess Plateau.
It should be pointed out that in this study area, different types of shrublands showed obvious differences in the MAP gradient, especially for the desert shrublands, which grow in their entirety in the area with a MAP < 450 mm. Although there was no obvious boundary between the alpine shrubland and secondary shrubland along the MAP gradient, the amplitude of variation in the MAP of the secondary shrubland was higher than that of the alpine shrubland, and the MAP level of certain shrubland sites was higher (>800 mm). This indicates that the alpine shrublands have a clear range of suitable precipitation, and their habitat is more consistent; that is, alpine shrubland communities are formed under specific precipitation conditions (the alpine shrubland community is the climax community). Secondary shrubland is a mid-successional community formed by human activities, and the range of variation in the MAP in its habitat is large; however, the successional outcome of secondary shrubland is inconsistent under different MAP values, and the climax community in the late-successional period may be determined by the type of zonal vegetation in the area. Therefore, a difference in precipitation can determine the type of shrubland formed in the western region of the Loess Plateau and can affect the diversity level of shrubland communities.
Interestingly, for the same type of shrubland, diversity did not show a change trend with MAP. However, this trend appeared when the three shrubland types were brought together (Figure 6, Figure 7 and Figure 8). In fact, this phenomenon has been mentioned in a study of ecological processes at the local and regional scales [25]. In this study, because different types of shrublands have relatively similar combinations of environmental factors, their geographical distribution areas may also be relatively concentrated. We can assume that different types of shrublands are distributed in different areas. Due to differences in the ecological processes at the local to regional scales, different results will appear when considering the diversity of all types of shrubland across MAP gradients at the regional scale.

4.3. The Assembly Process of Different Types of Shrubland Communities Is Similar

Exploring and revealing the spatial patterns and maintenance mechanisms of biodiversity is a hot topic in ecology and biogeography [32]. In this study, the phylogenetic structures of the different types of shrublands were similar, indicating that shrubland communities whose formation had different causes in the western region of the Loess Plateau may have been affected by similar assembly processes. In this region, due to the harshness of the habitat, some species will have been excluded by environmental filtering during the assembly process [33]. The selected species may have similar traits, and the interspecific competition may be more intense. However, due to the limitations of the natural conditions in the region, the available niche within the community is limited, which may further exacerbate the degree of interspecific competition. In this study, the three types of communities all showed a weak pattern of phylogenetic convergence, but there was no significant difference from the null model, indicating that the shrubland communities may have been superimposed by environmental filtering and similarity limitations, which masked the independent influence of two opposite ecological processes on community formation and maintenance [34].
It should be noted in particular that in this study, species pools of different sizes were constructed when the phylogenetic structure was analyzed for the three different types of shrublands. The species pools of the three types of shrubland were the sum total shrub species in the community of that type. In previous studies, selecting species pools of different sizes had a great impact on the phylogenetic structure of a region [25]. In general, a species pool should be constructed according to all species and potential species in the study area. However, the three types of shrubland in this study have significant environmental differences (the PCA analysis showed that the three types of shrubland have obvious environmental differences), and the selection of a larger regional species pool will inevitably exaggerate the role of many species that are unlikely to appear in certain types of shrubland, resulting in a convergence pattern when calculating the phylogenetic structure that is inconsistent with reality. Therefore, in this study, only species that had appeared in the specific type of shrubland were selected for constructing the phylogenetic supertree. Relatively speaking, the selection of smaller supertrees may reduce the detection sensitivity of the patterns of phylogenetic convergence, further resulting in most of the communities in this study exhibiting insignificant convergence or even random community assembly patterns.

5. Conclusions

This study combined the diversity levels and phylogenetic structures of desert, alpine, and secondary shrublands in the western region of the Loess Plateau with the environmental factors affecting these communities to explore the effect of the former on the latter. It was found that the environmental factors in different types of shrubland communities were significantly different, and the diversity level of the desert shrublands was significantly lower than that of the alpine and secondary shrublands. The precipitation gradient can distinguish different types of shrublands and seemingly influences their diversity levels. In the western part of the Loess Plateau, the community assembly mechanism of different types of shrublands is similar, and they are all affected by higher levels of environmental filtering and competitive exclusion processes.

Author Contributions

Conceptualization, X.L.; writing—original draft preparation, J.X.; formal analysis, H.D.; supervision, D.H. and P.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the National Natural Science Foundation of China (41871036); the National Basic Science and Technology Projects of China (2015FY110300); the General Survey of Herbaceous Plant Germplasm Resources in Shandong Province (LCYZ [2021] No. 1); the Key R&D Plan of Shandong Province (Major Scientific and Technological Innovation Project, 2021LZGC023); the Shandong Provincial Agricultural Elite Varieties Project (2019LZGC01805); and the Open Research Fund of Guizhou Provincial Key Laboratory for Biodiversity Conservation and Utilization in the Fanjing Mountain Region ([2020]2003).

Data Availability Statement

The original contributions presented in this study are included in the article; further inquiries can be directed to the corresponding authors.

Acknowledgments

The authors wish to thank Yongfu Chai, Mao Wang, Pengcheng Wan, Jiangchao Guo, Guoxun Wang, Yu Chen, Maolin Lei, Lei Wang, Peng Zhao, and Chengcheng Xiang for their assistance in the field sampling. We also thank Ming Yue and Peiliang Liu for their guidance on the manuscript, especially on species identification.

Conflicts of Interest

Han Dang was employed by the China ENFI Engineering Co., Ltd. The remaining 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. The China ENFI Engineering Co., Ltd. had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Climate and soil environment information of the sampling sites.
Table A1. Climate and soil environment information of the sampling sites.
Shrubland TypeMean
Annual
Temperature (MAT)/°C
Mean Annual Precipitation (MAP)/mmSoil Organic Carbon (SOC)/mg·g−1Soil Total
Nitrogen (TN)/mg·g−1
Soil Total
Phosphorus (TP)/mg·g−1
pH of Soil
DNRS8.3293.06.60450.58650.38067.7956
DNAC8.3293.03.33250.27180.36277.6611
DNCJ7.4306.08.10560.77670.65868.3144
DNSP7.4262.09.25820.70910.53268.3078
DNOA7.4262.010.58790.85660.51928.0156
DNCT7.4262.013.29801.32530.56538.1833
DNAO8.6260.02.02430.15700.23078.8378
DNKG8.6275.03.96580.32740.41608.5856
DNNT8.6275.04.45280.33970.29268.7189
DNAM6.7228.021.77941.92200.57477.8400
DNLO6.7228.027.91032.51420.62407.8356
DNEL6.8226.024.50272.13160.52237.9311
DNCK8.6142.00.83310.08540.12398.7622
DNMM8.6168.014.73070.35280.19888.2056
DGAB8.4201.02.65180.21040.14809.0356
DGSX8.3202.01.79350.06550.13719.2689
DGNS8.4200.01.13540.03080.21219.1067
DGCM8.1127.01.14070.02870.19588.2256
DGEP8.0118.02.69560.12300.14249.0167
DGNP8.1119.01.70130.08030.27047.9667
DGNR8.2116.04.21690.23770.99468.2111
DGTC8.3146.01.62100.07200.21169.5744
DGHA8.3146.0119.140511.42170.82795.9922
DSSC7.2436.00.38000.01860.18959.3847
AGRC1.0513.080.01247.73461.10165.4000
AGPF1.3520.080.96456.24580.69806.3489
AGRP0.4504.0130.204314.37741.19395.5011
AGSG0.0495.018.70895.16081.30716.8867
AGBT0.0495.051.668911.42481.09946.9111
AGSA0.0495.062.26864.56670.74577.3756
AGSC4.4681.054.56815.94401.19256.0133
AGSS4.4681.0106.396813.07981.45157.2644
AGHR4.4681.070.70298.22891.16645.8967
AGPP2.6605.079.15627.07740.91846.4967
AGBK1.9603.048.47804.76311.10457.2233
AGSQ1.9603.084.77388.35911.04117.2789
SSOD10.9600.01.50000.21000.57249.0615
SSVN10.2608.00.63120.16550.41068.7353
SSRH12.9550.01.17930.06550.51819.1766
SSSV13.8540.00.99380.08600.58569.0032
SSSO10.6564.00.85520.15960.57089.2357
SSZJ10.7616.00.78670.09410.62359.3220
SSRP14.2606.00.84660.10590.24077.1597
SSGB14.2606.00.80460.17040.60326.2511
SSLF14.2606.00.99860.23590.52275.7449
SSCG13.7658.01.83280.17240.73469.0609
SSCN13.1748.00.24720.04500.32879.5806
SSMA16.0719.02.15330.25510.58779.4934
SSFS16.0719.02.87700.29750.47679.1644
SSSS14.8763.02.35910.30610.78698.6214
SSLG13.9937.00.91970.32860.60796.2379
SSPP13.9937.00.85080.57630.75627.2977
SSCC10.3729.018.29051.55750.38028.3170
SSQB14.3832.019.19170.93010.15734.8156
Table A2. Pearon correlation between environmental factors and diversity indicators.
Table A2. Pearon correlation between environmental factors and diversity indicators.
FactorsRichnessPD
Elevation−0.125−0.099
Slope0.2280.281
mean annual temperature (MAT)0.2680.254
mean annual precipitation (MAP)0.475 **0.530 **
organic carbon (TOC)−0.0320.01
total nitrogen (TN)−0.0350.006
total phosphorus (TP)0.1290.171
pH0.022−0.062
** means significant correlation.

References

  1. Liang, H.; Xue, Y.; Li, Z.; Wang, S.; Wu, X.; Gao, G.; Liu, G.; Fu, B. Soil moisture decline following the plantation of Robinia pseudoacacia forests: Evidence from the Loess Plateau. For. Ecol. Manag. 2018, 412, 62–69. [Google Scholar] [CrossRef]
  2. Fu, B.; Wang, S.; Liu, Y.; Liu, J.; Liang, W.; Miao, C. Hydrogeomorphic Ecosystem Responses to Natural and Anthropogenic Changes in the Loess Plateau of China. Annu. Rev. Earth Planet. Sci. 2017, 45, 223–243. [Google Scholar] [CrossRef]
  3. Jia, X.; Shao, M.; Zhu, Y.; Luo, Y. Soil moisture decline due to afforestation across the Loess Plateau, China. J. Hydrol. 2017, 546, 113–122. [Google Scholar] [CrossRef]
  4. Chu, C.; Lutz, J.A.; Král, K.; Vrška, T.; Yin, X.; Myers, J.A.; Abiem, I.; Alonso, A.; Bourg, N.; Burslem, D.F.; et al. Direct and indirect effects of climate on richness drive the latitudinal diversity gradient in forest trees. Ecol. Lett. 2018, 22, 245–255. [Google Scholar] [CrossRef]
  5. Liu, J.; Yunhong, T.; Slik, J.F. Topography related habitat associations of tree species traits, composition and diversity in a Chinese tropical forest. For. Ecol. Manag. 2014, 330, 75–81. [Google Scholar] [CrossRef]
  6. Fang, S.; Munoz, F.; Ye, J.; Lin, F.; Yuan, Z.; Kuang, X.; Hao, Z.; Wang, X. Deterministic processes drive functional and phylogenetic temporal changes of woody species in temperate forests in Northeast China. Ann. For. Sci. 2019, 76, 42. [Google Scholar] [CrossRef]
  7. Park, D.S.; Potter, D. Why close relatives make bad neighbours: Phylogenetic conservatism in niche preferences and dispersal disproves Darwin’s naturalization hypothesis in the thistle tribe. Mol. Ecol. 2015, 24, 3181–3193. [Google Scholar] [CrossRef] [PubMed]
  8. Wiens, J.J.; Graham, H. Niche conservatism: Integrating evolution, ecology, and conservation biology. Ann. Rev. Ecol. Evol. Syst. 2005, 36, 519–539. [Google Scholar] [CrossRef]
  9. Swenson, N.G.; Erickson, D.L.; Mi, X.C.; Bourg, N.A.; Forero-Montaña, J.; Ge, X.; Howe, R.; Lake, J.K.; Liu, X.; Ma, K.; et al. Phylogenetic and functional alpha and beta diversity in temperate and tropical tree communities. Ecology 2012, 93, S112–S125. [Google Scholar] [CrossRef]
  10. McGill, B.J.; Etienne, R.S.; Gray, J.S.; Alonso, D.; Anderson, M.J.; Benecha, H.K.; Dornelas, M.; Enquist, B.J.; Green, J.L.; He, F. Species abundance distributions: Moving beyond single prediction theories to integration within an ecological framework. Ecol. Lett. 2007, 10, 995–1015. [Google Scholar] [CrossRef]
  11. Xu, J.; Dang, H.; Wang, M.; Chai, Y.; Guo, Y.; Chen, Y.; Zhang, C.; Yue, M. Is Phylogeny More Useful than Functional Traits for Assessing Diversity Patterns under Community Assembly Processes? Forests 2019, 10, 1159. [Google Scholar] [CrossRef]
  12. Cornelissen, J.H.C.; Lavorel, S.; Garnier, E.; Díaz, S.; Buchmann, N.; Gurvich, D.E.; Reich, P.B.; ter Steege, H.; Morgan, H.D.; van der Heijden, M.G.A.; et al. A handbook of protocols for standardised and easy measurement of plant functional traits worldwide. Aust. J. Bot. 2003, 51, 335–380. [Google Scholar] [CrossRef]
  13. Xu, J.; Zhang, Y.; Liu, S.; Chai, Y.; Dang, H.; Yue, M.; Liu, X.; Guo, Y. Patterns of diversity and community assembly change across local to regional scales: An evidence of deterministic assembly processes along resource availability gradient at temperate forest. Ecol. Indic. 2021, 132, 108261. [Google Scholar] [CrossRef]
  14. Gentry, A.H. Patterns of Neotropical Plant Species Diversity. In Evolutionary Biology; Springer: Boston, MA, USA, 1982; pp. 1–84. [Google Scholar]
  15. Faith, D.P. Conservation evaluation and phylogenetic diversity. Biolog. Conser. 1992, 61, 1–10. [Google Scholar] [CrossRef]
  16. Zanne, A.E.; Tank, D.C.; Cornwell, W.K.; Eastman, J.M.; Smith, S.A.; FitzJohn, R.G.; McGlinn, D.J.; O’meara, B.C.; Moles, A.T.; Reich, P.B.; et al. Three keys to the radiation of angiosperms into freezing environments. Nature 2013, 506, 89–92. [Google Scholar] [CrossRef] [PubMed]
  17. Kraft, N.J.; Cornwell, W.K.; Webb, C.O.; Ackerly, D.D.; Heard, S.B.; Whitlock, M.C. Trait evolution, community assembly, and the phylogenetic structure of ecological communities. Am. Nat. 2007, 170, 271–283. [Google Scholar] [CrossRef] [PubMed]
  18. Webb, C.O.; Ackerly, D.D.; McPeek, M.A.; Donoghue, M.J. Phylogenies and Community Ecology. Annu. Rev. Ecol. Syst. 2002, 33, 475–505. [Google Scholar] [CrossRef]
  19. Safi, K.; Cianciaruso, M.V.; Loyola, R.D.; Brito, D.; Armour-Marshall, K.; Diniz-Filho, J.A.F. Understanding global patterns of mammalian functional and phylogenetic diversity. Phil. Trans. R. Soc. B Biol. Sci. 2011, 366, 2536–2544. [Google Scholar] [CrossRef]
  20. Brym, Z.T.; Lake, J.K.; Allen, D.; Ostling, A. Plant functional traitssuggest novel ecological strategy for an invasive shrub in an understorey woody plant community. J. Appl. Ecol. 2011, 48, 1098–1106. [Google Scholar] [CrossRef]
  21. Xu, J.S.; Dang, H.; Tian, T.; Liu, S.; Chai, Y.; Liu, X.; Yue, M.; Xiang, C.; Chang, J. Human disturbance rather than habitat factors drives plant community assembly and diversity patterns in a semi-arid region. Land Degrad. Dev. 2020, 31, 1803–1811. [Google Scholar] [CrossRef]
  22. Laliberté, E.; Zemunik, G.; Turner, B.L. Environmental filtering explains variation in plant diversity along resource gradients. Science 2014, 345, 1602–1605. [Google Scholar] [CrossRef] [PubMed]
  23. Chase, J.M. Community assembly: When should history matter? Oecologia 2003, 136, 489–498. [Google Scholar] [CrossRef]
  24. Punchi-Manage, R.; Wiegand, T.; Wiegand, K.; Getzin, S.; Gunatilleke, C.V.S.; Gunatilleke, I.A.U.N. Effect of spatial processes and topography on structuring species assemblages in a Sri Lankan dipterocarp forest. Ecology 2014, 95, 376–386. [Google Scholar] [CrossRef] [PubMed]
  25. Xu, J.; Chen, Y.; Zhang, L.; Chai, Y.; Wang, M.; Guo, Y.; Li, T.; Yue, M. Using phylogeny and functional traits for assessing community assembly along environmental gradients: A deterministic process driven by elevation. Ecol. Evol. 2017, 7, 5056–5069. [Google Scholar] [CrossRef] [PubMed]
  26. Shiono, T.; Kusumoto, B.; Maeshiro, R.; Fujii, S.J.; Götzenberger, L.; de Bello, F.; Kubota, Y. Climatic drivers of trait assembly in woody plants in Japan. J. Biogeogr. 2015, 42, 1176–1186. [Google Scholar] [CrossRef]
  27. Dang, H.; Li, J.; Xu, J.; Chu, G.; Zhang, J.; Yu, Y.; Jin, Z. Differences in soil water and nutrients under catchment afforestation and natural restoration shape herbaceous communities on the Chinese Loess Plateau. For. Ecol. Manag. 2022, 505, 119925. [Google Scholar] [CrossRef]
  28. Hulshof, C.M.; Violle, C.; Spasojevic, M.J.; McGill, B.; Damschen, E.; Harrison, S.; Enquist, B.J. Intra-specific and inter-specific variation in specific leaf area reveal the importance of abiotic and biotic drivers of species diversity across elevation and latitude. J. Veg. Sci. 2013, 24, 921–931. [Google Scholar] [CrossRef]
  29. Qian, H.; Hao, Z.Q.; Zhang, J. Phylogenetic structure and phylogenetic diversity of angiosperm assemblages in forests along an elevational gradient in Changbaishan, China. J. Plant Ecol. 2014, 7, 154–165. [Google Scholar] [CrossRef]
  30. Zhang, H.; Gilbert, B.; Zhang, X.; Zhou, S. Community assembly along a successional gradient in sub-alpine meadows of the Qinghai-Tibetan Plateau, China. Oikos 2013, 122, 952–960. [Google Scholar] [CrossRef]
  31. Eiserhardt, W.L.; Svenning, J.-C.; Kissling, W.D.; Balslev, H. Geographical ecology of the palms (Arecaceae): Determinants of diversity and distributions across spatial scales. Ann. Bot. 2011, 108, 1391–1416. [Google Scholar] [CrossRef]
  32. Cavender-Bares, J.; Kozak, K.H.; Fine, P.V.A.; Kembel, S.W. The merging of community ecology and phylogenetic biology. Ecol. Lett. 2009, 12, 693–715. [Google Scholar] [CrossRef] [PubMed]
  33. Gianuca, A.T.; Dias, R.A.; Debastiani, V.J.; Duarte, L.D.S. Habitat filtering influences the phylogenetic structure of avian communities across a coastal gradient in southern Brazil. Austral Ecol. 2013, 39, 29–38. [Google Scholar] [CrossRef]
  34. Mayfield, M.M.; Levine, J.M. Opposing effects of competitive exclusion on the phylogenetic structure of communities. Ecol. Lett. 2010, 13, 1085–1093. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The PCA analysis of different types of shrublands. Red circles refer to alpine shrublands, green circles refer to desert shrublands, and blue circles refer to secondary shrublands. The ellipses represent aggregations of samples.
Figure 1. The PCA analysis of different types of shrublands. Red circles refer to alpine shrublands, green circles refer to desert shrublands, and blue circles refer to secondary shrublands. The ellipses represent aggregations of samples.
Forests 15 00222 g001
Figure 2. The richness of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Figure 2. The richness of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Forests 15 00222 g002
Figure 3. The Shannon–Wiener index of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Figure 3. The Shannon–Wiener index of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Forests 15 00222 g003
Figure 4. The phylogenetic diversity (PD) of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Figure 4. The phylogenetic diversity (PD) of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Forests 15 00222 g004
Figure 5. The net relatedness index (NRI) of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Figure 5. The net relatedness index (NRI) of different types of shrublands. Different letters indicate significant differences. A refers to alpine shrublands, D refers to desert shrublands, and S refers to secondary shrublands.
Forests 15 00222 g005
Figure 6. The species richness of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Figure 6. The species richness of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Forests 15 00222 g006
Figure 7. The Shannon–Wiener index of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Figure 7. The Shannon–Wiener index of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Forests 15 00222 g007
Figure 8. The phylogenetic diversity (PD) of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Figure 8. The phylogenetic diversity (PD) of shrubland community along mean annual precipitation (MAP) gradient across different types of shrublands. The green circles represent desert shrublands, the red circles represent alpine shrublands, and the blue circles represent secondary shrublands, while the black line represents line fitting.
Forests 15 00222 g008
Table 1. Site information.
Table 1. Site information.
Community TypeSite CodeDominant SpeciesLocationLatitude/°NLongitude/°EElevation/m
Desert
shrubland
DNRSReaumuria songarica (Pall.) Maxim.Tongxin, Ningxia36.915375105.8105811544
DNACAsterothamnus centraliasiaticus Novopokr.Tongxin, Ningxia36.915222105.8105891530
DNCJCaragana jubata (Pall.) Poir.Wuzhong, Ningxia36.719008105.5406221779
DNSPCaroxylon passerinum (Bunge) Akhani and RoalsonZhongning, Ningxia36.999789105.2574141761
DNOAOxytropis aciphylla Ledeb.Zhongning, Ningxia36.999381105.2581141780
DNCTCaragana tibetica Kom.Zhongning, Ningxia36.999097105.2579421784
DNAOArtemisia ordosica Krasch.Zhongwei, Ningxia37.448169106.1733311316
DNKGKalidium gracile FenzlTongxin, Ningxia37.437297106.4498081313
DNNTNitraria tangutorum Bobr.Tongxin, Ningxia37.437444106.4480361320
DNAMPrunus mongolica (Maxim.) RickerHelan, Ningxia38.609011105.9270331591
DNLOLeptodermis ordosica H. C. Fu et E. W. MaHelan, Ningxia38.609169105.9272081621
DNELEphedra rhytidosperma PachomovaHelan, Ningxia38.606356105.9280501655
DNCKCaragana korshinskii Kom.Zhongning, Ningxia39.348111105.6349031290
DNMMAmmopiptanthus mongolicus (Maxim. ex Kom.) Cheng f.Shizuishan, Ningxia39.405397106.7298781097
DGABAnabasis brevifolia C. A. Mey.Jingtai, Gansu37.363667104.1379001592
DGSXZygophyllum xanthoxylum (Bunge) Maxim.Jingtai, Gansu37.374603104.1632691614
DGNSNitraria sibirica Pall.Jingtai, Gansu37.375506104.1598001613
DGCMCalligonum mongolicum Turcz.Wuwei, Gansu38.443858102.9322361388
DGEPEphedra przewalskii StapfMinqin, Gansu38.566175102.8452191373
DGNPNitraria sphaerocarpa Maxim.Yongchang, Gansu38.729933102.3347811372
DGNRNitraria roborowskii Kom.Yongchang, Gansu38.758275102.3506001348
DGTCTamarix chinensis Lour.Minqin, Gansu38.203958102.7651561443
DGHAHaloxylon ammodendron (C. A. Mey.) BungeMinqin, Gansu38.204842102.7663001449
DSSCSalix cheilophila Schneid.Shenmu, Shaanxi38.787222110.2430561223
Alpine
shrubland
AGRCRhododendron capitatum MaximTianzhu, Gansu36.959111102.8161693326
AGPFDasiphora fruticosa (L.) Rydb.Tianzhu, Gansu36.958769102.8168643313
AGRPRhododendron przewalskii Maxim.Tianzhu, Gansu36.955908102.8289063320
AGSGSalix gilashanica C. Wang et P. Y. FuTianzhu, Gansu36.944772102.8312283244
AGBTBerberis thunbergii DC.Tianzhu, Gansu36.945528102.8319583251
AGSASpiraea alpina Pall.Tianzhu, Gansu36.945497102.8319973249
AGSCSalix cupularis Rehd.Gannan, Gansu33.841131104.2642502739
AGSSSorbaria sorbifolia (L.) A. Br.Gannan, Gansu33.840844104.2628692809
AGHRHippophae rhamnoides L.Gannan, Gansu33.841431104.2646332740
AGPPDasiphora parvifolia (Fisch. ex Lehm.) Juz.Lintan, Gansu34.757606103.6095283187
AGBKBerberis kansuensis Schneid.Kangle, Gansu34.939761103.7522563212
AGSQJuniperus squamata Buchanan-Hamilton ex D. DonKangle, Gansu34.939764103.7522583300
Secondary shrublandSSODOstryopsis davidiana DecaisneYaozhou, Shaanxi35.152300108.8369441150
SSVNVitex negundo var. heterophylla (Franch.) Rehd.Xunyi, Shaanxi35.156047108.2811111138
SSRHRosa hugonis Hemsl.Chengcheng, Shaanxi35.185314109.8686111160
SSSVSophora davidii (Franch.) SkeelsChengcheng, Shaanxi35.070753109.8697221000
SSSOSyringa oblata Lindl.Huanglong, Shaanxi35.583330109.8166701212
SSZJZiziphus jujuba var. spinosa (Bunge) Hu ex H.F.Chow.Fuping, Shaanxi35.062778109.278611947
SSRPRubus piluliferus FockeChang’an, Shaanxi34.042500108.836667566
SSGBGrewia biloba G. DonChang’an, Shaanxi34.042222108.836667650
SSLFLespedeza thunbergii subsp. formosa (Vogel) H. OhashiChang’an, Shaanxi34.042222108.836944687
SSCGCelastrus gemmatus Loes.Zhouzhi, Shaanxi34.058333108.281111606
SSCNCoriaria nepalensis Wall.Shanyang, Shaanxi33.520278109.869722739
SSMAMyrsine africana L.Baihe, Shaanxi32.859167110.010833249
SSFSForsythia suspensa (Thunb.) VahlBaihe, Shaanxi32.859167110.011111295
SSSSSpiraea sericea Turcz.Baihe, Shaanxi32.846944110.009444344
SSLGLindera glauca (Siebold and Zucc.) BlumeNanzheng, Shaanxi32.591667107.157222891
SSPPPittosporum podocarpum Gagnep.Nanzheng, Shaanxi32.593611107.155000858
SSCCCotinus coggygria Scop.Fengxian, Shaanxi34.169419106.7737331250
SSQBQuercus baronii SkanChenggu, Shaanxi32.855461107.213786682
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, J.; Dang, H.; Hu, D.; Zhang, P.; Liu, X. Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau. Forests 2024, 15, 222. https://doi.org/10.3390/f15020222

AMA Style

Xu J, Dang H, Hu D, Zhang P, Liu X. Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau. Forests. 2024; 15(2):222. https://doi.org/10.3390/f15020222

Chicago/Turabian Style

Xu, Jinshi, Han Dang, Dechang Hu, Ping Zhang, and Xiao Liu. 2024. "Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau" Forests 15, no. 2: 222. https://doi.org/10.3390/f15020222

APA Style

Xu, J., Dang, H., Hu, D., Zhang, P., & Liu, X. (2024). Patterns of Diversity and Community Assembly and Their Environmental Explanation across Different Types of Shrublands in the Western Loess Plateau. Forests, 15(2), 222. https://doi.org/10.3390/f15020222

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop