Stand Structure, Allometric Equations, Biomass and Carbon Sequestration Capacity of Acacia mangium Wild. (Mimosaceae) in Côte d’Ivoire ()
1. Introduction
Climate change has become a major concern for the world. Poor countries suffer the most of the adverse effects of global warming, given their vulnerability, lack or insufficient resilience and capacity to adapt (Makundi, 2014) . The land use change that represents between 10% and 20% of total carbon emissions is exacerbated by the conversion of tropical forests into cultivated and overexploited land, which considerably disrupts global carbon cycles (IPCC, 2007; Werf, 2009) . In Côte d’Ivoire, over 85% of the forest cover was lost between 1880 and 2008, under the combined effect of agricultural expansion and logging (SOFRECO, 2009) . This deforestation could continue further by 2030, if suitable measures are not taken.
Faced with these problems, the country has been committed since 2011 to the International Mechanism for Reducing Emissions from Deforestation and Forest Degradation (REDD+) implemented in the United Nations Framework Convention on Climate Change (UNFCCC). This mechanism should enable developing countries to benefit from financial compensation for their efforts of reducing deforestation, forest degradation and increase forest carbon stocks and forest conservation (UNFCCC, 2009) . Assessments of forest carbon stocks are therefore critical for countries planning to contribute on climate change mitigation impacts through their forest activities (Picard et al., 2012) . One of the strategic options is to use fast-growing species such as Acacias, Albizias, etc. to increase carbon stocks and also to meet the high demand of fuel for domestic energy.
Many studies have been undertaken since 1980 by Ivorian forest research on Australian Acacias and other fast-growing forest species (Dupuy & N’Guessan, 1991; Krisnawati et al., 2011) . These plants today offer agro-ecological and socio-economic potentialities through their ability to improve soil fertility, produce wood energy, timber and serve as windbreaks (Peltier & Ballé, 1993) . Acacia sp. plantations created by agroforestry systems or pure plantations can shoot one year after final tree cuttings (Gnahoua et al., 2014; Krisnawati et al., 2011; N’Guessan et al., 2013) through natural regeneration. In addition to bioenergy production, Acacias also play a major role in climate change mitigation through sequestration of carbon emitted in the atmosphere (Bakayoko, 2009; Bakayoko et al., 2012; Benbrahim et al., 2014; Doumbia, 2014; Krisnawati et al., 2011) . However, huge efforts are still needed to determine the best methodologies for CO2 assessment with reliable results.
Several methods are often used to determine aboveground biomass. Among these methods are biomass expansion factors and the destructive method for allometric equations. The biomass expansion factors (BAFs) are used to convert volumes from direct inventory data or direct field measurements (diameter, height) to total aboveground biomass (IPCC, 2006) . The most commonly used method for estimating carbon stocks is the destructive method and the development of allometric equations (Djomo et al., 2010; Fayolle et al., 2013; Ngomanda et al., 2014; Djomo & Chimi, 2017) . In Côte d’Ivoire from 73 allometric equations developed for 32 species, there is only one equation for Acacia mangium. Among these equations, fifty-four (54) allowed to evaluate only commercial volume, 15 trunk, stump and branches, and two for branches only (Henry, 2010) . None of these equations comply fully for estimation of biomass (Genet et al. 2011) , including trunk, branches, twigs, leaves and fruits. Henry et al. (2010) showed that improper selection of models for biomass estimation can lead up to 20% - 40% uncertainty on carbon stock estimates. In this context, the development of allometric equations for biomass estimation that takes into account all compartments including stump, trunk, branches, leaves, and fruits is becoming a national priority. Therefore, the aim of this study was to develop allometric equations that take into account all the compartments of the aboveground biomass. The specific objectives of the study were to:
- characterize three Acacia mangium stands of 3, 7 and 11 years old obtained from natural regeneration;
- develop allometric equations for estimating total aboveground biomass and also biomass of other tree compartments such as trunk, branches and leaves;
- evaluate the carbon sequestration potential of Acacia mangium at different ages.
2. Material and Methods
2.1. Study Site
The study was carried out in 2014-2015 at CNRA (National Center for Agronomic Research) forestry research station located within Anguédédou Forest (Figure 1). The station was created on March 25, 1965 and the management done by the former Centre Technique Forestier Tropical-Côte d’Ivoire (CTFT-CI) for forestry research. It has an area of 176 ha and includes many mono-specific stands or multi-species forest and agroforestry stands from former agricultural fallows with manioc, maize, etc. (Adou, 1986) . This station allows CNRA to carry out research activities on silviculture of local and exotic species, agroforestry, improve forest resources management, reforestation and study on energy properties of wood. Anguédédou Forest is located in a moist evergreen forest in the Southern part of Côte d’Ivoire, about twenty (20) kilometers on the Bouaké/Abidjan road (SODEFOR, 1996) . It lies between latitude north 5˚22' and 5˚26' and longitude west 4˚04' and 4˚13'. Its altitude varies from 40 m to 100 m asl (Adou, 1986) . This forest belongs to the rainforest area linked to sand and sandstone of the Continental Terminal characterized by the abundance of Thurraeanthus africana (Avodiré), Heisteria parvifolia (Amimimon), Xylopia acutiflora (Elo with small leaves), Monodora myristica (Adou, 1986; Guillaumet & Adjanohoun, 1971) .
Figure 1. Map of the location of the study area in Côte d’Ivoire.
The climate is hot and humid of subequatorial type with uneven rainy seasons (Eldin, 1971) , characterized by low temperatures of 24˚C - 34˚C, high humidity of 80% - 90%, and precipitation that reaches 1766 mm per year (Pega, 2008) . The ombrothermic diagram of the Abidjan region of the period between 2002 and 2012 shows 4 seasons: 1) a major rainy season from March to July, 2) a small rainy season from October to November, 3) a large dry season from December to February, and 4) a small dry season from August to September. The relief of the area consists of valleys and lowland plateaus varying from 40 to 60 m in the South and 70 or 100 m in the North. It belongs to the area which extends from the Fresco region in the west to Ghana in the east, over more than 440 km. The width of this plateau is about 25 to 30 km, except to the west of Bandama where it gradually decreases (SODEFOR, 1996) . The area is drained by Banco in east, Niangon (the main river of the research station) at the Center and Anguédédou in the west, which are three small rivers (Adou, 1986) . This area is located on a plateau belonging to the Lower Côte d’Ivoire sedimentary basin. This basin consists essentially of continental neogene sediments called tertiary or continental terminal sand, generally constituted by coarse detrital deposits, the most common of which is clayey and ferruginous sand (10 to 30 m thick) (Perraud, 1971) . Other facies such as clays, lateritic sandy clays, exceptionally non-clayey sand can also be encountered. Potentially fertile, these soils are subject to leaching and/or podzolisation (Adou, 1986) .
Acacia mangium is native from Australia (North Queensland), western Papua New Guinea and the provinces of Irian and Maluku in Indonesia (Benbrahim et al., 2014; Gnahoua et al., 2014; Dupuy & N’Guessan, 1991; N’Guessan & N’Goran, 1999; N’Guessan & Soro, 2006) . Its natural range is between 1˚ and 18˚ South latitude. Beyond its natural range, this tree can be found in Asian countries such as Cambodia, China, Japan, Philippines, Thailand and Vietnam. Acacia mangium plantations are estimated to be over 600,000 ha throughout the humid intertropical zone, especially in Southeast Asia (Gnahoua & Louppe, 2003) . Outside the Asian continent, Acacia mangium is used in reforestation in several African countries including Côte d’Ivoire, Senegal, Guinea, Benin, Cameroon, Democratic Republic of Congo, etc. In Côte d’Ivoire, experimental private plantations (PHCI, PALMINDUSTRIE, SODEFOR, etc.) of the species were created in Port Bouët, Oume, Dabou, La Mé, Assinie-France, Mopri, Yapo, San Pédro, Korhogo and Anguédédou (N’Guessan, 1991) .
2.2. Data Collection
2.2.1. Sampling Technique
The stands were subdivided into 3, 7 and 11 years old of Acacia mangium corresponding to the different stages of development. These stands were regenerated naturally one year after the final cut of the stands planted at 3 m × 3 m spacing after burning of the litter which stimulates spontaneous germination of the fallen seeds (N’Guessan et al., 2013) . A circular plot with a radius of 14 m (615 m2) was installed in each Acacia Magium stand. Within each plot, all trees were counted, and the height and diameter (circumference) taken at 1.30 m aboveground measured (Ekoungoulou et al., 2014) . To enable appropriate size distribution, tree sampling approach (CTFT, 1989) based on diametric or circumference structure enabled tree selection for destruction for allometric equations development. Hence, at least one tree was selected per diameter class and per stand to be felled. A total of 24 trees between of cbh between 31 and116 cm were felled and weighed entirely corresponding to 13 trees in the 3-year stand, 6 trees in the 7-year stand, and 5 trees in the 11-year stand. The inconsistency in number of trees in the three stands was due to the fact that it was laborious to felled and process bigger trees.
2.2.2. Biomass of Tree Compartments
After felling, each tree was separated into trunk, branches, leaves and fruits and the fresh biomass of each compartment weighed using a 500 kg scale. To obtain the dry weight of the trunk, three samples each of 4 cm thick were collected at the collar, at 1.30 m aboveground and at the top of the trunk (22 cm of circumference or 7 cm of diameter) of each tree (Girard, 1992) . Samples of leaves and branches (three samples of 4 cm thick) were also collected. In the laboratory samples of stems and branches were oven-dried at a constant temperature of 105˚C and leaves at 75˚C to a constant weight after 72 hours.
The amount of water (H%) in the various compartments (stem, branches and leaves) was determined after drying of the samples using the formula by
, with Wf the fresh weight of the sample and the Wo weight of the sample after drying. The dry weight of each compartment was determined using the relationship
with Wo the dry weight of the compartment under consideration, Wf the fresh weight of this compartment and H% the amount of water in the sample of the compartment (Girard, 1992) . The total dry weight of each tree was estimated by adding the dry weight of the various compartments of the tree.
3. Data analysis Structure of Stands
The stand structure showing the distribution and the variability of the number of stems per diameter class was used with the class interval of 4 cm. Other stand parameters such as density (number of trees per hectare), mean circumference, average height, mean annual increment, basal area were analyzed.
The mean circumference correlates strongly with stand volume. Moreover, for even stands, the distribution of trees by circumference class at breast height follows a normal distribution (Snowdon et al., 2002; Rondeux, 1993) . The formula
used for determining the average diameter (circumference) was
with Dm the mean diameter (circumference) of the considered stand in cm, fi the frequency of the class i and xi the circumference in cm of the tree at the center of the class i. The average height was obtained using the following formula (CTFT, 1989) :
with Hm the mean total stand height, Hi the total height of the tree i and N the number of trees sampled. The average annual increment was determined from the formula
with Cm the average diameter (circumference) of the stand and A the stand age. The basal area is the sum of the area per hectare, of the cross-sections of the trees taken at 1.30 m above ground, expressed in m2∙ha−1. When the basal area of the average tree was equal to the average stand area (Pardé & Bouchon, 1988) , the basal area G of the stand was estimated using the following relationship (CTFT, 1989)
where G is the basal area in m2∙ha−1, Cm the average circumference in m and d the stand density in N∙ha−1.
3.1. Development of Allometric Equations
Allometric equations are relations that link biomass with one or two independent variables such as diameter and height (Lotfi, 2008) . Most allometric equations use diameter at breast height taken at 1.30 m (dbh). The models most frequently found in the literature for biomass estimation are of two types: the power model and the polynomial model (Henry et al., 2011) . In this study, the power model was used because polynomial models frequently exhibit abnormal behavior outside their validity domain (e.g., negative predicted values). The most frequently used mathematical formula for adjusting a biomass is B = aDb (Picard et al., 2012) , where B is the total dry aboveground biomass, D the diameter at breast height (ddh) and a, b the regression coefficients. To take into account the heteroscedasticity of the data (Djomo et al., 2010; Djomo et al., 2016) , this formula is often modified using the logarithmic transformation through the relation
(Niklas, 1993; Djomo et al., 2010; Picard et al., 2012; Djomo et al., 2016) . The data from this study were adjusted by the following three models:
; (1)
; (2)
(3)
where AGB is the total or compartment aboveground biomass in kg, D the circumference taken at 1.30 m aboveground in cm and H the total height of the tree in m. The three models were repeated with the circumference to allow reporting regression equations with circumferences as input variables.
The first step in the analysis of the data was the graphical exploration which allowed visualizing the relationships between variables from a point cloud in order to get an idea of the type of model to be adjusted. Although the graphs of the non-transformed and log-transformed data looked similar, the log-transformed data had a correlation coefficient higher in all the scenarios tested and the log- transformation was therefore selected (Figure 2). The logarithmic transformation of the data generally leads to a bias in the estimation of biomass (Parresol, 1999; Chave et al., 2005) , which was corrected by multiplying the estimated biomass by a correction factor (CF) determined by the relation
(Djomo et al., 2016; Djomo & Chimi, 2017) .
Two criteria were used to select the best predictive model, the residual standard error (RSE) and the Akaike information criterion (AIC). RSE represents the standard deviation between the observed value and its prediction, estimated by the relation
where SSE is the sum of squared errors between the observed value and the prediction and
the number of degrees of freedom. The Akaike Information Criterion (AIC) is a measure of the quality
Figure 2. Representation of the scattered plot of the total aboveground biomass without log-transformation (left) and with log-transformation (right).
of the model used for the set of data considered. It allows to compare several models and to make the selection of the best model.
where p is the number of parameters in the model and L the maximum likelihood. These criteria make possible to judge the goodness of the model’s fit; the lower the criteria, the better the model (Chave et al., 2005; Saint-André et al., 2015) .
3.2. Determination of Carbon Stock
The total biomass per hectare Bh (t∙ha−1) was determined using the formula
with Bc the biomass in compartment c (trunk, branches, leaves), i the number of trees used and S the surface area considered. The carbon content in wood of 50% (IPCC, 2003) allowed estimating the carbon stock CS per hectare (tC∙ha−1) using the relation
. The annual carbon sequestration capacity C of aboveground biomass per hectare per year (tC∙ha−1∙year−1) was estimated using the relationship
in which A is the age of the stand. Since the CO2 equivalent is the measure of greenhouse gas emissions, the conversion of the carbon equivalent (eqC) into carbon dioxide equivalent (eqCO2) was done by multiplying C by 44/12, where 1 eqC = 3.66 eqCO2.
The biomass expansion factor (BEF) which is a factor used to convert trunk biomass to total biomass (Mavouroulou, 2012) was estimated in this study using the formula BEF = Total Biomass/Trunk Biomass.
The data of this study were compiled with Excel spreadsheet and the linear regressions analyzed with R software version 3.1.1 (2014-07-10).
4. Results
4.1. Stands Structures
Figure 3 shows the diameter distribution of the stands studied. For the stands of 3 to 11 years studied, this distribution was unimodal characteristic of regular forests. Within the 3-year-old stand there were four diameter classes ranging from 9 to 24 cm with the highest proportion of individuals (40% trees) in the 12 - 16
Figure 3. Diameter distribution of the studied stands.
cm class. Within the 7-year stand, there were five diameter classes ranging between 12 and 32 cm with the highest proportion of individuals (32% trees) in the 20 - 24 cm class. Within the 11-year-old stand, there were five diameter classes ranging between 20 and 40 cm with the highest proportion of individuals (41% of trees) in the 28 - 32 cm class.
In Table 1, the main characteristics of the stands showed that the stand structure varies greatly with age. The mean diameter was 12.5 cm for the 3-year stand, 23.4 cm for the 7-year stand and 29.1 cm for the 11-year stand. The average total height was 12.5 m for the 3-year stand, 19.0 m for the 7-year stand and 23.3 m for the 11-year stand. The average basal area was 10.42 m2∙ha−1 for the 3-year stand, 30.75 m2∙ha−1 for the 7-year-old stand, and 36.51 m2∙ha−1 for the 11-year stand. The highest density 845 stems∙ha−1 was obtained in the 3-year stand and the lowest 553 stems∙ha−1 in the 11-year stand. The mean annual increase on dbh and height was inversely proportional to the age of the stand. For dbh, it was 2.64 cm∙year−1 for 11-year stand, 3.34 cm∙year−1 for 7-year stand and 4.18 cm∙year−1 for 3-year stand. For height, it was 2.11 m∙year−1 for 11-year stand, 2.71 m∙year−1 for 7-year stand and 4.19 m∙year−1 for 3-year stand.
4.2. Distribution of Tree Biomass
Table 2 presents the distribution of the aboveground biomass by dbh and height for trunk, branches, leaves and total biomass. In the 3-year stand, the total tree biomass varied from 40 kg to 235 kg; the 7-year stand was relatively homogeneous varying from 190 to 327 kg. In the 11-year stand the total biomass varied from 395 to 568 kg. The largest biomass was found in the trunk followed by branches and leaves (Figure 4). The results of this study showed that the trunk stored on average 60%, 72% and 80% of the total biomass in 3-year, 7-year and 11-year stands respectively. The largest leaves was observed in the 3-year stand
Table 1. Stand characteristics by age group.
Table 2. Biomass (kg) measured for the 24 felled trees.
(19%) and the lowest in the 11-year stand (8%). The distribution of branches biomass in the different stand age classes was similar to that of leaves. It was more important in the 3-year stand (21%) and lower in the 11-year stand (12%).
4.3. Allometric Equations
For the estimation of the total biomass and also of the compartments trunk, branches and leaves, this study made it possible to develop allometric equations. The best equation for estimating total biomass was total
with Adj∙r2 = 0.97 and AIC = −89 which equivalent to
if circumference was used as input variable (Table 3). For trunk biomass, model II had the same value of Adj∙r2 = 0.98 with model III and the lowest AIC = −92. However, validation with
Figure 4. Biomass (%) in the various tree compartments.
Table 3. Results of regression analysis with models I. II. and III., a, b, c adjusted parameters of the model, D diameter at breast height, C circumference at breast height, H total height and B biomass. The best fit parameters are reported for each model. N: number of trees measured. CF: correction factor. AdjR2: Adjusted square regression coefficient. RSE: Residual standard error and AIC: Akaike information criterion.
*: selected equation.
original data showed that model III estimated the total and mean with an error of −3.81%, (model II: −4.22%, model I: −4.15%) and selected as the best fit. Therefore, the best equation to estimate the trunk biomass was
which corresponded to
if circumference was used as input variable. For branch biomass, the highest Adj∙r2 and lowest AIC was obtained with model III; however, validation with original data showed that equation II estimated the total and the mean with an error of −2.46% (model I: −3.02%; model III: −2.98%) and was selected. Therefore, the best equation for estimating branch biomass was
which corresponded to
if circumference was used as input variable. For leaf biomass, model III had the highest Adj∙r2 and lowest AIC, followed by model I. However, validation with original data showed that model I was the best equation with an error of −2.14% (model II: −2.78%; model III: −2.33%). Therefore, the best equation for estimating leaf biomass was
which corresponded to
if circumference was used as input variable. These equations were used to determine carbon stocks per hectare in the study stands.
4.4. Carbon Stock, Carbon Sequestration and BEF
The aboveground biomass stored in average tree increased with the age varying from 66 kg in 3-year stand to 241 kg in 7-year stand and 380 kg for 11-year stand. It was the same case for biomass stored in the leaves (10 kg in 3-year stand to 30 kg in 11-year stand), branches (15 kg in 3-year stand to 49 kg in 11-year stand) and trunk (40 kg in 3-year stand to 313 kg in 11-year stand) (Table 4). The corresponding total aboveground carbon stored was 28 tC∙ha−1 in 3-year stand; 86 tC∙ha−1 in 7-year stand and 105 tC∙ha−1 to 11 years. The equivalent carbon dioxide sequestration was higher for 7-year stand with 45 teqCO2 ha−1∙year−1 and lower in 3-year stand with 34 teqCO2∙ha−1∙year−1. The biomass expansion factor (BEF) decreased with stand age and was 1.66 in 3-year stand, 1.37 in 7-year stand and 1.21 in 11-year stand (Table 4).
5. Discussion
The diameter structure of the different stands showed a normal distribution of the stems with a maximum frequency of 40% in the class of 12 - 16 cm in 3-year stand, 32% in the class 20 - 24 cm in 7-year stand and 41% in class 28 - 32 cm in 11-year stand. This diameter structure is typical of even forest stands with high regeneration capacity (Rollet, 1974) . This diameter distribution structure was consistent with the study of Deb et al. (2014) that found normal distribution in a Shorea robusta plantation in Bangladesh. It was also consistent with the study of Ribeiro et al. (2014) which found the same structure in a plantation of Eremanthus erythropappus. The mean diameter of Acacia mangium 3-year, 7-year and 11-year stands of this study were 12.53 cm, 23.41 cm and 29.1 cm respectively.
Table 4. Biomass stock and carbon sequestration capacity.
The mean total height (m) varied between 12.58 m and 23.26 m. It was also the case with the stands basal area which grew with age from 10.42 m2∙ha−1 for 3-year stand to 36.51 m2∙ha−1 for 11-year stand. This growth in mean diameter, mean height and basal area with age was close to the results of other studies in forest plantations (Dupuy & N’Guessan, 1990; Lehtonen et al., 2004) .
Stand Densities ranged from 845, 715 and 553 stems∙ha−1 respectively for 3-year, 7-year and 11-year forest. These values were lower compared to the growth model prediction of Acacia mangium (Dupuy & N’Guessan, 1990) with stand density varying from 1350 stems∙ha−1 in 3-year stand to 970 stems∙ha−1 in 9 year stand. The difference could be due to the fact these values are results of model prediction and not direct measurements; also human activities such as tree harvest for fire wood, logging, bush fires and other factors such as natural disturbance (e.g. Chablis) and intraspecific competition (light, water and nutrient availability) can also explain this difference with the model (Oliver & Larson, 1996; Guillaume, 2006) .
The biomass stored in the trunk (up to 22 cm circumference or 7 cm diameter) represented 60%, 72% and 80% of the total biomass in 3-year, 7-year and 11-year stands respectively. For leaves, the highest biomass (19%) was stored in the Acacia mangium of 3-year stand. This proportion decreased with age to reach 8% in the 11-year stand. The proportions of biomass stored in the trunk were consistent with the study of Heriansyah et al. (2007) in Acacia mangium stands of 3 to 10 years in West Java in Indonesia which showed that the biomass stored in the trunk was between 63 and 71% of the total biomass. At the branch level, the highest proportion representing 21% of total biomass was observed in the Acacia mangium stand of 3 years. This proportion decreases with age, reaching 12% in 11-year stand. These results were close to those of the 7-year stand of Acacia mangium provenance trials in Ibi carbon sink project in Democratic Republic of Congo (Yusufu, 2014) . The study of Henry et al. (2010) in a dense moist forest of Ghana showed a significant difference with this study in leaf biomass that represented the smallest proportion with about 1% of the total aboveground biomass. In three native mixed species stand of a Hyrcanian forest in Iran, constituted by Fraxinus excelsior, Pterocarya fraxinifolia, Populus caspica, Quercus castaneifolia and Ficus carica, the crown biomass (branches and leaves) varied between 20%, 21% and 25% (Vahedi et al., 2016) which was close to this study ranging between 20%, 28% and 40% in 11-year, 7-year and 3-year stands respectively.
Data collected from 24 trees in 3 stands of Acacia mangium allowed developing allometric equations using two independent variables the diameter or circumference at 1.30 m aboveground, and the total height. To measure the goodness of each of the models the residual standard error (RSE) and Akaike Information Criteria (AIC) was assessed and the error made on the estimation of total and mean biomass used to validate the final selection. For total biomass, Model I selected with only diameter (circumference) measured perfectly the total and mean biomass and had also the lowest AIC and highest Adj∙r2. This showed that the total height did not influence significantly the total biomass as the total height is a function of D (Chave et al., 2005; Fayolle et al., 2013) . The best equation for estimating trunk biomass was obtained with model III (RSE = 0.138, AIC = −89, error = −2.98% on total or mean) and for branches model II (RSE = 0.314, AIC = −52, error = −2.46% on total or mean) showing that the total biomass was needed for precise estimate of these compartments. As with total biomass, the best fit for leaf biomass was model I (error −2.14%) which is not influenced by the total height.
The aboveground biomass of carbon stored was 27.8 tC∙ha−1 in 3-year, 86.3 tC∙ha−1 in 7-year and 105.0 tC∙ha−1 in 11-year stands respectively. These values were close to those obtained by Yusufu (2014) in the 46 provenances trial of Acacia mangium in the Democratic Republic of Congo. Bakayoko et al. (2007) study in an Australian Acacia experimental plantation of 14 and 15 year old in Côte d’Ivoire found aboveground biomass varying from 35.7 tC∙ha−1 to 121.3 tC∙ha-1 which was close to the results of this study.
The BEF of this study converted the trunk biomass to the total biomass (trunk, branches and leaves). In case the volume of the trunk is the data available, the wood density is needed to convert this volume data to trunk biomass (biomass = volume × wood density) before applying the BEF to obtain the total biomass. The BEF was 1.66 in 3-year, 1.37 in 7-year and 1.21 in 11-year stands respectively. Henry et al. (2010) found BEF ranging between 1.13 and 2.20 while Mavouroulou (2012) obtained BEF ranging between 1.04 and 3.88. For tropical forests, IPCC (2003) suggests using locally developed BEF for reliable estimates or 3.4 BEF in case there is no value in the study area.
6. Conclusion
This study allowed determining the diameter structure of Acacia mangium stands of 3-year, 7-year and 11-year stands respectively obtained from natural regeneration. Destructive data from 24 trees allowed developing allometric equations for the estimation of total aboveground biomass and the biomass of trunk, branches and leaves. Three models were tested and the residual standard error (RSE), Adj∙r2 and AIC were used to determine the goodness of fit and the error of total or mean biomass used to confirm and make the final selection of the best allometric equation. The total height did not influence the estimate of total aboveground biomass and the biomass of leaves. However for precise estimate of leaves and branches, total height needed to be added in the model and independent variable. The allometric developed were used to estimate carbon stock, the CO2eq and the carbon sequestration capacity of the three forest stand of Acacia magium studied. For each stand, the BEF was also developed to convert the trunk data available for example from forest inventory data into total aboveground biomass.
The results of this study contribute to improve the estimation of the aboveground biomass of even stands resulting from natural regeneration. However, this research did not include the root biomass and recommends further studies providing precise information on this important compartment for the estimation of carbon stocks.
Acknowledgements
The authors would like to address their heartfelt thanks to CNRA who provided the financial support for this study and most especially to Dr. DEA Goué Bernard, Regional Director of CNRA of Asbidjan who made the financial arrangement. The authors extend their thanks to Dr YTE Wongbé, General Director of CNRA and in particular Dr ZAKRA Nicodème, Director of Scientific Research and Development Support who provided the training facilities needed for this research.