Can Mangrove Silviculture Be Carbon Neutral?

: Matang Mangrove Forest Reserve (MMFR) in peninsular Malaysia has been managed for pole and charcoal production from Rhizophora stands with a 30-year rotation cycle since 1902. The aim of this study is to estimate the carbon budget of the MMFR by considering the carbon stock of the forest, evaluated from remote sensing data (Landsat TM and ETM+, JERS-1 SAR, ALOS PALSAR, ALOS-2 PALSAR-2, SRTM, TANDEM-X, and WorldView-2) for aboveground carbon and ﬁeld data for belowground carbon. This was investigated in combination with the emissions from the silvicultural activities in the production chain, plus the distribution and consumer-related activities covering the supply chain, estimated with appropriate emission factors. The aboveground biomass carbon stock of the productive forest was of 1.4 TgC, while for the protective forest (not used for silviculture) it was at least equal to 1.2 TgC. The total soil carbon of ca. 32 TgC shows the potential of the MMFR as a carbon sink. However, the commercial exploitation of mangroves also generates greenhouse gasses with an estimate of nearly 152.80 Mg C ha − 1 during charcoal production and up to 0.53 Mg C ha − 1 during pole production, for a total emission of 1.8 TgC. Consequently, if the productive forest alone is considered, then the carbon budget is negative, and the ongoing silvicultural management seems to be an unsustainable practice that needs a reduction in the exploited area of at least 20% to achieve carbon neutrality. However, even with the current management, and considering the protective forest together with the productive zones, the MMFR carbon budget is slightly positive, thus showing the importance of mangrove conservation as part of the management for the preservation of the carbon stock.


Introduction
In recent years the importance of sustainable conservation and management of natural ecosystems has grown exponentially due to increased visibility of the negative impacts of greenhouse gas emissions and climate change [1,2]. One of the major effects is the increased vulnerability of coastal systems to sea-level rise and storm surges [3][4][5]. In tropical and subtropical countries this has increased the importance of the conservation and adequate management of mangrove forests, an essential ecosystem for coastal protection and carbon sequestration [6][7][8]. Recent literature on mangroves has been focused on finding strategies to stabilize and offset carbon emissions [9][10][11], including through programs such as Reducing Emissions from Deforestation and forest Degradation'(REDD+) that aims at maintaining economic growth while maintaining carbon neutrality in the use of ecosystem resources [12].
Estimates of global mangrove carbon stocks have been improved over the years [13][14][15][16]. Due to the peculiarity of having forest stands of known age, the MMFR has been the centre of many studies: Ong et al. [28] were the first to investigate the aboveground biomass of 5-, 10-, 15-, and 25-year-old trees using allometric equations (ranging from 16 to 300 Mg ha −1 ), followed by Putz and Chan [29] (62-462 Mg ha −1 ). More recently Goessens et al. [30] and Otero et al. [31] estimated the aboveground biomass to range (on average) from 216 to 372 Mg ha −1 and 217 to 238 Mg ha −1 , respectively. Subsequently, by combining site-and species-specific biomass observations from the MMFR [32] with those from Thailand and Indonesia [33], Komiyama et al. [34] provided a review on the allometric equations that have proved essential for mangrove research. Recently, Adame et al. [35] claimed that carbon loss due to the clear-felling activity at the MMFR recovers to carbon gain over a period of time, with the greatest recovery rates (9.5 Mg C ha −1 yr −1 ) being in the first ten years. This can also be seen in changes in forest cover in remote sensing data [36][37][38] with a regeneration time of 5.9 ± 2.7 years [36] associated with aboveground biomass accumulation rates of 8. 36 Mg ha −1 [37].
Although MMFR management is keen to support sustainable practices and ensure the ecological functions (such as habitat and coastal protection) are incessant [26], there is no sufficient information on the impact of the carbon emissions from the charcoal production yet. In terms of timber extraction [39,40], the charcoal production can have negative ecological consequences due to the loss of the natural and/or planted vegetation. In order to determine the environmental impact of charcoal production in the MMFR, it is imperative to quantify emissions released by different production activities. For mangroves, the studies conducted by Lacaux et al. [41] in the African tropics, Pennise et al. [42] in Brazil and Kenya, Kridiborworn et al. [43] in Thailand, and Bailis et al. [44] in Brazil, estimated the carbon emissions from wood burning inside charcoal kilns can be up to 3027 g CO2 per Kg of charcoal depending on the kiln type; however, they did not consider all carbon emission sources related to tree cutting and the transport of logs and charcoal. Until now, we did not find any scientific records on the estimated total carbon impact of charcoal production for any mangrove location in the world, including Malaysia.
The present study, centred on MMFR, aimed to (i) estimate carbon stock in the entire MMFR based on remote sensing data; (ii) estimate carbon emissions of the entire pole and Although MMFR management is keen to support sustainable practices and ensure the ecological functions (such as habitat and coastal protection) are incessant [26], there is no sufficient information on the impact of the carbon emissions from the charcoal production yet. In terms of timber extraction [39,40], the charcoal production can have negative ecological consequences due to the loss of the natural and/or planted vegetation. In order to determine the environmental impact of charcoal production in the MMFR, it is imperative to quantify emissions released by different production activities. For mangroves, the studies conducted by Lacaux et al. [41] in the African tropics, Pennise et al. [42] in Brazil and Kenya, Kridiborworn et al. [43] in Thailand, and Bailis et al. [44] in Brazil, estimated the carbon emissions from wood burning inside charcoal kilns can be up to 3027 g CO 2 per Kg of charcoal depending on the kiln type; however, they did not consider all carbon emission sources related to tree cutting and the transport of logs and charcoal. Until now, we did not find any scientific records on the estimated total carbon impact of charcoal production for any mangrove location in the world, including Malaysia.
The present study, centred on MMFR, aimed to (i) estimate carbon stock in the entire MMFR based on remote sensing data; (ii) estimate carbon emissions of the entire pole and charcoal production chains, i.e., including the production (silvicultural practices), distribution (over land and water), and consumer chain (end-user utilisation); and (iii) assess whether or not the area under management is carbon neutral.

Study Area
The Matang Mangrove Forest Reserve ( Figure 1) covers about 40,288 ha and has been monitored by the State Forestry Department of Perak for timber production since 1902 [25]. In 1906, it was declared a permanent forest reserve with the aim of charcoal and pole production [25]. The introduction of a new management zoning in the 2000s (management plan 2000-2009 [45]) divided the MMFR into four distinct categories: productive (74.8%), protective (17.4%), restrictive productive (6.8%), and unproductive forest (1%) [26,45,46]. While timber is extracted from the productive zone, both protective zones and restrictive zones are retained for erosion control, biodiversity preservation, preservation of the mangrove seed bank, etc. The productive zone is largely composed of mangrove stands dominated Remote Sens. 2022, 14, 2920 4 of 17 by Rhizophora species of different ages, located in 108 compartments or coupes [26,30]. According to the (2010-2019) management plan, nearly 27,430 ha were allocated for timber production; of which, 8159 ha were for the first thinning, 7678 ha were for the second thinning, and 11,593 ha were for the final felling [26]. In total, 144 charcoal contractors are present and each one is authorized to exploit 2.2 ha y −1 for a kiln (maximum 3 kilns allowed for each contractor). Similarly, 70 pole contractors are present with each allowed to exploit 22.3 ha y −1 [26].

Primary Source Information
This study is a follow-up from the study by Lucas et al. [37] who utilized remote sensing imagery from the MMFR from 1988 to 2016 to evaluate the structural characteristics of the mangrove forest. These data were integrated with data from Goessens et al. [30] and Adame et al. [35] to evaluate the total carbon stock of the MMFR.
Furthermore, in order to evaluate the carbon emissions from the charcoal and pole production activities, a systematic review of primary sources (books, theses, web pages, and articles) was conducted to collect the information on carbon analytical techniques, procedures estimating different emission factors, etc. In this context, a wealth of information was extracted from the works of Satyanarayana et al. [27] along with the 2010-2019 MMFR management plan [26].

Vegetation Carbon Stock
To evaluate the carbon stock for the aboveground biomass (AGB) the most recent remote sensing biomass data on the productive zones in the MMFR [37] were analysed. Lucas et al. [37] estimated forest age on an annual basis for the period 1988 to 2016 by combining the time series of the Landsat-derived Normalised Difference Moisture Index (NDMI) and Japanese L-band Synthetic Aperture Radar (SAR) data. The NDMI was further used to retrieve canopy cover (%). Furthermore, interferometric Shuttle Radar Topographic Mission (SRTM) X/C-band (2000), TanDEM-X-band (2010-2016), and stereo WorldView-2 stereo (2016) data were used to estimate canopy height (CH), from which above ground biomass (AGB, Mg ha −1 ), consisting of stem, branch, leaf biomass, and prop roots, was derived using pre-established allometry [16,34,[47][48][49].
Starting from the AGB data for the year 2016 obtained from the above methods, the total AGB values (Mg) were calculated for the different zones of the MMFR. To achieve these results the age map and the biomass map were combined using QGIS to be able to differentiate the different age categories within the MMFR. The AGB data of each forest patch were subsequently multiplied by the respective area to obtain values in Mg. Finally, the biomass data were converted to carbon data with a conversion factor of 0.47 [50].
To estimate the belowground biomass (BGB), i.e., the carbon pool associated with roots, the data from Goessens et al. [30] were analysed for forest age of 15, 20, and 30 years and aboveground and belowground carbon pools were computed using the equations by Komiyama et al. [34]. The biomass values were then converted to carbon using a conversion factor of 0.39 [50]. The BGB carbon values were equivalent to ca. 35% of the aboveground values for all ages considered and hence a conversion factor of 0.35 was used to integrate the BGB carbon pool with the AGB carbon from the remote sensing data. Finally, after summing the AGB and BGB carbon values, the biomass carbon stock in the productive zones within the MMFR was mapped in QGIS according to the age of the different forest coupes.
To calculate the total carbon stock of the MMFR according to the forest zones (productive, protective, and restrictive productive), the biomass values of each single forest coupe were summed. For the protective forest zones, i.e., areas that are not used for silviculture, the AGB data were not available for most forest compartments due to the limits of the canopy height model above a certain height/age of the tree and the higher heterogeneity of the forest structure in these areas. Therefore, the biomass of a 30-year-old stand was used for the calculation of the carbon stock of the protective forest zones and the results must be considered as an underestimation of the real biomass carbon stock.

Soil Carbon
Soil carbon data were estimated from Adame et al. [35] who measured the carbon content in 1 m soil cores at different forest ages (5,15,30, and 40 years and a clear-felled site). Mean carbon values were extrapolated for ages ranging 0 to 30 years and each forest patch of a certain age was associated with its soil carbon value, using the same procedure described above for the AGB. As for the biomass carbon, a map for soil carbon was generated for the year 2016.

Carbon Emission through Charcoal and Pole Production
To estimate the carbon emissions from the silvicultural practices in the Matang Mangrove Forest Reserve, the main activities involved in mangrove charcoal and pole production systems were identified and divided into direct (i.e., activities with a direct contact with the product) and indirect tasks (i.e., out of the activity chain).
The direct tasks in charcoal production begin with transportation of tree cutters from the charcoal factory to the allocated clear-felling areas by boat. While the cutters clear fell an area of 2.2 ha of vegetation in 14 to 24 days (using chainsaw), the boat drivers transport the logs (cut into 1.6 m logs) daily to the charcoal factory. The logs are arranged inside the kilns where greenwood is converted into charcoal in 40-45 days. Nearly 80% of the charcoal from the MMFR is sent to the exporting units (by lorry) for shipment to Japan; the remaining 20% is distributed to the local markets in peninsular Malaysia [27]. Indirect tasks were represented by mode of transportation (moped or car) by the local workers to reach charcoal factories and the amount of charcoal burned for barbecuing in Japan.
In contrast to the charcoal system, mangrove pole production follows a simpler activity chain. Poles are used in peninsular Malaysia as piling material for house and road construction. Workers reach the allocated thinning areas by boat and cut down the Rhizophora trees (using an axe) by following a stick method of 1.2 m distance between the trees in 15 years and 1.8 m in 20-year-old stands. Poles are transferred to the charcoal factory by boat and from there to the local markets by lorry [27].
To estimate the carbon emissions per task/activity we referred to average emission factors that were selected based on types of equipment used and brand specifications (details can be found in the Supplementary Materials of this article). In this context, all exhaust gasses were considered and converted to carbon following the IPCC criteria [51,52]. In particular, nitrous oxide (N 2 O) and methane (CH 4 ), though not released in large quantities, were considered after conversion to CO 2 equivalents and subsequently to carbon values. Conversion in CO 2 equivalents was conducted using the conversion factor of 21 for CH 4 and of 289.3 for N 2 O and subsequently CO 2 equivalents were converted in carbon values using the conversion factor of 3.67 [51]. Emission factors for each activity, i.e., how much carbon is released during one iteration of the activity, were calculated by summing the emissions of each gas converted into carbon (Equation (1)).
where C () represents the conversion of the different gasses into carbon and n 1 , n 2 , n 3 , and n 4 represent the total emissions per each gas. The actual carbon emission per activity (gC) was calculated by multiplying the emission factor by the corresponding distances for travelling (for boat, truck, moped, and car) and quantities of greenwood for conversion into charcoal in the kilns, and quantities of charcoal for export and use for barbecue (Table 1). All these were calculated firstly for one iteration of each activity using Equations (2)- (9).
Remote Sens. 2022, 14, 2920 where a is the carbon emission factor value, b is the unit of power from the selected chainsaw, and c is the time of use in hours.
where a is the carbon emission factor value and b is the total distance in kilometres.
where a is the carbon emission factor value and b is the amount of charcoal in grams.
where a is the carbon emission factor value and b is the total distance in kilometres.
where a is the carbon emission factor value, b is the amount of charcoal transported in tonnes, and c is the total distance in kilometres.
where a is the carbon emission factor value and b is the total distance in kilometres.
where a is the carbon emission factor value and b is the total distance in kilometres.
where a is the carbon emission factor value and b is the total amount of charcoal in kilograms. Table 1. Equations used for the calculation of carbon emissions for each activity in the pole and charcoal production chain. The emission factors were computed using Equation (1).

Activity Equation
Tree cutting Equation (2) Logs transport by small boat Equation (3) Conversion of wood into charcoal in kiln Equation (4) Charcoal transport by truck Equation (5) Charcoal transport by Cargo ship for exportation Equation (6) Tree cutters transport by boat Equation (3) Workers and contractors transport by moped Equation (7) Workers and contractors transport by car Equation (8) Barbecue Equation (9) Due to uncertain details of the transportation used by different stakeholders in both charcoal and pole production chains, some carbon emission estimates were formed based on best professional judgment.
To calculate the total emissions for the charcoal and pole production chains, the results from the single activities were combined after multiplication by the number of iterations of each task needed to complete the production for the licensed annual area allocated to one Remote Sens. 2022, 14, 2920 7 of 17 contractor (i.e., 2.2 ha for charcoal and 22.4 ha for poles). Afterwards, the emissions were also computed for the total allocated area for charcoal and pole production reported in the latest management plan.
Additionally, carbon loss from the soil after clear-felling was also considered as a side effect of the charcoal production by referring to the findings by Adame et al. [35]. In this context, the carbon content of the upper first meter of soil from a 30-year-old stand was duly compared with a clear-felled site.

Comparison between Carbon Content and Carbon Emission
To evaluate whether or not the MMFR area under silvicultural management can be considered carbon neutral, the estimates of carbon sequestration and carbon emissions were duly computed. First, the emissions from pole and charcoal production were compared with the carbon stock of the productive forest and then the protective forest carbon stock was added to calculate the total carbon budget of the MMFR.
The comparison between carbon emissions and carbon stock was possible because we considered the emissions for a rotation cycle of 30 years based on the area information provided by the latest management plan [26] and we compared it with the carbon stock accumulated in the productive areas over 30 years.

Carbon Stock in the MMFR
To classify the areas based on the rotation cycle and silvicultural thinning practice at 15 and 20 years, respectively, the productive zone was divided into three categories: 0-15 years, 16-20 years, and 21-30 years (Table 2, Figure 2). The carbon stock ranged from 0 (clear-felled area) to 219.4 Mg ha −1 (30-year-old stand) for the AGB and from 0 to 36.1 Mg ha −1 for the BGB. Soil carbon in the top 1 m of soil had a range of 385.2 Mg ha −1 to 497.5 Mg ha −1 from a clear-felled area to a 30-year-old site (Figure 3). Considering these values, the total carbon stock in the biomass for the productive forest was equal to 1,893,571.1 Mg C, while the total for the protective forest was at least equal to 1,623,427.7 Mg C. Regarding the soil carbon pool, this was equal to 12,683,428.4 Mg C and 5,801,745.5 Mg C for the productive forest and the protective forest, respectively. Table 2. Total carbon stock (C) in AGB, BGB, and top 1 m of soil for the productive, restrictive productive, and protective forest zones in the MMFR. For the protective forest the biomass and carbon per ha was considered at least equal to a 30-year-old stand, as these forest areas are in many cases older than that; hence, the values reported for the protective forest are an underestimation of the real carbon stock.

Area (ha)
Total   Biomass data used to produce the map were derived from Lucas et al. [38] and Goessens et al. [30] and soil data were estimated from Adame et al. [35]. The white areas on the map are mostly protective forest zones for which biomass values were not directly available.    Biomass data used to produce the map were derived from Lucas et al. [38] and Goessen et al. [30] and soil data were estimated from Adame et al. [35]. The white areas on the map are mostly protective forest zones for which biomass values were not directly available.

Carbon Emissions Factors from Charcoal and Pole Production
For the production of charcoal, transportation of people/logs in the boats had the minimum average carbon emission factors (1.14 g C km −1 ), while greenwood conversion into charcoal showed the maximum (702 g C kg −1 ) ( Table 3). In the case of mangrove pole production, transport by boats had the lowest carbon emission factor (1.16 g C km −1 whereas transport by lorries distributing poles to local markets had the highest (167.45 g C km −1 ) ( Table 3). Table 3. Average carbon emission factors per each activity in the two production chains (char Biomass data used to produce the map were derived from Lucas et al. [38] and Goessens et al. [30] and soil data were estimated from Adame et al. [35]. The white areas on the map are mostly protective forest zones for which biomass values were not directly available. Considering the total areas occupied by different forest zones (productive, restrictive productive, and protective), the total C stock for the MMFR (AGB, BGB, and soil combined) would be equal to 36,212,677.1 Mg C.

Carbon Emissions Factors from Charcoal and Pole Production
For the production of charcoal, transportation of people/logs in the boats had the minimum average carbon emission factors (1.14 g C km −1 ), while greenwood conversion into Remote Sens. 2022, 14, 2920 9 of 17 charcoal showed the maximum (702 g C kg −1 ) ( Table 3). In the case of mangrove pole production, transport by boats had the lowest carbon emission factor (1.16 g C km −1 ) whereas transport by lorries distributing poles to local markets had the highest (167.45 g C km −1 ) ( Table 3). Table 3. Average carbon emission factors per each activity in the two production chains (charcoal and poles). These values were calculated considering the brands and models of the equipment used, number of workers, and mean distances travelled.

Carbon Emission from the Areas Assigned to One Contractor
The area of 2.2 ha allocated for charcoal production in a year provides an average yield of 396 Mg (i.e., 180 Mg ha −1 ). According to the management plan, 40.8 Mg of greenwood (with additional 12.2 Mg as firing material) is required per kiln to produce 11 Mg of charcoal in 40-45 days. On this basis, 304.85 Mg out of 396 Mg greenwood in 2.2 ha is used for charcoal production and 91. 15 Mg as firewood which ultimately delivers 82.08 Mg of charcoal. After considering the amount of wood used in charcoal production and its yield, along with the number of workers involved and the kilometres travelled for transport, the carbon emission estimates were found to be ranging from 0.001 Mg C for boat transport to 278 Mg C for greenwood burning inside the kilns for a total of 336.17 Mg C ( Figure 4). Likewise, the pole production from a 22.3 ha area of intermediate thinning, together with average number of workers involved in the different activities, amounted to carbon emission ranging from 0.04 Mg C for boat/lorry transports to 9.31 Mg C for car transport, for a total of 11.82 Mg C (Figure 4).

Carbon Emission from the Total Areas of Pole and Charcoal Production
Dividing the emissions for the area assigned to one contractor by the area in hectares, the carbon emission factor per ha would be equal to 152.8 Mg C/ha for charcoal production and 0.53 Mg C ha −1 for pole production.
Considering the total areas for thinning I and II and clear-felling reported by the latest management plan of 2010-2019 [26] equal to 8159 ha, 7678 ha, and 11,593 ha, respectively (named as allocated areas in Table 4), the total emissions are of 4.32 Gg C for pole production from thinning I, 8.13 Gg C for pole production from thinning I and 1.77 Tg C for charcoal production.

Carbon Emission from the Total Areas of Pole and Charcoal Production
Dividing the emissions for the area assigned to one contractor by the area in hectares, the carbon emission factor per ha would be equal to 152.8 Mg C/ha for charcoal production and 0.53 Mg C ha −1 for pole production.
Considering the total areas for thinning I and II and clear-felling reported by the latest management plan of 2010-2019 [26] equal to 8159 ha, 7678 ha, and 11,593 ha, respectively (named as allocated areas in Table 4), the total emissions are of 4.32 Gg C for pole production from thinning I, 8.13 Gg C for pole production from thinning I and 1.77 Tg C for charcoal production.

Carbon Loss from Soil after Clear-Felling
Based on the data of Adame et al. [35], the carbon loss from the soil due to clear-felling activities was also calculated. On average, the carbon stock in the soil of a clear-felled site is 385.2 Mg C ha −1 , while a 30-year-old site has a carbon stock of 497.5 Mg C ha −1 . Consequently, there is a loss of 112.3 Mg C ha −1 each time an area is cleared for charcoal production. By considering the total area assigned for clear-felling in the 2010-2019 management plan (11,593 ha), the total carbon loss from the soil would be equal to 1,301,893.9 Mg C.

Carbon Stock Versus Carbon Emission
The biomass carbon stock was compared with the carbon emissions considering the emissions due to the activities practiced during the first and second thinning (pole production) and clear-felling (charcoal production), as well as the carbon loss from the soil after clear-felling (Table 4). For a licensed area assigned to one contractor, i.e., 22.3 ha for thinnings and 2.2 ha for clear-felling, the emissions are much lower for the pole production compared with the charcoal production. On the other hand, in the case of a clear-felling area, the emissions are higher compared with the biomass carbon stock.
The carbon budget of the productive forest was computed taking into account the areas allocated for poles and charcoal reported in the latest management plan and their derived emission values together with the vegetation carbon stock calculated from the remote sensing data. The AGB carbon stock of the productive forest alone, equal to 1,402,645.3 Mg C (Table 1) was compared with the emissions of 1,771,463.09 Mg C (Table 4). Hence, the carbon budget of the silvicultural activities were negative (−381,228.1 Mg C) ( Figure 5). Considering the total area of the MMFR, i.e., including protective and restrictive productive zones, the total carbon stock was equal to 2,698,331.7 Mg C (Table 1), which compared with the emission values reported above, resulted in a positive carbon budget (+341,113.4 Mg C) ( Figure 5).
Based on the data of Adame et al. [35], the carbon loss from the soil due to clearfelling activities was also calculated. On average, the carbon stock in the soil of a clearfelled site is 385.2 Mg C ha −1 , while a 30-year-old site has a carbon stock of 497.5 Mg C ha −1 . Consequently, there is a loss of 112.3 Mg C ha −1 each time an area is cleared for charcoal production. By considering the total area assigned for clear-felling in the 2010-2019 management plan (11,593 ha), the total carbon loss from the soil would be equal to 1,301,893.9 Mg C.

Carbon Stock Versus Carbon Emission
The biomass carbon stock was compared with the carbon emissions considering the emissions due to the activities practiced during the first and second thinning (pole production) and clear-felling (charcoal production), as well as the carbon loss from the soil after clear-felling (Table 4). For a licensed area assigned to one contractor, i.e., 22.3 ha for thinnings and 2.2 ha for clear-felling, the emissions are much lower for the pole production compared with the charcoal production. On the other hand, in the case of a clearfelling area, the emissions are higher compared with the biomass carbon stock.
The carbon budget of the productive forest was computed taking into account the areas allocated for poles and charcoal reported in the latest management plan and their derived emission values together with the vegetation carbon stock calculated from the remote sensing data. The AGB carbon stock of the productive forest alone, equal to 1,402,645.3 Mg C (Table 1) was compared with the emissions of 1,771,463.09 Mg C (Table  4). Hence, the carbon budget of the silvicultural activities were negative (−381,228.1 Mg C) ( Figure 5). Considering the total area of the MMFR, i.e., including protective and restrictive productive zones, the total carbon stock was equal to 2,698,331.7 Mg C (Table 1), which compared with the emission values reported above, resulted in a positive carbon budget (+341,113.4 Mg C) ( Figure 5).
Moreover, if we also consider the soil carbon, the total carbon stock of the productive zones of the MMFR would be equal to 14,576,999.5 Mg C. On the same line, the total emissions, including the carbon loss from the soil, would be equal to 3,073,357.0 Mg C. This was not included in the budget calculation as the carbon in the soil accumulates over time and cannot be compared with the emissions of a thirty-year cycle.  Moreover, if we also consider the soil carbon, the total carbon stock of the productive zones of the MMFR would be equal to 14,576,999.5 Mg C. On the same line, the total emissions, including the carbon loss from the soil, would be equal to 3,073,357.0 Mg C. This was not included in the budget calculation as the carbon in the soil accumulates over time and cannot be compared with the emissions of a thirty-year cycle.

Vegetation Biomass and Soil Carbon Stock
Given the key role of mangroves as carbon stock, we measured the carbon sequestration potential of MMFR, and subsequently compared this with the C emissions originating from pole and charcoal production. The total carbon in the soil was also evaluated to provide an overview of the total carbon stock of the MMFR.
Regarding the vegetation biomass, remote sensing data are a valuable source for biomass estimations in mangroves and can provide a very detailed view of an area over time [37,38]. Furthermore, we integrated data from field measurements [30,35] with the age maps obtained through remote sensing [37], producing for the first time an accurate representation of the carbon storage in the entire MMFR that includes not only vegetation but also the belowground component (i.e., roots and soil).
In this study the carbon stock of the MMFR was firstly calculated based on AGB data derived from remote sensing analysis. The main advantage of the use of remote sensing data in this study is that each forest coupe was considered with its own age and biomass values whereas previous biomass estimations were based on field data collection alone (e.g., Goessens et al. [30]) and could not provide an estimate for the entire forest as only specific forest ages were selected for the study (namely 15-, 20-, 30-year-old and protective forest). Thus, our results should provide a more reliable overview of the carbon stored in the MMFR. Furthermore, to have a complete understanding of the carbon stored in the ecosystem, other components such as belowground and soil carbon, which are often overlooked when only remote sensing data are used, were considered. Previous studies [8,35,[53][54][55][56][57] all reported the importance of belowground carbon which is considered to store from 50 to 90% [50] of the total mangrove ecosystem carbon stock. Until now the only study in the MMFR that measured soil carbon in the soil was Adame et al. [35] and their estimates ranged (approximately) from 380 Mg C ha −1 for a 1-year-old forest to 480 Mg C ha −1 for a 30-year-old forest and reached 540 Mg C ha −1 in the protective forest (for the first meter of soil). By extrapolating soil carbon values for all forest ages and linking them to the areal extent of each forest zone, we estimated that the total carbon storage in the soil can reach 31,637,347.8 Mg C. Hence the carbon storage in the soil seems to be much greater than the carbon stock in the vegetation biomass, with this estimated to be 2,698,331.7 Mg C.
However, the above estimates only consider the carbon in the soil in the first 1 m while it is known that mangrove soils can be much deeper than 1 m [50]. Thus, further research should be focused on deeper soil layers to obtain a more accurate representation of the carbon stock in the mangrove forest.

Charcoal and Pole Production Chains: Carbon Emissions
The different emissions in the two production chains were expected considering the difference in the number of activities and workers involved in each type of production, with total carbon emission values of 11.82 Mg C on average for pole production and a total average of 336.17 Mg C for charcoal production per licensed area. In other words, pole production in MMFR was found to release on average 84% less carbon than charcoal production.
Nevertheless, during our calculations of emissions in the various stages of the production chains, we noticed a lack of easily available specific information (e.g., regarding the mode of transport used, chainsaw brand, and specific number of users). For wood burning, other research studies such as Lacaux et al. [41], Kridiborwnorn et al. [43] and Sparrevik et al. [58] used a complex experimental design by considering more detailed values (e.g., wood carbon (%), wood calorific value (kJ/g), brands carbon (%)) with more comprehensive and complete results. At the same time, information about exhaust emission released by transports was limited or difficult to follow because of the need for specific data that were unavailable for our study. Moreover, the other factors such as the local guards' role, which involves use of the boats to monitor the productive area, and the activities linked to mangroves replanting with seedlings were not considered in the present study due to lack of information.

Vegetation Carbon Stock Versus Carbon Emissions
A comparison between the values for carbon sequestration and carbon emissions at MMFR was conducted taking into account the different areas of exploitation. Considering the biomass carbon stock in the productive area compared with the emissions from the pole and charcoal production for a thirty-year cycle, the emissions were estimated to be greater than the carbon stock. This is because the activities included in the production chain also use fuels other than mangrove wood, making the budget negative. If the forest was used entirely for charcoal production, the management would be unsustainable from a carbon point of view. Furthermore, if we also include the emissions from the soil after clear-felling calculated by Adame et al. [35] who estimated a drop of 43% in soil carbon after clearing, the emissions would be even higher.
In order to achieve carbon neutrality for the productive forest, it is necessary to reduce the area allocated for clear-felling by at least 20.8%, implying a reduction in the productive forest from 74.8% to 45.6% of the entire MMFR. Fortunately, even in the current situation, the MMFR also includes areas that are not exploited and are left unmanaged (i.e., protective forest zones) that can compensate for the loss of carbon from the productive zones. When the carbon stock of these areas was considered into the carbon budget, it was positive, providing the management a more sustainable outlook. However, the current management does not take into account possible disruptions caused by climate change (e.g., rising sea levels and increased storm intensity), that can affect the yield and functionality of the mangrove forest.
Given the importance of carbon sequestration for climate change mitigation, it would be better for the MMFR to increase the protective zone areas and REDD+ or similar projects can be considered as an alternative source of income, as already suggested by other studies [35,59].
On the contrary, in recent years an increase in exploitation of the restrictive productive zones has been observed (Harry Yong, personal communication, June 2019), probably due to the decrease in yield caused by mangrove degradation [60] and the parallel high demand for charcoal, leading to more clear-felling activities and subsequently higher emissions. This is concerning as the MMFR seems to slowly become more unsustainable in its silvicultural practices.

Conclusions
The present study provided valuable information regarding the emissions released by pole and charcoal production in the MMFR. By using the most recent biomass and carbon data we provided an accurate overview of the detailed carbon stock in MMFR and for the first time this was integrated with emission estimates caused by charcoal and pole production. Therefore, this study can be used as a pilot project to evaluate similar silvicultural mangrove management schemes worldwide.
At the carbon sequestration side, many unknowns remain, particularly regarding the depth of soils under mangroves and the actual carbon stock belowground together with how the thinning process, the clear cutting, and the replanting stages influence the carbon cycle dynamics. Furthermore, data on other ecosystem components such as leaf litter, dead wood, and fauna is very limited.
In general, considering the values obtained in this study, we can conclude that the total area of the MMFR can be considered carbon positive, despite the emissions caused by the silvicultural activities. However, we suggest a reduction in the productive zones of at least 20% to achieve carbon neutrality in the productive forest coupes.
This exploratory study helps to strengthen the assumption that the protection of mangrove forest must be a priority as they are among some of the carbon-richest ecosystems in the world. We also show that a sustainable exploitation of mangroves for pole and charcoal production is possible from a carbon point of view if the exploited area is managed properly. In this way mangroves can continue to be an excellent carbon stock and the livelihood of local people that depend on this ecosystem can be sustained.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/rs14122920/s1, Table S1: Sources used to calculate carbon emission factors from charcoal and pole production; Table S2: Parameters taken into account to estimate carbon emission per contractor for each type of production; Table S3: Emission factors for the Husqvarna chainsaws models; Table S4: Emission factors for the logging boats; Table S5: Distance (km) from the respective Charcoal Factory to the coupes available for Felling; Table S6: Emission factors, grams of pollutant per kilogram charcoal produced in kiln; Table S7: Emission factors for heavy truck transport (g.km −1 ); Table S8: Distance (km) from the charcoal factories to the Companies in charge of the charcoal exportation.; Table S9: Distance (km) from the Companies to the Port terminals in Malaysia for charcoal exportation; Table S10: Distance (km) from the charcoal factories to the registered workers' residence; Table S11: Estimates of exhaust emission factors for cargo ships (g/ tonne-km); Table S12: Distance from Malaysia port terminals (Port Klang and Penang) to Japanese terminals for charcoal exportation; Table S13: Car Standard Emissions (g/km); Table S14: Moped Standard Emissions (gC.km −1 ); Table S15: Distance in km from the charcoal factory to the places where pole are sold; Table S16: Emission factors for light and heavy truck transport (gC.km −1 ) used for pole transport. The references in SM can refer to [26,41,42,44,[61][62][63][64][65][66][67][68][69][70][71].