Next Article in Journal
C Stocks in Forest Floor and Mineral Soil of Two Mediterranean Beech Forests
Previous Article in Journal
Analysis of the Mechanism Hindering Sustainable Forestry Operations: A Case Study of Japanese Forest Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Allometric Equations for Estimating Tree Aboveground Biomass in Tropical Dipterocarp Forests of Vietnam

1
Department of Forest Resources and Environment Management, Tay Nguyen University, 567 Le Duan, Buon Ma Thuot, Dak Lak 630000, Vietnam
2
Department of Forest Engineering, Resources and Management, Oregon State University, Corvallis, OR 97333, USA
3
Forest Inventory and Planning Institute, Vinh Quynh, Thanh Tri, Ha Noi 100000, Vietnam
4
Vietnam National University of Forestry, Xuan Mai, Chuong My, Ha Noi 100000, Vietnam
5
Vietnamese Academy of Forest Sciences, Đuc Thang, Bac Tu Liem, Ha Noi 100000, Vietnam
*
Author to whom correspondence should be addressed.
Forests 2016, 7(8), 180; https://doi.org/10.3390/f7080180
Submission received: 23 May 2016 / Revised: 15 August 2016 / Accepted: 16 August 2016 / Published: 22 August 2016

Abstract

:
There are few allometric equations available for dipterocarp forests, despite the fact that this forest type covers extensive areas in tropical Southeast Asia. This study aims to develop a set of equations to estimate tree aboveground biomass (AGB) in dipterocarp forests in Vietnam and to validate and compare their predictive performance with allometric equations used for dipterocarps in Indonesia and pantropical areas. Diameter at breast height (DBH), total tree height (H), and wood density (WD) were used as input variables of the nonlinear weighted least square models. Akaike information criterion (AIC) and residual plots were used to select the best models; while percent bias, root mean square percentage error, and mean absolute percent error were used to compare their performance to published models. For mixed-species, the best equation was AGB = 0.06203 × DBH 2.26430 × H 0.51415 × WD 0.79456 . When applied to a random independent validation dataset, the predicted values from the generic equations and the dipterocarp equations in Indonesia overestimated the AGB for different sites, indicating the need for region-specific equations. At the genus level, the selected equations were AGB = 0.03713 × DBH 2.73813 and AGB = 0.07483 × DBH 2.54496 for two genera, Dipterocarpus and Shorea, respectively, in Vietnam. Compared to the mixed-species equations, the genus-specific equations improved the accuracy of the AGB estimates. Additionally, the genus-specific equations showed no significant differences in predictive performance in different regions (e.g., Indonesia, Vietnam) of Southeast Asia.

1. Introduction

There is a significant need globally and in Vietnam to develop equations that estimate forest biomass and carbon for national measuring, reporting, and verification systems. Under the United Nations Framework Convention on Climate Change, countries must report the state of their forest resources through initiatives such as the United Nations collaborative programme on Reducing Emissions from Deforestation and Forest Degradation Plus UN-REDD+ programme [1] in developing countries [2]. For high accuracy, equations specific to dominant forest types and ecological regions are necessary.
Dipterocarpaceae is one of the most well-known tree families in the tropics for its species’ iconic form and value as timber [3]. Timber species of Dipterocarpaceae are distributed widely throughout tropical Asia, spanning multiple climatic zones and geographic regions [4]. In Vietnam, this forest type covers approximately 650,000 ha and is distributed in two main ecological regions or ecoregions, namely the Central Highlands (CH) and the Southeast (SE) [5]. Differences in historical management of dipterocarps have resulted in an uneven balance of knowledge across their range, with much known about dipterocarps on the Indian subcontinent and comparatively little known about the dipterocarps of Southeast Asia [3]. In Vietnam, the ability to accurately estimate biomass and carbon in dipterocarp forests is complicated further by the mixed-species nature of this forest type [4].
Jenkins et al. [6,7] performed a modified national-level meta-analysis of allometric regression equations for total aboveground biomass (AGB) based on tree diameter and developed consistent equations for tree species in the United States using a logarithm function as the general form of these equations. Power models relating tree AGB to diameter at breast height (DBH) or more complex associations of input variables have also been widely developed and used for pantropic [8,9]. For example, Ter-Mikaelian and Korzukhin [10] used power models to estimate biomass for 65 North American tree species. In general, numerous publications suggest power models as a suitable option for building allometric equations based on one or more variables [8,11]. In developing biomass equations for pantropic or dipterocarp forests in the Southeast Asia, Brown [12], and Basuki et al. [13] applied logarithmic transformation of the power model. Picard et al. [14] pointed out that although being a good model for biomass prediction, the power model was outperformed by the allometric model.
The literature on allometric equations developed to estimate biomass for tropical forests is sizable [12,13,15,16,17,18,19,20,21]. However, there is limited literature and allometric equations for tropical dipterocarp forests [13] despite covering extensive areas in tropical Southeast Asia. To date, there are only a few general pantropical equations, but due to the range of diverse forest types with different characteristics in this region, there is a need to develop more specific biomass equations. Using data collected from 122 destructively sampled trees, Basuki et al. [13] developed allometric equations for dipterocarp forests, including one for each of the four main genera (Dipterocarpus, Hopea, Polaquium, and Shorea), one for commercial forests, and one for mixed-species of dipterocarps in East Kalimantan, Indonesia. The input variables of these equations were DBH, wood density (WD), and commercial bole height. Despite the development of equations for the main genera of the dipterocarp forest type, these authors focused on only one model form: a power model linearized by log transformation. Furthermore, they did not compare the predictive performances of genus or mixed-species allometric equations.
In Vietnam, there were a few allometric equations developed before the beginning of the UN-REDD programme [8,9]. Recently, Huy et al. [22] developed aboveground biomass equations for the evergreen broadleaf forests of the South Central Coastal ecoregion. The authors of this study report on the development of national-scale allometric equations to estimate tree biomass for the dipterocarp forest type in Vietnam, which was summarized and introduced by Sola et al. [8,9]. The objectives of this study were to: (i) develop equations to estimate AGB of dipterocarp forests in Vietnam; (ii) examine the performance of nonlinear mixed models using ecoregion as a random effect and develop separate equations for both mixed-species and dominant genera; and (iii) validate the accuracy of selected equations and compare them with both allometric equations for dipterocarps in Indonesia [13] and general pantropic equations [12,16,17,19].

2. Materials and Methods

2.1. Study Sites

The data for this study were collected in two of the eight agro-ecological zones or ecoregions of Vietnam [5,8]: CH and SE. (Binh Thuan Province). The ecoregion classification accounts for environmental variability with respect to climate, soil, and altitude, and is therefore relevant to forest stratification. Most of the country’s dipterocarp forests exist within the CH and SE ecoregions. Figure 1 shows the location of sample plots of dipterocarps of two agro-ecological zones.
The CH has an area of approximately 600,000 ha of the dipterocarp forest located at an elevation of 197 to 417 m with less than 10° slope. Mean annual rainfall is 1600 mm with a dry season lasting four months and mean annual temperature of 25.3 °C [23]. The SE ecoregion has an area of approximately 50,000 ha of the dipterocarp forest located at an elevation of 230 m with less than 3° slope. Mean annual rainfall for this ecoregion is 1003 mm with a dry season lasting three months and mean annual temperature of 25.5 °C [23]. Both ecoregions have a mixed-dipterocarp forest that is distributed primarily on a soil type of igneous rocks [24] with stand density ranging from 256 to 1292 trees ha−1 and basal area ranging from 10 to 25 m2 ha−1.

2.2. Sampling Design and Data Collection

Forest inventory data for the CH and SE ecoregions was collected with the support of Vietnam UN-REDD+ Phase I Program. Most of the data for the CH ecoregion was collected by the main author of this study (Huy et al. [25]) with the support of the Vietnam Ministry of Education and Training. Thirteen 0.25 ha (50 m × 50 m) plots from the CH, home of the dipterocarp forests in Vietnam, and one 1 ha plot from the SE ecoregion were surveyed. Attributes measured were plot location, forest status, basal area, slope, soil type, and standing tree measurements: species, DBH (cm), and tree height (H, m) of all trees with DBH ≥ 5 cm. Within the sample plots, the selection of sample trees was proportional to the number of trees in each of the 10 cm DBH class. The DBH of sample trees ranged from 5 to 50 cm and heights ranged from 2 to 24 m. A complete list of the species, families, and genera of the sample trees used in modeling is provided in Table 1. A total of 222 trees were destructively sampled for data collection between the CH (168 trees, 10–15 trees per 0.25-ha plot) and SE (54 trees) ecoregion.
Species name was recorded and DBH, H were measured on sample trees prior to felling. Total tree height (H) and green component mass for stems with bark, branches, and foliage was measured post felling. To establish the green-to-dry ratio of each tree, five samples for stem, three for branches, and two for old and new leaves were brought to the laboratory in sealed plastic bags to determine moisture content and reweighed to determine their mass in grams. Samples were then dried at 105 °C until a constant weight was reached [13,21,26]. The dry biomass of stem, branches, and leaves was estimated by multiplying the green biomass with the averaged green-to-dry ratio of each component. The total dry weight of a tree (AGB) was obtained by summing the dry biomass of the stem, branches, and foliage.
For wood density (WD) analysis, samples were taken from every one-fifth of the stem length. The water displacement method was used to measure the volume, and the averaged WD (g/cm3) per sample tree was calculated as oven-dry weight divided by volume at saturation.
Table 2 shows a summary for each of the predictors and the response variables of the destructively sampled trees for mixed-species and two main genera: Dipterocarpus and Shorea.

2.3. Equation Development and Selection

The dataset was randomly split 200 times into model development dataset (70% or 156 trees) and validation dataset (30% or 66 trees). After selecting, validating, and comparing the equations, the final parameters of the selected equations were obtained using the entire dataset (222 trees).
Equations were developed at two levels: mixed-species (all trees) and dominant genera of the dipterocarps (Dipterocarpus and Shorea). Exploratory analysis of destructively sampled tree data by ecoregion, family, and genus suggested the use of a power model with random effects (Figure 2).
The nonlinear mixed-effects models used in this study had the following form:
Y i j = ( a o + a i ) × X i j ( b o +   b i ) + ε i j
ε i j   ~   i i d   N ( 0 , σ 2 )
where Y i j is the observed ABG (kg) for the jth tree associated with the ith ecoregion; a o and b o are parameters of the model; a i and b i are parameters associated with the random effect of the ith ecoregion; and X i j is the DBH (cm), H (m), WD (g/cm3), DBH2H (m3), or DBH2HWD (kg) for the jth tree in ith ecoregion; and ε i j is the random error associated with the jth tree in ith ecoregion. The input variables DBH2H and DBH2HWD for X i j are approximations of volume and AGB, respectively, and were calculated as follows:
DBH 2 H = ( DBH 100 ) 2 × H
DBH 2 HWD = DBH 2 H × WD × 1000
Preliminary analysis also indicated that heterogeneous variance of residuals existed in all AGB models, as the variance tended to increase with increasing diameters. A power variance function as defined in Equation (5) was used to take into account the heterogeneity of residuals.
Var ( ε i j ) = σ 2 ( X i j ) k
where σ2 is the residual sum of squares, X i j is the weighting variable (DBH, DBH2H, or DBH2HWD in this study), and k is the variance function coefficient. The weighted nonlinear mixed-effects models with a specified variance function were fit using the NLME package [27] in R statistical software [28].
Allometric equations were first developed for mixed-species (all trees) and the best model forms for each input variable were selected. Best models were determined as those with significant parameters, no obvious issues with diagnostic plot, lower Akaike information criterion (AIC), and higher adjusted R2. The random effect of ecoregion was then tested for significance. Finally, the equations were built for two dominant genera of dipterocarp forests, Dipterocarpus and Shorea.

2.4. Equation Validation and Comparison

Performance of selected equations were validated and compared to each other and to equations developed for the dipterocarp and pantropic forests in terms of percent bias, root mean square percentage error (RMSPE), and mean absolute percent error (MAPE) [29]. These indicators were calculated for the validation dataset and smaller values were preferred. The cross-validation statistics were computed for each realization of randomly selected data and averaged over the 200 realizations [30].
Percent Bias = 100 R r = 1 R i = 1 n r [ y r i y ^ r i y r i ] / n r
RMSPE = 100 R r = 1 R i = 1 n r ( y r i   y ^ r i y r i ) 2 / n r
MAPE =   100 R r = 1 R i = 1 n r [ | y r i   y ^ r i | y r i ] / n r
where, R is the number of realizations (200);   n r is the number of trees per realization r; y r i and y ^ r i are the observed and predicted AGB (kg) for the ith tree in realization r, respectively.
The validation dataset was used to evaluate the performance of the following allometric equations developed for pantropic and dipterocarp forests of Indonesia and compared with the selected equations for mixed-species and genera developed in this study:
Pan-tropic equations for mixed-species:
Brown (1997) [12]:
AGB = exp ( 2.134   +   2.530 × log ( DBH ) )
IPCC (2003) [19]:
AGB = exp ( 2.289   +   2.649 × log ( DBH ) 0.021 × ( log ( DBH ) ) 2 )
Chave et al. (2005) [16] for dry forests:
AGB = WD × exp ( 0.667   + 1.784 × log ( DBH ) +   0.207 × ( log ( DBH ) ) 2 0.0281 × ( log ( DBH ) ) 3 )
Chave et al. (2014) [17]:
AGB =   0.0673 × ( WD × DBH 2 × H ) 0.976
Dipterocarp forest equations of Basuki et al. [13] for Indonesia:
Basuki et al. I (2009) for mixed-species:
AGB = exp ( 1.201 + 2.196 × log ( DBH ) )
Basuki et al. II (2009) for mixed-species:
AGB = exp ( 0.744 + 2.188 × log ( DBH ) + 0.832 × log ( WD ) )  
Basuki et al. I (2009) for Dipterocarpus genus:
AGB = exp ( 1.232 + 2.178 × log ( DBH ) )
Basuki et al. I (2009) for Shorea genus:
AGB = exp ( 2.193   +   2.371 × log ( DBH ) )
where, log(·) is natural logarithm and all other variables are as previously defined.
Finally, after selecting, validating, and comparing the equations with the pantropic and dipterocarp equations in Indonesia, the parameters of the selected equations in this study were obtained using the entire dataset (222 trees).

3. Results

Sizable differences were observed in WD among tree species, genus, and family. Average and standard deviation of WD for each plant species, genus, and family are presented in Table 1 and Figure 3.

3.1. Equations for Mixed-Species

The summary of the model fit for fixed- and mixed-effects models using different covariates to describe AGB in mixed-species forests are given in Table 3. The adjusted R2 of all models was more than 0.9 and the models that used separate input variables produced smaller AIC and RMSPE values compared to the models that used combined variables as the covariates. Therefore, equations with separate input variables were selected for mixed-species AGB estimation. Parameter estimates for these equations are presented in Table 4 and all parameters were statistically significant at the 0.001 level. Figure 4 shows the plots of fitted values and weighted residuals of those equations. The AIC value of the three-variable (DBH, H, and WD) model with random effect of ecoregion was smaller compared to the same model without random effect. However, there was no substantial changes in the parameter estimates for different ecoregions. Therefore, we selected the fixed-effects model as our final model.
The effect of ecoregion was then tested on the selected mixed-species models. The inclusion of ecoregion as a random effect did not improve the model fit (Table 3), nor did it change the parameter estimates substantially. After examining the between-ecoregion variability of the coefficient of the power equations, we found that there were no significant differences between the two ecoregions. Thus, the selected models with different input variables for mixed-species (equations in Table 4) can be applied to either CH or SE without considering the effect of different ecoregions. Among these models, the model with three input variables (DBH, H, and WD) performed the best and had the following form:
AGB = 0.06203 × DBH 2.26430 × H 0.51415 × WD 0.79456  

3.2. Equations for Dominant Genera

The summary of the model fit for fixed-effects models using different covariates to describe AGB for these genera are given in Table 5. AIC and R2 were very close between models for the same genus. Therefore, models with the fewest parameters were selected for parsimony and their parameter estimates are presented in Table 6.
However, as increasing the number of input variables from one (DBH) to three (DBH, H, and WD) did not improve model fit for either genus, the following simple AGB models with DBH as the only independent variable were selected for Dipterocarpus and Shorea, respectively.
Dipterocarpus :   AGB =   0.03713 × DBH 2.73813
Shorea :   AGB =   0.07483 × DBH 2.54496
When differences between observed and predicted AGB values were compared over tree size, no evidence of lack of fit was found. The models predicted AGB across tree sizes well for both Dipterocarpus and Shorea. The plots of fitted and weighted residuals of the selected models showed no trend or patterns for both genera (Figure 5).

3.3. Equation Validation and Comparison

Mixed-species equations:
Cross-validation was employed for the four selected mixed-species equations (M1–M4) and compared to the mixed-species equations for dipterocarp of Indonesia (Equations (13) and (14)) and pantropic forests (Equations (9)–(12)). Averaged validation statistics over the 200 realizations for mixed-species equations are presented in Table 7 for all three DBH classes.
The MAPE of the four selected equations (M1–M4) ranged from 24.0% to 25.7%, while the range was from 36.6% to 81.9% for the pantropic equations. Compared to the generic pantropic equations, the equations regionally specific to Vietnam reduced MAPE by 13% to 56% for dipterocarp forests in Vietnam.
The AGB estimates obtained from the four selected mixed-species equations (M1–M4), closely match the observed values of the validation dataset (Figure 6) and produced the lowest values for percent bias, RMSPE, and MAPE in all three DBH classes (Table 7). The pantropic equations and the equations developed for dipterocarp forests of Indonesia (Equations (9)–(14)) all overestimated AGB for the validation dataset (Figure 6) and had significantly higher percent bias, RMSPE, and MAPE than produced by equations M1–M4 in all three DBH classes (Table 7).
Genus-level equations:
The selected genus-level models for Dipterocarpus and Shorea (D1 and S1, respectively) were cross-validated and compared to those of Basuki et al. [13] that were developed for the same genera, though located in Indonesia (Equations (15) and (16), respectively). Predictions from Equations (D1), (S1) and (16) closely matched the random validation data, with low values for percent bias, RMSPE, and MAPE in all three DBH classes (Table 8, Figure 7). The validation statistics (percent bias, RMSPE, and MAPE) were higher for Equation 15, Basuki et al.’s Dipterocarpus equation (overall bias −50.9% vs. 7.9%), but the biases were similar for Shorea (10.4% vs. 9.6%).

4. Discussion

We obtained a large dataset consisting of 222 destructively sampled trees from two ecoregions, the original location of dipterocarp forests, to develop biomass equations for the tropical dipterocarp forests of Vietnam. Since the selection of sample trees was based on the diameter distribution, the number of trees sampled in the larger diameter class is still small. Thus, having additional sample trees in the larger diameter classes would further improve the reliability of these equations. Additionally, biomass allometry is affected by stand structure, developmental stage and other stand attributes. Thus, having such attributes in the biomass equations might further improve the predictive ability of such equations. However, this could complicate the model forms.

4.1. Input Variables for AGB Mixed-Species Equations for Dipterocarp Forests

Aboveground tree biomass is known to be affected by tree diameter, height, and wood density; diameter and height yield volume with wood density act as a scaling factor for biomass. The mixed-species equations available in this region have used different explanatory variables to estimate aboveground biomass e.g., [12,13,16,17,19]. For modeling AGB in mixed-species dipterocarp forests, we found that an equation with DBH as the only covariate performed poorly (highest AIC) when compared to equations that also incorporated H, WD, or both H and WD variables. Furthermore, the mixed-species equation developed in this study that incorporated all three variables resulted in the lowest uncertainty and lowest AIC. This could be due to the fact that at the mixed-species level, the wood density component helps in adjusting for the differences in allometric relationship between DBH and H among species. Thus, the best option for mixed-species was the AGB equation with three variables, DBH, H, and WD, in the form of M4. We would like to point out the use of height in fitting biomass equations. Commonly, the height measured after felling the sample trees is used in model development but the standing tree heights are used in applying these equations. This discrepancy in the variable used in model fitting and application can result in biased estimate of the AGB and the magnitude of bias depends on the sign and magnitude of bias in height measurement as well as the model form and coefficients.

4.2. Equations at the Mixed-Species Level in the Dipterocarp Forests of Southeast Asia

While generic equations for estimating AGB in tropical forests exist [12,16,17,19], the dipterocarp forest is a specific forest type dominated by the Dipterocarpaceae plant family and has specific site conditions that differ from those of other forest types. Therefore, the application of pantropic equations for dipterocarp forests must be evaluated.
The mixed-species pantropic equations of Brown [12], IPCC [19], and Chave et al. [16,17] for tropical forests all overestimate the AGB of mixed-species dipterocarp forests in Vietnam (Figure 6). Based on the Vietnam mixed-species validation dataset, the pantropic equations produced higher percent bias, RMSPE, and MAPE than the equations developed in this study for mixed-species dipterocarp forests of Vietnam (Table 7).
Basuki et al.’s [13] equations for mixed-species dipterocarp forests in Indonesia (Equations (13) and (14)) also overestimated the AGB of mixed-species dipterocarp forests in Vietnam (Figure 6). These mixed-species equations developed based on dipterocarp forests of Indonesia resulted in substantially higher MAPE (94% to 108%) than the selected mixed-species equations developed in this study (24% to 26%). The substantial differences in performance between these sets of equations can partially be explained by the different species used in the studies. There was no overlap in species used in equation development between the 26 species reported by Basuki et al. [14] and the 26 species from this study, despite the fact that they include some of the same main families and genera, such as Dipterocarpus and Shorea. Our evaluation shows that while equations developed for specific ecoregions (e.g., CH or SE) do not increase the reliability and accuracy of the AGB estimates when applied at the mixed-species level of the dipterocarps, regionally specific equations (e.g., Vietnam) offer considerable improvement over pantropic equations or equations developed for other regions (e.g., Indonesia). This finding is supported by the work of Basuki et al. [14], Nelson et al. (1999) [31], and Cairns et al. (2003) [32], which indicated that site-specific equations improve the accuracy of biomass estimates [33].

4.3. Equations at the Genus Level in the Dipterocarp Forests of Southeast Asia

The results demonstrate that the allometric equations developed at the genus-level in different regions (e.g., Vietnam vs. Indonesia) were not significantly different in accuracy (Figure 7 and Table 8). As a subset of the species list for mixed-species equations, the species list for the Dipterocarpus and Shorea genera used in genus-level equation development by Basuki et al. [13] differed from that of this study (Table 1). However, the prediction of AGB values from the genus equations of Basuki et al. [13] (Equations 16 for Shorea genus) and from this study (S1) were not substantially different. This could be because the species in the same genus have similar of tree foliage structure, stem form, and in particular the WD, compared to trees in a different genus (Figure 3). Such differences are not accounted for in the generic equations. These evaluations support the application of common equations for each plant genus of the dipterocarp forests throughout Southeast Asia. Future research using additional data would further validate our findings.

5. Conclusions

All variables of diameter at breast height (DBH), height (H), and wood density (WD) affected aboveground biomass (AGB). Increasing number of input variables from one to three reduced the uncertainty of the estimates. However, from an application perspective, a simpler model might be preferred as obtaining the wood density variable is very difficult. In general, for tropical forests with a large number of tree species, using the mixed-species model with multi-tree input variables was appropriate to improve the accuracy of the AGB estimates.
The predicted values from generic equations overestimated AGB in mixed-species dipterocarp forests. Because the dipterocarp forest is a specific forest type, the equations specific to different regions of Southeast Asia are appropriate and increase the reliability and accuracy of the AGB estimates.
The power models with DBH as the only input variable for the Dipterocarpus and Shorea genera were developed and validated for application to the dipterocarps. We speculate that using genus-specific equations improve the reliability of the AGB estimates in the dipterocarp forests over generic mixed-species equations. There were no substantial differences in accuracy among the developed equations for the same genus in different regions or countries. Therefore, it seems reasonable to develop and apply the common equations for each genus of the dipterocarps throughout Southeast Asia.

Acknowledgments

This work was undertaken in 2014, with the support of the UN-REDD+ Vietnam Phase II Program and built on an extensive field measurement campaign supported by the UN-REDD+ Vietnam Phase I Program in 2012 and 2013, and a research was carried out in 2010–2013 in the Central Highlands by Tay Nguyen University with funding from Ministry of Education and Training. Four institutions collaborated on the field work and analysis: Forest Inventory and Planning Institute, Vietnamese Academy of Forest Sciences, Vietnam National University of Forestry and Tay Nguyen University, with technical assistance from FAO. The authors also would like to acknowledge technical support provided by Gael Sola (FAO). We acknowledge Caryn Davis for her edits on an earlier manuscript.

Author Contributions

Bao Huy, Nguyen Dinh Hung, Phung Van Khoa, Vu Tan Phuong collected data set and analyzed the raw data; Bao Huy, Hailemariam Temesgen, Krishna P. Poudel developed methodology and analyzed data for equation development and validation; Bao Huy, Krishna P. Poudel, Karin Kralicek, Hailemariam Temesgen, Nguyen Dinh Hung, Phung Van Khoa and Vu Tan Phuong wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Measurement, Reporting & Verification (MRV) Framework Document. Available online: http://theredddesk.org/resources/measurement-reporting-and-verification-framework-document (accessed on 17 August 2016).
  2. UNFCCC (United Nations Framework Convention on Climate Change). Report of the Conference of the Parties on its Thirteenth Session; Proceedings of UNFCCC: Bonn, Germany, 2008. [Google Scholar]
  3. Appanah, S. Management of Natural Forests. In A Review of Dipterocarp: Taxonomy, Ecology and Silviculture; Appanah, S., Turnbull, J.M., Eds.; Center for International Forestry Research (CIFOR): Bogor, Indonesia, 1998; pp. 130–149. [Google Scholar]
  4. Maury-Lechon, G.; Curtet, L. Biogeography and Evolutionary Systematics of Dipterocarpaceae. In A Review of Dipterocarp: Taxonomy, Ecology and Silviculture; Appanah, S., Turnbull, J.M., Eds.; Center for International Forestry Research (CIFOR): Bogor, Indonesia, 1998; pp. 5–44. [Google Scholar]
  5. Huy, B. Part B1: Equations for Biomass of Aboveground Trees, Branches and Leaves Biomass in Evergreen Broadleaved Forests, and for Aboveground Biomass of Six Tree Families in Evergreen and Deciduous Forests. In Allometric Equations at National Scale for Estimating Tree and Forest Biomass in Vietnam; Sola, G., Inoguchi, A., Garcia-Perez, J., Donegan, E., Birigazzi, L., Henry, M., Eds.; UN-REDD Programme: Ha Noi, Vietnam, 2014. [Google Scholar]
  6. Jenkins, J.C.; Chojnacky, D.C.; Heath, L.S.; Birdsey, R.A. Comprehensive Database of Diameter-based Biomass Regressions for North American Tree Species; General Technical Report NE-319; USDA Forest Service: Newtown Square, PA, USA, 2004; p. 45. [Google Scholar]
  7. Jenkins, J.C.; Chojnacky, D.C.; Heath, L.S.; Birdsey, R.A. National-scale biomass estimators for USA tree species. J. For. Sci. 2003, 49, 12–35. [Google Scholar]
  8. Sola, G.; Inoguchi, A.; Garcia-Perez, J.; Donegan, E.; Birigazzi, L.; Henry, M. Allometric Equations at National Scale for Tree Biomass Assessment in Vietnam, Context, methodology and summary of the results; UN-REDD Programme: Ha Noi, Vietnam, 2014. [Google Scholar]
  9. Sola, G.; Inoguchi, A.; Phuong, V.T.; Huy, B.; Khoa, P.V.; Hung, N.D.; Xuan, N.V.; Henry, M. Allometric Equations at National Scale for Tree Biomass Assessment in Vietnam; UN-REDD Programme: Ha Noi, Vietnam, 2014. [Google Scholar]
  10. Ter-Mikaelian, M.T.; Korzukhin, M.D. Biomass equations for sixty-five North American tree species. J. For. Ecol. Manag. 1997, 97, 1–24. [Google Scholar] [CrossRef]
  11. Pearson, T.R.H.; Brown, S.; Birdsey, R.A. Measurement Guidelines for the Sequestration of Forest Carbon; General Technical Report NRS-18; USDA Forest Service: Newton Square, PA, USA, 2007. [Google Scholar]
  12. Estimating Biomass and Biomass Change of Tropical Forests: A Primer. (FAO Forestry Paper-134). Available online: http://www.fao.org/docrep/W4095E/w4095e00.htm#Contents (accessed on 16 August 2016).
  13. Basuki, T.M.; Van Lake, P.E.; Skidmore, A.K.; Hussin, Y.A. Allometric equations for estimating the above-ground biomass in the tropical lowland Dipterocarp forests. For. Ecol. Manag. 2009, 257, 1684–1694. [Google Scholar] [CrossRef]
  14. Picard, N.; Rutishauser, E.; Ploton, P.; Ngomanda, A.; Henry, M. Should tree biomass allometry be restricted to power models? For. Ecol. Manag. 2015, 353, 156–163. [Google Scholar] [CrossRef]
  15. Chambers, J.Q.; dos Santos, J.; Ribeiro, R.J.; Higuchi, N. Tree damage allometric relationship and above-ground net primary production in central Amazon forest. For. Ecol. Manag. 2001, 152, 73–84. [Google Scholar] [CrossRef]
  16. Chave, J.; Andalo, C.; Brown, S.; Cairns, M.A.; Chambers, J.Q.; Eamus, D.; Folster, H.; Fromard, F.; Higuchi, N.; Kira, T.; et al. Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 2005, 145, 87–99. [Google Scholar] [CrossRef] [PubMed]
  17. Chave, J.; Mechain, M.R.; Burquez, A.; Chidumayo, E.; Colgan, M.S.; Delitti, W.B.C.; Duque, A.; Eid, T.; Fearnside, P.M.; Goodman, R.C.; et al. Improved allometric models to estimate the aboveground biomass of tropical trees. Glob. Change Biol. 2014, 20, 3177–3190. [Google Scholar] [CrossRef] [PubMed]
  18. IPCC National Greenhouse Gas Inventories Programme. Good Practice Guidance for Land Use, Land-Use Change and Forestry; IPCC: Hayama, Japan, 2003; p. 295. [Google Scholar]
  19. National Greenhouse Gas Inventories Programme. IPCC Guidelines for National Greenhouse Gas Inventories; IGES: Hayama, Kanagawa, Japan, 2006. [Google Scholar]
  20. Keller, M.; Palace, M.; Hurtt, G. Biomass estimation in the Tapajos National Forest. Brazil examination of sampling and allometric uncertainties. For. Ecol. Manag. 2001, 154, 371–382. [Google Scholar] [CrossRef]
  21. Ketterings, Q.M.; Coe, R.; van Noordwijk, M.; Ambagau, Y.; Palm, C.A. Reducing uncertainty in the use of allometric biomass equations for predicting above ground tree biomass in mixed secondary forests. For. Ecol. Manag. 2001, 146, 199–209. [Google Scholar] [CrossRef]
  22. Huy, B.; Poudel, K.P.; Temegen, H. Aboveground biomass equations for evergreen broadleaf forests in South Central Coastal of Vietnam: Selection of eco-regional or pantropical models. For. Ecol. Manag. 2016, 376, 276–283. [Google Scholar] [CrossRef]
  23. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  24. Nachtergaele, F.O.; van Velthuizen, H.T.; Verelst, L.; Wiberg, D.; Batjes, N.H.; Dijkshoorn, J.A.; van Engelen, V.W.P.; Fischer, G.; Jones, A.; Montanarella, L.; et al. Harmonized World Soil Database (Version 1.2); FAO: Rome, Italy; IIASA: Laxenburg, Austria, 2012; p. 42. [Google Scholar]
  25. Huy, B.; Ly, C.T.; Hung, V.; Nguyen, T.T.H.; Dinh, N.D. Define CO2 Absorbed in the Broad-Leaf Ever-Green Forests in the Central Highlands of Vietnam to Join the REDD+ Program; Technical Report; Ministry of Education and Training: Ha Noi, Vietnam, 2013.
  26. Stewart, J.L.; Dunsdon, A.J.; Hellin, J.J.; Hughes, C.E. Wood Biomass Estimation of the Central American Dry Zone Species; Tropical Forestry Papers 26; Oxford Forestry Institute, Department of Plant Sciences, University of Oxford: Oxford, UK, 1992. [Google Scholar]
  27. Nlme: Linear and nonlinear mixed effects models. Available online: http://CRAN.R-project.org/package=nlme. (accessed on 19 August 2016).
  28. The R Project for Statistical Computing. Available online: https://www.R-project.org/ (accessed on 16 August 2016).
  29. Swanson, D.A.; Tayman, J.; Bryan, T.M. MAPE-R: A rescaled measure of accuracy for cross-sectional subnational population forecasts. J. Pop. Res. 2011, 28, 225–243. [Google Scholar] [CrossRef]
  30. Temesgen, H.; Zhang, C.H.; Zhao, X.H. Modelling tree height-diameter relationships in multi-species and multi-layered forests: A large observational study from Northeast China. For. Ecol. Manag. 2014, 316, 78–89. [Google Scholar] [CrossRef]
  31. Nelson, B.W.; Mesquita, R.; Pereira, J.L.G.; de Souza, S.G.A.; Batista, G.T.; Couta, L.B. Allometric regressions for improved estimate of secondary forest biomass in the Central Amazon. For. Ecol. Manag. 1999, 177, 149–167. [Google Scholar] [CrossRef]
  32. Cairns, M.A.; Olmsted, I.; Granados, J.; Argaez, J. Composition and aboveground tree biomass of the dry semi-evergreen forest on Mexico’s Yucatan Peninsula. For. Ecol. Manag. 2003, 186, 125–132. [Google Scholar] [CrossRef]
  33. Temesgen, H.; Affleck, D.; Poudel, K.P.; Gray, A.; Sessions, J. A review of the challenges and opportunities in estimating above ground forest biomass using tree-level models. Scand. J. For. Res. 2015, 30, 326–335. [Google Scholar]
Figure 1. The location of sample plots in two ecoregions where the dipterocarp forests of Vietnam are distributed.
Figure 1. The location of sample plots in two ecoregions where the dipterocarp forests of Vietnam are distributed.
Forests 07 00180 g001
Figure 2. AGB vs. DBH of destructively sampled trees by ecoregion, plant family, and genus. The trend of relationship between AGB and DBH follows power model.
Figure 2. AGB vs. DBH of destructively sampled trees by ecoregion, plant family, and genus. The trend of relationship between AGB and DBH follows power model.
Forests 07 00180 g002
Figure 3. Box plot of WD by sampled genera in alphabetical order. Numbers inside the boxplots represent the sample size in that genus.
Figure 3. Box plot of WD by sampled genera in alphabetical order. Numbers inside the boxplots represent the sample size in that genus.
Forests 07 00180 g003
Figure 4. Plots of selected AGB models for mixed-species: fitted vs. observed (left column) and residuals vs. fitted (right column) using the entire dataset of 222 destructively sampled trees. (a) Model AGB = a × DBHb, (b) Model AGB = a × DBHb × Hc, (c) Model AGB = a × DBHb × Hc × WDd.
Figure 4. Plots of selected AGB models for mixed-species: fitted vs. observed (left column) and residuals vs. fitted (right column) using the entire dataset of 222 destructively sampled trees. (a) Model AGB = a × DBHb, (b) Model AGB = a × DBHb × Hc, (c) Model AGB = a × DBHb × Hc × WDd.
Forests 07 00180 g004
Figure 5. Plots of fitted values vs. DBH (left) and fitted values vs. weighted residuals (right) of the selected model AGB = a × DBHb for two main genera: (a) Dipterocarpus (the entire dataset, n = 94 trees), and (b) for Shorea (the entire dataset, n = 36 trees).
Figure 5. Plots of fitted values vs. DBH (left) and fitted values vs. weighted residuals (right) of the selected model AGB = a × DBHb for two main genera: (a) Dipterocarpus (the entire dataset, n = 94 trees), and (b) for Shorea (the entire dataset, n = 36 trees).
Forests 07 00180 g005aForests 07 00180 g005b
Figure 6. Observed vs. predicted AGB obtained from the mixed-species equations developed in this study, pantropical equations, and equations for dipterocarp forest of Indonesia.
Figure 6. Observed vs. predicted AGB obtained from the mixed-species equations developed in this study, pantropical equations, and equations for dipterocarp forest of Indonesia.
Forests 07 00180 g006
Figure 7. Predictions of AGB by the DBH-based genus-specific models developed in this study and those developed for dipterocarp forest of Indonesia for (a) Dipterocarpus and (b) Shorea genera.
Figure 7. Predictions of AGB by the DBH-based genus-specific models developed in this study and those developed for dipterocarp forest of Indonesia for (a) Dipterocarpus and (b) Shorea genera.
Forests 07 00180 g007
Table 1. Number of destructively sampled trees (N), average WD and standard deviation (SD) of WD by family, genus, and species. Family and genus-level averages are in bold.
Table 1. Number of destructively sampled trees (N), average WD and standard deviation (SD) of WD by family, genus, and species. Family and genus-level averages are in bold.
Family/Genus/SpeciesNAverage WD (g/cm3)SD of WD (g/cm3)
Chrysobalanaceae30.5356650.016343
Parinari30.5356650.016343
Parinari anamensis Hance30.5356650.016343
Combretaceae230.7470570.094749
Terminalia230.7470570.094749
Terminalia alata Roth160.7300930.100805
Terminalia corticosa Pierre ex Laness70.7858300.070653
Dilleniaceae20.7004230.246027
Dillenia20.7004230.246027
Dillenia sp.20.7004230.246027
Dipterocarpaceae1300.6575990.098560
Dipterocarpus940.6272940.080587
Dipterocarpus intricatus Dyer210.5896670.037682
Dipterocarpus obtusifolius Teijsm. ex Miq.100.6407040.050715
Dipterocarpus tuberculatus Roxb.630.6377070.091056
Shorea360.7367290.098185
Shorea obtusa Wall. ex Blume280.7724000.074928
Shorea siamensis Miq.80.6118810.060359
Phyllanthaceae70.6124780.074246
Aporosa70.6124780.074246
Aporosa sp.70.6124780.074246
Hypericaceae10.625000-
Cratoxylum10.625000-
Cratoxylum formosum (Jacq.) Benth. & Hook.f. ex Dyer10.625000-
Irvingiaceae60.7048330.017702
Irvingia60.7048330.017702
Irvingia malayana Oliv. ex A.W.Benn.60.7048330.017702
Lauraceae20.5855000.024749
Cryptocarya20.5855000.024749
Cryptocarya sp.20.5855000.024749
Lecythidaceae50.5925860.134872
Careya50.5925860.134872
Careya arborea Roxb.50.5925860.134872
Leguminosae190.6836480.082150
Dalbergia40.6060630.049240
Dalbergia cochinchinensis Pierre30.5876670.040079
Dalbergia nigrescens var. saigonensis (Pierre) Gagnep.10.661252-
Sindora10.818989-
Sindora siamensis Miq.10.818989-
Xylia140.6961480.073180
Xylia xylocarpa (Roxb.) Taub.140.6961480.073180
Loganiaceae20.6006130.006983
Strychnos20.6006130.006983
Strychnos nux-blanda A.W. Hill20.6006130.006983
Myrtaceae60.5212440.062018
Syzygium60.5212440.062018
Syzygium cumini (L.) Skeels20.5897320.062332
Syzygium sp.40.4870000.020607
Rubiaceae150.5890790.051010
Adina20.5710000.000000
Adina pilulifera (Lam.) Franch. ex Drake20.5710000.000000
Haldina20.6112930.031190
Haldina cordifolia (Roxb.) Ridsdale20.6112930.031190
Morinda10.463000-
Morinda citrifolia L.10.463000-
Nauclea100.6008600.042858
Nauclea orientalis (L.) L.100.6008600.042858
Sapindaceae10.630000-
Sapindus10.630000-
Sapindus saponaria L.10.630000-
Grand Total2220.6564750.101823
Table 2. Summary for each of the predictors and the response variables of the destructively sampled trees for mixed-species and two main genera.
Table 2. Summary for each of the predictors and the response variables of the destructively sampled trees for mixed-species and two main genera.
LevelSummaryDBH (cm)H (m)WD (g/cm3)AGB (kg)
Mixed-SpeciesMin3.42.50.3791.3
Average17.610.30.656184.9
Max48.823.50.9531710.8
Standard Deviation10.34.50.102315.0
n222222222222
Dipterocarpus GenusMin4.93.80.3791.5
Average20.411.20.627285.6
Max48.823.50.8581710.8
Standard Deviation12.45.30.081431.8
n94949494
Shorea GenusMin5.64.40.5072.9
Average15.29.40.73786.4
Max23.014.10.917205.9
Standard Deviation4.22.40.09852.5
n36363636
Table 3. Comparison of different AGB models for mixed-species with and without the random effect of ecoregion (156 trees).
Table 3. Comparison of different AGB models for mixed-species with and without the random effect of ecoregion (156 trees).
Model SelectedModel FormRandom Effect Weight VariableAICAdj. R2RMSPE
Covariate: DBH
M1AGB = a × DBHbNone 1/DBH k13350.96659.0
AGB = a × DBHbEcoregion1/DBH k13350.966 59.0
M2Covariates: DBH and H
AGB = a × (DBH2H)bNone1/(DBH2H) k13740.96748.1
AGB = a × DBHb × HcNone1/DBH k13240.97551.1
AGB = a × DBHb × HcEcoregion1/DBH k13240.975 51.1
M3Covariates: DBH and WD
AGB = a × (DBH2WD)bNone1/(DBH2WD) k13370.90153.7
AGB = a × DBHb × WDcNone1/DBH k13220.93354.8
AGB = a × DBHb × WDcEcoregion1/DBH k13220.93354.8
M4Covariates: DBH, H, and WD
AGB = a × (DBH2HWD)bNone1/(DBH2HWD) k13230.94539.3
AGB = a × DBHb × Hc × WDdNone1/DBH k12980.94043.7
AGB = a × DBHb × Hc × WDdEcoregion1/DBH k12800.94043.7
k is the variance function coefficient.
Table 4. Selected AGB models with different input variables, their parameter estimates, and standard errors for mixed-species of the dipterocarp forests (based on the entire dataset, 222 trees).
Table 4. Selected AGB models with different input variables, their parameter estimates, and standard errors for mixed-species of the dipterocarp forests (based on the entire dataset, 222 trees).
Model SelectedCovariate (s)Model FormParameter EstimatesStandard Error
M1DBH AGB = a × DBH b a0.047420.00434
b2.666630.03151
M2DBH, H AGB =   a × DBH b × H c a0.038440.00377
b2.407560.05861
c0.404080.07988
M3DBH, WD AGB = a × DBH b × WD c a0.071260.00843
b2.609400.03290
c0.592750.59275
M4DBH, H, WD AGB = a × DBH b × H c × WD d a0.062030.00679
b2.264300.05754
c0.514150.07497
d0.794560.11420
Table 5. Comparison of different AGB models with different input variables for Dipterocarpus (67 trees) and for Shorea (26 trees).
Table 5. Comparison of different AGB models with different input variables for Dipterocarpus (67 trees) and for Shorea (26 trees).
Model SelectedCovariate(s)Model FormWeight VariableAICAdj. R2RMSPE
Genus: Dipterocarpus
D1DBHAGB = a × DBHb1/DBH k6000.97060.7
D2DBH, HAGB = a × (DBH2H)b1/DBH k6160.96546.5
AGB = a × DBHb × Hc1/DBH k6000.97954.1
D3DBH, WDAGB = a × (DBH2WD)b1/(DBH2WD) k6030.92061.8
AGB = a × DBHb × WDc(#)1/DBH k6050.95259.9
D4DBH, H, WDAGB = a × (DBH2HWD)b1/(DBH2HWD) k6020.97045.2
AGB = a × DBHb × Hc × WDd1/DBH k6010.97049.8
Genus: Shorea
S1DBHAGB = a × DBHb1/DBH k2200.91925.4
S2DBH, HAGB = a × (DBH2H)b1/(DBH2H) k2310.86133.6
AGB = a × DBHb × Hc(#)1/DBH k2270.92125.9
S3DBH, WDAGB = a × (DBH2WD)b1/(DBH2WD) k2240.88122.3
AGB = a × DBHb × WDc(#)1/DBH k2260.91222.6
S4DBH, H, WDAGB = a × (DBH2HWD)b1/(DBH2HWD) k2200.89424.0
AGB = a × DBHb × Hc(#) × WDd1/DBH k2310.92021.0
(#): p-value of the parameter > 0.05 and k is the variance function coefficient.
Table 6. Parameter estimates and their standard errors for the selected AGB models with different input variables for Dipterocarpus (the entire dataset, 94 trees) and for Shorea (the entire dataset, 36 trees).
Table 6. Parameter estimates and their standard errors for the selected AGB models with different input variables for Dipterocarpus (the entire dataset, 94 trees) and for Shorea (the entire dataset, 36 trees).
Model SelectedCovariate (s)Model FormParameter EstimatesStandard Error
abab
Genus: Dipterocarpus
D1DBH AGB =   a × DBH b 0.037132.738130.004880.04310
D2DBH, H AGB =   a × ( DBH 2 H ) b 290.3701.039139.303990.01900
D3DBH, WD AGB =   a × ( DBH 2 WD ) b 0.093871.315390.010470.02101
D4DBH, H, WD AGB =   a × ( DBH 2 HWD ) b 0.458121.006730.047290.01707
Genus: Shorea
S1DBH AGB =   a × DBH b 0.074832.544960.028970.13722
S2DBH, H AGB =   a × ( DBH 2 H ) b 325.2640.9522024.20810.04909
S3DBH, WD AGB =   a × ( DBH 2 WD ) b 0.191141.163900.061810.06166
S4DBH, H, WD AGB =   a × ( DBH 2 HWD ) b 0.698380.921730.138070.03772
Table 7. Percent bias, RMSPE, and MAPE for the selected mixed-species models, mixed-species models for dipterocarp forest of Indonesia, and pantropical models for mixed-species by covariate (s).
Table 7. Percent bias, RMSPE, and MAPE for the selected mixed-species models, mixed-species models for dipterocarp forest of Indonesia, and pantropical models for mixed-species by covariate (s).
Covariate (s)ModelDBH Class (cm)Percent BiasRMSPEMAPE
DBHBrown (1997) [12]
Equation (9)
<15−86.5106.786.5
15–30−74.080.974.0
>30−79.092.779.0
All−80.9 101.5 80.9
IPCC (2003) [18]
Equation (10)
<15−86.8106.486.8
15–30−76.283.176.2
> 30−79.293.079.2
All−81.9102.081.9
Basuki et al. I (2009) [13]
Equation (13)
<15−132.2163.5132.2
15–30−64.672.064.6
> 30−37.757.647.8
All−92.3126.2 94.2
M1<152.532.024.3
15–30−3.421.919.4
>30−17.243.738.1
All−3.1 34.6 25.2
DBH, HM2<150.131.624.2
15–30−2.119.916.9
>30−18.041.836.2
All−3.9 33.524.0
DBH, WDChave et al. (2005) [16]
Equation (11)
<15−97.2117.397.3
15–30−73.883.873.8
>30−46.569.756.0
All−79.2105.081.1
Basuki et al. II (2009) [13]
Equation (14)
<15−139.2164.2139.2
15–30−90.399.190.3
>30−56.476.461.5
All−107.6136.9108.7
M3< 152.828.421.1
15–30−8.327.423.6
>30−15.547.239.6
All−4.136.525.7
DBH, H, WDChave et al. (2014) [17]
Equation (12)
<15−27.144.533.0
15–30−30.941.833.9
>30−42.960.348.2
All-31.253.536.6
M4< 150.126.520.6
15–30−7.423.720.2
>30−16.646.138.5
All−5.234.724.2
Table 8. Percent bias, RMSPE, and MAPE for the selected genus-specific models and genus-specific equations for dipterocarp forest of Indonesia.
Table 8. Percent bias, RMSPE, and MAPE for the selected genus-specific models and genus-specific equations for dipterocarp forest of Indonesia.
GenusCovariate (s)EquationDBH Class (cm)Percent BiasRMSPEMAPE
DipterocarpusDBHBasuki et al. I (2009) [13]
Equation (15)
<15−91.794.991.7
15–30−61.861.561.4
>30−5.730.527.6
All−50.970.659.7
D1<1518.419.518.4
15–30−5.217.217.2
>300.228.225.9
All7.924.221.8
ShoreaDBHBasuki et al. (2009) [13]
Equation (16)
<1513.315.715.8
15–306.819.417.9
>30---
All10.418.016.7
S1<15 16.818.218.4
15–30−1.223.522.7
>30---
All9.621.219.9

Share and Cite

MDPI and ACS Style

Huy, B.; Poudel, K.P.; Kralicek, K.; Hung, N.D.; Khoa, P.V.; Phương, V.T.; Temesgen, H. Allometric Equations for Estimating Tree Aboveground Biomass in Tropical Dipterocarp Forests of Vietnam. Forests 2016, 7, 180. https://doi.org/10.3390/f7080180

AMA Style

Huy B, Poudel KP, Kralicek K, Hung ND, Khoa PV, Phương VT, Temesgen H. Allometric Equations for Estimating Tree Aboveground Biomass in Tropical Dipterocarp Forests of Vietnam. Forests. 2016; 7(8):180. https://doi.org/10.3390/f7080180

Chicago/Turabian Style

Huy, Bao, Krishna P. Poudel, Karin Kralicek, Nguyen Dinh Hung, Phung Van Khoa, Vu Tan Phương, and Hailemariam Temesgen. 2016. "Allometric Equations for Estimating Tree Aboveground Biomass in Tropical Dipterocarp Forests of Vietnam" Forests 7, no. 8: 180. https://doi.org/10.3390/f7080180

APA Style

Huy, B., Poudel, K. P., Kralicek, K., Hung, N. D., Khoa, P. V., Phương, V. T., & Temesgen, H. (2016). Allometric Equations for Estimating Tree Aboveground Biomass in Tropical Dipterocarp Forests of Vietnam. Forests, 7(8), 180. https://doi.org/10.3390/f7080180

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

Article Metrics

Back to TopTop