Abstract Aims The community succession theory is much debated in ecology. We studied succession on Zokor rodent mounds on the Tibetan Plateau to address several fundamental questions, among them: (i) During secondary succession, does the community composition converge towards one community state or multiple states depending on the initial colonization? (ii) Do mound communities located in different background communities exhibit different assembly trajectories? Methods In a sub-alpine meadow, we investigated a total of 80 mound communities at several successional stages in three different background communities resulting from different management histories and compared their changes in species composition. The distribution of plant communities over time was analyzed with quantitative classification and ordination methods. The co-occurrence patterns of species were evaluated at each successional stage, and the degree of convergence/divergence among communities was obtained by calculating two beta-diversity indices. Important Findings During secondary succession, species richness of mound communities changed over time, and this change was dependent on the background community. Five life-form groups exhibited different dynamic patterns in species richness and plant cover. Community composition and the degree of species co-occurrence between communities increased over time since disturbance. There was much variation in species composition at earlier stages of succession, but communities on older mounds became more similar to each other and to their surrounding vegetation over the course of secondary succession. Post-disturbance succession of Zokor mound communities transitioned from ‘multiple alternative states’ to ‘background-based deterministic community assembly’ over time. Tradeoffs between competition and colonization, as well as the characteristics of different life-forms and mass effects within a limited species pool are the mechanisms responsible for convergence of mound communities. multiple alternative states, background-based deterministic community assembly, sub-alpine meadow, disturbance, plant communities INTRODUCTION Climatic determinism and historical contingency are alternative theories that have been proposed to explain patterns of succession: plant community change over time (Lepš and Rejmánek 1991; Samuels and Drake 1997; Suding et al. 2004; Suding and Hobbs 2009; Young et al. 2001). Deterministic climax refers to successional convergence of plant communities within a region on one state determined by environmental factors (Clements 1916), primarily climatic conditions and soil parent material (Lanta and Lepš 2009; Lepš and Rejmánek 1991). Historical contingency, on the other hand, suggests that plant communities can diverge towards multiple states due to stochastic factors (Egler 1954), such as the order and timing of species arriving at different localities, even if environmental conditions and species pool are similar (Chase 2010; Fukami 2010; Fukami et al. 2005). When an ecosystem is severely degraded, the original community can move to an alternative state after disturbance that causes niche release (Erfanzadeh et al. 2010; Fukami and Lee 2006; Saccone et al. 2014), after which priority effects affect community assembly process through niche preemption and modification. In some productive ecosystems, regional diverse species pool and local rapid revegetation dynamics together can also provide the conditions for divergent community development (Chase 2003a, 2003b). In contrast to alternative equilibrium states, increasing similarity during the course of community development has been reported in long-term primary succession on lava flows (Cutler 2010) and in secondary succession of tropical forests (Lebrija-Trejos et al. 2010). To date, most experiments testing deterministic versus alternative trajectories in succession have focused on restoration of abandoned lands (Fukami et al. 2005; Lebrija-Trejos et al. 2010), damaged or degraded ecosystems (Beals et al. 2014; Erfanzadeh et al. 2010; Larios et al. 2013) or other specific situations (Fridley 2013; do Nascimento et al. 2014; Yang et al. 2013), but few studies have investigated small-scale disturbance in a degraded system (Lanta and Lepš 2009). In the present study, we investigated the patterns of small-scale secondary succession under different land-use backgrounds using quantitative methods. On the Tibetan Plateau, one quarter of the natural grasslands are seriously damaged by rodents, which consume much vegetation and alter the soil physically (Wang et al. 2008). Plateau zokor (Eospalax baileyi syn. Myospalax baileyi), is one of the two burrowing mammals in the high-elevation interior of Qinghai-Tibetan Plateau. It feeds primarily on plant roots at a depth of 3–20 cm below the soil surface (Kang et al. 2007) and builds soil mounds on which succession is initiated. It has been estimated that one plateau zokor can produce ~240 mounds each year with a total mound surface area of around 22.5 m2 of grassland (Kang et al. 2007). These mounds are ideal systems to investigate changes of species composition during early succession. Sub-alpine meadow ecosystems on the Tibetan Plateau have been experiencing serious degradation due to overgrazing and planting of grasses or crops. Several attempts at ecological restoration are changing the community structure. For example, enclosure fencing has been established to prevent intensive grazing. These local interventions mean that secondary succession of zokor mounds is not occurring within a single type of climax or sub-climax ‘background’ or ‘matrix community’ vegetation (Zhang and Liu 2003). We sampled the local communities within the fenced area to investigate the influence of the local community on succession on mounds. In the present study, our aims are to (i) describe the dynamics of species replacement, co-occurrence patterns and species diversity during succession on Zokor mounds, (ii) determine effects of different background communities on the development of local communities, (iii) test the alternative models of succession and (iv) ask which potential mechanism are consistent with the patterns observed. MATERIALS AND METHODS Study area Our study was conducted at the Research Station of Alpine Meadow and Wetland Ecosystems of Lanzhou University (34°55ʹN, 102°53ʹE; 2900 m.a.s.l.), located on the northeast part of the Tibetan Plateau, China. The region has a continental dry climate, with a mean annual temperature of ~2°C, mean annual precipitation of 550 mm, and a plant growth period of 120–140 days from May to September. The topography is flat, and the vegetation is classified as sub-alpine meadow, dominated by Artemisia tangutica, Trigonella ruthenica and Potentilla anserina. The soil is classified as a chestnut soil or a typical calcareous soil (Table 1). The vegetation has deteriorated due to digging by rodents and overgrazing by antelope and yaks. Table 1: the physical and chemical characteristics of soil from the three background communities in areas A, B and C Parameters Areas A B C pH 7.2 ± 0.03 7.22 ± 0.07 7.43 ± 0.06 Organic carbon mg*g−1 36.79 ± 0.75 33.80 ± 0.9 32.49 ± 1.29 Total nitrogen mg*g−1 4.02 ± 0.02 3.88 ± 0.06 3.65 ± 0.06 Total phosphorus mg*g−1 0.41 ± 0.01 0.4 ± 0.01 0.38 ± 0.01 Soil moisture % 27.52 ± 0.55 25.4 ± 0.54 20.58 ± 0.54 Bulk density g*cm−3 1.41 ± 0.01 1.45 ± 0.02 1.51 ± 0.05 Parameters Areas A B C pH 7.2 ± 0.03 7.22 ± 0.07 7.43 ± 0.06 Organic carbon mg*g−1 36.79 ± 0.75 33.80 ± 0.9 32.49 ± 1.29 Total nitrogen mg*g−1 4.02 ± 0.02 3.88 ± 0.06 3.65 ± 0.06 Total phosphorus mg*g−1 0.41 ± 0.01 0.4 ± 0.01 0.38 ± 0.01 Soil moisture % 27.52 ± 0.55 25.4 ± 0.54 20.58 ± 0.54 Bulk density g*cm−3 1.41 ± 0.01 1.45 ± 0.02 1.51 ± 0.05 View Large Experimental design and data collection A 300 × 300 m field has been fenced to control grazing since 2001, but 15 days of grazing by antelope and yaks each winter has been permitted. The field was equally divided into three 300 × 100 m areas (A, B and C) representing 3 different vegetation management histories (Fig. 1). Areas A and B were managed for domination by native perennial grasses until 2002 according to government guidelines, after which the natural succession was permitted. Before abandonment, Festuca sinensis was planted on area A and Leymus secalinus on area B. F. sinensis is common in later native primary succession, while L. secalinus is considered a mid-successional species. In contrast, area C has been allowed to undergo natural succession since the enclosures were established in 2001. Fifty-eight plant species were identified on area A, 71 on B and 80 on C. Figure 1: View largeDownload slide the three experimental areas on the Tibetan Plateau, China. Areas A was abandoned and dominated by a later-successional native species Festuca sinensis, area B was abandoned and dominated by a mid-successional species Leymus secalinus and area C was experienced natural succession since 2001. Figure 1: View largeDownload slide the three experimental areas on the Tibetan Plateau, China. Areas A was abandoned and dominated by a later-successional native species Festuca sinensis, area B was abandoned and dominated by a mid-successional species Leymus secalinus and area C was experienced natural succession since 2001. Fieldwork was carried out in the summer of 2009. A 50 × 50 m plot was placed in the center of each area. To investigate the severity of zokor activity under different management regimes, we have followed mound-making activity in the three plots since 2003. Zokor mounds were labelled with the year they were created. In July 2009, 4–6 mounds of each of 5 age groups were selected and labeled: (a) 1 week to 3 months old, (b) 1-, (c) 2-, (d) 3- and (e) 4-year-old mounds. It was not always possible to ﬁnd six zokor mounds at each stage in all three areas. We selected and marked 27 mounds in area A, 26 in B and 27 in C. All selected mounds were of similar size (~2500 cm2) and intact since the time the rodents built them. In addition, six 50 × 50 cm2 permanent quadrats (f) that did not include any mounds were randomly placed in each plot to sample the background community. All vascular plants rooted in the quadrats were recorded, and the ground cover of each species estimated visually. Due to logistical limitations, each management regime was applied in only one field, so, like many studies on succession, our experiment design has no replication at the background community treatment level. Statistical analyses We analyzed species composition, relative cover and species richness and five life forms: annual forbs (AF), biennial forbs (BF), perennial forbs (PF), leguminous plants (L) and gramineous plants (G) on each mound. Regression analysis was used to test effects of the successional stages on species richness, and the relative cover of mound communities. We used quantitative classification and ordination methods to analyze the overall structure and distribution pattern of mound communities during secondary succession. TWINSPAN (two-way indicator species analysis) in program WinTwins (Hill and Šmilauer 2005) was used to identify major vegetation types from these data. Because the relative cover values of plant species ranged from 0.1% to 58.8%, the cutoff levels were defined as 0, 5%, 20% and 50% in our TWINSPAN analysis. Detrended correspondence analysis (DCA) was used to compare the overall changes of species composition at each successional stage of mound vegetation and the reference community (Oksanen and Minchin 1997). Principal component analysis was not used because it exhibits a strong ‘horseshoe effects’ (Legendre and Legendre 2012). Ordination was performed with the ‘vegan’ package (Oksanen et al. 2011) in the R statistical program (version 2.13, R Development Core Team 2011). The species dataset was down-weighted to reduce the influence of rare species. The checkerboard score (C-score; Gotelli and McCabe 2002; Oksanen et al. 2011) was used to estimate co-occurrence patterns among colonizing species within each mound age group. The C-score was obtained by calculating the average number of C units in a species presence/absence matrix. The C unit (cu) for each species pair was calculated as cui j=(Oi−S)(Oj−S), where Oi and Oj are the total number of quadrats occupied by species i and j, and S is the number of quadrats occupied by both species. A null model was established to ask whether the observed C-score differs from random expectations according to the method of Cutler (2010). A C-score greater than null model expectations is support for negative co-occurrence patterns (Diamond 1975), whereas a C-score less than that of the null model is support for positive co-occurrence. All calculations used functions loaded from the ‘bipartite’ package (Poisot et al. 2011) in the R environment. The degree of divergence/convergence of sequentially assembled plant communities (stages a to f) was our measure of community dissimilarity. We calculated community dissimilarity for all possible pairs of mounds within each age group as well as all possible pairs of a mound and its reference community with two beta-diversity indices. The Jaccard dissimilarity index is based on presence/absence data only, while Bray–Curtis distance index uses quantitative relative coverage values for each species. Both indices are widely applied into community studies because they are good at evaluating differences among communities (Faith et al. 1987). Distance analyses were computed using the “Vegan” package in R. RESULTS Species richness and plant cover Many plant species were recorded in the study area: 90 species belonging to 26 families were observed in the 80 mounds + 18 background vegetation quadrats (see online supplementary Table S1). Of these, 52 species were perennial herbs, 10 biennial herbs, 26 annual herbs and 2 shrubs. Most species belonged to nine plant families: Cyperaceae (3 species), Poaceae (7 species), Asteraceae (18 species), Fabaceae (6 species), Lamiaceae (8 species), Polygonaceae (6 species), Ranunculaceae (5 species), Gentianaceae (5 species) and Rosaceae (7 species). In 2009, species richness on mounds that we labeled in our study area ranged from 2 to 29, and total plant cover from 1% to 99%. The changes trend of total plant cover on a mound generally followed patterns of species richness within each plot (Fig. 2). Both species richness and plant cover of mound communities were low at stage a but then increased greatly within each plot (Fig. 2). Species richness and plant cover of mound communities in plots B and C increased continuously over time, but that in plot A increased and then leveled off (Fig. 2, see online supplementary Table S2). From stage b onwards, species richness of mound communities was higher in plot C than in the other two plots (Fig. 2). The total cover of mound communities was higher in plot A than the other two plots from stage b to d (Fig. 2). Species richness and total plant cover of the background community were higher than that of any stage in the mound communities. Figure 2: View largeDownload slide changes in species richness (above) and total cover (below) of mounds/quadrats during succession within three background communities (A, B and C). The dotted reference line on the y-axis was the mean of species richness at stage b (above) and cover at stage c (below) in plot A. a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. Figure 2: View largeDownload slide changes in species richness (above) and total cover (below) of mounds/quadrats during succession within three background communities (A, B and C). The dotted reference line on the y-axis was the mean of species richness at stage b (above) and cover at stage c (below) in plot A. a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. Species composition of the different life forms showed different patterns over time (Figs. 3 and 4, see online supplementary Table S2). Early stages of mound communities had high relative cover of annual forbs with low species richness. Over time, relative cover and species richness of mound communities increased for perennial forbs and legumes. In addition, species richness but relative cover of biennial forbs and gramineous plants of mound communities increased with succession stage. In background community, relative cover of perennial groups and species richness of annual forbs were lower than older mound communities. Figure 3: View largeDownload slide changes in species richness of different life-forms during succession on Zokor mounds in three background communities (A, B and C). a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. Figure 3: View largeDownload slide changes in species richness of different life-forms during succession on Zokor mounds in three background communities (A, B and C). a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. Figure 4: View largeDownload slide changes in relative cover of life forms during succession on Zokor rodent mounds in three background communities (A, B and C). a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. Figure 4: View largeDownload slide changes in relative cover of life forms during succession on Zokor rodent mounds in three background communities (A, B and C). a–e are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and f is the background community. Lines show significant regressions over the succession period a–e. Parameters and statistics for these are in supplementary Table S2. TWINSPAN classification In the TWINSPAN cluster analysis, younger communities were the first to segregate and older communities were the last (Fig. 5A). Similarly, the species growing on fresh mounds clustered separately from those on older mounds and the background communities (Fig. 5B, see online supplementary Fig. S1 and Table S1). Early-successional species, such as Elsholtzia densa, Polygonum nepalense, Pedicularis chinensis, Chenopodium glaucum, Geranium sibiricum, Polygonum convolvulus and Hypecoum erectum, were abundant on open patches but distributed sporadically on older mounds. The later-successional group was classified into two sub-groups: (i) species found on newly established communities and abundant on older mounds, such as Artemisia tangutica, Potentilla anserina and Lancea tibetica and (ii) species that occurred mainly on older mounds and in the background communities, such as Elymus nutans. Figure 5: View largeDownload slide TWINSPAN analysis for 98 communities in the study site. A. Dendrogram of TWINSPAN classification. Community number of each group figure into textbox on each branch. B. Stepwise indicator species for the groups at each branch together with pseudospecies values. Underline parts indicate different groups. I: fresh mounds (younger than 1 year old); II: 2- to 4-year-old mounds; III: background communities. The sequence of branching is given with numbers. (+) denoting left arm preferential species, (-) denoting right arm preferential species. Figure 5: View largeDownload slide TWINSPAN analysis for 98 communities in the study site. A. Dendrogram of TWINSPAN classification. Community number of each group figure into textbox on each branch. B. Stepwise indicator species for the groups at each branch together with pseudospecies values. Underline parts indicate different groups. I: fresh mounds (younger than 1 year old); II: 2- to 4-year-old mounds; III: background communities. The sequence of branching is given with numbers. (+) denoting left arm preferential species, (-) denoting right arm preferential species. Detrended correspondence analysis DCA ordination diagrams of mound communities with increasing age illustrate relationships among plant communities during succession (Fig. 6). Community composition changed over time after mound formation. Almost all new mounds are scattered in ordination space, but with increasing age mounds showed higher aggregation in all plots (Fig. 6, see online supplementary Table S3). Over the sequence of successional stages, the DCA of species composition produced eigenvalues (λ) of 0.849, 0.618, 0.525, 0.437, 0.417 and 0.284, and gradient lengths (GL) of 5.122, 3.083, 3.261, 2.785, 2.488 and 1.851 SD units from axis I of each diagram. These community species scores exhibited notable changes from one succession stage to the next. In the first stage, the high-cover species, E. densa, was located close to the center of DCA ordination graph (Fig. 6A). Most species with high scores at the second stage belonged to seed bank species and bud bank colonizers, such as E. densa and P. anserina. As time progresses, the scores of annuals gradually decreased, while the scores of perennial plants increased, and several perennial species coexist in the community (Fig. 6B–E). For example, E. densa gradually decreased from stage c, while P. anserina gradually increased. Perennial grasses gradually dominated after stage d. Figure 6: View largeDownload slide ordination diagrams from detrended correspondence analyses of data on Zokor mound communities at several successional stages. A–E are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and F is the background community. The scientific names of the species are abbreviated (see online supplementary Table S1). Figure 6: View largeDownload slide ordination diagrams from detrended correspondence analyses of data on Zokor mound communities at several successional stages. A–E are 1 week to 3 month-, 1-, 2-, 3- and 4-year-old mounds, respectively, and F is the background community. The scientific names of the species are abbreviated (see online supplementary Table S1). Species co-occurrence patterns The C-scores for the consecutive communities showed clear and progressive changes with increasing age within each plot (Table 2). The C-scores for new mounds in plots B and C were significantly higher than random expectations. However, most C-scores of older mounds and ‘background’ communities had lower values than expected. Table 2: observed and expected [mean (95% confidence interval)] values for C-scores for co-occurrence analyses of successional stages on Zokor rodent mounds within three background communities C-score observed expected A a 6.7 6.32 (6.29, 6.35) b 19.4 24(23.9, 24.1) c 68.8 90.3 (90.04, 90.56) d 22.3 43.13 (43.01, 43.25) e 35.4 57.68 (57.49, 57.88) f 32.1 64.67 (64.5, 64.84) B a 13.8 9.63 (9.58, 9.67) b 36.1 38.38 (38.24, 38.52) c 26.2 48.31 (48.14, 48.47) d 56.1 91.96 (91.75, 92.17) e 97.2 130.28 (129.94, 130.63) f 50.8 125.5 (125.23, 125.78) C a 7.9 6.38 (6.35, 6.41) b 39.2 72.02 (71.85, 72.2) c 85.1 133.73 (133.38, 134.08) d 59.3 102.13 (101.9, 102.36) e 72.1 109.24 (108.94, 109.54) f 74.3 163.86 (163.53, 164.2) C-score observed expected A a 6.7 6.32 (6.29, 6.35) b 19.4 24(23.9, 24.1) c 68.8 90.3 (90.04, 90.56) d 22.3 43.13 (43.01, 43.25) e 35.4 57.68 (57.49, 57.88) f 32.1 64.67 (64.5, 64.84) B a 13.8 9.63 (9.58, 9.67) b 36.1 38.38 (38.24, 38.52) c 26.2 48.31 (48.14, 48.47) d 56.1 91.96 (91.75, 92.17) e 97.2 130.28 (129.94, 130.63) f 50.8 125.5 (125.23, 125.78) C a 7.9 6.38 (6.35, 6.41) b 39.2 72.02 (71.85, 72.2) c 85.1 133.73 (133.38, 134.08) d 59.3 102.13 (101.9, 102.36) e 72.1 109.24 (108.94, 109.54) f 74.3 163.86 (163.53, 164.2) A, B and C are the three areas, and a–f are 1 week to 3 months, 1-, 2-, 3-, 4-year-old mounds and the background communities, respectively. View Large Dissimilarity of mound communities Based on presence-absence and cover data for species, four levels of community convergence during succession in the mound communities were identified (Fig. 7): (i) convergence of communities at the same stage within each area, (ii) convergence in species composition of zokor mounds and their background community within each area, (iii) convergence of communities at the same stage across the three sample areas and (iv) convergence in species composition of zokor mounds and their background community across the three sample areas. Figure 7: View largeDownload slide changes in mean plant community dissimilarity among plant communities pairs of Zokor mound and/or background quadrat along 6 successional stages in three background communities, A, B and C. All refers to the meta-community of all three background communities. SSJ (black dotted line): Jaccard dissimilarity between community pairs at the same stage; SSB (black solid line): Bray–Curtis dissimilarity between community pairs at the same stage; SBJ (grey dotted line): Jaccard dissimilarity between mound community and their background community; SBB (grey solid line): Bray–Curtis dissimilarity index between mound community and their background community. Mean ± SE, n = 6~153 pairs of communities. Figure 7: View largeDownload slide changes in mean plant community dissimilarity among plant communities pairs of Zokor mound and/or background quadrat along 6 successional stages in three background communities, A, B and C. All refers to the meta-community of all three background communities. SSJ (black dotted line): Jaccard dissimilarity between community pairs at the same stage; SSB (black solid line): Bray–Curtis dissimilarity between community pairs at the same stage; SBJ (grey dotted line): Jaccard dissimilarity between mound community and their background community; SBB (grey solid line): Bray–Curtis dissimilarity index between mound community and their background community. Mean ± SE, n = 6~153 pairs of communities. DISCUSSION The results from TWINSPAN classification and life-form groups analysis showed that recently disturbed mounds were closely associated with annual forbs (e.g. E. densa, P. convolvulus, C. glaucum, H. erectum, P. chinensis, P. nepalense) or perennial forbs (e.g. G. sibiricum, L. tibetica, A. tangutica, P. anserine). These early-arriving and disturbance-sensitive colonists with the characteristics of r-strategists are short-lived annuals emerging from the seed bank or longer-lived perennials sporadically resprouting from root fragments in the soil of fresh mounds, or dispersing from adjacent habitats, establishing high-abundance populations with rapid growth rates because they can benefit from the resource-rich patches immediately after the disturbance (Houseman et al. 2008; Latzel et al. 2008). It was in line with some results from a seed addition experiment, which indicated that, in the case of unrestricted dispersal, superior colonizers profit from non-equilibrium conditions caused by small disturbances (Inouye and Tilman 1995; Lanta and Lepš 2009). Also, it provided evidence to Sheil and Burslem (2003)’s view, which pointed out that the significant increase in species richness following gap creation has been shown to arise primarily from regrowth of existing propagules and immigration. However, patterns at the life-form level showed trends different for early-colonizing annual versus perennial forbs. We found evidence that annuals are best at colonizing open habitats but have poor competitive ability, and are therefore not able to compete against more vigorous perennial species when the latter arrive, while perennial forbs have an advantage in competition (Jones et al. 2008). This could account for the continuous increase in species richness and coverage of plant communities in the later period. This result was consistent with many studies on succession, in which pioneer species are replaced by more competitive taxa (Cutler 2010; Lanta and Lepš 2009; Lebrija-Trejos et al. 2010). However, it contrasts with the pattern seen in three meadows in the western Cascade Range of Oregon, where few annuals and little evidence of competitive exclusion of subordinate species appeared on gopher mounds (Jones et al. 2008). Several factors may help to explain why competitive exclusion occurred on zokor mounds in our meadow, including a competition–colonization tradeoff, sensitivity of different life-forms and the presence of a seed bank and/or below-ground organs (Cadotte 2007; Jones et al. 2008; Questad and Foster 2007). For example, early-appearing annual forbs (such as E. densa) declined before the slower growing superior competitors (such as P. anserina) became dominant or reached their maximum abundances, probably in part because the crust formed on the mound surface as time passes prevented seedlings from germinating (Collins and Smith 2006). Three observations support the conclusion that disturbed systems can contain species with both good colonization and strong competitive ability (Limberger and Wickham 2012). First, E. densa, an annual forb with high dispersal ability, dominates on younger mounds but persists in older mounds and in background communities. The relevant mechanism may be a mass effect, which helps regulate and maintain annual forbs appearing in mound and non-mound communities (Chase et al. 2005). This is consistent with previous reports in which ephemeral plants with wind-dispersed seeds and fast growth can grow in unoccupied habitats in a wide variety of locations (Questad and Foster 2007; Schiffers et al. 2010), but contrary to the study by Jones et al. (2008), in which forbs were virtually absent from communities with little gopher disturbance. Second, our results also show that some perennial forbs that are good at colonizing empty mounds also dominated later, outcompeting other species. This could be explained by regenerative strategies by which perennial rhizomatous forbs with rhizomes, roots and tubers recolonize empty habitats (Latzel et al. 2008, 2009). Third, we also observe that the perennial grass L. secalinus regenerated sporadically on young mounds in plot B. This could be because (i) L. secalinus was the dominant species in the background community and was not limited by dispersal there, (ii) the soil is loose in newly formed mounds or (iii) some buried individuals of L. secalinus regenerated. Taken together, our observations imply either (i) that there is not always a tradeoff between these two abilities, or (ii) that dispersal does not always limit colonization by more competitive species, even if they lack long-distance dispersal ability. The C-scores within the community-level matrix increased with successional stage, implying that species interactions become more obvious over time, similar to the results of Cutler (2010). The negative co-occurrence patterns in the first stage may occur due to low population densities, resulting in unoccupied niches and weak interspecific interactions. In contrast, positive co-occurrence patterns were formed in later stages, indicating greater aggregation of species than expected by chance due to a change in the balance between competition and localized dispersal. In our study, two different diversity patterns were observed during post-disturbance succession on zokor mounds. Species richness and plant cover increased continuously in plots B and C. This is consistent with previous findings in other small-scale disturbed grassland ecosystems (Cutler 2010; Jones et al. 2008; Wang et al. 2008). Unlike plots B and C, however, species richness and plant cover in plot A quickly increased to a peak, followed by a slight subsequent decline. This could be due to a limited species pool in area A where the total species richness was much lower than in the other two areas. Our study showed mound communities accumulate species immediately following disturbance: both species richness and cover increased quickly, consistent with results of Wang et al. (2008). High species accumulation rates without competition exclusion could account for this observation (e.g. Hubbell et al. 1999; Tilman 1994). On the other hand, species accumulation on empty patches during secondary succession can also be slow (Cutler 2010; Sheil and Burslem 2003). After individuals fill the mound community, it is no longer in a ‘gap’ state. It has been hypothesized that if mound communities provided intermittent empty niches for colonization, diversity might rise continuously (Debra et al. 2014; Sheil and Burslem 2003). Alternatively, species diversity can decrease through competitive exclusion (Chesson 2000; Limberger and Wickham 2012). Thus, we conclude that two patterns in our study may reflect transient changes in species accumulation rates. None of our results support the hypothesis of historical contingency and climatic determinism (Clements 1916; Egler 1954). Initial mound communities are highly variable (Lepš and Rejmánek 1991; Schiffers et al. 2010). Only a small number of individuals occurred on the mounds in the first few days, and the identities of the first colonists to arrive on newly formed mound communities are highly variable, so stochastic factors may dominate. Alternatively or in addition, small-scale disturbance can lead to fluctuations in micro-environmental conditions, resulting in colonization by fugitive species (Collins and Smith 2006; Questad and Foster 2007). In our study area, both annual and perennial forbs have the chance to occur on new zokor mounds, but species richness was low on each mound, resulting in communities with negative species co-occurrence patterns at this stage. This result contrasts with recent study in California grasslands (Larios et al. 2013), in which observed annual exotic dominant species appeared in the first year after an extreme disturbance. Zokor mound communities at the same stage became more similar to each other and to their background communities dominated by gramineous species over time. Previous studies of long-term primary succession on lava flows (Cutler 2010) and secondary succession in tropical forests (Lebrija-Trejos et al. 2010) have produced similar results. We also observe a decline in eigenvalues and GLs of DCA ordination with increasing age, suggesting that community structure of all the quadrats tended to converge over time. Such community patterns resulted from the combination of a sharp decline of few dominant pioneer species (annual forbs) and a gradual increase of species from the background community (perennial forbs) from a limited species pool. This is not consistent with the historical contingency hypothesis, which holds that community is determined by the order and timing that individual species entered the system, even if environmental conditions and species pool are similar. Our results suggest that regeneration from broken roots and germination from the seed bank drive early succession, but background-based species pool and competition exclusion become more important over time. Our results provide one of the first examples of zokor mound community convergence over time, whereas mound communities among the three plots at the same stage apparently diverged later in succession. In other words, older mound communities in different background communities diverged due to background-based assembly, while mounds in the same ‘background community’ and even the whole study area converged. This is partially consistent with results from a 5-year container experiment in meadow and floodplain communities, in which final composition converged in each site but diverged between the two habitats (Lanta and Lepš 2009). Different background vegetation resulting from three restoration treatments is one of the potential explanations for our result. As is most studies of succession, the absence of replication at the background-community level limits the strength of our inferences. Future studies should use randomized blocks for background-community treatment, and include a broader spectrum of grassland types. Also, our study has focused on vegetation colonization within the first four years after disturbance in a sub-alpine meadow, but longer term studies are needed. Ongoing and future studies on zokor disturbance communities can help us test theories of succession, as well as helping us to develop strategies for the management of grassland ecosystems. CONCLUSIONS In our 4-year field experiment on plant succession on zokor mounds, we observed two different patterns of plant community diversity over time: (i) continuously increasing and (ii) an initial increase followed by a later decline. The relative cover changes in life forms do not correspond with patterns in species richness. Succession of zokor mound communities was initially consistent with ‘multiple alternative states’, whereas later succession moved towards ‘non-random community structuring by background-based assembly’. Species composition on newly formed mound communities was primarily determined by species’ dispersal colonization abilities and regenerative strategies, while species co-occurring later were primarily the outcome of species interactions. Tradeoffs between competition and colonization, characteristics of different life-forms and mass effects in a limited species pool are mechanisms that could explain the convergent development pattern under low disturbance frequency. Small-scale disturbances by rodents provide unique opportunities to investigate and test theories of secondary succession and increase our ability to manage the degraded grasslands. SUPPLEMENTARY MATERIAL Supplementary data are available at Journal of Plant Ecology online. FUNDING Fund for Zokor Experiment: National Natural Science Foundation of China (No. 31601984, 41101524), Fundamental Research Funds for the Central Universities (lzujbky-2012–108, lzujbky-2015–191) and the Special Fund for Agro-scientific Research in the Public Interest (No. 201203006). ACKNOWLEDGEMENTS We thank Xiaojuan Xin, Renyi Zhang, Rui Qi, Rui Hou, Zhiyong Yang, Yuanzhi Li, Peng Jia and Xiaoliang Qin for help during the study, and two anonymous reviewers for helpful comments on the manuscript. Conflict of interest statement. None declared. REFERENCES Beals SC Hartley LM Prevéy JSet al. ( 2014) The effects of black-tailed prairie dogs on plant communities within a complex urban landscape: an ecological surprise? Ecology 95: 1349– 59. Google Scholar CrossRef Search ADS PubMed Cadotte MW( 2007) Competition-colonization trade-offs and disturbance effects at multiple scales. Ecology 88: 823– 9. Google Scholar CrossRef Search ADS PubMed Chase JM ( 2003a) Community assembly: when should history matter? Oecologia 136: 489– 98. Google Scholar CrossRef Search ADS Chase JM ( 2003b) Experimental evidence for alternative stable equilibria in a benthic pond food web. Ecol Lett 6: 733– 41. Google Scholar CrossRef Search ADS Chase JM( 2010) Stochastic community assembly causes higher biodiversity in more productive environments. Science 328: 1388– 91. Google Scholar CrossRef Search ADS PubMed Chase JM Amarasekare P Cottenie Ket al. ( 2005) Competing theories for competitive metacommunities. In Holyoak M Leibold MA Holt RD (eds). Metacommunities: Spatial Dynamics and Ecological Communities . Chicago/London: University of Chicago Press, 335– 54. Chesson P( 2000) Mechanisms of maintenance of species diversity. Annu Rev Ecol Syst 31: 343– 66. Google Scholar CrossRef Search ADS Clements FE( 1916) Plant succession: an analysis of the development of vegetation. Estuar Coast Shelf Sci 73: 223– 35. Collins SL Smith MD( 2006) Scale-dependent interaction of fire and grazing on community heterogeneity in tallgrass prairie. Ecology 87: 2058– 67. Google Scholar CrossRef Search ADS PubMed Cutler N( 2010) Long-term primary succession: a comparison of non-spatial and spatially explicit inferential techniques. Plant Ecol 208: 123– 36. Google Scholar CrossRef Search ADS Debra ZD Bernhard S Petermann JSet al. ( 2014) Selection for niche differentiation in plant communities increases biodiversity effects. Nature 515: 108– 11. Google Scholar CrossRef Search ADS PubMed Diamond JM( 1975) Assembly of species communities. Ecol Evol Commun 342: 444. Egler FE( 1954) Vegetation science concepts I. Initial floristic composition, a factor in old-field vegetation development with 2 figs. Vegetatio 4: 412– 7. Google Scholar CrossRef Search ADS Erfanzadeh R Pétillon J Maelfait J-Pet al. ( 2010) Environmental determinism versus biotic stochasticity in the appearance of plant species in salt-marsh succession. Plant Ecol Evol 143: 43– 50. Google Scholar CrossRef Search ADS Faith DP Minchin PR Belbin L( 1987) Compositional dissimilarity as a robust measure of ecological distance. Vegetatio 69: 57– 68. Google Scholar CrossRef Search ADS Fridley JD( 2013) Successional convergence, stochastic assembly and the future of tropical forests. J Veg Sci 24: 415– 6. Google Scholar CrossRef Search ADS Fukami T Lee W( 2006) Alternative stable states, trait dispersion and ecological restoration. Oikos 113: 353– 6. Google Scholar CrossRef Search ADS Fukami T Martijn Bezemer T Mortimer SRet al. ( 2005) Species divergence and trait convergence in experimental plant community assembly. Ecol Lett 8: 1283– 90. Google Scholar CrossRef Search ADS Fukami T( 2010) Community assembly dynamics in space. In Verhoef HA Morin PJ (eds). Community Ecology: Processes, Models, and Applications . Oxford: Oxford University Press, 45– 54. Google Scholar CrossRef Search ADS Gotelli NJ McCabe DJ( 2002) Species co-occurrence: a meta-analysis of JM Diamond’s assembly rules model. Ecology 83: 2091– 6. Google Scholar CrossRef Search ADS Hill MO Šmilauer P( 2005) TWINSPAN for Windows. Version 2.3 . Huntingdon & Ceské Budejovice: Centre for Ecology and Hydrology, University of South Bohemia. Houseman GR Mittelbach GG Reynolds HLet al. ( 2008) Perturbations alter community convergence, divergence, and formation of multiple community states. Ecology 89: 2172– 80. Google Scholar CrossRef Search ADS PubMed Hubbell SP Foster RB O’Brien STet al. ( 1999) Light-gap disturbances, recruitment limitation, and tree diversity in a neotropical forest. Science 283: 554– 7. Google Scholar CrossRef Search ADS PubMed Inouye RS Tilman D( 1995) Convergence and divergence of old-field vegetation after 11 yr of nitrogen addition. Ecology 76: 1872– 87. Google Scholar CrossRef Search ADS Jones CC Halpern CB Niederer J( 2008) Plant succession on gopher mounds in western Cascade meadows: consequences for species diversity and heterogeneity. Am Midl Nat 159: 275– 86. Google Scholar CrossRef Search ADS Kang L Han X Zhang Zet al. ( 2007) Grassland ecosystems in China: review of current knowledge and research advancement. Phil Trans R Soc B Biol Sci 362: 997– 1008. Google Scholar CrossRef Search ADS Lanta V Lepš J( 2009) How does surrounding vegetation affect the course of succession: a five‐year container experiment. J Veg Sci 20: 686– 94. Google Scholar CrossRef Search ADS Larios L Aicher RJ Suding KN( 2013) Effect of propagule pressure on recovery of a California grassland after an extreme disturbance. J Veg Sci 24: 1043– 52. Google Scholar CrossRef Search ADS Latzel V Hájek T Klimešová Jet al. ( 2009) Nutrients and disturbance history in two Plantago species: maternal effects as a clue for observed dichotomy between resprouting and seeding strategies. Oikos 118: 1669– 78. Google Scholar CrossRef Search ADS Latzel V Mihulka S Klimešová J( 2008) Plant traits and regeneration of urban plant communities after disturbance: does the bud bank play any role? Appl Veg Sci 11: 387– 94. Google Scholar CrossRef Search ADS Lebrija-Trejos E Pérez-García EA Meave JAet al. ( 2010) Functional traits and environmental filtering drive community assembly in a species-rich tropical system. Ecology 91: 386– 98. Google Scholar CrossRef Search ADS PubMed Legendre P Legendre L( 2012) Numerical Ecology . Oxford: Elsevier. Lepš J Rejmánek M( 1991) Convergence or divergence: what should we expect from vegetation succession? Oikos 62: 261– 4. Google Scholar CrossRef Search ADS Limberger R Wickham SA( 2012) Disturbance and diversity at two spatial scales. Oecologia 168: 785– 95. Google Scholar CrossRef Search ADS PubMed do Nascimento LM Sampaio EVdSB Rodal MJNet al. ( 2014) Secondary succession in a fragmented Atlantic Forest landscape: evidence of structural and diversity convergence along a chronosequence. J For Res 19: 1– 13. Google Scholar CrossRef Search ADS Oksanen J Minchin PR( 1997) Instability of ordination results under changes in input data order: explanations and remedies. J Veg Sci 8: 447– 54. Google Scholar CrossRef Search ADS Oksanen J Blanchet F Kindt Ret al. ( 2011) vegan: Community Ecology Package. R package version 2.0-2. http://CRAN.R-project.org/package=vegan. Poisot T Bever JD Nemri Aet al. ( 2011) A conceptual framework for the evolution of ecological specialisation. Ecol Lett 14: 841– 51. Google Scholar CrossRef Search ADS PubMed Questad EJ Foster BL( 2007) Vole disturbances and plant diversity in a grassland metacommunity. Oecologia 153: 341– 51. Google Scholar CrossRef Search ADS PubMed Saccone P Pyykkonen T Eskelinen Aet al. ( 2014) Environmental perturbation, grazing pressure and soil wetness jointly drive mountain tundra toward divergent alternative states. J Ecol 102: 1661– 72. Google Scholar CrossRef Search ADS Samuels CL Drake JA( 1997) Divergent perspectives on community convergence. Trends Ecol Evol 12: 427– 32. Google Scholar CrossRef Search ADS PubMed Schiffers K Tielbörger K Jeltsch F( 2010) Changing importance of environmental factors driving secondary succession on molehills. J Veg Sci 21: 500– 6. Google Scholar CrossRef Search ADS Sheil D Burslem DFRP( 2003) Disturbing hypotheses in tropical forests. Trends Ecol Evol 18: 18– 26. Google Scholar CrossRef Search ADS Suding KN Gross KL Houseman GR( 2004) Alternative states and positive feedbacks in restoration ecology. Trends Ecol Evol 19: 46– 53. Google Scholar CrossRef Search ADS PubMed Suding KN Hobbs RJ( 2009) Threshold models in restoration and conservation: a developing framework. Trends Ecol Evol 24: 271– 9. Google Scholar CrossRef Search ADS PubMed Tilman D( 1994) Competition and biodiversity in spatially structured habitats. Ecol 75: 2– 16. Google Scholar CrossRef Search ADS Wang TC Xiong YC Ge JPet al. ( 2008) Four-year dynamic of vegetation on mounds created by zokors (Myospalax baileyi) in a subalpine meadow of the Qinghai-Tibet Plateau. J Arid Environ 72: 84– 96. Google Scholar CrossRef Search ADS Yang Z Guo H Zhang Jet al. ( 2013) Stochastic and deterministic processes together determine alpine meadow plant community composition on the Tibetan Plateau. Oecologia 171: 495– 504. Google Scholar CrossRef Search ADS PubMed Young TP Chase JM Huddleston RT( 2001) Community succession and assembly comparing, contrasting and combining paradigms in the context of ecological restoration. Ecol Restor 19: 5– 18. Google Scholar CrossRef Search ADS Zhang YM Liu JK( 2003) Effects of plateau zokors (Myospalax fontanierii) on plant community and soil in an alpine meadow. J Mammal 84: 644– 51. Google Scholar CrossRef Search ADS © The Author(s) 2017. Published by Oxford University Press on behalf of the Institute of Botany, Chinese Academy of Sciences and the Botanical Society of China. All rights reserved. For permissions, please email: email@example.com
Journal of Plant Ecology – Oxford University Press
Published: Jun 1, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
All the latest content is available, no embargo periods.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud