Stand structure and species diversity regulate biomass carbon stock under major Central Himalayan forest types of India

Data on the impact of species diversity on biomass in the Central Himalayas, along with stand structural attributes is sparse and inconsistent. Moreover, few studies in the region have related population structure and the influence of large trees on biomass. Such data is crucial for maintaining Himalayan biodiversity and carbon stock. Therefore, we investigated these relationships in major Central Himalayan forest types using non-destructive methodologies to determine key factors and underlying mechanisms. Tropical Shorea robusta dominant forest has the highest total biomass density (1280.79 Mg ha−1) and total carbon density (577.77 Mg C ha−1) along with the highest total species richness (21 species). The stem density ranged between 153 and 457 trees ha−1 with large trees (> 70 cm diameter) contributing 0–22%. Conifer dominant forest types had higher median diameter and Cedrus deodara forest had the highest growing stock (718.87 m3 ha−1); furthermore, C. deodara contributed maximally toward total carbon density (14.6%) among all the 53 species combined. Quercus semecarpifolia–Rhododendron arboreum association forest had the highest total basal area (94.75 m2 ha−1). We found large trees to contribute up to 65% of the growing stock. Nine percent of the species contributed more than 50% of the carbon stock. Species dominance regulated the growing stock significantly (R2 = 0.707, p < 0.001). Temperate forest types had heterogeneous biomass distribution within the forest stands. We found total basal area, large tree density, maximum diameter, species richness, and species diversity as the predominant variables with a significant positive influence on biomass carbon stock. Both structural attributes and diversity influenced the ordination of study sites under PCA analysis. Elevation showed no significant correlation with either biomass or species diversity components. The results suggest biomass hyperdominance with both selection effects and niche complementarity to play a complex mechanism in enhancing Central Himalayan biomass carbon stock. Major climax forests are in an alarming state regarding future carbon security. Large trees and selective species act as key regulators of biomass stocks; however, species diversity also has a positive influence and should also reflect under management implications.


Background
For adequate ecosystem functioning and management, biomass as well as species diversity together play a critical role and they also address the two most predominant challenges viz. climate change and biodiversity loss. Experimental models and long-term studies under grassland ecosystems have established a positive relationship between species diversity and productivity (Tilman et al. 1997;Tilman et al. 2001). These studies regarded niche complementarity as the underlying principle behind the effects of diversity on productivity with a rationale that resource use and availability is greater under diverse species. However, the effects become more prominent in the long term. In short-term studies, the single species dominant systems appear the most productive (Cardinale et al. 2007). Several studies (Szwagrzyk and Gazda 2007;Day et al. 2014;Mensah et al. 2016;Zhang et al. 2017;Li et al. 2018) under forest ecosystem implicate a positive relationship between tree species diversity and biomass; however, most of them reported weak positive correlation. Furthermore, the tree strata diversity affects the diversity and biomass relation of the understory vegetation (shrubs and herbs). The trend, however, is not ubiquitous. Jacob et al. (2010) reported higher biomass under monocultures than under diverse oldgrowth forest counterparts. While, mature to old-growth primary forests are indispensable for their carbon and biodiversity repository, and it is primarily the presence of large diameter trees that hold major carbon stock in these forests (Gibson et al. 2011). From a global perspective, it is evident that irrespective of their density, the large diameter trees (diameter > 60 cm) contribute significantly to biomass and their loss could yield a reduction in structural heterogeneity and the carbon capture potential of forests (Lutz et al. 2018). These mature forests, therefore are more resilient and provide long-term carbon pool reservoirs (Day et al. 2014). Thus, apart from diversity, the forest structure is also critical in deciphering management implications. Estimation of biomass carbon has other imperatives as well. Consumption of fossil fuels leads to greenhouse gas emissions which have increased substantially (Chu 2009). Around 80% of the total energy consumption in 2018-2019 was from coal, lignite, and crude oil as a major fuel source in India (Energy Statistics 2020). Consequently, the global atmospheric annual CO 2 concentration has increased from 338.8 ppm in 1980 to 409.8 ppm in 2019, gaining 71 ppm in only 39 years (Dlugokencky and Tans 2020). Loading of this labile CO 2 into the atmospheric carbon (C) pool needs to be managed. One route to achieve this is through mitigation options which involve reducing the emissions and sequestering the labile emissions present into long-term pools (Lal 2008). A practical climate stabilization measure among other sequestration alternatives is terrestrial biotic carbon sequestration. Forest ecosystems are the largest of all terrestrial ecosystems, globally covering 4.06 billion hectares (30.8%) of area and India contributes 2% to this global forest area (FAO and UNEP 2020). Forests have a huge repository of C, storing 360 Pg C (1 Pg = 10 15 g) in live and dead biomass components along with 398 Pg C stored in forest soils globally. Along with this, the forest ecosystems sequester around 1.7 ± 0.5 Pg C year −1 . However, destructive anthropogenic activities like deforestation make forests a net source, e.g., deforestation of tropical forests alone release 0.3 Pg C year −1 ; therefore, the forest ecosystems are not only viable sequestration option but they need efficient management to preserve their C stock (Tans et al. 1990;Lal 2008;Lorenz and Lal 2010). India's forest cover amount to 71.2 million hectares which constitute 21.67% of the total geographic area (GA) of the country (FSI 2019). Situated in the foothills of the Himalayan mountain ranges lie the Indian state of Uttarakhand; it constitutes the central portion of the Indian Himalayan Region (IHR) (Kafaltia and Kafaltia 2019). Uttarakhand has the second-highest forest growing stock (406.8 million cubic meters) and the fourth-highest forest carbon stock per hectare in aboveground biomass (62.77 Mg ha −1 ) in India (FSI 2019). In the Central Himalayan grasslands, the species diversity and productivity relationships are well established (Bhattarai et al. 2004;Singh et al. 2005). The information on tree species, however, is inconsistent. Previous studies in the region determined associations between species diversity and biomass; however, elevation attributed the predominant focus, and the impact of diversity or stand structural attributes alone were not elucidated (Singh et al. 1994;Sharma et al. 2010;. Moreover, the reports revealed contrasting associations between diversity and biomass in Central Himalayas.  reported no significant correlation between total tree carbon density (TCD) and diversity in twelve major temperate Himalayan forest types, while Sharma et al. (2010) showed a significant negative correlation between TCD and species diversity in their study involving twenty major Garhwal Himalayan forest types. With limited studies under the topic, this creates a need for understanding this relation in Central Himalayas, since apart from carbon management, the conservation of biodiversity is also crucial for optimum ecosystem functioning. Therefore, the present study was undertaken with two broad objectives: (i) determination of growing stock, biomass, and carbon stock of major forest types in Central Himalaya along with their comparison with similar forest types to access their status, and (ii) understanding the effects of stand structural attributes and species diversity on biomass and carbon stocks.

Study area description
The Garhwal Himalayas of Uttarakhand state in India ( Fig. 1) which lies in the central portion of the Indian Himalayan Region covering an area of around 32,449 km 2 represented our study area. Situated in the foothills of the Himalayas, it gives rise to the two prominent river systems the Ganges and the Yamuna-Tons basin (Nandy et al. 2006). The study sites are situated in two districts viz. Dehradun and Rudraprayag lying between 30°00′ N-30°46′ N latitude and 77°48′ E-79°12′ E longitude covering an elevational range of around 3000 m. To select areas that represent major forest types of Uttarakhand, we selected the forest types having the maximum forest cover in the state as per FSI (2019) report. The forest types thus selected represents 63.4% of Uttarakhand's forest cover (FSI 2019). We conducted a preliminary reconnaissance study in 2017 to select the study areas that comprise of these forest types. The selected areas represent large stands of a particular forest type, constitute naturally regenerating forests without plantations, and have a protected legal status to minimize anthropogenic influence. Based on these criteria along with the ease of access and availability of field research permission, we selected eight forest types as per Champion and Seth (1968) forest type classification (detailed in Table 1 and Supplementary Table S1). Physiographically, the region rises from the Bhabar-Siwalik zone of the outer Himalaya to the Himadri of the inner Himalayas (Rajwar 1993). F1 (RNP) and F2 (ASA) had generally flat topography with a < 10% slope. All the other sites (F3-F8) had undulating hilly terrain with highly variable slope angles in a range of 20°to 35°and with low (< 30%) to medium (30-80%) rockiness. The region has a seasonal climate with three major seasons, i.e., winter season (November to February), summer season (April to June), and rainy (July to September). A small period of spring (March) and autumn (October) are also present between the seasons (Singh and Singh 1987).
The overall climate follows the monsoon pattern of rainfall having most of the annual rainfall from mid-June to mid-September. Therefore, the summers represent the warm dry months and the rainy season represents the warm and wet months. The maximum temperature in 2019 ranged between 24 and 46°C and minimum between 3 and − 30°C (AccuWeather©, Inc. 2020). Snowfall in higher elevations begins as early as October with heavy snowfall during January and February in temperate and sub-alpine elevations. In the sub-alpine elevation, F8 (TUN) the snowfall can continue up to April (Rai et al. 2020).
Furthermore, the microclimate within an area also varies due to slope aspects, topographic relief, windward, and leeward sides, and even position on the ridge (Kafaltia and Kafaltia 2019). Geologically, the study sites lie in three major tectonic belts viz. (a) the Sub-Himalayan (SH) belt in the lower elevations comprising of clays, sandstones, and conglomerates; (b) Lesser Himalayan (LH) belt which predominantly comprises of sedimentary rocks, limestones, grits, shales, and schists; and (c) The Higher Himalayan Crystalline (HHC) belt mainly comprising of gneisses and granite (Rajwar 1993;Mukherjee et al. 2019).

Study design and vegetation sampling
Based on FSI field inventory manual, we classified the general properties of the area like legal status, topography, and rockiness (FSI 2002). The geographic coordinates and elevation were estimated using a handheld GPS logger (Garmin ® GPS 72™) and the slope angle for the sites was determined using the topographic relief map from Google Maps (https://www.google.com/maps). For forest inventory, we randomly laid square quadrats of dimensions 31.6 m × 31.6 m (~0.1 ha) in each forest type. The minimum distance between the quadrats was not less than 300 m while in some forest types it was as large as 5000 m. The plot size is based on FSI field inventory manual (FSI 2002), and is adopted by other investigators in Himalayan forest types (Thokchom and Yadava 2017;Banik et al. 2018;Sharma et al. 2018). In each forest type F1-F6, we laid ten quadrats (10 × 0.1 = 1 ha); however, for upper temperate and sub-alpine elevations F7 and F8 due to harsh terrain and climatic conditions, we placed only five quadrats each (5 × 0.1 = 0.5 ha). Overall, 70 quadrats (7 ha) were laid to cover the forest types in this study. Within each quadrat, we identified all the tree species, i.e., species that classify as trees in their growth form by recognized floras, and have a diameter ≥ 10 cm at breast height (137 cm from the ground). The circumference over bark was measured at breast height (CBH) using a measuring tape as per standard rules (Singh et al. 1994;UNFCCC 2015). For taxonomic identification, we consulted local floras viz. Flora of the district Garhwal (Gaur 1999), Flora of Jaunsar Bawar (Agarwal 2017), Flora of Rajaji National Park (Singh and Prakash 2002), Plants of Kedarnath Wildlife Sanctuary, Western Himalaya: a field guide (Rai et al. 2017) as well as the taxonomic knowledge from the local community.

Stand structure and species diversity
The stand structure or the population structure was analyzed using density-diameter curves which were plotted along with the growing stock volume for a comparative Forest type groups and classes are based upon Champion and Seth (1968) forest type classification S visualization. The distribution and skewness of the DBH distribution are represented using a box-whisker plot indicating the overall girth distribution along with the distribution pattern of large trees (DBH > 70 cm). Various studies use different criteria to delineate large trees ranging from > 60 cm DBH to > 100 cm, and even associated biomass values (Lutz et al. 2013;Sist et al. 2014). In the absence of any unified classification system for the Himalayan region, we regard trees with DBH > 70 cm as large trees since this parameter is used by various investigators (Brown and Lugo 1992;Slik et al. 2013;Bradford and Murphy 2019) thus facilitating easier comparisons.
To determine the species dominance, we calculated the Importance Value Index (IVI) for each species. Summation of relative density, relative frequency, and relative basal area of the species constituted the IVI (Misra 1968

Growing stock, biomass, and carbon stock estimation
Since the climatic, edaphic, and topographic factors along with inherent anthropogenic influence modulate the biomass of a region, therefore, the use of general biomass allometric equations might be inaccurate for region-specific studies (Brown and Iverson 1992;Zianis 2008). Moreover, since destructive biomass sampling is prohibited in protected forest areas, and species-specific or general allometric equations for biomass estimation are not available for most of the species in the Indian Himalayan region, therefore we used a volume-based method. Other investigators (Brown et al. 1999;Chhabra et al. 2002;Dimri et al. 2017a) also used a similar approach. For each individual, wood volume (m 3 ) was therefore determined using species and region-specific volume equations (detailed in Table 2) as per the Forest Survey of India publications (FSI 1996;FSI 2011). For species where equations were not available, we used either genus-specific or general equations for the region ( Table 2). The growing stock volume density (GSVD, m 3 ha −1 ) for each forest type and their constituent species was determined through the summation of volume for each tree individual. Since for forest types, F7 and F8 only 0.5 ha area were sampled, the summation of growing stock volume (species wise) was extrapolated to 1 ha. To convert GSVD into aboveground biomass density (AGBD, Mg ha −1 ), we multiplied the GSVD with appropriate biomass expansion factor (BEF, Mg m −3 ) which converts wood volume into biomass and accounts for non-commercial tree components like branches, twigs, bark, foliage, and reproductive parts thus, providing total tree AGBD (Schroeder et al. 1997). Table 3 details the regression equations of BEFs and BEF values (Brown et al. 1999). To determine the belowground biomass density (BGBD, Mg ha −1 ) which represents the biomass of root components, we used the regression equation (Table 3) given by Cairns et al. (1997) (Brown et al. 1999). The equation determines the BGBD based upon AGBD and latitudinal zone (temperate zone in our case). Summation of AGBD and BGBD yields the total biomass density (TBD, Mg ha −1 ) for trees. For aboveground carbon density (AGCD, Mg C ha −1 ) and total carbon density (TCD, Mg C ha −1 ), we multiplied the AGBD and TBD, respectively with species-specific carbon factor.

Statistical analysis
To understand the impact of species dominance over growing stock distribution, we ran a linear regression between Log e IVI and Log e GSVD (at species level) to quantify the dependence between the variables. For determining the evenness in growing stock distribution within the forest types, we used one-way analysis of variance (ANOVA) to compare means of GSVD at quadrat level within forest types. A two-tailed Pearson's correlation between the topographic (elevation) variable, total stem density, large tree stem density, maximum tree diameter, TBA, GSVD, TBD, TCD, TSR, MR, and H′ using IBM SPSS software (Ver. 23.0.0.0) determined the degree and direction of the association between the variables. Mean values for the entire forest types were used to run the correlation to avoid error bias due to plot variances because of the smaller plot size (0.1 ha). The dimension reduction was done by principal component analysis (PCA) using Paleontological Statistics software (PAST, Version 4.03) of biomass variable (TBD), species richness variables (TSR and MR), diversity variable (H′), and stand structure variable (large tree stem density). PCA was run using a correlation matrix, and components with cumulative eigenvalue (as %) 80% or greater were retained. Furthermore, a scree plot of eigenvalues (as % variation) is also considered, and components at the elbow point (below the broken stick model) were considered less significant (Zuur et al. 2007). Loading coefficients and PCA biplot were used to interpret the association between the variables and principal axes. Abies spectabilis (D. Don) Mirb.

Association of biomass carbon stock with stand structure and large trees
The biomass distribution between pure and mixed forest stands did not show any trend upon comparison with similar forest types in this study. For instance, Sal mixed forest (F1) had higher biomass density than pure Sal counterpart (F2); however, pure C. deodara forest (F6) had higher biomass density than mixed oak-deodar forest (F4) ( Table 4). Under this aspect, population structure and species composition played a significant role. a b  The density-diameter curves (Fig. 2) not only represent the population structure of a forest but also give interpretation over age structure, successional stage, and regeneration status (Rao et al. 1990). Forest types F1, F2, F4, and F8 show lower stem density in the lowest diameter classes (Figs. 2 and 3). This indicates poor regeneration potential in these forests or selective harvest of smaller diameter individuals. Being forest types with climax species viz. S. robusta, Q. oblongata, and A. spectabilis, this loss of regeneration could influence the future carbon stocks of these forests drastically. Previous studies (Singh and Singh 1986;Singh and Singh 1987;Rai et al. 2012) in the region also highlight the regeneration loss. F2 and F4 indicate unimodal stacking of individuals (Fig. 2) in the intermediate DBH classes along with a lower concentration of large trees (Fig. 3). This is further corroborated with the least spread (smaller box) in the box plot analysis (Fig. 3), suggesting that these forests are still maturing and represent even-aged forests having maximum GSVD stacking in intermediate diameter classes. These forests require stringent management since natural or anthropogenic disturbance could diminish the carbon stock of these forests especially F2 wherein S. robusta is highly susceptible to pest infestation (Singh and Thapa 1988). Species-specific pest infestations impact the carbon stock of a region immensely, causing the loss of standing biomass. Furthermore, the snags and deadwood thus produced affect the regeneration potential, species composition, carbon pool, and fire dynamics in the region (Lutz et al. 2021). F3, F5, F6, and F7 showed efficient regeneration with maximum GSVD allocation to the larger diameter classes (Figs. 2 and 3). The most prominent observation from stand structural analysis yields large trees as the chief contributors toward growing stock (Fig. 2). For the even-aged forest stands, the diameter classes with the maximum stem density did not have maximum GSVD (Fig. 3). Our study, therefore, signifies that overall stem density had no relation with either growing stock or biomass carbon stock. Furthermore, a significant positive association of total basal area, stem density of large trees with GSVD, and maximum DBH with TBD and TCD signify the importance of large trees in regulating the growing stock. Thus, the occurrence of trees in the larger diameter classes instead of total stem density is a major stand structural factor regulating the growing stock in these major Central Himalayan forest types. Sharma et al. (2010) also observed no significant correlation between stem density and total carbon density in twenty major forest types of Garhwal Himalayas. In this study, the large trees though constituting only 0-22.3% of stem density contributed 0-65.1% of the total GSVD. In Eastern Himalayas, also large trees contributed significantly to biomass carbon stock (Sundriyal and Sharma 1996;Baishya and Barik 2011). The enormous carbon stock of large trees has global recognition. Brown and Lugo (1992) concluded that large trees contributed around 40% to aboveground biomass in tropical Amazon forest with stem density no more than 3% of the total. Slik et al. (2013) also reported 25.1-44.5% biomass allocation to large trees in the Neotropical and Paleotropical forest where such trees occupied no more than 3.8% of the total stem density. Large trees contributed 33% of the total biomass in the Australian wet tropical rainforest though being only 2.4% of total stem density (Bradford and Murphy 2019). This indicates that these large trees are rare in the forest structure, the box-plot also place such individuals as outliers from the general population distribution (Fig. 3). Therefore, the loss of even a single large tree though may not give any pronounced effect over stem density, and may go unnoticed but would reduce the carbon stock of the forest substantially. Furthermore, large trees are regarded as important keystone species in various terrestrial ecosystems, and prevent the establishment of invasive and ruderal species (Gandhi and Sundarapandian 2017;Bradford and Murphy 2019). Luyssaert et al. (2008) and Körner (2017) aptly justified the significance of oldgrowth large trees in maintaining the carbon pool globally as well as acting as efficient global carbon sinks refuting the claim that large trees are not efficient carbon sequesters. There is significance of small and intermediate diameter trees since these individuals would enhance the carbon stock of the forest types with time on account of their higher carbon sequestration potential (Gandhi and Sundarapandian 2017). In our study, we saw the preponderance of intermediate DBH classes in all forest types except F5 and F7 (Fig. 2); this explains lower recruitment and exploitation of higher diameter classes in maturing forests. A similar trend was also observed by Saxena et al. (1984) in Kumaun Himalayas. F5 and F7 indicate decreasing tree density with increasing DBH classes (Fig. 2), and Sundriyal and Sharma (1996) regarded such a trend as a sign of healthy forest structure.
ANOVA results show between groups (quadrats) variability among all sub-tropical (F3) and temperate (F4-F7) forest types, also visualized by broader IQ range for these zones (Fig. 3). Uneven distribution of large trees predominantly causes such heterogeneity. Similar observations were also reported by Gandhi and Sundarapandian (2017) in the tropical dry deciduous forest of Eastern Ghats, India. Thus, even the slightest destruction on account of anthropogenic activity could lead to greater carbon stock loss for such forests.

Association of biomass carbon stock with species diversity
Species composition plays a vital role in modulating forest carbon stocks. Negi et al. (2003) reported maximum carbon allocation in the conifer woods followed by dicot deciduous and dicot evergreen species. This pattern is evident since the top ten species in order of TCD (Supplementary Table S3) are either conifers or dicot deciduous. Oaks being evergreen still have a higher carbon content in their wood along with a longer rotation period as compared to other evergreen species leading to higher TCD (Sharma et al. 2010). The total species richness recorded in this study (4-21 species) is within the range reported (7-22 species) for Garhwal Himalayas (Sharma et al. 2018). Our diversity (H′) value (0.77-2.63) was higher than the values reported (0.28-1.75) for major forest types of Garhwal Himalayas (Sharma et al. 2010). Since the study sites fall under reserved forests, wildlife sanctuary, and national park, therefore, the anthropogenic disturbance is quite limited thereby realizing the true diversity potential of the region. The regression analysis (Fig. 4) indicates the importance of species dominance as a predictor in growing stock contribution. IVI incorporates the three major components of a species viz. relative density, relative frequency, and relative basal area, therefore IVI gives a true sense of species dominance to a species relative to its population. This relation justifies the organization of species in decreasing order of TCD contribution (Supplementary Table S3) since the top-level species are either dominant or co-dominant species for different forest types. This relation is significant under forest management since basic phytosociology could identify prominent species that could be earmarked for conservation and afforestation plans.
In our study, just 9.4% of the total species accounted for 50.6% of the TCD in the study (Supplementary Table  S3). Corollary to this in the Amazonian forest just 1% of the total species account for 50% of the total carbon stock and productivity (Fauset et al. 2015). This indicates that there is biomass hyperdominance in these major Central Himalayan forest types. Biomass hyperdominance is a concept wherein a small number of species have a disproportionate effect on total biomass and productivity (Fauset et al. 2015). Such species thus become indispensable under management implications. Our results also indicate a strong significant positive correlation between species richness, diversity, and biomass density. Gogoi et al. (2020) in Eastern Himalayas reported a weak positive association between species diversity (H′) and TCD (R 2 = 0.25, p < 0.05). Behera et al. (2017) reported strong positive correlation between species richness, species diversity, and above-ground biomass in Indian tropical deciduous forests with S. robusta and Tectona grandis as major elements. Several investigators have also shown a positive relationship between species richness and biomass globally (Day et al. 2014;Mensah et al. 2016;Zhang et al. 2017;Sintayehu et al. 2020). This association of biomass with species diversity is mutually beneficial for conservation, and plays a key role under management implication. Any management plan must not only focus on biomass enrichment but also aim to maintain if not enrich the biodiversity. In this regard, monospecific plantations would harm the ecosystem functions of such forests. Furthermore, with appropriate management strategies, the forest types studied provide co-benefit of biomass and biodiversity enhancement. It is recommended to consider the prominent native species for each forest type along with endemic and threatened Himalayan species for management practices. In recent times, exotic and commercial species such as T. grandis, Eucalyptus sp., P. roxburghii are being given more preference over the native species for management purposes (Singh and Singh 1987). These species have a smaller rotation period than their counterparts like S. robusta and oaks (Sharma et al. 2010), thus such management efforts would yield only short-term carbon gains on the cost of long-term carbon storage and biodiversity.
The results from PCA analysis ( Fig. 5 and Table 5) signify the role of both large tree density and species diversity in influencing the ordination of study sites. F1 and F5 cluster together in the positive direction of PC 1 since they have maximum diversity, species richness, and TBD (Fig. 5b). Therefore, in these sites, management of diversity is crucial for securing carbon stock. Sites F6, F7, and F3 ordinate in the positive plane of PC2 having a maximum correlation with large tree stem density (Fig. 5b). These sites thus require stringent management to secure large trees to preserve the carbon stock; moreover, these sites represent poor species richness and diversity. Management strategies could focus on this aspect too. Sites F2, F8, and F4 club in the bottom left quadrant indicating that these sites have the lowest values for all the variables as all the variables point in the positive and increasing values of components (Fig. 5b). These sites require special attention, the management strategy for such regions require focus for the preservation of natural regeneration, enrichment of diversity through plantations, and preservation of existing large diameter trees. Sintayehu et al. (2020) reported that functional diversity traits (maximum DBH, maximum tree height, wood density) are more strongly correlated with Ethiopian woodland aboveground biomass than species diversity components. Based on this, they considered that the natural ecosystems might have achieved their species richness potential. In our study, both the functional traits (maximum DBH) and species richness traits correlate significantly with biomass, indicating that our forests can still accommodate more species. Therefore, management implications can aim in enhancing the biodiversity of forest stands which would concurrently enhance the carbon stocks. Our PCA findings were comparable to Day et al. (2014) in the Central African rainforest; however, their second PCA component had a significant correlation between stem density and biomass. This indicates that different sites have different factors that yield maximum biomass production.
The two most predominant mechanisms elucidated for the effect of species richness over biomass are complementarity effects (diverse species have different niche requirements leading to resource partitioning thus leading to maximum resource utilization and enhanced productivity) and selection effects (enhanced productivity on account of the dominance of species with specific traits) (Loreau and Hector 2001;Cardinale et al. 2007). Under our observations, not only a few dominant species and species with larger diameter contribute to biomass carbon stocks but the species diversity and richness are also positively associated. Thus, these forest types have a complex interplay of both mechanisms. Similar results were reported by Cavanaugh et al. (2014) on a global scale. We, however, do not generalize our findings to Central Himalayas since we did not consider all the forest types present in the region, and on account of the smaller sample size. Zhang et al. (2017) also emphasized that positive species diversity and biomass relationships though ubiquitous across all vegetation strata; however, the overstorey diversity and biomass profoundly impact the understorey biomass predominantly through resource filtering. This study would help researchers focusing on understorey diversity and biomass relations in the region. Furthermore, the relationship between species diversity and biomass is inconsistent; there were reports of negative to no significant association between species diversity and biomass within the region (Sharma et al. 2010;), in India (Sahu et al. 2016Behera et al. 2017;Gandhi and Sundarapandian 2017) and the world (Jacob et al. 2010;Szwagrzyk and Gazda 2007). Therefore, on a regional scale, such studies are imperative for deciphering adequate management strategies since the associations between variables might get diluted in large-scale studies thereafter hindering forest management and carbon stock enhancement at the regional scale.

Variations in biomass carbon stocks
The segregation of TBD into above and belowground components showed a distinct pattern in all the forest types. The ccontribution of AGBD (79.6%) and BGBD (20.4%) to the TBD in the present study is comparable to the value reported (78.54% AGBD and 21.46% BGBD) by Sharma et al. (2018) for Garhwal Himalayas, (79% AGBD and 21% BGBD) for pan-India biomass (Chhabra et al. 2002), and Brazilian Amazon tropical moist forest (10-50% by BGBD) (Brown and Lugo 1992). Lower contribution by belowground biomass is normal since plants invest in belowground biomass significantly only under stressful conditions. In our study, elevation had a negative association with stem density, species richness, species diversity, TBD, and TCD; however, the associations were not significant. Singh et al. (1994) reported that elevation had a significant negative association with stem density, species richness, species diversity, and biomass in the Central Himalayan forest types. Since we selected predominant Central Himalayan forest types and not all within the elevational range therefore, our trends with elevation were not significant. At elevations around 2800-2900 m, we saw a spike in stem density, basal area, and biomass due to much-branched Krumholtz like growth of R. arboreum in these elevations which forms a major under canopy tree layer contributing significantly to stem density, basal area, and biomass (Singh et al. 1994;Rai et al. 2018). The treeline ecotone forest F8 depicts lower values across all parameters since these forests form a limit of closed-canopy forests, and only a few species acclimated to the harsh alpine conditions thrive at such elevations. Furthermore, new recruitment in A. spectabilis dominant Tungnath (F8) treeline ecotone is sparse which is reflected in the population structure also (Fig.  2) (Rai et al. 2012). Stringent management approach is deemed necessary in such regions to facilitate and conserve natural regeneration to preserve treeline carbon stock. Sub-tropical pine forest (F2) recorded the lowest carbon stock in our elevational range since these forests are inflicted with forest fires which limit their productivity and affect regeneration (Sharma et al. 2010). On average, temperate forest types had higher carbon stock due to the presence of mature to old-growth large-diameter hardwood and conifer trees (Dimri et al. 2017a).
The estimates for growing stock, biomass, and carbon stock though give a quantitative measure to compare the forest types; however, the state of the forests could only be visualized upon comparing the results with similar forest types. Our biomass and carbon stocks were higher or within the range as reported by previous investigators (Table 6). Our estimates were most comparable with either mature old-growth forests, community-managed forests, or sacred grooves (Table 6). These forests are pristine undisturbed forests representing the true potential of biomass and carbon stock in the region. Such forests are used as control sites in biomass modeling studies (Osuri et al. 2014). This study, therefore, elucidates that the selected predominant Central Himalayan forest types have huge biomass carbon stock and thus, require protection from anthropogenic activities.

Conclusions
The knowledge under growing stock provides quantitative value to the forest land helping forest managers to delineate regions as carbon-rich or deficient, and proceed accordingly. Results from quantitative stand structural attributes implicate predominant Himalayan climax species and forest types have a preponderance of intermediate DBH class along with a significant decline in new recruitment. The total basal area, the density of large trees, and maximum DBH are key stand structural factors that regulate the biomass carbon stock. Large trees although low in abundance (up to 22%) were key contributors to the total biomass (up to 65%). Though species richness and diversity correlated significantly with biomass carbon density only a handful of species (9.4%) contributed more than 50% of the total biomass indicating biomass hyperdominance. The biomass distribution was also heterogeneous within the forest types, and species dominance regulate the growing stock directly. Elevation did not show significant correlations indicating that the entire elevational range has huge carbon stock. Regulation of biomass carbon by both structural and diversity components suggest a complex interplay of complementarity effects and selection effect that enhance the carbon stock of the region. The comparative review illustrates that these Central Himalayan forest types are carbon-rich and represent stock equivalent to sacred grooves and communitymanaged forests where communities preserve the forest resources.
Under management implications, we recommend that for any timber requirements the small diameter trees are a viable alternative since they have greater stem density and lower contribution towards growing stock. However, sustainable utilization is necessary since these individuals represent future carbon stock security. For further enhancing the carbon security, the preservation as well as supplementing the natural regeneration is pertinent. For conservation of present carbon stock, the large trees require special focus in the region. Tagging and enumeration of these individuals will secure their conservation. The deadwood large trees are also vital ecosystem components and support essential ecological processes. Uniform implementation of forest management within an area is essential since there is heterogeneity in biomass distribution thus, creating biomass hotspots, and any disturbance especially in the peripheral areas of protected regions close to habitations could lead to significant loss. All management strategies require focus over both biomass and diversity enhancement, monospecific strategies are not a viable option. The management of these regions by forest authorities alone is not sufficient; the co-operation from local communities is a key necessity. Therefore, sensitization of local communities for sustainable use of forest resources along with enhancement of forest biodiversity is crucial. Further research is required for generalization of results, especially incorporating all the forest types along entire elevational gradients.
Additional file 1 Supplementary Table S1. GPS coordinates and elevational range of the quadrats laid under different forest types. Supplementary Table S2. Percentage distribution of stem density and total carbon density of different species in studied forest types. Supplementary Table S3