Skip to main content

Improving soil carbon estimates of Philippine mangroves using localized organic matter to organic carbon equations

Abstract

Background

Southeast Asian (SEA) mangroves are globally recognized as blue carbon hotspots. Methodologies that measure mangrove soil carbon stock (SCS) are either accurate but costly (i.e., elemental analyzers), or economical but less accurate (i.e., loss-on-ignition [LOI]). Most SEA countries estimate SCS by measuring soil organic matter (OM) through the LOI method then converting it into organic carbon (OC) using a conventional conversion equation (%Corg = 0.415 * % LOI + 2.89, R2 = 0.59, n = 78) developed from Palau mangroves. The local site conditions in Palau does not reflect the wide range of environmental settings and disturbances in the Philippines. Consequently, the conventional conversion equation possibly compounds the inaccuracies of converting OM to OC causing over- or under-estimated SCS. Here, we generated a localized OM-OC conversion equation and tested its accuracy in computing SCS against the conventional equation. The localized equation was generated by plotting % OC (from elemental analyzer) against the % OM (from LOI). The study was conducted in different mangrove stands (natural, restored, and mangrove-recolonized fishponds) in Oriental Mindoro and Sorsogon, Philippines from the West and North Philippine Sea biogeographic regions, respectively. The OM:OC ratios were also statistically tested based on (a) stand types, (b) among natural stands, and (c) across different ages of the restored and recolonized stands. Increasing the accuracy of OM-OC conversion equations will improve SCS estimates that will yield reasonable C emission reduction targets for the country’s commitments on Nationally Determined Contributions (NDC) under the Paris Agreement.

Results

The localized conversion equation is %OC = 0.36 * % LOI + 2.40 (R2 = 0.67; n = 458). The SOM:OC ratios showed significant differences based on stand types (x2 = 19.24; P = 6.63 × 10–05), among natural stands (F = 23.22; p = 1.17 × 10–08), and among ages of restored (F = 5.14; P = 0.03) and recolonized stands (F = 3.4; P = 0.02). SCS estimates using the localized (5%) and stand-specific equations (7%) were similar with the values derived from an elemental analyzer. In contrast, the conventional equation overestimates SCS by 20%.

Conclusions

The calculated SCS improves as the conversion equation becomes more reflective of localized site conditions. Both localized and stand-specific conversion equations yielded more accurate SCS compared to the conventional equation. While our study explored only two out of the six marine biogeographic regions in the Philippines, we proved that having a localized conversion equation leads to improved SCS measurements. Using our proposed equations will make more realistic SCS targets (and therefore GHG reductions) in designing mangrove restoration programs to achieve the country’s NDC commitments.

Background

Blue carbon ecosystems (BCEs) are vegetated coastal ecosystems (e.g., mangroves, seagrasses, and tidal marshes) that sequester disproportionately large amounts of organic carbon (OC) per area [1]. Their anoxic and waterlogged soils capture and sequester OC by slowing down decay rates [2]. Mangrove forests are one of the important BCEs that rose into prominence for their unique role in the global climate crisis by sequestering atmospheric carbon while providing many ecosystem services [3]. The role of mangroves in climate change adaptation and mitigation (CCAM) through conservation and restoration programs has been emphasized in many countries’ Nationally Determined Contributions (NDCs) under the Paris Agreement [4]. However, when mangroves are disturbed, their stored carbon are released as greenhouse gases (GHG) back into the atmosphere [5].

Mangroves contain the largest carbon stock (CS) per unit area when compared to any other ecosystems in the tropical ocean (6.17 Pg Corg; 17% of the total tropical marine CS, [6]). However, wide variability in CS estimates exists due to local macroscale ecosystem differences (e.g., forest age, geomorphology, among others, [7]). Mangrove soil organic carbon (SOC) locks in 49–98% of the total ecosystem carbon stock (TECS), but is the least studied carbon pool [8, 9] particularly in tropical developing countries. Elemental analyzers that yield direct OC measurements are costly, making them impractical to many tropical developing countries. OC concentration is often estimated semi-quantitatively by converting organic matter (OM) derived from the more economical and more accessible loss-on-ignition (LOI) method using a conversion equation. The usual conversion equation applied in most sediment carbon stock (SCS) estimates in the Philippines is derived from Palau mangroves, %Corg = 0.415 * % LOI + 2.89 ([9]; referred from hereon as the conventional equation). Inherent inaccuracies are expected when determining OC through the LOI method [10]. Site-specific differences, especially on local climate and geomorphology, compound these inherent inaccuracies [6] leading to less accurate SCS measurements [1, 6]. Ultimately, the accuracy of the conversion equation relies heavily on the degree of similarity between the study site and the reference site where the conventional equation was derived [11].

The Southeast Asian (SEA) region is a “blue carbon hotspot” hosting the largest and most diverse mangroves in the world. Among the SEA countries, the Philippines has the third largest mangrove CS (102–576 MgC ha−1 [12]) and the fourth largest mangrove extent (284,798 ha [13]). Philippine mangroves exist in a variety of stand conditions (natural, colonized, and restored) based on vegetation structure and the developmental stage with different climatic and geomorphological conditions [14]. These mangroves are also frequently affected by natural and anthropogenic disturbances such as typhoons and conversion to aquaculture ponds [15, 16]. Applying a general OM-OC equation for the Philippines may be prone to misuse as it may discard site-specific mangrove conditions (based on stand types, occurrences/types of disturbances, etc.). These country-specific conditions are less prominent in Palau mangroves hence may not necessarily be accounted for by the conventional OM-OC conversion equation [9, 17].

In this study, we generated localized OM-OC site-specific conversion equations representative of the different mangrove stand types (e.g., natural, colonized, restored, and no vegetation) in the Philippines. The soil OM:OC ratios and the slopes of the conversion equations for the (a) natural stands and the (b) restored and recolonized stands based on chronosequence of ages were compared to infer the influence of mangrove age (cf. [18]). The accuracy of the generated equations was then compared to the conventional equation.

Methods

Site description

The study was conducted in two areas in the Philippines: Prieto Diaz in the province of Sorsogon (13.0179300°, 124.1867900°), and the municipalities of Bongabong (12.6991487°, 121.5307829°), Mansalay (12.4338929°, 121.4114457°), and Roxas (12.6046100°, 121.5443800°) in the province of Oriental Mindoro (Fig. 1). The municipality of Prieto Diaz is located on the easternmost point of the Bicol Peninsula on the southern portion of Luzon Island and belongs to the North Philippine Sea biogeographic region. The sites in Oriental Mindoro are on the southeastern portion of Mindoro Island and belong to the West Philippine Sea biogeographic region. Prieto Diaz has a tropical rainforest climate while Oriental Mindoro has a tropical monsoon climate based on the Koppen Classification (Table 1 [19]). Both sites are categorized as open coast (based on typology, cf. [20]). Natural stands are mangroves that are not planted nor recolonized. Restored mangroves are “planted” stands while recolonized mangroves are mangroves that undergo succession post-fishpond abandonment. The estimated ages of the stands were based on local accounts/records complemented with local maps. We also added a site without mangroves (referred to as “no vegetation”) to depict clear-cut vegetation and represent disturbed mangroves.

Fig. 1
figure 1

Location of sampling sites and mangroves within the study areas (A Oriental Mindoro: Bongabong Recolonized, BC; Roxas Restored, RR; Roxas Natural, RN; Mansalay Natural, MN; Mansalay Recolonized, MC; and B Prieto Diaz: Natural, PN; Restored, PR; Young recolonized, PCY; Mature recolonized, PCM). Inset map shows the location of study areas within the Philippine marine biogeographic regions

Table 1 Biophysical profile of the study sites

Field sampling

Sampling was conducted in nine sites across different mangrove stand types and stand ages (Table 1). The “no vegetation” samples were collected adjacent to the colonized mangroves in Mansalay. Within each site, core samples were collected across a zonation gradient, if applicable. Sediment cores of up to one meter depth were collected using a 6.5 cm diameter open-faced auger. The core samples were cut every two cm for the upper 0–20 cm surface, and every four cm for the 20–100 cm depths. Each subsample was individually placed in airtight polyethylene bags. All soil samples from Sorsogon were collected in February 2022 while samples from Bongabong, Mansalay, and Roxas were collected in July 2022.

Laboratory analysis

In the laboratory, each sub-sample was air-dried then weighed. The samples were then placed in a drying oven at 60 °C until constant mass was attained. Bulk density was calculated using the formula below:

$$Bulk\, Density\, (g/{cm}^{3})=\frac{weight\, of\, subsample\, (g)}{volume\, of\, subsample\, ({cm}^{3})}$$

The OM content was determined through the LOI method (cf. [9]). Dried samples were pulverized with an agate mortar and pestle. Stones, twigs, roots, and other large materials were removed during pulverization (cf [11]). The ground samples were then oven-dried at 60 °C for 24 h. Five grams per subsample were weighed, then placed in a muffle furnace at 550 °C. We standardized combustion for four hours. The change in mass after combustion represented the OM concentration.

$$\%OM=\frac{initial\, mass\, \left(g\right)-final\, mass(g)}{initial\, mass (g)}\times 100$$

A separate set of samples was sent to Jakarta, Indonesia for the analysis of OC content using an elemental analyzer (LECO CHN 628 Series: CHN Analyzer, LECO Corp., St. Joseph, MI). In total, 458 soil samples from different sites, stand types, and zonation were analyzed separately for OM and OC concentrations.

Data analyses

Using simple linear regression, an allometric equation was derived to determine the relationships between the OM and OC concentrations across sites and stand types (referred from hereon as general equation). Allometric equations were also derived for each stand type and age (referred from hereon as stand-specific and age-specific equation, respectively). The SCS of each sampling site was calculated using the generated general OM-OC equation and then compared with the OC derived from the elemental analyzer and the conventional equation derived from Palau mangroves [9].

The OM derived from the LOI method and the OC from the elemental analyzer were used to create OM:OC ratios. The OM:OC ratios were compared across stand types, among natural stands, and among ages for the recolonized and restored stands. The OM:OC ratios were first filtered for outliers (categorized as values greater than or equal to 1.0, following [10]). Shapiro–Wilk normality test and Levene’s homogeneity test were then applied to the filtered data then analyzed using one-way ANOVA (followed by Tukey’s Highly Significant Difference for post-hoc comparisons). Non-normal data was analyzed using the Kruskal Wallis Test then post-hoc Dunn test. Data uncertainties were expressed in a 95% confidence interval (CI). All data analyses and visualizations were performed using R version 4.0.3 and RStudio version 1.4.1103.

Results

Organic carbon strongly correlated with organic matter content

The OM% (derived from the LOI) was significantly correlated with the OC% (derived from the elemental analyzer), %Corg  =  0.36 * % LOI  +  2.40 (R2 = 0.67, x2 = 3; P = 0.28; Fig. 2). There were no significant differences in regression slopes across stand types. However, the strength of relationship varied with stand types (Restored, R2 = 0.59; Recolonized, R2 = 0.69; Natural, R2 = 0.60; Fig. 3; Table 2a). The natural stands had the steepest slope (0.39), closely followed by the recolonized stands (0.38). Both the natural and recolonized stands have 23% steeper slopes than that of the restored stands (0.30; Fig. 3).

Fig. 2
figure 2

Scatter plot of organic matter (OM) and organic carbon (OC) contents used to derive the general OM-OC conversion equation for Philippine mangroves. The OM was calculated using the loss-on-ignition (LOI) method, while OC was calculated using the conversion equation, %Corg = 0.415 * % LOI + 2.89 [9]

Fig. 3
figure 3

Stand-specific OM-OC conversion equations derived for natural stands (A), restored stands (B), recolonized stands (C), and no vegetation (D)

Table 2 Summary of slopes, intercepts and R2 values of the OM-OC conversion equations (A) and the results of the analysis of variance (ANOVA) tests on the OM:OC ratios (mean ± CI) according to stand-specific and site-specific parameters (B)

Variability of OM:OC ratios within and between stand types

The OM:OC ratios of the restored stands (2.25 ± 0.05) were 11% and 26% significantly lower than the natural (2.53 ± 0.09) and recolonized (3.03 ± 0.14) stands, respectively (Table 2b). Variabilities within groups were also observed. The natural mangroves in Mansalay (1.93 ± 0.06) had 29–31% lower ratio than Prieto Diaz (2.81 ± 0.22) and Roxas (2.74 ± 0.09). In the restored stands, the 10-yr-old plantations (2.12 ± 0.06) had 11% lower ratio than the 30-yr-old plantations (2.38 ± 0.08). In the recolonized stands, the ratios varied with stand ages although pairwise comparisons didn’t yield significant differences between ages (8-yrs = 2.25 ± 0.07; 10-yrs = 2.77 ± 0.19, 17-yrs = 2.23 ± 0.08, 35-yrs = 3.85 ± 0.30).

Locally-derived equation has higher accuracy than the conventional equation

The mean SCS estimated using %OC from the elemental analyzer was 585.61 ± 72.16 Mg/ha (Fig. 4A). Our general equation was 5% higher across stands (616.04 ± 64.09 Mg/ha; Fig. 4A) and 7% higher with stand-specific calculations (625.56 ± 61.79 Mg/ha; Fig. 4B) than the values derived from the elemental analyzer (Fig. 4C). In contrast, the conventional equation overestimates SCS by 20% (718.69 ± 71.85 Mg/ha; Fig. 4C). Our derived general and stand-specific equations yielded higher SCS for the natural stands (%Corg = 0.39 * % LOI + 0.94 [R2 = 0.60, x2 = 2; P = 0.37; Fig. 3A]), but lower SCS for both the restored (%Corg = 0.38 * % LOI + 2.3 [R2 = 0.69, x2 = 1; P = 0.32; Fig. 3B; Fig. 5]) and recolonized stands (%Corg = 0.30*% LOI + 3.6 [R2 = 0.59, x2 = 3; P = 0.32; Fig. 3C; Fig. 5]).

Fig. 4
figure 4

(Mean ±  CI) SCS across stand types using the stand-specific equation (A), across ages in the restored and recolonized stands using the age-specific equations (B), and pooled comparison with an elemental analyzer relative to the general, stand-specific and the conventional equations (C; [9])

Fig. 5
figure 5

Stand-specific scatter plots of the OM to OC ratios per stand type (A), among natural stands (B), across ages of the restored stands (C), and across ages of the recolonized stands (D)

Discussion

Our study generated general and stand-specific OM-OC conversion equations for Philippine mangroves. The equations we developed for Philippine mangroves, specifically for open coast typology [20], are within the range of the slopes (0.35 vs 0.23–0.50) and correlation strengths (R2 = 0.67 vs. 0.25–0.79) reported in SE Asia and Micronesia (Table 3). The regression slopes in the stand-specific and age-specific equations, although not statistically significant between and within groups, proved that localized equations will be more appropriate to calculate SCS compared to the conventional equation [9]. To our knowledge, these equations were the first account that established OM-OC relationship that was used in the measurements of SCS in Philippine mangroves.

Table 3 Comparison of OC-OM relationships from SCS studies in Southeast Asia and Micronesia

The need for a localized OM-OC equation

Blue carbon studies in Philippine mangroves have increasingly become popular since 2009 [21]. With the need to utilize blue carbon as a nature-based solution (NbS) and as an integral part of CCAM strategy, the topic is anticipated to become even more prominent as a priority research agenda in the country [22]. Empirical studies are therefore needed to provide science-based policy support in the conservation and restoration of mangroves.

In the Philippines, empirical blue carbon studies are limited and concentrated in few places (Table 4). Out of the published studies, most are limited to the assessment of aboveground carbon stocks (AGCS, n = 62) as compared to the SCS (n = 38; [14]). This could be due to the labor-intensive nature of sampling, lack of adequate funding, and lack of instrument needed for the direct measurement of OC (i.e., through an elemental analyzer; [10]). In determining soil OC, the most used is the LOI method [23, 24] because of cheaper cost despite being considered as semi-quantitative [25]. Some studies used the Walkley–Black method [26,27,28]; see also Table 4. Although elemental analyzer is used in some studies (n = 3; [29,30,31]), it is less used primarily because of the expensive analytical costs.

Table 4 SCS estimation methods of mangrove carbon stock studies in the Philippines

The calculated OM from LOI method is typically converted into OC for SE Asian mangroves as %Corg = 0.415 * % LOI + 2.89 (R2 = 0.59 [9]) and then used in the calculation of SCS. We acknowledge that some of our samples were below the < 18% OM threshold for the OM to OC conversion [9]. But, this OM to OC conversion equation has been widely used, not just in the Philippines [14], but in most SE Asian countries [32], despite it being derived from limited sampling points and having low correlation coefficient. It has already been criticized (see for example [33]) because of the overestimation of SCS. The overestimation is likely due to the small sample sizes (n = 78) and the differences in biophysical conditions between the site that the equation was developed from (Palau and Yap Islands [9, 17]) and the site it was applied to (for example this study). When compared to the mangroves in Palau and Yap Islands, our mangroves are visited more frequently by stronger catastrophic typhoons. Typhoons have profound effects on the wash out (i.e., decreased SCS) and possible post-typhoon recovery (i.e., increased SCS) especially in the upper sediment strata [34].

Our study supports the assertion that site-specific OM-OC equations (with larger sample sizes) are needed for more accurate calculations of OC [10]. The localized equations will be more sensitive to account the different mangrove vegetation conditions and local geomorphologies. And in the case of the Philippines, a localized equation will be needed to account for the effects of disturbances (e.g., typhoons, conversion of mangroves to aquaculture ponds) and post-disturbance recovery of carbon. Having a locally derived OM-OC conversion equation can enhance the SCS estimation by improving the accuracy of the LOI method and reducing reliance on costlier and less accessible elemental analyzers.

Limitations, opportunities, and further improvements of localized OM-OC conversion equations

Estimating SCS using the LOI method requires intensive labor in the laboratory to determine the OM. Any inaccuracy in OM determination and the subsequent conversion to OC will also lead to inaccuracy in SCS estimates. The conventional equation [9] overestimate SCS by 20%. In comparison, our proposed equation for the open coast Philippine mangroves (%OC = 0.36 * % LOI + 2.40 (R2 = 0.67, n = 458) is within 5%–7% when compared with the values derived from the elemental analyzer. Despite having a relatively stronger relationship, our proposed equation will still produce uncertainties, as with other conversion equations due to the semi-quantitative nature of the LOI method [10].

The sources of carbon are generally categorized as autochthonous or allochthonous [35], hence the stability of the accumulated carbon largely depend on the state of ecosystem health and geomorphological conditions of mangroves. In general, the factors that contribute to SCS are biomass, latitude, precipitation, tidal range, soil pH, and soil depth [36]. In the Philippines, out of the available empirical datasets, the factors that contribute to SCS are related to the “ecosystem health” and geomorphological conditions, e.g., temperature, redox, biomass, mangrove area, stand age, and latitude [10]. Our sites are composed of wide ecosystem status and environmental settings, from natural to recolonized, and restored stands of different ages. These stands have also likely undergone periods of disturbance and regeneration (e.g., primarily typhoons). Although our proposed equation is already relatively better than the conventional equation, it can still be improved by applying stand-specific equations (Fig. 3) consistent with the proposition of Breithaupt et al. [10]. Otherwise, SCS will still either be over- or under-estimated because it does not integrate yet the chronosequence of changes in carbon linked from the bare vegetation (synonymous to a situation when mangroves are degraded; Fig. 3d) to the vegetation growth and development in the restored (Fig. 3b) and recolonized stands (Fig. 3c; see also the differences in OM-OC ratios across stands; Table 2b). In fact, even for natural mangroves, there were some variabilities in regression slopes (although not statistically significant; Table 2). Furthermore, although our equation encompasses a large variety of site conditions, it only represents two (West Philippine Sea and North Philippine Sea) out of the six marine biogeographic regions in the Philippines [37]. Based on our results, it is likely that a biogeographic region-specific conversion equations will be more realistic.

Implications for conservation and restoration

The Philippines experienced a massive reduction of about 46% in mangrove cover from about 450,000 ha in 1918 to about 264,818 ha in 2020 [38] caused by various anthropogenic and climate-related factors [39]. Mangrove losses contribute to CO2 emission (globally estimated at 0.09–0.45 Pg CO2/year, [40]). In accordance with the Paris Agreement, the country committed to reduce and avoid its projected GHG emissions to 75% by 2030 in its first NDC submitted to the United Nations Framework Convention on Climate Change (UNFCCC) in 2021 [41]. As successive NDC’s call for adaptive climate change mitigation program every five years, preparing for future NDCs is a critical opportunity to revisit the current NDC and consider how to meet the country’s commitment in the Paris Agreement.

The Philippines is one of the three countries in Southeast Asia, alongside Indonesia and Malaysia, that incorporates conservation and restoration of blue carbon ecosystems in their NDCs [42]. In Indonesia, optimal mangrove restoration sites have been prioritized nationwide with the aim to contribute to achieving the country’s rehabilitation targets and to reduce national GHG emissions [43]. Increased accuracy of CS estimates of Philippine mangroves through the use of our proposed localized conversion equation can help provide insights on the extent of restoration required to achieve the country’s NDC by 2030. These sites include mangrove-recolonized abandoned, undeveloped, or underutilized (AUU) fishponds which have exhibited similar levels of CS with natural mangroves, decades after restoration [44]. The rehabilitation of recolonized AUU fishponds has been noted for its potential in contribution to mangrove conservation [45,46,47]. Our proposed stand-specific equations can provide further empirical evidence on the carbon storage potential of different mangrove stands, particularly of recolonized ponds, further incentivizing their rehabilitation.

The TECS estimation has been used to determine the success of restoration in forest ecosystems [44, 48,49,50,51]. However, only a limited number of studies in the Philippines has assessed CS as a restoration indicator [44, 52]. Our proposed equation could make SCS estimation more accurate. Additional studies could provide more evidence on the effectiveness of various rehabilitation and restoration strategies and determine the influence of disturbances on SCS. Strategic mangrove restoration combined with focused conservation using more accurate OM to OC conversion equation will be key in achieving the national commitment on Paris Agreement.

Conclusions

The proposed localized OM-OC conversion equations were more accurate in estimating SCS of Philippine mangroves than the conventional conversion equation. The proposed equations are also comparable to published equations with other Southeast Asian mangroves in terms of slopes and relationship strengths. The existence of localized OM-OC conversion equations can promote further research on mangrove SCS in the Philippines by increasing the accuracy of the LOI method and reducing reliance on costlier and less accessible elemental analyzers. Increased accuracy in carbon stock estimates can have implications in conservation and restoration by contributing to the identification of optimal and priority restoration sites, particularly the mangrove-recolonized fishponds. Elemental analyzers will still provide the most accurate OC values for SCS estimation. But in the absence of elemental analyzers, using the LOI method with our proposed equations can be the most suitable and accurate alternative in estimating the SCS of Philippine mangroves.

Data availability

Data will be made available upon request.

Abbreviations

OC:

Organic carbon

CS:

Carbon stock

SOC:

Soil organic carbon

TECS:

Total ecosystem carbon stock

OM:

Organic matter

SCS:

Soil carbon stock

References

  1. Macreadie P, Anton A, Raven J, Beaumont N, Connolly R, Friess D, et al. The future of blue carbon science. Nat Commun. 2019. https://doi.org/10.1038/s41467-019-11693-w.

    Article  Google Scholar 

  2. Murdiyarso D, Purbopuspito J, Kauffman J, Warren MW, Sasmito S, Donato D, et al. The potential of Indonesian mangrove forests for global climate change mitigation. Nat Clim Chang. 2015;5:1089–92. https://doi.org/10.1038/nclimate2734.

    Article  CAS  Google Scholar 

  3. Azman M, Sharma S, Hamzah M, Zakaria R, Palaniveloo K, MacKenzie R. Total ecosystem blue carbon stocks and sequestration potential along a naturally regenerated mangrove forest chronosequence. For Ecol Manage. 2023. https://doi.org/10.1016/j.foreco.2022.120611.

    Article  Google Scholar 

  4. Global Mangrove Alliance. Integrating Mangrove Ecosystems into NDCs; 2022. https://www.mangrovealliance.org/wp-content/uploads/2022/11/Global-Mangrove-Watch_NDC-Policy-Brief_v5.pdf. Accessed Feb 2024.

  5. Sasmito S, Taillardat P, Clendenning J, Cameron C, Friess D, Murdiyarso D, et al. Effect of land-use and land-cover change on mangrove blue carbon: a systematic review. Glob Change Biol. 2019;25:4291–302. https://doi.org/10.1111/gcb.14774.

    Article  Google Scholar 

  6. Alongi D. Global significance of mangrove blue carbon in climate change mitigation. Sci. 2020;2(3):67. https://doi.org/10.3390/sci2030067.

    Article  Google Scholar 

  7. Ahmed S, Sarker S, Kamruzzaman M, Ema J, Naabeh C, Cudjoe E, et al. How biotic, abiotic, and functional variables drive belowground soil carbon stocks along stress gradient in the Sundarbans Mangrove Forest? J Environ Manage. 2023. https://doi.org/10.1016/j.jenvman.2023.117772.

    Article  Google Scholar 

  8. Donato D, Kauffman J, Murdiyarso D, Kurnianto S, Stidham M, Kanninen M. Mangroves among the most carbon-rich forests in the tropics. Nat Geosci. 2011;4(5):293–7. https://doi.org/10.1038/ngeo1123.

    Article  CAS  Google Scholar 

  9. Kauffman J, Donato D. Protocols for the measurement, monitoring and reporting of structure, biomass and carbon stocks in mangrove forests. CIFOR. 2012. https://doi.org/10.1752/cifor/003749.

    Article  Google Scholar 

  10. Breithaupt J, Steinmuller H, Rovai A, Engelbert K, Smoak J, Chambers L, et al. An improved framework for estimating organic carbon content of mangrove soils using loss-on-ignition and coastal environmental setting. Wetlands. 2023. https://doi.org/10.1007/s13157-023-01698-z.

    Article  Google Scholar 

  11. Howard J, Hoyt S, Isensee K, Pidgeon E, Telszewski M, editors. Coastal Blue Carbon: Methods for assessing carbon stocks and emissions factors in mangroves, tidal salt marshes, and seagrass meadows. Arlington, Virginia, USA: Conservation International, Intergovernmental Oceanographic Commission of UNESCO, International Union for Conservation of Nature; 2019.

  12. Thorhaug A, Gallagher J, Kiswara W, Prathep A, Huang X, Yap T, et al. Coastal and estuarine blue carbon stocks in the greater Southeast Asia region: seagrasses and mangroves per nation and sum of total. Marine Pollut Bull. 2020. https://doi.org/10.1016/j.marpolbul.2020.111168.

    Article  Google Scholar 

  13. Bunting P, Rosenqvist A, Hilarides L, Lucas R, Thomas N, Tadono T, et al. Global mangrove extent change 1996–2020: global mangrove watch version 3.0. Remote Sens. 2022;14(15):3657. https://doi.org/10.3390/rs14153657.

    Article  Google Scholar 

  14. Corcino R, Gerona-Daga M, Samoza S, Fraga J, Salmo S III. Status, limitations, and challenges of blue carbon studies in the Philippines: a bibliographic analysis. Reg Stud Mar Sci. 2023. https://doi.org/10.1016/j.rsma.2023.102916.

    Article  Google Scholar 

  15. Salmo S, Lovelock C, Duke N. Assessment of vegetation and soil conditions in restored mangroves interrupted by severe tropical typhoon ‘Chan-hom’ in the Philippines. Hydrobiologia. 2014;733(1):85–102. https://doi.org/10.1007/s10750-013-1766-4.

    Article  Google Scholar 

  16. Samson M, Rollon R. Growth performance of planted mangroves in the Philippines: revisiting forest management strategies. PubMed. 2008;37(4):234–40. https://doi.org/10.1579/0044-7447(2008)37[234:gpopmi]2.0.co;2.

    Article  Google Scholar 

  17. Kauffman J, Heider C, Cole T, Dwire K, Donato D. Ecosystem carbon stocks of Micronesian mangrove forests. Wetlands. 2011;31(2):343–52. https://doi.org/10.1007/s13157-011-0148-9.

    Article  Google Scholar 

  18. Osland M, Hughes A, Armitage A, et al. The impacts of mangrove range expansion on wetland ecosystem services in the southeastern United States: current understanding, knowledge gaps, and emerging research needs. Glob Change Biol. 2022;28:3163–87.

    Article  CAS  Google Scholar 

  19. Beck H, Zimmermann N, McVicar T, Vergopolan N, Berg A, Wood E. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci Data. 2018. https://doi.org/10.1038/sdata.2018.214.

    Article  Google Scholar 

  20. Worthington T, Ermgassen P, Friess D, et al. A global biophysical typology of mangroves and its relevance for ecosystem structure and deforestation. Sci Rep. 2020. https://doi.org/10.1038/s41598-020-71194-5.

    Article  Google Scholar 

  21. Lasco R, Pulhin F, Sales R. Analysis of leakage in carbon sequestration projects in forestry: a case study of upper Magat watershed, Philippines. Mitig Adapt Strat Glob Change. 2007;12:1189–221. https://doi.org/10.1007/s11027-006-9059-0.

    Article  Google Scholar 

  22. Quevedo J, Uchiyama Y, Kohsaka R. Progress of blue carbon research: 12 years of global trends based on content analysis of peer-reviewed and ‘gray literature’ documents. Ocean Coast Manag. 2023. https://doi.org/10.1016/j.ocecoaman.2023.106495.

    Article  Google Scholar 

  23. Nesperos V, Visco R, Malabrigo P Jr, Paelmo R, Pleto J. Blue carbon stock of aquasilviculture-based agroforestry system in Infanta, Quezon, Philippines. Int J Agric Technol. 2023;19(4):1967–1714.

    Google Scholar 

  24. Salmo S III, Malapit V, Garcia M, Pagkalinawan H. Establishing rates of carbon sequestration in mangroves from an earthquake uplift event. Biol Lett. 2019. https://doi.org/10.1098/rsbl.2018.0799.

    Article  Google Scholar 

  25. Nelson D, Sommers L. Total carbon, organic carbon, and organic matter. In: Sparks DL, Page AL, Helmke PA, Loeppert RH, Soltanpour PN, Tabatabai MA, Johnston CT, Sumner ME, editors. Methods of soil analysis part 3. Chemical methods. Madison: Soil Science Society of America, American Society of Agronomy; 1996.

    Google Scholar 

  26. Barrientos K, Apolonio J. Species diversity and soil carbon sequestration potential of mangrove species at Katunggan It Ibajay (KII) eco-park in Aklan, Philippines. Prism. 2017;22(1). https://ssrn.com/abstract=3824652.

  27. Gevaña D, Camacho L, Camacho S. Stand density management and blue carbon stock of monospecific mangrove plantation in Bohol, Philippines. Metsanduslikud Uurimused. 2017;66(1):75–83. https://doi.org/10.1515/fsmu-2017-0008.

    Article  Google Scholar 

  28. Abino A, Castillo J, Lee Y. Species diversity, biomass, and carbon stock assessments of a natural mangrove forest in Palawan, Philippines. Pak J Bot. 2014;46(6):1955–62.

    Google Scholar 

  29. Lomoljo-Bantayan N, Tatil W, Dagoc F, Tampus A, Amparado R. Carbon stock assessment of mangrove forests along Macajalar Bay, Misamis Oriental, Philippines. Int J Adv Appl Sci. 2023;10(10):36–45. https://doi.org/10.2183/ijaas.2023.10.004.

    Article  Google Scholar 

  30. Castillo J, Apan A, Maraseni T, Salmo S III. Soil C quantities of mangrove forests, their competing land uses, and their spatial distribution in the coast of Honda Bay, Philippines. Geoderma. 2017;293:82–90. https://doi.org/10.1016/j.geoderma.2017.01.025.

    Article  CAS  Google Scholar 

  31. Castillo J, Breva L. Carbon stock assessment of four mangrove reforestations/plantation stands in the Philippines. Proceedings of the 1st ASEAN Congress on Mangrove Research and Development, Manila. 2012.

  32. Stankovic M, Mishra A, Rahayu Y, Lefcheck J, Murdiyarso D, Friess D, et al. Blue carbon assessments of seagrass and mangrove ecosystems in South and Southeast Asia: current progress and knowledge gaps. Sci Total Environ. 2023. https://doi.org/10.1016/j.scitotenv.2023.166618.

    Article  Google Scholar 

  33. Ouyang X, Lee J. Improved estimates on global carbon stock and carbon pools in tidal wetlands. Nat Commun. 2020. https://doi.org/10.1038/s41467-019-14120-2.

    Article  Google Scholar 

  34. Salmo S III. Assessment of typhoon impacts and post-typhoon recovery in Philippine mangroves lessons and challenges for adaptive management. In: Sidik F and Friess D, editors. Dynamic sedimentary environments of mangrove coasts. Amsterdam: Elsevier; 2021. p. 539–62.

    Chapter  Google Scholar 

  35. Kristensen E, Bouillon S, Dittmar T, Marchand C. Organic carbon dynamics in mangrove ecosystems: a review. Aquat Bot. 2008;89(2):201–19. https://doi.org/10.1016/j.aquabot.2007.12.005.

    Article  CAS  Google Scholar 

  36. Kauffman J, Adame M, Arifanti V, Schile-Beers L, Bernardino A, Bhomia R, Donato D, et al. Total ecosystem carbon stocks of mangroves across broad global environmental and physical gradients. Ecol Monogr. 2020. https://doi.org/10.1002/ecm.1405.

    Article  Google Scholar 

  37. Aliño P, Gomez E. Philippine coral reef conservation: Its significance to the South China Sea. In: Yamazato K, Ishijima S, Sakihara S, Taira H, Shimabukuro Z, Teruya F, Nishihira F. eds. Development and conservation in the Asia-Pacific Region. Proceedings of the Regional Conference of East-West Center Association: Regional Development in the 21st Century, Think Globally, Act Locally: 1993. Okinawa, Japan: The East-West Center Association Okinawa Chapter. p. 222–229

  38. Baloloy A, Martinez K, Blanco A, Neri M, Ticman K, Burgos D, et al. Mapping multi-decadal mangrove forest change in the Philippines: vegetation extent and impacts of anthropogenic and climate-related factors. In: Filho WL, Kovaleva M, Alves F, Abubakar IR, editors., et al., Climate change strategies: handling the challenges of adapting to a changing climate. Cham: Springer Nature Switzerland; 2023. p. 217–48.

    Chapter  Google Scholar 

  39. Garcia K, Malabrigo P Jr, Gevaña D. Philippines’ mangrove ecosystem: status, threats and conservation. In: Faridah-Hanum I, Latiff A, Hakeem K, Ozturk M, editors. Mangrove ecosystems of Asia. New York: Springer; 2013. p. 81–94.

    Google Scholar 

  40. Pendleton L, Donato D, Murray B, Crooks S, Jenkins W, Sifleet S, et al. Estimating global “Blue Carbon” emissions from conversion and degradation of vegetated coastal ecosystems. PLoS ONE. 2012. https://doi.org/10.1371/journal.pone.0043542.

    Article  Google Scholar 

  41. United Nations Framework Convention on Climate Change (UNFCCC). Philippines - NDC. 2021 Apr. https://unfccc.int/sites/default/files/NDC/2022-06/Philippines%20-%20NDC.pdf. Accessed February 2024.

  42. Pham T, Thi T. Incorporating blue carbon into nationally determined contributions: Current status, opportunities and challenges of 13 Asia-Pacific countries. 2019. https://doi.org/10.1752/cifor/007554.

  43. Sasmito S, Basyuni M, Kridalaksana A, Saragi-Sasmito M, Lovelock C, Murdiyarso D. Challenges and opportunities for achieving sustainable development goals through restoration of Indonesia’s mangroves. Nat Ecol Evol. 2023;7(1):62–70. https://doi.org/10.1038/s41559-022-01926-5.

    Article  Google Scholar 

  44. Mariano H, Aguilos M, Dagoc F, Sumalinab B, Amparado R Jr. Abandoned fishpond reversal to mangrove forest: will the carbon storage potential match the natural stand 30 years after reforestation? Forests. 2022;13(6):847. https://doi.org/10.3390/f13060847.

    Article  Google Scholar 

  45. Matsui N, Suekuni J, Nogami M, Havanond S, Salikul P. Mangrove rehabilitation dynamics and soil organic carbon changes as a result of full hydraulic restoration and re-grading of a previously intensively managed shrimp pond. Wetlands Ecol Manage. 2010;18(2):233–42. https://doi.org/10.1007/s11273-009-9162-6.

    Article  Google Scholar 

  46. Primavera J, Esteban J. A review of mangrove rehabilitation in the Philippines: successes, failures and future prospects. Wetlands Ecol Manage. 2008;16(5):345–58. https://doi.org/10.1007/s11273-008-9101-y.

    Article  Google Scholar 

  47. Brown B, Fadillah R, Nurdin Y, Soulsby I, Ahmad R. Case study: Community-based ecological mangrove rehabilitation (CBEMR) in Indonesia. Sapiens. 2014; 7. http://sapiens.revues.org/pdf/1589.

  48. Soeprobowati T, Sularto R, Hadiyanto H, Puryono S, Rahim A, Jumari J, et al. The carbon stock potential of the restored mangrove ecosystem of Pasarbanggi, Rembang, Central Java. Marine Environ Res. 2024. https://doi.org/10.1016/j.marenvres.2023.106257.

    Article  Google Scholar 

  49. Matatula J, Wirabuana P, Yasin E, Mulyana B. Species composition and carbon stock of rehabilitated mangrove forest in Kupang District, East Nusa Tenggara, Indonesia. Environ Res Eng Manag. 2023;79(3):24–34. https://doi.org/10.5755/j01.erem.79.3.33123.

    Article  Google Scholar 

  50. Soboka D, Yimer F. Restoration of degraded lands for carbon stock enhancement and climate change mitigation: the case of Rebu watershed, Woliso Woreda, Southwest Shoa, Ethiopia. J Degrade Min Lands Manag. 2022;9(2):3387–96. https://doi.org/10.1524/jdmlm.2022.092.3387.

    Article  Google Scholar 

  51. Zanini A, Mayrinck R, Vieira S, De Camargo P, Rodrigues R. The effect of ecological restoration methods on carbon stocks in the Brazilian Atlantic Forest. For Ecol Manag. 2021. https://doi.org/10.1016/j.foreco.2020.118734.

    Article  Google Scholar 

  52. Duncan C, Primavera J, Pettorelli N, Thompson J, Loma R, Koldewey H. Rehabilitating mangrove ecosystem services: a case study on the relative benefits of abandoned pond reversion from Panay Island, Philippines. Marine Pollut Bull. 2016;109(2):772–82. https://doi.org/10.1016/j.marpolbul.2016.05.049.

    Article  CAS  Google Scholar 

  53. Adame MF, Zakaria RM, Fry B, Chong VC, Then YHA, Brown CJ, Lee SY. Loss and recovery of carbon and nitrogen after mangrove clearing. Ocean Coast Manag. 2018;161:117–26.

    Article  Google Scholar 

  54. Van Hieu P, Dung LV, Tue NT, Omori K. Will restored mangrove forests enhance sediment organic carbon and ecosystem carbon storage? Region Stud Marine Sci. 2017;14:43–52.

    Article  Google Scholar 

  55. Chaikaew P, Chavanich S. Spatial variability and relationship of mangrove soil organic matter to organic carbon. Appl Environ Soil Sci. 2017. https://doi.org/10.1155/2017/4010381.

    Article  Google Scholar 

  56. Dung LV, Tue NT, Nhuan MT, Omori K. Carbon storage in a restored mangrove forest in Can Gio mangrove forest park, Mekong Delta, Vietnam. For Ecol Manag. 2016;380:31–40.

    Article  Google Scholar 

  57. Nguyen PT, Tue NT, Quy TD, Thai ND. Quantifying organic carbon storage and sources in sediments of Dong Rui mangrove forests, Tien Yen district, Quang Ninh province using carbon stable isotope. Vietnam J Earth Sci. 2016;38(4):317–26.

    Google Scholar 

  58. Phang VXH, Chou LM, Friess DA. Ecosystem carbon stocks across a tropical intertidal habitat mosaic of mangrove forest, seagrass meadow, mudflat and sandbar. Earth Surf Proc Land. 2015;40(10):1387–400.

    Article  CAS  Google Scholar 

  59. Tue NT, Dung LV, Nhuan MT, Omori K. Carbon storage of a tropical mangrove forest in Mui Ca Mau National Park. Vietnam Catena. 2014;121:119–26.

    Article  Google Scholar 

  60. Castillo J, MacKenzie R, Manahan J, Castillo J. Monitoring the sediment surface elevation change across a chronosequence of restored stands of tropical mangroves and their contemporary carbon sequestration in soil pool. Forests. 2022;13(2):241. https://doi.org/10.3390/f13020241.

    Article  Google Scholar 

  61. Dimalen F, Rojo M. Carbon stock assessment of a mangrove forest in Cotabato City, Philippines. Journal of Biodiversity and Environmental Science. 2019. https://ssrn.com/abstract=3588114.

  62. Salmo S III, Gianan E. Post-disturbance carbon stocks and rates of sequestration: implications on “blue carbon” estimates in Philippine mangroves. Philipp Sci Lett. 2019;12(2):122–32.

    Google Scholar 

  63. Gevaña D, Pampolina N. Plant diversity and carbon storage of a Rhizophora stand in Verde Passage, San Juan, Batangas, Philippines. J Environ Sci Manag. 2009;12(2):1–10.

    Google Scholar 

Download references

Acknowledgements

We extend our deep gratitude to the local officials, communities, and field guides from the municipal governments of Bongabong, Mansalay, and Roxas (in Oriental Mindoro), the provincial government of Oriental Mindoro, and the municipal government of Prieto Diaz (in Sorsogon). John Ryan Navidad and Jayson Del Rosario assisted in the field sampling and preliminary lab analyses.

Funding

This study was supported by the “Assessment and Comparison of Recovery of Biodiversity and Carbon Sequestration in Philippine Mangroves Among Natural, Replanted, and Naturally-recolonized Mangrove Stands” (PEER ManCoRe) Project funded by the National Academy of Sciences (NAS) and the United States Agency for International Development (USAID) under USAID cooperative agreement number AID-OAA-A-11-00012 through the Partnerships of Enhanced Engagement in Research (PEER) with Grant Number 9–379. This study was also partly funded by the United States Forest Service (USFS) under International Cooperative Agreement number 21-IC-11132762-290 and the Ministry of Education and Culture Republic of Indonesia under the World Class University Program of Universitas Sumatera Utara Year 2022–2023 (Grant no. 20084.1/UN 5.4.17/ TPM/2022).

Author information

Authors and Affiliations

Authors

Contributions

SGS: Conceptualization, writing, revision, analysis, supervision, project administration, funding acquisition SBM: Lab analysis, writing, data collection and organization, data/statistical analysis, visualization PBJ: Lab analysis, writing, data collection and organization, data/statistical analysis, visualization MGD: Sampling, writing CFP: Sampling, lab analysis, writing MWM: Lab analysis, writing MB: Writing, funding acquisition FS: Writing RM: Writing, funding acquisition.

Corresponding author

Correspondence to Severino G. Salmo III.

Ethics declarations

Competing interests

The authors declare 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-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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 http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Salmo, S.G., Manalo, S.P.B., Jacob, P.B. et al. Improving soil carbon estimates of Philippine mangroves using localized organic matter to organic carbon equations. Carbon Balance Manage 19, 31 (2024). https://doi.org/10.1186/s13021-024-00276-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13021-024-00276-y

Keywords