International Journal of Forestry Research

International Journal of Forestry Research / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 1876329 | https://doi.org/10.1155/2019/1876329

Bernardol John Manyanda, Wilson Ancelm Mugasha, Emannuel F. Nzunda, Rogers Ernest Malimbwi, "Biomass and Volume Models Based on Stump Diameter for Assessing Degradation of Miombo Woodlands in Tanzania", International Journal of Forestry Research, vol. 2019, Article ID 1876329, 15 pages, 2019. https://doi.org/10.1155/2019/1876329

Biomass and Volume Models Based on Stump Diameter for Assessing Degradation of Miombo Woodlands in Tanzania

Academic Editor: Scott D. Roberts
Received29 May 2018
Revised11 Nov 2018
Accepted27 Dec 2018
Published15 Jan 2019

Abstract

Models to estimate forest degradation in terms of removed volume and biomass from the extraction of wood fuel and logging using stump diameter (SD) are lacking. The common method of estimating removals is through estimating diameter at breast height (D) by applying equations relating measured D and SD. The estimated D is then used to estimate biomass and volume by means of allometric equations, which utilize D. Through this sequence of procedures, it is apparent that there is an accumulation of errors. This study developed equations for estimating volume, aboveground biomass (ABG), and belowground biomass (BGB) using SD in miombo woodlands of mainland Tanzania. Volume models were developed from 114 sample trees while AGB and BGB models were developed from 127 and 57 sample trees, respectively. Both site specific and regional models were developed. Over 70% of the variations in BGB, AGB, and volume were explained by SD. It was apparent that SD is inferior compared to measured D in explaining variation in volume and BGB but not AGB. However, the accuracy of BGB and volume estimates emanating directly from SD were far better than those obtained indirectly, i.e., volume or BGB estimates obtained from estimated D from SD, since the latter is affected by accumulation of regression equation errors. For improved accuracy of ABG, BGB, and volume estimates, we recommend the use of site specific models. However, for areas with no site specific models, application of regional models is recommended. The developed models will facilitate the addition of forest degradation as a REDD+ activity into the forthcoming FREL.

1. Introduction

Miombo woodlands are lands dominated by deciduous trees of the genera Brachystegia, Julbernadia, and Isoberlinia [1, 2]. They cover an area of approximately 2.7 million km2 equivalent to 9% of the African land area spanning across ten countries in southern and eastern Africa including Tanzania [1, 36]. In Tanzania, woodlands cover about 93% of the total forest area of 48.1 million ha [7, 8]. Miombo woodlands can be divided into dry (annual precipitation <1000 mm) and wet miombo (annual precipitation>1000 mm). Dry miombo are found in Zimbabwe, central Tanzania, southern areas of Mozambique, Malawi, and Zambia [6]. Wet Miombo are found in eastern Angola, northern Zambia, central Malawi, and south western Tanzania [6].

Miombo woodlands offer both direct and indirect benefits. Direct benefits include energy (fuel wood and charcoal), construction and craft materials, medicines, employment, income, food (fruits, mushrooms, honey, and edible insect), and fodder for animals. These benefits are either for domestic consumption or for local sale [3, 5]. On the other hand, indirect benefits encompass soil nutrient inputs through nutrient cycling and through nitrogen fixation and environmental services such as soil and water conservation, biodiversity, and carbon sequestration [912].

Miombo woodlands are under pressure from increasing demands for land-based products and services. This has led to forest degradation and deforestation [13, 14]. Forest degradation jeopardizes the capacity of forests to function as regulators of the environment. Consequently, increasing flooding, erosion, reduced soil fertility, and loss of plant and animal diversity have been common [13]. Sustainable provision of goods and services from the miombo woodlands requires effective forest management efforts, which ultimately may make a significant contribution to national goals for Reduction of Emission from Deforestation and Forest Degradation “plus,” the role of conservation, sustainable management of forests, and enhanced carbon stock (REDD+). Implementation of REDD+ as a carbon credit market approvals among others requires measurement and monitoring of carbon emissions from forest degradation and deforestation.

While deforestation is relatively easy to estimate, forest degradation is more challenging [8]. Recently, Tanzania has engaged in developing Forest Reference Emission Level (FREL) standards [14]. Deforestation and conservation are the only REDD+ activities among the five activities, which have been included in the Tanzania FREL. Degradation has not been included due to inadequate data for establishing baselines and monitoring. Forest degradation is taking place throughout the country. Monitoring degradation by means of remote sensing techniques poses a significant challenge since degraded forests frequently maintain a closed canopy. The main drivers of forest degradation are extraction of wood fuel (charcoal and firewood), logging, grazing, and wildfire. Techniques to estimate forest degradation need to be developed following a stepwise approach. One of the techniques is the use of Stump diameter (SD) which may be included during forest inventories as a means of assessing forest degradation.

Diameter at breast height (D) and total tree height have often been used as standard predictors of biomass (both above- and belowground biomass) and volume [35, 15]. This is because these variables are highly correlated with biomass and volume. The common method of estimating removals is through estimating D by developing equations relating D and SD from sample trees measured for both D and SD [1618]. The estimated D is then used to estimate biomass and volume by means of allometric equations. Through this sequence of procedures, it is apparent that there is an accumulation of errors from the estimation of D to the estimation of biomass or volume [16, 1921]. In this case, it is considered more precise to estimate biomass and volume directly from SD. In addition, the recent National Forest Monitoring and Assessment (NAFORMA) included SD for monitoring forest degradation. Therefore, having biomass/volume – SD allometric equations will facilitate the estimation and monitoring of forest degradation in Tanzania associated with extraction of wood fuel and logging.

To date, there is only one biomass/volume – SD allometric equation that was developed from limited sample trees (30 trees) from only a single site in Tanzania [22]. Topography, soil, climate, and species are the main factors, which affect stand architectural variability in miombo woodlands [22]. Therefore, adequate sample trees collected from different sites are imperative to cover the variability in SD, D, and H of trees. Thus, this study developed aboveground biomass (AGB), belowground biomass (BGB), and volume-SD allometric equations that utilize SD as the sole predictor for estimating removal and emission from forest degradation associated with extraction of fuel wood and timber in miombo woodlands.

2. Materials and Methods

2.1. Study Area Description

The data for developing biomass and volume models were collected in Tanzanian miombo woodlands of Manyara, Lindi, and Tabora regions (Figure 1). In Manyara region, data were collected in Ayasanda and Duru Haitemba Village Land Forest Reserves. The dominant miombo woodlands tree species in these reserves include Albizia versicolor, Brachystegia microphylla, Julbernardia globiflora Brachystegia spiciformis, Brachystegia boehmii, Combretum collinum, Parinari curatellifolia, Markhamia obtusifolia, Tamarindus indica., and Senegalia nigrescens (Acacia nigrescens). Likewise, Brachystegia spp., Julbernadia spp., and Pterocarpus angolensis are common miombo woodlands trees found in Angai Village Forest Reserve in the Lindi region. Furthermore, in Tabora region, the data were collected in Nyahua Forest Reserve in which Pterocarpus angolensis, Afzelia quanzensis, Dalbergia melanoxylon, Burkea africana,Pterocarpus tinctorius, and Swartzia madagascariensis are the common tree species. Details on locations and conditions of the study sites are described in Table 1.


RegionForestLocationDominant soil TypeAltitude (m)Mean annual temperature (°C)Mean annual rainfall (mm)
MinMax

ManyaraAyasanda and Duru Haitemba4°20′S, 35°47′EClay alluvial soils1 300–1 8001526854

LindiAngai Villages Forest Reserve9°47′S, 37°55′ESandy loam soils330–6002031873

TaboraNyahua Forest Reserve5°18′S, 32°58′ESandy clay loam soils1 096–1 1031730771

The weather conditions for all sites may be divided into three categories, i.e., a hot dry season from mid-August to the end of October, a hot wet season from November to the beginning of April, and a relatively cool dry season from April to mid-August. Furthermore, two rainfall regimes exist. In the southern, southwestern, central, and western parts of the country, including Lindi and Tabora, the rainy season starts in mid-November and ends in mid-May. In the north and in the northern coastal zone, including Manyara, the rain is distributed over two shorter periods (October–December and March–May).

2.2. Data Collection
2.2.1. Sampling Design

This study utilized a data set that was collected to develop biomass and volume models utilizing D and H [3, 5]. Additionally, SD of the sampled trees was measured. The SD was measured over bark immediately under the cutting point (felling cut). When a tree was irregular, the SD was measured at a higher point where a regular shape commenced but must be below D. Circular sample plots of 15 m radius were distributed within the study sites where 40 plots in each site were systematically established making a total of 120 plots across all three sites. Tree size and species distribution of this inventory were used as the target distribution for the selection of sample trees. One or two trees in each plot were selected for destructive sampling in order to match the target distribution. Some of the sampled trees were selected outside the plot area, i.e., large trees that were not present inside the plot.

2.2.2. Destructive Sampling

Before felling, species of each sample tree was recorded and the tree was measured for D and SD. The SD was measured at least 30 cm from the ground. In case of irregularity of the stump, we moved further above up to a point where irregularity ends. However, the measured point was always breast height. The selected sample trees were divided into two main parts, i.e., aboveground and belowground components. The aboveground component included all biomass above a point of SD measurement.

Aboveground Component. The aboveground component (Table 2) was divided into three sections, merchantable stem and branches including tops (up to a minimum diameter of 2.5 cm) and twigs (with diameter less than 2.5 cm). For small trees with D < 10 cm, no merchantable stem part was considered. For trees with D ≥ 10 cm, no specific minimum diameter was set to distinguish between merchantable stem biomass and branch biomass, but the decision was based on a subjective judgment of the length of the stem that could be used to produce timber. Leaves were excluded from twigs and thus not included in the modelling. Stems and branches were trimmed and cross cut into manageable sections ranging from 1 to 2.5 m in length and weighed for green weight. Mid-diameter and length were measured for each log section. At least two sample disks (depending on the length) from stem and branches were extracted and weighed for determination of the dry to green weight ratio (DG-ratio). Twigs were collected into separate bundles and the green weight of each twig was recorded. Small sample disks were collected from each bundle, labelled, and measured for green weight ready for drying in laboratory.


ComponentSiteStump diameter (cm)
nMeanMinMax

VolumeManyara3637.08580
Lindi3934.254.582
Tabora3937.051.8102
All11436.101.8102

AGBManyara4038.312.480
Lindi4738.481.9114
Tabora4037.271.8102
All12738.051.8114

BGBManyara1933.25555
Lindi1936.6110.272
Tabora1948.9413102
All5739.605102

Belowground Component. For the belowground component (57 sample trees) (Table 2), areas around the sample tree were excavated to expose all of the roots emanating from the root crown. Three main roots (largest, medium, and smallest in diameter) were selected and excavated in full, measured for diameter at the branching point from the root crown, and weighed. Up to three side roots were selected from the excavated main roots, measured for diameter at the branching points from the main root, and weighed. Other side roots from the excavated main roots were measured for diameter at the branching point. For trees with a tap root, excavation was done until reaching a diameter corresponding to the largest selected main roots (down to 4 m depth). Diameter was measured and the remaining portion was treated as a side root (Figure 2). All main roots that were not excavated were measured for diameter at the branching point of the root crown. The root crown was also recorded for green weight. Details on excavation and sampling procedures for BGB are described in [3].

In order to obtain estimates of the dry weights of the belowground components, at least two samples were taken from all main and side roots and at least two from the root crown. All were weighed for green weight, labelled, and prepared for the laboratory procedures. In the laboratory, all above- and belowground samples were oven-dried to a constant weight at 105°C for stems, branches, and twigs for at least 48 hours and weight was recorded after every 6 hours until they maintain constant weight.

2.3. Data Analysis
2.3.1. Determination of Tree Volume

Each log section volume was calculated by multiplying the cross-sectional area at the midpoint of each log by its length. The volume of the merchantable stem and branches for a tree was obtained by summing the volumes of the respective sections for that specific tree. Total tree volume was finally obtained through summation of the volumes of the merchantable stem and the branches. Summary statistics for total observed volume over sites are shown in Table 3.


SiteAGB (Kg)BGB (kg)Volume(m3)
nMeanMinMaxnMeanMinMaxnMeanMinMax

Manyara4011260.775143191992.90643361.260.00157.87
Lindi4714630.12104181930817.201130391.210.00165.79
Tabora4011490.1987901948618.722355391.420.00027.53

All12712580.1110418573312.8923561141.300.00027.87

2.3.2. Determination of Aboveground Biomass

Dry to green weight (DG) ratios were determined for each tree component, i.e., stem, branches, and twigs [3]. The mean DG for each tree component was computed. The dry weight of each component (merchantable stem, branches, and twigs) was obtained as a product of mean DG-ratio and the green weight of the respective tree component. Finally, total AGB was computed as the sum of the dry weights of merchantable stem, branches, and twigs. Summary statistics for the AGB for each site are presented in Table 3.

2.3.3. Determination of Belowground Biomass

Dry to green weight (DG) ratios were determined for each belowground components, i.e., root crown, main roots, and side roots [3]. Since few sample roots were measured for green weight in the belowground section, it was necessary to develop a green weight-root diameter relationship to be able to establish the green weight of unexcavated roots, which were only measured for diameter. Once the green weight of each root was determined, green weight of each section (main roots, side roots, and root crowns) was converted to biomass by multiplying its green weight and its respective DG. The BGB of a tree was obtained by summing the dry weights of the roots and that of the root crown. Details on the procedures and uncertainty associated with these equations are described in [3]. Summary statistics for BGB for each site are shown in Table 3.

2.3.4. Model Development

Prior to model fitting, and to avoid blind fittings, response variables were plotted against the explanatory variable to examine the patterns and extent of variance for each site and for the combination of sites. The scatter plots (Figures 3, 4, and 5) displayed a positive nonlinear relationship between volume/biomass and SD. As a rule of thumb, a best model emanates from fitting several models since tree allometry differs due to biotic and abiotic factors [35, 19, 23, 24]. Since all trees during data collection were cut to a height of at least 30 cm above ground level, the stump height was not included in the models as an independent variable. We fitted three model forms, i.e., one nonlinear and two linear models. where Y is volume (m3 tree−1), AGB, or BGB (kg tree−1); SD is stump diameter (cm); Ɛ is the error term (or residual); and a, b, and c are the regression parameters to be estimated.

Modelling data used in this study represent a hierarchical structure, i.e., three sites and tree species nested into sites. In this case, a nonlinear mixed effect modelling approach was considered ideal for developing predictive models that would account for dependence of the species within sites and for preservation of original scale. The list of species is presented in the Appendix. Preliminary findings showed a significant improvement of model with mixed effect compared to models without mixed effect. In addition, to account for variation (i.e., heteroscedasticity due to SD) not accounted for by the random effects, we also included a power variance function structure, i.e., varPower, implemented in the nlme packages of R software for nonlinear and linear models [25, 26]. Model (1) (nonlinear) was fitted with an nlme function while models (2) and (3) (linear models) were fitted with the lme function in the nlme package [27] in R software (R Development Core Team 2018). We developed regional as well as site specific AGB, BGB, and volume models. Models were fitted by allowing random effects on all model parameters (see (1)). With the inclusion of random effects, (4) emerged. where represents model parameters a, b, and c and expresses the difference in parameter of site i and tree species j from the mean value obtained from (1), (2), and (3) or typical site or tree species and Ri is the random effects.

2.3.5. Model Selection and Evaluation

To select the best models for volume, AGB, and BGB, we computed Root Mean Square Error (SE%), Coefficient of Determination (R2), Mean Prediction Error percentage (E%), and Akaike Information Criterion (AIC). Models with small values of SE%, E%, AIC, and high R2 were considered to have a good fit. Mean prediction error percentage (E%) was computed using (5) based on fixed effects parameters only. However, with an increase in model parameters, SE% tends to be smaller and R2 becomes higher regardless of the contribution of the added parameters. Therefore, to address this problem we used AIC to select the best model (see (6)). The selected models were further examined by using residual plots. where Y is observed AGB, BGB, or volume; is predicted AGB, BGB, or volume; and n is the number of observations (i.e., number of trees). where L is the value of the likelihood and k is the number of parameters.

In addition, we tested the previously developed volume, BGB, and AGB models from similar data sets except that their independent variables were D and total tree height. The aim was to gauge the strength of SD relative to D. We further evaluated the performance of regional models to the studied sites using E%.

3. Results

3.1. Volume Models

Table 4 presents parameter estimates and performance criteria of the fitted volume models. Parameter estimates of all fitted models were significantly different from zero (p < 0.05). Model (1) was consistently the best in all cases, i.e., smallest AIC and E%, and therefore was selected for further evaluation. For the selected models, the coefficient of determination (R2) was found to be greater than 70% except for Tabora (62%). With the exception of Manyara (absolute E% = 0.85%), E% were found to range from 9.6% to 11.8%. Residual plots of selected models are shown in Figure 6. Residual plots did not show any pattern that indicates model bias.


SiteModel formsModel coefficientsAICSESE%E%R2
abc

Tabora0.0000312.78805-41.00.7245.7-11.790.62
-0.002340.00054-1.2571.0163.93-67.20.79
0.00443-0.00340.00059-24.660.9861.96-14.50.80

Manyara0.0000352.774061.3480.5142.10-0.850.89
-0.01290.0005023.851.0688.12-13.550.536
0.0205-0.00770.000799.750.7965.04-9.510.75

Lindi0.0000262.8908-39.20.8156.2-9.60.789
-0.01370.000592.1320.8861.39-17.00.749
0.0385-0.012370.00098-27.20.5336.6-9.800.91

Regional0.0000322.7992-90.161.0171.7-10.50.709
-0.002940.0005630.70.9970.4-48.30.719
0.005570.00410.00066-41.80.9164.6-15.30.76

Note: Bold type indicates selected model.
3.2. Above- and Belowground Biomass Models

The parameter estimates and performance criteria of the AGB and BGB models are presented in Tables 5 and 6, respectively. Model (1) was consistently the best for AGB in all cases based on the lowest E% and AIC (Table 7). Since selection was based on models with smallest AIC, model (1) was selected for all cases. Coefficient of determination (R2) of the selected regional AGB model was 0.92, while R2 for the selected site specific models ranged from 0.88 to 0.94. The E% of the selected models ranged from 9.7% to 11.1%. Residual plots of selected AGB models are shown in Figure 7. Residual plots did not show any pattern to indicate model bias.


SiteModelModel coefficientsAICSESE%E%R2
abc

Tabora0.03912.6491485.9483.2242.9-11.10.925
-1.36260.3755524.9986.7287.6-39.60.687
3.0632-2.40730.4559501.8848.1775.3-13.40.76

Manyara0.06052.5490511.7649.6144.49.30.88
-1.45510.35343544.31319.7090.2-10.30.52
5.5304-3.18940.5021516.31035.8670.8-9.70.70

Lindi0.028212.7508552.8482.8442.0-9.70.94
-0.78730.2337620.51640.49142.7-10.30.28
1.9440-1.66670.3871598.91300.21113.18.30.55

Regional0.037852.67001541.1551.0043.8-7.90.92
-1.10670.31941679.61419.02112.8-14.90.49
2.9345-2.27790.45511614.61124.6589.4-8.960.68

Note: Bold type indicates selected model.

SiteModel formsModel coefficientsAICSESE%E%R2
abc

Tabora0.08822.1386222.887.6944.05-13.590.86
-10.17270.1590228.290.8845.65-13.310.85
-5.8520-0.37810.1651230.289.9845.2-13.310.85

Manyara0.14342.0144199.6122.5239.72-8.360.78
0.70400.1546199.0122.1239.59-9.270.79
-8.95451.90070.0945196.2141.9646.02-3.180.71

Lindi0.081472.1896211.2180.7137.18-13.420.89
0.15530.1454212.6259.7853.456.150.77
41.0525-4.93790.25975212.6157.4732.4-10.840.91

Regional0.10562.1035623.5154.0146.5-12.20.86
-0.80400.1526624.9164.6749.72-11.410.84
0.9476-0.36880.1631626.5158.9848.0-12.040.85

Note: Bold type indicates selected model.

SiteComponentReferenceEquationE

TaboraAGBMugasha et al. (2013)-3.9
BGBMugasha et al. (2013)-2.7
VolumeMauya et al. (2014)2.0

ManyaraAGBMugasha et al. (2013)-12.5
BGBMugasha et al. (2013)-11.8
VolumeMauya et al. (2014)1.1

LindiAGBMugasha et al. (2013)-10.1
BGBMugasha et al. (2013)-8.8
VolumeMauya et al. (2014)0.2

RegionalAGBMugasha et al. (2013)-8.9
BGBMugasha et al. (2013)-7.8
VolumeMauya et al. (2014)-0.6

All parameter estimates for BGB models were significantly different from zero. Except for Manyara where model (3) had a good fit, model (1) performed well in Lindi and Tabora and for the regional model (Table 6). These models had lower AIC and lower E% values and therefore were selected for further evaluation. The selected BGB model for Tabora, Manyara, and Lindi had R2 of 0.86, 0.71, and 0.86, respectively. The E% of the selected models ranged from 3.2% to 13.6%. Residual plots of selected models are shown in Figure 8. Residual plots did not show any pattern indicating model bias.

3.3. Evaluation of Previously Developed Models

The prediction power of the previously developed volume and biomass models for Miombo woodlands which utilize D as sole independent variable was also tested on the current data set (Table 7). The tested volume and biomass models were developed from similar data [3, 5]. The E% for AGB, BGB, and volume model over the sites ranged from -2.7 to 11.8, -3.9 to 12.5, and 1.1 to 2.0, respectively. These values were slightly higher than those produced by the current developed models.

3.4. Performance of Developed Regional Models to Study Sites

Performance of the regional models relative to study sites results is presented in Table 8. The findings show that regional models E% were not significantly different from zero (p > 0.05) except for Tabora region.


SiteE%
BGB AGBVolume

Manyara-8.0-7.7-9.1
Lindi-9.5-10.5-0.7
Tabora-19.4-21.231.4

4. Discussion

4.1. Modelling Data

This study developed volume and biomass models using a comprehensive data collected from three regions rich in miombo woodlands in Tanzania, i.e., Tabora, Lindi, and Manyara. Sampling was guided by the observed species distributions from previous systematic sample plot inventories carried out in each site so that the most frequently occurring genera of miombo woodlands, such as Brachystegia, Pterocarpus, Julbernadia, and Combretum, were as representative as possible (see the Appendix). Such considerations were not made in the development of models [18, 22] that are currently applied for removed volume and biomass estimation in the miombo woodlands of Tanzania. The data covered a wide range of growth conditions such as climate, topography, and soils. In addition, the data also covered a large number of observations (AGB = 127; BGB = 57; and volume = 114) and wide range of tree sizes (minimum SD of 1.8cm to maximum SD of 102 cm) for modelling which ensured that the majority of volume and biomass variation is explained by SD. Trees of large sizes were also included to avoid extrapolation beyond the data ranges. The inclusion of larger trees is of particular importance for mature forests because large trees account for the largest part of the volume and biomass and drives the model fits [3, 5, 6]. With such data set characteristics, we are confident that the developed models are superior to the previously developed volume and biomass model from a single site and limited data in eastern Tanzania that also utilized SD as a sole predictor [21, 22].

4.2. Volume Models

Model (1) was selected in all cases (site specific and for regional models). Similar model forms had also fitted well to the same data utilizing D as the explanatory variable [3, 5]. When comparing R2 with other studies, [21] found R2 ranging from 0.70 to 0.81 which are in line with those obtained in the present study. This also suggests that SD may explain over 70% of variation in tree volume.

The lower E% which was also not significantly different from zero indicates that the developed models were unbiased (Table 4). Generally, based on E%, models for individual sites were superior to regional models as expected since by combining all data sets there is an increase of volume, AGB, and BGB variations which were not able to be explained by only SD [3, 5]. Evaluation of previously developed volume, AGB, and BGB models shows that E% was slightly lower than those of earlier models except for Manyara. This trend suggests that D is a slightly better explanatory variable than SD when estimating standing trees volume. However, when estimating removed volume, the most viable option is to utilize a model with SD as sole predictor. Furthermore, the findings show that application of regional models for all components was limited to Manyara and Lindi where prediction errors were not significantly different from zero. It is not clear why prediction was poor in Tabora. However, differences in tree allometry, especially at root collar, may explain the difference.

4.3. Aboveground Biomass Models

Model (1) was selected in all cases with lower AIC and E%. This goodness of fit and flexibility of this model form have been reported by other scholars [3, 5, 6, 22, 24]. It is apparent in most cases that models are applied outside their modelling data. However, models which are simple (e.g., model (1) with few coefficients) are flexible and mostly unbiased relative to more complex models (with large number of coefficients) which are more specific and therefore biased when applied outside their range [28].

The selected site specific and regional models for AGB had R2 values ranging from 0.91 to 0.93. These values are similar to previously reported values for models developed for miombo woodlands that utilize SD as the predictor [22] and lowland forest in Dindili forest reserve in Morogoro that utilize D as the predictor [29] with R2 ranging from 83% to 89%. This implies that, for miombo woodlands, SD is adequate to explain majority of variation in AGB, as does D. This is further supported by residual plots, which did not show any indication of model bias. In addition, E% obtained from evaluation of previous model to the current data set were similar to the E% obtained for the models developed in the current study (Tables 5 and 7). Therefore, the conventional approach of estimating degradation (trees removal) by first estimating D from SD and then applying estimated D to estimate AGB may jeopardize the estimates accuracy more than using SD directly as suggested in this study.

4.4. Belowground Biomass Models

Model (1) for all sites and for the regional models was selected as the preferred model except for Manyara region. The inconsistency of model (1) performance for Manyara region was also reported by [3] for BGB. In comparison with R2 values reported by [3], i.e., R2 ranging from 0.89 to 0.94, R2 reported in this study was relatively lower. Similarly, E% values also followed that trend. For example, the current regional model had an E% of 12.2% while the previous regional model using D as the predictor had an E% of -7.8%. This suggests that not all variation in BGB explained by D can adequately be explained by SD. In addition, the accuracy of BGB estimates dropped when using SD as the predictor rather than D.

There are several studies, which have modelled BGB with SD as the explanatory variable to make a concrete comparison. However, [30] reported findings that are contrary to the current findings. They found in Pinus densiflora in Samcheok that root collar diameter had a higher correlation with BGB than D. The fact that this study dealt with multiple tree species growing in different environmental conditions, i.e., soil depth, different moisture stress (climate) [30], and topography, may explain the inconsistency between the two studies.

5. Conclusion

Predicting volume and biomass directly from SD is useful in situations where D is not available. This study developed robust volume, AGB, and BGB models utilizing SD as the predictor that can be used to estimate forest degradation carried out through tree cutting in miombo woodlands. These models will facilitate the addition of forest degradation as a REDD+ activity in the forthcoming FREL. It was apparent that SD was inferior to D in explaining variation in volume and BGB but not AGB. However, the accuracy of BGB and volume estimates emanating directly from SD are far better than those obtained indirectly, i.e., volume or BGB estimates obtained from estimated D from SD, since the latter are affected by the accumulation of regression equation errors. For improved accuracy of ABG and BGB and volume estimates, it is recommended that site specific models be used. However, for areas where no site specific models exist, application of regional models is recommended.

Appendix

See Table 9.


SiteTree speciesSD (cm)

ManyaraBrachystegia microphylla55
ManyaraBrachystegia microphylla45
ManyaraLannea schimperi32
ManyaraJulbernadia globiflora77
ManyaraStrychnos spinosa34
ManyaraBrachystegia microphylla55
ManyaraBrachystegia microphylla51
ManyaraJulbernadia globiflora30.9
ManyaraBrachystegia microphylla53
ManyaraGrewia sp.8
ManyaraJulbernadia globiflora46
ManyaraBrachystegia microphylla10
ManyaraBrachystegia microphylla37
ManyaraVernonia exserstiflora5
ManyaraJulbernadia globiflora36
ManyaraLannea schimperi32
ManyaraBrachystegia spiciformis29
ManyaraBrachystegia spiciformis9.5
ManyaraBrachystegia microphylla22.7
ManyaraBrucea antidysenterica34.5
ManyaraBrachystegia microphylla78
ManyaraRhus vulgaris2.4
ManyaraBrachystegia spiciformis40
ManyaraBrachystegia microphylla68
ManyaraJulbernadia globiflora80
ManyaraBrachystegia spiciformis32
ManyaraJulbernadia globiflora28
ManyaraJulbernadia globiflora54
ManyaraBrachystegia microphylla69
ManyaraBrachystegia microphylla67
ManyaraAlbizia antunesiana63
ManyaraJulbernadia globiflora40
ManyaraGrewia sp.31.5
ManyaraJulbernadia globiflora29
ManyaraBrachystegia microphylla26.3
ManyaraBrachystegia spiciformis10.9
ManyaraBrachystegia microphylla17
ManyaraJulbernadia globiflora52
ManyaraLannea schimperi25.7
ManyaraBrucea antidysenterica16
LindiPterocarpus angolensis72
LindiDichrostachys sp.18
LindiLonchocarpus laxiflorus11
LindiPterocarpus angolensis37
LindiDeinbollia borbonica18
LindiPterocarpus angolensis32
LindiPterocarpus angolensis41.7
LindiPterocarpus angolensis60
LindiPterocarpus angolensis28.1
LindiPterocarpus angolensis54
LindiPterocarpus angolensis29.8
LindiBrachystegia boehmii29.2
LindiJulbernadia globiflora46.5
LindiZanha africana27.5
LindiPterocarpus angolensis16.7
LindiBrachystegia spiciformis71
LindiFicus natalensis11.5
LindiSterculia appendiculata65.2
LindiBrachystegia boehmii29
LindiBrachystegia spiciformis82
LindiPseudolachnostylis maprouneifolia41
LindiVangueria acutiloba12
LindiLonchocarpus bussei6.2
LindiSterculia appendiculata33
LindiPseudolachnostylis maprouneifolia24.2
LindiBridelia scleroneura65
LindiBrachystegia spiciformis25
LindiCordyla africana54
LindiAfzelia quanzensis48
LindiDiplorhynchus condylocarpon29
LindiDeinbollia borbonica32
LindiFicus natalensis34.2
LindiLonchocarpus laxiflorus2.4
LindiDeinbollia borbonica1.9
LindiHugonia castaneifolia30.2
LindiJulbernadia globiflora54
LindiDeinbollia borbonica86
LindiAfzelia quanzensis39
LindiAnnona senegalensis18.2
LindiBrachystegia spiciformis114
LindiBrachystegia spiciformis98
LindiDeinbollia borbonica19
LindiCombretum apiculatum24.2
LindiDodonaea viscosa80
LindiLonchocarpus bussei10.2
LindiBrachystegia spiciformis43
LindiJulbernadia globiflora4.5
TaboraAcacia robusta41.5
TaboraTerminalia sericea34
TaboraAfzelia quanzensis31
TaboraCombretum zeyheri25
TaboraCombretum molle12.5
TaboraPericopsis angolensis92
TaboraIsoberlinia globiflora82.8
TaboraBridelia scleroneura51
TaboraSchrebera koiloneura22
TaboraCombretum zeyheri13
TaboraBrachystegia boehmii55.3
TaboraLannea amaniensis27
TaboraLannea schimperi28
TaboraUnknown29
TaboraBrachystegia boehmii70
TaboraPterocarpus angolensis32.7
TaboraBrachystegia boehmii51.7
TaboraAlbizia antunesiana45.9
TaboraBrachystegia spiciformis102
TaboraBrachystegia spiciformis89
TaboraDalbergia melanoxylon17.2
TaboraScrerocaria birrea36
TaboraPericopsis angolensis36.5
TaboraPsorospermum febrifugum3.8
TaboraAntiaris toxicaria5.5
TaboraTamarindus indica4.3
TaboraBrachystegia spiciformis77
TaboraLonchocarpus bussei17
TaboraCordyla africana1.8
TaboraGarcinia sp.19
TaboraDalbergia melanoxylon12
TaboraSwartzia madagascariensis46
TaboraErythrophleum africanum25
TaboraBridelia scleroneura55
TaboraKigelia africana39
TaboraCombretum zeyheri25
TaboraStrychnos pungens6.5
TaboraBrachystegia spiciformis60
TaboraXeroderris stuhlmannii35
TaboraDichrostachys glomerata34

Data Availability

The data regarding this research work will be readily available when requested.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully thank the Tanzania Forest Services Agency for giving them a permit to conduct destructive sampling in all the study sites.

References

  1. P. Frost, “The ecology of Miombo woodlands,” in The Miombo in Transition: Woodlands and Welfare in Africa, B. M. Campbell, Ed., pp. 11–57, Center for International Forestry Research, Bogor, Indonesia, 1996. View at: Google Scholar
  2. I. Backéus, B. Pettersson, L. Strömquist, and C. Ruffo, “Tree communities and structural dynamics in miombo (Brachystegia-Julbernardia) woodland, Tanzania,” Forest Ecology and Management, vol. 230, no. 1–3, pp. 171–178, 2006. View at: Publisher Site | Google Scholar
  3. W. A. Mugasha, O. M. Bollandsås, R. E. Malimbwi et al., “Allometric models for prediction of above- and belowground biomass of trees in the miombo woodlands of Tanzania,” Forest Ecology and Management, vol. 310, pp. 87–101, 2013. View at: Publisher Site | Google Scholar
  4. E. Mwakalukwa, E. Meilby H, and T. Treue, “Volume and aboveground biomass models for dry Miombo woodland in Tanzania,” International Journal of Forestry Research, vol. 11, 2014. View at: Google Scholar
  5. E. W. Mauya, W. A. Mugasha, E. Zahabu, O. M. Bollandsås, and T. Eid, “Models for estimation of tree volume in the miombo woodlands of Tanzania,” Southern Forests: Journal of Forest Science, vol. 10, pp. 1–11, 2014. View at: Publisher Site | Google Scholar
  6. D. J. Kachamba, T. Eid, and T. Gobakken, “Above- and belowground biomass models for trees in the miombo woodlands of Malawi,” Forests, vol. 7, no. 38, pp. 1–22, 2016. View at: Google Scholar
  7. J. M. Abdallah and G. C. Monela, “Overview of Miombo woodlands in Tanzania. MITMIOMBO–Management of indigenous tree species for ecosystem restoration and wood production in semi-arid Miombo woodlands in eastern Africa,” Working papers of the Finnish Forest Research Institute, vol. 50, pp. 9–23, 2007. View at: Google Scholar
  8. United Republic of Tanzania, National Forest Resources Monitoring and Assessment of Tanzania Mainland, (NAFORMA): Main Results, Ministry of Natural Resources & Tourism, Tanzania Forest Services Agency, The Government of Finland and Food and Agriculture Organization(FAO) of the United Nations, 2015.
  9. P. K. T. Munishi and T. H. Shear, “Carbon storage in afromontane rain forests of the eastern arc mountains of Tanzania: their net contribution to atmospheric carbon,” Journal of Tropical Forest Science, vol. 16, no. 1, pp. 78–93, 2004. View at: Google Scholar
  10. E. Zahabu, Sinks and Sources: A Strategy to Involve Forest Communities in Tanzania in Global Climate Policy [Thesis for Award of PhD Degree at University of Twente], 2008.
  11. N. D. Burgess, B. Bahane, T. Clairs et al., “Getting ready for REDD+ in Tanzania: a case study of progress and challenges,” Oryx, vol. 44, no. 3, pp. 339–351, 2010. View at: Publisher Site | Google Scholar
  12. P. A. Dewees, B. M. Campbell, Y. Katerere et al., “Managing the miombo woodlands of southern Africa: policies, incentives and options for the rural poor,” Journal of Natural Resources Policy Research, vol. 2, no. 1, pp. 57–73, 2010. View at: Publisher Site | Google Scholar
  13. Food Agriculture Organization, Global Forest Resources Asssessment, vol. 163, Food and Agriculture Organization of the United Nations, Rome, Italy, 2010.
  14. United Republic of Tanzania, Tanzania’s Forest Reference Emission Level Submission to the United Nations Framework Convention on Climate Change (FREL), Ministry of Natural Resources and Tourism, Dar es Salaam, Tanzania, 2017.
  15. J. S. Makero, R. E. Malimbwi, T. Eid, and E. Zahabu, “Models Predicting Above- and Belowground Biomass of Thicket and Associate Tree Species in Itigi Thicket Vegetation of Tanzania,” Agriculture, Forestry and Fisheries, vol. 5, no. 4, pp. 115–125, 2016. View at: Google Scholar
  16. R. B. Heiligmann, M. Golitz, and M. E. Dale, “Predicting board-foot tree volume from stump diameter for eight hardwood species in Ohio,” Journal of Science, vol. 84, no. 5, pp. 259–263, 1984. View at: Google Scholar
  17. J. J. Corral-Rivas, M. Barrio-Anta, O. A. Aguirre-Calderón, and U. Diéguez-Aranda, “Use of stump diameter to estimate diameter at breast height and tree volume for major pine species in El Salto, Durango (Mexico),” Forestry, vol. 80, no. 1, pp. 29–40, 2007. View at: Publisher Site | Google Scholar
  18. T. C. Sawe, P. K. T. Munishi, and S. M. Maliondo, “Miombo of Tanzania: Implications on conservation and carbon stocks,” International Journal of Biodiversity and Conservation, vol. 6, no. 3, pp. 230–237, 2014. View at: Google Scholar
  19. H. I. Aigbe, W. W. Modugu, and B. A. Oyebade, “Modeling volume from stump diameter of Terminalia ivorensis (a. chev) in sokponba forest reserve, edo state, Nigeria,” Journal of Agricultural and Biological Science, vol. 7, no. 3, pp. 146–151, 2012. View at: Google Scholar
  20. S. B. Shamaki and S. O. Akindele, “Volume estimation models from stump diameter for Teak (Tectona grandis Linn f) plantation in Nimbia forest reserve Nigeria,” Journal of Environmental Science and Water Resources, vol. 2, no. 3, pp. 89–94, 2013. View at: Google Scholar
  21. V. B. Carl, “Volume prediction from stump diameter and stump height of selected species in Louisiana, New Orleans, United States of department of agriculture,” Research Paper, no. SO-182, 1982. View at: Google Scholar
  22. S. A. O. Chamshama, A. G. Mugasha, and E. Zahabu, “Stand biomass and volume estimation for Miombo woodlands at Kitulangalo, Morogoro, Tanzania,” Southern African Forestry Journal, vol. 200, pp. 59–70, 2004. View at: Publisher Site | Google Scholar
  23. A. M. Masota, E. Zahabu, R. E. Malimbwi, O. M. Bollandsås, and T. H. Eid, “Volume models for single trees in tropical rainforests in Tanzania,” Journal of Energy and Natural Resources, vol. 3, no. 5, pp. 66–76, 2014. View at: Publisher Site | Google Scholar
  24. M. Henry, N. Picard, C. Trotta et al., “Estimating tree biomass of sub-Saharan African forests: A review of available allometric equations,” Silva Fennica, vol. 45, no. 3, pp. 477–569, 2011. View at: Google Scholar
  25. N. Picard, L. Saint-André, and M. Henry, “Manual for building tree volume and biomass allometric equations: from field measurement to prediction,” Food and Agriculture Organization of the United Nations, vol. 215, 2012. View at: Google Scholar
  26. B. R. Parresol, “Modeling multiplicative error variance - An example predicting tree diameter from stump dimensions in baldcypress,” Forest Science, vol. 39, no. 4, pp. 670–679, 1993. View at: Google Scholar
  27. F. Baty and M. L. Delignette-Muller, Tools for Nonlinear Regression Diagnostics, 2018, http://cran.r-project.org/web/packages/nlstools/index.html.
  28. W. A. Mugasha, T. Eid, O. M. Bollandsås, and L. Mbwambo, “Modelling diameter growth, mortality and recruitment of trees in miombo woodlands of Tanzania,” Southern Forests, vol. 79, no. 1, pp. 51–64, 2017. View at: Publisher Site | Google Scholar
  29. W. A. Mugasha, E. E. Mwakalukwa, E. Luoga et al., “Allometric models for estimating tree volume and aboveground biomass in lowland forests of Tanzania,” International Journal of Forestry Research, vol. 2016, Article ID 8076271, 13 pages, 2016. View at: Google Scholar
  30. A. Ledo, K. I. Paul, D. F. R. P. Burslem et al., “Tree size and climatic water deficit control root to shoot ratio in individual trees globally,” New Phytologist, vol. 217, no. 1, pp. 8–11, 2018. View at: Publisher Site | Google Scholar

Copyright © 2019 Bernardol John Manyanda et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views9041
Downloads632
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.