Effect of Restoration Actions on Organic Carbon Pools in the Lagoon—Delta Ci é naga Grande de Santa Marta, Colombian Caribbean

: Mangroves provide multiple ecosystem services and are essential for mitigating global warming owing to their capacity to store large carbon (C) stocks. Due to widespread mangrove degradation, actions have been implemented to restore them worldwide. An important representative case in Colombia is the Ci é naga Grande de Santa Marta’s restoration plan. This management intervention focused on restoring the natural hydrological functioning after massive mangrove mortality (~25,000 ha) due to soil hyper-salinization after river water input from the Magdalena River was eliminated. A partial recovery occurred during subsequent years, and hydrological management is still being implemented today. To understand how the degradation and subsequent management have affected mangrove C stocks, we compared C stocks in stands with different intervention levels reﬂected in their current forest structure. We found that the total C stock (398–1160 Mg C ha − 1 ) was within the range measured in other neotropical mangroves without vegetation deterioration. The aboveground C was signiﬁcantly higher in the stands where hydraulic connectivity was restored. By contrast, the belowground C was higher in the stands with low hydraulic connectivity due to channel clogging and a lack of sufﬁcient maintenance. Our results show that hydrological management measures inﬂuenced above- and belowground C stocks, even at a 2 m depth. In addition, a strong indirect relationship useful for estimating carbon content from organic matter content was found.


Introduction
Mangroves provide numerous ecosystem services, including their extraordinary capacity to store and retain carbon (C) [1,2]. These ecosystems can store up to five times the amount of C present in tropical forests [3]. Kauffman and Bhomia [4] and Bindoff et al. [5] showed that, on a global scale, mangrove C stocks vary between 50 and 2200 Mg C ha −1 . Mangroves store large amounts of C owing to several mechanisms. The high productivity of mangroves is a source of autochthonous organic debris, and their ability to reduce water flow energy favors the accumulation of autochthonous and allochthonous particles by increasing sedimentation and reducing resuspension [6,7]. Additionally, their soils saturated with water remain in an anaerobic state, preventing the degradation of organic matter, and promoting its accumulation [1,8,9]. La Ciénaga Grande de Santa Marta (CGSM) is an estuarine lagoon system located on the central Caribbean coast of Colombia (10 • 30 -11 • 15 N, 74 • 15 -74 • 45 W). Geomorphologically, it is part of the deltaic plain of the Magdalena River [23], the largest river in Colombia with an annual average water discharge of 7200 m 3 s −1 and 144 × 106 t yr −1 of sediment [24]. Limited to the east by the Sierra Nevada de Santa Marta (the highest coastal mountain in the world), to the north by the Caribbean Sea, from which it is separated by the Salamanca barrier island, and to the west and south west by the Magdalena River floodplain [18,25,26] (Figure 1). The delta and the CGSM rests on a coastal plain, which was formed by marine and fluvial sedimentary deposition during the Holocene transgression and regression phases (beginning about 10,000 B.P.), while the surrounding formations date back to the Cenozoic era. The lagoon complex (1321 km 2 ) comprises the Ciénaga Grande (450 km 2 ), the Ciénaga de Pajarales (120 km 2 ), several smaller lagoons, creeks, and channels (150 km 2 ), and mangrove swamps [22,23].
The climate of the area is tropical semi-arid, with a dry period from December to March and a rainy period for the rest of the year, periods greatly influenced by years with El Niño or La Niña events [27]. The maximum precipitation occurs in October, and in general, during the year, the region has a water deficit due to evaporation (1953.8 mm yr −1 ) exceeding precipitation values (mean = 807 mm yr −1 ) [26,27]. Given the wide variety of environments in the delta-lagoon system covering approximately 130,000 ha [22], vegetation in the area ranges from tropical thorn bush and wetland forest stands growing near the rivers in soils that are periodically flooded, to wetland and halophytes herbaceous vegetation and mangroves; the latter considered the most important vegetation type in the area [28,29]. The climate of the area is tropical semi-arid, with a dry period from December to March and a rainy period for the rest of the year, periods greatly influenced by years with El Niño or La Niña events [27]. The maximum precipitation occurs in October, and in general, during the year, the region has a water deficit due to evaporation (1953.8 mm yr −1 ) exceeding precipitation values (mean = 807 mm yr −1 ) [26,27]. Given the wide variety of environments in the delta-lagoon system covering approximately 130,000 ha [22], vegetation in the area ranges from tropical thorn bush and wetland forest stands growing near the rivers in soils that are periodically flooded, to wetland and halophytes herbaceous vegetation and mangroves; the latter considered the most important vegetation type in the area [28,29].
The CGSM, with a primary production value of 5.76 g C m −2 d −1 [19], is one of the most productive coastal lagoons in the world [20]. It is the economic support of seven fishing villages with a total population of close to 20,000 people. Due to its ecological and socioeconomic relevance, the area includes several special management zones such as the Parque Isla de Salamanca, the CGSM Fauna and Flora Sanctuary, and the Exclusive Reserve Zone (the CGSM main coastal lagoon and the Complejo de Pajarales lagoons). It was also declared a Ramsar Wetland, a Biosphere Reserve, and an area of international importance for bird conservation (IAICA) [30,31].
Mangrove forests in CGSM covered 51,150 ha until 1956, but a continuous die-off that lasted 40 years reduced the forests by approximately 24,624 ha. This massive death of mangroves was due to a series of physical alterations made worse by climatic effects. The physical alterations were caused by (1) the construction of the Barranquilla-Ciénaga The CGSM, with a primary production value of 5.76 g C m −2 d −1 [19], is one of the most productive coastal lagoons in the world [20]. It is the economic support of seven fishing villages with a total population of close to 20,000 people. Due to its ecological and socioeconomic relevance, the area includes several special management zones such as the Parque Isla de Salamanca, the CGSM Fauna and Flora Sanctuary, and the Exclusive Reserve Zone (the CGSM main coastal lagoon and the Complejo de Pajarales lagoons). It was also declared a Ramsar Wetland, a Biosphere Reserve, and an area of international importance for bird conservation (IAICA) [30,31].
Mangrove forests in CGSM covered 51,150 ha until 1956, but a continuous die-off that lasted 40 years reduced the forests by approximately 24,624 ha. This massive death of mangroves was due to a series of physical alterations made worse by climatic effects. The physical alterations were caused by (1) the construction of the Barranquilla-Ciénaga highways, on the Salamanca Bar between 1956 and 1960, and Palermo-Sitio Nuevo, marginal to the Magdalena River, on the eastern bank, in the 1970s, without considering the communication between the river and the CGSM; (2) the construction of dikes and embankments to prevent the overflowing of the Magdalena River and the flow of freshwater from the river to the lagoon complex; (3) sedimentation/obstruction of the culverts coming from the Magdalena River and the mouths of the Sierra Nevada de Santa Marta (SNSM) rivers as a consequence of erosion due to deforestation in the basins; (4) a decrease in the flow of the SNSM rivers due to the collection of water for the irrigation districts; and (5) soil compaction for agriculture and livestock [32]. The main direct consequence of these physical alterations was the breakdown of the hydrological balance of the entire system due to the interruption of freshwater flows from the rivers to the lagoon system and of marine and estuarine water from and to the Caribbean Sea and the swamps across the Salamanca Bar. On the other hand, El Niño southern oscillation (ENSO) triggers lower precipitation along the Magdalena River catchment during El Niño years and higher precipitation during La Niña years that influence CGSM salinity [24,27,33,34]. The loss of hydrologic connectivity in CGSM occurred during an average neutral ENSO condition; this combination generated a high mangrove coverage loss that accumulated for many years in this dry climate region [7].
In an area as semi-arid as CGSM, where evapotranspiration far exceeds precipitation, mangrove trees with 20 m heights and 65 cm in diameters [35], indicate a freshwater subsidy. This case was mainly from the Magdalena River, on which the function of the system largely depends. The hydrological imbalance in the entire system is mainly reflected in the extreme hyper-salinization of the mangrove soils [18,26,36] and in the progressive increase in the salinity of the water bodies. Much of the soils exhibited interstitial water salinities greater than 100 for 8 months of the year. The soil hyper-salinization was identified as the main cause of the massive mortality of the mangroves, and it was shown that there is an inverse relationship between the salinity of the soil and the basal area of the mangrove in some sites [36]. The rate of mangrove loss gradually increased from 1.75 km 2 yr −1 from 1956 to 1968, to 9.8 km 2 yr −1 from 1968 to 1987, to 13.3 km 2 yr −1 from 1987 to 1993, and to 18,4 km 2 yr −1 between 1993 and 1995 [37,38]. The loss of mangroves caused dramatic consequences in the CGSM, including the loss of biodiversity associated with the forest. While in conserved mangrove areas an average of 29 species of birds, 32 species of fish and 50 species of invertebrates were recorded, in dead mangrove only 19 species of birds, 14 species of fish and 5 species of invertebrates were found [39]. In addition to extensive mangrove mortality, several fish kills caused by low oxygen concentrations and eutrophication have occurred during recent years [40].
With the intention of restoring the historical hydrological conditions, in 1992 the Project for the Rehabilitation of the CGSM (ProCiénaga) was initiated, with the reopening of 5 pre-existing channels that connected the Magdalena River with the mangrove areas, 3 of which have a direct effect on the areas of high mangrove mortality: Clarín, Renegado and Aguas Negras [18] (Figure 1). With the reestablishment of water inflows, the mangrove began to regenerate naturally, and by 2013 it was estimated that approx. 14,000 ha. had recovered [41]. After the ProCiénaga Project, the local environmental authority took over the maintenance and cleaning of canals, which has been done sporadically. Since 2014, mangrove dieback episodes have occurred due to the lack of canal maintenance and a prolonged drought associated with the strong El Niño event of 2015 and 2016 [31]. Within the measures taken to mitigate the new mangrove deterioration, a hydraulic recovery campaign was executed between 2016 and 2019 [42]. The campaign began in 2016 with the dredging of pre-existing channels, including those that had been reopened with the ProCiénga Project. In the case of the Clarín channel, the recovery included, in addition to dredging, the clearing of weeds (aquatic macrophytes) and the construction of lateral channels of the fishbone type. In 2017, channels were reopened in the northwestern sector of the Ciénaga (Caleta del Tambor, Bristol, Covado, Machetico, Covado 2); in 2018, the Clarín channel had its first two kilometers newly dredged in the area of Vía Parque Isla de Salamanca, and in 2019, maintenance was performed at the mouth of rivers that reach the eastern part of the Ciénaga Grande [42].
For this study, we selected four sites with different intervention levels located on the west side of the CGSM. RIN is in the northeastern corner of the main water body of the CGSM and close to the mouth of the Clarín channel ( Figure 1). The forest here is composed of Avicennia germinans (L.) L., Laguncularia racemosa (L.) C.F. Gaertn and Rhizophora mangle L., with the former being the dominant species. RIN was not affected by the massive mangrove mortality of the 1970s. Although in this site some large trees have died, the forest structure has been relatively stable over the years, with a basal area (BA) = 34.1 m 2 ha −1 [41]. K22 is located around the middle of the Clarín channel course. This site lost approximately 70% of its forest cover during the massive mortality (BA = 2 m 2 ha −1 by 1990) and recovered after the opening of the freshwater channels (BA = 27.2 m 2 ha −1 by 2013) [41,43]. All three common mangrove species are present. K22 has been directly influenced by the opening of the Clarín channel during the ProCiénaga project and by maintenance dredging operations. Since 2014, it has again had massive mangrove mortality. AN site is located to the southwest of CGSM system, 200 m to the south of the Aguas Negras channel mouth ( Figure 1). This area, which was named AN according to the initials of the channel, lost 90% of its mangrove cover (BA = 0.1 m 2 ha −1 ; [35]). After the ProCiénaga project, the mangrove forest recovered and by 2014 had a BA of 27.3 m 2 ha −1 [41], with A. germinans, L. racemosa, and R. mangle trees. Finally, a site (AHU) located in the northern margin of the Ciénaga La Ahuyama (AHU) that experienced 100% tree mortality is indirectly influenced by the canals opened by the ProCiénaga project. No forest regeneration has been observed in AHU after the mangrove dieback [44,45].

Ecosystem C Stock
This study estimated the organic carbon (OC) stock in trees, living pneumatophores, living roots, and soil. The tree component is considered highly vulnerable, while the soil component has been described as the main carbon stock present in mangroves [46]. Direct field measurements of pneumatophores and belowground roots were also carried out because their contribution to total C stocks is usually estimated from allometric equations developed in sites with different conditions [47].

Carbon Stocks in Trees
In each study site, we delimited 20 × 30 m plots located 30 m from water bodies (coastal lagoons or channels). In RIN, where A. germinans predominates, we established two plots. In K22 and AN, there are stands where A. germinans or L. racemosa predominates, and therefore we established two plots in each of these stands. C in trees was not considered in AHU because it still does not have adult mangrove vegetation. In each plot, we counted and measured the stem diameter at breast height (DBH) of individuals with DBH equal to or greater than 2.5 cm during February of 2017. Based on these data, we estimated tree density and basal area (BA) according to Cintrón and Schaeffer-Novelli [48]. Aboveground biomass was estimated from the allometric equations developed for A. germinans and L. racemosa in Bahia de Cispatá and San Andres island, respectively [49,50]. The Bahia de Cispatá (Sinú River Delta) and San Andres Island are in the Colombian Caribbean. Carbon content in trees was estimated by multiplying the aboveground biomass by 0.48 as recommended by Kauffman and Donato [46].

Root Carbon Stocks
Root sampling was carried out during December 2016 and January 2017 in five 25 × 25 cm subplots located in the center and corners of one of the 20 × 30 m plots used to estimate the tree structures per site (RIN) or per species and site (K22 and AN). C in roots was not considered in AHU because it remains without adult mangrove vegetation. Inside subplots, pneumatophores were collected with adventitious roots attached to them ("feeding roots" according to Angeles et al. [51]). Belowground roots were collected close to each subplot using a PVC coring device [52][53][54], obtaining samples 10.2 cm in diameter and 1 m in length. Samples were stored in labeled plastic bags and transported on ice to the laboratory. Collected pneumatophores were separated into living and dead according to texture, appearance, and floatability. On the other hand, the core samples were washed through a 1 mm sieve to separate roots from sediments and other materials. Then, living and dead roots were differentiated by floatability [54], texture, and color. Living pneumatophores and belowground roots were packed separately in paper bags and then dried at 60 • to constant weight [46]. Soil coring is impaired in zones where large roots are present, and therefore roots above 20 mm were not considered in the study.
Ten and fifteen samples for each root structure (pneumatophores and belowground roots) were extracted in L. racemosa and A. germinans, respectively, for the percentage of C content estimation. Samples were dried at 60 • C up to a constant weight, homogenized with a tissue pulverizer (Retsch, GmBH & Co., Haan, Germany), and processed in an elemental analyzer (Eurovector EA 3000, Pavia, Italy). C stocks (Mg C ha −1 ) were estimated by multiplying belowground root biomass and pneumatophore biomass by the percentage of C content of each root structure.

Soil Carbon Stocks
Three soil cores were extracted in one 20 × 30 m plot per site (RIN and AHU) or per species and site (K22 and AN) up to a depth of 2 m with a 5 cm-diameter Russian peat borer during December 2016 (RIN, K22 and AN) and August 2017 (AHU). Cores were extracted in opposite corners and the center of plots. The cores were divided into 10 cm segments in the first meter, while from 1 to 2 m depth, they were divided into four segments (100-115, 115-130, 130-150, and 150-200 cm). From these segments, a subsample of 10 cm was taken in the middle point of the segment. Samples were stored in labeled plastic bags and transported on ice to the laboratory for organic matter, bulk density, and carbon content estimation.
In the laboratory, samples were dried at 60 • C up to a constant weight. Bulk density was estimated by dividing the sample dry weight by the sample volume according to Kauffman and Donato [46]. Organic matter (OM) content was estimated by the combustion of 1 g of soil sample at 500 • C for two hours [55]. Soil samples were previously pulverized using a FRITSCH Analisette 3 SPARTAN-Pulverisette 0, after the removal of wood debris, rocks, and shells. C content was estimated in each of the ten segments in the first meter of soil with an elemental analyzer (Eurovector EA 3000., Pavia, Italy) using samples of 10 to 20 mg. Prior to OC analyses, samples were acidified with HCL 1 M for carbonate removal. C content analyses were carried out at Leibniz Centre for Tropical Marine Research-ZMT in Bremen-Germany. C from 1 to 2 m depth was estimated by exploring significant relationship between soil properties variables and % of C. A detailed description of these analyses is shown in the statistical analysis section. The results of % C were converted into Mg C ha −1 as recommended by Kauffman and Donato [46], using bulk density and soil segment length. C stocks in the two first meters of soil were estimated by summing C content of all soil segments. Then, C stocks per site were estimated by averaging the C stocks values of different cores extracted (RIN = 3; K22 = 6; AN = 6; AHU = 3). Soil depth in each site was measured with a stainless still rod of 0.25" and 6 m length that was inserted into the soil four times within each plot to parent rock.
The mangrove C stocks per site were estimated according to Kauffman and Donato [46] by summing the value of each component (trees, pneumatophores, belowground roots, and soil) and by estimating the standard error of the sum.

Statistical Analyses
To estimate difference in the C stock among sites, we ran one-way ANOVAs. Data normality was verified with the Shapiro-Wilk test using ANOVA residuals, and the homoscedasticity was verified with the Bartlett test. When significant differences were found in the ANOVA test, we used the post hoc Tukey test to identify which treatment comparisons showed significant differences. The non-parametric Kruskal-Wallis test was used to compare C content between sites of belowground roots and pneumatophores owing to the violation of the homoscedasticity assumption. We used a significance level of α < 0.05. Analyzes were run with R-3.5.1 for Windows. In addition, we ran simple regression anal-ysis between the organic matter content (%) and bulk density as predictors of the C (%) obtained from 0 to 1 m soil samples.

Carbon Stocks in Trees
We found a mature forest in RIN, with low tree density and dominated by A. germinans with individuals up to 40 cm DBH (Table 1). This site had the highest BA (27.1 ± 0.5 m 2 ha −1 ) and tree biomass (212.6 ± 3.8 Mg ha −1 ). In the restored sites of K22 and AN, tree density was higher and dominated by A. germinans or L. racemosa, and DBH was lower in comparison to RIN (up to 35 cm of DBH). Tree biomass in K22 was 94.3 ± 12.1 Mg ha −1 , while in AN it was 153.2 ± 11.8 Mg ha −1 . As described in the study area section, the AHU site has not recovered its mangrove coverage since the massive forest mortality of the 1970s and 1980s. The carbon stocks in trees were significantly different between sites (F = 40.9; df = 3; p < 0.05) and were between 45.2 ± 5.8 (K22) and 102.1 ± 1.8 Mg C ha −1 (RIN). RIN had the largest tree biomass and therefore the largest aboveground C stock, followed by AN, K22, and AHU (Table 1).  [49] and L. racemosa: [50]). ** Estimated using a factor of 0.48 [46]. Plots of 20 × 30 m (RIN: n = 2; K22 and AN: n = 4). Different letters represent significant differences among sites.

Root Carbon Stocks
The percentage of C in A. germinans and L. racemosa pneumatophores was 45 ± 0.6% and 43 ± 0.8%, respectively. On the other hand, belowground roots had a lower C percentage: 37 ± 1.4% for A. germinans and 42 ± 0.1% for L. racemosa. Considering these values for the percentage of C and the root biomass structure, the C content of pneumatophores was significantly higher in RIN, the conserved site, compared to restored sites (F = 7.43; df = 2; p < 0.05) ( Table 2). In addition, belowground roots did not demonstrate significant differences among sites (χ 2 = 5.47; df = 2; p < 0.05).

Soil Carbon Stocks
In the first meter, the average soil bulk density (BD) was between 0.18 ± 0.01 (g cm −3 ) in AHU and 0.57 ± 0.08 (g cm −3 ) in RIN. As for organic matter (OM) content, the lowest average value occurred in RIN (15.3%) and the highest occurred in AHU (47.46%). A strong and positive relationship was observed between % OM and % OC in the samples collected at the study sites in the first meter of soil (r 2 = 0.96; p = 0.001; Figure 2). According to the linear model, 96.14% of OC variability can be explained by the OM content (%). On the other hand, an equally strong but negative exponential relationship was observed between BD and % OC concentration (Figure 3). Since the relationship between the OM and % OC values of the soil samples at one meter had a higher coefficient of determination than the BD vs. % OC relationship, the first relationship was used to estimate the OC concentration of the second meter of soil.

Soil Carbon Stocks
In the first meter, the average soil bulk density (BD) was between 0.18 ± 0.01 (g cm −3 ) in AHU and 0.57 ± 0.08 (g cm −3 ) in RIN. As for organic matter (OM) content, the lowest average value occurred in RIN (15.3%) and the highest occurred in AHU (47.46%). A strong and positive relationship was observed between % OM and % OC in the samples collected at the study sites in the first meter of soil (r 2 = 0.96; p = 0.001; Figure 2). According to the linear model, 96.14% of OC variability can be explained by the OM content (%). On the other hand, an equally strong but negative exponential relationship was observed between BD and % OC concentration (Figure 3). Since the relationship between the OM and % OC values of the soil samples at one meter had a higher coefficient of determination than the BD vs. % OC relationship, the first relationship was used to estimate the OC concentration of the second meter of soil.   The elemental analyzer results for the first meter showed that the OC ranged between 1.6% (at a 70 cm depth in RIN) and 26% (the first 10 cm in AHU). In addition, the inorganic carbon of the soil segments was lower than 1% in all the sites. On the other hand, the estimation of the C content in the second meter of soil had values between 0.2% (at a 150-200 cm depth in RIN) and 36% (at 115-130 cm in K22). The soil OC at 0-1 m depth was significantly different among all the sites, while at 0-2 m, it was only significantly higher in K22 in relation to other sites (F = 262.9; df = 3; p < 0.05) ( Table 3). The soil depth to bedrock ranged between 2.2 m in AN and 4.4 m in RIN. The elemental analyzer results for the first meter showed that the OC ranged between 1.6% (at a 70 cm depth in RIN) and 26% (the first 10 cm in AHU). In addition, the inorganic carbon of the soil segments was lower than 1% in all the sites. On the other hand, the estimation of the C content in the second meter of soil had values between 0.2% (at a 150-200 cm depth in RIN) and 36% (at 115-130 cm in K22). The soil OC at 0-1 m depth was significantly different among all the sites, while at 0-2 m, it was only significantly higher in K22 in relation to other sites (F = 262.9; df = 3; p < 0.05) ( Table 3). The soil depth to bedrock ranged between 2.2 m in AN and 4.4 m in RIN.

Mangrove Total Carbon
The total C stock considering trees, pneumatophores, belowground roots, and soil up to a 2 m depth ranged between 398.2 ± 76.7 and 1160.8 ± 130.9 Mg C ha −1 in AHU and K22, respectively. The aboveground C was significantly different among all the sites and ranged from 0 Mg C ha −1 in AHU, the deforested and unvegetated site, up to 102 Mg C ha −1 in RIN, the conserved site. On the other hand, the belowground C up to a 2 m depth was 2.7-3.8 times significantly higher in K22, a restored site, than the other sites ( Figure 4). The largest C stocks were those of soils (0-2 m depth) that represented more than 70% of total C stocks. C stocks in trees represented 25%, 4%, and 16% of total mangrove carbon stocks in RIN, K22, and AN, respectively, and were the second-highest component in importance.
Regarding the contribution of roots and pneumatophores, at all the sites, their combined contribution was less than 3%.

Indirect Estimations of Soil Carbon
The use of the relationship between OM, determined by the loss on ignition method

Indirect Estimations of Soil Carbon
The use of the relationship between OM, determined by the loss on ignition method (LOI), and the concentration of C (%) in soil samples from mangroves, marshes, and seagrasses, is recommended for determining the C concentration when it is not possible to access an elemental analyzer [1,46]. Kauffman and Donato [46] reported a positive association (r 2 = 0.59) between these two variables in mangrove soil samples of the Republic of Palau. In the present study, we also found a positive and stronger relationship (r 2 = 0.96, p = 0.00, Table 4), which can offer a cost-effective solution in future studies to predict the concentration of soil C based on the values of OM. When comparing the results between the models, the percentage of C obtained by the model of Kauffman and Donato [46] had significant differences compared with the values measured in the laboratory and with those obtained with the regression model found in this study (χ 2 = 44.907, df = 2, p < 0.05). A higher relationship in the model determined in this study may be due to the fact that the soil samples in this case were acidified to remove inorganic C prior to determination in the elemental analyzer.
Bhomia et al. [57] presented another relationship that could be used to estimate the soil C content, between the bulk density of the soils and the C percentage. These authors found a strong and inverse relationship between these variables in mangroves of the Caribbean and Pacific coasts of Honduras (r 2 = 0.83, p = 0.00). In CGSM soils, we observed a similar strong relationship between the two variables (r 2 = 0.80, Table 5). This relationship is useful in studies where there is no elemental analyzer or muffle furnace. There were no statistically significant differences between the values calculated with these models and the % C measured in the laboratory (χ 2 = 8.3894, df = 2, p < 0.05).

Mangrove Total Carbon
The carbon stocks found in the CGSM mangroves are within the range recorded for neotropical mangroves including mangroves on the Caribbean and Pacific coasts of Colombia (Table 6). However, the average (604.2 ± 180 Mg C ha −1 up to 2 m depth) is lower than the global average (738 ± 28 Mg C ha −1 up to 1 m depth) estimated by Alongi [58], but nevertheless the mangroves of the CGSM store significant C reserves, despite continued disturbances in the area.
Spatial variations in hydroperiod and physicochemical conditions drive soil biochemistry together with soil anoxia levels and mangrove regulators and resources. As a result, there is a broad variability in C stocks in different sites (Table 6). Tree structure is related to mangrove C stocks in some sites. For example, Adame et al. [65] found that tall forest had the highest C stocks in comparison to medium and dwarf forest in Mexico (987, 623, 381 Mg C ha −1 , respectively). However, other authors, such as Kauffman et al. [66], Bhomia et al. [57], and Kauffman et al. [64], did not find a correlation between forest height or structure and total C stocks in mangroves of the Dominican Republic, Honduras, and Mexico, respectively. Our results agree with the findings of the latter authors, as the highest C stocks were found in K22 (deforested/restored), the site that had the lowest structural development among the three vegetated sites (Table 1). C stored in tree biomass of the CGSM sites (45-102 Mg C ha −1 ) is comparable to the values recorded in other mangroves of the Colombian Caribbean, such as Cispatá (65 ± 10 Mg C ha −1 ) and Golfo de Urabá (38-83 Mg C ha −1 ) [49,67]. The C stocks in trees were slightly low compared to the global average of 109.3 Mg C ha −1 [58] and were as expected: higher in the conserved site (RIN = 102 Mg C ha −1 ) than the sites where mangrove vegetation reestablished after the massive mangrove loss and subsequent hydrologic restoration (AN = 73.5 Mg C ha −1 and K22 = 45 Mg C ha −1 ). 1 Aboveground biomass + belowground biomass + soil up to 2 m depth. 2,4 Aboveground biomass + below ground biomass + necromass+ soil up to 0.45 m depth. 3 Aboveground biomass + belowground biomass + downed wood + litterfall + soil up to 1 m depth. 5 Aboveground biomass + belowground biomass + soil up to 0.3 m depth. 6 Aboveground biomass + downed wood + soil up to 2.7 m depth. 7 Aboveground biomass + belowground biomass + downed wood + dead trees + soil up to 1.5 m depth. 8 Aboveground biomass + belowground biomass + downed wood + dead trees + soil up to 3 m depth. 9 Aboveground biomass + belowground biomass + downed wood + soil up to depth of parent material. 10 Aboveground biomass + belowground biomass + downed wood + soil up to 2.5 m depth. * Graph estimated values.
There was a strong increase in K22 soil C stocks when the second meter of soil depth is included. This is explained by a recurrent delivery of sediments that favor large C preservation from 1 to 2 m depths. Since the ProCiénaga restoration project in 1996, the opened channels have been dredged sporadically, and the sediments extracted were deposited along the channel's borders. Monitoring of CGSM mangroves carried out between 2014 and 2015 indicated that this approach negatively impacted the restoration process, especially around the middle of the Clarín channel's course, where K22 was located [68]. This was due to two impacts. First, sediment disposal from channel dredging created spoil banks that functioned as dikes and limited water flux from the channels to the mangrove forest. Second, channel currents and rains during wet seasons relocated sediments from the channel borders to the forest interior or back to the channel bed. Dredged material continued to be deposited along the Clarín channel borders in the hydrodynamic recovery campaign carried out between 2016 and 2019 [42]. Soil C stock up to 2 m depth in AN and AHU (deforested/restored and deforested/unvegetated sites), compared to RIN (conserved), did not have significant differences, which shows that C loss in disturbed sites can be very slow and can take years or decades [69]. In our case, the massive forest mortality already happened 40 to 50 years ago.
The localization of the sampling sites in relation to the Magdalena River could help to explain the current situation regarding carbon concentrations in soils and vegetation. K22, being located further from the Magdalena River discharges than AN (Figure 1), receives lower nutrient inputs. Before the hydraulic disturbances, the carbon vegetation:soil ratio was probably lower in K22 than in AN. Additionally, the formation of dikes decreased the flow of the Clarín channel through the mangroves, affecting the tree's biomass in K22. This could explain why this site has the lowest aerial biomass compared to AN, the other deforested/restored site. This situation was aggravated during 2015 and 2016, when, due to high salinities (greater than 90 in the interstitial water at 50 cm depth), a new mangrove mortality occurred at the K22 site [21]. The hyper-salinization of the soil accelerated because of the El Niño event, which in 2015 was very intense (average ONI index 1.5), and the water flow limitation by debris banks.
The soil depth to the base rock in all the study sites was greater than the 2 m depth limit used in this study to estimate C stocks directly. Thus, C stocks of study sites are greater than those shown in Figure 4 and Table 6. If C stocks found at 0-2 m are extrapolated to the total soil depth, carbon stocks would increase to 749.2 Mg C ha −1 ; 1738.6 Mg C ha −1 ; 489.8 Mg C ha −1 and 796.5 Mg C ha −1 in RIN, K22, AN, and AHU, respectively. C stocks extrapolated to the total soil depth have two important error sources: The C content variability along the soil profile and the channel maintain influence in some sites. However, the estimation of C up to base rock highlights the considerable C stocks present in CGSM, even in non-recovered sites.
Our results showing that the main C stocks of mangroves are soils (>71%) (Figure 4) agree with other studies of neotropical mangroves [57,64,65]. In contrast, pneumatophores and belowground roots had low contributions to total C stocks (<3%), a low value compared to reports in Mexico for the contribution of these structures to C stocks (1-16%) [63]. Technique differences may explain this difference. In this study, soil cores were used to estimate the belowground root biomass, while Adame et al. [63] used the general equation of Komiyama et al. [70]. While the soil core technique does not allow the sampling of coarse roots (>20 mm), because these would limit the penetration of the core into the soil, the general equation was developed from the biomass found in a trench 15.5 m long by 20 cm wide and 1 m deep, where all root sizes were sampled.
Given the sensitivity of the vegetation of the CGSM lagoon system to hydrological variability [29], conditions leading to forest dieback will have consequences on mangrove C stocks. A useful tool to raise awareness of the need to preserve or restore an ecosystem is to carry out ecosystem service valuation in sites with different degrees of degradation including preserved sites and sites where ecosystems have been lost or damaged. Although the suggested carbon value to achieve the Paris Accord goals (limit temperature rise to below 2 • C) is USD 40-80/t CO 2 -equivalents for 2020, an analysis by the World Bank indicates that 51% of the emissions covered by C pricing initiatives have been below USD 10/t CO 2 -equivalents [71]. In Colombia, C emissions had a government-regulated value of USD 5.08/t CO 2 -equivalents in 2017, a year of C biomass measurements in CGSM. Based on this valuation, the C in the CGSM sites, expressed as potential CO 2 emissions or CO 2 -equivalents, would have an economic value between USD 7283 and 21,230/ha −1 (AHU and K22, respectively) ( Table 7). Our estimations of the value of C in CGSM sites are higher than those reported by De la Peña et al. [72]. These authors only considered carbon stocks in trees, while we considered above and belowground C stocks. In addition, these authors used a different reference value of t CO 2 -equivalents, while we used the government-regulated value in Colombia.

Conclusions
This study showed that, although the massive mortality of the mangrove forests of CGSM occurred several decades ago, these mangrove forests have important reserves of C both in the aboveground and belowground biomass and soils. In addition, these reservoirs are constantly influenced by channel dredging and respond to these actions in a contrasting way: the ProCienaga channels open greater hydrological connectivity in some sites that favor the recovery of mangroves (AN) or have not managed to reduce the salinity of the soil and recover the mangrove vegetation in places far from its influence (AHU). At other sites (K22), inadequate canal maintenance operations have affected tree biomass and increased underground carbon stocks by reducing hydraulic connectivity and burying mangrove soils with dredging materials that are relocated from the edges of the channels towards the interior of the forest due to heavy rains and flood events. Therefore, the future management of the CGSM should increase the water connectivity of the sites away from the channels and improve the management of the dredged sediments to consolidate the mangrove recovery trend. Finally, together with data on the economic valuation of C stocks in forests, these results will be useful for decision makers when considering management actions in ecosystems where rapid environmental changes are taking place in mangrove areas.