Open Access
Translator Disclaimer
1 May 2017 Allometric Bark Biomass Model for Daphne bholua in the Mid-Hills of Nepal
Ram P. Sharma, Shes K. Bhandari, Ram Bahadur BK
Author Affiliations +

Bark of Daphne bholua is an important non-timber forest product and makes a substantial contribution to the Nepalese economy. A precise estimate of the amount of D. bholua bark in mountain forests is possible using a biomass model. We developed an allometric bark biomass model for naturally grown D. bholua in Baglung District in the mid-hills of Nepal. The model was based on data from 101 destructively sampled D. bholua on 20 sample plots representing different growth stages (regeneration, established, and matured), site qualities, and stand densities, and we used diameter and height–diameter ratio as predictors. Among 9 functions evaluated, a simple power function showed the best fit to the data. This model described most of the variations in bark biomass with no substantial trends in the residuals. Leave-one-out cross-validation also confirmed the high precision of this model, because it described most of the variations in bark biomass with no substantial trends in the prediction errors. The model can be applied for a precise prediction of bark biomass for individuals of D. bholua with diameters and height–diameter ratios similar to those used in this study. It is site-specific, and its application should therefore be limited to sites with growth stage, site quality, stand density, and species distribution similar to those that formed the basis of this study. Further validation and verification of this model, with a larger dataset collected from sites with a wider range of these characteristics, is recommended.


Plant biomass can be estimated by using direct or indirect methods (Vogt et al 1998). Direct methods, traditionally preferred, require destructive sampling—that is, felling an individual and separating and weighing its parts. This method is the most accurate, but it is costly and time consuming. Generally, it is more suitable for small plants and sample sizes, which require fewer resources than similar testing for big trees and large sample sizes. Since forest area is often large, estimation of tree biomass for an entire forest using destructive sampling is almost impossible. An indirect approach using allometric models is the most appropriate option in such a situation (Chaturvedi et al 2010, 2012; Kuyah et al 2012; Subedi and Sharma 2012; Chaturvedi and Raghubanshi 2015). Allometric models serve as tools to relate variables such as diameter, height, height–diameter ratio, crown dimension, and wood density to biomass of the whole individual or its part. Felling and measurement of various components, including bark, of a number of sampled individuals of a plant species in any population is needed to establish a precise allometric biomass model, which can then be applied to predict biomass for the rest of the individuals in the population.

Because of diverse chemical and physical properties, tree bark could emerge as a potential raw material for various biomass-based industries. Bark offers a substantial added value to the biomass assortments and utilities. Because of the higher concentrations of several elements in the bark, it has higher metabolic activity than that of wood (Kenney et al 1990; Adler et al 2005). This is the reason that bark can be used for various purposes. For example, bark fibers from coppiced branches of Tilia cordata are used in Europe to make cordage for fishing and trap construction (Myking et al 2005). Bark of Betula species can be used as a covering in the making of canoes, as a drainage layer in roofs, and in the manufacture of shoes and backpacks (Adney and Howard 1964). Many commercial products are made from bark, including cork, cinnamon, and quinine; it is also used in tanning. Bark chips are used as mulch to protect soil moisture and nutrients (Adney and Howard 1964; Aronson et al 2009). The inner bark (phloem) of some species is also edible; for example, the Sami people of northern Europe use the bark of Scots pine (Pinus sylvestris) as a staple food (Zackrisson et al 2000). Similarly, bark of Cinnamomum tamala is used as a food flavoring (Subedi and Sharma 2012).

Non-timber forest products including bark play an important role in efforts to reduce the poverty of forest-dependent people, contributing to livelihoods and improvements in food security, income, and health (FAO 1995; Falconer 1997; Ahenkan and Boon 2008; Uprety et al 2010). Most of the people residing in rural parts of developing countries like Nepal depend heavily on non-timber forest products for food, medicine, and income (Olsen 1998; Olsen and Larsen 2003; Chhetri and Gupta 2006; Kunwar and Bussmann 2008; Piya et al 2011). These products also play an important role in the development of entrepreneurship and contribute substantially to local and national economies, especially in mountain areas (Shrestha et al 2003; Gauli and Hauser 2009; Uprety et al 2010).

Daphne bholua (family Thymelaeaceae), locally known as lokta, is an important non-timber forest product that grows naturally in the mid-hills of Nepal. It is an erect 1 to 3 m tall shrub found at elevations between 1600 and 4000 m (Khadgi et al 2013). Two species of Daphne are common in Nepal: D. bholua and D. papyracea (Khadgi et al 2013). Both species have strong fibers in their bark, which is used to make handmade Nepali paper; because of its higher fiber density, D. bholua is preferred (Paudyal 2004; Kharal et al 2011; Khadgi et al 2013). The handmade paper is used to prepare various value-added products, which have strong markets in Europe, the United States, and Japan (Banjara 2007). D. bholua has also been identified in Nepal's Master Plan for the Forestry Sector as 1 of the 7 non-timber forest products available in Nepal (MPFS 1989). D. bholua covers 2.9 million ha in 55 districts in Nepal, and total bark stock has been estimated at 110,481 metric tons, which could support paper production of more than 950 metric tons every year (FSRO 1984). Several enterprises that produce high-quality Nepali paper from the bark of D. bholua provide employment to local people (Paudyal 2004).

Realizing its economic importance, the Nepal government and some private enterprises have attempted to promote sustainable management of D. bholua (Kharal et al 2011). However, the total amount of available D. bholua bark remains unknown, because of lack of bark biomass models and scientific forest inventory methods suited to local conditions. Even though several studies have been conducted on Daphne species (Paudyal 2004; Biggs and Messerschmidt 2005; Banjara 2007; Kharal et al 2011; Khadgi et al 2013), none of them developed a bark biomass model. This study developed such a model using data from 101 destructively sampled D. bholua in Baglung District in Nepal's mid-hills. The proposed model can be useful for predicting the total amount of D. bholua bark available in a forest and can thus support effective management of Daphne species.

Material and methods

Study area

This study was conducted in the Dekhothulojanti Community Forest of Lekhani village development committee (VDC), Baglung District, Nepal (Figure 1). This district experiences subtropical to alpine climates, with mean annual rainfall ranging from 1500 to 3000 mm and mean annual temperature from 19 to 27°C. Most of the forests in this district, including Dekhothulojanti, are managed by local user groups. Dekhothulojanti covers 129.12 ha and is a natural, uneven-aged, mixed-species forest. The main tree species are Abies spp, Rhododendron spp, Quercus lamellosa, Quercus semecarpifolia, Aesculus indica, Alnus nepalensis, Saurauia nepalensis, Semecarpus anacardium, and Daphniphyllum himalense (DFO 2012). Daphne bholua or Nepalese paper plant, the species on which this study focused, grows throughout this forest.


Map of the study area and photo of the Dekhothulojanti Community Forest (Map of Nepal by Andreas Brodbeck; map of study area by Shes K. Bhandari; photo by Ram Bahadur BK).


Sampling and measurement

The forest was divided into 6 strata based on stand conditions, growth stages, and management modalities. Stand condition could be even-aged or uneven aged, monospecies or mixed species, and sparse or dense. We took 3 growth stages into account (defined in the next paragraph). Management modality was defined as management practices performed by user groups—ie, priority given either to timber or fodder production or biodiversity conservation. We applied stratified systematic line plot sampling (sampling intensity about 0.1%), but only in the forest area where D. bholua grows (effective area). The sampling design applied in this study is presented in Table 1. D. bholua was absent in the first stratum, therefore only 5 strata were used for allocation of sample plots. In total, 20 sample plots (5 × 5 m) were used, allocated to each stratum using the proportional allocation method (Chaturvedi and Khanna 2011). We avoided sampling damaged individuals. We identified 3 growth stages based on the diameter (at 30 cm above ground) and height of D. bholua individuals:

  • 1.  Regeneration was defined as a seedling with diameter <1 cm and height up to 1 m.

  • 2.  Established individuals were defined as those with diameter 1–3 cm and height >2 m or diameter >3 cm and height 1–2 m.

  • 3.  Matured individuals were defined as those with diameter > 3 cm and height > 2 m.


Sampling design used in the study.


For each sample plot, regeneration and established individuals were counted and 2 from each category were felled; matured individuals were counted and all were felled. The total height of each felled individual was measured. We felled 101 individuals. Because of the greater abundance of smaller (regeneration and established) individuals in each sample plot, the number of felled individuals in those categories was higher than that for matured plants, even though all matured individuals were felled (Table 2).


Summary statistics of the modeling data.


The stems of all felled individuals were debarked. Where applicable, large branches (>1 cm diameter) from established and matured individuals were also debarked and the bark added to the stem-bark samples. The total green weight of the bark of each felled individual was measured using a digital weighing machine with a precision of 10 mg. Of the 101 felled individuals, we randomly selected 75 and took 75% of the total bark sample of each to the laboratory, where it was dried in an oven at 80°C until a constant weight was achieved. The weight of the sample was recorded after 24 hours and again at 6-hour intervals. It took 30 hours to get constant weights for most of the samples, and only a few samples needed to be kept in the oven longer (30–36 hours) in order to get constant weights. From this, we derived an average dry-to-fresh-weight ratio, which we used to estimate the dry weight of the remaining 26 individuals. Summary statistics for these data are presented in Table 2, and the relation of bark biomass to diameter is displayed in Figure 2.


Relationship between diameter and bark biomass for D. bholua.


Model development

Plant biomass models commonly treat biomass as a function of

We plotted height–diameter ratio (HDR) against diameter (Figure 3) and used the results to inform our modeling. The HDR is also known as a tree slenderness quotient and is defined as total height divided either by diameter at breast height or diameter at the base of tree or other smaller individual. As HDR is a measure of slenderness, it is an appropriate characteristic to describe the form of an individual plant (Opio et al 2000; Sharma et al 2016). We developed and compared models based on the first 2 approaches, with the second approach adjusted to use HDR instead of height.


Relationship between HDR and diameter for D. bholua.


In preliminary analysis, we fitted 9 different functions (Table 3) applying the first approach (modeling bark biomass as a function of diameter alone). However, the resulting models performed poorly, and this approach was dropped from further evaluation. Then we applied the second approach using the parameter prediction method (Clutter 1983; Chapagain et al 2014), in which parameter b1 of each function in Table 3 was modeled as a nonlinear function of HDR. A simple power function (F1 in Table 3) with b1 modeled as a nonlinear function of HDR best fitted the data and therefore was selected (Equation 1) for further evaluation.

where Wi = bark biomass of individual i (g); Di = diameter of individual i (cm); HDRi = height-diameter ratio of individual i (m/cm); b1, b2, α1, α2 = parameters; and ℰi = residual error for individual i, which was assumed to be independent with zero mean and variance σ2.


Candidate functions used to fit data.a)


Model estimation and evaluation

We estimated models using PROC MODEL in SAS (SAS Institute 2012). The models were evaluated using root mean squared error (RMSE) and adjusted coefficient of determination (R2adj) (Montgomery et al 2001). We examined graphs of residuals against the fitted biomass and predictors. We used a 1% level of significance in our analyses.

Cross-validation increases the reliability and credibility of models; a common method of doing this is use of split data (Vanclay 1994; Kozak and Kozak 2003; Yang et al 2004). However, our dataset was too small for this to be feasible. Instead, we evaluated our model by applying the leave-one-out cross-validation (LOOCV) method, which is considered suitable for small datasets (Molinaro et al 2005; Nord-Larsen et al 2009; Timilsina and Staudhammer 2013; Rushing et al 2015). This involved exclusion of 1 observation from the full dataset, fitting the model to the remaining data, and predicting the biomass of the left-out observation. This was replicated for all 101 observations. The RMSE and R2 were calculated using the difference between predicted and observed biomasses; graphs of the mean prediction errors were also examined.


Among 9 functions evaluated, a simple power function (Equation 1) emerged as the best, describing a larger part of the variations in bark biomass of D. bholua (R2adj = 0.8873, RMSE = 20.98) than the other functions (Table 3). Inclusion of HDR significantly improved the fit of statistics over that of the model with diameter as a single predictor. Without inclusion of HDR, the model fitted relatively poorly (R2adj = 0.8516, RMSE = 25.39). All other functions fitted with diameter as single predictor also showed poorer fit statistics than those of a simple power function (F1 in Table 3). All parameter estimates of each function in Table 3 were significant (p < 0.01) and biologically plausible. The allometric bark biomass model with estimated parameters is


All symbols and abbreviations are the same as in Equation 1.

No heteroscedasticity problem in the residuals was observed. Also, no serious systematic trend in the residuals was observed within the ranges of observed diameters and heights of the individuals. The LOOCV did not show substantial trends in the prediction errors for the observed ranges of diameter and height (Figure 4). The model described a large part of the variations in the observed biomass of the left-out individuals in the LOOCV (R2 = 0.8653, RMSE = 23.03), and all parameter estimates from each iteration in the LOOCV were significant (p < 0.01). The model predicted the biomass data for larger individuals less precisely than that for smaller individuals. The bark biomass curves simulated using Equation 2 adequately covered the observed data with a few exceptions (Figure 5).


Mean residuals and mean prediction errors of bark biomass model by size class. Mean residuals were produced with the model, fitted to the full dataset. Mean prediction errors were produced with the model and fitted to the left-out data in the leave-one-out cross-validation.



Bark biomass curves produced with the model (Equation 2), overlaid on the observed biomass of D. bholua. The curves were produced using HDRs at 0.2 intervals, with the highest and lowest curves belonging to biomass for individuals with HDR of 0.2 and 1.6, respectively. The dots represent the same observed values shown in Figure 2.



Even though the model described a large part of the bark biomass variations, some residuals seemed larger, but were distributed randomly. The absence of serious trends in the residuals and prediction errors within the observed range confirms the high precision of the model (Figure 4). As in other studies (Subedi and Sharma 2012; Chapagain et al 2014), the model we developed using diameter as a single predictor performed poorly compared to the model that combined diameter and HDR and therefore could not be applied for precise prediction of all individuals in a stand. This is because D. bholua individuals with similar diameters within the same stand may have different heights, HDRs, and bark densities and consequently different bark biomasses. Variation in HDR was much greater for small diameters than for large diameters (Figure 3). This might be due to the sampling design, whereby individuals were selected from stands with large variations in competition. Greater competition results in thinner and taller individuals—that is, larger HDRs (Nykänen et al 1997; Opio et al 2000; Sharma et al 2016). To describe more variation in the bark biomass, inclusion of HDR as an additional predictor is thus justifiable.

A clear differentiation of the curves within the observed data range, even for individuals with similar diameters, is due to the large effect of HDR on the model. Because of differing individual heights (Table 2), their predicted biomasses are expected to be different, even for individuals with similar diameters, as shown in Figure 5. Adequate covering of the observed data by simulated curves also suggests that our model is precise for all possible sizes of individuals. Smaller variations of the residuals and prediction errors for smaller individuals (Figure 4) suggest that our model could be more precise for smaller individuals. Model precision depends largely on sample size, the predictors used, and application of the model to growth conditions of individuals. It may not be relevant to compare biomass models developed for different species, components, and locations, because several factors determine the accuracy of the model. Our model is as precise as a recent bark biomass model for Cinnamomum tamala (Subedi and Sharma 2012), which included both diameter and height as predictors. Our sample size (101 individuals) is larger than those used by other recent biomass studies (Ajit et al 2011; Sharma 2011; Subedi and Sharma 2012; Chaturvedi and Raghubanshi 2013; Chapagain et al 2014), which ranged from 27 to 46 individuals.

Total bark biomass of D. bholua in the forest can be estimated using our model. This information can be useful for sustainable management and effective commercial utilization of this economically important species. Several enterprises are producing high-quality handmade paper from the bark of D. bholua and providing employment to people residing in the mid-hills of Nepal (Paudyal 2004). Our model can also serve as a tool for predicting the total quantity of the raw material (bark) in the forest, which will directly contribute to feasibility studies for new enterprises to be established to meet growing demand for handmade paper, which is used in the preparation of various value-added products that have strong markets in Europe, the United States, and Japan (Banjara 2007).

Our model is site-specific and therefore may not be applicable to a wider geographical range, where growth conditions for D. bholua may differ. To make the model more comprehensive, accurate, and broadly applicable, it needs to incorporate measurements of site quality (eg site index), stand density (eg competition index), topographic features such as aspect, slope, and altitude; and climate and soil features such as temperature, precipitation, and soil properties—elements for which we did not have data.

Furthermore, we did not validate our model with data from other sites and regions, which would be necessary to confirm whether the model can be used for precise prediction under different growing conditions. Several studies have found that validating a model by applying LOOCV can yield less reliable results than validating it against external data (Molinaro et al 2005; Nord-Larsen et al 2009; Timilsina and Staudhammer 2013; Rushing et al 2015). Collecting a large amount of biomass data or data external to the main study site is often difficult due to lack of resources and time (Vanclay 1994). With a few exceptions (eg Brown et al 1993; Chave et al 2005), biomass modeling studies requiring destructive sampling often use small sample sizes from a small geographic area.


An allometric bark biomass model was developed using HDR and diameter of D. bholua individuals. The model described most of the variations in bark biomass (R2adj = 0.8873, RMSE = 20.98) with no substantial trends in the residuals. LOOCV confirmed the high precision of the model, because it described most of the variations in bark biomass (R2adj = 0.8653, RMSE = 23.03) with no substantial trends in the prediction errors.

Since site quality measures (eg site index), stand density measures (eg competition index), and measures of other factors that influence growth and development of D. bholua were not included in the model, it is site-specific. For high precision in prediction, application of the model should therefore be limited to site, growth stage, and stand conditions similar to those that formed the basis of this study.

The advantage of this model is that its application only requires measurement of outside bark diameter (0.3 m above ground) and total height of individuals to predict the total bark biomass of D. bholua in a forest, information which can then be used in formulation of a management plan. Further research is recommended to validate and verify our model using a larger dataset with a wider range of values for site quality, stand density, growth stage, and species distribution.


This publication was funded by Faculty of Forestry and Wood Sciences in Czech University of Life Sciences Prague through the funding scheme for research and publication (CZ60460709; IGA B-03/17). The District Forest Office, Baglung, provided financial support for inventory work. We are also thankful to the user groups of the Dekhothulojanti Community Forest for their support for inventory work. We are grateful to 2 anonymous reviewers, the MRD Associate Editor, and the external editor, who provided important comments and suggestions to improve this article.



Adler A, Verwijst T, Aronsson P. 2005. Estimation and relevance of bark proportion in a willow stand. Biomass and Bioenergy 29(2):102–113. Google Scholar


Adney T, Howard IC. 1964. The Bark Canoes and Skin Boats of North America. Washington, DC: Smithsonian Institution. Google Scholar


Ahenkan A, Boon EK. 2008. Enhancing Food Security and Poverty Reduction in Ghana through Non-timber Forest Products Farming: Case Study of Sefwi Wiawso District. Munich, Germany: GRIN Publishers. Google Scholar


Ajit, Das DK, Chaturvedi OP, Jabeend N, Dhyani SK. 2011. Predictive models for dry weight estimation of above and below ground biomass components of Populus deltoides in India: Development and comparative diagnosis. Biomass Bioenergy 35:1145–1152. Google Scholar


Alvarez E, Duque A, Saldarriaga J, Cabrera K, Salas GD, Valle I, Lema A, Morenob F, Orregob S, Rodrıguez L. 2012. Tree above-ground biomass allometries for carbon stocks estimation in the natural forests of Colombia. Forest Ecology and Management 267:297–308. Google Scholar


Aronson J, Pereira JS, Pausas JG. 2009. Cork Oak Woodlands on the Edge: Conservation, Adaptive Management, and Restoration. Washington, DC: Island Press. Google Scholar


Banjara GB. 2007. Handmade Paper in Nepal: Upgrading with Value Chain Network. Kathmandu, Nepal: German Technical Cooperation/Private Sector Promotion. Google Scholar


Basuki TM, van Laake PE, Skidmore AK, Hussin YA. 2009. Allometric equations for estimating the above-ground biomass in tropical lowland Dipterocarp forests. Forest Ecology and Management 257:1684–1694. Google Scholar


Bertalanffy LV. 1957. Quantitative laws in metabolism and growth. Quarterly Review of Biology 32:217–231. Google Scholar


Biggs S, Messerschmidt D. 2005. Social responsibility in the growing handmade paper industry of Nepal. World Development 33:1821–1843. Google Scholar


Brown S, Iverson LR, Prasad A, Liu D. 1993. Geographical distributions of carbon in biomass and soils of tropical Asian forests. Geocarto International 8:45–59. Google Scholar


Chapagain TR, Sharma RP, Bhandari SK. 2014. Modeling above-ground biomass for three tropical tree species at their juvenile stage. Forest Science and Technology 10(2):51–60. Google Scholar


Chaturvedi AN, Khanna LS. 2011. Forest Mensuration and Biometry. 5th edition. Dheradun, India: Khanna Bandhu. Google Scholar


Chaturvedi RK, Raghubanshi AS. 2013. Above-ground biomass estimation of small diameter woody species of tropical dry forest. New Forest 44(4):509–519. Google Scholar


Chaturvedi RK, Raghubanshi AS. 2015. Allometric model for accurate estimation of above ground biomass of teak in tropical dry forest of India. Forest Science 61(5):938–949 Google Scholar


Chaturvedi RK, Raghubansi AS, Singh JS. 2010. Non-destructive estimation of tree biomass by using wood specific gravity in the estimator. National Academy of Science Letters 33(5–6):133–138. Google Scholar


Chaturvedi RK, Raghubanshi AS, Singh JS. 2012. Biomass estimation of dry tropical woody species at juvenile stage. Scientific World Journal1–5.  10.1100/2012/790219Google Scholar


Chave J, Andalo C, Brown S, Cairns MA, Chambers JQ, Eamus D, Folster H, Fromard F, Higuchi N, Kira T, Lescure JP et al. 2005. Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 145(1):87–99. Google Scholar


Chhetri HB, Gupta VNP. 2006. NTFP potential of Upper Mustang: A trans-Himalayan region in western Nepal. Scientific World 4(4):38–43. Google Scholar


Clutter JL, Fortson JC, Pienaar LV, Brister GH, Bailey RL. 1983. Timber Management: A Quantitative Approach. New York, NY: Wiley. Google Scholar


DFO [District Forest Office]. 2012. Program glimpse of district forest office, Baglung. Baglung, Nepal: DFO. Google Scholar


Falconer J. 1997. Non-timber Forest Products in Southern Ghana: A Summary Report. ODA Forestry Series No 2. Chatham Maritime, United Kingdom: Natural Resources Institute. Google Scholar


FAO [Food and Agricultural Organization]. 1995. Non-wood Forest Products for Rural Income and Sustainable Forestry. Rome, Italy: FAO. Google Scholar


FSRO [Forest Survey and Research Office]. 1984. Preliminary Survey Report of Daphne spp. Publication No 41. Kathmandu, Nepal: Forest Survey and Research Office, Department of Forests. Google Scholar


Gauli K, Hauser M. 2009. Pro-poor commercial management of non-timber forest products in Nepal's community forest user groups: Factors and success. Mountain Research and Development 29(4):298–307. Google Scholar


Hosoda K, Iehara T. 2010. Above-ground biomass equations for individual trees of Cryptomeria japonica, Chamaecyparis obtusa and Larix kaempferi in Japan. Journal of Forest Research 15:299–306. Google Scholar


Huxley JS, Teissier G. 1936. Terminology of relative growth. Nature 137:780–781. Google Scholar


Kenney WA, Sennerby-Forsse L, Layton P. 1990. A review of biomass quality research relevant to the use of poplar and willow for energy conversion. Biomass 21:163–188. Google Scholar


Ketterings QM, Coe R, van Noordwijk M, Ambagau Y, Palm CA. 2001. Reducing uncertainty in the use of allometric biomass equations for predicting above-ground tree biomass in mixed secondary forests. Forest Ecology and Management 146:199–209. Google Scholar


Khadgi N, Shrestha BB, Siwakoti M. 2013. Resource assessment and habitat analysis of Daphne bholua in Bhujung of Annapurna Conservation Area, central Nepal. Research Journal of Agriculture and Forestry Science 1(4):1–9. Google Scholar


Kharal DK, Oli BN, Poudel I. 2011. Distribution and availability of raw materials for production of Nepali handmade paper from Daphne species in Darchula district, Nepal. Banko Jankari 21(2):24–30. Google Scholar


Kozak A, Kozak R. 2003. Does cross validation provide additional information in the evaluation of regression models? Canadian Journal of Forest Research 33(6):976–987. Google Scholar


Kunwar RK, Bussmann RW. 2008. Ethnobotany in the Nepal Himalaya. Journal of Ethnobiology and Ethnomedicine 4(24).; accessed 25 February 2016. Google Scholar


Kuyah S, Dietz J, Muthuri C, Jamnadass R, Mwangi P, Coe R, Neufeldt H. 2012. Allometric equations for estimating biomass in agricultural landscapes: I. Above ground biomass. Agriculture, Ecosystem and Environment 158:216–224. Google Scholar


Molinaro AM, Simon R, Pfeiffer RM. 2005. Prediction error estimation: A comparison of resampling methods. Bioinformatics 21:3301–3307. Google Scholar


Montgomery DC, Peck EA, Vining GG. 2001. Introduction to Linear Regression Analysis. New York, NY: Wiley. Google Scholar


MPFS [Master Plan for the Forestry Sector]. 1989. Master Plan for the Forestry Sector. Kathmandu, Nepal: Government of Nepal. Google Scholar


Myking, T, Hertzberg A, Skrøppa T. 2005. History, manufacture and properties of lime bast cordage in northern Europe. Forestry 78(1):65–71. Google Scholar


Nord-Larsen T, Meilby H, Skovsgaard JP. 2009. Site-specific height growth models for six common tree species in Denmark. Scandinavian Journal of Forest Research 24:194–204. Google Scholar


Nykänen ML, Peltola H, Quine C, Kellomäki S, Broadgate M. 1997. Factors affecting snow damage of trees with particular reference to European conditions. Silva Fennica 31:193–213. Google Scholar


Olsen CS. 1998. The trade in medicinal and aromatic plants from central Nepal to northern India. Economic Botany 52(3):279–292. Google Scholar


Olsen CS, Larsen HO. 2003. Alpine medicinal plant trade and Himalayan mountain livelihood strategies. Geographical Journal 169:243–254. Google Scholar


Opio C, Jacob N, Coopersmith D. 2000. Height to diameter ratio as a competition index for young conifer plantations in northern British Columbia, Canada. Forest Ecology and Management 137:245–252. Google Scholar


Paudyal AS. 2004. Sustainability of local handmade paper enterprises. Journal of Forest and Livelihood 4(1):64–69. Google Scholar


Piya L, Maharjan KL, Joshi NP, Dangol DR. 2011. Collection and marketing of non-timber forest products by Chepang community in Nepal. Journal of Agricultural Environment 12:10–21. Google Scholar


Rizvi RH, Gupta VK, Ajit . 2008. Comparison of various linear and non-linear functions for estimating biomass and volume of Dalbergia sissoo grown under rainfed conditions. Indian Journal of Agriculture Science 78:138–141. Google Scholar


Rushing C, Bulusu A, Hurwitz HI, Nixon AB, Pang H. 2015. A leave-one-out cross-validation SAS macro for the identification of markers associated with survival. Computers in Biology and Medicine 57:123–129. Google Scholar


SAS Institute. 2012. SAS/ETS1 9.3 User's Guide. Cary, NC: SAS Institute. Google Scholar


Schumacher FX. 1939. A new growth curve and its application to timber yield studies. Journal of Forestry Research 37:819–820. Google Scholar


Senelwa K, Sims REH. 1997. Tree biomass equations for short rotation eucalypts grown in New Zealand. Biomass and Bioenergy 13(3):133–140. Google Scholar


Sharma RP. 2011. Allometric models for total-tree and component-tree biomass of Alnus nepalensis D. Don in Nepal. Indian Forester 137:1386–1390. Google Scholar


Sharma RP, Vacek Z, Vacek S. 2016. Modeling individual tree height to diameter ratio for Norway spruce and European beech in Czech Republic. Trees 30:1969–1982. Google Scholar


Shrestha KK, Tiwari NN, Rajbhandari S, Shrestha S, Uprety Y, Poudel RC. 2003. Non-timber Forest Products in the Critical Bottlenecks and Corridors of Terai Arc Landscape: Documentation, Utilization, Trade and People's Livelihood. Kathmandu, Nepal: WWF-Nepal. Google Scholar


Singh V, Tewari A, Kushwaha SPS, Dadhwal V. 2011. Formulating allometric equations for estimating biomass and carbon stock in small diameter trees. Forest Ecology and Management 261:1945–1949. Google Scholar


Subedi M, Sharma RP. 2012. Allometric biomass models for bark of Cinnamomum tamala in mid-hill of Nepal. Biomass and Bioenergy 47:44–49. Google Scholar


Timilsina N, Staudhammer CL. 2013. Individual tree-based diameter growth model of Slash pine in Florida using nonlinear mixed modeling. Forest Science 59:27–37. Google Scholar


Uprety Y, Boon EK, Paudel RC, Shrestha KK, Rajbhandary S, Ahenkan A, Tiwari NN. 2010. Non-timber forest products in Bardiya district of Nepal: Indigenous use, trade and conservation. Journal of Human Ecology 30(3):143–158. Google Scholar


Vanclay JK. 1994. Modelling Forest Growth and Yield: Applications to Mixed Tropical Forests. Wallingford, United Kingdom: CAB International. Google Scholar


Vogt KA, Vogt DJ, Bloomfield J. 1998. Analysis of some direct and indirect methods for estimating root biomass and production of forests at an ecosystem level. Plant Soil 200:71–89. Google Scholar


Yang Y, Monserud RA, Huang S. 2004. An evaluation of diagnostic tests and their roles in validating forest biometric models. Canadian Journal of Forest Research 34(3):619–629. Google Scholar


Zackrisson O, Östlund L, Korhonen O, Bergman I. 2000. The ancient use of Pinus sylvestris L. (Scots pine) inner bark by Sami people in northern Sweden, related to cultural and ecological factors. Vegetation History and Archaeobotany 9(2):99–109. Google Scholar
© 2017 Sharma et al. This open access article is licensed under a Creative Commons Attribution 4.0 International License ( Please credit the authors and the full source.
Ram P. Sharma, Shes K. Bhandari, and Ram Bahadur BK "Allometric Bark Biomass Model for Daphne bholua in the Mid-Hills of Nepal," Mountain Research and Development 37(2), 206-215, (1 May 2017).
Received: 1 February 2017; Accepted: 1 March 2017; Published: 1 May 2017

Allometric modeling
destructive sampling
height–diameter ratio
leave-one-out cross-validation
non-timber forest product
Get copyright permission
Back to Top