Tropical forests sheltered most of the world’s biodiversity, accounting for 2/3 of the 1.3 million recorded animal and plant species  . They contain 40% - 50% of terrestrial carbon and play a major role in the global carbon cycle  . At the beginning of the 19th century, the world’s tropical forests covered around 16 million km2. Now, they cover only around 3.8 billion hectares  . Since the 1992 Earth Summit, it has been established that biodiversity is seriously threatened by humans  . The loss of forest cover resulting from deforestation and forest degradation results in the loss of more than 13 million hectares of forests worldwide each year  and damages biodiversity  . This activity contributes actively to about 10% - 15% of global annual greenhouse gas emissions  .
In Republic of Congo, forests cover 22.4 million hectares divided into three major groups: the North Congo massif with 17.4 million hectares, and the Mayombe and Chaillu massifs in southern Congo with around 4.7 million hectares. These last two were the first exploited since colonial times by the reason of the Pointe-Noire carriage proximity and the presence of the CFCO (Railways Congo Ocean Company) railways and the Ogooué Mining Company  . This led foresters to ask themselves very early on the availability of biodiversity and the restoration of forest heritage. The creation in 1934 of the Mbuku Nsitu arboretum in the Mayombe Massif, has historically testified to the interest that some forestry services or researchers have made early in the question concerning the sustainability of Congolese forests  . Thus, an initiative of reforestation of the degraded zones of Mayombe began in the season 1949-1950 in order to contribute to the maintenance and the restoration of the biological diversity and productive functions of the southern Congolese forests. After independence, the Congolese government supported in 1970 through cooperation with UNESCO through the Intergovernmental Program on Man and the Biosphere (MAB) and the United Nations Development Program launched a vast development project for the Mayombe Massif  .
Today, the Mayombe afforestation and reforestation cover an area of around 10,000 ha  , including the reforestation of Terminalia superba Engl. & Diels from Bilala. Covering an area of 421 hectares and set up in 1954 in open forest, these reforestations were originally intended to offset the lack of timber. However, new opportunities have emerged over time; the evolution of Bilala’s reforestation has resulted in the rebuilding of the natural forest, which has led to the resettlement of biodiversity and has attracted the attention of many scientists and local populations. In a context of mitigation to the effects of climate change  , a new potential valorization of Bilala’s reforestation could be linked to its carbon sink function. Indeed, under the Kyoto Protocol ratified in December 1997, carbon credits can be granted to better conserve its forests   . In accordance with the REDD+ (Reducing emissions from deforestation, forest degradation, and forest conservation, sustainable management of forest, and enhancement of forest carbon stocks) process set up by the Clean Development Mechanism of the Kyoto Protocol, the signing of a Carbon credits with the World Bank, through the Carbon Fund set up since 2011, could help the National Reforestation Service (SNR) to valorize forest carbon and boost the management of its plantations  . In order to benefit from the carbon-credit, it must be possible to demonstrate that reforestation contributes to increased carbon sequestration and to quantify this increase.
Studies on floristic diversity and carbon stock estimation  in southern Republic of Congo remains fragmented and mostly located in protected areas (Dimonika and Conkouati-Douli Reserve). They focus on general knowledge of flora  . They have sometimes focused on undergrowth vegetation in forest plantations  , while others have focused on estimating the carbon stock   . A study on the evolution of tired out after Cassava plantations  . As the issues are different, the results are difficult to compare  . There is existence of a few data on reforestation of Terminalia superba Engl. & Diels in Republic of Congo. Those that are available focus on cuttings and spacing trials. No studies on undergrowth vegetation of Terminalia superba Engl. & Diels reforestation has not been conducted yet. As forestry research is concerned, this current study on biodiversity, carbon stock assessment and forest management of the Bilala reforestation is the first to contribute to the knowledge of Mayombe flora in its diversity. In the context of the fight against climate change and preservation of biodiversity, How can reforestation contribute to the implementation of the REDD + process in Bilala?
The general objective of this study was to assess the biodiversity of Terminalia superba Engl. & Diels plantations, and its undergrowth and also to quantify carbon stocks about the above-and below-ground biomass of trees in order to assess the impact of reforestation on forest regeneration and the generation of carbon stocks in order to provide the SNR (National Reforestation Service) with a tool to assist the management of plantations and local populations. The specific objectives of this study was (I) to identify, by floristic surveys, trees with DBH ≥ 5 cm at 1.30 m from the soil, as well as the herbaceous present; (II) Quantify carbon stocks sequestered by trees with DBH ≥ 5 cm; and (III) identify bottlenecks that hinder the successful management of Terminalia superba Engl. & Diels reforestation to propose a long-term management strategy.
2. Materials and Methods
2.1. Study Area
The study sites were located at Bilala (4˚30'S, 12˚13'E) in southeastern Republic of Congo, in Kouilou Department (Mayombe), close to Mvouti District with an altitude of 30m. Bilala is bounded at the North by Bilinga, at the South by the Paris forest and the Congo-Angola international border, at the East by Tchyvala village and at the West by the coastal plain and Yanga village  . Thus, the research work was carried out in Bilala and its surroundings in three plantation sites under the supervision of the Mayombe station. Specifically, the three sites studied are: 1) the Matete1 reforestation site; 2) the Bilala site installed for the protection of the “Gutter” watershed against erosion; 3) the coastal plain site.
Bilala is subject to a warm and humid climate combining equatorial, tropical and oceanic traits. According to the  study, the climatic characteristics of Bilala are as follows: The climate is low-Congolese with a variant of the “Guinean” type; the average annual rainfall is between 1200 and 1500 mm  . The hottest month is March with an annual temperature of 27.5˚C. A rainy season from mid-October to mid-May with the highest temperature and illumination maxima of the year; it is interrupted by a small dry period which is not recognized as a dry season. A dry season from mid-May to mid-October with a lowest temperature and illumination minimums of the year. Relative humidity is always around 85%; relative humidity is greater during the dry season than during the rainy season. The sky is often permanently covered with fog and clouds (Figure 1 and Figure 2). Bilala is flat and at the same time rugged. The rocks are of mainly schistose or sandstone origin, with intrusions of various crystalline rocks: granites, quartz diorites, green rocks  . The soils of Bilala are mostly ferralitic, poor in exchangeable bases, highly saturated and very acidic  . The Bilala river system is organized around the Loémé River. The village is contorted by several swamps  .
2.2. Data Collection
The floristic surveys were carried out on the one hand on the monospecific stands planted with Terminalia superba Engl. & Diels aged 64, 31 and 20, and at
Figure 1. Study area location.
Figure 2. Climograph of Bilala provided by the meteorological station around the study area from January to December 2018.
the same time in the undergrowth of its stands. All data were collected during the period from July 16, 2018 to November 1, 2018. In the first place, the surveys were carried out on the block of Terminalia superba Engl. & Diels planted in Matete in 1954 with an area of 421 ha and then in two stands of 31 years and 20 years of respective areas of 0.3 ha and 1 ha. Thus, in the undergrowth of its three blocks, woody, herbaceous and lianascent species served as identification and measurement support (case of woody plants).
Data collection was carried out in nine rectangular sub-plots of 20 × 25 m each, installed in three blocks of Terminalia superba Engl. & Diels still called Limba and in a fallow field, a total area of 0.5 ha. The shape of the plot was inspired by the work of  where the estimation and monitoring of carbon stock could be done in a temporary plot. Each sub-plot was divided into two strips of 10 m × 25 m each. The floristic surveys consisted of recording each vegetative species present in the sampling area band by band and at the same time measuring all the planted or spontaneous wood specimens with a DBH ≥ 5 cm at 1.30 m from the soil. For herbaceous species, counting of individuals was not done in full. However, band-by-band notes ranging from 25 to 100 individuals representing the numerical abundance of the species within the band were assigned to each herbaceous species inventoried according to the Ushuaia multiple inventory scale which is as follows:
25: slightly abundant; 50: moderately abundant; 75: abundant; 100: very abundant.
The plots were dispersed in the three blocks studied. This operation began with a careful observation of the place. The planting line was then found using the compass that gave an East-West direction. After that, the first base milestone (milestone 1) was placed next to the first tree of the planting line. Using the double decameter of 20 m, 20 m were drawn starting from milestone 1 following the South-North direction and 25 m following the planting line using the double decameter of 30 m. Four basic milestones were marked at the ends of the sub-plot. Finally, the sub-plot was divided into two equal rectangular strips (10 m × 25 m) where two milestones were placed at the ends of each strip. At the end of the operation, six milestones were marked with red paint and a GPS coordinate was recorded on the first tree of the sub-plot.
Concerning woody trees, the circumference was first measured using a tape measure at 1.3 m from the ground and the diameter calculated on site. The diameter of the trees was measured above the buttresses or stilt roots. The height was then measured using the vertex and the species was identified. The height of considered tree was which of the highest branch. To finalized, a number was assigned in ascending order and marked with red paint to each tree measured and identified. Species identification was made with the naked eye using the following botanical characteristics: bark shape, latex color, leaf arrangement, crown shape, trunk base shape, and its termination. A manual has been used as a basis for the identification of species: Botanical Manual and Inventory of Biodiversity of the Association of International Tropical Timber Technicians  . Regarding the case of Ricinodendron heudelotii (Baill.) Pierre ex Pax which loss its leaves, the bark allowed the identification of the species. Thus, all trees with DBH ≥ 5 cm were measured, identified and counted from band 1 to band 2 and at the same time all herbaceous and lianascent species were identified. Tree species recorded in study area have been identified using earned scientific knowledge during the biosystematics courses at university. Field work experiences during the university studies period also were a help to identify flora species of study area.
2.3. Data Analysis
The general checklist of species composing the flora procession has been established after digital processing of nine sample sub-plots. On the basis of The African plants database (version 3.4.0) of Conservatory and botanical garden of Geneva, Switzerland & South African National Biodiversity Institute, Pretoria (accessed 1 April 2019 at http://www.ville-ge.ch/musinfo/bd/cjb/africa/recherche.php), The Global plants database (Accessed on 7 April 2019 at http://plants.jstor.org, The working list of all plant species database (Retrieved April 1, 2019 from http://www.theplantlist.org), and The Xycol database (The list of scientific and vernacular woods names: accessed on April 7, 2019 at http://www.xycol.net/inde x.php?categorie=0&sess_langue=430), all flora species were checked and confirmed. Using the above databases and websites, all identified flora species were confirmed. A diversity index is a quantitative measure that reflects how many different species there are in a community, and simultaneously take into account how evenly the basic entities (such as individuals) are distributed among those types. Vernacular and local nomenclatures of the species recorded in the undergrowth of the stands of Terminalia superba Engl. & Diels were verified in the South-Congo Plant Database of the National Center for Forest Resources and Wildlife Inventories and Management  and tough database of “Global plants database” (Accessed April 1, 2019 at http://plants.jstor.org and http://plantes-botanique.org). Wood density for each species has been founded from the database “global wood density database” available at http://datadryad.org/handle/10255/dryad.235 (Retrieved 10 October 2018).
2.3.1. Analysis of Floristic Diversity
Floristic composition and spontaneous floristic diversity were assessed using four main parameters:
1) Species richness (R) which represents the total number of species recorded;
2) Absolute Abundance (A) which represents the total number of individuals of a species  ;
3) The Shannon Weaver (H') floristic diversity index is an index expressed in bits between 0 - 5  (2). This index was obtained from the relative abundance (Pi) of each species, which is the ratio between the total number of individuals in a survey and absolute abundance  . Relative abundance is given by Equation (1).
represent the total number of species individuals i and N the total number of individuals.
4) The evenness index of Pielou (E) is between 0 and 1 (3). It is the ratio between the diversity of a species and the total number of species “S” of a sub-plot  . This index makes it possible to measure the distribution of individuals within species.
2.3.2. Distribution of Recorded Species by Diameter Class
This part concerns all spontaneous species of diameter ≥ 5 cm inventoried. The diameters of the trees were obtained in the field by dividing the circumference by 3.14  . These diameters were used to divide the studied stands by amplitude class.
2.3.3. Survival Rate
The number of plants installed being known, after inventory and measurement of the trees, we obtained the number of trees for each plot, which made it possible to calculate the survival rate given by Equation (4):
2.3.4. Biomass Estimation and Carbon Stock
1) Aboveground biomass (AGB)
The estimation of aboveground biomass was performed using the standard allometric equation developed by  (5). The choice of this equation was made based on the type of ecosystem in the study area and the threshold of DBH ≥ 5 cm. This equation has been recommended by the REDD+ process, successfully used and repeated in many studies    .
: wood density (g・cm−3);
H: height of tree (m);
D: diameter at breast height (cm);
AGB: aboveground biomass (t・ha−1).
2) Belowground biomass
Belowground biomass was obtained from the  equation for tropical forests  .
So “y” is belowground biomass (BGB) and AGB is aboveground biomass.
3) Total biomass
The total biomass of standing woody was obtained by adding the aboveground biomass and belowground biomass.
4) Carbon stock and ecological service value
The sequestered carbon stock was obtained by multiplying the total biomass by 0.47  . The amount of carbon dioxide (CO2) emissions was then determined by multiplying the sequestered carbon stock by the ratio 44/12 corresponding to the CO2/C ratio. This is the monetary value deducted from the amount of CO2 emissions captured from the atmosphere by the identified trees. This value was obtained by evaluating the amount of CO2 in very fluctuating monetary value which is 5 USD/t eq.CO2  .
Data analysis has been performed using PAST v.3.05 and QGIS v.2 softwares.
3.1. Floristic Diversity Overview
The inventory carried out in the blocks of Terminalia superba Engl. & Diels of 64, 31 and 20 years old has enabled us to record 9, 27 and 15 alive trees respectively (Table 1). Table 1 shows important differences depending on the surveys and blocks. The largest value of the diameter is observed at the level of the sub-plot 3 (P3), while that of the height is obtained at the level of the sub-plot1
Table 1. Trees recorded in study area by sub-plot. n: Number of individuals for each species recorded in study area; H: Average height of each tree in Bilala (m); DBH: Average of diameter at breast height for each species in study area (cm).
(P1). The 64-year-old blocks have an average diameter of 65.28 cm and an average height of 44.5 m. Floristic surveys carried out in the block of Terminalia superba Engl. & Diels aged 64, 31 and 20 years; and in the fallow field, the total of 62 species belonging to 36 families, around 3586 individuals. These species consist of woody, herbaceous and lianascent plants. Table 2 & Table 3 present the richness and floristic diversity of the studied blocks.
Table 2 shows that the number of individuals varies from 40 to 596. Species richness varies from one survey to another. In the fallow field, there are around 579 individuals in 13 families and 18 species. The large value was obtained in P3 and P8 (18 species). The Shannon diversity index and Pielou evenness ranged from 1.30 to 2.44 bits and 0.35 to 0.77, respectively, in all Terminalia superba Engl. & Diels surveys. In the fallow field, diversity is very important with a value of 3.02. Table 3, however, shows that the number of individuals varies between 680 and 1329. The family number varies between 10 and 24. Species richness in turn varies from 12 to 33 species. The greatest species richness is observed in 20-year-old and 64-year-old blocks with 27 species in 22 families. The lowest species richness is observed in 31-year-old blocks (12 species in 10 families). Thus, in the undergrowth of Terminalia superba Engl. & Diels, the floristic inventory yielded around 3007 individuals in 33 families and 52 species. However, Shannon indices ranged from 1.81 to 1.97 bits. The highest value is observed in the 31-year-old sub-plot. The Pielou index varies from 0.47 to 0.61. The highest average value of the Pielou index goes back to the 31-year-old sub-plot. The high values of the Shannon diversity indices and Pielou evenness shows that its stands are diverse and the individuals are well distributed among the species.
Figure 3 shows the percentages of species per family for the 64, 31 and 20 year old blocks. Figure 3(a) shows that 9 out of 22 families represent more than 67% of species. The most represented families are Burseraceae with 12% of species, followed by Euphorbiaceae with 10%, Moraceae with 8%, Rutaceae and Sapotaceae with 7% each. The Arecaceae family with 5% of species ranks last. The term
Table 2. Specific richness and diversity indices of the different surveys.
Table 3. Specific richness and diversity index of studies blocks.
Figure 3. Proportion of the most represented families in number of species in the 20-year-old block (a). Proportion of most represented families in number of species in the 31-year-old block (b). Proportion of most represented families in number of species in 64-year-old block (c).
“other” represents the 15 families that each has less than 5% of species. Figure 3(b) shows that 7 families account for 92% of the species recorded. The first rank is occupied by the Poaceae family with 27%. In the last row is the family Cucurbitaceae with 5% of species. Five other families represent 8% of the listed species. Five out of 24 families represent 65% of the species found in 20-year-old blocks (Figure 3(c)). The most represented families are: Fabaceae-Mimosoïdeae (24%); Euphorbiaceae and Fabaceae-Caesalpinoideae with each 13%; Annonaceae (9%); Meliaceae (6%). Nineteen other families represent 35% of recorded species.
3.2. Tree Species Distribution by Diameter Class
Seventy-eight woody plants identified in the undergrowth of Limba (Terminalia superba Engl. & Diels) blocks studied and distributed in 28 species (Figure 4) have different physiognomies depending on the diameter. Figure 4 shows the distribution of woody plants in the 64, 31 and 20-year sub-plots by amplitude diameter class 10.
Figure 4(a) shows that the total number of individuals is 44 trees. The largest number of individuals is represented by class 10 - 20 (13) followed by class 0 - 10 with 6 individuals. Class 80 - 90 was not represented. This Figure 4(a) shows that the diametric distribution of individuals has a decreasing exponential form (equation y = 12.464e−0.25x). Figure 4(b) shows that the largest number of individuals is represented by the 0 - 10 class with 2 individuals followed by the 10 - 20 class with only 1 individual. Classes from 20 to 120 were not represented. The distribution of individuals by diameter class once again shows a decreasing exponential form (equation y = 4e−0.693x). Figure 4(c) shows that the number of individuals is 32. The largest number of individuals is represented by the 0 - 10 class with 9 individuals followed by the 20 - 30 class with 7 individuals. Classes from 70 to 120 were not represented unfortunately. The diametric structure of individuals, however, has a decreasing exponential form of equation y = 10.60e−0.247x. Figure 4 shows in the different diameter classes that the number of individuals decreases with increasing diameter.
3.3. Biomass and Carbon Stock
The floristic inventory conducted in the blocks of Terminalia superba Engl. & Diels of 64, 31 and 20 years of age identified a total of 129 trees of DBH ≥ 5 cm consisting of 51 trees of Terminalia superba Engl. & Diels and 78 feet of various species, which allowed the calculation of biomass and carbon stock.
Table 4 presents the distribution of biomass by survey and pool. Table 4 presents the distribution of biomass by survey and pool. Table 4 shows that the largest biomass is represented by the P1 with 55.21 tons. From P2 to P9, a considerable variation of the biomass is observed. Table 5 presents the distribution of biomass and carbon stock by pool, block, sub-plot and species. However, biomass ranged from 46.88 to 131.60 tons in all blocks. The total aboveground
Figure 4. Distribution of species founded in 64-year-old blocks by diameter class (a); Distribution of identified species in 31-year-old blocks by diameter class (b); Distribution of species recorded in 20-year-old block by diameter class (c).
Table 4. Biomass distribution by sub-plot (in t・ha−1); P: Sub-plot; AGB: aboveground biomass (t・ha−1); BGB: belowground biomass (t・ha−1).
Table 5. Biomass and carbon stock distribution in study area by forest age. AGB: Trees aboveground biomass (t・ha−1); BGB: Trees belowground biomass (t・ha−1); Biomass: Total trees biomass (t・ha−1); Carbon: Total carbon stock (t・ha−1); Block: Block of planted forest by age (in year).
The largest biomass is represented by the oldest block (64 years) followed by the block of 31 years (60.82 t・ha−1) and lastly the stand of 20 years (56.49 t・ha−1). This Table 5 shows that the biomass varies with the age of the block. This distribution of biomass is also observed for the carbon stock.
With regard to the distribution of carbon stock per sub-plot, P1 sequestered more than 25 tons of carbon, followed by the P3 with a stock of 25.35 tons. In this ecosystem, for Limba (Terminalia superba Engl. & Diels) we recorded 181.42 t・ha−1 of biomass, or 85.27 t・ha−1 of carbon stock. For other species, we recorded 96.53 t・ha−1 of biomass, or 19.79 t・ha−1 of carbon stock. The total biomass of the study area was 275.89 tons and the carbon stock 129.67 tons of carbon. The quantity of CO2 captured in the atmosphere by the identified plant procession deducted from this carbon stock has been 475.89 tons of CO2, an economic value equal to US$ 2379.45 or XAF 1,413,278.
For distribution by diameter class of the carbon stock recorded at the 64-year block level, the largest value is represented by the class 50 - 60 cm (36.97 t・ha−1) followed by the 70 - 80 cm class (36.51 t・ha−1). The class of 100 to 110 was unfortunately not represented. For Distribution by diameter class of carbon stock sequestered by the 31-year-old block, the 60 - 70 cm (18.78 t・ha−1) class sequestered the largest carbon stock, followed by class 30 - 40 cm (12.97 t・ha−1). Classes 80 - 90, 90 - 100, 100 - 110 and 110 - 120 have unfortunately not been represented. For the diameter class distribution of carbon sequestered by the 20-year-old block, the largest carbon stock was sequestered by the 40 - 50 cm diameter class (23.02 t・ha−1). Diameter classes of 0 - 10 and 70 - 120 were not represented. However, for the height class distribution of the carbon stock recorded in this ecosystem, the distribution pattern of the carbon stock has been the same for all three blocks. The largest value of the carbon stock was represented by the class 30 - 40 m, followed by class classes 20 - 30 m, 10 - 20 m and 0 - 10 m. Classes 40 - 50 m and 50 - 60 m were not represented in the blocks of 31 and 20 years.
3.4. Survival Rate
This study shows different survival rates (4). The sub-plot1 has a survival rate of 100% followed by the second row P2 with a rate of 75%. P7, P9 and P8 with respectively 66.67%, 44.44% and 44.44% survival rate, high values compared to P4, P5 and P6 with respectively 30%, 26.67% and 36.67%. The monospecific stands of 64, 31 and 20 years represent the respective survival rates of 75%, 30% and 55.55%.
4.1. Dendrology and Survival Rate State
This study carried out within the Limba (Terminalia superba Engl. & Diels) plantations of 64, 31 and 20 years, allowed us to record respectively 9, 27 and 15 alive floras. These alive floras counts correspond to a survival rate of 75%, 30% and 55.55% respectively for blocks of 64, 31 and 20 years. The study conducted by  in the plot of Aucoumea klaineana Pierre with an area of 1.4 ha shows a survival rate of 44.05%. Thus,  stated that this rate is explained by a poor monitoring of the plot once the plants were installed. Regarding our study, the different survival rates recorded in the blocks of Terminalia superba Engl. & Diels highlighted is explained in various ways.
In fact, the 64-year-old block installed by a silvicultural method of artificial regeneration consisting in destroying the existing vegetation and then reforesting it, has given very good results as testifies Mr. Bakala, former worker of the OCF Company. He added that regular talks were given to his reforestation. Given this testimony, it’s clear to reject the hypothesis of poor monitoring of the block. However, this survival rate of 75% can be granted to large-scale anthropogenic actions in these blocks. This plantation is the most disturbed of our study, 3 strains found on the spot in the sub-plots 2 and sub-plot3 confirm the slaughter of 3 trees. As for the 31-year-old block, it was set up in a swampy area and at a distance of 5 m by 5 m by a method of total destruction of the wild gallery forest. This block has the lowest survival rate of the study (30%). This plot located in Bilala 100 m from the administrative building of the station to protect the soil of the creek “gutter” against erosion, remains immune to human disturbance. However, the factors that may cause this survival rate may be soil type and planting spacing.
4.2. Floristic Diversity Assessment
In terms of biodiversity conservation, the indicator commonly used is the specific richness that determines the state of degradation of an ecosystem. Overall, the plant procession studied reported 62 species, or around 3586 individuals belonging to 36 families. Species richness varies from one survey to another from 6 to 19 species. The Shannon diversity index and Pielou evenness ranged from 1.30 to 2.44 bits and 0.35 to 0.77, respectively, in all Terminalia superba Engl. & Diels surveys. In the fallow field, the diversity has been very important with a value of 3.02 compared to the records of Terminalia superba Engl. & Diels.
In all blocks of Terminalia superba Engl. & Diels, species richness ranges from 12 to 33 species. The greatest species richness was observed in 20-year-old blocks followed by 64-year-old blocks with 27 species. These two blocks of 64 and 20 years have a high number in terms of species compared to fallow field which has 18 species. This richness of 18 species was greater than that obtained in block of 31 years (12 species). These results do not corroborate those of  where the number of species increases with block age.  further reveals that floristic wealth is conditioned mainly by ecological conditions and anthropogenic disturbances  .
It is reasonable to say, however, that activities within the 64-year-old block have had an impact on forest recovery. These differences can also be explained by the silvicultural method applied for each block. Indeed, the 20-year-old block has been set up right in the heart of a gallery forest by Professor Aubreville’s technique of opening strips, which partially destroys existing vegetation. This block is very stable; no activity has been identified within this block. The 31-year-old and 64-year-old block were established by the destruction of the existing forest, which had a negative impact on forest recovery. The 31-year-old stand was established in a flooded area while  noted that wetlands have relatively low species richness.
Shannon indices range from 1.81 to 1.97 bits in all three stands. The Pielou index varies from 0.47 to 0.61. The highest average value of the Pielou index goes back to the 31-year-old plot. These strong values reflect a great diversity and a possible reconstitution of forest. These results are superior to those found by  in monocultures of Eucalyptus and Pines (0.38 to 0.89 bits). Differences in silvicultural crops and the age of their blocks could explain the differences between their indices.
Resettlement can also be characterized by the structural distribution of individuals by diameter classes  . It appears from this study that the individuals surveyed are highly concentrated in numbers in small diameter classes. The largest number of individuals is represented by the 10 - 20 class for the 64 year old block and 0 - 10 in the 31 and 20-year-old blocks. In all blocks of Terminalia superba Engl. & Diels, the distribution of individuals by diameter class has a decreasing exponential form. The number of individuals per diameter class decreases with the increase in diameter reminiscent of a recent block and good regeneration of vegetation. This distribution corroborates with the works from    ; where the general trend shows a distribution with more important numbers for trees of small diameters.
4.3. Trees Biomass and Carbon Stock
One hundred and twenty-nine floras of DBH ≥ 5 cm, estimated the biomass of the block at 275.89 tons divided into two pools: aboveground pool was 228.96 tons and that of the belowground pool was 46.94 tons. In all the plots, this biomass varies from 0.051 to 55.21. This is explained by the fact that the plots studied do not have the same physiognomy.
However, biomass varies from 56.49 to 158.58 tons in all the blocks studied. The largest biomass was represented by the oldest block (64 years) followed by the block of 31 years and lastly the block of 20 years. These differences are due to the diversity of each block and the specific density of the inventoried trees. These results corroborate with the study conducted by  , in the agroforestry systems of cocoa plantation from 5 to 15 years, where biomass increases with the age of planting  .
The carbon stock sequestered by all studied block was 129.67 tons of carbon. It ranges from 26.55 to 74.53 tons of carbon throughout the block. The largest carbon stock was sequestered by the Terminalia superba Engl. & Diels block with a value of 85.27 tons. The other 27 species represent a carbon stock of 19.27 tons of carbon. This difference is likely related to the diametric structure, sequestration capacity, and age of species that have been identified. The amount of CO2 captured in the atmosphere by the plant procession is 475.89 tons CO2 equivalent, an economic value equal to USD 2379.45 or XAF 1,413,278. This economic value obtained in 2018 is very fluctuating. The price of a ton of carbon in the carbon market is oscillating.
This work has shown that other forest species can exist in a monodominant forest of Terminalia superba Engl. & Diels. The distribution of exhausted flora according to diameter classes shows that stems with a diameter of less than 10 cm accumulate a total of 43 individuals while large diameter exhausted flora becomes rare. This presence of small diameter stems shows that the sub-plots studied have several individuals with a future to ensure regeneration. It is therefore appropriate to state that the vegetation of the undergrowth of Terminalia superba Engl. & Diels evolves towards the characteristics of the natural forest. Recolonization of exhausted flora by the forest is rapid. In the long term, this will probably lead to an installation of the climax forest; particularly the heliophilous species such as Ricinodendron heudelotii (Baill.) Pierre ex Heckel, Musanga cecropioides R.Br., etc., already reflect this phenomenon. The total biomass has been 275.89 tons and the carbon stock has been 129.67 tons of carbon corresponding to an amount of 475.89 tons of carbon dioxide equivalent, which is an economic value equal to XAF 1,413,278. This study has shown that strong anthropogenic disturbances play a very detrimental role in forest recovery, which prevents the development of trees and the generation of carbon stocks. This study highlighted the reforestation of Terminalia superba Engl. & Diels in its floristic diversity, which contributes to fight against the climate change. Finally, it proposes socio-economic integration of the local population of Bilala into agroforestry projects in order to diversify the activities of SNR, which generates spinoffs that can overcome the economic crisis to which the structure and population of Bilala are confronted with.
The authors would like to thank National Superior School for Agronomy and Forestry (Ecole Nationale Superieure d’Agronomie et de Foresterie―ENSAF) from Marien Ngouabi University (Université Marien Ngouabi) for supporting this study. Reforestation National Service (SNR) from Congolese Ministry of Forest Economy also supported this work. We greatly acknowledge Joseph Faustin Dembi, Guy Julien Kazotti and Gad Kinouani from Reforestation National Service (Service National de Reboisement―SNR) for providing facilities about field measurements in Bilala. Different anonymous referees have provided substantial contribution and the authors address to them their heartfelt thanks.
Supplementary material related to this paper is available online at:
 Pan, Y., Birdsey, R.A., Fang, J., Houghton, R., Kauppi, P.E., Kurz, W.A. and Phillips, O.L. (2011) A Large and Persistent Carbon Sink in the World. Science, 333, 988-993. https://doi.org/10.1126/science.1201609
 PEFC (2017) European Program for Forest Certification.
 Van Der Werf, G.R., Morton, D.C., Defries, R.S., Olivier, J.G.J., Kasibhatla, P.S. and Jackson (2009) CO2 Emissions from Forest Loss. Nature Geoscience, 2, 737-738.
 Ngoumba, B.Y. (2019) Biodiversity, Estimation of Carbon Stock and Forest Management for Terminalia superba Engl. & Diels Reforestation in Bilala (Kouilou Department). Master Memory. National Superior School for Agronomy and Forestry. Marien Ngouabi University, Brazzaville.
 Ekoungoulou, R., Liu, X.D., Loumeto, J.J., Ifo, S.A., Bocko, Y.E., Koula, F.E. and Niu, S.K. (2014) Tree Allometry in Tropical Forest of Congo for Carbon Stocks Estimation in Above-Ground Biomass. Open Journal of Forestry, 4, 481-491.
 Ekoungoulou, R., Niu, S.K., Folega, F., Nzala, D. and Liu, X.D. (2018) Carbon Stocks of Coarse Woody Debris in Central African Tropical Forests. Sustainability in Environment, 3, 142-160. https://doi.org/10.22158/se.v3n2p142
 Jonard, M., Colmant, R. and Heylen, C. (2014) Impact of Coniferous Forestation on Carbon Sequestration of Peruvian Andes: Case of Pinus patula Plantations in the Cooperative of Granja Porcon (Cajamarca), Perou. Bois et Forêts des Tropiques, 322, 17-27. https://doi.org/10.19182/bft2014.322.a31226
 Ekoungoulou, R., Liu X.D., Loumeto, J.J. and Ifo, S.A. (2014) Tree Above-And Below-Ground Biomass Allometries for Carbon Stocks Estimation in Secondary Forest of Congo. Journal of Environmental Science, Toxicology and Food Technology, 8, 9-20. https://doi.org/10.9790/2402-08420920
 Ekoungoulou, R., Liu, X.D., Ifo, S.A., Loumeto, J.J. and Folega, F. (2014) Carbon Stock Estimation in Secondary Forest and Gallery Forest of Congo Using Allometric Equations. International Journal of Scientific and Technology Research, 3, 465-474.
 Kimpouni, V., Loumeto, J.L. and Mizingou, J. (2008) Floristic Diversity of Forest Facies of Aucoumea klaineana (Okoume) in the Congolese Littoral. Acta Botanica Gallic: Botany Letters, 153, 323-334.
 Boukono, D.R.V. (2018) Estimation of Carbon Stock in a Block of Okoumea klaineana from the Reforestation National Service. Bachelor Memory. National Superior School for Agronomy and Forestry, Brazzaville.
 Ekoungoulou, R., Niu, S.K., Loumeto, J.J., Ifo, S.A., Bocko, Y.E., Mikieleko, F.E.K., Guiekisse, E.D.M., Senou, H. and Liu, X.D. (2015) Evaluating the Carbon Stock in Above-And Below-Ground Biomass in a Moist Central African Forest. Applied Ecology and Environmental Sciences, 3, 51-59.
 Goma-Tchimbakala, J., Ndondou-Hockemba, M., Kokolo, A. and Mboussou-Kimbangou, A.N.S. (2005) Variations on Litter Bringing and Mineral Component in the Limba (Terminalia superba) Plantation in Congo. Tropicultura, 23, 53-59.
 Pedel, L. and Fabri, M.C. (2012) Art State on Existing Index Concerning Ecological State of Benthic Habitats of Deep Domain. Laboratory of Environment Resources (LER). MEDDTL Convention, Ifremer for DCSMM, BEE.
 Ekoungoulou, R., Nzala, D., Liu, X.D. and Niu, S.K. (2017) Ecological and Structural Analyses of Trees in an Evergreen Lowland Congo Basin Forest. International Journal of Biology, 10, 31-43. https://doi.org/10.5539/ijb.v10n1p31
 Ekoungoulou, R., Folega, F., Mukete, B., Ifo, S.A., Loumeto, J.J., Liu, X.D. and Niu, S.K. (2018) Assessing the Effectiveness of Protected Areas on Floristic Diversity in Tropical Forests. Applied Ecology and Environmental Research, 16, 837-853.
 Chave, J., Rejou-Mechain, M., Burquez, A., Chidumayo, E., Colgan, M.S., Delitti, W.B.C., Duque, A., Eid, T., Fearnside, P.M., Goodman, R.C., Henry, M., Martinez-Yrizar, A., Mugasha, W.A., Muller-Landau, H.C., Mencuccini, M., Nelson, B.W., Ngomanda, A., Nogueira, E.M., Ortiz-Malavassi, E., Pelissier, R., Ploton, P., Ryan, C.M., Saldarriaga, J.G. and Vieilledent, G. (2014) Improved Allometric Models to Estimate the Aboveground Biomass of Tropical Trees. Global Change Biology, 20, 3177-3190. https://doi.org/10.1111/gcb.12629
 Djomo, A.N., Picard, N., Fayolle, A., Henry, M., Ngomanda, A., Ploton, P., McLellan, J., Saborowski, J., Adamou, I. and Lejeune, P. (2016) Tree Allometry for Estimation of Carbon Stocks in African Tropical Forests. Forestry, 89, 446-455.
 Ekoungoulou, R., Nzala, D., Liu, X.D. and Niu, S.K. (2018) Tree Biomass Estimation in Central African Forests Using Allometric Models. Open Journal of Ecology, 8, 209-237. https://doi.org/10.4236/oje.2018.83014
 Fayolle, A., Loubota, P.G.J., Drouet, T., Swaine, M.D., Bauwens, S., Vleminckx, J., Biwole, A., Lejeune, P. and Doucet, J.L. (2016) Taller Trees, Denser Stands and Greater Biomass in Semi-Deciduous than in Evergreen Lowland Central African Forests. Forest Ecology and Management, 374, 42-50.
 Mokany, K., Raison, R.J. and Prokushkin, A.S. (2006) Critical Analysis of Root: Shoot Ratios in Terrestrial Biomes. Global Change Biology, 12, 84-96.
 IPCC (2006) The 2006 Directrix Line of International Panel on Climate Change (IPCC) for the Nationals Inventories of Greenhouse Gas. Agriculture, Forestry and Others Earth Affections. Chapter 3: Earth Coherent Representation.
 EMP (2018) State of the Voluntary Carbon Market. Ecosystem Services. Ecosystems Marketplace. http://www.ecosystemmarketplace.com
 Tsoumou, B.R., Lumande, K.J., Kampe, J.P. and Nzila, J.D. (2016) Estimation of Carbon Stock Sequestred by Dimonika Reference Forest (Southwestern Republic of Congo). Revue Scientifique et Technique Foret et Environnement du Bassin du Congo, 6, 39-45.
 Temgoua, L.F., Dongmo, W., Nguimdo, V. and Nguena, C. (2018) Ligneous Diversity and Carbon Stock of Agroforestry System Based on Cocoa Plantations at East of Cameroon. Journal of Applied Biosciences, 122, 12274-12286.