Chuanwu Chen, Marcel Holyoak, Yanping Wang, Xingfeng Si, Ping Ding. 2019: Spatiotemporal distribution of seasonal bird assemblages on land-bridge islands: linking dynamic and static views of metacommunities. Avian Research, 10(1): 25. DOI: 10.1186/s40657-019-0164-7
Citation: Chuanwu Chen, Marcel Holyoak, Yanping Wang, Xingfeng Si, Ping Ding. 2019: Spatiotemporal distribution of seasonal bird assemblages on land-bridge islands: linking dynamic and static views of metacommunities. Avian Research, 10(1): 25. DOI: 10.1186/s40657-019-0164-7

Spatiotemporal distribution of seasonal bird assemblages on land-bridge islands: linking dynamic and static views of metacommunities

More Information
  • Corresponding author:

    Ping Ding, dingping@zju.edu.cn

  • Received Date: 21 Mar 2019
  • Accepted Date: 13 Jun 2019
  • Available Online: 24 Apr 2022
  • Publish Date: 19 Jun 2019
  • Background 

    Although assessing temporal dynamics of populations is crucial for understanding metacommunities, empirical studies have primarily analyzed only static snapshots of communities. Here, we present a holistic view of how species traits and habitat characteristics relate to metacommunity dynamics and use it to test for differences in the spatiotemporal distribution of seasonal bird assemblages.

    Methods 

    We surveyed forest birds in breeding and winter seasons within 36 islands for 9 years. We then grouped birds into four landbird assemblages, selected on the basis of published differences in biology or ecology: winter residents, migratory winter visitors, breeding summer residents, and migratory summer visitors. We estimated dynamic species colonization and extirpation through the 9-year period, and evaluated the associations among island attributes, species attributes and community composition.

    Results 

    Overall, winter and summer residents showed strong associations between composition and habitat structure of the islands. In addition, winter and summer residents on large islands had lower extirpation and turnover than winter and summer visitors. Visitor assemblages showed no significant habitat associations, and in winter had high extirpation rates and small body sizes. By contrast, local extirpation of summer visitors was correlated with local species richness, indicating a likely effect of competition on extirpation.

    Conclusions 

    Our results demonstrated repeated patterns among species composition, bird traits, habitat/island characteristics and observed metacommunity dynamics. Winter and summer residents best matched species sorting and patch dynamics, respectively, due to differences in resource availability and requirements of overwinter survival versus breeding. Summer visitors were consistent with species sorting and winter visitors were randomly distributed, likely because of interactions with resident competitors. Our results highlight that coexisting seasonal migrant and resident assemblages differ in their spatial dynamics, with consequences for relevant conservation and management strategies.

  • Understanding how species composition varies among sites across space or time (i.e., beta diversity) is a key issue and fundamental topic in ecology (Svenning et al., 2011; da Silva et al., 2018; Carlos-Júnior et al., 2019; Zeng et al., 2022). The increased interest in beta diversity comes from the recognition of its important roles in revealing community assembly mechanisms (Zellweger et al., 2017; Soininen et al., 2018; Du et al., 2021), and helping guide conservation practices (reviewed in Socolar et al., 2016). Major progresses have been made in the field of beta diversity studies by decomposing beta diversity into its turnover and nestedness components (Baselga, 2010), and/or by incorporating species' functional traits and evolutionary histories (Swenson et al., 2011; Leprieur et al., 2012). Specifically, the former enables deeper insights into mechanisms driving the variation in biodiversity (Baselga et al., 2012; Si et al., 2015; Soininen et al., 2018), while the latter offers a complete perspective on biodiversity and better captures community assembly (Perez Rocha et al., 2018; Branco et al., 2020; González-Trujillo et al., 2020; Jiang et al., 2021). Thus, it has now been widely acknowledged that simultaneous consideration of taxonomic, functional, and phylogenetic dimensions of biodiversity, and adoption of beta diversity partitioning framework would help us better understand community structure and the underlying mechanism of community ecological processes.

    Beta diversity was found to be driven by niche-based deterministic or dispersal-based stochastic processes (Chase et al., 2011; Swenson et al., 2011; Segre et al., 2014; Du et al., 2021). These two processes are recognized as mutually inclusive and can co-occur simultaneously within communities or sequentially along environmental gradients (Soininen et al., 2018; Wang et al., 2020; Du et al., 2021). Assessing the relative contribution of deterministic and stochastic processes in driving community dissimilarities remains a central challenge and the research frontier in community ecology (Siefert et al., 2013; Vellend et al., 2014; Du et al., 2021). Among the many methods used in this effort, patterns of distance decay, the decrease in community similarity with increasing geographical distance, have been regarded as a primary tool for inferring community assembly mechanisms (Nekola and White, 1999; Nekola and McGill, 2014; Soininen et al., 2007; Steinbauer et al., 2012; Morlon et al., 2008; Siefert et al., 2013; Graco-Roza et al., 2022; Martínez-Santalla et al., 2022). Because both climatic constraint and dispersal limitation are the main drivers of distance-decay of similarity (Baselga and Gómez-Rodríguez, 2021), several studies have been conducted to compare the relative roles of climatic and spatial distances to infer which processes played a dominant role in shaping community assemblages. Specifically, higher correlations of assemblage similarity with spatial than with climatic distance indicate dominant role of dispersal limitation, while species sorting across climatic gradients via niche processes dominate community assemblages otherwise (Condit et al., 2002; Tuomisto et al., 2003; Baselga and Gómez-Rodríguez, 2021).

    Another promising approach which measures traits or lineage composition varying among sites (i.e., functional or phylogenetic beta diversity) could provide additional insights into community assembly (Swenson et al., 2011; Leprieur et al., 2012; Siefert et al., 2013). Specifically, if variation in trait or lineage composition is similar as random expectations given the observed species beta diversity, a product of neutral processes could be drawn. In contrast, when variation in trait or lineage composition was higher or lower than random expectations, niche-based deterministic processes were considered dominant in driving community dissimilarity (Swenson et al., 2011; Vellend et al., 2014; Du et al., 2021; Ge et al., 2021). However, integrating both above methods jointly remains scarce, which may hamper the opportunity of gaining new perspectives on community structure and assembly of communities.

    Mountains only occupy about 25% of the world's total continental land surface, but harbor roughly 87% of the world's species, and are thus believed to be critical for global biodiversity conservation (reviewed in Rahbek et al., 2019). Indeed, mountains are considered to provide the most powerful natural experiment for exploring spatial patterns of biodiversity and community structure as they contain dramatic geophysical changes over relatively short spatial distances (Körner, 2007). Yet since ecologists committed to mountain studies over two centuries (Rahbek et al., 2019), much less attention have been paid to elevational beta-diversity, especially multidimensional beta diversity, compared to the well-known elevational alpha-diversity (Wang et al., 2012, 2022; Mori et al., 2018; García-Navas et al., 2020). This lack of attention may hamper our understanding on mechanisms driving community dissimilarities (Wang et al., 2012; García-Navas et al., 2020, 2022), thereby limiting the opportunities for informing conservation practices and improving habitat management and restoration (Sreekar et al., 2020; Du et al., 2021).

    The Himalayan Mountains are the highest mountains in the world, containing a typical mountain environment that supports high biodiversity with many rare species. They are also popularly known for their high beta diversity (Gaston et al., 2007). This system provides unique opportunity to improve our understanding of niche-and dispersal-based processes, and then enrich knowledge bases of beta diversity studies. Here, we sampled bird communities along a 3600-m altitudinal gradient at Gyriong Valley of Gyirong County in Xizang, China, and simultaneously consider distance-decay models and multiple dimensions of beta diversity to address the following two questions: 1) What are the elevational patterns of taxonomic, functional, and phylogenetic bird beta diversity, and how do they differ from one another? 2) What are the taxonomic, functional, and phylogenetic structures of bird communities, and the relative importance of niche- and neutral-based processes underlying them? According to the distance-decay relationship (Morlon et al., 2008), we predicted an increase in beta diversity with increasing elevational distance among communities. Further, we expected functional and phylogenetic beta diversity should be lower than their taxonomic counterparts because stressful environments at high elevations select for certain traits/lineages that are suited to those conditions, resulting in functional/phylogenetic clustering between elevation bands. On the other hand, when species with different ecological performance (e.g., unique traits or lineages) frequently loss or gain, functional/phylogenetic beta diversity could be higher than its taxonomic counterparts (Siefert et al., 2013; Weinstein et al., 2014; Du et al., 2021). In addition, we expected multiple dimensions of beta diversity should be driven largely by turnover in mountain birds because species are replaced frequently across various ecosystems and taxonomic groups (Soininen et al., 2018). Finally, strong competition could be induced as the environment changes from benign to harsh when moving from low to high elevation. We then expect stronger competitive exclusion and a weaker environmental filtering in structuring the communities with increasing elevational distances.

    Field studies were conducted in Gyirong Valley (28°15′–29°0′ N, 85°6′–85°41′ E), the longest and westernmost canyon of Mt. Qomolangma National Nature Reserve on the southern slope of the central Himalaya, China. The Himalayas is the youngest and highest mountain chain on the Earth, and is also known as 'Asia's water tower'. The studied valley has a length of ca. 72 km from the lowest Resuo Village at 1800 m a.s.l. to the summit of Mt. Mala at 5400 m a.s.l. This valley includes five distinct vegetation types, which are evergreen broadleaf forests (1800–2500 m), coniferous and broadleaf mixed forests (2500–3300 m), subalpine coniferous forests (3300–3900 m), alpine bush and meadow (3900–4700 m), and alpine tundra with sparse herbs (4700–5400 m) (Ding et al., 2021).

    We collected bird data (richness and abundance) using line transect methods (Bibby et al., 2000) in the wet season from 2012 to 2013 (i.e., four repeated surveys, May and June in 2012, August in 2012, September and October in 2012, July and August in 2013). This dataset was published in Ding et al. (2021) and we here re-analyzed this dataset for beta diversity measures. More detailed descriptions about bird survey methods are given in Ding et al. (2021), which focused on alpha components of taxonomic, functional, and phylogenetic diversity in the same region.

    Following previous studies (McCain, 2009; Wu et al., 2013; Quintero and Jetz, 2018; Ding et al., 2021), we excluded migrant birds which likely cause potential biases from seasonal, long-distance movements, but instead only included breeding resident birds for subsequent analyses. Information of migratory status was extracted from the local literature (Zheng et al., 1983). In addition, we also excluded shorebirds and owls because of their highly specific habitats, and nocturnal behavior, respectively. As a result, a total of 151 breeding birds were recorded in the present study (Appendix Fig. S1).

    We used four traits in this study, including one continuous (body mass) and three categorical types (feeding guild, foraging behavior, or stratum). Presence or absence was scored for feeding guild (carnivore, granivore, insectivore, or omnivore), foraging behavior (glean, probe, leap, or sally), and foraging stratum (ground, understorey, midstorey, canopy, or air) (Ding et al., 2013). All these selected traits were related to resource use, energy requirements and trophic levels (Ding et al., 2013; Si et al., 2016), which were compiled from the local literature The Avifauna of Xizang (Zheng et al., 1983). Note that all traits used in this study exhibited significant phylogenetic signal (Ding et al., 2021).

    We constructed a phylogenetic tree for the study species by pruning the global phylogenetic tree of birds from BirdTree (http://birdtree.org) under the option of 'Hackett All Species: a set of 10,000 trees with 9993 OTUs each' (Jetz et al., 2012). We then sampled 5000 pseudo-posterior distributions and estimated the Maximum Clade Credibility tree using mean node heights by the software TreeAnnonator v1.10.4 of the BEAST package (Drummond and Rambaut, 2007). The resulting tree was used for all subsequent analyses on phylogenetic analyses (Appendix Fig. S1).

    We computed multiple-site and pairwise taxonomic, functional and phylogenetic beta diversity following Baselga's method (Baselga, 2010). For multiple-site (indicated in capital letters), overall dissimilarity (βSOR) was partitioned into the pure spatial turnover (βSIM) and nestedness components (βSNE). The ratio between βSIM and βSOR (i.e., βRATIO = βSIM/βSOR) indicated the relative contribution of the turnover component to β-diversity from taxonomic, functional, and phylogenetic dimensions (i.e., TβRATIO, FβRATIO and PβRATIO). For pairwise site (indicated in small letters), the overall β diversity (βsor, βfuncsor and βphylosor refer to taxonomic, functional, and phylogenetic beta diversity, respectively), which was indicated by the Sørensen dissimilarity index, was additively decomposed into two components accounting for pure spatial turnover (βsim, βfuncsim and βphylosim) and nestedness (βsne, βfuncsne and βphylosne). All the above analyses were performed in R using 'betapart' package (Baselga and Orme, 2012).

    We used null models to assess whether taxonomic, functional, phylogenetic dissimilarities and its partitioned components were deviated from expected values by chance. We run null models by maintaining species occupancy and species richness 1000 times (i.e., fixed-row fixed-column algorithm, or quasiswap algorithm; Jonsson, 2001; Miklós and Podani, 2004), and calculated the standardized effect size of functional, phylogenetic dissimilarities and its partitioned components as:

    SES=ObsMeannullSDnull

    where Obs is the observed value of β-diversity measures, Meannull and SDnull is the mean and standard deviation of the null expectations, respectively.

    SES values deviated from null expectations indicate the dominance of non-random assembly (deterministic processes; Webb et al., 2002). Specifically, if SES values were higher than 1.96, it indicates the observed beta diversity was significantly higher than expected, and competitive exclusion may drive community overdispersion; if SES values were below −1.96, it indicates that observed beta diversity was significantly lower than expected, and environmental filtering may drive community clustering (Fig. 1). In addition to the extent of the deviation from our null expectation, we also aimed to investigate directional trends in the SES values along the elevational gradient. To this end, we used the loess function in R to plot the relationship of each SES values along the elevation gradient (see Fig. 2).

    Figure  1.  Conceptual framework showing the analytical approach in this study. (A) Beta diversity have three complementary components including taxonomic, functional and phylogenetic dimensions; (B) Pairwise community dissimilarities and its turnover and nestedness components; (C) Distance decay of similarities in three dimensions, and its corresponding turnover and nestedness components; (D) The relative importance of niche-based and neutral community assembly mechanisms by applying null model procedures.
    Figure  2.  The multiple-site Sørensen dissimilarity (βSOR) and its spatial turnover (βSIM) and nestedness components (βSNE). Values outside the brackets are multiple-site Sørensen dissimilarity, and those inside the brackets represent the ratio of each component to total dissimilarity (βRATIO).
    Figure  3.  The dissimilarities driven by turnover, nestedness and their combinations at taxonomic, functional, and phylogenetic dimensions of breeding birds in the central Himalaya, China.

    We computed the similarity index, 1–βsor (of βsim and βsne), using the beta.pair function of the package 'betapart' (Baselga and Orme, 2012). In this study, distance decay of similarity was modelled by the negative exponential curve between the community similarity and elevational distances (Nekola and McGill, 2014; Graco-Roza et al., 2022). These non-linear fits were performed using decay.model of the package 'betapart' (Baselga and Orme, 2012), and specified the similarity index and elevational distances as the response and explanatory variables, respectively. Specifically, we conducted the distance decay model via generalized linear models (GLMs) following a quasi-binomial family with log-link (Millar et al., 2011). In this framework, the slope of the model represents the rate of similarity decay, and the r2 represents the strength of similarity-distance relationships. Climate data (i.e., precipitation and temperature) were obtained from WordClim database (http://www.worldclim.org, 1950–2000) at a 30 arc-s resolution. The value for each elevational band was then calculated as the average of all grid cells within it (Ding et al., 2021). Following Baselga and Gómez-Rodríguez (2021), we quantified climatic effects using mean annual temperature and annual precipitation. Given their high correlation between mean annual temperature and annual precipitation (r = 0.986*), we used mean annual temperature only in this study.

    For multiple-site beta diversity, the spatial turnover component was the larger component of total dissimilarity. Specifically, spatial turnover had the highest importance in functional beta diversity, followed by phylogenetic and taxonomic analogues. In contrast, nestedness-resultant component played a minor role in generating beta-diversity patterns (all βRATIO > 0.5).

    For pairwise beta diversity, taxonomic beta diversity was relatively higher than its functional and phylogenetic dimensions (βsor > βfuncsim/βphylosim). Moreover, taxonomic turnover was relatively higher than its functional and phylogenetic dimensions (βsim > βphylosim/βfuncsim), whereas nestedness-resultant dissimilarities of functional dimensions tended to be higher than its species and phylogenetic dimension (βfuncsne > βphylosne/βsne).

    Considering absolute value, spatial distance decays of taxonomic similarities were the strongest, followed by the distance decays of phylogenetic similarities and functional similarities (Fig. 4). Similar results were found for turnover components of multidimensional beta diversity, whereas contrasting patterns were observed in its nestedness components. Specifically, slopes of the spatial distance decay were largely negative, and ranged from −0.00095 to 0.000034. Furthermore, we found that all correlations of assemblage similarity with climatic distance were higher than that with spatial distances, and these slopes ranged from −0.18 to 0.0065 (Appendix Fig. S2).

    Figure  4.  Distance-decay models for breeding birds along the central Himalayan elevation gradient. The negative values were highlighted in red. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

    Generally, the standardized effect size of overall taxonomic, functional, and phylogenetic beta diversity and their turnover components increased with increasing elevational distance (Fig. 5), indicating a general tendency that SES values vary from less than 0 to greater than 0. In contrast, the SES values of taxonomic and phylogenetic nestedness decreased with increasing elevational distance.

    Figure  5.  The relationships between the standardized effect size (SES) of turnover, nestedness and their combinations at taxonomic, functional, phylogenetic dimensions and elevational distance. The pink shaded area indicated the α = 0.05 threshold of ±1.96 SES for significantly non-random values. The blue line represents the fit based on Loess smoothing. Each circle represents a comparison between the lowest elevation (1800 m) and the subsequent higher elevations. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

    In this study, we analyzed taxonomic, functional, and phylogenetic beta diversity and its turnover and nestedness components of breeding birds along a 3600-m elevation gradient in the central Himalaya. We showed that spatial turnover components were the main drivers of multiple dimensions of beta diversity, and significant distance-decay relationships were found in all overall beta diversity and turnover components, and most of nestedness components. In addition, higher correlations of assemblage similarity with climatic than with spatial distance were detected. Furthermore, we found a central role for deterministic processes in shaping multiple dimensions of beta diversity of breeding birds in the central Himalayan elevational gradient.

    For multiple-site dissimilarities, turnover likely dominate multiple dimensions of beta diversity compared to nestedness across a wide range of taxa and realms (Soininen et al., 2018), indicating that turnover mainly measures compositional dissimilarities of communities between bands, whereas nestedness quantifies more about richness differences (e.g., Soininen et al., 2018; Zhao et al., 2021; Wang et al., 2022). This also suggests that species turnover caused by species sorting or dispersal processes quantifies mostly the same phenomenon as total dissimilarity (Baselga, 2010; Soininen et al., 2018).

    When considering pairwise beta diversity, the results of this study indicated that taxonomic turnover was higher than their functional and phylogenetic dimensions, whereas functional nestedness was higher than their taxonomic and phylogenetic analogues (Fig. 3). This is probably because species replacement could occur between phylogenetically related species with similar functional traits (Graham and Fine, 2008; Siefert et al., 2013; Villéger et al., 2013; Wang et al., 2022), but species loss/gain occurred more commonly among distantly related species with similar traits (Zhao et al., 2021). These results meant that different but interdependent outcomes could be produced by the assembling processes at different dimensions of beta diversity, highlighting the necessity of simultaneously integrating multiple dimensions when inferring assembly processes (Du et al., 2021; Wang et al., 2022).

    At the functional dimensions, we found an important role of nestedness in total dissimilarities. This result is consistent with previous work on ants (Du et al., 2021; Zhao et al., 2021), spiders (Wu et al., 2017), amphibians (Wang et al., 2022), passerines, and rodents, suggesting an ordered loss of birds' functional roles along elevational gradient to some extent (e.g., Bishop et al., 2015; Si et al., 2016; Du et al., 2021). In addition, variations in functional turnover reveal functional redundancy among bird communities (Ricotta et al., 2020), and this is also true for alpha functional diversity (Ding et al., 2021). In contrast, phylogenetic turnover remains dominating the total beta diversity, suggesting that turnover occurs in species with distant phylogenies (e.g., Du et al., 2021; Zhao et al., 2021). For instance, ground-feeding birds like quails and pheasants occupy a relatively distant lineage, but inhabit mainly at high elevations (Ding et al., 2019), and our previous findings showed that bird phylogenetic structure tended towards clustering at mid-elevations and overdispersion towards low and high elevations (Ding et al., 2021). Taken together, this might indicate that distantly related lineages occupy at mid-elevations, whereas closed related lineages occupy at both ends of gradient.

    However, we noted that our surveys spanned several months from May to early October, largely overlapping the breeding season in this region that are from May to the end of September (Zheng et al., 1983), so it may result in high proportion of turnover component than that if we had only collected bird data during a short period of the breeding season (Liang et al., 2021; He et al., 2022). In this region, some species (e.g., Chaimarrornis leucocephalus, Heterophasia capistrata) have the records of seasonal migration, the estimates of turnover component will be higher than that if we surveyed birds in the non-breeding seasons because of more constrained ranges of breeding birds. While such overestimation is certainly present, our results are still reliable, given the turnover is the dominant process of multiple dimensions of β diversity (Soininen et al., 2018). At the meantime, more studies are needed to test the seasonal movements of breeding birds in this area when the data are available.

    We found that elevational distances negatively influenced multiple dimensions of community similarity (except taxonomic nestedness), and distance decay of total taxonomic similarities was stronger than its functional and phylogenetic dimensions (Fig. 4). This result is consistent with previous global syntheses revealing that taxonomic distance decay is generally stronger than its functional dimensions (reviewed in Graco-Roza et al., 2022), likely reflecting the fact that species variation occurred in different traits/phylogeny among assemblages (Fig. 1). The major reason lies in that taxonomic composition was more influenced by dispersal limitation or species pool effects than functional composition (Soininen et al., 2016). Indeed, the distance decay of functional similarity was believed to have limited ability to reflect the spatial differences in biotic communities, probably because deterministic and stochastic factors have different roles in shaping community composition (Graco-Roza et al., 2022).

    In this study, the distance decay of turnover component was similar to that of total similarities. This is because turnover dominated the total beta diversity, and to some extent, turnover seemed to be a good surrogate of total beta diversity (Soininen et al., 2018). However, the distance decay of nestedness component showed contrasting pattern across three dimensions. Specifically, the distance decay of taxonomic nestedness increased with elevational distance, probably driven by the loss or gain of small-ranged species or rare species. Although counterintuitive, this relationship has been found previously (e.g., Fontana et al., 2020; Kusumoto et al., 2021; Wang et al., 2022); however, given the relatively low statistic power (Fig. 4), more comprehensive and more systematic studies would be very necessary.

    Further, our finding that strong relationship exists between assemblage similarity and climatic distance suggests that ecological niches shaped breeding bird assemblages in central Himalayan elevation gradient. This is especially true for birds because they are often considered as good dispersers, and seem not to be affected by dispersal limitation especially at small scales in this study (e.g., Araújo et al., 2005; Baselga and Gómez-Rodríguez, 2021), as indicated by previous studies that the degree to which the environment controls distance-decay patterns via niche processes varies with dispersal ability (e.g., Astorga et al., 2012; Saito et al., 2015; Gómez-Rodríguez and Baselga, 2018; Baselaga and Gómez-Rodríguez, 2021).

    Generally, there is a shift from clustered to overdispersed taxonomic, functional, and phylogenetic community structure (total dissimilarities and its turnover components) (Fig. 5). It should be noted that a relatively high proportion of SES values fall between −1.96 and 1.96, implying a product of random structure (Du et al., 2021; Wang et al., 2022). However, we could not tease apart randomness from deterministic (the balance of environmental filtering and limiting similarity) or neutral process (Fig. 1). Specifically, the SES of taxonomic, functional and phylogenetic total dissimilarities and turnover increased with increasing elevational distances, showing a gradual weakening trend of environmental filtering but strengthened trend of limiting similarity. This pattern is consistent with previous work on plants in the same study system (unpublished data), and is commonly reported in beta diversity along the elevational gradient in a wide range of taxa (e.g., Swenson et al., 2011; Dehling et al., 2014; Bishop et al., 2015; Tello et al., 2015; Tolmos et al., 2022). It seems clear that additional studies are needed to test the universality of our findings.

    In contrast, the SES of functional, and taxonomic and phylogenetic nestedness show no clear pattern, and a decreasing trend with the increase in elevational distances, respectively. Given that most SES values of functional and phylogenetic nestedness were between −1.96 and 1.96, we could infer that these communities were assembled by random process. However, nestedness components of phylogenetic dissimilarity still showed downward trend, indicating a gradual strengthened of clustering. This pattern could be jointly driven by the weaker competitive exclusion and stronger environmental filtering. Combined with the pattern of taxonomic nestedness, we could infer that phylogenetic losses/gains occurred frequently among distantly related species.

    We acknowledged that our line-transect surveys may be affected by imperfect detection of rare or elusive species residing in dense forests, leading to the underestimation and potential impacts on our observed beta diversity patterns (Si et al., 2018). In fact, the problem of imperfect detection is more likely to occur in areas with low and mid-elevations due to the prevalence of dense forests (Ding et al., 2019). Therefore, to obtain reliable results, continuous efforts on data collection are required (Feeley and Silman, 2011), and we should also improve the sampling design to better record rare species (Specht et al., 2017).

    Distance decay of similarity at multiple dimensions is a useful tool reflecting dispersal-related factors, and also for understanding species' response to climatic and spatially structured environmental variables (Soininen et al., 2007; Graco-Roza et al., 2022). When combining this method and beta diversity partitioning approach across multiple dimensions, additional insights on elevational patterns and assembly processes could be obtained (e.g., Du et al., 2021; Zhao et al., 2021; Wang et al., 2022). By comprehensively considering many aspects of beta diversity across a wide elevation range varying from 1800 to 5400 m, we found that turnover dominated taxonomic and phylogenetic beta diversity, whereas nestedness played an important role in functional dimensions. Additionally, distance decay of similarity decreased with the increase in elevational distance except taxonomic nestedness. Elevational beta diversity was primarily driven by deterministic process (from decreasing filtering to increasing limiting similarity), although we could not disentangle the effects of environmental filtering and limiting similarity (e.g., Du et al., 2021; Wang et al., 2022). Our results highlighted the niche-based deterministic processes in shaping elevational bird diversity patterns that were determined by the relative roles of decreasing trend of environmental filtering and increasing trend of limiting similarity along elevation distances.

    Ethical approval was not required since we did not manipulate any living animals. All data were obtained from field observation.

    Zhifeng Ding: Writing – review & editing, Writing – original draft, Visualization, Software, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. Jianchao Liang: Writing – review & editing, Investigation. Le Yang: Writing – review & editing, Funding acquisition. Cong Wei: Writing – review & editing, Funding acquisition. Huijian Hu: Writing – review & editing, Investigation, Funding acquisition, Conceptualization. Xingfeng Si: Writing – review & editing, Writing – original draft, Visualization, Software, Methodology, Formal analysis, Data curation, Conceptualization.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    We are grateful to Jingjing Li, Hongfen Cao, Li Xie and many other graduate students in our group for the field surveys, and thank the Qomolangma National Nature Reserve for the permits to conduct the research in the Gyirong Valley. We also thank Dr. Yuhao Zhao and Dr. Di Zeng at East China Normal University for constructive suggestions and patient helps on data analyses.

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.avrs.2024.100170.

  • Alatalo RV, Eriksson D, Gustafsson L, Larsson K. Exploitation competition influences the use of foraging sites by tits: experimental evidence. Ecology. 1987;68:284-90.
    Amarasekare P, Hoopes MF, Mouquet N, Holyoak M. Mechanisms of coexistence in competitive metacommunities. Am Nat. 2004;164:310-26.
    Barbaro L, van Halder I. Linking bird, carabid beetle and butterfly life-history traits to habitat fragmentation in mosaic landscapes. Ecography. 2009;32:321-33.
    Beven G. Changes in breeding bird populations of an oak-wood on Bookham Common, Surrey, over twenty-seven years. Lond Nat. 1976;55:23-42.
    Bibby CJ, Burgess ND, Hill DA, Mustoe S. Bird census techniques. Amsterdam: Elsevier; 2000.
    Bolnick DI, Amarasekare P, Araújo MS, Bürger R, Levine JM, Novak M, Rudolf VHW, Schreiber SJ, Urban MC, Vasseur DA. Why intraspecific trait variation matters in community ecology. Trends Ecol Evol. 2011;26:183-92.
    Brook BW, Buettel JC. Emigration is costly, but immigration has benefits in human-altered landscapes. Funct Ecol. 2016;30:1478-9.
    Brown JH, Kodric-Brown A. Turnover rates in insular biogeography: effect of immigration on extinction. Ecology. 1977;58:445-9.
    Brown WL, Wilson EO. Character displacement. Syst Zool. 1956;5:49-64.
    Brown BL, Sokol ER, Skelton J, Tornwall B. Making sense of metacommunities: dispelling the mythology of a metacommunity typology. Oecologia. 2017;183:643-52.
    Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2003.
    Chave J. Neutral theory and community ecology. Ecol Lett. 2004;7:241-53.
    Chen C, Holyoak M, Si X, Wang Y, Ding P. Do seasonal species assemblages differ in their biogeography? Evidence from the spatial structure of bird communities on land-bridge islands. J Biogeogr. 2018a;45:473-83.
    Chen C, Holyoak M, Si X, Wang Y, Ding P. Species traits, habitat characteristics and spatial dynamics of birds on land-bridge islands. China: Thousand Island Lake. Pangaea; 2018b. .
    Cook RR, Hanski I. On expected lifetimes of small-bodied and large-bodied species of birds on islands. Am Nat. 1995;145:307-15.
    Core Team R. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.
    Cottenie K. Integrating environmental and spatial processes in ecological community dynamics. Ecol Lett. 2005;8:1175-82.
    Dolédec S, Chessel D, ter Braak CJF, Champely S. Matching species traits to environmental variables: a new three-table ordination method. Environ Ecol Stat. 1996;3:143-66.
    Dray S, Dufour AB. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22:1-20.
    Dray S, Choler P, Dolédec S, Peres-Neto PR, Thuiller W, Pavoine S, ter Braak CJF. Combining the fourth-corner and the RLQ methods for assessing trait responses to environmental variation. Ecology. 2014;95:14-21.
    Driscoll DA, Lindenmayer DB. Empirical tests of metacommunity theory using an isolation gradient. Ecol Monogr. 2009;79:485-501.
    Ellis AM, Lounibos LP, Holyoak M. Evaluating the long-term metacommunity dynamics of tree hole mosquitoes. Ecology. 2006;87:2582-90.
    Farneda FZ, Rocha R, López-Baucells A, Groenenberg M, Silva I, Palmeirim JM, Bobrowiec PED, Meyer CFJ. Trait-related responses to habitat fragmentation in Amazonian bats. J Appl Ecol. 2015;52:1381-91.
    Fattorini S. The influence of geographical and ecological factors on island beta diversity patterns. J Biogeogr. 2010;37:1061-70.
    Fischer J, Lindenmayer DB. Nestedness in fragmented landscapes: a case study on birds, arboreal marsupials and lizards. J Biogeogr. 2005;32:1737-50.
    Fukami T. Historical contingency in community assembly: integrating niches, species pools, and priority effects. Annu Rev Ecol Evol Syst. 2015;46:1-23.
    Hanski I. Population dynamics of shrews on small islands accord with the equilibrium model. Biol J Linn Soc. 1986;28:23-36.
    Henle K, Davies KF, Kleyer M, Margules C, Settele J. Predictors of species sensitivity to fragmentation. Biodivers Conserv. 2004;13:207-51.
    Herrera CM. Ecological correlates of residence and non-residence in a mediterranean passerine bird community. J Anim Ecol. 1978;47:871-90.
    Holyoak M, Leibold MA, Holt RD. Metacommunities. Chicago: University of Chicago Press; 2005.
    Hubbell SP. The unified neutral theory of species abundance and diversity. Princeton: Princeton University Press; 2001.
    Hutto RL. Habitat selection by nonbreeding, migratory land birds. In: Cody ML, editor. Habitat selection in birds. London: Academic Press; 1985. p. 455-73.
    Klingbeil BT, Willig MR. Community assembly in temperate forest birds: habitat filtering, interspecific interactions and priority effects. Evol Ecol. 2016;30:703-22.
    Leck CF. Seasonal changes in feeding pressures of fruit- and nectar-eating birds in Panama. Condor. 1972;74:54-60.
    Legendre P, Galzin R, Harmelin-Vivien ML. Relating behavior to habitat: solutions to the fourth-corner problem. Ecology. 1997;78:547-62.
    Leibold MA, Chase JM. Metacommunity ecology. Princeton: Princeton University Press; 2017.
    Leibold MA, Holyoak M, Mouquet N, Amarasekare P, Chase JM, Hoopes MF, Holt RD, Shurin JB, Law R, Tilman D, Loreau M, Gonzalez A. The metacommunity concept: a framework for multi-scale community ecology. Ecol Lett. 2004;7:601-13.
    Leisler B. Habitat selection and coexistence of migrants and Afrotropical residents. Ibis. 1992;134:77-82.
    Levey DJ, Stiles FG. Evolutionary precursors of long-distance migration: resource availability and movement patterns in neotropical landbirds. Am Nat. 1992;140:447-76.
    Logue JB, Mouquet N, Peter H, Hillebrand H. The metacommunity working group. Empirical approaches to metacommunities: a review and comparison with theory. Trends Ecol Evol. 2011;26:482-91.
    MacKenzie DI, Nichols JD, Seamans ME, Gutiérrez RJ. Modeling species occurrence dynamics with multiple states and imperfect detection. Ecology. 2009;90:823-35.
    MacKinnon JR, MacKinnon J, Phillipps K, He FQ. A field guide to the birds of China. Oxford: Oxford University Press; 2000.
    Marra PP, Hobson KA, Holmes RT. Linking winter and summer events in a migratory bird by using stable-carbon isotopes. Science. 1998;282:1884-6.
    Meynard CN, Quinn JF. Bird metacommunities in temperate South American forest: vegetation structure, area, and climate effects. Ecology. 2008;89:981-90.
    Murcia C. Edge effects in fragmented forests: implications for conservation. Trends Ecol Evol. 1995;10:58-62.
    Murgui E. Effects of seasonality on the species-area relationship: a case study with birds in urban parks. Glob Ecol Biogeogr. 2007;16:319-29.
    Murgui E. Seasonality and nestedness of bird communities in urban parks in Valencia, Spain. Ecography. 2010;33:979-84.
    Nice MM. The role of territory in bird life. Am Midl Nat. 1941;26:441-87.
    Norris DR, Marra PP. Seasonal interactions, habitat quality, and population dynamics in migratory birds. Condor. 2007;109:535-47.
    Nour N, van Damme R, Matthysen E, Dhondt AA. Forest birds in forest fragments: are fragmentation effects independent of season? Bird Study. 2010;46:279-88.
    Özkan K, Svenning JC, Jeppesen E. Environmental species sorting dominates forest-bird community assembly across scales. J Anim Ecol. 2013;82:266-74.
    Pavoine S, Baguette M, Stevens VM, Leibold MA, Turlure C, Bonsall MB. Life history traits, but not phylogeny, drive compositional patterns in a butterfly metacommunity. Ecology. 2014;95:3304-13.
    Peltonen A, Hanski I. Patterns of island occupancy explained by colonization and extinction rates in shrews. Ecology. 1991;72:1698-708.
    Pimm SL. The balance of nature? Ecological issues in the conservation of species and communities. Chicago: University of Chicago Press; 1991.
    Rayner L, Lindenmayer DB, Wood JT, Gibbons P, Manning AD. Are protected areas maintaining bird diversity? Ecography. 2014;37:43-53.
    Russell GJ, Diamond JM, Reed TM, Pimm SL. Breeding birds on small islands: island biogeography or optimal foraging? J Anim Ecol. 2006;75:324-39.
    Schoener TW. Sizes of feeding territories among birds. Ecology. 1968;49:123-41.
    Schoereder JH, Galbiati C, Ribas CR, Sobrinho TG, Sperber CF, DeSouza O, Lopes-Andrade C. Should we use proportional sampling for species-area studies? J Biogeogr. 2004;31:1219-26.
    Si X, Pimm SL, Russell GJ, Ding P. Turnover of breeding bird communities on islands in an inundated lake. J Biogeogr. 2014;41:2283-92.
    Si X, Cadotte MW, Zeng D, Baselga A, Zhao Y, Li J, Wu Y, Wang S, Ding P. Functional and phylogenetic structure of island bird communities. J Anim Ecol. 2017;86:532-42.
    Song X, Holt RD, Si X, Christman MC, Ding P. When the species-time-area relationship meets island biogeography: diversity patterns of avian communities over time and space in a subtropical archipelago. J Biogeogr. 2018;45:664-75.
    Stracey CM, Pimm SL. Testing island biogeography theory with visitation rates of birds to British islands. J Biogeogr. 2009;36:1532-9.
    Tellería JL, Santos T. Seasonal and interannual occupation of a forest archipelago by insectivorous passerines. Oikos. 1997;78:239-48.
    Tramer EJ. Proportions of wintering North American birds in disturbed and undisturbed dry tropical habitats. Condor. 1974;76:460-4.
    Wang Y, Bao Y, Yu M, Xu G, Ding P. Nestedness for different reasons: the distributions of birds, lizards and small mammals on islands of an inundated lake. Divers Distrib. 2010;16:862-73.
    Wang Y, Chen S, Ding P. Testing multiple assembly rule models in avian communities on islands of an inundated lake, Zhejiang Province, China. J Biogeogr. 2011;38:1330-44.
    Wang Y, Thornton DH, Ge D, Wang S, Ding P. Ecological correlates of vulnerability to fragmentation in forest birds on inundated subtropical land-bridge islands. Biol Conserv. 2015;191:251-7.
    White EP, Hurlbert AH. The combined influence of the local environment and regional enrichment on bird species richness. Am Nat. 2010;175:E35-43.
    Wiens JA. The ecology of bird communities. Cambridge: Cambridge University Press; 1992.
    Zhuge Y. Fauna of Zhejiang: aves. Zhejiang: Zhejiang Science and Technology Publishing House; 1990.
    Zuckerberg B, Fink D, La Sorte FA, Hochachka WM, Kelling S. Novel seasonal land cover associations for eastern North American forest birds identified through dynamic species distribution modelling. Divers Distrib. 2016;22:717-30.
  • Related Articles

  • Cited by

    Periodical cited type(8)

    1. Junjian Zhang, Xianghuang Li, Nyambayar Batbayar, et al. Exploring potential reasons for the increase in the East Asian Greylag Goose Anser anser population by assessing habitat use and use of protected areas. Bird Conservation International, 2023, 33 DOI:10.1017/S095927092200017X
    2. Baasansuren Erdenechimeg, Gankhuyag Purev-Ochir, Amarkhuu Gungaa, et al. Migration Pattern, Habitat Use, and Conservation Status of the Eastern Common Crane (Grus grus lilfordi) from Eastern Mongolia. Animals, 2023, 13(14): 2287. DOI:10.3390/ani13142287
    3. Rui Ma, Shuo Ma, Xu Wei, et al. Tracking Migration of Eastern Spot-Billed Ducks Anas zonorhyncha and Mallards Anas platyrhynchos Wintering in Shanghai, China. Waterbirds, 2023, 45(4) DOI:10.1675/063.045.0411
    4. Xin Yu, Ping Fan, Yun Wu, et al. GPS tracking data reveal the annual spatiotemporal movement patterns of Bridled Terns. Avian Research, 2022, 13: 100065. DOI:10.1016/j.avrs.2022.100065
    5. Iderbat Damba, Junjian Zhang, Kunpeng Yi, et al. Seasonal and regional differences in migration patterns and conservation status of Swan Geese (Anser cygnoides) in the East Asian Flyway. Avian Research, 2021, 12(1) DOI:10.1186/s40657-021-00308-y
    6. Fanjuan Meng, Xin Wang, Nyambayar Batbayar, et al. Consistent habitat preference underpins the geographically divergent autumn migration of individual Mongolian common shelducks. Current Zoology, 2020, 66(4): 355. DOI:10.1093/cz/zoz056
    7. Ye Wang, Chunrong Mi, Yumin Guo. Satellite tracking reveals a new migration route of black-necked cranes (Grus nigricollis) in Qinghai-Tibet Plateau. PeerJ, 2020, 8: e9715. DOI:10.7717/peerj.9715
    8. Ros M.W. Green, Niall H.K. Burton, Aonghais S.C.P. Cook. Migratory movements of British and Irish Common Shelduck Tadorna tadorna: a review of ringing data and a pilot tracking study to inform potential interactions with offshore wind farms in the North Sea. Ringing & Migration, 2019, 34(2): 71. DOI:10.1080/03078698.2019.1887670

    Other cited types(0)

Catalog

    Figures(4)  /  Tables(2)

    Article Metrics

    Article views (257) PDF downloads (9) Cited by(8)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return