Skip to main content

Effects of tree species diversity and stand structure on carbon stocks of homestead forests in Maheshkhali Island, Southern Bangladesh



The homestead forests of Bangladesh occupy 0.27 million hectares (10% of the total forested area) and have potential to store carbon (C) and conserve biodiversity. Small scale forestry practices, however, are lacking reliable estimation of C stocks and tree species diversity. This may hinder successful implementation of REDD + and similar mechanisms as they concentrate on large-scale forests. This study aimed to estimate the above- and below-ground carbon stocks in homestead forests of Maheshkhali Island in Bangladesh and how tree species diversity and stand structural variation affect these C stocks. We randomly surveyed a total of 239 homestead forests in the hillside, beachside, and inland in 2019.


Tree biomass C stocks were 48–67% greater in the inland and hillside forests than in the beachside due to significantly greater stand density, basal area, tree diameter. In total we found 52 tree species, but most abundant species in the inland and hillside forests, Mangifera indica, Samanea saman, and Artocarpus heterophyllus stored the most C in tree biomass. Greater tree species richness and diversity index in the inland and hillside forests indicated greater above- and below-ground tree biomass C stocks. An increase in tree species richness and diversity index by one unit was found to increase the tree biomass C stock by 22 and 30 Mg C ha−1, respectively. The total soil C stock was also affected by tree species diversity, stand density, and their interaction with soil properties. Total soil C stocks were greatest (51 Mg ha−1) in the inland forests, having also the greatest stand density and tree species richness. C stock in soil surface was greatest in the hillside forests due to the greatest litterfall, but the average share of litterfall from the total biomass C was only 0.1%.


Homestead forest ecosystems could store 96 Mg C ha−1 in total, which can contribute to climate change mitigation by generating C credits for small-scale homestead forests owners. Above- and below-ground tree biomass C stocks were found to correlate with tree species diversity, which may also contribute to biodiversity conservation in the REDD + in Bangladesh and countries alike.


A gradual increase in the global emissions of carbon dioxide (CO2) and consequent temperature increase has become a major concern to work on emissions mitigation [1,2,3]. Tropical forests play an important role in removing atmospheric CO2 as they store one fourth of the global terrestrial carbon [4,5,6]. In addition, tropical forests support at least two-thirds of the world’s biodiversity (e.g., [7]) and have 50% of all known plant species, while their coverage from the total land area of the Earth is about 12% [8].

In Bangladesh, CO2 emissions over the time are rapidly increasing, for example, 609% in 2017 (78 Mt of CO2) compared to that of 1990 (11 Mt of CO2) due to an increase in energy consumption [9, 10]. Bangladesh´s contribution to global emissions is very low; however, its carbon rich forest ecosystems are highly affected by land-use change and adverse impacts of climate change [11], such as changes in precipitation and global sea level rise. Bangladesh has forest areas of about 2.53 million hectares, representing 17.5% of the total land area [12, 13]. These tropical forests are consisted of hill, mangrove, sal forests, and coastal mangrove plantations, with semi-evergreen and deciduous tree coverage. Natural hill and sal forests are being degraded due to illicit felling, shifting cultivation, and conversion to other land uses [14, 15]. Apart from these natural and planted forests, tree outside forest (TOF) including homestead forests, roadside plantation is booming in Bangladesh [13].

The homestead forests of Bangladesh occupy 0.27 million hectares land area, representing 10% and 2% of the total forested land area and total land area, respectively [13, 16] and have potential to store carbon in biomass [17, 18]. The contribution of homestead forests to rural economy is second to agriculture and these forests provide people’s daily needs. The homestead forests of Bangladesh supply 70% of total timber and 90% of fuelwood and bamboo demand in the country [13, 19] and thus, release pressure on natural forests. The homestead forests can be characterized as well-established land use systems for sustenance and conservation of biodiversity [20], which are maintained by at least 20 million households. Since homestead forestry is practiced primarily for supplying daily necessities as a livelihood option, understanding the carbon stocks of the homestead forests is required to address their potential in climate change mitigation. The homestead forests in Bangladesh are in pressure due to fragmentation of landholdings [21] and not under the national forest management plan.

A few researches have estimated above-ground forest carbon stocks in Bangladesh. Carbon stocks were found to vary with land uses, including mangrove and coastal (99 Mg C ha−1; [22]), protected contiguous and fragmented (34–53 Mg C ha−1; [15]), bamboo (53 Mg C ha−1; [23]), hill (103 Mg C ha−1; [16]), and homestead (53 Mg C ha−1; [17]) forests. These carbon stocks have been found to be dependent on the stand structure (e.g., tree height, DBH, density, basal area) [24, 25] and tree species [26, 27], and stands with fast-growing tropical tree species having the highest forest carbon stocks (201 Mg C ha−1; [28]). Tree species diversity may increase above-ground biomass carbon stocks of tropical forests [29, 30].

The carbon stock in litterfall is only a small fraction relative to above-ground biomass carbon in forest ecosystems [31], but this needs to be studied when estimating carbon dynamics among pools [32]. A balance between accumulation and decay of litter controls the accumulation of organic matter in an ecosystem [33]. Within the same climate, forest and tree species types are the main drivers of the litterfall, e.g., in mixed species natural forests and planted forests [34]. Litterfall is also affected by the management such as tree harvest and pruning [35,36,37]. Research on carbon stock in litterfall has been very scarce in the homestead forest of Bangladesh, and generally in tropical and temperate forests [38].

In Bangladesh, soil carbon concentration has been found variable (1–16 mg g−1), mainly responsible to the site and depth in soil [39, 40]. Soil carbon stocks estimated in Bangladesh were 23 Mg C ha−1 in semievergreen [41], 34 Mg C ha−1 in mangrove [22], and 59 Mg C ha−1 in deciduous [42] forests. Earlier, changes in stand structure and litter quality have found to modify the soil carbon dynamics in agroforestry ecosystems and also in tropical homestead forests [43,44,45,46]. The dynamics are also influenced by microclimatic and edaphic conditions [47,48,49] and they vary in space and time [50]. For instance, tree size, stand density, and species richness positively affected soil carbon in tropical forests [51, 52]. While significant advances in estimating the carbon balance of forests have been attained, there are still critical uncertainties in the magnitude of soil carbon stocks [37].

The United Nations Framework Convention on Climate Change (UNFCCC) introduced mitigation instruments, clean development mechanism (CDM) and reducing emissions from deforestation and forest degradation and conservation and enhancement of forest carbon stock with sustainable management (REDD +) [2, 53, 54]. REDD aims to maintain carbon stock within tropical forests while conserving biodiversity [55]. However, there is uncertainty in biodiversity provision within REDD [56] and lack of clear understanding on interactions of carbon dynamics and biodiversity [57]. A positive corelation between estimated biomass carbon and biodiversity exits globally [58], but there is spatial variation [57], which makes REDD initiative complex at regional or sub-national scales. Moreover, inadequate data on carbon stock in local forests (Baccini et al. [59]) and the lack of reliable estimation of tropical forests carbon stocks may hinder the effective implementation of REDD + and similar mechanisms [60]. More importantly, REDD, or derivative REDD + , takes only large-scale forests into account, while ignoring the small-scale forests, such as homestead forests [61]. The evidence of carbon sequestration potential of TOF can be significant for small-scale forest landowners or households in developing countries, such as Bangladesh, and it can be also used to support international treaties of the Paris agreement 2015 and the Kyoto Protocol 1997. In addition to the natural and planted forests managed by Bangladesh Forest Department, the estimation of the carbon stocks in homestead forests is imperative for investigating their potentials for carbon enhancement and credits [61, 62].

Under this circumstance, the study about homestead forests for estimation of carbon stocks would be the scientific-based information for the policymakers and scientists with a view to support climate change mitigation efforts. Our study aims to estimate the carbon stocks in homestead forest ecosystems (trees, litterfalls, and soil) of Maheshkhali Island under Cox’s Bazar District in Bangladesh and how tree species diversity and stand structural variation affect these carbon stocks.

Material and methods

Study area

The study was conducted in homestead forests of Maheshkhali Island under Cox’s Bazar District, a coastal area of Bangladesh (Fig. 1), emerged as most vulnerable to climate change impacts [11]. Maheshkhali island is the only hilly island with complex geological system on the eastern cliff coast of Bangladesh [63], located between 21°28′ and 21°46′N latitude and 91°51′ and 91°59′E longitude [64]. It occupies an area of 362.18 km2, with a total of 33,287 households [65]. The island has a moist tropical climate with a long wet season (April–October) and a relatively short dry season (November–March). The mean annual precipitation, temperature, and relative humidity are 3627 mm, 25.7 °C, and 70–90%, respectively [66]. This region is prone to cyclonic storms, tidal surges, and flood due to proximity to the Bay of Bengal, a source of cyclones, usually occurring during April–May and October–November.

Fig. 1

Map of (a) Bangladesh and (b) Maheshkhali Island showing three categorized study sites with (c) sampling points of homestead forests

The island has four subdivisions including active, young, and old coastal plain, and hilly areas with piedmont plain. Geological deposition of sedimentation forms landmasses [63]. Maheshkhali with an accretion rate of 1.2 sq. km. per year since 1972 formed huge landmasses in the southwest coastal plain (e.g., Bara Maheshkhali) and western coastal plain (e.g., Gorakghata) [67], contributing to the land use and land cover changes. The major land uses include salt cultivation, agricultural land, hill forests, and coastal forests, which have been changed markedly since 1972. Expansion of salt fields caused a decline in the agricultural land at an average rate of 14.5 ha per year, and extensive and illegal hill cutting for settlement, betel leaf cultivation, and unpanned development caused reduced hill forests at an average rate of 90 ha per year [64]. Shrimp cultivation, another type of land use threatened the coastal forests, to some extent, especially old coastal zone [64].

Bangladesh Forest Department manages its hill and coastal forests consisting of mangrove plantations of 534.54 ha and 8129.9 ha and non-mangrove plantations of 2667.3 sq. km and 232.66 sq. km, under two range,Footnote 1 Maheshkhali and Gorokghata, respectively [68]. Forest covers several hills of up to 23 m and the low-laying valleys. Soils of the forest vary from clay to sandy loam and to some extent yellowish red sandy clay [69]. Peoples` dependence on hill forests were significant for collecting fuelwood, house and boat making materials, traditional medicine, and non-timber forest products such as bamboo, honey, fodder etc. However, from the time immemorial, this overexploitation of the resource declined in biodiversity in the area [70]. This situation made the individuals for taking care and managing homestead forests very carefully for their protection from coastal storms, surges, floods. Homestead forests are managed by owners themselves.

Reconnaissance survey

Before starting the data collection, three initial field visits for reconnaissance survey were made to get an overview of the study area in February 2019. This included observation of general conditions including geographical location, physiography, hill and coastal forests, and existence of homestead forests along the hillside, beachside, and in the inland of the area. We categorized three UnionsFootnote 2 Chhato Maheshkhali (recently converted as municipality), Gorokghata, and Bara Maheshkhali as the hillside, beachside, and inland, respectively, under Maheshkhali sadar Upazila according to the geographic location (Fig. 1). Settlement was assumed to be associated with the homestead forests in the three sites. Settlement is an important type of land use, and the population density was relatively higher in the hillside (Chhato Maheshkhali), beachside (Gorokghata), and inland (Bara Maheshkhali) compared to the northwestern part of the island [71]. From the key informant (KI) interview with the Chairman of the Union in February 2019 it was known that some of the settlements started in the foot of the hills under the hillside by human intervention in modifying slope of the hills. Since after the loss of lives in 1991 cyclone, people started migrating from Gorakghata to other places [67]. Settlement started earlier in the inland site and population was relatively higher [71]. The homestead forests of these three sites were also assumed to represent the same kind of ecotype. We hypothesized that carbon stocks of the homestead forests differ from each other among the three sites.

The researchers collected relevant data of villages, number of households, and homestead forests of each village from the respective administrative (Union Parishad) offices. Additionally, we interviewed the Chairman of these three Unions as key informants to gather knowledge about the study sites and as well as inform them about the purpose of the study.

Sample selection and woody vegetation measurement in homestead forests

The sampling procedure followed from Upazila to Union, Union to village, village to homestead forests of the households. From the lists of the number of households provided by the office, with a sampling intensity of 5% as accepted by the United Nations [72], a total of 239 homestead forests were determined. Then, based on the total number of households in each of the three sites, 67, 69, and 103 homestead forests from hillside, beachside, and inland, were randomly allocated for the study in Maheshkhali sadar Upazila in 2019 (from February to April). The mean area of the studied homestead forests in these three sites were 0.02, 0.01, and 0.02 ha per household.

Each of the homestead forests was divided into quadrats (5 m × 5 m) based on the area and the direction from the dwelling. The surveyed data were recorded which included all woody plants identification, with measurement of height (m), diameter at breast height (DBH, cm) and the area of the homestead forests. The owners of the homestead forests helped in identification with local name, and in few cases, herbarium was prepared to ensure the identification with scientific names. The height measurement was made by rangefinder and DBH by diameter tape. The coordinates of each point of sample collections was recorded by using GPS. Herbs and shrubs were not considered as 98% of total forest biomass consists of tree biomass; they may be ignored in estimation of carbon [73]. Homestead forests are well managed and therefore, are usually free from herbs and lianas.

Soil and litterfall sampling for estimating C stock

A sampling of the litterfall was made in 4 points wherever available for each of the three different sites in 2019 (from February to April), thus making a total of 12 (3 × 4 = 12) samples. All litterfalls at each point of an area of 1 m2 (1 m × 1 m) was collected using a metallic frame. A pit of 30 cm depth, under the litterfall layer sampling point, was dug by using a soil auger and mineral soil samples were collected at 10, 20, and 30 cm depths. This procedure was followed for four samples consisting of 12 (4 × 3 depths = 12) subsamples for each of the three different sites, thus making a total of (12 × 3 = 36) subsamples. Accordingly, following the same procedure, 36 unaltered soil subsamples were collected using a core (volume 100 cm3) to measure bulk density (BD) at the same three depths in each point, following Blake [74].

Data analyses

Estimation of tree (above- and below-ground) biomass, and density and basal area of stands

Above-ground biomass (AGB) was estimated by converting tree data into biomass using allometric Eqs. (1), (2), (3), and (4) for tropical trees, Cocos nucifera, Areca catechu, and Phoenix dactylifera, respectively (Table 1) [60, 75,76,77]. Below-ground tree biomass (BGB) was estimated as 15% of AGB [78]. Tree total biomass (TB) was the summing up of AGB and BGB. Finally, total carbon stock (Mg ha−1) was estimated as carbon content is assumed to be 50% of dry TB [79]. To estimate AGB, wood density (g cm−3), a required variable, which was collected from Bangladesh Forest Research Institute (BFRI) [80]. For those not found in BFRI publications we used global wood density database [81, 82]. Additionally, species level carbon was also estimated for most frequent tree species and expressed in kilogram (kg) carbon per individual across three sites. Stand density (individual ha−1) and basal area BA (m2 ha−1) were estimated (Eqs. 5–7). Mean values of tree biomass, density, and BA were compared among three different homestead forest sites.

Table 1 Equations used in analyses of data

Laboratory analysis and estimation of carbon of litterfall and mineral soil and bulk density (BD)

To estimate soil organic carbon (SOC), soils were oven-dried at 105 °C for 72 h. After cleaning, washed silica crucibles were dried in an oven at 105 °C for half an hour and cooled in desiccators, and then mass was taken. Oven-dried soils were ground by pestle mortar and then exactly 5 g of grind soils were kept in silica crucibles and reweighed by an electric balance. The crucibles with soil were then transferred to an electric muffle furnace for igniting at 850 °C for one and a half hour. Then crucibles with soils were cooled in the desiccator and reweighed to determine the percent loss of ignition LOI (%), from which, SOC (%) was calculated (Eqs. 8 and 9). C stocks in mineral soil at three depths were calculated using BD (g cm−3) (Eq. 10) and expressed in Mg ha−1 for three different sites (Table 1). We calculated soil BD as the quotient between the dry mass of the fine fraction in the core segment and volume of the cylinder [86].

Regarding the estimation of biomass of litterfall, after taking the fresh mass of the original samples collected from each point of litter collection, adequate subsamples from the weighed original sample were made and labelled. In each plot, the number of original samples was four and subsamples three to five, depending on the wet masses of the original samples. The wet masses of all the subsamples were measured and recorded. Subsamples were oven-dried at 65 °C until reaching a constant mass and dried masses were recorded. Then, the dry mass of the original sample from the wet to dry ratio of the subsamples was estimated (Table 1; Eqs. 11 and 12). The carbon concentration was considered to be 50% of the dry mass of litter [88]. The process was repeated for all 12 original samples collected from homestead forests across three different sites. Carbon stocks in litterfall were calculated and expressed in Mg C ha−1 for three different sites. These carbon concentrations and stocks of litterfall and soils were compared among three different homestead forest sites.

Estimation of tree species richness, diversity and relative frequency and relative density

Tree species richness (Margalef index) and diversity (Shannon-Weiner Index, H) were estimated according to Eqs. 13 and 14 (Table 1). The greater value of indices of diversity indicates greater species richness and diversity in an area. In addition, the relative frequency of occurrence (RF %) and relative density (RD %) for species were estimated (Eqs. 15–17). Mean values of tree height (m), DBH (cm), all indices, RF, and RD were compared among three different sites.

Statistical analyses and modelling work

We run Kolmogorov–Smirnov (K–S) Test and found the data were normally distributed. Therefore, one-way analysis of variance (ANOVA) and Tukey’s post hoc test were used to determine whether there are any statistically significant differences (p ≤ 0.05) between the three homestead forest sites and which site significantly differed from the other sites in tree biomass (Mg C ha−1), height (m), DBH (cm), density (individual ha−1), basal area, BA (m2 ha−1), Margalef richness index and Shannon–Wiener diversity index. Moreover, two-way analysis of variance (ANOVA) was performed to determine whether there are any statistically significant differences (p ≤ 0.05) of soil carbon stock (Mg C ha−1) against three sites and three soil depths.

Relationship between tree biomass carbon stock (Mg C ha−1) and tree (a) height (m), (b) DBH (cm), (c) density (individual ha−1), (d) basal area, BA (m2 ha−1), (e) Margalef richness index, and (f) Shannon–Wiener diversity index were modelled by using linear regression analysis. In addition, multiple regression analysis was used to model the effect of all variables (a–f) to tree biomass. The regression analyses were performed to determine the what independent variables contribute to the explanation of the biomass carbon stock and to what degree. All these statistical analyses were performed by using Statistical Package for the Social Sciences (SPSS) version 26.


Stand structure and tree species diversity in homestead forests

We found the greatest mean tree DBH and height in the hillside and lowest in the beachside homestead forests, with significant (p ≤ 0.05) differences among the three sites (Table 2). The greatest number of trees were in 16–20 cm DBH class, with the inland homestead forests site being dominated. The trees with large DBH (31–45 cm) were greater in the inland compared to the other two sites (Fig. 2a). However, there were only few trees with DBH < 15 cm due to the harvesting of those at pole stage to be used as fuels. Related to height, the greatest number of trees were in 6–9 m class. Taller trees (10–13 m) were greater in inland and hillside homestead forests (Fig. 2b). The tree species diversity and richness of homestead forests were significantly (p ≤ 0.05) greater in the hillside and inland, compared to those on the beachside (Table 2).

Table 2 Mean values of tree DBH, height, species diversity and richness indices in the homestead forests
Fig. 2

Tree (a) DBH (cm) and (b) height (m) classes in homestead forests

Mean stand density and BA of homestead forests were 601 individuals ha−1 and 27 m2 ha−1, respectively, across the study area. Regarding the site, stand density in the inland homestead forests was significantly (p ≤ 0.05) greatest, compared to that in the other two sites, while no significant difference in these two (Fig. 3). BA of homestead forests was significantly (p ≤ 0.05) greater in the inland and hillside, compared to that of other one (Fig. 3).

Fig. 3

Stand density (primary y axis) and basal area (secondary y axis) in the homestead forests. Bars represent standard error of mean

Among 52 tree species found in the homestead forests, the numbers of species in the hillside, beachside, and inland were 41, 42, 48, respectively (Table 3). The number of tree individuals sampled were 840, 540, and 1504, respectively. The most five frequent species across the area were Mangifera indica, Acacia auriculiformis, Cocos nucifera, Artocarpus heterophyllus, and Samanea saman and these also corresponded to the RD (Table 3).

Table 3 Relative frequency (RF, %) and relative density (RD, %) of species found in the homestead forests

Tree (above-and below-ground) biomass and litterfall carbon in homestead forests

Mean tree (above-and below-ground) biomass in the homestead forests was estimated to be 46.11 Mg C ha−1 across the study area. Tree biomass was significantly (p ≤ 0.05) greater in the hillside and inland, compared to that in the beachside homestead forests (Fig. 4). The mean dry biomass of the litterfall was 0.04 ±  0.01 Mg C ha−1 across the study area. These were 0.06 ± 0.02,  0.04 ± 0.01, and 0.03 ± 0.02  Mg C ha−1 in the hillside, beachside and inland, respectively.

Fig. 4

Tree (above-and below-ground) biomass in the homestead forests. Bars represent standard error of mean

Among the species, Samanea saman dominated in storing carbon, with Mangifera indica, Artocarpus heterophyllus, Diptercarpus turbinatus, and Albizia procera stored relatively greater amount of biomass carbon in the inland and hillside homestead forests, compared to the beachside forest (Fig. 5).

Fig. 5

Tree (above-and below-ground) biomass for most frequent species in the homestead forests. Bars represent standard error of mean

Carbon concentration and stocks in mineral soil

Carbon concentration and stocks diminished with increasing depth of soil in homestead forests across three sites (Table 4; Fig. 6). The greatest mean carbon stocks were found throughout the soil depths in the inland homestead forests. However, carbon stocks did not significantly (p ≤ 0.05) vary with sites and with soil depths (Fig. 6). The bulk density of soils increased with depth for all sites (Table 4).

Table 4 Soil carbon (C) concentration and bulk density (BD) at three different depths in the homestead forests
Fig. 6

Carbon stocks in mineral soil at three different depths in homestead forests. Bars represent standard error of mean

Relationship of tree biomass with structural compositions in homestead forests

Figure 7 shows the significant (p ≤ 0.05) relationship between tree biomass (Mg C ha−1) and height (m), DBH (cm), density (individual ha−1), basal area, BA (m2 ha−1), Margalef richness index, and Shannon–Wiener diversity index in the homestead forests across three sites. Multiple regression analysis reveald that 90% of the variability in biomass C was explained by these factors together.

Fig. 7

Relationship between tree biomass and (a) height, (b) DBH, (c) stand density, (d) basal area, (e) Margalef index, and (f) Shannon–Wiener index in the homestead forests


Taking urgent action to combat climate change and its impacts are amongst the Sustainable Development Goal 13, which provides us with a common plan and agenda to tackle climate change [92]. Storing carbon in a forest ecosystem helps in removing CO2 emissions from the atmosphere and thus, contributing to the climate change mitigation. This requires growing trees in and outside the large-scale forests, for instance in homestead forests. Realizing this potential, this study estimated the carbon stocks in the homestead forest ecosystems in an Upazila of Maheshkhali island, a hilly and coastal area in Southern Bangladesh, and estimated how tree species diversity and stand structural variation affect these carbon stocks.

Mean biomass carbon stock estimated (46 Mg C ha−1) in our case was close to that found (54 Mg C ha−1) in homestead forests in northern Bangladesh [17]. This, however, is lower than that found in the mangrove (99 Mg C ha−1; [22]) and total forests (49–121 Mg C ha−1; [12]) of Bangladesh. This disagreement could be explained by the lower overall species diversity and richness in our study, which indicates lower biomass C stock [15]. For example, Nath et al. [18] estimated tree biomass of 118 Mg C ha−1, with a diversity index of 2.21 in homestead forests, whereas in our case the index was 1.24. However, we found the significantly positive effects of tree species diversity and richness on biomass carbon stock. Greater species richness and diversity index in the inland and hillside homestead forests indicated higher above- and below-ground biomass carbon stocks compared to that in the beachside. An increase in species richness and diversity index by one unit increased the biomass carbon stock by 22 and 30 Mg C ha−1, respectively (Fig. 7e, f). Our findings agreed with earlier studies in tropical forests of Asia and Africa [29, 30, 93], implying the more tree species diversity and richness the more likely higher above-ground biomass carbon stock.

Our study revealed a greater contribution of some of the frequently occurred species to the total carbon stocks. The relatively greater number of individuals of Mangifera indica, Samanea saman, Artocarpus heterophyllus, and Dipterocarpus turbinatus in the inland and hillside homestead forests, contributed to the greater carbon stocks compared to the beachside. Alamgir and Al-Amin [26] also found a greater biomass carbon stock in these tree species in the hill forests of Bangladesh. The strongly positive effects of BA and stand density on biomass carbon stocks are also generally in line with the findings of [94] who reported that greater density (4258 trees ha−1) and BA (53 m2 ha−1) increased biomass carbon stocks in roadside plantations.

In this study, with greater DBH and height, above- and below-ground biomass carbon stocks were 48–67% greater in inland and hillside homestead forests than in beachside forests. The overall share of individuals with DBH of 31–40 cm was also greater in the hillside and inland homestead forests which contributed to the greater above- and below-ground biomass carbon stock in comparison with that of the beachside. In our case, when tree height and DBH increased by one unit each, the biomass carbon stock increased by 11 and 3 Mg C ha−1, respectively (Fig. 7a, b). The importance of contribution of larger trees to the biomass C stock is in line with [24] who depicted that individuals with DBH of 10–56 cm, constituting only 28% of stand density, contributed 84% of the total biomass carbon stock in mangrove forest.

The carbon stocks of tropical litterfall have not received much attention in research as it constitutes a small fraction of above-ground biomass [32, 95]. The overall carbon stored in litterfall was 0.1% of the total biomass C in this homestead forest, while it was 1.8% in the natural forests of Bangladesh [34]. The carbon stock in litterfall was greatest in the hillside, which was up to 53–83% greater than that in the beachside and inland homestead forests. Litter accumulates in natural forests as they are no longer under silvicultural management due to the harvesting restriction, while litters in homestead forests are used as cooking fuel [96, 97], which leads to lower C stocks.

There was a clear decline in soil carbon concentration and stocks across three homestead forest sites. Total soil carbon stock in the inland homestead forests across the three depths was greatest (51 Mg ha−1), with increased stand density and species richness. The greater carbon in soil was correlated with greater stand density, and species diversity and richness, which has been found earlier in tropical agroforestry systems [46, 52] and temperate forests [43, 98]. The surface soil had 5–38% and 29–75% greater carbon stock of that stored at a depth of 20 and 30 cm, respectively, depending on the site. Hillside forests with greater litterfall, had 8–9% greater surface soil carbon stock compared to beachside and inland forests. Acacia auriculiformis, Acacia mangium, and Swietenia mahagoni species were abundant in the hillside and inland forests, contributing to soil C stock. This was because, Acacia and Mahagani litters were not preferred as fodder or fuels due to being small leaflet and unpalatable, contrary to Mangifera indica and Artocarpus heterophyllus [35]. Acacia species planted site in Bangladesh and African Mahogani in Ghana were found enhancing soil carbon stock [35, 41, 99]. Earlier studies on homestead forest and agroforestry in India also reported that the slower decay rate of Acacia and S. mahagoni litters resulted in accumulation of organic matters in the soil, compared to M. indica, A. heterophyllus and Anacardium occidentale [44, 100].

The relationships found between both the stand structure and tree species diversity with the biomass carbon have some important implications for emission reduction under the REDD and its derivative REDD + programme. The homestead forests with high floral diversity and biomass carbon indicate their high conservation potential. Therefore, carbon storing in homestead forests can provide co-benefits of biodiversity conservation under the REDD programme as it aims to maintain carbon stock in tropical forests while protecting threatened tree species [55, 58]. In this study, Garcinia cowa and Vitex peduncularis appeared as rare species [101] and can be conserved by protecting them from further erosion. While state-owned forests decline in Bangladesh, reducing carbon stocks and biodiversity in forests, the biodiverse-rich homestead forests are in a crucial role in enhancing carbon sinks, reducing emissions from deforestation, and contributing to global carbon cycle. These forests are managed for variable household livelihood options, which could provide carbon credits under REDD + programme. An appropriate management in community forests in Nepal have contributed to REDD + and local livelihoods [102,103,104], for example. However, a regulatory framework would be required to take homestead or small-scale forests under REDD and REDD + programme for policy initiatives to safeguard carbon, biodiversity, and local livelihood.

At local level, stand structural traits (DBH, height, BA, density) can easily be measured in the field by homestead forest owners or local communities. Based on these field data from small spatial scale, one can also assess biomass or produce map over large area using remote sensing techniques, and estimate national carbon storage and deforestation in TOF for REDD + monitoring.

Conclusions and policy implications

The hillside and inland homestead forests stored remarkable amount of tree biomass carbon, which was significantly increased with increasing stand density, BA, DBH, height, species richness and diversity. A smaller carbon stocks in litter found in our study compared to earlier studies could be linked to the removal of litter for using it as fuels. Higher litterfall in the hillside homestead forests may have contributed to surface soil carbon stocks, but the overall soil carbon stock in the study was also affected by types of litter of species, the stand density, and species richness. However, the decay of litter and humus and underground process in tropical forests influencing soil carbon stock, depending on environmental factors, would need to be studied further.

In Bangladesh, total annual emissions of the energy sector are 78 Mt of CO2 [9]. Conversely, the land use and land use change and forestry (LULUCF) sector’s carbon sink was 81 Mt in 2010 [105]. According to our study, the homestead forest ecosystems (trees, litterfall and soil) store 96 Mg C ha−1, which is 73 and 62% of that found in mangrove and hill forests (133 and 154 Mg C ha−1, respectively [16, 22]. This can be upscaled to be 26 Mt considering total homestead forest area of the country, which thus can contribute to climate change mitigation through REDD + and CDM mechanisms as emphasized in UNFCCC’s mitigation strategies [2, 106].

This study reduced the gap with documentation and producing estimated carbon of homestead forests that would help for applying REDD + mechanism since the contribution of TOF in carbon sequestration is ignored due to the scarce documentation [62]. In addition, the empirical and analytical results of this study could be a source of carbon credits through Payment for Environmental Services for small-scale homestead forests owners or households in developing countries, such as Bangladesh, while generating livelihood options and biodiversity conservation [61, 62].

Availability of data and materials

All data generated or analysed during this study are included in this article.


  1. 1.

    The second lowest administrative and management unit at the field level of Bangladesh Forest Department, BFD.

  2. 2.

    The second lowest administrative and management unit at the local level of the governing system of Bangladesh.



Above ground biomass


Below ground biomass


Basal area


Bulk density




Diameter at breast height


Loss of ignition


Relative frequency


Relative density


Soil organic carbon


Total biomass


Tree outside forest


  1. 1.

    IPCC: Climate change. 2014. Mitigation of climate change. Contribution of working group III to the fifth assessment report of the intergovernmental panel on climate change. 2014. Accessed 20 Mar 2020.

  2. 2.

    IPCC: Summary for policymakers 2018. In: Zhai P, Pörtner HO, Roberts D, Skea J, Shukla PR, Pirani A, Moufouma-Okia W, Pean C, Pidcock R, Connors S, Matthews JBR, Chen Y, Zhou X, Gomis MI, Lonnoy E, Maycock T, Tignor M, Waterfield T, editors. Global Warming of 1.5 °C. An IPCC Special Report on the Impacts of Global Warming of 1.5 °C above Pre-industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty [V. Masson-Delmotte. World Meteorological Organization, Geneva: Switzerland; 2008. p. 32

  3. 3.

    UNFCCC. Adoption of the paris agreement. Geneva: Switzerland; 2015. p. 32.

    Google Scholar 

  4. 4.

    Adame MF, Kauffman JB, Medina I, Gamboa JN, Torres O, Caamal JP, Reza M, Herrera-Silveira JA. Carbon stocks of tropical coastal wetlands within the karstic landscape of the Mexican Caribbean. PLoS ONE. 2013;8:565–9.

    Article  CAS  Google Scholar 

  5. 5.

    Dixon RK, Brown S, Houghton RA, Solomon AM, Trexler MC, Wisniewski J. Carbon pools and flux of global forest ecosystems. Science. 1994;263:185–90.

    CAS  Article  Google Scholar 

  6. 6.

    Mitchard ETA. The tropical forest carbon cycle and climate change. Nature. 2018;559:527–34.

    CAS  Article  Google Scholar 

  7. 7.

    Wilson EO, Peter FM, Raven PH. Our diminishing tropical forests. In: Wilson EO, Peter FM, editors. Biodiversity. Washington, DC: National Academy Press; 1988.

    Google Scholar 

  8. 8.

    Ashton MS, Tyrrell ML, Spalding D, Gentry B. Managing forest carbon in a changing climate. The Netherlands: Springer Science & Business Media; 2012.

    Book  Google Scholar 

  9. 9.

    Sarkar MSK, Sadeka S, Sikdar MMH, Badiuzzaman. Energy Consumption and CO2 Emission in Bangladesh: trends and policy implications. Asia Pac J Energy Environ. 2015;2(3):175–182

  10. 10.

    IEA (International Energy Agency): (2020). Accessed 22 Jun 2020.

  11. 11.

    Mani M, Bandyopadhyay S, Chonabayashi S, Markandya A, Mosier T. South Asia’s hotspots: the impact of temperature and precipitation changes on living standards. 2018; Accessed 22 Jun 2020

  12. 12.

    Alamgir M, Turton SM. Climate change and organic carbon storage in Bangladesh forests. In: Tuteja N, Gill SS, editors. Climate change and plant abiotic stress tolerance. Germany: Wiley; 2014.

    Google Scholar 

  13. 13.

    BFD (Bangladesh Forest Department): (2020). Accessed 23 Mar 2020.

  14. 14.

    Hansen MC, Potapov PV, Moore R, Hancher M, Turubanova SA, Tyukavina A, et al. High-resolution global maps of 21st-century forest cover change. Science. 2013;342:850–3.

    CAS  Article  Google Scholar 

  15. 15.

    Islam M, Deb GP, Rahman M. Forest fragmentation reduced carbon storage in a moist tropical forest in Bangladesh: implications for policy development. Land Use Policy. 2017;65:15–25.

    Article  Google Scholar 

  16. 16.

    Mukul SA. Biodiversity conservation and ecosystem functions of indigenous agroforestry systems: case study from three tribal communities in and around Lawachara National Park. In: Chowdhury MSH, editor. Forest conservation in protected areas of Bangladesh: policy and community development perspective. Springer: Switzerland; 2014. p. 171–9.

    Google Scholar 

  17. 17.

    Jaman MS, Hossain MF, Shariful I, J Helal MG, Jamil M, Mizanur R. Quantification of carbon stock and tree diversity of Homegardens in Rangpur District, Bangladesh. Int J Agric For. 2016;6:169–80. Available from:

  18. 18.

    Nath TK, Aziz N, Inoue M. Contribution of homestead forests to rural economy and climate change mitigation: a study from the ecologically critical area of Cox’s Bazar—Teknaf Peninsula, Bangladesh. Small-scale For. 2015;14:1–18.

    Article  Google Scholar 

  19. 19.

    Singh CD. Valuation and evaluation of trees outside of the forest (TOF) of Bangladesh. Draft paper for a regional study for Asia and Pacific in contribution to the forest resource assessment (FRA). Rome: FAO; 2000.

  20. 20.

    Kabir ME, Webb EL. Can homegarden conserve biodiversity in Bangladesh? Biotropica. 2008;40(1):95–103.

    Google Scholar 

  21. 21.

    Motiur RM, Furukawa Y, Kawata I, Rahman MM, Alam M. Role of homestead forests in household economy and factors affecting forest production: a case study in southwest Bangladesh. J For Res. 2006;11:89–97.

    Article  Google Scholar 

  22. 22.

    Mizanur Rahman M, Nabiul Islam Khan M, Fazlul Hoque AK, Ahmed I. Carbon stock in the Sundarbans mangrove forest: spatial variations in vegetation types and salinity zones. Wetl Ecol Manag. 2015;23:269–83.

    Article  Google Scholar 

  23. 23.

    Sohel MSI, Alamgir M, Akhter S, Rahman M. Carbon storage in a bamboo (Bambusa vulgaris) plantation in the degraded tropical forests: implications for policy development. Land Use Policy. 2015;49:142–51.

    Article  Google Scholar 

  24. 24.

    Kamruzzaman M, Ahmed S, Paul S, Rahman MM, Osawa A. Stand structure and carbon storage in the oligohaline zone of the Sundarbans mangrove forest, Bangladesh. Forest Sci Technol. 2018;14:23–8.

    Article  Google Scholar 

  25. 25.

    Khan MNI, Islam MR, Rahman A, Azad MS, Mollick AS, Kamruzzaman M, et al. Allometric relationships of stand level carbon stocks to basal area, tree height and wood density of nine tree species in Bangladesh. Glob Ecol Conserv. 2020;22:e01025.

    Article  Google Scholar 

  26. 26.

    Alamgir M, Al-Amin M. Organic carbon storage in trees within different Geopositions of Chittagong (South) Forest Division, Bangladesh. J For Res. 2007;18:174–80.

    CAS  Article  Google Scholar 

  27. 27.

    Mandal RA, Jha PK, Dutta IC, Thapa U, Karmacharya SB. Carbon sequestration in tropical and subtropical plant species in collaborative and community forests of Nepal. Adv Ecol. 2016.

    Article  Google Scholar 

  28. 28.

    Ullah MR, Banik GR, Banik R. Developing allometric models for carbon stock estimation in eighteen year old plantation forests of Bangladesh. Jacobs J Microbiol Patholog. 2014;1:1–8.

    Google Scholar 

  29. 29.

    Day M, Baldauf C, Rutishauser E, Sunderland TCH. Relationships between tree species diversity and above-ground biomass in Central African rain forests: implications for REDD. Environ Conserv. 2013;41:64–72.

    Article  Google Scholar 

  30. 30.

    Shen Y, Yu S, Lian J, Shen H, Cao H, Lu H, Ye W. Tree aboveground carbon storage correlates with environmental gradients and functional diversity in a tropical forest. Sci Rep. 2016;6:1–10.

    CAS  Article  Google Scholar 

  31. 31.

    Dondini M, Hastings A, Saiz G, Jones M, Smith P. The potential of Miscanthus to sequester carbon in soils: comparing field measurements in Carlow, Ireland to model predictions. Glob Change Biol Bioenergy. 2009;1:413–25.

    CAS  Article  Google Scholar 

  32. 32.

    Ostertag R, Marín-Spiotta E, Silver WL, Schulten J. Litterfall and decomposition in relation to soil carbon pools along a secondary forest chronosequence in Puerto Rico. Ecosystems. 2008;11:701–14.

    CAS  Article  Google Scholar 

  33. 33.

    Houghton RA. Balancing the global carbon budget. Annu Rev Earth Planet Sci. 2007;35:313–47.

    CAS  Article  Google Scholar 

  34. 34.

    Ullah MR, Al-Amin M. Above- and below-ground carbon stock estimation in a natural forest of Bangladesh. J For Sci. 2012;58:372–9.

    Article  Google Scholar 

  35. 35.

    Danquah JA, Appiah M, Pappinen A. Effect of African mahogany species on soil chemical properties in degraded dry semi-deciduous forest ecosystems in Ghana. Int J Agric Biol. 2012;14:321–8.

    CAS  Google Scholar 

  36. 36.

    Sayer EJ, Lopez-Sangil L, Crawford JA, Bréchet LM, Birkett AJ, Baxendale C, et al. Tropical forest soil carbon stocks do not increase despite 15 years of doubled litter inputs. Sci Rep. 2019;9:1–9.

    Article  CAS  Google Scholar 

  37. 37.

    Miao R, Ma J, Liu Y, Liu Y, Yang Z, Guo M. Variability of aboveground litter inputs alters soil carbon and nitrogen in a coniferous-broadleaf mixed forest of central China. Forests. 2019;10:188.

    Article  Google Scholar 

  38. 38.

    Krishna MP, Mohan M. Litter decomposition in forest ecosystems: a review. Joint Center on Global Change and Earth System Science of the University of Maryland and Beijing Normal University. Energy Ecol Environ. 2017;2:236–49.

    Article  Google Scholar 

  39. 39.

    Akhtaruzzaman M, Osman KT, Sirajul Haque SM. Soil properties in two forest sites in Cox’s Bazar, Bangladesh. J For Environ Sci. 2015;31:280–7.

    Google Scholar 

  40. 40.

    Zaman MA, Osman KT, Sirajul Haque SM. Comparative study of some soil properties in forested and deforested areas in Cox’s Bazar and Rangamati Districts, Bangladesh. J For Res. 2010;21:319–22.

    CAS  Article  Google Scholar 

  41. 41.

    Barua SK, Haque SMS. Soil characteristics and carbon sequestration potentials of vegetation in degraded hills of Chittagong, Bangladesh. Land Degrad Develop. 2013;24:63–71.

    Article  Google Scholar 

  42. 42.

    Kibria MG, Saha N. Analysis of existing agroforestry practices in Madhupur Sal forest: an assessment based on ecological and economic perspectives. J For Res. 2011;22:533–42.

    CAS  Article  Google Scholar 

  43. 43.

    Chen S, Wang W, Xu W, Wang Y, Wan H, Chen D, et al. Plant diversity enhances productivity and soil carbon storage. Proc Natl Acad Sci USA. 2018;115:4027–32.

    CAS  Article  Google Scholar 

  44. 44.

    Isaac SR, Nair MA. Biodegradation of leaf litter in the warm humid tropics of Kerala, India. Soil Biol Biochem. 2005;37:1656–64.

    CAS  Article  Google Scholar 

  45. 45.

    Nair PKR, Nair VD, Kumar BM, Haile SG. Soil carbon sequestration in tropical agroforestry systems: a feasibility appraisal. Environ Sci Policy. 2009;12:1099–111.

    CAS  Article  Google Scholar 

  46. 46.

    Saha SK, Nair PKR, Nair VD, Kumar BM. Soil carbon stock in relation to plant diversity of homegardens in Kerala, India. Agroforest Syst. 2009;76:53–65.

    Article  Google Scholar 

  47. 47.

    Lal R. Forest soils and carbon sequestration. For Ecol Manage. 2005;220:242–58.

    Article  Google Scholar 

  48. 48.

    Lal R. Soil carbon stocks under present and future climate with specific reference to European ecoregions. Nutr Cycl Agroecosyst. 2008;81:113–27.

    Article  Google Scholar 

  49. 49.

    Nair PKR, Kumar BM, Nair VD. Agroforestry as a strategy for carbon sequestration. J Plant Nutr Soil Sci. 2009;172:10–23.

    CAS  Article  Google Scholar 

  50. 50.

    Homann PS, Bormann BT, Boyle JR, Darbyshire RL, Bigley R. Soil C and N minimum detectable changes and treatment differences in a multi-treatment forest experiment. For Ecol Manag. 2008;255:1724–34.

    Article  Google Scholar 

  51. 51.

    Casals P, Romero J, Rusch GM, Ibrahim M. Soil organic C and nutrient contents under trees with different functional characteristics in seasonally dry tropical silvopastures. Plant Soil. 2014;374:643–59.

    CAS  Article  Google Scholar 

  52. 52.

    Islam M, Dey A, Rahman M. Effect of tree diversity on soil organic carbon content in the homegarden agroforestry system of north-eastern Bangladesh. Small-scale For. 2015;14:91–101.

    Article  Google Scholar 

  53. 53.

    Gardner TA, Burgess ND, Aguilar-Amuchastegui N, Barlow J, Berenguer E, Clements T, et al. A framework for integrating biodiversity concerns into national REDD+ programmes. Biol Conserv. 2012;154:61–71.

    Article  Google Scholar 

  54. 54.

    Takimoto A, Nair PKR, Nair VD. Carbon stock and sequestration potential of traditional and improved agroforestry systems in the West African Sahel. Agric Ecosyst Environ. 2008;125:159–66.

    CAS  Article  Google Scholar 

  55. 55.

    Paoli GD, Wells PL, Meijaard E, Struebig MJ. et al. Biodiversity conservation in the REDD. Carbon Balance Manag. 2010;5:7.

  56. 56.

    Grainger A, Boucher DH, Frumhoff PC, Laurance WF, Lovejoy T, McNeely J, Niekisch M, Raven P, Sodhi NS, Venter O, Pimm SL. Biodiversity and REDD at Copenhagen. Curr Biol. 2009;19:R974–6.

    CAS  Article  Google Scholar 

  57. 57.

    UN-REDD. Carbon and biodiversity relationships in tropical forests. 2010. Multiple Benefits Series 4. UN-REDD programme.

  58. 58.

    Strassburg BBN, Kelly A, Balmford A, Davies RG, Gibbs HK, Lovett A, Miles L, Orme CDL, Price J, Turner RK, Rodrigues ASL. Global congruence of carbon storage and biodiversity in terrestrial ecosystems. Conserv Lett. 2010;3:98–105.

    Article  Google Scholar 

  59. 59.

    Baccini A, Goetz SJ, Walker WS, Laporte N, Sun M, Sulla-Menashe D, Hackler J, Beck PSA, Dubayah R, Friedl MA, Samanta S, Houghton RA. Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nat Clim Change. 2010;2:182–5.

    Article  CAS  Google Scholar 

  60. 60.

    Chave J, Réjou-Méchain M, Búrquez A, Chidumayo E, Colgan MS, Delitti WBC, et al. Improved allometric models to estimate the aboveground biomass of tropical trees. Glob Chang Biol. 2014;20:3177–90.

    Article  Google Scholar 

  61. 61.

    Nair PKR. Carbon sequestration studies in agroforestry systems: a reality check. Agroforest Syst. 2012;86(2):243–53.

    Article  Google Scholar 

  62. 62.

    Jashimuddin M, Inoue M. Management of village common forests in the Chittagong hill tracts of Bangladesh: historical background and current issues in terms of sustainability. Open J For. 2012;2(3):121–37.

    Google Scholar 

  63. 63.

    Majlis AB, Islam M, Khasru MA, Ahsan MK. Protectred to open basin depositional system: an appraisal for the quaternary evolution of the Moheshkhali-Kutubdia Coastal Plain, Bangladesh, Abstract volume. National seminar on Bangladesh Coast: Geology, Hazards and Resources. Dhaka, Bangladesh. 2011.

  64. 64.

    Islam MA, Maitra MK, Baquee A, Majlis K, Murshed S, Rahman S. Spatial changes of land use/land cover of Moheshkhali Island, Bangladesh: a fact finding approach by remote sensing analysis. Dhaka Univ J Earth Environ Sci. 2012;2:43–54.

    Google Scholar 

  65. 65.

    Upazila Parishad office. Maheshkhali Upazila, Cox´s Bazar District: Bangladesh; 2019.

  66. 66.

    Bangladesh Meteorological Department, Cox’s Bazar Region, Chittagong Division: Bangladesh; 2019.

  67. 67.

    Islam MA, Majlis ABK, Bazlar R. Changing face of Bangladesh Coast. Abstract volume. National seminar on Bangladesh Coast: Geology, Hazards and Resources. Dhaka, Bangladesh; 2011.

  68. 68.

    Papry RI. Status of coastal plantation in Chittagong coastal forest division. IOSR J Environ Sci Toxicol Food Technol. 2014;8:79–83.

    Article  Google Scholar 

  69. 69.

    DoE (Department of Environment). GIS and Cartographic Services—Final Report, Pre-Investment Facility Study: Coastal and Wetland Biodiversity Management Project (Project BGD/94/G41). Dhaka, Bangladesh; 1999.

  70. 70.

    Maheshkhali Forest Range Office. Chittagong Coastal Forest Division, Chittagong, Bangladesh Forest Department: Bangladesh; 2019.

  71. 71.

    Union Parishad Office. Bara Maheshkhali, Chatto Maheshkhali, Gorokghata Unions, Cox’s Bazar District, Chittagong Division: Bangladesh; 2019.

  72. 72.

    UN (United Nations). Designing household survey samples: practical guidelines. Publication No. ST/ESA/STAT/SER.F/98. Department of Economic and Social Affairs, Statistics Division: New York; 2005

  73. 73.

    Ensslin A, Rutten G, Pommer U, Zimmermann R, Hemp A, Fischer M. Effects of elevation and land use on the biomass of trees, shrubs and herbs at Mount Kilimanjaro. Ecosphere. 2015;6:art45.

    Article  Google Scholar 

  74. 74.

    Blake GR. 1965. Bulk density. In: Black CA, Evans DD, White JL, Ensminger IE, Clark FE, editors. Methods of soil analysis. Part 1. American Society of Agronomy, Wisconsin, USA: Inc. Publisher, Madison; pp. 894–895.

  75. 75.

    Hairiah K. Measuring carbon stocks: across land use systems: a manual. Published in close cooperation with Brawijaya University and ICALRRD (Indonesian Center for Agricultural Land Resources Research and Development); 2011.

  76. 76.

    Frangi JL, Lugo AE. Ecosystem dynamics of a subtropical floodplain forest. Ecol Monogr. 1985;55(3):351–69.

    Article  Google Scholar 

  77. 77.

    Issa S, Dahy B, Ksiksi T, Saleous N. Development of a new allometric equation correlated wth variables for the assessment of date palm biomass. Conference paper. UAE University, College of Science: United Arab Emirates; 2018.

  78. 78.

    MacDicken KG. A guide to monitoring carbon storage in forestry and agro-forestry projects. Arlington, USA: Winrock International; 1997.

    Google Scholar 

  79. 79.

    Brown S. Estimating biomass and biomass change of tropical forests: a primer. FAO Forestry Paper 134. Rome: FAO; 1997.

  80. 80.

    Sattar MA, Bhattacharje DK, Kabir MF. Physical and mechanical properties and uses of timbers of Bangladesh. Chittagong, Bangladesh: Bangladesh Forest Research Institute; 1999.

    Google Scholar 

  81. 81.

    Chave J, Coomes D, Jansen S, Lewis SL, Swenson NG, Zanne AE. Towards a worldwide wood economics spectrum. Ecol Lett. 2009;12(4):351–66.

    Article  Google Scholar 

  82. 82.

    Zanne AE, Lopez-Gonzalez G, Coomes DA, Ilic J, Jansen S, Lewis SL, et al. Data from: towards a worldwide wood economics spectrum, Dryad digital repository. 2009. Global Wood Density Database.

  83. 83.

    Shukla RS, Chandel PS. Plant ecology and soil science. 9th ed. India: New Delhi; 2000.

    Google Scholar 

  84. 84.

    Ball DF. Loss-on-ignition as an estimate of organic matter and organic carbon in noncalcareous soils. J Soil Sci. 1964;15:84–92.

    CAS  Article  Google Scholar 

  85. 85.

    Pearson T, Walker S, Brown S. Sourcebook for land use, land-use change and forestry projects 29 (Winrock International and the BioCarbon Fund of the World Bank). 2005.

  86. 86.

    Davidson DA. Soil physical measurement and interpretation for land evaluation. Aclep, Land Evaluation Program. 1986;2:1–16. Acessed 25 Mar 2020.

  87. 87.

    Pearson TRH, Brown SL, Birdsey RA. Measurement Guidelines for the Sequestration of Forest Carbon. General Technical Report NRS-18. USDA FOREST SERVICE, United States Department of Agriculture; 2007.

  88. 88.

    Coleman DC. Soil carbon balance in a successional grassland. Oikos. 1973;24:195–9.

    CAS  Article  Google Scholar 

  89. 89.

    Margalef R. Information theory in ecology. General Syst. 1958;3:36–71.

    Google Scholar 

  90. 90.

    Michael P. Ecological methods for field and laboratory investigation. New Delhi, India: Tata Mc Graw Hill; 1990.

    Google Scholar 

  91. 91.

    Dallmeier F, Kabel M, Rice R. Methods for long-term biodiversity inventory plots in protected tropical forests. In: Dallmeier F, editor. Long-term monitoring of biological diversity in tropical forest areas: methods for establishment and inventory of permanent plots. Paris, France: UNESCO; 1992. p. 11–46.

    Google Scholar 

  92. 92.

    Hasan MA, Abubakar IA, Rahman SM, Aina YA, Chowdhury MMI, Khondaker AN. The synergy between climate change policies and national development goals: implications for sustainability. J Clean Prod. 2019.

    Article  Google Scholar 

  93. 93.

    Poorter L, van der Sande MT, Thompson J, Arets EJMM, Alarcón A, Álvarez-Sánchez J, et al. Diversity enhances carbon storage in tropical forests. Glob Ecol Biogeogr. 2015;24:1314–28.

    Article  Google Scholar 

  94. 94.

    Rahman MM, Kabir ME, JahirUddinAkon ASM, Ando K. High carbon stocks in roadside plantations under participatory management in Bangladesh. Glob Ecol Conserv (Elsevier BV). 2015;3:412–23.

    Article  Google Scholar 

  95. 95.

    Ordóñez JAB, de Jong BHJ, Garcia-Oliva F, Avina FL, Perez JV, Guerrero G, et al. Carbon content in vegetation, litter, and soil under 10 different land-use and land-cover classes in the Central Highlands of Michoacan, Mexico. For Ecol Manage. 2008;255:2074–84.

    Article  Google Scholar 

  96. 96.

    Baul TK, Datta D, Alam A. A comparative study on household level energy consumption and related emissions from renewable (biomass) and non-renewable energy sources in Bangladesh. Energy Pol. 2018;114:598–608.

    Article  Google Scholar 

  97. 97.

    Hossen MM, Rahman AHMS, Kabir AS, Hasan MMF, Ahmed S. Systematic assessment of the availability and utilization potential of biomass in Bangladesh. Renew Sustain Energy Rev. 2017;67:94–105.

    Article  Google Scholar 

  98. 98.

    Schleuß P-M, Heitkamp F, Leuschner C, Fender A-C, Jungkunst HF. Higher subsoil carbon storage in species-rich than species-poor temperate forests. Environ Res Lett. 2014;9:1–10.

    CAS  Article  Google Scholar 

  99. 99.

    Islam KR, Weil RR. Land use effects on soil quality in a tropical forest ecosystem of Bangladesh. Agric Ecosyst Environ. 2000;79:9–16.

    Article  Google Scholar 

  100. 100.

    Jamaludheen V, Kumar BM. Litter of multipurpose trees in Kerala, India: variations in the amount, quality, decay rates and release of nutrients. For Ecol Manage. 1999;115:1–11.

    Article  Google Scholar 

  101. 101.

    Hossain MK, Alim A, Hossen S, et al. Diversity and conservation status of tree species in Hazarikhil Wildlife Sanctuary (HWS) of Chittagong, Bangladesh. Geol Ecol Landscapes. 2020;4:298–305.

    Article  Google Scholar 

  102. 102.

    Maraseni TN, Neupane PR, Lopez Casero F, Cadman T. An assessment of the impacts of the REDD+ pilot project on community forests user groups (CFUGs) and their community forests in Nepal. J Environ Manage. 2014;136:37–46.

    CAS  Article  Google Scholar 

  103. 103.

    Poudyal BH, Maraseni TK, Cockfield G. Impacts of forest management on tree species richness and composition: assessment of forest management regimes in Tarai landscape, Nepal. Appl Geogr. 2019;111:102078.

    Article  Google Scholar 

  104. 104.

    Luintel H, Scheller RM, Bluffstone RA. Assessments of biodiversity, carbon, and their relationships in Nepalese forest commons: implications for global climate initiatives. For Sci. 2018;64:418–28.

    Article  Google Scholar 

  105. 105.

    Tubiello FN, Cóndor-Golec RD, Salvatore M, Piersante A, Federici S, Ferrara A, Rossi S, et al. Estimating greenhouse gas emissions in agriculture, a manual to address data requirements for developing countries. Rome: Food and Agriculture Organization; 2015.

    Google Scholar 

  106. 106.

    Mehling MA, Metcalf GE, Stavins RN. Linking climate policies to advance global mitigation. Science. 2018;359:997–8.

    CAS  Article  Google Scholar 

Download references


Authors would like to acknowledge the Research Cell and Publication office of University of Chittagong, Bangladesh for providing research grant. Authors are highly thankful to the homestead forest owners for the permission to survey their forests and giving us valuable time during survey. In addition, concerned officers of the local government and Forest Department offices in the study area are greatly acknowledged for providing us with valuable information related to the study.


This research work was funded by the Research Cell and Publication office of University of Chittagong, Bangladesh.

Author information




All authors have contributed significantly to produce this manuscript. TKB has made research conception and design, analyses, and interpretation of data and been involved in preparing and revising the manuscript. AC has been in field survey and organizing field data and involved in analyses of soil and litter data in the lab. RN has been involved in research conception, design of data acquisition and been involved in revising the manuscript. MM has been involved in drafting or revising the manuscript critically. AK has made substantial contribution with editing and revising the manuscript rigorously. TS has been involved in analyses of soil and litter data in the lab. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Tarit Kumar Baul.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Baul, T.K., Chakraborty, A., Nandi, R. et al. Effects of tree species diversity and stand structure on carbon stocks of homestead forests in Maheshkhali Island, Southern Bangladesh. Carbon Balance Manage 16, 11 (2021).

Download citation


  • Biomass
  • Biodiversity
  • Carbon stock
  • Homestead forests
  • Litterfall
  • REDD
  • Soil
  • Tree species