Next Article in Journal
Hypocrealean Fungi Associated with Hylobius abietis in Slovakia, Their Virulence against Weevil Adults and Effect on Feeding Damage in Laboratory
Next Article in Special Issue
On the Management of Large-Diameter Trees in China’s Forests
Previous Article in Journal
Selection of Suitable Reference Genes in Pinus massoniana Lamb. Under Different Abiotic Stresses for qPCR Normalization
Previous Article in Special Issue
Forest Structure, Wood Standing Stock, and Tree Biomass in Different Restoration Systems in the Brazilian Atlantic Forest
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Forest Structure and Composition along the Altitudinal Gradient in the State of Sikkim, Eastern Himalayas, India

by
Yangchenla Bhutia
1,2,*,
Ravikanth Gudasalamani
1,
Rengaian Ganesan
1 and
Somidh Saha
3,4,*
1
Suri Sehgal Centre for Biodiversity and Conservation, Ashoka Trust for Research in Ecology and the Environment, Royal Enclave, Sriramapura, Jakkur Post, Bengaluru 560064, India
2
Manipal Academy of Higher Education (MAHE), Tiger Circle Road, Madhav Nagar, Manipal, Karnataka 576104, India
3
Institute for Technology Assessment and Systems Analysis, Karlsruhe Institute of Technology, Karlstr. 11, D-76133 Karlsruhe, Germany
4
Chair of Silviculture, Institute of Forest Sciences, University of Freiburg, Tennenbacherstr. 4, D-79085 Freiburg, Germany
*
Authors to whom correspondence should be addressed.
Forests 2019, 10(8), 633; https://doi.org/10.3390/f10080633
Submission received: 9 June 2019 / Revised: 21 July 2019 / Accepted: 24 July 2019 / Published: 27 July 2019
(This article belongs to the Special Issue Forest Structure and Sustainable Resource Management)

Abstract

:
Understanding the structure and composition of native forests is a prerequisite in developing an adaptive forest management plan for Himalayan forest ecosystems where climate change is rapid. However, basic information on forest structure and composition are still limited in many places of the Eastern Himalayas. In this study, we aimed to understand the diversity, structure, and composition of forests and their variations along an altitudinal gradient in Himalayan forests. The study was conducted in the Indian federal state of Sikkim, Eastern Himalayas. We carried out a comprehensive and comparative evaluation of species diversity, stand basal area, and stem density along the altitudinal gradient from 900 m a.s.l. to 3200 m a.s.l. We used stratified random sampling to survey eighty-three plots each 0.1 ha in forest communities that occurred along the altitudinal gradient: (a) lower (900–1700 m) altitude forest (N = 24), (b) mid (1700–2500 m) altitude forests (N = 37), and (c) higher (2500–3200 m)altitude forests (N = 22). We measured and identified all living trees with a >3 cm diameter at breast height in each plot. We counted 10,344 individual plants, representing 114 woody species belonging to 42 families and 75 genera. The family Fagaceae and its species Lithocarpus pachyphyllus (Kurz) Rehder. were reported as the most dominant forest trees with the highest Importance Value Index. The Shannon diversity index was recorded as being the highest for the low-altitude forests, whereas measures of structural diversity varied among forests along with altitude: the mid-altitude forests recorded the highest stem density and the high-altitude forests showed the highest mean stem DBH (diameter at 1.3 m height). One significant finding of our study was the disparity of the size class distribution among forests along the altitudinal gradient. Overall, we found a reverse J-shape distribution of tree diameter signifying the uneven-agedness. However, we showed, for the first time, a complete lack of large trees (>93 cm DBH) in the lower altitude forests. Our study highlights conservation concerns for the low-altitude forests that record high species diversity, although lacked large-diameter trees. We anticipate that our study will provide a comprehensive understanding of forest diversity, composition, and structure along the altitudinal gradient to design conservation and sustainable management strategies

1. Introduction

1.1. Background

Quantitative assessment of the diversity, structure, and composition of forest trees are essential, not only to understand forest biodiversity and health but also for designing conservation strategies towards climate change [1], which is a significant threat to biodiversity in the Himalayas [2,3,4,5]. The Eastern Himalaya region is a confluence of multiple biogeographic origins such as the Indo–Malayan realm, Palearctic, and the Sino–Japanese region, and it is characterized by remarkable biodiversity of major ecological and global significance and represents one of the 34 global biodiversity hotspots [6,7]. It covers an extensive area of 524,190 sqkm, ranging from east of Nepal to northwest Yunnan Province in China, including the northeastern region of India [2]. The Himalayan forests stores abundant carbon in forest soil and vegetation, and thus, are of prime importance for the regional and global carbon cycles [8,9]. It supplies constant water and is the prime regulator of life [10]. Further, the forest-based biological resources are the chief source of livelihood for people living in the Himalayas [11]. Simultaneously, the Himalayas are one of the most fragile systems on earth [12], geologically young (45 million years), tectonically active, and highly susceptible to natural hazards [4]. Moreover, the Himalayas are reported as the prime site for changing climate [13], with forest vegetation notably vulnerable to climate change [14,15,16]. The degrading Himalayan forests not only affect the surrounding region but have global implications [17]. While the ecological fragility and the significance of Himalayan forests are much realized, basic knowledge on the structure and composition in Himalayan forests are still limited in many regions particularly in the remote eastern parts [16].
Forest trees are the dominant structural and functional component of the forest ecosystem [18]. Tree species diversity, size class distribution, stem density, and basal area are the essential attributes that describes a forest’s ecosystem [19,20], and measuring these attributes are fundamental in designing conservation strategies. Species diversity is the most crucial descriptors, which not only captures information on species richness (number of species in a community) but the relative abundance of species in a forest. It also provides information on the rarity and the commonness of a species. The ability of species diversity to quantify forest composition has provided ecologist with the most prominent tool, often used as the scorecard to preserve and restore a forest. Likewise, the size class distribution of forest trees is the prime indicator of forest structure and dynamics, widely used to examine the forest’s health, including regeneration [21]. Tree diameter is a statistically proven parameter to measure forest carbon stocks [22]. Furthermore, stem density and basal area are an excellent surrogate to estimate forest biomass and carbon [23,24]. A proper understanding of forest composition and structure allow foresters, national park rangers, and landowners to maximize forest ecosystem’s goods and services by maintaining or conserving a desired structure and composition of the forest at stand and landscape level.

1.2. Past Studies in the Eastern Himalayas

Studies in the Eastern Himalayan forests outside of India were mainly carried out in Nepal [25,26], Bhutan [27,28,29], and China [30,31] to assess forest structure and composition. Within the Indian part of the Eastern Himalayas, studies were carried out in Arunachal Pradesh [32,33,34], Meghalaya [8,35,36], Darjeeling, West Bengal [37], and Assam [38,39]. The majority of these studies presented the structure and composition of different forests along altitudinal gradients with varying intensities of anthropogenic disturbance. However, we did not find similar studies from the Sikkim State of the Eastern Himalaya region in India. Studies on the forest of Sikkim have mostly been carried out as case studies with concentrated sampling on watershed [40], alpine forests [41,42], timberline forest [43], and trekking corridors [44], with narrow altitudinal range [42,43]. Albeit more recently, Acharya et al. [45] and Sharma et al. [46] performed research with samples from multiple locations but did not describe the composition and size class distribution of trees.

1.3. Aims of this Study

In this study, we aimed to fill the existing knowledge gap on the diversity, structure, and composition of forest tree species, using well-represented samples covering forests along the altitudinal gradient in Sikkim, a part of the Eastern Himalayas. Our study not only provides an overall assessment on the diversity, composition, and structure of forest tree communities but also evaluates how these attributes vary among forest communities which occur along the altitudinal gradient.

2. Materials and Methods

2.1. Study Area Description and Data Collection

Sikkim (27°07′ to 28°13′ N and 88°01′ to 88°92′ E), a northeastern region of India, falls within the Eastern Himalayan biodiversity hotspot (Figure 1). It has a remarkable altitudinal gradient with only an area of 7096 km2, ranging from 300 m a.s.l to 8586 m a.s.l. with Mount Kanchendzonga—the third highest peak in the world. The complex topography and enormous altitude are characterized by 12 major vegetation types, from a tropical warm broad-leaved forest at the lowest altitude to alpine meadow at the highest altitude [47]. It has a subtropical to temperate climate, with annual rainfall between 2700 mm to 3200 mm, and the mean annual temperature (MAT) varies from 8.4 °C to 23.2 °C [48].
The field study was conducted over 3 years (2013–2015), during the dry winter seasons from the beginning of October to late May to avoid challenges caused by understory growth of thick shrubs and herbs during the monsoon season. Stratified random sampling was used to sample plots along the altitudinal gradients, from 900 m a.s.l. to 3200 m a.s.l. The forests were selected because, forests within the range has undergone significant loss in forest cover [47]. Thus, quantitative assessment of these forests is imperative to plan conservation strategies. The forests were stratified into three altitudinal zones: low-altitude forest (LF), ranging from 900 m a.s.l to 1700 m a.s.l., the mid-altitude forest (MF), ranging from 1700 m a.s.l. to 2500 m a.s.l., and the high-altitude forest (HF) ranging from 2500 m a.s.l to 3200 m a.s.l. Plots were randomly installed at different locations for each altitude zone. Altogether, we installed 83 plots of 100 m × 10 m (each 0.1 ha), covering a total area of 8.3 ha. The number of plots varied across forests, and were distributed as LF = 24, MF = 37, and HF = 22 plots among three forest communities along the altitudinal gradient. Within each plot, all living woody trees with >3 cm diameter at breast height (DBH, or trunk diameter at 1.3 m above ground level) were measured and identified. Voucher specimens were collected for individuals challenging to identify in the field, and was used later for further identification.

2.2. Data Analysis

The analysis included quantitative data on the abundance, DBH, taxonomic, and geo-coordinate details for all sites. Prior to the quantification of forest attributes, we conducted cluster analysis and ordination to examine and validate the presence of ecologically meaningful clusters among our sampled plots designed for the study. First, we computed the Bray–Curtis distance matrix using the abundance data, with species in the columns and sites in ascending order of altitude in rows. Then, we performed hierarchical Ward’s minimum variance clustering on the Bray–Curtis distance matrix. This method is based on the least squares linear model criteria, and forms cluster with the least variance [49]. Secondly, we performed the Nonmetirc multidimensional scaling method (NMDS), an indirect gradient analysis approach to ordinate our sampled sites using the same distance matrix. Nonmetric multidimensional scaling represents the pairwise dissimilarity among sites in a low dimensional space as closely as possible [49]. We also estimated the goodness of fit measured as “stress” from the Shepard diagram to evaluate the appropriateness of the NMDS result [49].
We plotted the species accumulation curve to determine the sampling effectiveness for the entire sample installed in our study. Further, we plotted sample-based rarefaction curve for different forest communities. It computes the average richness by randomly drawing 1,2,3,....n samples representing a community, without replacement [50]. The richness estimated from the sample-based rarefaction curves was also used to compare the richness for different forest communities with distinct sampling effort.
Shannon’s diversity index (H’) was estimated by multiplying the proportion of each species to their natural log. Shannon Index (H’) = Σpilog(ln)pi, where, pi is the proportion (n/N) of individuals of a particular species found (n) divided by the total number of individuals recorded (N), ln is the natural log, and Σ is the sum of the calculations.
Shannon diversity provides information on both species richness and relative abundance among plots, and thus are sensitive to the sampling effort or the number of individuals sampled. To avoid bias, we estimated other indices, namely, the rarefied richness for forest communities and compared richness among subsets using a fixed sample size. Further, we estimated the Fisher alpha, a parametric richness estimate, which compares richness among samples with varying individuals. Pielou’s J was also estimated to check for differences in evenness among communities with different sample size. Apart from these, we estimated observed species richness, genus richness, and family richness for the overall forest community and three forest communities along the altitudinal gradient.
Additionally, we estimated four essential stand structural parameters: the median DBH, stem density (ha−1), basal area (m2 ha−1), and size class distribution of forest trees. We mostly used the median as a measure of central tendency because variables often did not follow a Gaussian distribution. We calculated the Important Values Index (IVI) both for species(SIVI) and families(FIVI). We followed the equations by Keel et al. [51] and Ganesh et al. [52] for IVI estimation. Diversity analysis, clustering, and ordination were conducted using the R package “vegan”and “BiodiversityR”V 2.9-2, and the built-in function “decostand”, “vegdist”, “hclust”, “ward.D2”and “metaNMDS”in R Studio (version 1.1383).

3. Results

3.1. Distinct Forest Communities along the Altitude Gradient from 900 m to 3200 m

The cluster analysis revealed that three forest communities could be formed based on the composition of tree species at different sites, although few plots deflected from its designated categories (Figure 2a). The NMDS result superimposed with a cluster dendrogram also showed similar groups (Figure 2b). The Shepard diagram and the goodness of fit also supported our result with high R2 value (Figure 3). It can be stated that both the cluster analysis and NMDS suggested the presence of three ecologically meaningful forest communities along the altitudinal gradient with distinct species composition. These three forest communities corresponded to the altitudinal gradient, hence, reiterated the importance of altitude as a covariate to determine the composition of tree species in the Himalayas. Based on the Species Importance Value Index and a proportion of stem density and basal area (see Table 3), all forests can be classified as mixed-broadleaf forests but with distinct communities at different altitudes. The lower altitude (LF) forests can be classified as the association of Schima wallichii Choisy., Castanopsis tribuloides Smith., Engelhardtia spicata Blume., and Gynocardia odorata R.Br. (Figure 4). The mid-altitude (MF) forests can be classified asthe association of Symplocos ramosissima Wall. ex. G. Don, Castanopsis hystrix Hook.f. & Thomson ex A.DC., Viburnum erubescens Wall., and Quercus lamellosa Sm. (Figure 5). Whereas, high-altitude (HF) forests can be classified as the association of Lithocarpus pachyphyllus, Rhododendron arboreum Sm., Magnolia campbellii Hook.f.& Thomson, and Symplocos heishanensis Hayata. (Figure 6).

3.2. Diversity and Composition of Tree Species in Three Forest Communities

Overall, the study recorded 10,344 individuals from 8.3 ha (83 sample plots of each 0.1 ha in size), including 87 unidentified individuals, with the observed richness of 114 species from 42 families and 75 genera (Table 1). The species accumulation curve for the entire tree community showed a gentle slope for the cumulative species richness after the sample size of 70 indicating an adequate number of plots (Figure 7a). However, for different forest communities, the species accumulation curve showed an increasing asymptote (Figure 7b). Nevertheless, the curve clearly recorded the highest species richness (77) for the MF, followed by LF (68) and HF(51). However, for the same sample size of 20, the LF recorded the highest average species richness of 65.36 (Figure 7b). The Shannon (H) diversity varied from 1.71 ± 0.40 to 2.10 ± 0.35, with the highest value for the LF (Table 1), and the Pielou’s evenness (J) varied from 0.46 ± 0.7 to 0.35 ± 0.10, with the highest value for the LF (Table 1).
Table 2 and Table 3 show the results for the Importance Value Index (IVI), estimated for both species (SIVI) and families (FIVI), along with species richness, families, stem density, basal area, their proportion to the contribution for the overall forest community, and different forest communities. We found Lithocarpus pachyphyllus (9.30%), Symplocos ramosissima (8.72%), Castanopsis hystrix (8.23%), Castanopsis tribuloides (4.99%), and Schima wallichii (4.46%) to be the five most important species with the highest Species Importance Value Index (SIVI). Similarly, we recorded Fagaceae (27.86%), Symplocaceae (12.26%), Theaceae (8.04%), Ericaceae (7.70%), and Lauraceae (6.62%), as the top five families. Lauraceae was also recorded as the most diverse family with 13 species (Table 2).
Species such as Schima wallichii (18.90%), Castanopsis tribuloides (17.91%), Engelhardtia spicata (6.72%), Gynocardia odorata (6.05%), and Eurya acuminata (4.83%) were recorded as the five most dominant species with higher species Important Value Index (SIVI) for LF. Similarly, Symplocos ramosissima (15.09%), Castanopsis hystrix (14.99%), Viburnum erubescens (7.09%), Quercus lamellosa (6.85%) and Eurya acuminata (4.42%) were reported as the most dominant species for the MF. The HF was dominated by species such as Lithocarpus pachyphyllus (24.37%), Rhododendron arboreum (16.78%), Magnolia campbellii (4.54%), Symplocos heishanensis (4.14%), and Symplocos ramosissima (3.96%) with the highest SIVI (Table 3).

3.3. Forest Structure

Overall, we recorded the median stem DBH of 15.42 ± 1.11 cm, stem density of 1246.27 ± 568.28 ha−1, and the basal area of 54.35 ± 4.44 m2 ha−1 (Figure 8). Among forest communities, the HF recorded the highest median stem DBH (20.89 ± 3.09 cm) but showed the lowest stem density. Similarly, MF recorded the highest stem density of about 1320 ± 119.72 ha−1, with the second highest value for the stem DBH and basal area. However, the LF recorded the second highest value for median stem density and the minimum value for the median DBH and basal area.
The overall size class distribution showed a reverse “J”-shaped pattern (Figure 9), with the highest number of individuals in the smallest DBH class between 3 to 13 cm (Figure 9a), and the least number of individuals in the largest DBH class between 293–303 cm (Figure 9b). A similar distribution was documented for different forest communities. However, they differed in the number of DBH classes: MF recorded the maximum number of DBH classes with the biggest tree of 295 cm DBH recorded in the study; LF and MF lacked comparable DBH trees. The LF recorded the least number of DBH classes and completely lacked trees above 93 cm DBH.

4. Discussion

4.1. Forest Communities and Species Composition

The cluster analysis highlighted the significance of altitude in driving compositional change of forest tree species recorded in the Himalayan forests, like Sikkim. The finding could be related to the close association of altitude with climatic variables: temperature and precipitation [53], and thus draws attention towards the susceptibility of Himalayan forests to climate change.
Additionally, cluster analysis also implied the presence of three distinct forest communities between 900 ma.s.l. and 3200 m a.s.l. Moreover, according to the Grierson and Long [54] classification, our study resembles three forest types, i.e., LF forest as the warm broad-leaved forest, MF as the evergreen oak forest, and HF as the temperate mixed forests [54]. Nevertheless, we recorded some differences in species recorded for high-elevation forests. For example, Grierson and Long’s [54] classification reported Tsuga dumosa (D.Don)Eichler., Larix graffithii Hook.f., and Picea smithiana (Wall.) Boiss. as the major trees for high-elevation forests; however, we did not report similar findings. A reason in the change in species composition in these three forest types compared to the 1960s and 1970s could be related to increasing anthropogenic pressure, such as open grazing, tree felling for its use as scaffolding, timber for construction of houses, firewood, and clearing for forest for agriculture practice, as well as climate change effects; however, warrants more investigation. Future research should focus on documentation of forest trees for a detailed description of the Himalayan forest and revision of forest type classifications if needed.

4.2. Tree Species Richness and Diversity

Species accumulation curves are especially useful to decide sampling completeness and also to compare richness for communities with different sample sizes [50]. In the current analysis, the accumulation curves did not record saturation when plotted for different communities, which implies finding more species with increasing sampling effort. Nevertheless, it is clear that low-altitude forests record the highest number of species richness and evenness, and the disparity between the observed and rarified richness in the MF was the mere result of the sampling size differences. The varied species richness between forests covering different altitudes can be associated with numerous factors. Broadly, we assumed the widely reported climatic variables, namely, precipitation, temperature, and their interaction, as the prime factor for varied richness along the altitudinal gradient of the Himalayan system [46,55,56]. Another crucial factor could be anthropogenic disturbance and its intensity [27,38]. Future research could explore these variables to explain the difference in richness recorded in our study. The LF recorded not only the highest number of species but recorded more uniformly distributed tree species. However, the region lacks comparable literature to discuss optimum richness and evenness. We recommend long-term monitoring of permanent plots, which is completely lacking in the Indian Eastern Himalayas.
One of the important findings of our study was the significance of the family Fagaceae. It was recorded as having the highest FIVI and SIVI and dominated the forest along the altitudinal gradient between 900 m and3200 m a.s.l. in Sikkim. The reasons should be further investigated because ecological and silvicultural knowledge of Fagaceae and its species are very limited for the Eastern Himalayan forests. Furthermore, we also showed for the first time the existence of huge trees of Lithocarpus pachyphylls (>1.5 m DBH) in high number, which implies that this species was the remnant of the old growth primary forests in the altitude from 2500 to 3200 in Sikkim Himalayas. These old-growth trees should be marked, registered, monitored, and conserved as habitat trees by the forest department.
However, not all forest communities were dominated by the Fagaceae species. For example, Schima wallichii and Castanopsis tribuloides dominated the lower altitude forests, whereas species like Symplocos ramosissima and Castanopsis hystrix dominated the mid-altitude forests. Different physiological and climatic adaptation of a species may have caused distinct species to dominate different altitude, resulting in three unique clusters of our samples. Moreover, the Himalayas are renowned for the pronounced variation in climatic, topographic and edaphic factors along the altitudinal gradient. Apart from these, the anthropogenic disturbance is a widely discussed factor for the dominance of fast-growing pioneer species like the Symplocos racemose [57].Future research could explore the role of anthropogenic disturbance and different climatic and topographic factors on species dominance and composition along the altitudinal gradient.

4.3. Forest Structure: Stem Diameter, Stem Density, and Basal Area

Size class distribution of trees provides the population structure of forest [58], and are extensively used to understand regeneration status [59]. The present assessment presented a reverse J-shaped distribution, suggesting uneven-aged forests for sustainable reproduction and regeneration [60], with a sufficient number of young individuals to replace the old mature stand. However, between forest types, the low-altitude forests completely lacked trees above 93 cm diameter, a scenario undesirable for a sustainable forest. The finding could be associated with forest degradation and deforestation reported by earlier studies [47,61]. One potential factor for the absence of large trees could be the practice of cardamom cultivation mainly in the low-altitude forest of Sikkim. Cardamom is a high-valued cash crop and it is intensively cultivated in the low- and mid-altitude forests, where forests are partially cleared for its cultivation [61]. Furthermore, the drying of cardamom demands a substantial supply of fuelwood, often generated from the nearby forests. Moreover, until 2004, Sikkim produced the highest among of large cardamom (Amomum subulatum Roxb.) in India, with the largest share in the world’s market [62], number of trees felled must been profound.

5. Conclusions

Our study provides a comprehensive and quantitative understanding of diversity, structure, and composition of forest trees along the altitudinal gradient, ranging from 900 m a.s.l. to 3200 m a.s.l. in the Sikkim Himalayas. The study estimated high-species richness in the low-altitude forest occurring at 900 m a.s.l. to 1700 m a.s.l., whereas the high-altitude forest covering 2500 m a.s.l. to 3200 m a.s.l. recorded the highest mean stem DBH and stand basal area. Furthermore, we found that Fagaceae trees are the most significant, with a dominant contribution towards the total basal area and stem density. Hence, Fagaceae trees could prove as having a potential for carbon storage and climate change mitigation in the Himalayas. At the same time, large-sized Fagaceae trees may have very high biodiversity values which are yet to be quantified. One of the significant findings of our study was the disproportionate size class distribution within forests at different altitudes. The low-altitude forests completely lacked trees above 93 cm DBH. Overall, the study highlights the need to conserve low-altitude forests to maintain diversity and improve forest structure. We suggest further research to understand the factors associated with varied diversity, composition, and uneven forest structure recorded in the study.

Author Contributions

Y.B. conceptualized the research concept, methodology, collected field data, conducted all the statistical analyses, and wrote the original draft. S.S., R.G. (Ravikanth Gudasalamani), and R.G. (Rengaian Ganesan) supervised and reviewed the article. R.G. (Rengaian Ganesan) helped with plant identification. S.S. acquired funding for the Article Processing Charge.

Funding

The Department of Biotechnology, Government of India funded the fieldwork and data collection (Grant No. BT/01/NE/PS/NCBS/09). The work was partially funded by the National Mission for Himalayan Studies (NMHS) under the Ministry of Environment, Forest, and Climate Change, Government of India (NMHS/2015–16/HF11/11).The International Union of Forest Research Organizations, Vienna (IUFRO) awarded the Young Scientist Award to Y.B. for a short-term visit to the Thünen Institute, Germany. The Karlsruhe Institute of Technology, Germany funded the Article Processing Charge.

Acknowledgments

Y.B. acknowledges the support of the Department of Forest Environment and Wildlife Management Department, Government of Sikkim, including various officials of the concerned Sanctuary for permits. Further, Y.B. is grateful to Janga Bir Subba, Director (retd.), Khangchendzonga National Park for providing permits, field crew, and a forest guest house during field data collection. Y.B. also gratefully acknowledges the Department of Science and Technology for sharing satellite images during the initial phase of the fieldwork. Y.B. also thanks Robert John Chandran for acquiring the grant from the Department of Biotechnology, and guidance during fieldwork. Y.B. gratefully acknowledges IUFRO, for awarding her the IUFRO-EFI Young Scientists Initiative award of 2019. Lastly, Y.B. thanks Roshan Subba, field assistant, for the genuine enthusiasm and support in the field. S.S. acknowledges his colleagues in IUFRO’s task force on “Forest Restoration and Adaptation under Global Change” for their helpful suggestions on the scope of forest restoration research in the Himalayas.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Condit, R.; Sukumar, R.; Hubbell, S.P.; Foster, R.B. Predicting population trends from size distributions: A direct test in a tropical tree community. Am. Nat. 1998, 152, 495–509. [Google Scholar] [CrossRef] [PubMed]
  2. Sharma, E.; Tse-ring, K.; Chettri, N.; Shrestha, A.; Kathmandu, N. Biodiversity in the Himalayas–trends, perception and impacts of climate change. In Proceedings of the International Mountain Biodiversity Conference, Kathmandu, Nepal, 16–18 November 2008. [Google Scholar]
  3. Chaudhary, P.; Bawa, K.S. Local perceptions of climate change validated by scientific evidence in the Himalayas. Biol. Lett. 2011, 7, 767–770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Pandit, M.K. The Himalayas must be protected. Nat. News 2013, 501, 283. [Google Scholar] [CrossRef] [PubMed]
  5. Gautam, M.R.; Timilsina, G.R.; Acharya, K. Climate Change in the Himalayas: Current State of Knowledge; The World Bank: Washington, DC, USA, 2013. [Google Scholar]
  6. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; Da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  7. Mittermeier, R.A.; Myers, N.; Mittermeier, C.G.; Robles, G. Hotspots: Earth’s Biologically Richest and Most Endangered Terrestrial Ecoregions; CEMEX, SA, Agrupación Sierra Madre, SC: Mexico City, Mexico, 1999. [Google Scholar]
  8. Upadhaya, K. Structure and floristic composition of subtropical broad-leaved humid forest of Cherapunjee in Meghalaya, Northeast India. J. Biodivers. Manag. For. 2015, 4, 2. [Google Scholar] [CrossRef]
  9. Brown, S.; Sathaye, J.; Cannell, M.; Kauppi, P.E. Mitigation of carbon emissions to the atmosphere by forest management. Commonw. For. Rev. 1996, 75, 80–91. [Google Scholar]
  10. Haripriya, G. Biomass carbon of truncated diameter classes in Indian forests. For. Ecol. Manag. 2002, 168, 1–13. [Google Scholar] [CrossRef]
  11. Saha, D.; Sundriyal, R.C. Utilization of non-timber forest products in humid tropics: Implications for management and livelihood. For. Policy Econ. 2012, 14, 28–40. [Google Scholar] [CrossRef]
  12. Diaz, H.F.; Grosjean, M.; Graumlich, L. Climate variability and change in high elevation regions: Past, present and future. Clim. Chang. 2003, 59, 1–4. [Google Scholar] [CrossRef]
  13. Williams, J.W.; Jackson, S.T.; Kutzbach, J.E. Projected distributions of novel and disappearing climates by 2100 AD. Proc. Natl. Acad. Sci. USA 2007, 104, 5738–5742. [Google Scholar] [CrossRef] [Green Version]
  14. Shrestha, U.B.; Gautam, S.; Bawa, K.S. Widespread climate change in the Himalayas and associated changes in local ecosystems. PLoS ONE 2012, 7, e36741. [Google Scholar] [CrossRef] [PubMed]
  15. Krishnaswamy, J.; John, R.; Joseph, S. Consistent response of vegetation dynamics to recent climate change in tropical mountain regions. Glob. Chang. Biol. 2014, 20, 203–215. [Google Scholar] [CrossRef] [PubMed]
  16. Chakraborty, A.; Saha, S.; Sachdeva, K.; Joshi, P.K. Vulnerability of forests in the Himalayan region to climate change impacts and anthropogenic disturbances: A systematic review. Reg. Environ. Chang. 2018, 18, 1783–1799. [Google Scholar] [CrossRef]
  17. Xu, J.; Grumbine, R.E.; Shrestha, A.; Eriksson, M.; Yang, X.; Wang, Y.; Wilkes, A. The melting Himalayas: Cascading effects of climate change on water, biodiversity, and livelihoods. Conserv. Biol. 2009, 23, 520–530. [Google Scholar] [CrossRef]
  18. Khan, M.L.; Menon, S.; Bawa, K.S. Effectiveness of the protected area network in biodiversity conservation: A case-study of Meghalaya state. Biodivers. Conserv. 1997, 6, 853–868. [Google Scholar] [CrossRef]
  19. Shaheen, H.; Ullah, Z.; Khan, S.M.; Harper, D.M. Species composition and community structure of western Himalayan moist temperate forests in Kashmir. For. Ecol. Manag. 2012, 278, 138–145. [Google Scholar] [CrossRef]
  20. Peet, R.K. The measurement of species diversity. Annu. Rev. Ecol. Syst. 1974, 5, 285–307. [Google Scholar] [CrossRef]
  21. Khan, M.; Rai, J.; Tripathi, R. Population structure of some tree species in disturbed and protected subtropical forests of north-east India. Acta Oecolog. 1987, 8, 247–255. [Google Scholar]
  22. Brown, S. Measuring carbon in forests: Current status and future challenges. Environ. Pollut. 2002, 116, 363–372. [Google Scholar] [CrossRef]
  23. Chave, J.; Andalo, C.; Brown, S.; Cairns, M.A.; Chambers, J.; Eamus, D.; Fölster, H.; Fromard, F.; Higuchi, N.; Kira, T.; et al. Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 2005, 145, 87–99. [Google Scholar] [CrossRef]
  24. Forrester, D.I.; Tachauer, I.H.H.; Annighoefer, P.; Barbeito, I.; Pretzsch, H.; Ruiz-Peinado, R.; Stark, H.; Vacchiano, G.; Zlatanov, T.; Chakraborty, T.; et al. Generalized biomass and leaf area allometric equations for European tree species incorporating stand structure, tree age and climate. For. Ecol. Manag. 2017, 396, 160–175. [Google Scholar] [CrossRef]
  25. Pandey, K.P. Structure, Composition and Diversity of Forest along the Altitudinal Gradient in the Himalayas, Nepal. Appl. Ecol. Environ. Res. 2016, 14, 235–251. [Google Scholar] [CrossRef]
  26. Panthi, M.P.; Chaudhary, R.P.; Vetaas, O.R. Plant species richness and composition in a trans-Himalayan inner valley of Manang district, central Nepal. Himal. J. Sci. 2007, 4, 57–64. [Google Scholar] [CrossRef]
  27. Tenzin, J.; Hasenauer, H. Tree species composition and diversity in relation to anthropogenic disturbances in broad-leaved forests of Bhutan. Int. J. Biodivers. Sci. Ecosyst. Serv. Manag. 2016, 12, 274–290. [Google Scholar] [CrossRef]
  28. Mukhia, P.K.; Wangyal, J.T.; Gurung, D. Floristic composition and species diversity of the chirpine forest ecosystem, Lobesa, Western Bhutan. For. Nepal 2011, 1–4. [Google Scholar]
  29. Wangda, P.; Ohsawa, M. Structure and regeneration dynamics of dominant tree species along altitudinal gradient in a dry valley slopes of the Bhutan Himalaya. For. Ecol. Manag. 2006, 230, 136–150. [Google Scholar] [CrossRef]
  30. Tang, C.Q.; Li, Y.-H.; Zhang, Z.-Y. Species Diversity Patterns in Natural Secondary Plant Communities and Man-made Forests in a Subtropical Mountainous Karst Area, Yunnan, SW China. Mt. Res. Dev. 2010, 30, 244–251. [Google Scholar] [CrossRef]
  31. Li, R.; Dao, Z.; Li, H. Seed Plant Species Diversity and Conservation in the Northern Gaoligong Mountains in Western Yunnan, China. Mt. Res. Dev. 2011, 31, 160–165. [Google Scholar] [CrossRef]
  32. Bharali, S.; Paul, A.; Khan, M.L.; Singha, L.B. Species diversity and community structure of a temperate mixed Rhododendron forest along an altitudinal gradient in West Siang District of Arunachal Pradesh, India. Nat. Sci. 2011, 9, 125–140. [Google Scholar]
  33. Nath, P.; Arunachalam, A.; Khan, M.; Arunachalam, K.; Barbhuiya, A. Vegetation analysis and tree population structure of tropical wet evergreen forests in and around Namdapha National Park, northeast India. Biodivers. Conserv. 2005, 14, 2109–2135. [Google Scholar] [CrossRef]
  34. Yam, G.; Tripathi, O.P. Tree diversity and community characteristics in Talle Wildlife Sanctuary, Arunachal Pradesh, Eastern Himalaya, India. J. Asia-Pac. Biodivers. 2016, 9, 160–165. [Google Scholar] [CrossRef] [Green Version]
  35. Tripathi, O.; Pandey, H.; Tripathi, R. Distribution, community characteristics and tree population structure of subtropical pine forest of Meghalaya, northeast India. Int. J. Ecol. Environ. Sci. 2004, 29, 207–214. [Google Scholar]
  36. Tripathi, O.P.; Tripathi, R.S. Community composition, structure and management of subtropical vegetation of forests in Meghalaya State, northeast India. Int. J. Biodivers. Sci. Ecosyst. Serv. Manag. 2011, 6, 157–163. [Google Scholar] [CrossRef]
  37. Sinha, S.; Badola, H.K.; Chhetri, B.; Gaira, K.S.; Lepcha, J.; Dhyani, P.P. Effect of altitude and climate in shaping the forest compositions of Singalila National Park in Khangchendzonga Landscape, Eastern Himalaya, India. J. Asia-Pac. Biodivers. 2018, 11, 267–275. [Google Scholar] [CrossRef]
  38. Gogoi, A.; Sahoo, U.K. Impact of anthropogenic disturbance on species diversity and vegetation structure of a lowland tropical rainforest of eastern Himalaya, India. J. Mt. Sci. 2018, 15, 2453–2465. [Google Scholar] [CrossRef]
  39. Dutta, G.; Devi, A. Plant diversity, population structure, and regeneration status in disturbed tropical forests in Assam, northeast India. J. For. Res. 2013, 24, 715–720. [Google Scholar] [CrossRef]
  40. Sundriyal, R.C.; Sharma, E.; Rai, L.K.; Rao, S.C. Tree structure, regeneration and woody biomass removal in a subtropical forest of Mamlay watershed in the Sikkim Himalaya. Vegetatio 1994, 113, 53–63. [Google Scholar] [CrossRef]
  41. Singh, H.B.; Sundriyal, R.C. Composition, Economic Use, and Nutrient Contents of Alpine Vegetation in the Khangchendzonga Biosphere Reserve, Sikkim Himalaya, India. Arct. Antarct. Alp. Res. 2005, 37, 591–601. [Google Scholar] [CrossRef] [Green Version]
  42. Tambe, S.; Rawat, G.S. The Alpine Vegetation of the Khangchendzonga Landscape, Sikkim Himalaya. Mt. Res. Dev. 2010, 30, 266–274. [Google Scholar] [CrossRef]
  43. Pandey, A.; Rai, S.; Kumar, D. Changes in vegetation attributes along an elevation gradient towards timberline in Khangchendzonga National Park, Sikkim. Trop. Ecol. 2018, 59, 259–271. [Google Scholar]
  44. Chettri, N.; Sharma, E.; Deb, D.C.; Sundriyal, R.C. Impact of Firewood Extraction on Tree Structure, Regeneration and Woody Biomass Productivity in a Trekking Corridor of the Sikkim Himalaya. Mt. Res. Dev. 2002, 22, 150–158. [Google Scholar] [CrossRef] [Green Version]
  45. Acharya, B.K.; Chettri, B.; Vijayan, L. Distribution pattern of trees along an elevation gradient of Eastern Himalaya, India. Acta Oecolog. 2011, 37, 329–336. [Google Scholar] [CrossRef]
  46. Sharma, N.; Behera, M.D.; Das, A.P.; Panda, R.M. Plant richness pattern in an elevation gradient in the Eastern Himalaya. Biodivers. Conserv. 2019, 28, 2085–2104. [Google Scholar] [CrossRef]
  47. Tambe, S.; Arrawatia, M.L.; Sharma, N. Assessing the Priorities for Sustainable Forest Management in the Sikkim Himalaya, India: A Remote Sensing Based Approach. J. Indian Soc. Remote Sens. 2011, 39, 555–564. [Google Scholar] [CrossRef]
  48. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  49. Legendre, P.; Legendre, L.F. Numerical Ecology; Elsevier: Amsterdam, The Netherlands, 2012; Volume 24. [Google Scholar]
  50. Gotelli, N.J.; Colwell, R.K. Estimating species richness. Biol. Divers. Front. Meas. Assess. 2011, 12, 39–54. [Google Scholar]
  51. Keel, S.; Gentry, A.H.; Spinzi, L. Using vegetation analysis to facilitate the selection of conservation sites in eastern Paraguay. Conserv. Biol. 1993, 7, 66–75. [Google Scholar] [CrossRef]
  52. Ganesh, T.; Ganesan, R.; Devy, M.S.; Davidar, P.; Bawa, K.S. Assessment of plant biodiversity at a mid-elevation evergreen forest of Kalakad-Mundanthurai Tiger Reserve, Western Ghats, India. Curr. Sci. 1996, 71, 379–392. [Google Scholar]
  53. Korner, C. The use of ‘altitude’ in ecological research. Trends Ecol. Evol. 2007, 22, 569–574. [Google Scholar] [CrossRef]
  54. Grierson, A.J.; Long, D.G. Flora of Bhutan, Including a Record of Plants from Sikkim; Royal Botanic Garden: Edinburgh, UK, 1983. [Google Scholar]
  55. Manish, K.; Telwala, Y.; Nautiyal, D.C.; Pandit, M.K. Modelling the impacts of future climate change on plant communities in the Himalaya: A case study from Eastern Himalaya, India. Model. Earth Syst. Environ. 2016, 2, 92. [Google Scholar] [CrossRef]
  56. Kharkwal, G.; Mehrotra, P.; Rawat, Y.S.; Pangtey, Y.P.S. Phytodiversity and growth form in relation to altitudinal gradient in the Central Himalayan (Kumaun) region of India. Curr. Sci. 2005, 89, 873–878. [Google Scholar]
  57. Mishra, B.; Tripathi, O.; Tripathi, R.; Pandey, H. Effects of anthropogenic disturbance on plant diversity and community structure of a sacred grove in Meghalaya, northeast India. Biodivers. Conserv. 2004, 13, 421–436. [Google Scholar] [CrossRef]
  58. Saxena, A.K.; Singh, J.S. Tree population-structure of certain Himalayan forest associations and implications concerning their future composition. Vegetatio 1984, 58, 61–69. [Google Scholar] [CrossRef]
  59. Veblen, T.T. Structure and dynamics of nothofagus forests near timberline in south-central Chile. Ecology 1979, 60, 937–945. [Google Scholar]
  60. Vetaas, O.R. The effect of environmental factors on the regeneration of Quercus semecarpifolia Sm. in Central Himalaya, Nepal. Plant Ecol. 2000, 146, 137–144. [Google Scholar] [CrossRef]
  61. Kanade, R.; John, R. Topographical influence on recent deforestation and degradation in the Sikkim Himalaya in India; Implications for conservation of East Himalayan broadleaf forest. Appl. Geogr. 2018, 92, 85–93. [Google Scholar] [CrossRef]
  62. Gudade, B.; Chhetri, P.; Gupta, U.; Deka, T.; Vijayan, A. Traditional practices of large cardamom cultivation in Sikkim and Darjeeling. Life Sci. Leafl. 2013, 9, 62–68. [Google Scholar]
Figure 1. Study area with sampling locations along an altitude gradient in Sikkim, India.
Figure 1. Study area with sampling locations along an altitude gradient in Sikkim, India.
Forests 10 00633 g001
Figure 2. Hierarchically clustered dendrogram and Nonmetic multidimensional scaling (NMDS) ordination on the Bray–Curtis distance estimated using abundance data. (a) Hierarchical clustering was done using Ward’s algorithm, the y-axis represents the Bray–Curtis distance shared among each cluster; (b) NMDS ordination superimposed with the result of the cluster analysis, more similar plots are grouped closer to another. The red, green and blue dots and lines correspond to the low-altitude forest (LF), mid-altitude forest (MF), and high-altitude forest (HF), forest categories in our study.
Figure 2. Hierarchically clustered dendrogram and Nonmetic multidimensional scaling (NMDS) ordination on the Bray–Curtis distance estimated using abundance data. (a) Hierarchical clustering was done using Ward’s algorithm, the y-axis represents the Bray–Curtis distance shared among each cluster; (b) NMDS ordination superimposed with the result of the cluster analysis, more similar plots are grouped closer to another. The red, green and blue dots and lines correspond to the low-altitude forest (LF), mid-altitude forest (MF), and high-altitude forest (HF), forest categories in our study.
Forests 10 00633 g002
Figure 3. Shepard diagram with the goodness of fit (R2) from the NMDS result presented in Figure 2.
Figure 3. Shepard diagram with the goodness of fit (R2) from the NMDS result presented in Figure 2.
Forests 10 00633 g003
Figure 4. The low-altitude forest (LF): (a) Engelhardtia spicata, one of the dominant trees in flowering; (b)Schima wallichii, Symplocos, and Castanopsis mixed forest; (cd) even sized Castanopsis and Schima wallichii trees recorded at 1400 m a.s.l.
Figure 4. The low-altitude forest (LF): (a) Engelhardtia spicata, one of the dominant trees in flowering; (b)Schima wallichii, Symplocos, and Castanopsis mixed forest; (cd) even sized Castanopsis and Schima wallichii trees recorded at 1400 m a.s.l.
Forests 10 00633 g004
Figure 5. Mid-altitude forest (MF): (a) Castanopsis hystrix with 82cm DBH recorded at 2100 m a.s.l.; (b) forests dominated by Castanopsis tribuloides at 2000 m a.s.l.; (cd) forests dominated by Symplocos species at 2300 m a.s.l.
Figure 5. Mid-altitude forest (MF): (a) Castanopsis hystrix with 82cm DBH recorded at 2100 m a.s.l.; (b) forests dominated by Castanopsis tribuloides at 2000 m a.s.l.; (cd) forests dominated by Symplocos species at 2300 m a.s.l.
Forests 10 00633 g005
Figure 6. High altitude forest (HF): (a) Lithocarpus pachyphyllus with 110 cm DBH recorded at 2860 m a.s.l.; (b) Lithocarpus pachyphyllus and Rhododendron mixed forest 2875 m a.s.l.; (c) Open forest at 2792 m a.s.l.,where Lithocarpus pachyphyllus, Magnolia campbellii and Acer campbellii Hoof.f.& Thomson ex Hiern., were reported as the major tree species; (d) a Rhododendron forest stand at 3200 m a.s.l.
Figure 6. High altitude forest (HF): (a) Lithocarpus pachyphyllus with 110 cm DBH recorded at 2860 m a.s.l.; (b) Lithocarpus pachyphyllus and Rhododendron mixed forest 2875 m a.s.l.; (c) Open forest at 2792 m a.s.l.,where Lithocarpus pachyphyllus, Magnolia campbellii and Acer campbellii Hoof.f.& Thomson ex Hiern., were reported as the major tree species; (d) a Rhododendron forest stand at 3200 m a.s.l.
Forests 10 00633 g006
Figure 7. Species accumulation curve based on the rarefaction method: (a) whole tree population; (b)three forests: LF = low-altitude forest, MF = mid-altitude forest, HF= high-altitude forest.
Figure 7. Species accumulation curve based on the rarefaction method: (a) whole tree population; (b)three forests: LF = low-altitude forest, MF = mid-altitude forest, HF= high-altitude forest.
Forests 10 00633 g007
Figure 8. Boxplot with jittered data points showing the distribution of (a) basal area, (b) stem density, and (c) DBH. The horizontal line crossing the box in the center represents the median, the box represents the 25th and the 75th percentiles. The vertical line outside the box represents the minimum and maximum values. ALL = total forest community, LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
Figure 8. Boxplot with jittered data points showing the distribution of (a) basal area, (b) stem density, and (c) DBH. The horizontal line crossing the box in the center represents the median, the box represents the 25th and the 75th percentiles. The vertical line outside the box represents the minimum and maximum values. ALL = total forest community, LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
Forests 10 00633 g008
Figure 9. Class distribution of tress with >3 cmDBH recorded in the study. (a) Overall stems recorded in the study, (b) DBH ranging from >3 to 93 cm among three forest communities, and (c) DBH ranging from >93 cm to 303 cm. LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
Figure 9. Class distribution of tress with >3 cmDBH recorded in the study. (a) Overall stems recorded in the study, (b) DBH ranging from >3 to 93 cm among three forest communities, and (c) DBH ranging from >93 cm to 303 cm. LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
Forests 10 00633 g009
Table 1. Richness, species diversity, stem densities, and basal area for different forest communities along an altitude gradient. The table provides the median and its standard error since richness showed a right-skewed distribution. LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
Table 1. Richness, species diversity, stem densities, and basal area for different forest communities along an altitude gradient. The table provides the median and its standard error since richness showed a right-skewed distribution. LF = low-altitude forest, MF = mid-altitude forest, HF = high-altitude forest.
ParametersALLLFMFHF
Altitude (m)900–3200900–1700>1700–2500>2500–3200
Total sampled area (ha)8.32.43.72.2
Total species richness 114687751
Total genus richness 75525437
Total family richness 42333525
Median of species’ richness ± SE14 ± 4.3416 ± 4.4314 ± 3.0410 ± 4.50
Median of genus’ richness ± SE12 ± 3.9615.50 ± 4.0612 ± 2.7210 ± 3.51
Median of families’ richness ± SE10 ± 2.8913 ± 2.949 ± 1.959 ± 2.81
Total individual counted10344259154632290
Shannon index (H’) ± SD1.87 ± 0.442.10 ± 0.351.71 ± 0.401.78 ± 0.48
Fisher’s alpha richness4.44 ± 1.845.68 ± 2.004.14 ± 1.353.51 ± 1.64
Pielou’s evenness (J)0.39 ± 0.390.46 ± 0.70.35 ± 0.100.38 ± 0.13
Dominant familyFagaceaeFagaceaeFagaceaeFagaceae
Dominant speciesLithocarpus pachyphyllusSchima wallichiiSymplocos ramosissimaLithocarpus pachyphyllus
Table 2. Species richness, density, basal area, Species Importance Value Index (SIVI), Families Importance Value Index (FIVI), and their proportion to the contribution for the top ten species and families for the overall forest community recorded in the study.
Table 2. Species richness, density, basal area, Species Importance Value Index (SIVI), Families Importance Value Index (FIVI), and their proportion to the contribution for the top ten species and families for the overall forest community recorded in the study.
SpeciesFamilyStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%SIVI%
Lithocarpus pachyphyllusFagaceae51.204.1112.3121.079.30
Symplocos ramosissimaSymplocaceae248.0719.191.612.768.72
Castanopsis hystrixFagaceae62.054.989.9116.978.23
Castanopsis tribuloidesFagaceae62.415.014.227.234.99
Schima wallichiiTheaceae78.196.272.404.114.46
Rhododendron arboreumEricaceae85.546.862.544.354.33
Viburnum erubescensAdoxaceae103.988.340.380.644.16
Eurya acuminata DC.Pentaphylacaceae60.964.890.811.393.91
Quercus lamellosaFagaceae9.880.794.958.493.86
Symplocos heishanensisSymplocaceae49.403.960.380.652.93
FamiliesRichnessStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%FIVI%
Fagaceae920.8816.7533.6157.1927.86
Symplocaceae433.5526.922.183.7112.26
Theaceae213.9011.163.205.458.04
Ericaceae714.7611.843.816.487.70
Lauraceae134.633.714.978.456.62
Adoxaceae210.418.350.380.644.64
Magnoliaceae41.120.902.093.642.70
Betulaceae41.931.550.991.692.19
Primulaceae32.281.830.410.232.10
Sapindaceae10.690.551.232.101.89
Table 3. Stem density, basal area, Family (FIVI) and Species Importance Value Indexes (SIVI) and their proportion to the contribution of the ten dominant species for three forest communities: LF, MF, and HF.
Table 3. Stem density, basal area, Family (FIVI) and Species Importance Value Indexes (SIVI) and their proportion to the contribution of the ten dominant species for three forest communities: LF, MF, and HF.
SpeciesFamilyStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%SIVI%
LF (900 m to 1700 m a.s.l.)
Schima wallichiiTheaceae259.1724.017.0826.7318.90
Castanopsis tribuloidesFagaceae202.5018.767.6829.0017.91
Engelhardtia spicataJuglandaceae60.425.602.358.866.72
Gynocardia odorataAchariaceae92.508.571.234.636.05
Euryaa cuminataPentaphylacaceae66.678.570.762.864.83
Betula alnoides Buch.-Ham.ex.D.DonBetulaceae21.256.180.993.752.57
Macaranga denticulata BlumeEuphorbiaceae18.081.970.391.482.43
Albizia lebbeck (L.) Benth.Leguminosae17.081.581.234.642.31
Alnus nepalensis D.DonBetulaceae12.080.661.003.772.21
Brassaiopsis hispida Seem.Araliaceae36.253.360.190.531.87
FamiliesRichnessStem density (ha−1)Stem density%Basal area (m2 ha−1)Basal area%FIVI%
Theaceae29.4230.187.8429.6022.54
Fagaceae66.7321.578.2331.0920.17
Juglandaceae31.805.752.388.997.42
Achariaceae12.678.571.234.636.58
Betulaceae31.043.322.097.895.26
Euphorbiaceae51.233.940.592.224.34
Araliaceae31.575.020.170.663.53
Leguminosae30.300.961.274.802.90
Anacardiaceae30.591.890.481.812.76
Lauraceae60.632.010.090.352.64
SpeciesFamilyStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%SIVI%
MF (>1700 m to 2500 m a.s.l.)
Symplocos ramosissimaSymplocaceae501.8933.993.364.6715.09
Castanopsis hystrixFagaceae132.438.9721.7230.1514.99
Viburnum erubescensAdoxaceae211.8914.350.751.057.09
Quercus lamellosaFagaceae20.541.3910.6714.816.85
Eurya acuminataPentaphylacaceae83.785.671.231.714.42
Lithocarpus pachyphyllusFagaceae48.923.315.177.184.32
Symplocos heishanensisSymplocaceae73.785.000.560.784.07
Castanopsis tribuloidesFagaceae8.650.594.496.232.78
Quercus glaucaFagaceae14.320.973.094.292.70
Elaeocarpus sikkimensis Mast.Elaeocarpaceae10.540.712.213.072.52
FamiliesRichnessStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%FIVI%
Fagaceae710.8216.4445.7045.7030.12
Symplocaceae428.9043.914.334.3320.14
Lauraceae113.084.699.609.609.22
Adoxaceae29.4614.370.750.758.07
Theaceae23.956.001.711.715.73
Elaeocarpaceae20.480.732.232.233.17
Ericaceae52.493.790.490.492.72
Sapindaceae10.400.601.501.502.60
Primulaceae31.241.890.160.162.60
Magnoliaceae40.590.901.751.752.43
SpeciesFamilyStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%SIVI%
HF (2500 m to 3200 m a.s.l.)
Lithocarpus pachyphyllusFagaceae110.9110.667.8854.5624.37
Rhododendron arboreumEricaceae310.0029.787.4713.0816.78
Magnolia campbelliiMagnoliaceae16.821.625.396.624.54
Symplocos heishanensisSymplocaceae61.825.945.810.694.14
Rhododendron hodgsonii Hook.f.Ericaceae54.555.242.903.153.76
Corylus ferox Wall.Betulaceae26.822.584.151.332.68
Viburnum erubescens Wall.Adoxaceae35.913.454.150.212.60
Acer campbelliiSapindaceae10.911.053.732.982.59
Rhododendron falconeri Hook.f.Ericaceae39.553.802.071.862.58
FamiliesRichnessStem Density (ha−1)Stem Density%Basal Area (m2 ha−1)Basal Area%FIVI%
Fagaceae43.3390.1140.9657.1679.21
Ericaceae711.8229.6313.4718.7972.16
Symplocaceae34.642.070.941.3126.02
Magnoliaceae20.4810.904.956.9116.03
Lauraceae60.925.512.513.5015.24
Betulaceae20.762.160.981.379.91
Theaceae20.531.460.670.939.69
Berberidaceae21.390.370.170.249.47
Adoxaceae10.950.330.150.219.92
Sapindaceae10.294.702.142.988.77

Share and Cite

MDPI and ACS Style

Bhutia, Y.; Gudasalamani, R.; Ganesan, R.; Saha, S. Assessing Forest Structure and Composition along the Altitudinal Gradient in the State of Sikkim, Eastern Himalayas, India. Forests 2019, 10, 633. https://doi.org/10.3390/f10080633

AMA Style

Bhutia Y, Gudasalamani R, Ganesan R, Saha S. Assessing Forest Structure and Composition along the Altitudinal Gradient in the State of Sikkim, Eastern Himalayas, India. Forests. 2019; 10(8):633. https://doi.org/10.3390/f10080633

Chicago/Turabian Style

Bhutia, Yangchenla, Ravikanth Gudasalamani, Rengaian Ganesan, and Somidh Saha. 2019. "Assessing Forest Structure and Composition along the Altitudinal Gradient in the State of Sikkim, Eastern Himalayas, India" Forests 10, no. 8: 633. https://doi.org/10.3390/f10080633

APA Style

Bhutia, Y., Gudasalamani, R., Ganesan, R., & Saha, S. (2019). Assessing Forest Structure and Composition along the Altitudinal Gradient in the State of Sikkim, Eastern Himalayas, India. Forests, 10(8), 633. https://doi.org/10.3390/f10080633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop