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

Background 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. Results 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%. Conclusions 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.


Background
A gradual increase in the global emissions of carbon dioxide (CO 2 ) 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 CO 2 as they store one fourth of the global terrestrial carbon [4][5][6]. In addition, tropical

Open Access
Carbon Balance and Management *Correspondence: tarit.ifes@cu.ac.bd 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, CO 2 emissions over the time are rapidly increasing, for example, 609% in 2017 (78 Mt of CO 2 ) compared to that of 1990 (11 Mt of CO 2 ) 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 landuse 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 semievergreen 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.
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 subnational 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.

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 km 2 , 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.
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, 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 Unions 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 m 2 (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 cm 3 ) 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 (m 2 ha −1 ) were estimated (Eqs. [5][6][7]. Mean values of tree biomass, density, and BA were compared among three different homestead forest sites.

Laboratory analysis and estimation of carbon of litterfall and mineral soil and bulk density (BD)
To estimate soil organic carbon (SOC), soils were ovendried 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  Dallmeier et al. [91] 17 Relative density, RD(%) = S n × 100 Dallmeier et al. [91] 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][16][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 (m 2 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 (m 2 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).
Mean stand density and BA of homestead forests were 601 individuals ha −1 and 27 m 2 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).
Among 52 tree species found in the homestead forests, the numbers of species in the hillside, beachside, and  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).

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.
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).

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). 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 (m 2 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.

Discussion
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 CO 2 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 m 2 ha −1 ) increased biomass carbon stocks in roadside plantations.
In this study, with greater DBH and height, aboveand 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    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 CO 2 [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].