Abstract

Degradation and destruction of mangrove forests in many regions have resulted in the alteration of carbon cycling. Objectives of this study were established to answer the question regarding how much soil organic carbon (SOC) is stored in wetland soils in part of the upper northeastern Gulf of Thailand and to what extent SOC is related to organic matter (OM). A total of 29 soil samples were collected in October 2015. Soil physiochemical analyses followed the standard protocol. Spatial distributions were estimated by a kriging method. Linear regression and coefficient were used to determine the suitable conversion factor for mangrove soils. The results showed that surface soil (0–5 cm) contained higher SOC content as compared to subsurface soil (5–10 cm). Considering a depth of 10 cm, this area had a high potential to sequester carbon with a mean ± standard deviation of %. The spatial variability of OM and SOC revealed that organic matter and carbon decreased with the distance from upstream areas toward the gulf. Based on the assumption that OM is 50% SOC, the conversion factor of 2 is recommended for more accuracy rather than the conventional factor of 1.724.

1. Introduction

Since the Industrial Revolution, there has been a speedy increase in atmospheric carbon dioxide (CO2) concentration. Global levels of CO2 already passed above 400 ppm mark in 2015 [1]. The natural movement of carbon across the atmosphere, vegetation, soils, and the oceans is the key to mitigate climate change due to elevated CO2. Globally, soils store more carbon, about 3.3 times, than the atmosphere and 4.5 times the bionic pool [2]. Among different ecosystems, wetlands often represent the largest carbon pool because of their anoxic conditions and thus play a vital role in carbon cycles. Of the 1,500 Pg C stored in the Earth’s soils, peat-forming wetlands are estimated to contain about 300–600 Pg C [3]. Southeast Asian peatlands account for the most important carbon sink, representing 68.5 Pg C or 77% of global tropical peatlands [4]. Despite accounting for the importance of wetlands to ecosystems and humans, one-third of global wetlands have been lost due to land use change, aquaculture inversion, and wetland degradation. Carbon emissions from mangrove loss are uncertain, but it is estimated that about 10% of all carbon emissions are released by deforestation worldwide [5]. Thailand alone is estimated to have lost approximately 82% of its wetlands [6]. Chonburi Province, Thailand, also has experienced wetland loss as a result of urban sprawl, industrial development, and deforestation. Between 1961 and 2007, its wetland coverage was reduced from 38.2 km2 to 7.8 km2, or by about 80% [7]. Wetland conversion directly affects the size of the soil carbon pool.

Soil organic carbon (SOC) is part of the carbon stored in soil organic matter (OM). The estimation of SOC has followed various techniques from the oldest to the simplest method. The assumption that OM contains 58% carbon was first established by Sprengel in 1826 [8]. Since then, the conversion factor of 1.724 has been repeatedly used as a rough estimate of organic carbon content for universal applications. Other early studies published in German scientific literature assumed that OM holds 58% carbon (conversion factor of 1.724) [9, 10] and 60% carbon (conversion factor of 1.667) [11, 12]. However, these numbers are considered too low in many areas. The application of 1.9 and 2.5 was recommended to convert OM to SOC for soils in the top and subsurface layers, respectively, worldwide [13]. Soil type and environment result in SOC-to-OM variations. In wetland ecosystems, a factor of 1.842 has been determined for most peats [14], 1.88 (advocated 2) for peat in Wales [15], and 2.07 and 2.34 for wetlands in England [16, 17]. A recent study has claimed that a factor of 2, based on the assumption that OM is 50% carbon, is more accurate in almost all cases than the conventional factor of 1.724 [18].

Soil management can manipulate the balance of carbon sink and source processes. When OM is broken down, some of the SOC can be quickly mineralized and converted to CO2. Carbon thus disappears from the soil. The high amount of soil organic matter tends to be limited to the soil surface, probably at a depth of 5 to 10 cm [19]. Organic matter often binds to fine particles, particularly clay, and the binding processes further prevent soil carbon loss [20]. The study of SOC and its relationship to OM is important; yet there are research gaps to be filled regarding mangrove soils in Thailand. It is still questionable how much SOC is stored in the wetland soils and to what extent SOC is related to OM. The aim of this study is to spatially investigate the SOC content in mangrove areas at surface and subsurface layers, as well as the association of SOC-to-OM conversion. The result of the conversion factor from this study is useful for estimating either SOC or OM in similar environments. The specific objectives of this study are threefold:(1)Investigating the concentrations of SOC in surface soil (0–5 cm) and subsurface soil (5–10 cm).(2)Assessing spatially explicit OM and SOC contents (0–10 cm).(3)Identifying the conversional function between OM and SOC in mangrove soils (0–10 cm).

2. Materials and Methods

2.1. Site Description

Estuarine mangrove forests along the coast of the Gulf of Thailand are influenced by the open-coastal environment during high tide and fresh water flow during low tide. The study area is part of the Nature Education Center for Mangrove Conservation and Ecotourism located in Chonburi Province (13°20′37.05′′N, 100°56′34.83′′E) on the eastern coast of the Gulf of Thailand. The center was established in 2001 for mangrove conservation and education purposes. The conservation area covers approximately 480,000 m2, but the focus area of interest is limited to the area within and around the wooden walkway due to the constraints of area accessibility. The area within the walkway covers 113,240 m2 and extends to 220,000 m2 for the study area. The water channel receives runoff from residential communities and industries upstream and then drains through the mangrove conservation area before being discharged into the Gulf of Thailand. Chonburi has a tropical savanna climate [21] and temperatures rise in April, with the average daily maximum at 35.2°C (95.4°F). The monsoon season, May through October, has heavy rain and somewhat cooler temperatures during the daytime but nights remain warm. Rich variations of fauna such as mud skipper, flying fox, fiddler crab, sesarma mederi, oyster, cockle, banana shrimp, and snake are commonly found. The main mangrove dominance is Avicennia alba Bl. due to their adaptive characteristics in tolerating a saline environment. Other plant species found across the conservation area include Rhizophora apiculata Blume, Rhizophora mucronata, Sonneratia griffithii, Sonneratia alba, Xylocarpus granatum, and Xylocarpus moluccensis. The presence of poor mangrove growth was near the entrance of the walkway and appears in the southeast center of the study area.

2.2. Sampling Design

The sampling site was marked by a handheld GPS (Garmin GPSMAP® 62sc model) on site. The boundary of the area was mapped using Google Earth imaginary layout [22]. The geographic coordinate system used was WGS84 before being projected to the UTM_Zone_47N system for consistent data analysis. The spatial sampling sites were planned on the grid-based sampling design with a cell size of . From an initial grid design providing 35 soil samples covering the entire area, a total of 29 soil samples from the surface (0–10 cm) across the study area were available for analysis due to limited accessibility (Figure 1). Soil samples were collected in October 2015. All the samples were collected during low tides; thus a hand collection technique was implemented. To collect samples, a 10 × 5 cm2 soil core was used from each designated location within a square grid. The collected sediments were cut into two layers: surface (0–5 cm) and subsurface (5–10 cm). All samples were kept at 4°C at the sampling moisture content for laboratory processing.

2.3. Analysis of Soil Characteristics

The physiochemical sediment properties, that is, texture, pH, salinity, electrical conductivity (EC), total dissolved solid (TDS), and soil texture, were analyzed. After air-dried samples were prepared, the percentages of sand, silt, and clay in the sediments were determined by hydrometer testing [23] and then compared with a soil texture triangle to identify soil texture [24]. The ratio of soil to water was 1 : 2.5 as suggested by Jackson [25]. The measurement of pH was determined by a pH meter (Vetus, pH-009 (III), USA). Salinity, EC, and TDS were analyzed by a portable handheld EC meter (Hach, sensION 156, USA) [25].

2.4. Assessment of Organic Matter

The loss on ignition (LOI) method was employed to analyze organic matter content at a 10 cm depth. The wet soil samples were completely oven-dried at a temperature of 105°C until the weight of dried soils remained constant. The dried soil samples were then burned in a furnace at 550°C for 4 hours before being weighed. Organic matter content was calculated from the soil sample weight before and after burning by using the following equation:

2.5. Assessment of Soil Organic Carbon

The soil samples were freeze-dried overnight in a freeze dry/shell system (FreeZone6, Labconco, USA) and then ground to a fine powder. The fractions passing a 2-mm sieve were stored in a polypropylene centrifuge tube. Total carbon (TC) and soil inorganic carbon (SIC) were measured by the Shimadzu gas analyzer (TOC-V CPH) using separate gas analysis procedures. The CO2 evolution using tiny ball-milled samples with gas combustion was applied at different temperatures. The temperature was set at 900°C for the combustion process to measure TC, while adding phosphoric acid (H3PO) reaction and setting the temperature at 200°C to measure IC [26, 27]. The organic carbon can then be derived by subtracting SIC from TC.

2.6. Statistical and Geostatistical Analyses

Descriptive statistics were performed to establish physiochemical properties information using SPSS 22.0. An ordinary kriging approach was used to estimate the spatial SOC distributions in both surface and subsurface layers. Three main steps in the kriging included exploring measured data to characterize spatial continuity, developing a semivariogram model for investigating spatial autocorrelation, and estimating unknown values based on optimized weighted averages of defined neighboring locations [28]. The relationship between OM and SOC was carried out by a Pearson’s correlation coefficient and a simple linear regression in statistical software (version 3.2.0).

3. Results

3.1. Soil Organic Carbon and Physiochemical Characteristics at 0–5 cm and 5–10 cm Depths

Sandy clay loam covered the majority of the surface soil, followed by sandy loam and loam. On the other hand, clay texture was dominant in the subsurface soil, followed by sandy clay loam and clay loam. Low salinity with mean values of 2.04 ppt in surface soil and 2.22 ppt in subsurface soil was influenced by strong water inflow from the estuary to the seaside because the data collected in this study was taken during low tides. These values were consistent with average EC values of 3.66 mS/cm and 3.96 mS/cm which indicated moderately saline condition and sea water was diluted by the fresh water. Its salinity, the mixture of saline and fresh water, ranged from 1.00 to 3.52 ppt. According to the salinity classification of Levinton [29], mangrove soil in this area was considered under brackish and slightly alkaline condition with an average pH above 7. Descriptive statistics of SOC content and related parameters are summarized in Table 1.

The measured values of SOC were equally alike to the amount of TC, indicating very low inorganic carbon forms in this area. The surface soil exhibited a higher SOC content than the subsurface layer. The results show that, vertically, SOC content decreased with increasing depth. The mean observation with standard deviation SOC was % and % for the surface and subsurface soils, respectively. Thus surface soils with higher productive mangroves in this area retain a greater amount of SOC input.

Local topography and soil drainage are considered important drivers in determining carbon storage. Since the mangrove area in this study was generally flat, the erosion process to transport carbon from topsoil to the bottom of the slopes was difficult. The water-logged condition in the mangrove forest could promote depositional rates and inhibit decomposition rates [30]. It can be observed that the disturbance of soil surface from tillage, dredging, or other aquaculture activities was extremely low in the conservation area compared to agricultural space. Less soil disturbance contributed to a net positive SOC storage of estuarine wetlands [31]. The different amount of carbon deposits along depths was largely influenced by vegetation. In an area with similar decomposition rates, carbon inputs greatly depend on primary productivity. This may be explained by the size of fraction SOC inputs and decomposition rates. Since a large fraction comes from root turnover, the variability of root inputs is an influential factor in carbon storage [32].

3.2. Carbon Storage (0–10 cm)

At a depth of 10 cm, SOC varied across vegetation species and areas. The Nature Education Center for Mangrove Conservation and Ecotourism in this study, heavily dominated by Avicennia alba Bl., stored SOC at about 5.59% and OM 11.12%. Observed organic matter ranged from 6.14 to 17.34% with a mean ± standard deviation of %. Soil organic matter ranged from 2.26 to 9.80% with a mean ± standard deviation of % (Table 2).

3.3. Estimated Spatial Distribution of Organic Matter and Soil Organic Carbon (0–10 cm)

The close values between mean and median with low skewness described the normal central tendency of the data. Log transformation was not required to perform krigged estimates. The total spatial distribution of OM showed a minimum of 9.19%, maximum of 12.47%, and mean ± standard deviation of %. The SOC variability demonstrated a minimum of 3.32%, maximum of 7.91%, and mean ± standard deviation of %. The estimated OM and SOC derived by ordinary kriging are shown in Figure 2. A spherical model was applied to fit the experimental semivariograms of OM and SOC with a cell size of 3 × 3 m2. The semivariograms of the data were found to vary between OM and SOC. Organic matter had longer correlation range (945.84 m) than SOC (322.49 m). The nugget and sill values of OM were 3.65 and 6.30, whereas nugget and sill of SOC were 2.61 and 6.34. Moderate spatial autocorrelation was found for both OM (0.58) and SOC (0.41) as indicated by the nugget-to-sill ratio. A small ratio in SOC informs the stronger spatial dependency as compared to OM. The large spatial autocorrelation range of 945.8 m of OM explained the broader homogenous scale compared to a shorter range of 322.5 m of SOC. Organic matter and SOC showed uniform spatial patterns, with high content near the land and low content toward the Gulf of Thailand. The tidal export of organic matter tended to be impeded by the mangrove root system and forest structure. The findings designated that plant debris and riverine influx of dissolved and particulate organic matter upstream supply a substantial degree of organic loading near the upper intertidal area.

3.4. SOC-to-OM Conversion

The average SOC content of OM in soils at a 10-cm depth was approximately 50%. The relationship between SOC and OM is illustrated by a scatter plot with a linear conversion in Figure 3, the regression formula beingwhere represents dried-weight soil organic carbon (%) measured by gas combustion analysis and represents dried-weight soil organic matter (%) measured by loss on ignition method.

In Figure 3(a), the relationship between OM and SOC of soils in this study area was weak but statistically significant (,   value < 0.001) and Pearson’s correlation coefficient was 0.50 for all the samples. When considering 17 samples near the ocean (Figure 3(b)), a model was moderately fitted with statistical significance (,   value < 0.001), and Pearson’s correlation coefficient was 0.68. When considering 12 samples near the land, the regression model showed a very weak relationship and was not statistically significant (,   value = 0.74) and Pearson’s correlation coefficient was as low as 0.10. The different relationships of OM and SOC detected between samples close to the ocean and far from the ocean were probably due to the natural settings such as inputs of OM, anaerobic environment, and influence of tides.

An average ratio of spatial OM : SOC was across the entire study area (Figure 4). Higher ratios were observed in the direction toward the ocean and findings from spatial analysis were compatible with results from our point-based data. Mangrove soils comprised of organic matter holding 50% carbon identified a conversion factor of 2. Even though the practical use of the 1.724 factor found at the beginning of nineteenth century has been established through authority and has been adopted in many publications, the results from this study were similar to those of other investigators who directly and indirectly indicated the value of 1.724 to be too low to estimate the amount of OM of mangrove and/or wetland soils [1518].

4. Discussion and Comparisons

Looking back at Table 2 we can see that the results were comparatively equal as compared to carbon content under dominance of Rhizophora mangle (mean SOC %) but were lower when compared to dominance of Avicennia schaueriana (mean SOC %) in Brazil (Table 3) [33]. In Indonesia, soils under Avicennia forest and Ceriops forest showed mean ± SD SOC contents of % and % at a depth of 20 cm, respectively. An overview of SOC content studies at different depths and mangrove forests are shown in Table 3.

Of note also is Howard’s study [16], where a regression model was expressed across various soil types as . However, the altered factors and amount of SOC in OM relied on type of soils being from 1.77 to 1.93 (51.92–56.39%) in moors, 1.90 to 1.95 (51.29–52.50%) in acid peats, 1.92 (52.19–52.20%) in fen peats, 1.81–1.83 (48.36–50.87%) in alluvial soils, and 1.97–2.07 in mulls. Published results early in the twentieth century showed that 1.8 was a suitable OM-to-SOC factor for the marine sediment [37]. Later, the simulated model of OM and nutrients in mangrove wetland soils located in the Shark River estuary of south Florida showed OM : SOC ratios from 1.81 to 2.10 with a mean value of 1.98 based on four mangrove sites [38]. Robinson et al. [15], Ponomareva and Plotnikova [39], and Pribyl [18] advocated 2 as the recommended factor for accurate conversion. The SOC-to-OM conversion factor of 2 is, thus, considered to provide better accuracy for estimating the carbon content in mangrove environments or similar types of wetland. Our results coincided with the study from Everglades National Park, Florida [38], in which OM concentrations decreased from 82% to 30% in upstream locations down to marine sites.

The recent study mentioned two sources of variability in any estimated factor: first, the different methods used to measure OM and SOC and, second, the diversity in soil composition [18]. In this study, the methods used to estimate OM and SOC were consistent; thus natural factors played a dominant role. The natural setting of the coastal area can reflect differences in carbon content of OM that, in this study, showed higher conversion factors of SOC to OM in the mangrove environment further from shore. In nature, mixed origins of organic matter in mangrove soil can be very diverse and complex. The near-ocean area showed a high density of mangrove plants that can gain soil OM from fallen detritus decomposition. On the contrary, near-shore areas potentially receive OM sources from land via human activity as it is closer to the community. Therefore, the rate and source of organic matter input affect microbial reactions as well as nutrient availability [40]. Continued production and slow decomposition can lead to very large OM content in soil with long periods of water saturation, but raised OM decomposition occurs in more aerated or less anoxic conditions. However, the speed and ease of carbon mineralization depend on the OM fraction in which it resides; for example, humus-carbon mineralizes slowly when compared to plant residue, particulate OM, and soil microbial biomass [41]. The tidal process is one factor driving a great influx of SOC content and turnover in mangrove forests [42]. Coastal mangroves act as a barrier to any fresh water and tidal flows; this buffer zone often develops suspended matter including OM and SOC in the directions parallel to the coast. Areas further from shore with a high density of mangroves have a vertical mixing of water caused by strong waves and tides which tend to circulate suspended matter; thus OM and SOC usually remain in this zone. Estimates of OM : SOC in soil vary from site to site and may be affected by other factors such as types of vegetation cover, clay in the soil, degree of decomposition, and soil composition, which are out of the study scope.

5. Conclusions

The Nature Education Center for Mangrove Conservation and Ecotourism under dominance of mixed natural and planted Avicennia alba Bl. showed a great potential to sequester carbon in soils with high amounts of SOC stored in the top 10 cm. Surface soil (0–5 cm) contained an average ± standard deviation in SOC of % which is greater when compared to subsurface soil (5–10 cm) with %. Considering the ecosystem in the study area, the origins of the carbon inputs were derived from primary productivity such as leaf debris and plant roots. The spatial variability of OM and SOC demonstrated the decreasing pattern from land to sea or upper to lower tidal zones. Besides the influence of the tides, we summarized that a substantial load of dissolved and particulate organic matter comes from canals and the community upstream. In terms of the strength of the spatial autocorrelation, the nugget-to-sill ratio of 0.41 for SOC was stronger than 0.58 for OM yet indicated moderate spatial dependency for both variables. Organic matter showed homogenous spatial variability due to the long spatial autocorrelation range. This study contradicts the assumption that has been repeatedly used for more than a century that OM is 58% organic carbon and that the conversion factor is 1.724 (van Bemmelen factor). The use of a low factor can cause sequential errors in estimating OM when multiplying carbon by 1.724. Organic matter in mangrove soils in this study contained about 50% organic carbon with the conversion factor of 2. The use of the factor of 2 in different types of wetlands was also recommended by considerable supporting evidence. However, this number is not for universal purposes but rather for rough estimation. The conversion factor may vary depending on vegetation cover, temperature, soil type, soil depth, and the technique used to analyze OM and SOC.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Funding research was provided by Grants for Development of New Faculty Staff, Ratchadaphiseksomphot Endowment Fund, Chulalongkorn University.