Predicting the upper stem diameters and volume of a tropical dominant tree species

Sunita Ulak · Keshav Ghimire · Rabindra Gautam · Shes Kanta Bhandari,3 ·Krishna Prasad Poudel · Yajna Prasad Timilsina · Dhirendra Pradhan ·

Thakur Subedi1

Abstract Predicting upper stem diameters and individual tree volumes is important for product quantif ication and can provide important information for the sustainable management of forests of important commercial tree species ( Shorea robusta) in Nepal. The aim of this study was to develop a taper equation for S. robusta. Fifty-four trees were selected and felled in the southern low land of Nepal.A mixed eff ect modelling approach was used to evaluate 17 diff erent taper functions. ‘Leave-one-out cross-validation’was used to validate the f itted taper functions. The variable exponent taper function best f itted our data and described more than 99% of the variation in upper stem diameters.Results also showed signif icant eff ects of stand density on tree taper. Individual tree volume prediction using the local volume model developed in this study was more accurate compared to the volume predicted through the taper function and existing volume model. The taper function developed in this study provides the benef it of predicting upper stem diameter and can be used for predicting volume to any merchantable height of individual trees. It will have implications in estimates of volume, biomass, and carbon and thus may be a potential supporting tool in carbon trade and revenue generation.

Keywords Mixed-model approach · Shorea robusta ·Taper · Volume · Nepal

Introduction

Measurement of individual tree attributes is important for estimating growth and yield as it provides guidelines for decision- making for sustainable forest management (Dau et al. 2015). However, the measurement of all individual tree variables is costly and time-consuming. Diameter at breast height (DBH) and the total height of individual trees are commonly measured variables (Sumida et al. 2013).Compared to DBH, the total height is less frequently measured. Generally, height of sample trees is measured and the height of the rest of the stand is predicted using heightdiameter allometry (Peng et al. 2001; Zhang et al. 2014;Bhandari et al. 2021a, 2021b, 2021c). Although DBH can be measured easily, measurement of the upper stem diameter at a particular height is diffi cult. Taper equations allow for the prediction of upper stem diameters using easily measured variables such as DBH and total height (Kozak and Smith 1966; Kozak et al. 1969; Kozak 2004; Poudel et al. 2018). A taper equation can directly, or through integration, estimate (1) diameter at any point along the stem,(2) the height at a given diameter, (3) total stem volume,(4) merchantable volume and merchantable height to any top diameter from any stump height, and (5) volumes for logs of any length at any height from the ground (Kozak 2004; Poudel et al. 2018).Taper in a tree stem is the rate of decrease in stem diameter from base upwards (Chaturvedi and Khanna 2011) and therefore depends on DBH, total height (H) and upper stem height (h). Taper functions vary in their mathematical form; however, a simple taper function can be represented by Eq. 1.

wherediis the stem diameter at heighthiabove the ground for a tree with a diameter at breast height (DBH) and total tree height (H).

The study of tree taper has a long history. Kozak and Smith ( 1966) and Kozak et al. ( 1969) recommended a simple quadratic model for describing tree taper. Bruce et al.( 1968) developed a polynomial model for predicting upper stem diameter of red alder (Alnus rubraBong.). Bennett and Swindel ( 1972) also used a polynomial equation for the upper stem diameter prediction of slash pine (Pinus elliottiiEngelm.). All these taper functions assumed that an entire stem has a single form that can be described with a single function. Consequently, these were referred to as single taper functions (Tang et al. 2017). However,the form of a tree stem varies even within a single tree.The basal portion of the stem resembles the frustum of a neiloid, the middle portion with a paraboloid, and the top portion with a conoid (Husch et al. 1972). Realizing the requirement for three diff erent models for three different portions of the stem of an individual tree, Max and Burkhart ( 1976) developed segmented taper functions for natural and planted loblolly pines (Pinus taedaL.). The segmented taper functions are more accurate than the single taper function because they consider the variation in tree form at diff erent heights.

Tree form varies among species, individual trees of the same species, and within and between stands as well.To account for such variations, Newnham ( 1988) and Kozak ( 1988) introduced a variable-exponent taper function which uses changing exponents to compensate for the changes in the form at diff erent sections of the stem(Perez et al. 1990). Several studies have shown that variable exponent taper functions predict upper stem diameters more accurately than the other two types of taper functions (Sharma and Zhang 2004; Rojo et al. 2005; Duan et al. 2016). However, other studies have also reported that better performance of segmented taper functions over variable exponent taper functions (Crecente-Campo et al. 2009; ?z?elik and Crecente-Campo 2016; Shahzad et al. 2021). Several studies have been carried out on tree tapers in diff erent areas of the world (Candy 1989; Fang and Bailey 1999; Kozak 2004; Sharma and Patron 2009;Arias-Rodil et al. 2015; Poudel et al. 2018; Silwal et al.2018; Alkan and ?z?elik 2020; Shahzad et al. 2021). Most of studies have taken place in North America (Max and Burkhart 1976; Shaw et al. 2003; Sharma and Patron 2009;Li and Weiskittel 2010; Poudel et al. 2018), in Europe(Arias-Rodil et al. 2015; Socha et al. 2020), Asia (Fang and Bailey 1999; Tang et al. 2016; Lumbres et al. 2017;Silwal et al. 2018; Shahzad et al. 2021; Koirala et al. 2021)and relatively a lower number in Australia (Gordon 1983;Bi and Turner 1994; Bi 2000), South America (Nunes and G?rgens 2016; Beltran et al. 2017) and Africa (Mabvurira and Eerik?inen 2002; Gomat et al. 2011). To our knowledge, only two taper functions have been developed in Nepal and these were developed forS. robustaC.F.Gaertn. (Silwal et al. 2018) andTectona grandisL. f. (Koirala et al. 2021).

Tree taper depends on species, site, age, and stand density (Sharma and Zhang 2004; Chaturvedi and Khanna 2011; Duan et al. 2016). These variations require diff erent taper equations and are more important in forests with high degree of species diversity, various management regimes, stand density, and human inf luence. The tropical forests in the western low land of Nepal are rich in diverse species and are managed under diff erent regimes(community forest, leasehold forest, collaborative forest,religious forest, national forest).S. robusta(family Dipterocarpaceae) is one of the dominant and more important timber species in the tropical forests of Nepal. Its distribution ranges from 120 m in the southern lowlands to 1200 m in the middle hills (Jackson 1994; Sah 2000).S. robustacontributes 19.28% (31.76 m 3 ha -1 ) in total standing volume and 15.27% in forest type coverage in Nepal (DFRS 2015).S. robustais a multipurpose species as its stem is used as timber, construction material, and fuelwood (Jackson 1994), leaves as fodder (Kibria et al. 1994) and plates(Kora 2019), and resin as medicine for dysentery and gonorrhea (Joshi 2003). S. robustaforests have the highest total standing volume in Nepal. However, studies in Nepal on the species are limited and mostly focused on volume models (Sharma and Pukkala 1990; Subedi 2017), biomass(Chapagain et al. 2014; Bhandari and Chhetri 2020), form factors (Baral et al. 2020; Chapagain and Sharma 2021;Subedi et al. 2021), regeneration dynamics and diversity(Awasthi et al. 2015, 2020). There is only one study by Silwal et al. ( 2018) on predicting upper stem diameters ofS. robusta. However, their function tested only one form of taper equation, did not consider the eff ect of stand density,and was based only on 33 individual trees.

The aim of this study was to develop models to predict upper stem diameter and volume of a dominant tree species(S. robusta). Our research questions were: (1) what is the best function to describe tapering in individual trees ofS.robusta? (2) does stem taper vary with stand density? (3) is the volume predicted using the taper function more accurate than that predicted by the local volume model? and, (4) is the volume predicted using the taper function more accurate than that predicted by volume models currently used in Nepal?

Materials and methods

Study area

This study used data from a forest of Jogikuti, Shankar Nagar (27°42’ N, 83°28’ E), one of the research sites maintained by the Forest Research and Training Centre, (previously known as Department of Forest Research and Survey)of the Government of Nepal, 25 km northeast of Bhairahawa Airport, Nepal (Fig. 1). The forest was naturally regenerated in 1986 and is dominated byS. robustaor sal. Field measurements were carried out in 2019 when the forest was 33 years old. After a decade of regeneration, stand density was an average of 5200 trees ha -1 (Hartz 1997). The soil is loamy and deep with adequate nutrients. The area has a tropical climate with a long-term average annual rainfall of 1649 mm, 85% from June to September (CBS 2019). Longterm rainfall data shows a declining trend. Mean monthly maximum temperature is 31.1 °C and the mean monthly minimum temperature is 19.2 °C (CBS 2019). The study site is mostly f lat with zero to 5% slope.

Experimental design

The research plot covers an area of 92 ha of natural sal forest and is intended to generate information on the growth potential and biomass production ofS. robustaunder different thinning regimes. The experiment was designed in a randomized complete block design (RCBD) with four treatments and three replications. One of the four treatments was a control; data from the control was not included because the trees were not harvested and taper measurements were unavailable. The remaining three treatments were “natural forest with the highest number of stems ha -1 (500 trees ha -1 , T1)”, “natural forest with average number of stems ha -1 (400 trees ha -1 , T2)” and“natural forest with the lowest number of stems ha -1 (300 trees ha -1 , T3)” (Table 1). The RCBD design minimized any existing unidirectional variability in the study area.Each plot was 100 m × 100 m (1 ha), whereas each block covered an area of 4.7 ha. Spacing between the treatments was a minimum of 5 m buff er on each side. The buff er was also used as a f ire line and to make the plot more accessible for silvicultural operations. Spacing between the blocks was at least 50 m.

Fig. 1 Map of Nepal showing the study site

Table 1 Variables used in this study; diameter at breast height (DBH), upper stem diameter (d), height, and volume; number of trees ha -1 is a stand level variable at 17, 22 and 33 years

Sampling and measurement

Fifty-four trees (18 from each treatment) were selected for destructive sampling. The individual trees from each treatment were selected, based on two main criteria. The f irst was the individual tree must be dominant or co-dominant,the second was the selected dominant or co-dominant trees must be free from disease, decay, and defects on a physical inspection. Once selected, DBH was measured and the tree felled at stump height (at 30 cm from the ground). Total height was measured with a linear tape, and diameter at 1,2.5, 5, 7.5, 10, 15, 20, 30, 40, 50, 60, 70, 80, and 90% of the total height.

Descriptive statistics of the data

The descriptive statistics of individual trees of all three treatments used in developing taper functions are presented in Table 1. Average DBH was 21.2 cm and minimum and maximum DBH were12.0 and 34.4 cm. Average height was 21.1 m, and minimum and maximum heights were 12.3 and 28.5 m. Average volume of the sample trees was 0.4 m 3 , and minimum and maximum volumes were 0.1 and 1.0 m 3 . The distribution of DBH and height showed that the stand was relatively young and had not reached its peak. The scatter plot (not included) of relative diameter and height showed a decreasing trend of diameter with an increase in height.

Model f itting and evaluation

Seventeen commonly used taper functions were f itted to our data and evaluated for their performance (Table 2). To account for the repeated measurement of diameter at different heights of the same tree, a non-linear mixed eff ect modelling approach was used. In this approach, DBH, total height (H), and height at which upper stem diameter was measured (h) were considered the f ixed eff ect parameters and the individual-tree variation was considered as random eff ects. Values of all the parameter and f it statistics were estimated using nlme package (Pinheiro et al. 2018) in Rversion4.0.5 (R Core Team 2021). Wherever available, the published values of parameters were used as the initial values in model f itting to facilitate the fast convergence. Orthogonal polynomials were used for predictor variables to minimize the multicollinearity (Bhandari et al. 2021a, 2021c).The f itted taper functions were evaluated using diff erent criteria, including the signif icance of estimated parameters (at 5% level of signif icance); coeffi cient of determination (R 2 ;higher values indicate better models) (Eq. 2); root mean squared error (RMSE; lower values indicate better models)(Eq. 3) (Montgomery et al. 2001); and Akaike informationcriterion (AIC; lower values indicate better models) (Eq. 4)(Akaike 1972; Burnham and Anderson 2002). The distribution of residuals was plotted against the predicted upper stem diameter to evaluate whether the taper function had a systematic bias. Residuals were also plotted against the predictor variables and examined for evidence of bias. The best three taper functions were selected using the above evaluation criteria.

Table 2 Diff erent taper models used in the study

Eff ect of stand density on taper

The method described by Bhandari et al. ( 2021b) was used to test the eff ect of stand density on stem tapering. A linear model of upper stem diameter was f itted as a function of height at which the upper stem diameter was measured (h),total height and DBH was allowed to vary with stand density by including an interaction: d ～ (h- + -H- + -DBH)*Tr,where Tr was the number of stems per ha (stand density). It was then tested whether the eff ect of stand density (Tr) was signif icant in the model.

Validation of taper models

The leave-one-out cross-validation (LOOCV) method was used to validate the selected taper function. In this method,one observation was left out for validation andn-1 observations were used for model f itting, wherenis the total number of observations used in data collection. This process was continued n times until all the observations were used as validation data. The predicted diameter on validation data was plotted against the observed diameter and the f it statistics including R 2 , RMSE, and mean absolute error (MAE)were estimated.

Volume prediction

Actual volume calculation

Diameter at diff erent proportions of total height and length between two consecutive diameter measurement points were recorded for each tree (See “ Sampling and measurement”section). From the diameter and length of each section, a separate volume was calculated. Formulae of a cylinder(Eq. 7) and cone (Eq. 9) were used to calculate the volume of the stump (30 cm height above the ground) and the top section of individual trees, respectively. Volumes of all other sections were calculated using Smalian’s formula (Eq. 8), as this is commonly used to estimate the volume of small sections accurately (Avery and Burkhart 1994; Chaturvedi and Khanna 2011; Beltran et al. 2017).

Taper based volume

Since all the taper equations used in this study do not have a closed form, numerical integration was used to calculate the stem volume of each tree. This method was also followed by Poudel et al. ( 2018). The total length of the stem was divided into 100 equal sections. The diameters at two ends of those sections were predicted using the f itted taper model, and the volume of each section was calculated using Smalian’s formula. Total volume of the tree was obtained by adding the volume of all 100 sections.

Local volume model developed in this study

The stem volume of each tree was modelled as a function of DBH and total height (Eq. 10). A commonly used power function in the logarithmic form was f itted as this is recommended by several authors for predicting the volume ofS.robustain Nepal (Sharma and Pukkala 1990; Subedi 2017).Logarithmic transformation was used because such a model produced smaller residuals than the nonlinear model. The local volume model developed in this study was also validated using LOOCV method.

Volume equation of Sharma and Pukkala ( 1990)

Sharma and Pukkala ( 1990) developed a model using DBH(D) and total height (H) as variables to predict the stem volume of individual trees ofS. robustain the Nepalese forest (Eq. 11).

For comparison, we predicted the volume using the local volume model developed in this study (Eq. 10), taper equation developed in this study (M1), and models of Sharma and Pukkala ( 1990) (Eq. 11). Bias (Eq. 5) and RMSE (Eq. 4)were estimated and compared using these equations. The predicted volume was then plotted against DBH.

Results

Taper functions

The parameter estimates and f it statistics (R 2 , RMSE and AIC) for the seventeen taper functions evaluated in this study are presented in Table 3. Except for β1of M16 and β2of M17,all other parameters of the 17 taper functions were signif icant at 95% conf idence level. Therefore, the taper functions M16 and M17were excluded from further analysis. Evaluation based on f it statistics showed that most of the remaining taper functions provided a strong f it except M15 which described only 80% of the variation in upper stem diameter with comparatively higher RMSE and AIC. Taper function M15 was also excluded from further analysis. The rest of the taper functions explained more than 96% of the variation in upper stem diameters. M1, M2, and M3 were selected as the best three as they provided higher values of R 2 and lower RMSE and AIC values than other remaining taper functions.The distribution of raw residuals against the predicted upper stem diameter and predictor variables (total height,height at which diameter was measured, and DBH) with a smooth superimposed curve for the three best taper functions (M1, M2 and M3) using tree-level random eff ect is presented in Fig. 2. Plot of residuals against the predicted upper stem diameter and predictor variables did not show any systematic trend in residual distribution. The taper function M1 performed the best with the highest R 2 (0.99),lowest RMSE (0.57), and lowest AIC (1829). The graph of taper function M1 showed a minimum error with negligible curvature and outliers than other taper functions, suggesting the absence of local bias in this taper function. Therefore,M1 was selected as the best taper function to predict the upper stem diameters of individual trees ofS. robusta. The predicted upper stem diameters at diff erent tree heights using three best taper functions are presented in Fig. 3 a. Stand density had a signif icant eff ect on the tapering of the stem ofS. robusta(P＜ 0.05).

Model validation

The f it statistics R 2 and RMSE of the model f itted with validation data was similar with that of the model f itting data(Fig. 3 b - d). The plot of predicted diameter against the observed diameter and predicted volume against observedvolume further demonstrated that the developed taper and volume models are validated and able to predict upper stem diameters and volume with an acceptable level of accuracy.

Table 3 Parameter estimates and f it statistics for diff erent taper models used in this study

Fig. 2 Distribution of residuals against the predicted upper stem diameter( a- c); total height( d- f); height at which diameter was measured ( g- i); diameter at breast height ( j- l) for the best three taper functions

Volume prediction

All parameter estimates of the f itted model were signif icant at a 99% conf idence level. The model described more than 99% variation in the total volume of the individual tree(Eq. 10). Validation of this model (Eq. 10) also showed a similar f it statistic (Fig. 4 a).

The volume calculated on destructive measurement was considered as a reference volume, and the volume predicted by the local volume model developed in this study, the taper equation developed in this study, and the model of Sharma and Pukkala ( 1990) were compared with the reference volume. Volume predicted by the local volume model and taper equation showed a positive bias, demonstrating that these two models underestimated the volume ofS. robusta. Volume predicted by the Sharma and Pukkala ( 1990) model had a negative bias, demonstrating an over prediction. The local volume model showed a signif icantly smaller bias and RMSE compared to the other models (Table 4) and indicates that this model is more accurate in volume prediction. The diff erence in predicted volume among the models is smaller for small-sized trees and increased with size (DBH) of thetree (Fig. 4 b).The predicted volume of individual trees differed signif icantly (P＜ 0.05) among three diff erent stand densities (300, 400 and 500 trees per ha) (Fig. 4 c, d).

Fig. 3 Predicted upper stem diameters at diff erent heights using best three taper functions for a Shorea robusta with DBH 32.7 cm and height 26.8 m ( a);leave-one-out cross-validation of taper model M1 ( b); M2 ( c);and M3 ( d). R 2 is coeffi cient of determination, RMSE is root mean squared error and MAE is mean absolute error. The dotted line in b- d is the 1:1 reference line

Discussion

Seventeen taper functions were f itted and evaluated using data collected from 54 felled trees and the best model was selected for predicting upper stem diameters at a given height ofS. robusta. The eff ect of stand density on taper was also examined. In addition, the actual volume of the individual sampled trees was compared with the volume predicted by the integration of the developed taper function, the local volume model and previously developed volume models forS. robustain Nepal. The main f indings indicate that the taper function ofS. robustadeveloped using the treelevel random eff ect in mixed-eff ect modelling performed well for upper stem diameter prediction. It was also found that stem taper varied signif icantly with stand density. The merchantable and/or total stem volume of an individual tree can be predicted by integration of the taper function, but the accuracy of the predicted volume is lower compared to the stem volume predicted by the local volume model.

This study showed that the variable-exponent taper function of Sharma and Zhang ( 2004) performed best in predicting upper stem diameters forS. robusta. This has also been reported as the best function in predicting upper stem diameters of Chinese f ir (Cunninghamia lanceolata(Lamb.)Hock.) in southern China (Duan et al. 2016). Several other studies have reported the better performance of variable exponent taper function in comparison to a single form and segmented taper functions (Sharma and Zhang 2004; Rojo et al. 2005; Duan et al. 2016). In contrast, other studies have shown that segmented taper functions are better than variable exponent taper functions (Crecente-Campo et al. 2009;?z?elik and Crecente-Campo 2016; Shahzad et al. 2021),suggesting that a particular form of the taper function may not be universally applicable.

Stand density has a signif icant eff ect on the tapering of individual trees ofS. robusta. Growth in DBH and height of individual trees varies with stand density (Bhandari et al.2021b). Diameter growth is more sensitive to a reduction in stand density than height. This indicates that diameter growth takes place more rapidly in a stand with low densitythan height. This results in variations in stem form and tapering of individual trees. Variation in stem form due to stand density may be explained by Metzger’s theory (Metzger 1893).Individual trees growing in an open stand allocate more resources towards the base of the tree to withstand external forces, including wind (Chaturvedi and Khanna 2011). Allocation of large amounts of resources towards the base of the tree results in a larger base and greater tapering.In contrast, trees growing in a dense stand face a relatively lower magnitude of external forces and avoid the unnecessary allocation of resources towards their base. Individual trees growing in dense stands do not receive suffi cient lightfor growth compared to trees growing in less denser stands.Because of the limited availability of light, trees growing in the dense stands have a tendency of growing upward to receive more light. This results in less tapering and a more cylindrical stem. This study also showed a signif icant diff erence in individual tree volumes based on stand density. The diff erence in individual tree volume between 400 trees ha -1 and 500 trees ha -1 was relatively small but the diff erence between 300 and 500 trees ha -1 was relatively large (Fig. 4 c). The diff erence between the average volume of trees in the stand with 300 trees ha -1 and with 500 trees ha -1 is statistically signif icant because the diff erence is ＞ 20% (Table 1). The diff erences observed is for stands at age 33 years. As the size and age of the stand increase, the diff erences in volumes are expected to increase. Trees growing in lower stand density yield higher timber than trees in higher density stands. This will have signif icant implications in decision- making for forest managers, especially in timber production industries.

Fig. 4 a Leave-one-out crossvalidation of local volume model; b comparison of the local volume model with that predicted from the taper equation and Sharma and Pukkala( 1990) model; c comparison of the predicted volume by the local volume model among three stand densities; and ( d)comparison of the predicted volume by the best taper equation among three stand densities; R 2 is coeffi cient of determination, RMSE is root mean squared error and MAE is mean absolute error; the dot line in a is the 1:1 reference line

Table 4 Mean bias and RMSE for predicting individual tree volume of Shorea robusta using volume model f itted in this study, taper equation, and Sharma and Pukkala ( 1990)

Taper functions for tree species of Nepal are limited.Only the taper function developed by Silwal et al. ( 2018)is available forS. robustawhich used 33 trees from two diff erent sites of the Terai region. Our study has used a relatively larger number of sample trees (n= 54). Data collected through destructive sampling is likely to improve accuracy more than data measured from the standing trees, although it increases the cost and time in data collection. Because of the higher accuracy of the data collected through destructive sampling, datasets generated from a small number of trees may be enough to develop the taper model and represent the whole population. Although the number of individual trees felled for data collection was 54, there were 756 data points to f it and evaluate the model. A mixed eff ect modelling approach was used to account for the correlation associated with multiple measurements taken from the same tree.The mixed eff ect modelling approach is more appropriate for repeatedly measured data (Lindstrom and Bates 1990).

This study showed that the volume predicted by the local volume model from the same forest and the general volume table developed for the whole of Nepal is more accurate than the volume predicted by the integration of the taper function. This indicates that the taper functions developed in this study are useful tools in predicting upper stem diameters but not in predicting the total stem volume of individual trees.Considering the poor performance of the taper function for total stem volume prediction, it is recommended using the local volume table for the prediction of total stem volume and taper function for the prediction of merchantable volume and/or volume up to a certain height for the stem and upper stem diameter. Poudel et al. ( 2018) also reported that the volume predicted by local volume models are more accurate than volumes predicted from taper functions developed for Douglas-f ir and Western hemlock in the Pacif ic Northwest of the USA. In a study carried out in Mongolian oak (Quercus mongolicaFisch. ex Ledeb.) in South Korea, Ko et al. ( 2019)found a signif icant diff erence between the volume predicted by existing volume models and taper functions. They recommended using new taper equations to predict the volume by assuming that the forest condition and structure had been signif icantly changed since the preparation of the existing volume model.

Taper functions can be used as a tool in quantifying timber and other forest products, either from an individual tree or from a stand. These functions also serve as decision support tools in carbon trading and other ecosystem service estimations. Predicting upper stem diameter and stem volume of individual trees provides crucial information to prepare management plans or to implement those plans for sustainable forest management. As the taper function varies with species, stand densities, and management regimes, it is suggested that the development of taper function for each species and management regimes would achieve more accurate prediction and useful information for the sustainable management of forests.

Conclusions

This study developed a taper function using a non-linear mixed eff ect model for predicting upper stem diameters ofS. robusta, the most important timber species in the southern lowlands of Nepal. The taper function of Sharma and Zhang ( 2004) performed the best in predicting upper stem diameter among the seventeen candidate functions studied.Stem taper ofS. robustavaried with stand density. Although the taper function facilitated the prediction of upper stem diameter, total tree volume prediction was not as accurate as that obtained from the local volume equation. The developed taper functions can be used as a tool for estimating volume,biomass and carbon and can provide important information for the sustainable management of the most important commercial tree species (S. robusta) in Nepal. The taper model can be further improved by using data from many individual trees and a wider range of geographical regions and management regimes.

Acknowledgements This study was part of a regular program of the Forest Research and Training Centre (FRTC), Nepal. We would like to acknowledge the FRTC for providing essential resources to accomplish this study and we extend our special thanks to Director General, Dr.Dipak Kumar Kharal and Deputy Director General, Dr. Buddhi Sagar Paudel for constructive suggestions and feedback. We would also thank Mr. Amrit Ghimire and Mr. Rohit Bhusal for their assistance during the f ield work.

Author contributions The concept of this study was designed by SU and KG. Data were by collected by SU, KG, TS and RG. SKB, KPP and YPT analysed the data. SKB and DP prepared the draft manuscript.All authors thoroughly read and agreed to the published version of the manuscript.

Funding Open Access funding enabled and organized by CAUL and its Member Institutions.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.