Species richness and Phylogenetic structure
Mean species richness was 16.4 species m− 2 in the control plots. DLP and DP reduced species richness by an average of 3.2 (19.5%; p < 0.01) and 3.8 species m− 2 (23.2%; p < 0.001), respectively. In contrast, DEP had no significant effect on species richness (Fig. 1a). ILP caused an increase in species richness of 3.8 species m− 2 (23.2%; p < 0.05), whereas neither IEP nor IP affected species richness.
NRI was lower than zero, indicating significant phylogenetic overdispersion (Fig. 1b). Changing precipitation had no effect on NRI, and average phylogenetic relatedness was stable among treatments (Fig. 1b). Functional NRI was greater than zero, indicating significant functional clustering (Fig. 1c). IEP and ILP decreased functional NRI and shifted phylogenetic relatedness from clustering to randomness. In contrast, DEP and DP significantly reduced functional NRI but did not alter the pattern of phylogenetic clustering. DLP and IP did not alter the functional NRI and average phylogenetic relatedness among species (Fig. 1c).
Species Colonization And Loss
An average of 4.6 (SE = 0.5) species lost in the control plots from 2015 to 2020 (Fig. 2). Some of the grasses (Agropyron cristatum and Leymus chinensis), forbs (e.g., Artemisia pubescens and Chamaerhodos erecta), and legumes (Medicago ruthenica and Gueldenstaedtia stenophylla) with lower initial abundances, sizes, and leaf dry matter content disappeared under DP (Fig. 3a, b, c, d), leading to the loss of 6.8 (SE = 0.4) species. In contrast, the number of species loss under DEP, DLP, IEP, ILP, and IP did not differ from the control. An average of 2.8 (SE = 0.5) species colonized the control plots from 2015 to 2020 (Fig. 2). DLP decreased the colonization of drought sensitive forbs (e.g., Chenopodium glaucum) by 1.4 (SE = 0.4) species m− 2 relative to the control from 2015 to 2020 (Fig. 2). In contrast, ILP increased forbs (Artemisia scoparia and Potentilla acaulis) and legumes (Gueldenstaedtia stenophylla and Melissitus ruthenica) colonization by 2.2 species m− 2 relative to the control from 2015 to 2020. DEP, DP, IEP, and IP had no effect on species colonization. Species with lower initial coverage in 2015 had greater increased likelihoods in 2020 (Fig. 3g), whereas traits, such as species height, leaf nitrogen content, leaf dry matter content, SLA, and sizes, were not good predictors of increased probability in species abundance (Fig. 3h, i, j, k, l).
βNRI between locally lost species and other resident species was significantly lower than zero in the control, IEP, ILP, and IP plots, indicating that lost species in this treatment were more distantly related to the other residents than expected by chance (Fig. 4a). DEP, DLP, and DP significantly increased βNRI but did not alter phylogenetic relatedness between lost species and other resident species. βNRI between colonists and residents in the control, DLP, ILP, and IP treatments was significantly lower than zero, indicating that colonists in these treatments were more distantly related to the residents than expected by chance. βNRI between colonists and residents in DEP, DP, and IEP did not differ from zero, indicating the phylogenetic randomness of species colonization (Fig. 4a).
Functional βNRI between lost species and other resident species was significantly greater than zero in the control and IEP, indicating that the traits of lost species in this treatment were more closely related to other residents than expected by chance (Fig. 4b). DP, ILP, and IP significantly and DLP marginally reduced the functional βNRI of lost species, which caused functional βNRI to be significantly lower than zero in DP, ILP, and IP; however, functional βNRI did not significantly differ from zero in DLP, indicating that the traits of lost species were more distantly related to other residents in DP, ILP, and IP but were random in DLP (Fig. 4b). Functional βNRI between colonists and residents in the control, DLP, DP, and IP did not differ from zero, indicating that the traits of colonists in these treatments were random. Functional βNRI in DEP, IEP, and ILP was significantly lower than zero, indicating that the traits of colonists in these treatments were more distantly related to the residents than expected by chance (Fig. 4b).
Factors Affecting Phylogenetic And Functional Community Structure
Structural equation modeling (SEM) showed that changing precipitation affected species richness directly by altering soil moisture and indirectly by altering the number of colonizing species. Changing precipitation affected functional NRI by regulating soil moisture and thus the number of colonizing species and species richness (Fig. 6). Functional phylogenetic relatedness was independent of structural phylogenetic relatedness, and changes in species richness did not affect phylogenetic relatedness (Fig. 6).