Abstract

The application of biomass models for quantifying forests’ above-ground biomass is essential for sustainable forest management. However, lack of knowledge in modelig biomass of individual tree growth hinders the sustainable management of Dry Afromontane forests. In this study, models to estimate above-ground biomass were developed for Rhus ruspolii, Ekebergia capensis, and Nuxia congesta. To develop the models, a total of 45 trees from different diameter classes were selected, felled, and divided into different biomass compartments. For the model’s development, diameter at breast height (DBH), total height (TH), diameter at stump height (DSH), and wood density (WD) were used as independent variables. Models’ performances were evaluated using RSE, adjusted coefficient of determination, and AIC. Also, model validations were done by using rRMSE, mean absolute deviation, bias, and coefficient of variation. Models 5 (Adj-R2 = 0.92), 6 (Adj-R2 = 0.97), and 8 (Adj-R22 = 0.82) were the best fitted models for Nuxia congesta, Ekebergia capensis, and Rhus ruspolii, respectively. The average wood densities of Ekebergia capensis, Nuxia congesta, and Rhus ruspolii were 0.59, 0.50, and 0.69, respectively. The variation between observed biomass and estimated biomass using new models was statistically not significant (). Thus, the biomass models developed here can be important tools to accurately estimate above-ground biomass in the Menagesha Suba forest and can be integrated into decision support tools.

1. Introduction

To understand the role of forest ecosystems in the global carbon cycle and to manage forest systems sustainably, quantifying biomass and carbon stocks in tropical forest species is very essential. Biomass estimation is specifically important to understand the role of forests as a carbon sink and for sustainable forest management [1]. The availability of tree biometric data like wood density (WD), total tree height (TH), diameter at breast height (DBH), volume, and biomass are critical for supporting forest management decision-making and reducing costs for ground-based forest resource assessment [2, 3].

Above-ground biomass (AGB), which is vital to measure forest stocking, is considered within the framework of sustainable forest management and carbon accounting purposes [4]. Also, it is a useful tool for comparing structural and functional attributes of the forest across a wide range of environmental conditions [5]. Furthermore, forest biomass and wood basic density (WBD) are crucial to effectively implementing climate change mitigation policies and reducing emissions from deforestation and forest degradation (REDD+) under the United Nations Framework Convention on Climate Change (UNFCCC) [6]. Hence, wbd of standing trees has become a fundamental component of biomass determinations in ecosystem studies.

Destructive sampling is the most accurate and direct method of estimating tree biomass [7]. Nonetheless, it is limited to small areas or small tree sample sizes and not suggested for endangered tree species [8]. In such cases, using generic models becomes preferable. But generic models may underestimate or overestimate the biomass in the tropical forest due to the high diversity of the species. Using generic models at the site level may lead to the production of systematic errors [911]. This error could have occurred due to environmental variation among different forests [9, 11, 12]. Besides, the variation in wood anatomy in the tropical forest is also another source of error. To minimize the error produced due to environmental variability and variation in wood anatomy, site-species-specific biomass becomes imperative. However, biomass estimation models specific to tropical Africa remain limited [1315]. Even though some efforts have been made to develop biomass equations in tropical Africa [11, 16, 17], attempts to develop biomass models for sub-Saharan tropical species still remain scarce [18]. These led to biomass estimates in Africa relying on generic models [1921].

Ethiopia is one of the sub-Saharan countries endowed with diverse vegetation types. The Ethiopian highlands constitute large parts of the Afromontane regions in Africa [22, 23]. These regions are dominated by dry afromontane forests, which are a complex ecosystem characterized by a variable rainfall pattern and a prolonged dry season [24]. It is the second-richest vegetation type, with about 460 species, next to Acacia-Commiphora woodland [25]. However, deforestation and forest degradation in the dry Afromontane forest ecosystem due to landscape fragmentation and land-use pressure are serious challenges in the forest system [22, 26]. Thus, Ethiopia has proposed and implemented a Climate Resilient Green Economy (CRGE) strategy to reduce GHG emissions [27]. Accordingly, to achieve the CRGE strategies, numerous forest resource assessment campaigns have been undertaken to understand the AGB status and management options of the forest resources [28]. Measurement, reporting, and verification (MRV) is one of the forest resource assessment methods established to quantify greenhouse gases (GHGs) emissions [29]. Biomass models are an important tool to facilitate the implementation of MRV in different forest systems across the country. To specifically, site-and species-specific biomass model is a crucial to accurately quantify biomass and carbon stocks by compromising the species diversity and environmental variability. Though efforts have been made to develop species-specific biomass models in dry Afromontane forests e.g. [12, 3034]; models and data which are basic for decision-making and facilitating forest management in dry Afromontane forests in Ethiopia, however, are very scarce [35].

The Menagesha Suba forest is one of the few remaining natural high forests in central Ethiopia [36]. This forest is dominated by multispecies stands, and about 97 plant species grow in this forest area [36, 37]. Of which, 46 species are woody plants [38]. They have been subjected to selective harvesting for more than a century to provide a mix of services to local communities [36, 39]. To preserve the forest and its environmental services, the Ethiopian government designated it a National Forest Priority Area in the 1980s. In addition to timber production, emphasis is being given to the concepts of ecological sustainability and biodiversity [36, 40, 41]. Even though the forest is considered a forest priority area, the current rate of deforestation is extremely high due to clearing for fuelwood, agricultural land expansion, and lumber. As a result, the forest cover was reduced from 7,360 ha to 2,720 ha [37]. Thus, developing a species specific biomass model is imperative to developing the management plan to sustainably conserve this forest. The overall aim of this study was therefore to develop species specific biomass models to estimate biomass and carbon stock of three codominant native species (i.e., Rhus ruspolii Engl, Ekebergia capensis Sparrm, and Nuxia congesta (R. Br. ex Fresen.) in Menagesha Suba dry Afromontane Forest.

2. Materials and Methods

2.1. Description of Study Area

The study was conducted in Menagesha Suba forest located in Oromia regional state, central highlands of Ethiopia, which is about 45 km southwest of Addis Ababa. Menagesha Suba Forest is one of the 58 National Forest Priority Areas (NFPAs) in Ethiopia which need management intervention to maintain the service provided by the forest. The forest supplies a wide range of ecosystem services, including carbon storage. It is also home to many bird species. This forest is also identified as one of the candidate areas for implementing the international climate change mitigation policy, Reducing Emission from Deforestation and Forest Degradation Plus Sustainable Forestry (REDD+). It is situated between 38°31′30″E and 38°34′30″E and 8°57′0″N and 9°0′N (Figure 1) with an altitude of 2200–3385 m a.s.l. [42]. The annual rainfall ranged between 900 mm and 1500 mm, while the mean average temperature of the forest is estimated at 16.5°C.

According to Zewdie’s [43] report, a total of 82 vascular plant species representing 44 families were recorded in the Menagesha Suba forest. Moreover, Feleke [44] reported that 135 vascular plant species belonging to 67 families were recorded in the Menagesha Suba forest. Oromia Wildlife and Forest Enterprise’s (OWFE) [45] report also indicated that about 46 woody species were recorded in Menagesha Suba’s natural forest. Based on the Important Value Index (IVI), Juniperus procera (Hochst. Ex. Endlicher), Olea europea subsp cuspidata (Wall. Ex. DC.), Prunus africana (Hook. f.) Kalkman, Podocarpus falcatus (Thunb.) Mirb, Nuxia congesta, Rhus ruspolii, Ficus sur Forssk., Acacia lahai Steud. and Hochst. Ex Benth, Ficus sycomorus L, and Ekebergia capensis are the top ten dominant woody tree species in the natural forest of Menagesha Suba (Table 1). The understory of the forest is dominated by Olea europea subsp cuspidata, Allophylus abyssinicus (Hochst. Radlk.), Maytenus sp., and Euphorbia ampliphylla Pax, and at higher altitudes, the smaller Juniperus procera Hochst. ex Endlicher is mixed with Erica arborea L., Rosa abyssinica R. Br., and the endemic Asminum stans Pax [43].

The harvest of economically important and endangered indigenous woody tree species such as Juniperus procera, Podocarpus falcatus (Thunb.) Mirb, Hagenia abyssinica, Prunus africana, and Olea europea, which were legally protected by the Forest Proclamation, was prohibited from the cut. Ficus sur, Acacia lahai, and Ficus sycomorus are not accessible to cut. Nevertheless, biomass equations for Olinia aequipetala and Scolopia theifolia Gilg were developed by Tesfaye et al. [33]. Consequently, among the dominant woody tree species, three codominant woody tree species (i.e., Ekebergia capensis, Nuxia congesta, and Rhus ruspolii) were selected for the aforementioned purposes based on their IVI and accessibility to cut down. These species have ecological and socio-economic importance for the surrounding community.

2.2. Sampling Design and Techniques

Menagesha Suba natural forest was stratified to avoid heterogeneity of the natural forest based on diameter distribution, accessibility, and representativeness. The stratification of the forest was done following the distribution of selected woody tree species in the forest. The selected species were not evenly distributed across the forest. Moreover, the number of species in each diameter range also varies from species to species. To accommodate the variation and get the targeted species, systematic random sampling was applied to select the sampled trees from the forest. A total of 45 sample trees and 15 trees for each tree species were selected following the study of Tesfaye et al. [33]. The diameter classes were formulated with 10 cm-diameter intervals. The basal area was used to determine the number of sample trees in each diameter class. Accordingly, 6, 6, and 3 trees were cut from diameter ranges of 5–15 cm, 16–25 cm, and 25–35 cm, respectively for E. capensis, whereas 6, 7, and 2 trees were cut from diameter classes of 5–15 cm, 16–25 cm, and 25–35 cm, respectively, for R. ruspolii, and 5, 8, and 2 trees were selected to be cut from diameter ranges of 5–15 cm, 16–25 cm, and 25–35 cm, respectively, for N. congesta. Representative sample trees were systematically marked [7, 46]. In each sampled tree across the forest, all trees with DBH >5 cm were recorded. Trees with DBH <5 cm were excluded from our data set [10] because such trees hold a small fraction of AGB in forests. Big trees with DBH >35 cm were excluded due to the OWFE prohibition om felling trees from natural forests by OWFE.

2.3. Field Data Collection

Diameter at breast height (DBH at 1.30 m) and diameter at stump height (DSH at 30 cm from the ground) of all tree species were measured before tree felling, whereas total height (TH) was measured after a tree was cut down as used in Tesfaye et al. [33]. The selected trees were felled at a height of 30 cm from the ground. The minimum distance between sampled tree species to be cut was 100 m from each other. The first tree was cut down randomly from one corner of the forest, and the other trees were cut down at a minimum distance of 100 m from each other. The felled tree species were separated into stems (from the stump to the top ≥7 cm diameter), branches (diameter >2 cm), and foliage (diameter < 2 cm + leaves) [33]. The stem and branches were crosscut into manageable pieces (logs), whose length is 2 m, to have a uniform shape to facilitate the volume and weight measurements, respectively. The real volume was calculated by using the Smalian formula. The total green weight of branches and foliage components was measured in the field by a beam balance of 100 kg (±0.01 kg) [47]. About 200 g of foliage and branches were collected for laboratory analysis. Their fresh weight was measured in the field. Three discs with 5 cm thickness were taken at the base, mid, and top of each stem. The fresh weight of each disc was measured at the field. The volume of the discs was measured via water displacement methods.

2.4. Laboratory Measurement

To determine the dry matter contents of the sampled trees, representative subsamples of the three components were taken to the laboratory for dry weight determination. In the laboratory, foliage samples were dried at 70°C, whereas woody parts or collected discs and branches were dried at a temperature of 105°C for all sample categories until a constant weight is reached, as used in [48, 49]. The weight of each dried disc, branch, and piece of foliage was measured using a beam balance weighing 100 kg.

2.5. Data Analysis

Due to several effect variables, multiple linear regression models were used [49] to develop biomass equations from independent variables DSH, DBH, H, and wood density. The summary of the variables used for biomass models is summarized in (Table 2). The general model is written as the following:where Y is the response variable, X1, …, Xp the p effect variables, a0, …, ap the coefficients to be estimated, and ε the residual error. General biomass models used for developing species-specific biomass equations were summarized in Table 2.

Fitting biomass models by using simple linear regression as the variable is generally rarely used in the case of biological data, such as tree volume or biomass because heteroscedasticity is the rule and homoscedasticity the exemption [49]. Therefore, the log transformation was applied to resolve this problem. The best-fitting models for each component were selected by evaluating their goodness-of-fit based on the residual standard error (RSE), their adjusted coefficient of determination (Adj. R2), and their akaike information criterion (AIC) based on the work of [19, 50, 51]. Thus, the model with the lowest value of RSE and AIC and the higher values of the adjusted coefficient of determination (adj. R2) were selected as the best-fitted biomass and volume regression models. All statistical analyses were carried out with the R3.5.3 software package. The AIC is as follows:where n: number of samples, RSS: the residual sums of squares, and K: parameters of the model including the p the parameter for error estimates (for example, the model y = a + bx, then k = 3).

Relative bias (RB), average deviation (S%), relative root mean square error (rRMSE), and a paired t-test were used to compare the predictive accuracy of the general equations developed for tropical dry forests by [4, 10, 19] with species-specific fitted models. Average deviation S% is used to evaluate variation between observations and predictions. The smaller S% value is preferredwhere Yilt: the predicted biomass; Yi: the observed biomass/carbon; n = number of observations

Then, the wood density was calculated as dry oven weight divided by the green volume of the samplewhere WD represents wood density, Wods oven-dry weight of the sample, and Vfresh.s fresh volume of the sample.

Besides, an ANOVA was conducted to know whether there is variation in wood density among species and between trees’ components or not.

2.6. Model Comparison

Model comparison was made using allometric models that were developed by [4, 10, 19]. Moreover, the model’s adequacy and precision were evaluated and tested using sampled trees. These models are frequently used in the tropical forest. For this reason, the models were selected for comparison.

3. Results

3.1. Wood Density

The stem wood densities of the targeted species were summarized in Table 3. The mean wood density of R. ruspolii ranged from 0.45–0.90, that of E. capensis ranged from 0.43–0.92, and that of N. congesta ranged from 0.33–0.70. Thus, there was statistically significant variation among selected woody tree species (i.e., R. ruspolii, E. capensis, and N. congesta).

3.2. Tree Component Biomass

Total AGB of each selected species is the summation of stem biomass (SBM), branches biomass (BBM), and foliage biomass (FBM). About 75.67, 69.6, and 64.96% of total above-ground biomass was accumulated in the stem component of E. capensis, R. ruspolii, and N. congesta, respectively, whereas 3.61, 6.47 and 4.44% of biomass was stored in the leaf part of E. capensis, R. ruspolii and N. congesta, respectively, and 20.74, 23.90, and 30.60% were stored in the branches of E. capensis, R. ruspolii, and N. congesta, respectively. Significant variation was observed between all components (), except between branch and foliage biomass of R. ruspolii and E. capensis (Table 4).

3.3. Correlation between Tree Components’ Biomass and Dendrometric Variables

When correlating the variables, we found strong correlation between AGB and DBH, AGB and DSH, DBH and DSH, and TH and DBH (; Corr. > 0.73; Figures 2(a)2(c)) for all species, whereas the correlation between DSH and TH was medium in the cases of N. congesta and R. ruspolii and weak for E. capensis. Wood density was not correlated with all variables, including AGB.

3.4. Biomass Regression Models

A total of eight fitted equations were selected for R. ruspolii, E. capensis, and N. congesta (Table 5). Out of tested models, model 8 was found the best-fitted model for R. ruspolii (Adj. R2 = 0.91, ) and model 5 for N. congesta (Adj. R2 = 0.97, ) whereas model 6 was the best fitted model for E. capensis (Adj. R2 = 0.95, ). All tested models were statistically significant (). In most cases, the coefficients β1, β2, β3, and β4 were positive, which indicates that the aboveground biomass of each tree species would increase for every 1 unit increase in the value of other tree parameters used in the model.

Residual plots of the best-fitted height equations for our studied species generally indicate an even spread of residuals above and below the zero line with no systematic trend (Figure 3). This suggests that natural log-transformed multiple linear regressions are effective in stabilizing error variance, and this study appeared to be appropriate for reducing heteroscedasticity.

The best-fit models from all tested models were also tested for accuracy based on observed and predicted data. Observed and predicted above-ground biomass values are close to the 1 : 1 line (Figure 4).

3.5. Model Comparison

The mean observed values of R. ruspolii, E. capensis, and N. congesta were 111.9, 125.94, and 191.53 kg, respectively. Whereas the average predicted biomass of R. ruspolii through the species-specific model, Brown [4]; Chave et al. [10, 19] were 125.71, 111.40, 155.94, and 185.53 kg, respectively. The average predicted biomass of E. capensis via a species-specific model Brown [4]; Chave et al. [10, 19] was 190.69, 153.51, 247.30, and 309.77 kg, respectively. Consequently, the average predicted biomass of N. congesta through the species-specific model of Brown [4]; Chave et al. [10, 19] was 110.80, 134.54, 146.18, and 142.26 kg, respectively. Thus, the variation between observed and predicted biomass through a species-specific allometric equation was statistically not significant for all targeted species, while the variation between observed and predicted biomass via Chave et al. [10, 19] was statistically significant (Table 6). On the other hand, the estimated biomass using the biomass equation developed by Brown [4] underestimated the biomass of R. ruspolii and E. capensis but overestimated that of N. congesta. Except for R. ruspolii, the variation between the predicted biomass using Brown [4] and the observed one was statistically significant ().

4. Discussions

Estimating wood density (WD) for these co-dominant native tree species in Ethiopian dry afromontane mixed forests will contribute to developing biomass models and sustainable management of these forests. Wood densities were determined for each targeted species. The mean wood density of each species ranged from 0.5–0.69 g/cm3. The WD of E. capensis reported by Desalegn et al. [52] is lower by 1% than that of this study. Consequently, the average WD of Rhus ruspolii was 0.51 [28], which is 18% lower than this study. The average WD of Nuxia congesta reported by MEFCC [28] was 12% higher than that of this study. The wood density of tropical woody species is between 0 and 1.5 [53] and 0.262–1.040 g/cm3 in Ethiopia [28]. Thus, the mean average WD of each targeted species is found within this range. Furthermore, significant variation in WD was observed among selected tree species. This indicated that the anatomical structure of these species could be varied. Thus, determining the WD of targeted species has an indispensable role in estimating the above-ground biomass of Ethiopian’s dry Afromontane forest. Tesfaye et al. [54] stated that determining WD for each species is essential rather than using a default value during the development of the biomass model. Henry et al. [55] also stated that tree species and stem position in a given forest influence wood density. For instance, Olmos and Lorena [56] found higher values of wood density in the bottom parts of a tree than in the top parts of a tree. Species with very thick-walled fibers had a significantly higher mean wood density than species with very thin-walled fibers [53].

We estimated the biomass of the stem, branches, and foliage of E. capensis, R. ruspolii, and N. congesta. For all targeted species, a larger proportion of carbon and biomass was stored in the stem than in the foliage and branches. Previous studies also confirmed that a large portion of biomass accumulated in the stem compared to other tree compartments [33, 55, 5762]. For instance, larger biomass was accumulated in the stem part of Casuarina equisetifolia than in the branches, leaves, and roots [63]. Thus, to increase the above-ground biomass, different silvicultural activities which enhance stem growth become imperative. In contrast to this study, the highest biomass fraction was observed in branches of Afzelia africana [58], Scolopia theifolia [33], Combretum glutinosum, and Terminalia laxiflora [64]. Thus, the variation of biomass among components of the woody plants could be due to plant architecture and morphology, climatic and edaphic factors, and forest management practices [7, 54, 58, 65, 66], availability of nutrients [67], age [68, 69], and spacing [70]. For example, annual maximum temperature and precipitation had a positive impact on stem biomass [71].

Developing biomass equations for tropical dry forest species is essential to quantifying biomass and carbon stock, which encourage sustainable forest management of dry Afromontane forests [33]. The biomass models developed in this study included DBH, DSH, TH, and WD as independent variables for all selected species. Above-ground biomass (AGB) of all targeted species was highly correlated with the aforementioned independent variables, except wood density. In line with this finding, the AGB of five native tree species was correlated with DBH, CH, and TH [54]. Previous findings also supported the idea that combinations of different dendrometric variables fit better than individual variables [12, 18, 33, 72]. For instance, DBH, crown diameter, and basal diameter were important variables to explain the biomass of Vitellaria paradoxa [73, 74]. Thus, combining these variables provides better fit results than using DBH and TH alone to accurately estimate AGB.

Despite some studies suggesting the importance of generic equation to estimate AGB of tropical African forests [5, 10, 19], others argued generic models are inappropriate for African tropical forests [11, 18, 30, 33]. Moreover, Goodman et al. [73] also emphasized that pan-tropical allometric models may underestimate the AGB of a very large tree by 20% due to sampling bias in the harvest data set. A comparison of generalized models [4, 10, and 19] to the fitted models for the species studied revealed that generalized models were not appropriate for all species to estimate biomass as compared to a newly developed model. All generalized models tested showed a high bias and were unsuitable for biomass estimation of selected tree species. The variation between observed biomass and predicted biomass through a generalized equation is due to unfair representation or a lack of sample trees from Ethiopia in their data set. For example, Tesfaye et al. [33] indicated that the generalized allometric models by Brown et al. [5] showed the poorest results, with 32–59% average deviation for AGB predictions of five tree species in Ethiopia. Correspondingly, the model developed by Chave et al. [19] was inappropriate for three species in Ethiopia, including Allophylus abyssinicus, Olinia rochetiana, and Rhus glutinosa [33]. Thus, compared to models developed by Brown [4] and Chave et al. [10, 19], the species-specific models can generate higher reliability biomass estimation for the studied species.

Overall, generalized models produce a larger bias for all species when estimating above-ground biomass. Due to the species diversity in tropical forests, species-specific biomass models are required to accurately estimate the biomass. Furthermore, the equations developed in this study can be used for estimating forest carbon stocks, identifying carbon sink capacity, establishing carbon trade value, and informing management policies. This study is the pioneer in the study area, and the models developed represent useful tools to support decision-making for selected ecologically and economically important tree species in Menagesha Suba natural forests. Thus, the biomass models developed here may help with the sustainable management of the dry forest of central Ethiopia.

5. Conclusions

We explored the applicability of species-specific biomass models to estimate the aboveground biomass of three economically and ecologically important indigenous woody species in Menagesha Suba, the central highlands of Ethiopia, which helped us understand the strategies to sustainably manage the dry Afromontane Forest. Furthermore, this study was used to indicate the gap in generic biomass equations used to estimate above-ground biomass in the study area. Moreover, the information generated by this study can be used for sustainable forest management, including REDD+ and MRV practices, in the dry Afromontane Forest of Ethiopia. The predictor variables used to develop biomass models of Ekebergia capensis, Nuxia congesta, and Rhus ruspolii were DBH, DSH, H, and WD. Despite the development of species-specific biomass being complex, the developed models are accurate in estimating the biomass of studied species in the study area. Including more variables improves the accuracy of the models to estimate the biomass. As compared to the general equation, the species-specific biomass equations measured biomass and carbon stocks accurately. The outputs of the model comparison showed that the application of generalized models for estimating AGB of selected species produced a very big bias. Generic equations highly deviated from the observed value, whereas species-specific equations deviated less from the actual value. Therefore, species-specific equations can accurately measure biomass as compared to general equations.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest on this work.

Acknowledgments

The authors acknowledge Ethiopian Environment and Forest Research Institute for its financial support. Shimelis Telila, Tewodros Andualem, Ashenafi Ayalew, and Kasahun Guluma were also acknowledged for their contribution to this work.