Allometric models for estimating aboveground biomass and carbon stock for Diospyros mespiliformis in West Africa

Accurate estimates of aboveground biomass (AGB) strongly depend on the suitability and precision of allometric models. Diospyros mespiliformis Hochst. ex A. DC. is a key component of most sub-Sahara agroforestry systems and, one of the most economically important trees in Africa. Despite its importance, very few scientific information exists regarding its biomass and carbon storage potential. In this study direct method was used to develop site-specific biomass models for D. mespiliformis tree components in Burkina Faso. Allometric models were developed for stem, branch and leaf biomass using data from 39 tree harvested in Sudanian savannas of Burkina Faso. Diameter at breast height (DBH), tree height, crown diameter (CD) and basal diameter (D20) were regressed on biomass component using non-linear models with DBH alone, and DBH in combination with height and/or CD as predictor variables. Carbon content was estimated for each tree component using the ash method. Allometric models differed between the experimental sites, except for branch biomass models. Site-specific models developed in this study exhibited good model fit and performance, with explained variance of 81–98%. Using models developed from other areas would have underestimated or overestimated biomass by between –72% and +98%. Carbon content in aboveground components of D. mespiliformis in Tiogo, Boulon and Tapoa-Boopo was 55.40% ± 1.50, 55.52% ± 1.06 and 55.63% ± 1.00, respectively, and did not vary significantly (P-value = 0.909). Site-specific models developed in this study are useful tool for estimating carbon stocks and can be used to accurately estimate tree components biomass in vegetation growing under similar conditions.


Introduction
Biomass estimates remain essential for understanding the role of forests as carbon sinks and sustainable forest management (IPCC 2006;Pan et al. 2011). Tropical forests play an important role for the climate change as a main carbon sink for the whole world. Biomass and carbon accounting in natural and semi-natural vegetation types has rapidly become the focus over the last decades (Brown 2002;Chave et al. 2005;Chauhan et al. 2010Chauhan et al. , 2012, as being paramount step towards the implementation of the emerging carbon credit market mechanism such as Reducing Emission from Deforestation and Forest Degradation (REDD+) (Bernard et al. 2014). Aboveground biomass (AGB) is a useful measure for comparing structural and functional attributes of forest and savanna ecosystems across a wide range of environmental conditions (Brown et al. 1999). AGB also serves as an indicator for the distribution and abundance of vegetation above the ground.
Tree biomass can be estimated by: (1) destructive methods (Balima et al. 2019;Dimobe et al. 2018a,b;Bayen et al. 2015), (2) Non-destructive methods (Dimobe et al. 2018c(Dimobe et al. ,d, 2019, and (3) remote sensing techniques (Qureshi et al. 2012;Galidaki et al. 2017). Destructive methods which involve harvesting of trees, drying and weighing different components (Brown 1997) are time consuming and costly. They are also not applicable to large areas and are restricted in protected areas or for species that are endangered compared to non-destructive methods, while remote sensing is limited by access to technology and cloud cover. Allometric biomass models are statistical formulas relating tree biometric variables (e.g. diameter at breast height, tree height, basal diameter or crown diameter) to tree dry biomass. These models can be developed by destructive sampling (Kuyah et al. 2013;Dimobe et al. 2018a,b), semi-destructive (Mensah et al. 2016(Mensah et al. , 2017 and non-destructive (MacFarlane et al. 2014;Dimobe et al. 2019), and used based on predictors such as diameter at breast height (DBH), tree height, wood density, crown diameter or basal diameter, or the combination of these variables. Although allometric models are developed by destructive sampling, once developed, they allow non-destructive estimation of biomass over a large area at relatively low cost. Allometric models for estimating AGB can be linear, logarithmic, exponential or power law (Návar 2009;Bayen et al. 2015;Mensah et al. 2016). Power law however remain the most common used model because it is supported by growth that assumes a constant scaling rate across ontogenies (Mensah et al. 2017).
A limited number of studies have developed allometric models for species in dryland environments. These are based on data derived from wider geographical range; for example Brown (1997) developed allometric models from field measurements in savanna ecosystems of South America and dry forests of India; Chave et al. (2014) reported allometric models for pantropical forests. Recent studies (Shirima et al. 2011;Kuyah et al. 2012;Chave et al. 2014;Mensah et al. 2017) have also reported allometric models specific to African ecosystems. However, few of them refer directly to Sudanian savanna ecosystems of West Africa (Sawadogo et al. 2010;Bayen et al. 2015;Dimobe et al. 2018a,b). Accurate estimates of biomass in many tropical dryland regions are still lacking because of lack of allometric models. Challenges in obtaining measured field data in these species-rich ecosystems have made it difficult to develop models for predicting biomass of individual trees (Dimobe et al. 2018a,b). This is a major constraint given that accurate estimates of biomass depend on the appropriateness of the allometric models that are used (Sileshi 2014;Picard et al. 2015). There is mounting evidence that application of models not suitable to target vegetation results in large systematic deviations from observed data.
Diospyros mespiliformis Hochst. ex A. DC. is a species of the family Ebenaceae providing important ecosystem services to local population in sub-Saharan Africa. The tree serves as wind breaks, controls soil erosion and provides a habitat for other organisms. It also has potential for climate mitigation through carbon sequestration -a function that has not been quantified. The role of D. mespiliformis in climate change mitigation is therefore poorly understood. Carbon storage potential of the species can be underestimated or overestimated depending on the allometric model chosen for biomass estimation. To fill this gap, we developed allometric model for prediction of aboveground biomass and the biomass of tree components for D. mespiliformis in Burkina Faso. The specific objectives were to (i) develop species-specific allometric models for estimating total biomass and biomass of tree components (stem, branch, leaf) using basal diameter, DBH, total height and crown diameter, (ii) compare the performance of the selected model to predict the AGB with previously published and well-established global models, and (iii) predict the carbon content in leaf, branch, stem components of D. mespiliformis trees.

Study area
The study was conducted in three protected areas located in the Eastern and Southwestern Regions of Burkina Faso, namely, Boulon, Tapoa-Boopo and Tiogo (Fig. 1). Description of the biophysical contexts of the three forests areas are presented in Table 1. Boulon and Tiogo forests were classified by the colonial administration in 1955 and 1940, while that of Tapoa-Boopo was classified in the 80s (FAO 2004). In terms of climate, the forest of Boulon is located in the Sudanian zone while Tapoa-Boopo and Tiogo are located in the Sudano-Sahelian zone (Guinko 1984). All three locations experience a single rain season from May to October with annual precipitation ranging between 600-1000 mm; the rest of the year is dry. The mean annual temperature across the three sites is 35 °C; temperatures above 35 °C occur during peak period in the month of April (Bognounou 2009). Vegetation is dominated by shrub savannas, tree savannas and woodlands.

Description of the species
Diospyros mespiliformis is a mesopharenophyte that belongs to the family of Ebenaceae with a maximum height of 25 m (Arbonnier 2002). The species is widely distributed in sub-Saharan Africa, occurring from Senegal east to Eritrea, Ethiopia and Kenya, and south to Namibia, northern South Africa and Swaziland, but it is nearly absent in the more humid forest zones of West and Central Africa. It is also found in Yemen. D. mespiliformis is widely distributed in agroforestry parklands, woodlands, savannas and along riverbanks. Also known as African ebony or jackalberry, D. mespiliformis is a multipurpose species. For instance, the wood of D. mespiliformis is used for posts in construction, making musical instruments (e.g. drums) and household utensils (e.g. cups, spoons, pestles and mortars). The fruit are directly consumed for their vitamins and energy contents. The seeds are also eaten. The leaves are occasionally eaten as vegetable and serve as fodder for livestock. D. mespiliformis is commonly managed as a crop, but also grows in the wild and parklands across the sub-Saharan African savannas. It contributes ecosystem service such as wind breaks, controls soil erosion, and climate regulation.

Forest inventory and biomass sampling
Data were collected through two phases of sampling. For the first phase, forest inventory was conducted using 180 plots (60 plots for each site) of 1000 m² (50 m × 20 m) to collect data on DBH, total tree height, crown diameter and basal diameter on all individual trees of D. mespiliformis with DBH ≥ 5 cm inside each plot. The plots were established based on a stratified random sampling design. For the second phase, a destructive sampling of the species was carried out for quantification of tree biomass in 39 plots (13 plots for each site). Thirteen (13) individual trees of the species (one felled tree per selected plot) were harvested at each of the three sites (Boulon, Tapoa-Boopo and Tiogo), giving a total of 39 trees for measurements of biomass and subsequent development of allometric models (Table 2). To facilitate sampling, the individuals in each site were grouped into seven stem diameter classes: 10-20 (4 individual trees in Boulon, 6 in Tapoa-Boopo and 5 in Tiogo), 20-30 (5; 1 and 3), 30-40 (2; 3 and 1), 40-50 (1; 1 and 1), 50-60 (1; 0 and 2), 60-70 (0; 1 and 0) and 70-80 (0; 1 and 1). Individual trees were selected during the rainy season in June 2018 when all vegetation was in leaf. Prior to the harvesting, biometric variables, including  (Guinko 1984 DBH, total height, crown diameter and basal diameter were measured for all individual trees of D. mespiliformis. Crown diameter, DBH and basal diameter were measured using a diameter tape graduated in cm. Total tree height was measured from the base of the trunk to the tip of the tree using a clinometer. To avoid measurement error and to account for structural variety, the crown diameter was measured twice, along the east-west direction and the north-south one (Meyer et al. 2014). The diameter of all ramification was measured for tree forking below 1.30 m and quadratic mean diameter was computed. Selected trees were cut at ground level using a chain saw. Felled trees were separated into stem, branches and leaves. Branches were removed from the stem and cut into weighable sections. Leaves were separated from branches and collected into bags for weighing. Fresh weights of stem, branches and leaves for each individual tree were separately determined in the field using a 100 kg scale balance. Subsamples of leaves, stem and branches were taken, and their fresh weight determined using a 5 kg electronic balance (precision 0.5 g). The subsamples of branches and stem were taken as discs of 5 cm thickness depending on the diameter of the trunk or branch. The fresh weight of each disc was recorded in the field immediately after harvesting the trees. Two hundred grams (200 g) of leaf subsample for each tree was collected into bags to facilitate weighing.
The discs and leaf samples were transported to the laboratory where they were oven dried to a constant weight at 105 °C (samples of branches and stems) or 75 °C (samples of leaves). The dry weight of the samples was recorded immediately after removal from the oven. The sample dry weight was divided by corresponding fresh weight to obtain the dry-to-green weight ratio for each tree component. Component dry weight (biomass) was computed by multiplying subsample dry-to-fresh weight ratio by the fresh weight of the respective tree component. Total aboveground biomass of the tree was calculated as the sum of the biomass of all components in kilograms.

Estimation of carbon content and carbon stocks
Carbon content in selected individuals of Diospyros mespiliformis was estimated using ash method (Allen et al. 1986;Jana et al. 2009; Chavan and Rasal 2011) as follows. Composite samples were obtained from the dry matter samples of leaves, branches and stems for determination of their total carbon content. The samples were crushed in a cutting mill. Five 2-g samples of each component were then collected and submitted for analysis to the Laboratory of Plants and Soils at the University Ouaga1 Joseph Ki-Zerbo (Burkina Faso). Each 2-g sample was placed in a lidless porcelain crucible and placed for 2 h inside a muffle furnace set at 550 °C, until calcination was completed. The samples were then removed and cooled in a desiccator to be weighed later. After cooling, the crucible with ash was weighed and the percentage of carbon was calculated according to the formula by Allen et al. (1986): where W 1 is weight of crucibles, W 2 is weight of oven-dried grind samples with crucibles, and W 3 the weight of ash with crucibles; 0.58 is the content of carbon in dry organic matter (Allen et al. 1986). The carbon stock in leaves, branches and stems of each individual tree was calculated separately and summed up to determine the total organic carbon stock for each individual tree.

Data analysis
Biomass data were screened for outliers in scatter plots that assisted also in visually assessing the relationships between dependent variables (aboveground biomass and the biomass of components) and independent variables (basal diameter, crown diameter, DBH and total tree height). We examined the relationship between dependent variables and independent variables using linear, logarithmic, exponential and power models. Scatter plots were set to identify the theoretical distribution that fitted well with the data. As a result, the power law model (Eq. 3) outperformed the linear, logarithmic and exponential models; and therefore, was used to fit the biomass data. We, therefore, developed site-specific allometric models for each of the tree biomass components (leaf, branch, and stem) using the following linearized form of the power-law model (Eq. 4).
where Y denotes the leaf biomass, branch biomass or stem biomass; X stands for D 20 , DBH or CD; β 0 and β 1 the parameters, and ε is the random error. This model was logarithmically transformed into its linear equivalent form defined as followed: where ε′ is the additive error. Based on the selected aboveground biomass model, we developed allometric models for stem, branch, leaves and aboveground components based on (i) DBH or basal diameter as lone predictor variable, (ii) DBH 2 × H, and (iii) DBH 2 × H and CD, as independent variables for equation (5), (6) and (7), respectively: where Yi represents stem biomass, branch biomass and leaf biomass. We tested whether additional use of total height or/and crown diameter as predictor would improve the precision of the prediction. As separate use of potentially correlated independent variables in regression models could give rise to collinearity (Sileshi 2014), we only considered the allometric model that combined DBH (or basal diameter) with tree height and crown diameter as single predictor variable. DBH and height were used as combined predictor to account for variation of height for the same value of diameter while solving collinearity issue in linear regression (Dimobe et al. 2018a).
The use of natural logarithmic transformations helped to normalize the response variable, and thus to meet the prerequisites of linear regression for simple estimation of the parameters (Dimobe et al. 2018a). However, the use of natural logarithmic transformations, induces a systematic bias in the final estimation of the response variable, which can generally be corrected (predicted values were back-transformed into the original values) by applying the correction factor (Baskerville 1972): where RSE is the residual standard error of the regression.
Normality and homoscedasticity of the residuals were tested using the Shapiro-Wilks normality and Breusch-Pagan tests, respectively. The performance of the fitted models was compared using adjusted determination coefficient (R 2 ), residual standard error (RSE), Akaike information criterion (AIC), a likelihood criterion that penalizes the number of parameters (Burnham and Anderson, 2002) and root mean squared error (RMSE), calculated as: where AGB obs is the observed AGB and AGB est is the estimated AGB.
The site-specific AGB models were compared to four published global models in predicting the individual AGB. Published models selected were: 1 00673 Chave et al. (2014) for pantropical forests, where AGB is the aboveground biomass per tree in kg per tree, DBH is the diameter at breast height, ρ is wood density in g cm -3 , and H is the total height of the tree.  (1989), using DBH as the only predictor.
Differences between the selected model developed in this study and published models were assessed graphically, and the models were also compared using adjusted R 2 , RMSE and mean prediction error. The average relative systematic error was calculated as follow: where AGB est and AGB obs are the estimated and observed aboveground biomass, respectively. The significance of the estimated parameters at α < 0.05 was considered. Allometric models with greater adjusted R 2 , and smaller RMSE and AIC were considered to be the best fitting.
All statistical analyses were performed using R statistical software version 3.4.2 (R Core Team 2018).

Allometric models for estimating stem, branch and leaf biomass
Allometric models fitted to stem, branch and leaf biomass of Diospyros mespiliformis showed AIC values between 7.35 and 34.79, adjusted R 2 between 64.93 and 96.70, and RMSE between 0.256 and 0.677 (Table 3). Candidate allometric models and their regression coefficients are presented in Table 3. The coefficients of log-transformed allometric biomass models were significant for all biomass components (p < 0.001, Table 3). In general, the fitted allometric models for the biomass components were more accurate for branch and stem biomass than for foliage biomass (Table 3). For instance, the fitted models explained more than 90% of the total variations, except for leaf. Stem biomass was more accurate predicted in Tiogo and Tapoa-Boopo when DBH and height were used Table 3. Allometric models developed for estimation of biomass of stem, branches and leaves of Diospyros mespiliformis in Burkina Faso. ln(β 0 ), β 1 and γ represent the intercept and regression coefficients of the models, and their respective standard errors. The indicators of performance for the models include the adjusted coefficient of determination (Adj.R 2 ), residual standard error (RSE), root mean squared error (RMSE), Akaike information criterion (AIC) and correction factor (CF).

Regression coefficients
Performance  (DBH 2 × H) as compound variable (Table 3) than in Boulon. In Boulon, stem biomass was more accurately predicted using basal diameter (D20) alone than in Tiogo and Tapoa-Boopo. Similarly, leaf biomass was more accurately predicted when using DBH and crown diameter (DBH 2 × CD) as single variable in Tiogo and Tapoa-Boopo while in Boulon addition of crown diameter to basal diameter as compound variable (D20 2 × CD) provided accurate estimation of leaf biomass (Table 3). However, in the three sites, branch biomass was more accurately predicted when DBH and crown diameter were used as compound variable (Table 3). The predictive ability and accuracy of these models are demonstrated by the homoscedastic trend and good coincidence to the y = x linear equation shown in the diagnostic graphs of residuals (and observed values) versus predicted values of stem biomass, branch biomass and leaf biomass (Fig. 2).

Allometric models for estimating aboveground biomass
For each study site, nine allometric models were developed for aboveground biomass (AGB). Metrics for evaluation the performance of candidate allometric models developed for estimating AGB of D. mespiliformis in Burkina Faso are presented in Table 4. The adjusted coefficient of determination for the selected models was greater, explaining between 92.26-98.92% of the variance observed in the biomass. The AIC and RMSE values ranged from -9.02 to 21.76 and 0.119 to 0.425, respectively. In the three sites, DBH alone was the most suitable predictor variable than basal diameter and crown diameter for AGB (Table 3). Compared with model based on DBH as single predictor, adding tree total height to DBH as compound variable (DBH 2 × H) did not significantly improve the statistical fits (Table 4). However, addition of DBH to crown diameter and height as compound variable (DBH 2 × H × CD) in Tiogo and Tapoa-Boopo reduced the RMSE, AIC and RSE and made further significant change except in Boulon where DBH provided greater adjusted R 2 , smaller RSE, smaller RMSE and smaller AIC (Table 3 and 4). Graphs of residuals (and observed values) versus predicted values of aboveground biomass for the selected models (Fig. 3) suggested homogeneity of residuals and y = x linear trend. These results indicate that the selected models were the most suitable and can be used to predict aboveground biomass of D. mespiliformis in the study area.

Performance of published models
Models developed in this study performed more accurately than published models for general purpose (Figs. 4a,4b,4c). The performance of models developed in this study was comparable to that of Brown (1997) in terms of RMSE. However, it had less mean prediction error and greater adjusted R 2 (Fig. 4a). The model by Brown (1989) had the less RMSE but a greater mean prediction error (+97.18%; +80.34% ; and +60.92% for Tapoa-Boopo; Tiogo and Boulon respectively), while the models developed in this study presented comparatively a greater adjusted R 2 and the lowest average systematic bias (+2.28%; +3.10% and -3.28% for Tiogo; Tapoa-Boopo and Boulon respectively). The comparison of model developed by Chave et al. (2014) with Brown's model showed that Chave et al.'s model fitted the data more accurately than the model of Brown (1989), which leads to an overestimation of observed AGB (Figs. 4a,4b,4c). The model proposed by Jucker et al. (2017)

Allometric biomass components models
Several allometric models for estimating tree components' biomass are used to describe the relationship between tree parameters and biomass Mensah et al. 2016;Dimobe et al. 2018b). Power law model was used in this study to predict tree components' biomass. This result is consistent with previous reports (Xiang et al. 2016;Dimobe et al. 2018a). The models tested to estimate the relationships between biomass and explanatory variables are of the loglog form, which have been shown to predict the biomass of woody species with great accuracy (Chaturvedi et al. 2013). Chave et al. (2005) pointed out that biological data are heteroskedastic and that it is necessary to log-transform the variables to comply with the conditions of normality and homoscedasticity. The explained variance of the fitted allometric models at stem, branch and leaf level, ranging from 81.41% to 98.92% and was smaller for leaf biomass models than for branch and stem biomass models, as also reported by Mensah et al. (2017) and Morote et al. (2012) in previous studies. These results indicate that foliage allometries are less responsive to tree size than are branch and stem allometries (Dimobe et al. 2018a). This is probably because leaves are more sensitive to light exposure than branches and stem (Antin et al. 2013). These results also showed that the predictors for biomass components were not the same in the three sites. DBH and height used as compound variable appears to be good predictor for stem biomass in Tapoa-Boopo and Tiogo (Sahelo-Sudanian zone) whereas in Boulon (Sudanian zone), the reliable predictor is basal diameter. This finding accords with that of Balima et al. (2019) who also found in Burkina Faso that basal diameter is a good predictor for the stem biomass of Afzelia Africana in the Sudanian zone whereas in Sahelo-Sudanian zone, it was basal diameter and height (D20 2 × h), indicating that stem biomass vary according to climatic gradient. The somewhat smaller adjusted R 2 values for component biomass in Sudanian zone compared to Sahelo-Sudanian zone seen in this study have also been reported by Balima et al (2019). The additional use of height and crown diameter in this study is expected to account for variation in biomass of trees having the same diameter. The inclusion of tree height as additional predictor in allometric models was studied by Picard et al. (2015) and Dimobe et al. (2018a). It helps accounting for variation in AGB among trees with same value of DBH (Picard et al. 2015), thus reducing the estimate errors (Chave et al. 2005). For instance, in this study, inclusion of tree height as additional predictor improved the models for stem biomass for both Tapoa-Boopo and Tiogo. Dimobe et al. (2018a) reported improved biomass models due to inclusion of tree height. Moreover, according to Mensah et al. (2017), the use of DBH 2 × H as compound variable helps avoiding collinearity issues while accounting for within species variation of height for a given value of DBH. Other authors recommend the use of wood density, crown diameter, and site as additional predictor variables to improve biomass prediction (Chave et al. 20005;Fayole et al. 2013;Ngomanda et al. 2014;Chave et al. 2014;). We found in this study that inclusion of crown diameter reduced the RMSE of branch and leaf biomass estimates and improved the goodness of fit for the three sites.

Aboveground biomass model
In this study, nine allometric models were tested for their relative performance in predicting aboveground biomass. These models were satisfactory for predicting the AGB with explained variance ranging from 87 to 95%. Among the candidate allometric models, only the one with DBH as single predictor in Boulon (Sudanian zone) and the compound variable (CD × DBH 2 × h) in Tapoa-Boopo and Tiogo (Sahelo-Sudanian zone) appeared to be the most suitable for estimating the total aboveground biomass of D. mespiliformis. This concurs with for instance, Balima et al (2019) who found that DBH is the reliable predictor of AGB of Afzelia africana in the Sudanian zone of Burkina Faso while in the Sahelo-Sudanian zone it is DBH 2 × CD. The difference in terms of predictors used in the two climatic zones might be due to rainfall impact on growth and morphology on individual trees of D. mespiliformis. In Boulon, the model with DBH as single predictor produced greater adjusted coefficient of determination (R 2 = 97.11%) compared to the other predictors. This finding is in line with previous study that showed that DBH is the most commonly used independent variable (Henry et al. 2011;Mbow et al. 2013;Vahedi et al. 2014) since it is the easiest tree variable to measure. However, according to other authors such as Alvarez et al. (2012), models based on DBH only may underestimate the aboveground biomass and may show uncertainty. In Tapoa-Boopo and Tiogo, it was expected that the addition of tree height and crown diameter to DBH will improve the estimates of the allometric models. To account for that, the inclusion of tree height, crown diameter and DBH (CD × DBH 2 × h) in the allometric models as a compound predictor reduced the RMSE, RSE and AIC, and produced greater adjusted R² in Tapoa-Boopo and Tiogo. Although the differences were small, the selected biomass models with CD × DBH 2 × h as compound variable performed slightly more accurately in terms of adjusted R 2 , RSE, RMSE and AIC than the model with DBH alone.

Comparison with previously published models
We compared the models developed in this study with four previously published models developed by Brown et al. (1989), Chave et al. (2005), Chave et al. (2014), andJucker et al. (2017). We choose these models for comparison with our own model because they are commonly used for biomass estimation in Burkina Faso by previous studies (Qasim et al. 2016;Dimobe et al. 2018d;Dimobe et al. 2019). The general lack of agreement between the models developed for D. mespiliformis in this study and the generic models, can only be viewed as discouraging. Actually, the previously published models of Brown et al. (1989), Chave et al. (2005), Chave et al. (2014), and Jucker et al. (2017) are not wrong; they are simply not right for this species or the study sites. Ducey et al. (2009) made similar observation in the eastern Amazon and found errors from -33% to +29% that occurred when using off-site relationships, indicating that site-and species specific allometric models differ with species, tree status, climate and soil (Zianis and Mencuccini 2004), and are generally preferred over generic allometric models (Montagu et al. 2005). Kim et al. (2011), in their study, emphasize that the site-specific allometric models are more accurate in predicting the forest biomass estimates on the local level as they take into account the site effects. Among the previous published models used in this study, Brown's model greatly overestimated the observed AGB when applied to our dataset, while the Jucker et al.'s model underestimated it. This suggests that the models developed in this study are adequate for AGB estimation for D. mespiliformis in Burkina Faso. The previous published models were widely used because of the great goodness of fit to the sample data used to develop them (Nath et al. 2019). The limitation of these site-specific models inherent to the cost of biomass measurement and they are generally based on a small sample size.

Carbon stock
In this study, the predicted values of carbon content are greater compared to the one (0.5) which was recognized as an acceptable average to be used as conversion factor (Sarmiento et al. 2005;Redondo-Brenes 2007). Hence, the application of the recommended value would underestimate the value of carbon content in leaf biomass, branch biomass, stem biomass and AGB. Our findings are consistent to the observation in Central Panama (Elias and Potvin 2003), Costa Rica (Fonseca et al. 2012), and in Burkina Faso (Bayen et al. 2015;Dimobe et al. 2018a) where authors detected underestimations of aboveground carbon when the recommended value of carbon, 0.5, is used or applied. Furthermore, our results revealed that the average carbon content varied among the different parts of the tree. It has been reported that carbon content of tree components was dependent on the ash content which in turn, depends on the amount of structural components (Negi et al. 2003). The carbon content in stems of D. mespiliformis, was slightly greater than the amount in branches and leaves. This can be explained by the fact that woody tissues of trunk, roots, branches and twigs were greater carbon content pools than soft tissues of leaves, flowers and fine roots (Kraenzel et al. 2003). The results are also in conformity with findings by Fonseca et al. (2012), Wani and Qaisar (2014), Bayen et al. (2015) and Dimobe et al. (2018a).

Conclusion
To our knowledge, this is the first attempt to develop site-specific allometric models for Diospyros mespiliformis in Burkina Faso. Diameter at breast height alone is a good predictor of AGB in Boulon forest (Sudanian climatic zone), while in Tapoa-Boopo and Tiogo forests (Sahelo-Sudanian climatic zone), the inclusion of crown diameter and height are required to improve the AGB predictions. The models developed in this study can be used accurately to estimate the aboveground and tree components biomass in similar environments (climatic zones). The applicability of these models is therefore restricted to the climatic zones and range of data covered by trees used in this study. The results revealed also that D. mespiliformis contains greater amount of carbon than the reference value suggested by the IPCC, indicating a greater potential for carbon sequestration for mitigating adverse effects of climate change over large areas of West African savannas.