Yun Chen
Shuai Niu2- 1College of Forestry, Henan Agricultural University, Zhengzhou, China
- 2College of Life Sciences, Henan Agricultural University, Zhengzhou, China
- 3Educational Administration Department, Henan Finance and Taxation College, Zhengzhou, China
- 4Xiaoqinling National Nature Reserve, Lingbao, China
Elucidating the major drivers of bryophyte distribution is the first step to protecting bryophyte diversity. Topography, forest, substrates (ground, tree trunks, roots, rocks, and rotten wood), and spatial factor, which factors are the major drivers of bryophyte distribution? In this study, 53 plots were set in 400 m2 along the elevation gradient in Xiaoqinling, China. All bryophytes in the plots were collected and identified. Regression analysis was used to examine the relationship between bryophyte and substrate diversity. We compared the patterns of overall bryophyte diversity and diversity of bryophytes found on the ground, tree, and rock along elevational gradients. Canonical correspondence analysis was applied to relate species composition to selected environmental variables. The importance of topography, forest, substrates, and spatial factors was determined by variance partitioning. A total of 1378 bryophyte specimens were collected, and 240 species were identified. Bryophyte diversity was closely related to substrate diversity. The overall bryophyte diversity significantly increased with elevation; however, the response varied among ground, tree, and rock bryophytes. Tree diversity and herb layer were considered important environmental factors in determining bryophyte distribution. Species abundance was best explained by stand structure (17%), and species diversity was best explained by stand structure (35%) and substrate (40%). Results directly indicated that substrate diversity can improve bryophyte species diversity. The effects of micro-habitat formed by stand structure and substrate diversity were higher than those of spatial processes and topography factors on bryophyte distribution. This study proved that the determinant factors influencing bryophyte diversity reflect the trends in recent forest management, providing a real opportunity to improve forest biodiversity conservation.
Introduction
Bryophytes are undoubtedly one of the earliest and most diversified groups of land plant, constituting an important part in vegetation diversity (Song et al., 2015). Bryophytes play important roles in the ecosystem. In particular, these species are involved in water balance and nutrient accumulation (Ódor et al., 2013). However, bryophytes are smaller than vascular plants and are often ignored in ecological studies (Stehn et al., 2010).
Elucidating the major drivers of bryophyte distribution is the firststep to conserving its diversity (Smith and Stark, 2014). Ecologists have conducted numerous explorations regarding the major drivers of vascular plant distribution. Most of the studies concluded that the coexistence of vascular plants is the joint result of niche and neutral theories (Jia et al., 2015; Chen et al., 2016). In recent years, neutral and niche theories have attracted considerable attention to explain this mechanism (Barot, 2004; Jia et al., 2015). Niche theory states that different species adapt to different habitats and are restricted by various factors (Hutchinson, 1957; Valladares et al., 2015; Conradi et al., 2017). Meanwhile, the neutral theory states that species have equal competitive abilities, and thus species patterns are only generated by dispersal limitation (Hubbell, 2001; Chen et al., 2016; Hidalgo et al., 2017). Bryophytes are reproduced by spores spreading, and dispersal limitation is one of the key processes determining the spatial structure of bryophyte distribution. Assuming that dispersal is limited in space (Hubbell, 2001), and environmental heterogeneity has been well sampled (Yuan et al., 2011), the existence of spatially structured variations of bryophyte distribution independent of the environmental fraction can indicate the importance of dispersal processes. However, the studies on the relative importance of environmental heterogeneity and spatially structured processes in explaining bryophyte distribution are relatively few. Although some ecologists conducted a few positive exploration (Mota de Oliveira et al., 2009; Silva et al., 2014; Smith and Stark, 2014), the major drivers of bryophyte distribution remain unclear.
What are the major drivers of bryophyte distribution? Smith and Stark (2014) investigated an arid desert that showed that species distribution of the bryophyte community is attributed to the two theories. Silva et al. (2014) studied the β-diversity of bryophytes on rocky outcrops and reported that randomization plays an important role in their distribution. Mota de Oliveira et al. (2009) studied the coexistence of trunk-epiphytic bryophytes and found that the role of the niche process is greater than that of the neutral process. However, these studies did not consider the roles of micro-habitats formed by stand structure and substrate diversity.
The spatial pattern of species along the elevational gradients is a basic pattern in biogeography, ecology, and conservation biology (Gaston, 2000), which also provides fundamental insights into the mechanisms of species distribution (Anderson et al., 2011). The species diversity of bryophytes along the elevation gradients has four major patterns: increasing richness with increasing elevation (Gradstein et al., 1989; Stehn et al., 2010), decreasing diversity with increasing elevation (Tusiime et al., 2007), hump-shaped distribution (Wolf, 1993; Ah-Peng et al., 2012; Song et al., 2015), or absence of any clear altitude trend (Andrew et al., 2003; Bruun et al., 2006; Sun et al., 2013). However, the growth and spread of bryophytes may be different among species that grow on different substrates, such as ground, tree, and rock bryophytes (During, 1979; Spitale, 2016). The patterns of overall bryophyte diversity may mask the true patterns of bryophyte diversity on different substrates. The study of the pattern comparison of bryophyte diversity on ground, tree, and rock is still rare.
Bryophytes inhabit a number of different substrates, such as trees, dead wood, rock, and ground in forests (Scott, 1994). The roles of substrate diversity in forest ecosystems in maintaining bryophyte diversity are frequently acknowledged (Turner and Pharo, 2005). Based on the substrate types, researchers conduct many positive explorations regarding the diversity and distribution of bryophytes (Soderstrom, 1993; Rambo and Muir, 1998; Pharo and Beattie, 2002; Szövényi et al., 2004). However, the difference of major drivers of bryophyte distribution among bryophytes on ground, tree, and rock remains unclear.
Forest management aims to harmonize timber production with biodiversity conservation to accomplish the sustainable forest management (Ódor et al., 2013). Biodiversity conservation is one of the important goals of sustainable forest management. Bryophytes are important part of biodiversity. Bryophytes can sensitively respond to environmental changes, and it is also an indicator for forest management (Nascimbene et al., 2014). The aim of this study is to improve the ecological understanding of bryophyte species distribution in forests in the context of sustainable forest management. We hope that this study can provide a reference for the protection of species diversity.
Xiaoqinling is a national nature reserve, a typical temperate deciduous broad-leaved forest. The communities have high diversity, clear layers of trees, shrubs, and herbs (Ye et al., 2004). Human disturbance is little, and bryophytes in Xiaoqinling are relatively rich (Ye et al., 2004). The highest point in the Nature Reserve is 2413.8 m, which is also the highest in the Henan province (Ye et al., 2004). So it is a natural platform for studying the factors driving the bryophyte species distribution. In this study, we set up 53 plots in 400 m2 along the altitudinal gradient in Xiaoqinling National Nature Reserve. We collected and identified bryophytes in all substrates within the plots. We focused on three issues: (1) How many bryophyte species exist and the difference of bryophyte diversity among different substrates; (2) What is the distribution pattern of bryophyte diversity along elevational gradients and the differences in the patterns of bryophytes among different substrates; (3) Topography, forest, substrates, and spatial factor, which factors are the major drivers of bryophyte distribution and the conservation implications for forest management in the forest ecosystem?
Materials and Methods
Study Area
Xiaoqinling National Nature Reserve is located in Lingbao (34° 23′–34° 31′ N, 110° 23′–110° 44′ E in Henan Province, China). The nature reserve covers 31 km from east to west and 12 km from north to south, and the total area is 15160 ha. The elevation is 308–2413.8 m, the annual average temperature is 5.9–14.0°C, the frost-free period is 170–215 days, and the annual precipitation is 506–719.2 mm (Ye et al., 2004). The most common soil types are alluvial soils and brown soil (pH of the 0–30 cm layer is 6.2–7.0) (Ye et al., 2004).
The area has rich biological resources with a diverse and high vegetation composition. Four vegetation belts exist along with elevation: scrub meadow zone (less than 1000 m), deciduous broad-leaved forest zone (1100–1700 m), mixed coniferous broad-leaved forest zone (1700–2000 m), and elfin-wood zone (more than 2000 m) (Ye et al., 2004). The forest coverage rate is 81.2%. The dominant species in forest were Quercus variabilis Bl., Quercus aliena Bl., Quercus glandulifera var. brevipetiolata Nakai, Betula luminifera H. Winkl., Betula albosinensis Burk., and Pteroceltis tatarinowii Maxim.
Sampling Design
We established 53 20 m × 20 m plots along the elevation gradient from the foothills (1020 m) to the peak of the mountain (2413 m) in August 2012. And data collected included species name of trees, number of individual trees, diameter at breast height (DBH), and height of all trees with DBH ≥ 1 cm. A ruler was used to measure the plant canopy length of east–west (EW) and north–south (NS) directions. Tree crown was calculated as EW(m) × NS(m), where EW is the plant canopy length of the EW direction and NS is the plant canopy length of the NS direction (Martens et al., 1991). A shrub community of 5 m × 5 m subplot was randomly selected within each 20 m × 20 m plot to record all shrub species, the number of individuals, and coverage. The height of each shrub species also was recorded. Four herb plots in 1 m × 1 m plot in each corner of the tree plot were also chosen, in which the data collected includes herbs species name, abundance, coverage, and height of each herb species. The longitude, latitude, and elevation of the plots were recorded by GPS. A geological compass and a clinometer were used to record aspect and slope, respectively.
We collected all bryophyte species in 53 plots of 20 m × 20 m and recorded substrate types, including ground, dead trees, rotten woods, dead leaves, trunks, roots, and rocks. All bryophyte species, including the different species and the same species at different points in the plot, were collected. We carefully checked the plot from every corner of a plot to avoid the omissions of bryophyte collections. Specimen collectors have had professional training before specimen collection. The training included basic characteristics of bryophytes, field identification of bryophytes, and field specimen collection methods. A total of 1378 specimens were collected. All specimens were identified by the microscope. Voucher specimens were deposited in the bryophyte herbarium of the Institute of Henan Agricultural University.
Stand Structure, Topographical, and Spatial Data
A total of 12 variables, namely, tree stand density, tree basal area, tree crown, tree diversity, shrub density, shrub diversity, shrub height, shrub crown, herb density, herb diversity, herb cover, and herb height, were used to represent the stand structure. Tree stand density is the number of individual trees in a plot. Tree basal area was calculated as π × R2, where R is the radius at the height of 1.3 m. Tree crown is the sum of crowns in a plot. Simpson’s diversity index (Hunter and Gaston, 1988) was used to measure tree diversity. Indexes for shrub and herb were calculated using methods similar to those of a tree. Shrub and herb height are species average height. Shrub and herb density is the number of individual in a plot. Moreover, the elevation, slope, and aspect of every plot were used to represent topographical variables. Slope and aspect were converted into sinusoidal form during analysis.
Spatial variables were obtained by principal coordinates of neighbor matrices (PCNM) according to Borcard and Legendre (2002) and Legendre et al. (2009). PCNM variables represent the spatial relationship between the plots. These variables more precisely represent spatial patterns than the Euclidean distance matrix, geographic coordinates, and cubic trend surface equation (Borcard and Legendre, 2002). This method is based on the calculation of a principal coordinate analysis from a truncated matrix of Euclidean distances among plots. The PCNM characteristic function was obtained by “spacemakeR” package in R language software (Blanchet et al., 2008).
Statistical Procedures
Regression analysis is one of the most robust methods for the relations between the variables. We examined how species diversity of all bryophytes was related to elevation using regression analysis, with significance at p < 0.05. Simpson’s diversity index measured species diversity. This index is one of the most commonly used methods for determining species diversity (Hunter and Gaston, 1988). We also examined the relation between species diversity of all bryophytes and substrate diversity using regression analysis. Based on the number of substrates in a plot, such as rock, ground, tree trunks, roots, and rotten wood, Simpson’s diversity index also was used to calculate the substrate diversity, so that the differences in different calculation methods can be avoided.
All bryophytes in the plots are classified into three main groups based on substrate: ground bryophytes, tree bryophytes, and rock bryophytes. The bryophytes found on the dead trees, rotten woods, dead leaves, trunks and roots are relatively few in number, and closely related to trees. So the bryophytes are classified as tree bryophytes. If there are no bryophytes in some plots, the plots will not be considered in the analysis process. Because of the substrate diversity have no method to calculate if there are no bryophytes in a plot. We examined the relation between substrate diversity and species diversity of ground bryophytes, tree bryophytes, and rock bryophytes using regression analysis. Linear regression was conducted in R 2.15.3 with Base Package (R Development Core Team).
Canonical correspondence analysis (CCA) is a simple method for arranging species along environmental variables. CCA can be used both for detecting species–environment relations and investigating specific questions regarding the response of species to environmental variables (Ter Braak, 1987). Prior to CCA analyses, detrended correspondence analysis (DCA) was conducted to detect the length of the first two axes and determine whether the linear or unimodal model should be selected. Unimodal distribution model, i.e., CCA, was selected because DCA showed a maximum axis length of more than four (Legendre and Legendre, 1998). In this study, CCA was used to determine the effects of stand structure, topographical factors, spatial factors, and substrate diversity on the bryophyte community. In CCA, species matrix is the species abundant in plots (number of individuals is more than one). Bryophytes are presence or absence at some point in one plot, abundant data of bryophytes can be obtained. The environmental matrix comprises stand structure, topographic variables, spatial variables, and substrate diversity. Stand structure includes tree stand density, tree basal area, tree crown, tree diversity, shrub density, shrub diversity, shrub height, shrub crown, herb density, herb diversity, herb cover, and herb height. Topographic variable includes elevation, slope, and aspect. Spatial variables represented by PCNM variables. Fifteen spatial PCNM variables were created in this study and used to model spatial structure at different scales. PCNM1 represents the spatial information of the entire research scale which belongs to a broad scale. PCNM2 to PCNM15 describe the fine-scale information. Stand structure, topographic variables, and spatial variables were forward selected by the “forward.se” function in R language software with “packfor” packages. Thus, a few environmental factors were removed, and the model was simplified. Moreover, Monte Carlo permutation test was used to analyze whether the model reached a significant level (P < 0.05). The “envfit” function in the R language software with “vegan” package (Oksanen et al., 2007) was used to test the significance of each environmental factor and bryophyte species distribution. The same methodology was used to determine the effects of stand structure, topographical factors, and spatial factors on the ground bryophytes, tree bryophytes, and rock bryophytes, respectively. PCNM variables were created by “PCNM” packages (Blanchet et al., 2008). CCA and DCA analyses were conducted in R 2.15 with “vegan” package (Oksanen et al., 2007) (R Development Core Team1).
Variance partitioning is useful in evaluating the major drivers of different species in community ecology (Violle et al., 2012; Coyle, 2017). Variation partitioning was performed to divide the variation in bryophytes species distribution and diversity, respectively, and explore the influence of spatial processes and topography factors, stand structure, and substrate diversity on bryophyte species distribution and species diversity. The dependent variable in variance partitioning is bryophyte species abundance matrix, which was transformed by “hellinger” (Legendre and Gallagher, 2001), or species diversity, which was calculated by Simpson’s diversity index. Hellinger transformations are appropriate alternatives giving low weights to rare species (Legendre and Gallagher, 2001). And this transformation might provide better resolutions for species by making them behave more like a Gaussian distribution (Griffith and Peres-Neto, 2006). Variance partitioning was also used to determine the effects of stand structure, topographical factors, and spatial factors on the ground bryophytes, tree bryophytes, and rock bryophytes, respectively. Variance partitioning was computed using the “vegan” library (Oksanen et al., 2007) of the R statistical language.
Results
Bryophyte Diversity in Xiaoqinling
Based on the 1378 bryophyte specimens collected in 53 plots, a total of 240 bryophyte species were collected, which belong to 33 families and 89 genera (Supplementary Table S1). Bryophyte abundance in 53 plots is provided in Supplementary Table S2. A total of 126 trees, 130 shrubs, and 263 herbs species were recorded. Statistics of the bryophyte diversity and stand structure are provided in Table 1. The bryophyte diversity (0.89) is second only to herbs diversity (0.92) and is far higher than the diversity of trees (0.75) and shrubs (0.77) (Table 1). Therefore, bryophytes have important significance for species diversity in this region.
TABLE 1. Overview of bryophytes, trees, shrubs, and herbaceous plants in 53 plots.
Bryophyte diversities remarkably differed among different substrates, and 57.5% species tended to be distributed in one kind of substrate. The numbers of bryophyte species in different substrate types are illustrated in Figure 1A. Species of saxicolous bryophytes are the richest species. A total of 72, 26, 19, 18, and 3 bryophyte species were recorded in the rock, ground, tree trunks, roots, and dead trees, respectively. Rock bryophytes are the richest. The regression analysis indicated that the overall bryophyte diversity significantly increased with the substrate diversity (Figure 1B). For overall bryophyte diversity, substrate diversity explained 17.20% of the total variance (p = 0.0020).
FIGURE 1. The number of bryophyte species in different substrates (A). (B) Shows the relationship between the bryophytes species diversity (log-transformed) and substrate diversity. Simpson’s diversity index was used to measure species diversity and substrate diversity, respectively. Substrate diversity was calculated based on the number of substrates in a plot.
Elevation Patterns of Species Diversity
The overall bryophyte diversity significantly increased with the elevation (R2 = 0.1738, p = 0.0023); however, bryophyte diversity varied among ground, tree, and rock bryophytes (Figures 2A–D). The tree bryophyte diversity exhibited highly similar trends along the elevation gradient with the overall bryophyte diversity (R2 = 0.2900, p = 0.0002). The diversity of ground (p = 0.6290) and rock bryophytes (p = 0.4200) did not have a significant correlation to elevation.
FIGURE 2. Species diversity patterns of (A) overall bryophytes, (B) ground bryophytes, (C) tree bryophytes, and (D) rock bryophytes along the elevation. Simpson’s diversity index was used to measure species diversity.
Influence of Environmental Factors on Bryophyte Communities
The CCA ordination of all plots and bryophyte species indicated that variation in the species composition of bryophytes was mainly explained by the stand structure (Trees crown, herb diversity, herb density, and herb height) and partly by spatial factors (PCNM1, PCNM2, and PCNM12) and elevation (Figure 3). In the ordination diagram, environmental variables explained 30.93% of the variation in the structure of the assemblages. The first CCA axis (eigenvalue 0.30714) was positively correlated with elevation and negatively correlated with large-scale spatial patterns (PCNM1) and herb density. The second CCA axis (eigenvalue 0.23170) represented a substrate diversity gradient. In the ordination analysis, the significant variables in the permutation test were elevation, trees crown, herb density, herb diversity, herb height, PCNM1, PCNM2, and PCNM12 (Table 2); however, elevation, trees crown, herb diversity, and height are the most important, whereas the effects of aspect and shrub height were less influential.
FIGURE 3. CCA ordination diagram of the bryophytes data. Yellow dots and blue numbers represent species and plots, respectively. Arrows represent environmental factors. The longer the arrows represent the correlation is larger or smaller. We used the PCNM approach to introduce space as an explanatory variable. The analysis was conducted using bryophyte species abundant. Numbers on the CCA1 or CCA2 is the eigenvalue value. Abbreviations: TSD, tree stand density; TBA, tree basal area; TC, trees crown; TD, tree diversity; SDE, shrub density; SDI, shrub diversity; SH, shrub height; SC, shrub crown; HDE, herb density; HDI, herb diversity; HH, herb height; HC, herb cover; SubD, substrate diversity; pcnm, spatial PCNM variables.
TABLE 2. Significant explanatory variables in the canonical correspondence analysis (CCA).
For ground, tree, and rock bryophytes, environmental variables explained 35, 35, and 27% of the variation in the structure of the bryophytes assemblages, respectively. The CCA ordination indicated that spatial factors are the most important for ground, tree, and rock bryophytes. Meanwhile, the effects of elevation are also important factors for tree bryophytes, and the effects of shrub and herb are more important factors for ground and rock bryophytes, respectively (Figure 3).
The Major Drivers of Bryophyte Distribution
For species abundance, variation partitioning results show that topographical factors, stand structure, spatial factors, and substrate diversity explained 30.93% of the variation in the overall bryophyte distribution. Of this, 16% was attributed to pure stand structure, 6% was attributed to pure spatial factors, 4% was attributed to pure topographical variables, and 2% was attributed to pure substrate diversity. The effect of stand structure was relatively remarkable (Figure 4A). For ground and tree bryophytes, the effects of pure spatial factors on bryophyte distribution (16 and 19%, respectively) were slightly higher than those of pure stand structure (14 and 13%), and the effects of pure topographical factors on bryophyte distribution were the lowest (2 and 2%, respectively). For rock bryophytes, the effects of pure stand structure on bryophyte distribution were the highest (15%).
FIGURE 4. Variance partitioning for the effects of topographical factors, spatial factors, stands structure, and substrate diversity on bryophyte species abundance (A) and species diversity (B). Values less than zero are not shown. Abundant data was calculated based on the presence or absence of bryophyte species at some point in one plot. Species diversity was measured by Simpson’s diversity index. Spatial variables represented by PCNM variables.
For species diversity, variation partitioning results show that topographical factors, stand structure, spatial factors, and substrate diversity explained 54% of the variation in the overall bryophyte diversity. Of this, 16% was solely attributed to substrate diversity, 10% was the joint effect of topographic and stand structure, and 8% was the joint effect of topography, spatial factors, and substrate diversity (Figure 4B). For ground bryophytes, the effects of pure stand structure on bryophyte diversity (18%) were higher than those of pure spatial factors (12%). For tree bryophytes, the effects of pure stand structure on bryophyte diversity were the lowest (2%). For rock bryophytes, the effects of pure spatial factors on bryophyte diversity (9%) were higher than those of pure stand structure (7%).
Discussion
Substrate Diversity
Our results indicated that substrate diversity can improve bryophyte species diversity. The diversity of bryophytes in forest ecosystems is related to forest management practices (Hofmeister et al., 2015). As forest management leads to insufficiency of old-growth forest attributes, bryophyte diversity has declined in most managed forests (Söderström, 1988; Friedel et al., 2006). However, Lindenmayer and Laurance (2012) argued that the effects of management practices are spatially and temporarily restricted. Hofmeister et al. (2015) studied European forests and found that the bryophyte species richness in managed forests was less than in unmanaged forests. Moreover, their study showed that only unmanaged forest stands that provided suitable substrates were able to support rare bryophyte species. Hofmeister et al. (2015) indirectly showed the importance of maintaining the substrate diversity by comparing the managed and unmanaged forests.
Species Spatial Patterns along the Elevation Gradient
Temperature, precipitation, and air humidity remarkably differ along elevation gradients (Lieberman et al., 1996). Our results showed that bryophyte diversity exhibited an increasing richness pattern with increasing elevation. Our results are also consistent with those of Gradstein et al. (1989) and Stehn et al. (2010). This is probably because bryophyte diversity increases with elevation as a result of increased moisture availability (Lee and Roi, 1979) and the ability of bryophytes to tolerate extreme conditions (Bruun et al., 2006) and partly because the degree of human disturbance is small at high elevation.
Our results show that ground, tree, and rock bryophytes exhibited different elevation patterns. The study of Söderström (1988) showed that many bryophyte species are therefore absent in forests with an uneven supply of substrate. If dispersal is limited, bryophytes will not, or rarely, appear (Söderström, 1988). Tree bryophyte diversity exhibited an increasing richness pattern with increasing elevation. This may because tree bryophytes mainly grow on the trunk and was closely related to the trees (Benítez et al., 2015). Trees are the main part of the forest ecosystem and are continuously distributed with increasing elevation. Accordingly, the growth and spread of tree bryophytes provide continuous space. However, the rocks are randomly distributed with increasing elevation, and dispersal limitation perhaps influence the spread of the bryophyte species. Therefore, the diversity of rock bryophytes did not have a significant correlation to elevation. For ground bryophytes, although it has a continuous space with species spread, it is strongly influenced by environmental spatial heterogeneity. For example, a few bryophytes prefer more shaded condition and high air humidity (Hosokawa and Odani, 1957; Lesica et al., 1991; Humphrey et al., 2002; Ranius et al., 2008), whereas some bryophytes prefer more open conditions as forest edges (Moe and Botnen, 2000; Vanderpoorten et al., 2004). Hence, the routes of ground bryophytes spread may are still restricted. Therefore, the diversity of ground bryophytes did not have a significant correlation to elevation. In this study, we discussed the distribution pattern of bryophyte only from the angle of the substrate. Other possible factors, such as the spread of bryophyte spores, were not considered here.
Bryophyte Composition and Environmental Variables
The CCA ordination of all bryophytes indicated that variation in the species composition of bryophytes was mainly explained by the stand structure, including the diversity and crown of the tree as well as the coverage, height, density, and diversity of the herb layer. The diversity of trees can provide a diverse habitat for bryophytes (Rosabal et al., 2013), and previous research has shown that tree diversity has a close relation with bryophyte diversity (Pharo et al., 1999, 2000; Chiarucci et al., 2007). Coverage, height, density, and diversity of the herb layer could elevate local air humidity (Aude and Poulsen, 2000; Brunialti et al., 2010). Therefore, tree diversity and herb layer were considered important environmental factors in determining the bryophyte distribution. Although the diversity and crown of the tree are important factors for overall bryophytes, our results show that no tree structure factors significantly affect tree bryophyte distribution in this region. Chiarucci et al. (2007) concluded their study on vascular plants and bryophytes which have similar results with us. The relations between tree structure factors and tree bryophyte diversity needs to study and investigate further.
Furthermore, the CCA ordination indicated that spatial variables on different scales significantly affected the distribution of overall bryophytes, signifying that local dispersal ability is a key determinant of the bryophyte community structure. These findings agree with the observations of local dispersal limitation in deserts (Smith and Stark, 2014), bogs (Andersen et al., 2011), forests (Löbel et al., 2006), and land (Chen et al., 2015) bryophyte communities. However, spatial variables measured are notably not equal to the dispersal capabilities of species. Many scholars have conducted active exploration in the dispersal capabilities of bryophyte (Miles and Longton, 1992; Nathan and Muller-Landau, 2000). This will help us further understand the effects of the spatial process on bryophyte distribution.
Community Assembly
Our results indicated that the effects of stand structure and substrate diversity on bryophyte distribution were higher than those of spatial processes and topography factors. Stand structure exhibited high explanatory power because trees, shrubs, and herbs constitute different micro-habitats (Mežaka et al., 2012). Smith and Stark (2014) showed that species distribution in the bryophyte community is attributed to the joint result of niche and neutral theories. Mota de Oliveira et al. (2009) studied the coexistence of trunk-epiphytic bryophytes and showed that the role of the niche process is greater than that of the neutral process. Our research is based on the niche and neutral theories, and we also considered the stand structure and substrate diversity. Our results are also inconsistent with those of Mota de Oliveira et al. (2009) and Smith and Stark (2014). Our results showed that the effects of micro-habitat formed by stand structure and substrate diversity were higher than those of spatial processes and topography factors on bryophyte community assembly. Vitt et al. (1995) studied the relative importance of the micro-habitat and climate in the Peatlands of Continental Western Canada, which showed that bryophyte richness is most highly correlated with habitat heterogeneity whereas climatic factors are insignificant. Thus, our results are consistent with those of Vitt et al. (1995). Micro-habitat formed by stand structure and substrate diversity are important factors for bryophyte species distribution.
When bryophytes in plots were classified into three main groups, variation partitioning results showed that the effects of spatial variables had increased. So stand structure and spatial processes are all important factors for ground, tree, and rock bryophytes. The mechanisms of community assembly may differ among ground, tree, and rock bryophytes (Halpern et al., 2014; Spitale, 2017). Our results showed that the effects of spatial processes on tree and ground bryophytes community were higher than rock bryophytes community. This may because the distribution of rock is more random and uneven. And trees are continuously distributed with increasing elevation. The distribution of tree bryophytes is affected by the spatially structured variations of tree community. Ground bryophytes also have a continuous space with species spread although it is influenced by environmental spatial heterogeneity. Therefore, the effects of spatial processes on tree and ground bryophytes community were higher than rock bryophytes community.
The study of Chang et al. (2013) showed that better environmental data could reverse the conclusions regarding community assembly processes. Substrate availability (Mills and Macdonald, 2005), climate and soil factors (Heilmann-Clausen et al., 2014), and the height of the host tree (Oliveira and Steege, 2015), as the important environmental factors for bryophyte distribution, were not considered here. Therefore, more environmental factors must be considered in future studies.
Implications for Species Conservation
In the perspective of sustainable forest management, this study proved that the determinant factors influencing bryophyte diversity reflect on trends in recent forest management (trees crown and diversity, coverage, height, density and diversity of herb layer, and substrate diversity), providing a real opportunity to improve forest biodiversity conservation. Based on our results, the main strategy of management focusing on diversity conservation should include the following: (1) maintenance of tree species diversity. Because of the diversity of trees can provide a diverse habitat for bryophytes. (2) Presence of herb and regeneration layer. This can elevate local air humidity and suitable for bryophyte growth. (3) The maintenance of substrate diversity, rock, and dead wood in the forest must be retained. More substrates can provide continuous and diversified spaces for bryophyte growth, and reduce the dispersal limitation. (4) Ground, tree, and rock bryophytes perhaps should receive different protection strategies. Because of the major drivers of bryophytes distribution have some difference among ground, tree, and rock bryophytes. (5) Protection of bryophyte species diversity may require propagule supplementation for individual species when dispersal limitation exists.
Author Contributions
ZY originally formulated the idea, YC developed methodology, SN, PL, HJ, HW, and YY conducted fieldwork, YC performed statistical analyses and wrote the manuscript.
Funding
This work was supported by Scientific and technological projects in Henan Province (132102110133).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgment
The insightful comments of two reviewers were very helpful in improving the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00874/full#supplementary-material
Footnotes
References
Ah-Peng, C., Wilding, N., Kluge, J., Descamps-Julien, B., Bardat, J., Chuah-Petiot, M., et al. (2012). Bryophyte diversity and range size distribution along two altitudinal gradients: continent vs. island. Acta Oecol. 42, 58–65. doi: 10.1016/j.actao.2012.04.010
Andersen, R., Poulin, M., Borcard, D., Laiho, R., Laine, J., Vasander, H., et al. (2011). Environmental control and spatial structures in peatland vegetation. J. Veg. Sci. 22, 878–890. doi: 10.1111/j.1654-1103.2011.01295.x
Anderson, M. J., Crist, T. O., Chase, J. M., Vellend, M., Inouye, B. D., Freestone, A. L., et al. (2011). Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecol. Lett. 14, 19–28. doi: 10.1111/j.1461-0248.2010.01552.x
Andrew, N. R., Rodgerson, L., and Dunlop, M. (2003). Variation in invertebrate–bryophyte community structure at different spatial scales along altitudinal gradients. J. Biogeogr. 30, 731–746. doi: 10.1046/j.1365-2699.2003.00849.x
Aude, E., and Poulsen, R. S. (2000). Influence of management on the species composition of epiphytic cryptogams in Danish Fagus forests. Appl. Veg. Sci. 3, 81–88. doi: 10.2307/1478921
Barot, S. (2004). Mechanisms promoting plant coexistence: can all the proposed processes be reconciled? Oikos 106, 185–192. doi: 10.1111/j.0030-1299.2004.13038.x
Benítez,Á., Prieto, M., and Aragón, G. (2015). Large trees and dense canopies: key factors for maintaining high epiphytic diversity on trunk bases (bryophytes and lichens) in tropical montane forests. Forestry 88, 521–527. doi: 10.1093/forestry/cpv022
Blanchet, F. G., Legendre, P., and Borcard, D. (2008). Forward selection of explanatory variables. Ecology 89, 2623–2632. doi: 10.1890/07-0986.1
Borcard, D., and Legendre, P. (2002). All-scale spatial analysis of ecological data by means of principal coordinates of neighbour matrices. Ecol. Model. 153, 51–68. doi: 10.1016/S0304-3800(01)00501-4
Brunialti, G., Frati, L., Aleffi, M., Marignani, M., Rosati, L., Burrascano, S., et al. (2010). Lichens and bryophytes as indicators of old-growth features in Mediterranean forests. Plant Biosyst. 144, 221–233. doi: 10.1080/11263500903560959
Bruun, H. H., Moen, J., Virtanen, R., Grytnes, J. A., Oksanen, L., and Angerbjörn, A. (2006). Effects of altitude and topography on species richness of vascular plants, bryophytes and lichens in alpine communities. J. Veg. Sci. 17, 37–46. doi: 10.1111/j.1654-1103.2006.tb02421.x
Chang, L. W., Zelený, D., Li, C. F., Chiu, S. T., and Hsieh, C. F. (2013). Better environmental data may reverse conclusions about niche-and dispersal-based processes in community assembly. Ecology 94, 2145–2151. doi: 10.1890/12-2053.1
Chen, S., Slik, J. F., Mao, L., Zhang, J., Sa, R., Zhou, K., et al. (2015). Spatial patterns and environmental correlates of bryophyte richness: sampling effort matters. Biodivers. Conserv. 24, 593–607. doi: 10.1007/s10531-014-0838-8
Chen, Y., Yuan, Z., Li, P., Cao, R., Jia, H., and Ye, Y. (2016). Effects of environment and space on species turnover of woody plants across multiple forest dynamic plots in East Asia. Front. Plant Sci. 7:1533. doi: 10.3389/fpls.2016.01533
Chiarucci, A., D’auria, F., and Bonini, I. (2007). Is vascular plant species diversity a predictor of bryophyte species diversity in Mediterranean forests? Biodivers. Conserv. 16, 525–545. doi: 10.1007/s10531-006-6733-1
Conradi, T., Temperton, V. M., and Kollmann, J. (2017). Resource availability determines the importance of niche-based versus stochastic community assembly in grasslands. Oikos doi: 10.1111/oik.03969
Coyle, J. R. (2017). Intraspecific variation in epiphyte functional traits reveals limited effects of microclimate on community assembly in temperate deciduous oak canopies. Oikos 126, 111–120. doi: 10.1111/oik.03239
Friedel, A., Dengler, J., and Härdtle, W. (2006). Species diversity and species composition of epiphytic bryophytes and lichens–a comparison of managed and unmanaged beech forests in NE Germany. Feddes Repert. 117, 172–185. doi: 10.1002/fedr.200511084
Gradstein, S. R., Van Reenen, G. B. A., and Griffin, D. (1989). Species richness and origin of the bryophyte flora of the Colombian Andes. Plant Biol. 38, 439–448. doi: 10.1111/j.1438-8677.1989.tb01375.x
Griffith, D. A., and Peres-Neto, P. R. (2006). Spatial modeling in ecology: the flexibility of eigenfunction spatial analyses. Ecology 87, 2603–2613. doi: 10.1890/0012-9658(2006)87[2603:SMIETF]2.0.CO;2
Halpern, C. B., Dovèiak, M., Urgenson, L. S., and Evans, S. A. (2014). Substrates mediate responses of forest bryophytes to a gradient in overstory retention. Can. J. For. Res. 44, 855–866. doi: 10.1139/cjfr-2014-0059
Heilmann-Clausen, J., Aude, E., Dort, K., Christensen, M., Piltaver, A., Veerkamp, M., et al. (2014). Communities of wood-inhabiting bryophytes and fungi on dead beech logs in Europe–reflecting substrate quality or shaped by climate and forest conditions? J. Biogeogr. 41, 2269–2282. doi: 10.1111/jbi.12388
Hidalgo, J., Suweis, S., and Maritan, A. (2017). Species coexistence in a neutral dynamics with environmental noise. J. Theor. Biol. 413, 1–10. doi: 10.1016/j.jtbi.2016.11.002
Hofmeister, J., Hošek, J., Holá, E., and Novozámská, E. (2015). Decline in bryophyte diversity in predominant types of central European managed forests. Biodivers. Conserv. 24, 1391–1402. doi: 10.1007/s10531-015-0863-2
Hosokawa, T., and Odani, N. (1957). The daily compensation period and vertical ranges of epiphytes in a beech forest. J. Ecol. 45, 901–915. doi: 10.2307/2256963
Hubbell, S. P. (2001). The Unified Neutral Theory of Biodiversity and Biogeography. Princeton, NJ: Princeton University Press.
Humphrey, J. W., Davey, S., Peace, A. J., Ferris, R., and Harding, K. (2002). Lichens and bryophyte communities of planted and semi-natural forests in Britain: the influence of site type, stand structure and deadwood. Biol. Conserv. 107, 165–180. doi: 10.1016/S0006-3207(02)00057-5
Hunter, P. R., and Gaston, M. A. (1988). Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J. Clin. Microbiol. 26, 2465–2466.
Hutchinson, G. E. (1957). Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 22, 415–427. doi: 10.1101/SQB.1957.022.01.039
Jia, H. R., Chen, Y., Yuan, Z. L., Ye, Y. Z., and Huang, Q. C. (2015). Effects of environmental and spatial heterogeneity on tree community assembly in Baotianman National Nature Reserve. Henan, China. Pol. J. Ecol. 63, 175–183. doi: 10.3161/15052249pje2015.63.2.002
Lee, T. D., and Roi, G. H. (1979). Bryophyte and understory vascular plant beta diversity in relation to moisture and elevation gradients. Plant Ecol. 40, 29–38. doi: 10.1007/BF00052012
Legendre, P., and Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280. doi: 10.1007/s004420100716
Legendre, P., and Legendre, L. (1998). Numerical Ecology: Developments in Environmental Modelling, Vol. 24. Amsterdam: Elsevier Science.
Legendre, P., Mi, X., Ren, H., Ma, K., Yu, M., Sun, I., et al. (2009). Partitioning beta diversity in a subtropical broad-leaved forest of China. Ecology 90, 663–674. doi: 10.1890/07-1880.1
Lesica, P., McCune, B., Cooper, S. V., and Hong, W. S. (1991). Differences in lichen and bryophyte communities between old-growth and managed second-growth forests in the Swan Valley, Montana. Can. J. Bot. 69, 1745–1755. doi: 10.1139/b91-222
Lieberman, D., Lieberman, M., Peralta, R., and Hartshorn, G. S. (1996). Tropical forest structure and composition on a large-scale altitudinal gradient in Costa Rica. J. Ecol. 84, 137–152. doi: 10.2307/2261350
Lindenmayer, D. B., and Laurance, W. F. (2012). A history of hubris–Cautionary lessons in ecologically sustainable forest management. Biol. Conserv. 151, 11–16. doi: 10.1016/j.biocon.2011.10.032
Löbel, S., Snäll, T., and Rydin, H. (2006). Metapopulation processes in epiphytes inferred from patterns of regional distribution and local abundance in fragmented forest landscapes. J. Ecol. 94, 856–868. doi: 10.1111/j.1365-2745.2006.01114.x
Martens, S. N., Ustin, S. L., and Norman, J. M. (1991). Measurement of tree canopy architecture. Int. J. Remote Sens. 12, 1525–1545. doi: 10.1080/01431169108955187
Mežaka, A., Brûmelis, G., and Piterâns, A. (2012). Tree and stand-scale factors affecting richness and composition of epiphytic bryophytes and lichens in deciduous woodland key habitats. Biodivers. Conserv. 21, 3221–3241. doi: 10.1007/s10531-012-0361-8
Miles, C. J., and Longton, R. E. (1992). Deposition of moss spores in relation to distance from parent gametophytes. J. Bryol. 17, 355–368. doi: 10.1179/jbr.1992.17.2.355
Mills, S. E., and Macdonald, S. E. (2005). Factors influencing bryophyte assemblage at different scales in the western Canadian boreal forest. Bryologist 108, 86–100. doi: 10.1639/0007-2745(2005)108[86:FIBAAD]2.0.CO;2
Moe, B., and Botnen, A. (2000). Epiphytic vegetation on pollarded trunks of Fraxinus excelsior in four different habitats at Grinde, Leikanger, western Norway. Plant Ecol. 151, 143–159. doi: 10.1023/A:1026585911823
Mota de Oliveira, S., Ter Steege, H., Cornelissen, J. H., and Robbert Gradstein, S. (2009). Niche assembly of epiphytic bryophyte communities in the Guianas: a regional approach. J. Biogeogr. 36, 2076–2084. doi: 10.1111/j.1365-2699.2009.02144.x
Nascimbene, J., Fontana, V., and Spitale, D. (2014). A multi-taxon approach reveals the effect of management intensity on biodiversity in Alpine larch grasslands. Sci. Total Environ. 487, 110–116. doi: 10.1016/j.scitotenv.2014.04.013
Nathan, R., and Muller-Landau, H. C. (2000). Spatial patterns of seed dispersal, their determinants and consequences for recruitment. Trends Ecol. Evol. 15, 278–285. doi: 10.1016/S0169-5347(00)01874-7
Ódor, P., Király, I., Tinya, F., Bortignon, F., and Nascimbene, J. (2013). Patterns and drivers of species composition of epiphytic bryophytes and lichens in managed temperate forests. For. Ecol. Manag. 306, 256–265. doi: 10.1016/j.foreco.2013.07.001
Oksanen, J., Kindt, R., Legendre, P., O’Hara, B., Stevens, M. H. H., Oksanen, M. J., et al. (2007). The vegan package. Community Ecol. Package 10, 631–637.
Oliveira, S. M., and Steege, H. (2015). Bryophyte communities in the Amazon forest are regulated by height on the host tree and site elevation. J. Ecol. 103, 441–450. doi: 10.1111/1365-2745.12359
Pharo, E. J., and Beattie, A. J. (2002). The association between substrate variability and bryophyte and lichen diversity in eastern Australian forests. Bryologist 105, 11–26. doi: 10.1639/0007-2745(2002)105[0011:TABSVA]2.0.CO;2
Pharo, E. J., Beattie, A. J., and Binns, D. (1999). Vascular plant diversity as a surrogate for bryophyte and lichen diversity. Conserv. Biol. 13, 282–292. doi: 10.1046/j.1523-1739.1999.013002282.x
Pharo, E. J., Beattie, A. J., and Pressey, R. L. (2000). Effectiveness of using vascular plants to select reserves for bryophytes and lichens. Biol. Conserv. 96, 371–378. doi: 10.1016/S0006-3207(00)00080-X
Rambo, T. R., and Muir, P. S. (1998). Forest floor bryophytes of Pseudotsuga menziesii-Tsuga heterophylla stands in Oregon: influences of substrate and overstory. Bryologist 101, 116–130. doi: 10.1639/0007-2745(1998)101[116:FFBOPM]2.0.CO;2
Ranius, T., Johansson, P., Berg, N., and Niklasson, M. (2008). The influence of tree age and microhabitat quality on the occurrence of crustose lichens associated with old oaks. J. Veg. Sci. 19, 653–662. doi: 10.3170/2008-8-18433
Rosabal, D., Burgaz, A. R., and Reyes, O. J. (2013). Substrate preferences and phorophyte specificity of corticolous lichens on five tree species of the montane rainforest of Gran Piedra, Santiago de Cuba. Bryologist 116, 113–121. doi: 10.1639/0007-2745-116.2.113
Scott, G. A. M. (1994). Elementary reflections on the biology of bryophytes. Vic. Nat. 111, 112–115.
Silva, J. B., Santos, N. D. D., and Pôrto, K. C. (2014). Beta-diversity: effect of geographical distance and environmental gradients on the rocky outcrop bryophytes. Cryptogam. Bryol. 35, 133–163. doi: 10.7872/cryb.v35.iss2.2014.133
Smith, R. J., and Stark, L. R. (2014). Habitat vs. dispersal constraints on bryophyte diversity in the Mojave Desert. USA. J. Arid Environ. 102, 76–81. doi: 10.1016/j.jaridenv.2013.11.011
Söderström, L. (1988). The occurrence of epixylic bryophyte and lichen species in an old natural and a managed forest stand in northeast Sweden. Biol. Conserv. 45, 169–178. doi: 10.1016/0006-3207(88)90137-1
Soderstrom, L. (1993). Substrate preference in some forest bryophytes: a quantitative study. Lindbergia 18, 98–103.
Song, L., Ma, W. Z., Yao, Y. L., Liu, W. Y., Li, S., Chen, K., et al. (2015). Bole bryophyte diversity and distribution patterns along three altitudinal gradients in Yunnan. China. J. Veg. Sci. 26, 576–587. doi: 10.1111/jvs.12263
Spitale, D. (2016). The interaction between elevational gradient and substratum reveals how bryophytes respond to the climate. J. Veg. Sci. 27, 844–853. doi: 10.1111/jvs.12403
Spitale, D. (2017). Forest and substrate type drive bryophyte distribution in the Alps. J. Bryol. 39, 128–140. doi: 10.1080/03736687.2016.1274090
Stehn, S. E., Webster, C. R., Glime, J. M., and Jenkins, M. A. (2010). Elevational gradients of bryophyte diversity, life forms, and community assemblage in the southern Appalachian Mountains. Can. J. For. Res. 40, 2164–2174. doi: 10.1139/X10-156
Sun, S. Q., Wu, Y. H., Wang, G. X., Zhou, J., Yu, D., Bing, H. J., et al. (2013). Bryophyte species richness and composition along an altitudinal gradient in Gongga Mountain. China. PLoS ONE 8:e58131. doi: 10.1371/journal.pone.0058131
Szövényi, P., Hock, Z. S., and Tóth, Z. (2004). Phorophyte preferences of epiphytic bryophytes in a stream valley in the Carpathian Basin. J. Bryol. 26, 137–146. doi: 10.1179/037366804225021092
Ter Braak, C. J. (1987). “The analysis of vegetation-environment relationships by canonical correspondence analysis,” in Theory and Models in Vegetation Science, eds I. C. Prentice and E. Van der MaarelIn (Dordrecht: Springer). doi: 10.1007/978-94-009-4061-1_7
Turner, P. A., and Pharo, E. J. (2005). Influence of substrate type and forest age on bryophyte species distribution in Tasmanian mixed forest. Bryologist 108, 67–85. doi: 10.1639/0007-2745(2005)108[67:IOSTAF]2.0.CO;2
Tusiime, F. M., Byarujali, S. M., and Bates, J. W. (2007). Diversity and distribution of bryophytes in three forest types of Bwindi Impenetrable National Park, Uganda. Afr. J. Ecol. 45, 79–87. doi: 10.1111/j.1365-2028.2007.00862.x
Valladares, F., Bastias, C. C., Godoy, O., Granda, E., and Escudero, A. (2015). Species coexistence in a changing world. Front. Plant Sci. 6:866. doi: 10.3389/fpls.2015.00866
Vanderpoorten, A., Engels, P., and Sotiaux, A. (2004). Trends in diversity and abundance of obligate epiphytic bryophytes in a highly managed landscape. Ecography 27, 567–576. doi: 10.1111/j.0906-7590.2004.03890.x
Violle, C., Enquist, B. J., McGill, B. J., Jiang, L., Albert, C. H., Hulshof, C., et al. (2012). The return of the variance: intraspecific variability in community ecology. Trends Ecol. Evol. 27, 244–252. doi: 10.1016/j.tree.2011.11.014
Vitt, D. H., Li, Y., and Belland, R. J. (1995). Patterns of bryophyte diversity in peatlands of continental western Canada. Bryologist 98, 218–227. doi: 10.2307/3243306
Wolf, J. H. (1993). Diversity patterns and biomass of epiphytic bryophytes and lichens along an altitudinal gradient in the northern Andes. Ann. Mo. Bot. Gard. 80, 928–960. doi: 10.2307/2399938
Ye, Y. Z., Wang, W. S., and Li, H. S. (2004). Scientific Survey of the Xiaoqinling Nature Reserve in Henan Province. Beijing: Science Press.
Keywords: bryophytes, elevation gradient, species diversity, species protection, variation partitioning
Citation: Chen Y, Niu S, Li P, Jia H, Wang H, Ye Y and Yuan Z (2017) Stand Structure and Substrate Diversity as Two Major Drivers for Bryophyte Distribution in a Temperate Montane Ecosystem. Front. Plant Sci. 8:874. doi: 10.3389/fpls.2017.00874
Received: 04 January 2017; Accepted: 10 May 2017;
Published: 26 May 2017.
Edited by:
Sebastian Leuzinger, Auckland University of Technology, New ZealandReviewed by:
Nagore G. Medina, Universidad Autonoma de Madrid, SpainJoan Bruno Silva, Federal University of Pernambuco, Brazil
Copyright © 2017 Chen, Niu, Li, Jia, Wang, Ye and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhiliang Yuan, zlyuan163@163.com