Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 07 August 2024
Sec. Functional Plant Ecology
This article is part of the Research Topic Plant Diversity: The Key to Ecosystem Resilience in a Changing World View all 10 articles

Shallow tillage mitigates plant competition by increasing diversity and altering plant community assembly process

Zihao Li,&#x;Zihao Li1,2†Jirong Qiao&#x;Jirong Qiao3†Zhuofan LiZhuofan Li1Xiaowei GaoXiaowei Gao1Guangyu HongGuangyu Hong1Haifeng YangHaifeng Yang1Ercha Hu,Ercha Hu1,2Chenming LiuChenming Liu1Xiaojiang WangXiaojiang Wang1Guanzhi Liu*Guanzhi Liu3*Runhong Gao*Runhong Gao2*
  • 1Institute of Forest Ecology, Inner Mongolia Academy of Forestry Sciences, Hohhot, China
  • 2College of Forestry, Inner Mongolia Agricultural University, Hohhot, China
  • 3College of Grassland, Resources and Environment, Inner Mongolia Agricultural University, Hohhot, China

Introduction: Understanding how human activities affect biodiversity is needed to inform systemic policies and targets for achieving sustainable development goals. Shallow tillage to remove Artemisia ordosica is commonly conducted in the Mu Us Desert. However, the impacts of shallow tillage on plant community species diversity, phylogenetic structure, and community assembly processes remain poorly understood.

Methods: This study explores the effects of shallow tillage on species diversity including three a-diversity and two b-diversity indicators, as well as phylogenetic structure [phylogenetic diversity (PD), net relatedness index (NRI), and nearest taxon index (NTI)]. Additionally, this research analyzes the effects of shallow tillage on the community assembly process.

Results and discussion: The results showed that the a-diversity index, b-diversity index, and PD of the shallow tillage (ST) communities were significantly higher than those of the non-shallow tillage (NT) communities, and the phylogenetic structures of both the ST and NT communities tended to be differentiated, with competitive exclusion being the main mechanism of plant assembly. However, shallow tillage increased the relative importance of the stochastic processes dominated by dispersal limitation, mitigating plant competition in the communities. This conclusion was supported by the Raup–Crick difference index-based analysis.

Conclusion: Therefore, for the ecological restoration of the Mu Us Desert, species with adaptability and low niche overlap should be selected to increase the utilization efficiency of the environmental resources. The results of this study provide a foundation for policy development for ecosystem management and restoration in the Mu Us Desert.

1 Introduction

Biodiversity underpins ecosystem stability and social development, yet it is declining globally due to human activities (Zhang and Chen, 2022; Zhang et al., 2023). Desert ecosystems cover 12% of the Earth’s land surface (Durant et al., 2012; Chen et al., 2023) and, compared to other ecosystems, are environmentally fragile and sensitive to disturbance, which can have disproportionately severe impacts on biodiversity. Previous studies have shown that even the smallest soil disturbance can significantly alter plant community structure and biodiversity in arid ecosystems (Valone, 2003; Erfanian et al., 2019). Shallow, a common disturbance measure in agroecosystems, is also used for desert grassland improvement, whereby local people remove vegetation that does not provide economic benefits by plowing it to utilize desert grassland resources. This practice, although tailored to local needs, significantly risks community species composition biodiversity and the natural processes supporting desert ecosystems (Randriamalala et al., 2012; Cordeau et al., 2017). Consequently, understanding the ecological consequences of anthropogenic interventions is essential for developing more sustainable management practices (Jaureguiberry et al., 2022).

Community assembly is an intrinsic mechanism of biodiversity change and maintenance and can be used to assess the impacts of external drivers, such as human activities, on biodiversity (Chase, 2010). Ecological niche theory (based on deterministic processes) and neutral theory (based on stochastic processes) are the two hypotheses for the community assembly process, and it is widely accepted that deterministic and stochastic factors work together in the community assembly process. However, assessing the relative importance of these factors is challenging, especially in arid regions. Traditionally, species α-diversity and β-diversity are often used in community assembly process analyses (Mori et al., 2018; Jin et al., 2023). However, species diversity reflects only one dimension of biodiversity and does not consider evolutionary diversity (Capmourteres and Anand, 2016; Heino and Tolonen, 2017). Therefore, the use of phylogenetic structures to analyze community assembly processes has gained prominence (Cadotte et al., 2009; Mishler et al., 2014). This methodology not only assesses the evolution of biodiversity and evolutionary distances among taxa but also reflects functional diversity, as closely related species often exhibit similar traits (Cadotte et al., 2008). Although species diversity and phylogenetic structure are often highly correlated, they may indicate different patterns of community assembly (Tucker and Cadotte, 2013). Accordingly, it is vital to integrate various measures of biodiversity to achieve a more comprehensive understanding of community assembly. Uncovering these hidden patterns can provide more targeted guidance for restoring damaged ecosystems and deepen the understanding of diversity formation and maintenance.

The studies on the effects of tillage on community diversity and assembly processes have mainly focused on agroecosystems (Travlos et al., 2018). It has been shown that the removal of native vegetation by tillage creates a large number of gap ecological niches that favor plant species colonization (Czerwiński et al., 2018). Consequently, tillage alters the species composition of weed communities, in which perennial species shifted to annual species, and increases the community evenness index (Booth and Swanton, 2002) but does not change diversity indices such as species richness and the Shannon–Weiner index (Légère et al., 2005; Randriamalala et al., 2012; Barroso et al., 2015). Tillage imposed environmental filtering, an artificial deterministic process, on the weed community and therefore altered the drivers of weed community assembly (Alarcón et al., 2018). However, the relative importance of deterministic and stochastic processes is largely influenced by community type (Cordeau et al., 2017). Although phylogenetic structure is important for understanding the mechanisms of community diversity change and formation, little is known about the effects of tillage on plant phylogenetic structure.

Artemisia ordosica plays a crucial role in wind and sand control (Wei et al., 2016). However, its strong odor leads to its unpopularity among livestock. Consequently, local herders frequently remove it through shallow tillage, which is defined as soil disturbance up to a depth of 10–15 cm (Cooper et al., 2016; Laudicina et al., 2017) to improve grasslands. This practice promotes the growth of more palatable herbaceous plants to support production and livelihoods in the Mu Us Desert but also heightens the risk of desertification (Wang et al., 2017). However, less research has been conducted on shallow tillage on plant diversity and community assembly processes in desert grasslands. This study evaluates the short-term impacts of shallow tillage on plant communities in the Mu Us Desert, focusing on changes in species diversity and phylogenetic structure, as well as community assembly processes, before and after tillage. Specifically, we aim to answer the following: (1) “Does shallow tillage affect the species diversity and phylogenetic structure of the communities and how?” and (2) “Does shallow tillage affect the ecological processes that drive community assembly?” We propose the hypothesis that shallow tillage increases community species and phylogenetic diversity but does not alter its community assembly processes because studies have shown that shallow tillage leads to a large number of blank ecological niches, while multiple tillage alters weed community assembly processes. The study of the effects of shallow tillage on multilevel of biodiversity (species diversity and phylogenetic diversity) of desert ecosystems is an important reference for conservation strategies and land management practices.

2 Materials and methods

2.1 Study area

The study was conducted in a fixed monitoring site at the Ulan Tolgoi Desertification Control Station located in Wushen Banner (38°48′N, 109°19′E, at 1,390 m), Inner Mongolia, China (Figure 1). This region in the hinterland of the Mu Us desert has a semi-arid continental monsoon climate, which is dry and windy. The average annual temperature ranges from 6.0°C to 9.0°C (Zheng et al., 2019a). The annual precipitation in the area varies greatly, ranging from 250 to 400 mm, and the precipitation is mainly concentrated in July–September, accounting for up to 75% of the annual precipitation. The precipitation in August can account for up to 54% of the whole rainy season, and the average annual evaporation ranges from 2,100 to 2,600 mm (She et al., 2017). The predominant soil type is sandy (FAO soil classification), and the geomorphological types include fixed, semi-fixed, wandering dune, and inter-dune lowlands (Chen et al., 2022). Drought-tolerant plants, such as A. ordosica, Salix cheilophila, Caragana microphylla, and Corethrodendron fruticosum var. mongolicum, are widely distributed in the area (Zhang et al., 2018).

Figure 1
www.frontiersin.org

Figure 1 Geographic location of the study area.

2.2 Sample plots and data collection

In March 2022, shallow tillage was applied to the A. ordosica community (there is no history of any tillage). In mid-August of the same year, 10 shallow tillage (ST) strips and 10 non-shallow tillage (NT) strips, each approximately 8 m wide and 200 m long, were selected, in which the two were interleaved, and the distance between adjacent strips under the same treatment was greater than 200 m to include more heterogeneous habitats. Three 8 m × 20 m sample plots were randomly set up in each shallow-tillage and non-shallow-tillage strip, and a detailed survey of plant species and their abundance occurring in the plots was conducted (Supplementary Table S1) to assess species diversity and phylogenetic structure.

2.3 Species diversity and phylogenetic structure

2.3.1 α-Diversity

Based on the species matrix data that were obtained from the survey, we calculated the Shannon–Weiner index (H), Simpson’s index (D), Pielou’s evenness index (J), and the cumulative number of species as measures of species α-diversity, and the formulas are as follows:

Shannon diversity index:

H=i=1SPilnPi

Simpson’s index:

D=1i=1SPi  2

Pielou’s evenness index:

J=HlnS

where S is the number of species and Pi is the relative abundance of species i. Pi = ni/N, where ni represents the number of individuals of species i and N represents the total number of individuals of all the plant species in a particular sample.

2.3.2 β-diversity

The Bray–Curtis dissimilarity index (βBC) and the Jaccard dissimilarity index (βJ) were used to explore the β-diversity characteristics of the plant communities in the shallow and non-shallow tillage belts. The βBC is based on the species multiplicity matrix to measure the compositional dissimilarity between the plots based on differences in species abundance, and it ranges from 0 (when both plots have the same species with equal abundances) to 1 (when the plots do not share any species; Catalán et al., 2023). The βBC was calculated using the following equation:

βBC=i=1S|xijxik|i=1Sxij+xik

where S is the number of species in the sample-species matrix, and xij and xik are the abundance of species i in j and k plots, respectively.

The βJ is based on the incidence to measure the compositional dissimilarity between the plots based only on the differences in species occurrence (presence/absence), and it ranges from 0 (when both plots have the same species) to 1 (when the plots do not share any species; Catalán et al., 2023). Podani and Schmera (2011) proposed that the βJ consists of two components: the first component (Rrel = 2 min(b,c)/n) accounts for species replacement, and the second (Drel = |b − c|/n) accounts for richness differences. βJ was calculated using the following equation:

βJ=b+ca+b+c=number of species not sharedtotal number of species=2min(b,c)a+b+c+|bc|a+b+c

2.3.3 Phylogenetic structure

Based on the list of species from the survey, the Chinese Flora (http://www.iplant.cn/frps) and the Plant List website (http://www.theplantlist.org/) were used to standardize the descriptions (Kalwij, 2012). To construct the phylogenetic tree, the awk version of the “phylomatic” software was used, and the phylogenetic tree was visualized using the “ggtree” package in R (Zheng et al., 2019b; Figure 2).

Figure 2
www.frontiersin.org

Figure 2 The phylogenetic tree of the plant species in the survey. The branch is the distance between the species.

To characterize the phylogenetic diversity (PD), the PD index was used, which represents the sum of the lengths of the evolutionary branches of all the species in the community, and the phylogenetic structure was measured using the net relatedness index (NRI) and the net nearest taxon index (NTI). The PD index and phylogenetic structure were calculated using the “picante” package in R (Chai et al., 2016; Qian and Sandel, 2017).

The NRI and NTI were calculated using the following equations:

NRI=1×MPDsMPDmdsSD(MPDmds)
NTI=1×MNTDsMNTDrSD(MNTDr)

where MPD is the mean phylogenetic distance between all the species pairs in the community, and MPDs and MNTDs denote the average observed paired phylogenetic distance and the average observed nearest neighbor phylogenetic distance, respectively. MPDmds and MNTDmds denote the average paired phylogenetic distance and the average nearest neighbor phylogenetic distance, respectively, under 999 null model simulations. SD is the standard deviation.

2.4 Data analysis

The statistical analyses were performed in the R platform (version 3.5.3). The significant differences between the different communities for the diversity indicators were evaluated using the Kruskal–Wallis rank sum test and the “dplyr” library. The overall change in the species composition was measured using a similarity percentage analysis (SIMPER) and the “vegan” package to calculate the contribution of the individual species to the overall community dissimilarity. The β-diversity was separated into two components, species replacement [or turnover (Repl)] and richness differences [or nestedness (RichDiff)], to clarify their relative contribution using the “adespatial” package. To examine whether the ST and NT communities were spatially similar, non-metric multidimensional scaling and hierarchical clustering were performed based on the Raup–Crick distance using the “vegan” library. Then, a correlation analysis was conducted between the phylogenetic and species diversity indices using the “corrplot” package. The roles of the deterministic and stochastic processes in community assembly were assessed using the β-null model (Tucker et al., 2016).

3 Results

3.1 Species composition and diversity

The survey revealed that there were 20 species of seed plants from 12 families and 20 genera identified in the ST communities, with Asteraceae having the most species (20.0%), followed by Amaranthaceae (15.0%). The proportion of the remaining families was less than 10.0%. There were nine herbaceous annuals, nine herbaceous perennials, and two semi-shrub species, with Setaria arenaria and Corispermum candelabrum (Supplementary Table S1). A total of 15 species of seed plants from eight families and 14 genera were identified in the NT communities, among which Asteraceae had the highest number of species (33.3%), which was followed by Amaranthaceae (20.0%) and Gramineae (13.3%), and the proportion of the remaining families was less than 10.0%. There were nine herbaceous annuals, four herbaceous perennials, one semi-shrub species, and one shrub species, with A. ordosica being the dominant species (Supplementary Table S1).

3.1.1 α-diversity

The plant community in the shallow-tilled area had a high species diversity and evenness, with an H, D, J, and richness (S) of 1.68 ± 0.18, 0.72 ± 0.07, 0.66 ± 0.07, and 12.90 ± 2.85, respectively. In the control area, those measures were 0.98 ± 0.44, 0.46 ± 0.20, 0.54 ± 0.16, and 6.10 ± 2.38, respectively, and they were significantly higher (p < 0.05) in the shallow plowed area than in the control area (Table 1).

Table 1
www.frontiersin.org

Table 1 Species diversity indices of the different communities.

3.1.2 β-diversity

The β-diversity analysis that was based on the Bray–Curtis dissimilarity index showed that the mean values of the β-diversity of the plant communities in the ST and NT communities were 0.48 and 0.42, respectively, and there was a highly significant difference between them (p < 0.001, Figure 3). The β-diversity decomposition analyses showed that the species replacement process contributed more (60.82%) to the differences in the plant community composition in the ST, and the richness differences contributed less to the β-diversity (Figure 4A). The differences in the richness and species turnover processes contributed similarly to the differences in the plant community composition in the NT, being 54.68% and 45.32%, respectively (Figure 4B).

Figure 3
www.frontiersin.org

Figure 3 Differences in the Bray–Curtis dissimilarity index between the communities. The horizontal lines indicate the median values, the boxes indicate 25%–75% confidence intervals, the vertical lines indicate 10%–90% confidence intervals, the solid circles indicate outliers, and the different lowercase letters indicate significant differences among the communities. ST, shallow tillage area; NT, non-shallow tillage area.

Figure 4
www.frontiersin.org

Figure 4 Triangular plots of the β-diversity comparisons (using the Jaccard index) for the plant communities among the (A) shallow tillage area (ST) and (B) non-shallow tillage area (NT). Each point represents a pair of sites. Its position is determined by a triplet of values, namely, the similarity (S), replacement (Repl), and richness difference (RichDiff) matrices, with each triplet summing to 1. The mean values of S, Repl, and RichDiff are shown.

The Venn diagram (Supplementary Figure S1) shows the presence of 13 common species between the two communities. Of the reported species in this survey, eight were exclusively identified in the ST (Ferula bungeana, Oxytropis racemose, Polygala tenuifolia, Allium mongolicum, Dracocephalum moldavica, Tribulus terrestris, Panzeria alaschanica, and Cynanchum thesioides), and three were only found in the NT (Tragus mongolorum, Hypecoum erectum, and A. ordosica).

The results of SIMPER analysis showed that C. candelabrum, S. arenaria, and Chenopodium acuminatum contributed the most to the inter-community species differences and that the contribution of the three plants to the inter-community dissimilarity between the ST and the NT was more than 10% (C. candelabrum, 18.44%; S. arenaria, 17.52%; C. acuminatum, 10.57%). The three plants’ cumulative contribution to the inter-community variation amounted to 68.13%, while the contribution of the other four plants to the inter-community differences amounted to 10.57%. The cumulative contribution of the seven plants to the inter-community variation amounted to 86.87% (Table 2).

Table 2
www.frontiersin.org

Table 2 Contribution of the main species to community dissimilarity.

3.2 Phylogenetic structure

The mean values of the PD, NRI, and NTI were 1,251.06 ± 70.94, −0.74 ± 0.03, and −0.06 ± 0.18 in the ST, respectively; the mean values of the PD, NRI, and NTI were 683.26 ± 62.37, −0.89 ± 0.13, and −0.84 ± 0.51 in the NT, respectively. The PD was significantly higher in the ST than in the NT, and the NRI was significantly lower in the ST than in the NT (p < 0.001, Supplementary Table S2). No significant correlation was found between the PD values and species diversity indices in the ST and NT communities. In the ST communities, the NTI was significantly and positively correlated with the H, D, and J (p < 0.001; Figure 5A). The NTI and NRI were significantly and positively correlated with the SR in the NT communities (p < 0.05, Figure 5B).

Figure 5
www.frontiersin.org

Figure 5 Correlation analysis between the species diversity and phylogenetic indices. (A) Shallow tillage area (ST). (B) Non-shallow tillage area (NT). Blue indicates a positive correlation, and orange indicates a negative correlation; the darker the color, the stronger the correlation. PD, phylogenetic diversity; NRI, net relatedness index; NTI, net nearest taxon index; H, Shannon–Weiner index; D, Simpson’s index; J, Pielou’s evenness index; SR, species richness. *p < 0.05; **p < 0.01; ***p < 0.001.

3.3 Community assembly process

The results revealed that the NT communities were more similar, indicating that they deviated from the null hypothesis (Figure 6). However, the ST communities were farther apart and spread over half of the space, suggesting that the null hypothesis may be supported. These results suggest that stochastic processes play a greater role in the ST communities. A similar result was observed for the phylogenetic structure. In the ST and NT communities, the NRI and NTI were both less than 0 (Supplementary Table S2), indicating that the phylogenetic structure of the plant community diverged, there were more distantly related species in the community, the construction of the ST and NT communities was dominated by stochastic processes, and competitive exclusion was the primary mechanism of plant construction. However, when compared with the NT communities, the NRI and NTI were closer to 0 in the ST communities, indicating that the ST communities were subject to stronger stochastic processes.

Figure 6
www.frontiersin.org

Figure 6 Non-metric multidimensional scaling (NMDS) analysis of the species composition in the sample plots. ST, shallow tillage area (yellow circle); NT, non-shallow tillage area (blue circle).

4 Discussion

4.1 Species diversity and species composition

The plant community structure and function are determined by a variety of biotic and abiotic drivers, but anthropogenic environmental changes may alter these drivers and their effects on the plant communities (Brooker, 2006). The Mu Us Desert, as a typical agricultural and pastoral area, is ecologically fragile and frequently disturbed by human activities, and a full understanding of the biodiversity changes before and after disturbance is needed for effective conservation. Studies have confirmed that the SR of woody and herbaceous plants is strongly affected by environmental changes and anthropogenic disturbances (Zhou et al., 2018). Shallow tillage removes almost all the existing above-ground vegetation, increases the ecological niche space for communities (Tortorelli et al., 2022), and alters the physical properties, such as soil compactness, creating open and variable soil conditions and promoting soil seed germination and organic matter decomposition (Lu et al., 2018). Herbaceous plants, whose seeds are generally small and easy to disseminate, can quickly and effectively utilize temporary habitats for reproduction and respond to external disturbances than woody plants, which resulted in a significantly higher species diversity index in the ST communities than in the NT communities. Additionally, the allelopathic effect of A. ordosica significantly inhibits seed germination and seedling growth of herbaceous plants (Lopes et al., 2022). When a layer of up to 80% shrubs is formed, annual and biennial herbaceous plants are disadvantaged in resource competition and are mostly distributed in the unshaded edges of A. ordosica scrub. Thus, competition for available resources, such as light and soil nutrients, influences the maintenance of species diversity (Adler et al., 2018). At the same time, the stochastic nature of seed dispersal produced a preferential effect (Weidlich et al., 2021), with large differences in the species composition among the ST communities, which were clearly differentiated from the NT communities. Analyses of the β-diversity also revealed differences in the species composition between the ST and NT communities, with the βST being dominated by species replacement processes, indicating that the different communities tended to have more endemic species, and differences in abundance played a larger role than species replacement processes in the βNT. This showed that a few communities contained the vast majority of the species in that state.

Several studies have demonstrated a negative correlation between the community restoration time and the SR; with the extension of the restoration time, early colonizing herbaceous plants will largely disappear (Daws et al., 2021; Standish et al., 2021). Ecological communities are always dynamic and vary in space and time, and they often develop along relatively predictable successional trajectories (Mori et al., 2018). We speculate that without re-imposed anthropogenic disturbance, the ST community will become more similar to the NT and form a zonal sub-top-level community with A. ordosica as a single dominant species (Webb et al., 1987; Li et al., 2015), i.e., a convergent process of community succession.

4.2 Phylogenetic structure

The low PD values shown for the NT communities suggested that fewer clades successfully colonized the NT and that the NT has a high frequency of closely related species. This could reflect the fact that phylogenetic diversity is lower in environments with greater competitive pressures for resources, as these environments may represent ecophysiological barriers that are difficult to surmount evolutionarily (Honorio Coronado et al., 2015); however, due to the small sample size and the limited extent of plots, further research is needed.

The increase in the SR may imply more complex interspecific affinities, which in turn leads to increased PD values, but no significant positive correlation between the PD values and the species diversity index was found in both the ST and NT in this study. Most studies have demonstrated a correlation between the two, but the correlation decreases with unbalanced evolutionary trees or narrow species distributions (Gascuel et al., 2015). When the ST and NT communities were viewed as a whole, the PD, NTI, and NRI were significantly positively correlated with the H, D, and S (Supplementary Figure S2). This significant relationship indicates that the species composition and distribution have some influence on the phylogenetic structure.

4.3 Community assembly

The phylogenetic structure of both the ST and NT communities tended to diverge (NTI < 0 and NRI < 0), suggesting that competitive exclusion is a dominant mechanism influencing community assembly. This pattern may be attributed to the limited availability of soil water and nutrients in the Mu Us Desert, which likely prevents closely related species with similar ecological niches from coexisting due to competition for these scarce resources, ultimately leading to the strengthening of interspecific ecological niche differentiation and the dispersion of phylogenetic structure (Webb et al., 2002; Silvertown, 2004). It is important to note that shallow tillage did not alter the major ecological processes driving community assembly, but the relative importance of deterministic and stochastic processes differed between the two treatments. The NRI and NTI of the ST communities differed significantly from those of the NT communities (p < 0.05, Supplementary Table S2), and analysis based on the Raup–Crick dissimilarity index also showed that the βRC within the ST communities was closer to the zero expectation than in the NT communities (Supplementary Figure S4), suggesting that in addition to being dominated by competitive exclusion, the ST was also largely influenced by stochastic processes dominated by dispersal limitation, which resulted in lower species similarity within the communities. Studies have confirmed that disturbance may promote stochastic processes (Didham et al., 2005; Didham and Norton, 2006).

The desert maintains environmental consistency at large scales (Teramoto et al., 2022), and it can be assumed that the effective seed bank is the result of deterministic abiotic factors, such as temperature or precipitation, and stochastic factors, such as dispersal constraints. In contrast, habitat sieving at small scales mainly consists of biotic factors, such as interspecific interactions (Royo and Ristau, 2013; Zou et al., 2021). As A. ordosica requires a large amount of resources as a community-building species, including light, it results in a limited number of ecological niches being available in the community, which prevents certain species from persisting (Deng et al., 2022). In contrast, under shallow tillage conditions, where the community species are essentially removed and available ecological niches proliferate, herbaceous seeds can germinate rapidly and complete their life history cycles when the climatic conditions are favorable (Grime, 1977). Under these conditions, rapid reproductive processes may reduce the role of severe constraints on the hydrothermal conditions and inter- and intraspecific competition in shaping community membership. In future investigations, more attention should be paid to the role of annual and biennial herbaceous plants in the process of community assembly under human interference.

5 Conclusions

Shallow tillage increases species alpha, beta, and phylogenetic diversity and reduces the importance of competitive exclusion in community assembly. Therefore, we suggest that plant species with low ecological niche overlap should be selected for artificial vegetation restoration in the Mu Us Desert to reduce competition for the same resources.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

ZL: Conceptualization, Investigation, Data curation, Writing – original draft, Writing – review & editing. JQ: Writing – original draft, Writing – review & editing. ZL: Investigation, Writing – review & editing. XG: Investigation, Writing – review & editing. GH: Investigation, Writing – review & editing. HY: Investigation, Writing – review & editing. EH: Investigation, Writing – review & editing. CL: Investigation, Writing – review & editing. XW: Writing – review & editing. GL: Writing – review & editing, Funding acquisition. RG: Writing – review & editing, Supervision.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work is supported by the Key Special Project for the Action of “Revitalizing Inner Mongolia through Science and Technology” (Grant No. 2022EEDSKJXM003-1), Natural Science Foundation of Inner Mongolia Autonomous Region (Grant No. 2024QN03025), Inner Mongolia Autonomous Region Key R&D and Achievement Transformation Program Project (Grant No. 2022YFHH0106), National “14th Five-Year Plan” Key R&D Program (No. 2022YFF1302503-3), Scientific Research Capacity Enhancement Project of Inner Mongolia Academy of Forestry Sciences (Grant No. 2024NLTS01 and 2024NLTS03), and National Natural Science Foundation of China (Grant No. 32360422).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1409493/full#supplementary-material

References

Adler, P. B., Smull, D., Beard, K. H., Choi, R. T., Furniss, T., Kulmatiski, A., et al. (2018). Competition and coexistence in plant communities: intraspecific competition is stronger than interspecific competition. Ecol. Lett. 21, 1319–1329. doi: 10.1111/ele.13098

PubMed Abstract | CrossRef Full Text | Google Scholar

Alarcón, R., Hernández-Plaza, E., Navarrete, L., Sánchez, M. J., Escudero, A., Hernanz, J. L., et al. (2018). Effects of no-tillage and non-inversion tillage on weed community diversity and crop yield over nine years in a mediterranean cereal-legume cropland. Soil Tillage Res. 179, 54–62. doi: 10.1016/j.still.2018.01.014

CrossRef Full Text | Google Scholar

Barroso, J., Miller, Z. J., Lehnhoff, E. A., Hatfield, P. G., Menalled, F. D. (2015). Impacts of cropping system and management practices on the assembly of weed communities. Weed Research 55, 426–435. doi: 10.1111/wre.12155

CrossRef Full Text | Google Scholar

Booth, B. D., Swanton, C. J. (2002). Assembly theory applied to weed communities. Weed Science 50, 2–13. doi: 10.1614/0043-1745(2002)050[0002:AIATAT]2.0.CO;2

CrossRef Full Text | Google Scholar

Brooker, R. W. (2006). Plant–plant interactions and environmental change. New Phytol. 171, 271–284. doi: 10.1111/j.1469-8137.2006.01752.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cadotte, M. W., Cardinale, B. J., Oakley, T. H. (2008). Evolutionary history and the effect of biodiversity on plant productivity. Proc. Natl. Acad. Sci. U.S.A. 105, 17012–17017. doi: 10.1073/pnas.0805962105

PubMed Abstract | CrossRef Full Text | Google Scholar

Cadotte, M. W., Cavender-Bares, J., Tilman, D., Oakley, T. H. (2009). Using phylogenetic, functional and trait diversity to understand patterns of plant community productivity. PloS One 4, e5695. doi: 10.1371/journal.pone.0005695

PubMed Abstract | CrossRef Full Text | Google Scholar

Capmourteres, V., Anand, M. (2016). Conservation value”: a review of the concept and its quantification. Ecosphere 7, e01476. doi: 10.1002/ecs2.1476

CrossRef Full Text | Google Scholar

Catalán, A. M., López, D. N., Fica-Rojas, E., Broitman, B. R., Valdivia, N., Scrosati, R. A. (2023). Foundation species canopies affect understory beta diversity differently depending on species mobility. Ecology 104, e3999. doi: 10.1002/ecy.3999

PubMed Abstract | CrossRef Full Text | Google Scholar

Chai, Y., Yue, M., Liu, X., Guo, Y., Wang, M., Xu, J., et al. (2016). Patterns of taxonomic, phylogenetic diversity during a long-term succession of forest on the Loess Plateau, China: insights into assembly process. Sci. Rep. 6, 27087. doi: 10.1038/srep27087

PubMed Abstract | CrossRef Full Text | Google Scholar

Chase, J. M. (2010). Stochastic community assembly causes higher biodiversity in more productive environments. Science 328, 1388–1391. doi: 10.1126/science.1187820

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Lu, H., Wu, H., Wang, J., Lyu, N. (2023). Global desert variation under climatic impact during 1982–2020. Sci. China Earth Sci. 66, 1062–1071. doi: 10.1007/s11430-022-1052-1

CrossRef Full Text | Google Scholar

Chen, Z., Huang, M., Xiao, C., Qi, S., Du, W., Zhu, D., et al. (2022). Integrating remote sensing and spatiotemporal analysis to characterize artificial vegetation restoration suitability in desert areas: A Case Study of Mu Us Sandy Land. Remote Sens. 14, 4736. doi: 10.3390/rs14194736

CrossRef Full Text | Google Scholar

Cooper, J., Baranski, M., Stewart, G., Nobel-de Lange, M., Bàrberi, P., Fließbach, A., et al. (2016). Shallow non-inversion tillage in organic farming maintains crop yields and increases soil C stocks: a meta-analysis. Agron. Sustain. Dev. 36, 22. doi: 10.1007/s13593-016-0354-1

CrossRef Full Text | Google Scholar

Cordeau, S., Smith, R. G., Gallandt, E. R., Brown, B., Salon, P., DiTommaso, A., et al. (2017). Disentangling the effects of tillage timing and weather on weed community assembly. Agriculture 7, 66. doi: 10.3390/agriculture7080066

CrossRef Full Text | Google Scholar

Czerwiński, M., Woodcock, B. A., Golińska, B., Kotowski, W. (2018). The effect of tillage management and its interaction with site conditions and plant functional traits on plant species establishment during meadow restoration. Ecol. Eng. 117, 28–37. doi: 10.1016/j.ecoleng.2018.03.017

CrossRef Full Text | Google Scholar

Daws, M. I., Walters, S. J., Harris, R. J., Tibbett, M., Grigg, A. H., Morald, T. K., et al. (2021). Nutrient enrichment diminishes plant diversity and density, and alters long-term ecological trajectories, in a biodiverse forest restoration. Ecol. Eng. 165, 106222. doi: 10.1016/j.ecoleng.2021.106222

CrossRef Full Text | Google Scholar

Deng, Z., Zhao, J., Wang, Z., Li, R., Guo, Y., Luo, T., et al. (2022). Stochastic processes drive plant community assembly in alpine grassland during the restoration period. Diversity 14, 832. doi: 10.3390/d14100832

CrossRef Full Text | Google Scholar

Didham, R. K., Norton, D. A. (2006). When are alternative stable states more likely to occur?: A reply to T. Fukami and W. G. Lee. Oikos 113, 357–362. doi: 10.1111/j.0030-1299.2006.13508.x

CrossRef Full Text | Google Scholar

Didham, R. K., Watts, C. H., Norton, D. A. (2005). Are systems with strong underlying abiotic regimes more likely to exhibit alternative stable states? Oikos 110, 409–416. doi: 10.1111/j.0030-1299.2005.13883.x

CrossRef Full Text | Google Scholar

Durant, S. M., Pettorelli, N., Bashir, S., Woodroffe, R., Wacher, T., De Ornellas, P., et al. (2012). Forgotten biodiversity in desert ecosystems. Science 336, 1379–1380. doi: 10.1126/science.336.6087.1379

PubMed Abstract | CrossRef Full Text | Google Scholar

Erfanian, M. B., Ejtehadi, H., Vaezi, J., Moazzeni, H. (2019). Plant community responses to multiple disturbances in an arid region of northeast Iran. Land Degrad Dev. 30, 1554–1563. doi: 10.1002/ldr.3341

CrossRef Full Text | Google Scholar

Gascuel, F., Ferrière, R., Aguilée, R., Lambert, A. (2015). How ecology and landscape dynamics shape phylogenetic trees. Systematic Biol. 64, 590–607. doi: 10.1093/sysbio/syv014

CrossRef Full Text | Google Scholar

Grime, J. P. (1977). Evidence for the existence of three primary strategies in plants and its relevance to ecological and evolutionary theory. Am. Nat. 111, 1169–1194. doi: 10.1086/283244

CrossRef Full Text | Google Scholar

Heino, J., Tolonen, K. T. (2017). Untangling the assembly of littoral macroinvertebrate communities through measures of functional and phylogenetic alpha diversity. Freshw. Biol. 62, 1168–1179. doi: 10.1111/fwb.12934

CrossRef Full Text | Google Scholar

Honorio Coronado, E. N., Dexter, K. G., Pennington, R. T., Chave, J., Lewis, S. L., Alexiades, M. N., et al. (2015). Phylogenetic diversity of Amazonian tree communities. Diversity Distributions 21, 1295–1307. doi: 10.1111/ddi.12357

CrossRef Full Text | Google Scholar

Jaureguiberry, P., Titeux, N., Wiemers, M., Bowler, D. E., Coscieme, L., Golden, A. S., et al. (2022). The direct drivers of recent global anthropogenic biodiversity loss. Sci. Adv. 8, eabm9982. doi: 10.1126/sciadv.abm9982

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, H., Xu, J., Peng, Y., Xin, J., Peng, N., Li, Y., et al. (2023). Impacts of landscape patterns on plant species diversity at a global scale. Sci. Total Environ. 896, 165193. doi: 10.1016/j.scitotenv.2023.165193

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalwij, J. M. (2012). Review of ‘The Plant List, a working list of all plant species.’. J. Vegetation Sci. 23, 998–1002. doi: 10.1111/j.1654-1103.2012.01407.x

CrossRef Full Text | Google Scholar

Laudicina, V. A., Palazzolo, E., Catania, P., Vallone, M., García, A. D., Badalucco, L. (2017). Soil quality indicators as affected by shallow tillage in a vineyard grown in a semiarid mediterranean environment. Land Degradation Dev. 28, 1038–1046. doi: 10.1002/ldr.2581

CrossRef Full Text | Google Scholar

Légère, A., Stevenson, F. C., Benoit, D. L. (2005). Diversity and assembly of weed communities: contrasting responses across cropping systems. Weed Research 45, 303–315. doi: 10.1111/j.1365-3180.2005.00459.x

CrossRef Full Text | Google Scholar

Li, S.-L., Yu, F.-H., Werger, M. J., Dong, M., During, H. J., Zuidema, P. A. (2015). Mobile dune fixation by a fast-growing clonal plant: a full life-cycle analysis. Sci. Rep. 5, 8935. doi: 10.1038/srep08935

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes, A. D., Nunes, M. G. I. F., Francisco, J. P., dos Santos, E. H. (2022).“Potential allelopathic effect of species of the Asteraceae family and its use in agriculture,”. In: Vegetation dynamics, changing ecosystems and human responsibility (IntechOpen). Available online at: https://www.intechopen.com/chapters/84835 (Accessed May 29, 2024).

Google Scholar

Lu, D., Wang, G. G., Yu, L., Zhang, T., Zhu, J. (2018). Seedling survival within forest gaps: the effects of gap size, within-gap position and forest type on species of contrasting shade-tolerance in Northeast China. Forestry: Int. J. For. Res. 91, 470–479. doi: 10.1093/forestry/cpy007

CrossRef Full Text | Google Scholar

Mishler, B. D., Knerr, N., González-Orozco, C. E., Thornhill, A. H., Laffan, S. W., Miller, J. T. (2014). Phylogenetic measures of biodiversity and neo-and paleo-endemism in Australian Acacia. Nat. Commun. 5, 4473. doi: 10.1038/ncomms5473

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori, A. S., Isbell, F., Seidl, R. (2018). β-diversity, community assembly, and ecosystem functioning. Trends Ecol. Evol. 33, 549–564. doi: 10.1016/j.tree.2018.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Podani, J., Schmera, D. (2011). A new conceptual and methodological framework for exploring and explaining pattern in presence–absence data. Oikos 120, 1625–1638. doi: 10.1111/j.1600-0706.2011.19451.x

CrossRef Full Text | Google Scholar

Qian, H., Sandel, B. (2017). Phylogenetic relatedness of native and exotic plants along climate gradients in California, USA. Diversity Distributions 23, 1323–1333. doi: 10.1111/ddi.12620

CrossRef Full Text | Google Scholar

Randriamalala, J. R., Hervé, D., Randriamboavonjy, J.-C., Carrière, S. M. (2012). Effects of tillage regime, cropping duration and fallow age on diversity and structure of secondary vegetation in Madagascar. Agriculture Ecosyst. Environ. 155, 182–193.

Google Scholar

Royo, A. A., Ristau, T. E. (2013). Stochastic and deterministic processes regulate spatio-temporal variation in seed bank diversity. J. Vegetation Sci. 24, 724–734. doi: 10.1111/jvs.12011

CrossRef Full Text | Google Scholar

She, W., Bai, Y., Zhang, Y., Qin, S., Liu, Z., Wu, B. (2017). Plasticity in meristem allocation as an adaptive strategy of a desert shrub under contrasting environments. Front. Plant Sci. 8, 1933.

PubMed Abstract | Google Scholar

Silvertown, J. (2004). Plant coexistence and the niche. Trends Ecol. Evol. 19, 605–611. doi: 10.1016/j.tree.2004.09.003

CrossRef Full Text | Google Scholar

Standish, R. J., Gove, A. D., Grigg, A. H., Daws, M. I. (2021). Beyond species richness and community composition: Using plant functional diversity to measure restoration success in jarrah forest. Appl. Vegetation Sci. 24, e12607. doi: 10.1111/avsc.12607

CrossRef Full Text | Google Scholar

Teramoto, M., Hamamoto, T., Liang, N., Taniguchi, T., Ito, T. Y., Hu, R., et al. (2022). Abiotic and biotic factors controlling the dynamics of soil respiration in a coastal dune ecosystem in western Japan. Sci. Rep. 12, 14320. doi: 10.1038/s41598-022-17787-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Tortorelli, C. M., Kerns, B. K., Krawchuk, M. A. (2022). Community invasion resistance is influenced by interactions between plant traits and site productivity. Ecology 103, e3697. doi: 10.1002/ecy.3697

PubMed Abstract | CrossRef Full Text | Google Scholar

Travlos, I. S., Cheimona, N., Roussis, I., Bilalis, D. J. (2018). Weed-species abundance and diversity indices in relation to tillage systems and fertilization. Front. Environ. Sci. 6. doi: 10.3389/fenvs.2018.00011

CrossRef Full Text | Google Scholar

Tucker, C. M., Cadotte, M. W. (2013). Unifying measures of biodiversity: understanding when richness and phylogenetic diversity should be congruent. Diversity Distributions 19, 845–854. doi: 10.1111/ddi.12087

CrossRef Full Text | Google Scholar

Tucker, C. M., Shoemaker, L. G., Davies, K. F., Nemergut, D. R., Melbourne, B. A. (2016). Differentiating between niche and neutral assembly in metacommunities using null models of β-diversity. Oikos 125, 778–789. doi: 10.1111/oik.02803

CrossRef Full Text | Google Scholar

Valone, T. J. (2003). Examination of interaction effects of multiple disturbances on an arid plant community. Southwestern Nat. 48, 481–490. doi: 10.1894/0038-4909(2003)048<0481:eoieom>2.0.co;2

CrossRef Full Text | Google Scholar

Wang, X., Lou, J., Ma, W., Jiao, L., Li, H., Hua, T. (2017). The impact of reclamation on aeolian desertification of four species in the otindag desert, china. CATENA 157, 189–194. doi: 10.1016/j.catena.2017.05.024

CrossRef Full Text | Google Scholar

Webb, C. O., Ackerly, D. D., McPeek, M. A., Donoghue, M. J. (2002). Phylogenies and community ecology. Annu. Rev. Ecol. Syst. 33, 475–505. doi: 10.1146/annurev.ecolsys.33.010802.150448

CrossRef Full Text | Google Scholar

Webb, R. H., Steiger, J. W., Turner, R. M. (1987). Dynamics of Mojave Desert shrub assemblages in the panamint mountains, California: ecological archives E068-001. Ecology 68, 478–490. doi: 10.2307/1938453

CrossRef Full Text | Google Scholar

Wei, H., Luo, T., Wu, B. (2016). Optimal balance of water use efficiency and leaf construction cost with a link to the drought threshold of the desert steppe ecotone in northern China. Ann. Bot. 118, 541–553. doi: 10.1093/aob/mcw127

PubMed Abstract | CrossRef Full Text | Google Scholar

Weidlich, E. W. A., Nelson, C. R., Maron, J. L., Callaway, R. M., Delory, B. M., Temperton, V. M. (2021). Priority effects and ecological restoration. Restor. Ecol. 29, e13317. doi: 10.1111/rec.13317

CrossRef Full Text | Google Scholar

Zhang, L., Hong, G., Li, Z., Gao, X., Wu, Y., Wang, X., et al. (2018). Assessment of the ecosystem service function of sandy lands at different times following aerial seeding of an endemic species. Sustainability 10, 902. doi: 10.3390/su10040902

CrossRef Full Text | Google Scholar

Zhang, T., Chen, Y. (2022). The effects of landscape change on habitat quality in arid desert areas based on future scenarios: Tarim River Basin as a case study. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1031859

CrossRef Full Text | Google Scholar

Zhang, T., Chen, Y., Ali, S. (2023). Abiotic stress and human activities reduce plant diversity in desert riparian forests. Ecol. Indic. 152, 110340. doi: 10.1016/j.ecolind.2023.110340

CrossRef Full Text | Google Scholar

Zheng, Y., Dong, L., Li, Z., Zhang, J., Li, Z., Miao, B., et al. (2019b). Phylogenetic structure and formation mechanism of shrub communities in arid and semiarid areas of the Mongolian Plateau. Ecol. Evol. 9, 13320–13331. doi: 10.1002/ece3.5787

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, D., Wang, Y., Shao, Y., Wang, L. (2019a). The vegetation dynamics and climate change responses by leaf area index in the Mu Us Desert. Sustainability 11, 3151. doi: 10.3390/su11113151

CrossRef Full Text | Google Scholar

Zhou, Y., Chen, S., Hu, G., Mwachala, G., Yan, X., Wang, Q. (2018). Species richness and phylogenetic diversity of seed plants across vegetation zones of Mount Kenya, East Africa. Ecol. Evol. 8, 8930–8939. doi: 10.1002/ece3.4428

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, C., Martini, F., Xia, S.-W., Castillo-Diaz, D., Goodale, U. M. (2021). Elevation and micro environmental conditions directly and indirectly influence forests’ soil seed bank communities. Global Ecol. Conserv. 26, e01443. doi: 10.1016/j.gecco.2020.e01443

CrossRef Full Text | Google Scholar

Keywords: shallow tillage, species diversity, phylogenetic structure, community assembly, Mu Us Desert

Citation: Li Z, Qiao J, Li Z, Gao X, Hong G, Yang H, Hu E, Liu C, Wang X, Liu G and Gao R (2024) Shallow tillage mitigates plant competition by increasing diversity and altering plant community assembly process. Front. Plant Sci. 15:1409493. doi: 10.3389/fpls.2024.1409493

Received: 30 March 2024; Accepted: 15 July 2024;
Published: 07 August 2024.

Edited by:

Praveen Dhyani, University of Koblenz, Germany

Reviewed by:

Claudia Bita-Nicolae, Institute of Biology Bucharest of the Romanian Academy, Romania
Qicong Wu, Shandong Agricultural University, China

Copyright © 2024 Li, Qiao, Li, Gao, Hong, Yang, Hu, Liu, Wang, Liu and Gao. 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.

*Correspondence: Guanzhi Liu, Y3p5bGd6QGltYXUuZWR1LmNu; Runhong Gao, Z3JoendkbUAxNjMuY29t

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.