Effects of Invasive Spartina alterniflora Loisel . and Subsequent Ecological Replacement by Sonneratia apetala Buch .-Ham . on Soil Organic Carbon Fractions and Stock

Background and Objectives: The rapid spread of invasive Spartina alterniflora Loisel. in the mangrove ecosystems of China was reduced using Sonneratia apetala Buch.-Ham. as an ecological replacement. Here, we studied the effects of invasion and ecological replacement using S. apetala on soil organic carbon fractions and stock on Qi’ao Island. Materials and Methods: Seven sites, including unvegetated mudflat and S. alterniflora, rehabilitated mangroves with different ages (one, six, and 10 years) and mature native Kandelia obovata Sheue, Liu, and Yong areas were selected in this study. Samples in the top 50 cm of soil were collected and then different fractions of organic carbon, including the total organic carbon (TOC), particulate organic carbon (POC), soil water dissolved carbon (DOC) and microbial biomass carbon (MBC), and the total carbon stock were measured and calculated. Results: The growth of S. alterniflora and mangroves significantly increased the soil TOC, POC, and MBC levels when compared to the mudflat. S. alterniflora had the highest soil DOC contents at 0–10 cm and 20–30 cm and the one-year restored mangroves had the highest MBC content. S. alterniflora and mangroves both had higher soil total carbon pools than the mudflat. Conclusions: The invasive S. alterniflora and young S. apetala forests had significantly lower soil TOC and POC contents and total organic carbon than the mature K. obovata on Qi’ao Island. These results indicate that ecological replacement methods can enhance long term carbon storage in Spartina-invaded ecosystems and native mangrove species are recommended.


Introduction
Coastal wetland ecosystems are efficient "blue carbon" (carbon stored in coastal and marine ecosystems) reservoirs [1] and they play important roles in reducing man-made CO 2 emissions globally [2,3].Due to their high productivity and soil carbon storage, mangroves are among the most carbon-rich ecosystems [4][5][6].Mangroves only cover 0.5% of the coastal area and less than 0.1% of the total global continental surface area [4].However, they contribute 10%-15% to the total coastal Forests 2019, 10, 171 2 of 13 sediment carbon and 10%-11% of the carbon exported to the ocean [4][5][6].More than half of the total ecosystem carbon (49%-90%) in mangrove forests is stored in the soil [5,7].The soil organic carbon is mainly derived from autocthonous carbon sources [8], such as tree litter and root decomposition [4,6].Therefore, a vegetation composition shift would greatly affect the quality and quantity of plant material input to the soil and change the soil organic carbon pool and dynamics [9][10][11][12][13].
Spartina alterniflora Loisel., which is native to the eastern coast of North America, has invaded most of the mudflat areas that surround the mangrove forests of southeastern China [10,14,15].A shift of the vegetation community related to plant invasion or restoration can dramatically change soil organic matter content and its dynamics [9,16,17].Short term invasion of S. alterniflora increased the soil carbon content by 0.37-7.43times relative to the mudflat and native saltmarsh communities [13] in the Yangtz River Estuary.The highest carbon content in a S. alterniflora community was observed at 40-60 cm depth, owing to the below ground biomass in Jiuduansha Island [9].However, only the top 10 cm of soil, the organic carbon content of the mudflat was significantly increased due to invasive Spartina in mangrove ecosystems in Zhangjiang Estuary [10].The effect of invasive species on native ecosystems can be influenced by vegetation biomass, litter production, and quality and microbiological activities [12,13,18,19].Therefore, it is useful to conduct more field studies on the invasive effects of S. alterniflora on soil organic carbon dynamics in different coastal wetland ecosystems [17].
The restoration of native or natural habitat by removing invasive plants is an important method for ecosystem management [20][21][22].Owing to adverse influences on the native ecosystems [15,23], different measures have been used to control the spread of S. alterniflora, including ecological replacement using native or exotic mangrove species [21,24,25] after physical treatment.As a fastgrowing mangrove species, the exotic Sonneratia apetala Buch.-Ham.has been successfully used to eradicate S. alterniflora and to increase the mangrove areas in China [24,25].The primary goal of restoration is to reconstruct ecosystems with similar structure and function to those that were previously degraded or lost [22].The mangrove vegetation community can be re-established [24,25] and the associated macrofaunal community and trophic relationships can also be gradually restored [10,20,21] after the S. alterniflora is replaced with mangrove species.Short-term restoration of Kandelia obovata Sheue, Liu, and Yong, a mangrove species that is native to southeastern China, did not significantly enhance the soil organic carbon (SOC) stock after 15 years [21].However, the soil carbon accumulation and stock in mangrove ecosystems can be affected by factors, such as plant species, stand age, environmental conditions, and organic matter retention and export [5,8,14,26].S. apetala has shown great potential in biomass accumulation and litter production when compared to the native mangroves [27,28].Thus, it is necessary to determine whether the soil carbon stock capacity can be enhanced rapidly when the invasive S. alterniflora is replaced by exotic S. apetala to assess the restoration effectiveness.
The alteration of vegetation resulting from invasion or restoration can change the source of organic material entering the soil, influence the lability and preservation of soil organic matter, and affect the carbon transformation and biogeochemical cycle [17,29].The proportions of labile and recalcitrant carbon are important in the turnover and decomposition rate of total soil organic carbon and can influence long-term carbon sequestration and storage [16,17].Soil particulate organic carbon (POC), water-soluble carbon (DOC), and microbial biomass carbon (MBC) are important fractions of soil labile organic matter [17,30,31].DOC is the direct reservoir of nutrients that are available for soil microbial metabolism [17,30].The microbial biomass in the soil regulates the transformation of soil organic matter [17].Despite its significance, few studies have focused on the invasive effects of S. alterniflora on the dynamics of soil labile carbon.Inconsistent results have been reported between saltmarsh and mangrove ecosystems [13,17,32].It is necessary to conduct additional studies on the effects of invasion and restoration on soil labile organic carbon pools to determine the process that is responsible for carbon sequestration and preservation in coastal wetland ecosystems.
The Qi'ao Island mangrove ecosystem, which is located in southeastern China, was invaded by S. alterniflora in the 1990s and eradication efforts using exotic mangrove S. apetala have been conducted Forests 2019, 10, 171 3 of 13 since 1999 [24].Short-term restoration of native mangroves after the removal of S. alterniflora did not enhance the SOC stock as compared to mudflats [10].The rapidly growing S. apetala is superior in its growth characteristics to native mangrove species [33,34], and different results may occur when the S. alterniflora is replaced by S. apetala.In this study, we compared the content of different fractions of organic carbon and the total carbon stock of the top 50 cm soil to determine (1) whether the invasion of S. alterniflora increases soil organic carbon content; (2) whether S. apetala plantings can increase the soil carbon content; and, (3) which fractions of soil labile organic carbon are most impacted by the invasion of S. alterniflora and subsequent replacement by S. apetala.

Study Area
The Qi'ao Island Mangrove Nature Reserve is in Zhuhai City, Guangdong Province, southeastern China (22 • 26 N, 113 • 38 E, Figure 1).The area of the reserve is 5103.77ha and the area of mangroves is approximately 700 ha.It has a subtropical marine monsoon climate with an average temperature of 22.4 • C. Qi'ao Island has an irregular semidiurnal tide with a mean high tide level of 0.17 m and a mean low tide level of −0.14 m.The annual mean salinity of seawater is 18.4% [24,25].effects of invasion and restoration on soil labile organic carbon pools to determine the process that is responsible for carbon sequestration and preservation in coastal wetland ecosystems.
The Qi'ao Island mangrove ecosystem, which is located in southeastern China, was invaded by S. alterniflora in the 1990s and eradication efforts using exotic mangrove S. apetala have been conducted since 1999 [24].Short-term restoration of native mangroves after the removal of S. alterniflora did not enhance the SOC stock as compared to mudflats [10].The rapidly growing S. apetala is superior in its growth characteristics to native mangrove species [33,34], and different results may occur when the S. alterniflora is replaced by S. apetala.In this study, we compared the content of different fractions of organic carbon and the total carbon stock of the top 50 cm soil to determine (1) whether the invasion of S. alterniflora increases soil organic carbon content; (2) whether S. apetala plantings can increase the soil carbon content; and, (3) which fractions of soil labile organic carbon are most impacted by the invasion of S. alterniflora and subsequent replacement by S. apetala.

Study Area
The Qi'ao Island Mangrove Nature Reserve is in Zhuhai City, Guangdong Province, southeastern China (22°26′ N, 113°38′ E, Figures 1).The area of the reserve is 5103.77ha and the area of mangroves is approximately 700 ha.It has a subtropical marine monsoon climate with an average temperature of 22.4 °C.Qi'ao Island has an irregular semidiurnal tide with a mean high tide level of 0.17 m and a mean low tide level of −0.14 m.The annual mean salinity of seawater is 18.4% [24,25].S. alterniflora invaded the gaps of the natural mangroves and mudflats of the Qi'ao Island reserve in the early 1990s.It spread and formed monocultures rapidly and it covered 227 ha by 1995.Due to the construction of fish or shrimp ponds, some of the S. alterniflora was removed, but it remained dominant in the mudflats, with an area of 59 ha in 1999 [24].To control the expansion of S. alterniflora, S. apetala was introduced to Qi'ao Island in an attempt to replace the S. alterniflora stands after its aboveground biomass was cleared (without obvious disturbance on deep soil ), which started in 1999 [24].The S. apetala formed dense forests with a canopy about 60%-70% within two years and nearly and efficiently inhibited the growth of S. alterniflora.Only rare S. alterniflora individuals were observed four years after the eradication project [24].By 2011, the mangrove area of Qi'ao Island increased to more than 300 ha and the exotic S. apetala occupied about 95% of the restored mangroves.However, there are still 15 native true mangrove species that belong to 10 families and 13 genera on this island.Therefore, an area of natural mangroves that is dominated by native K. obovata was found as a reference site in this study.
A sampling design using space as a substitution for time was developed to study the soil carbon fractions and pools that are associated with the invasion and eradication of S. alterniflora and the subsequent establishment of mangroves.Seven sampling sites with different vegetation communities were chosen according to a chronosequence: (1) pre-invaded mudflat without any vegetation (abbreviated as MF); (2) invasive S. alterniflora stands (abbreviated as Sal); (3) exotic mangrove S. apetala planted in 2010, mixing with S. alterniflora (abbreviated as SA1); (4) mangrove forests that are exclusively composed of S. apetala planted after the aboveground biomass of S. alterniflora was removed manually in 2005 (abbreviated as SA6) and in 2002 (abbreviated as SA10); and, (5) an area of natural mangroves dominated by native K. obovata was selected as reference (abbreviated as KO) (Figure 1).

Sample Collection and Measurement
Soil core samples to 50 cm depth were collected from each site using a 5-cm diameter PVC pipe during August, 2011.Three 1 m × 1 m plots, separated from each other by a distance of more than 100 m, were established within each community.Within each plot, two or three soil cores were collected and divided into four sections (0-10 cm, 10-20 cm, 20-30 cm, and 30-50 cm).Samples from the same depth were pooled together as a replicate.Areas with flat soil surface without obvious root and crab burrows were selected to avoid the effects of biomass input from plant root material and animals.After obvious plant material and shells were removed from the soil, the redox potential (Eh), pH, and temperature were measured in situ using a Orion 5-star multiple function electrometer (Thermo Fisher Scientific, San Jose, CA, USA).All of the samples were divided into two subsamples.One subsample was placed into a 100 cm 2 aluminum cutting ring and the other was kept in a plastic bag.All of the samples were transported to the laboratory and kept at 4 • C until analysis.
The samples in the aluminum rings were freeze-dried using a vacuum freezer (Boyikang, Beijing, China) and then weighed to calculate bulk density.The dried samples were ground into a fine powder using a mortar and pestle to pass through a 250-µm sieve.Soil total organic carbon (TOC) was measured using the K 2 Cr 2 O 7 -H 2 SO 4 wet oxidation method.In brief, about 0.1-0.5gdried samples were added into the test tube, and 0.1g of Ag 2 SO 4 and 5.00 ml K 2 Cr 2 O 7 -H 2 SO 4 standard solution (0.80 mol/L) were added.The tubes were heated under oil-bathing (185-190 • C) and then kept boiling for 5 min.Two samples of quartz sand were used as blanks.All solution and residuals were transferred to a 250 ml beaker and kept at 60-70 ml.Subsequently, 5 ml H 3 PO 4 solution and 2-3 drops of o-phenanthroline indicator solution were added and titrated with FeSO 4 standard solution.
The POC contents of water stable aggregates were measured by the wet-sieving method using dried samples [35].In brief, about 10 g soil samples were submerged in water on a 53-µm screen that was contained in a shallow pan for 5 min.Afterwards, the sieve was moved vertically 50 times in 2 min in water.Materials retained on the sieve were transferred to a beaker and then dried in an oven at 60 • C. Subsequently, the carbon content of the dried materials (POC) was measured using K 2 Cr 2 O 7 -H 2 SO 4 wet oxidation method.
DOC was extracted from 10 g of fresh soil sample by adding 20 mL distilled water to 50-mL centrifuge tubes.The mixture was shaken for 30 min under 250 rpm and then centrifuged for 10 min at 6000 rpm [17].The supernatant liquids were filtered through 0.45 µm filters and the an automated TOC Analyser measured the carbon contents (Shimadzu TOC-VCPH, Kyoto, Japan).The MBC contents were determined by the fumigation-extraction method [36].In brief, subsamples of 10 g moist soil were fumigated with ethanol-free chloroform for 24 h under dark conditions, then quickly transferred to centrifuge tube, extracted with 40 mL 0.5 M K 2 SO 4 , shaken for 30 min under 250 r/min, and centrifuged for 5 min under 5000 r/min, and then the supernatant was vacuum filtered through a 0.45 µm filter.The unfumigated soil was also extracted using the same method.The carbon contents of the extracts were measured by an automated TOC Analyser (Shimadzu TOC-VCPH, Kyoto, Japan).

Statistical Analysis
One-way ANOVA was performed to test the significance of differences of all soil properties between different communities in each soil depth interval.Data were tested for normality and, when required, log-transformed to meet this assumption for ANOVA analysis.All of the data analyses were performed using SPSS 19.0 (IBM SPSS Statistics for Windows, IBM Corp., Armonk, NY, USA).A redundancy analysis (RDA) was used to explore the relationship between the soil carbon contents and environmental factors and ages of vegetation.The RDA was performed using the Vegan package [37] in R software (R 3.5.1).

Effects of Invasion and Restoration on Soil Physicochemical Parameters
The soil temperature decreased with depth at all sites and MF had the highest temperature in the top 30 cm (p < 0.05, Table 1).The soil water contents ranged from 0.4 to 0.57 and then decreased with depth in most of the sites except for Sal and SA6.KO had the highest water contents in the top 20 cm of soil (p < 0.05).For bulk density, negative relationships were observed between soil bulk density and water content (p < 0.05).The lowest bulk density in the surface 10 cm was found at the sites of KO and SA10 (p < 0.05), showing the effects of mangrove restoration.The pH showed no changes with depth in all sites except for MF.At each depth, MF had a higher pH than the other sites, except for 20 cm, where no significant differences were found.MF had the lowest Eh values and S. alterniflora invasion significantly increased the Eh values; Sal and SA1 sites showed the largest Eh values in the surface soil (0-10 cm).The restoration of mangroves following the eradication of S. alterniflora reduced the Eh values at all depths, except for 40-50 cm at SA6 and 20-40 cm at KO.

Effects of Invasion and Restoration on Soil Carbon Fractions and Carbon Stock
When compared with MF, the invasion of S. alterniflora and mangrove restoration significantly increased the soil TOC contents at all depths (Figure 2, p < 0.05).In the surface 0-10 cm, KO had much larger TOC contents than the other vegetation types (p < 0.05).Significantly higher TOC contents at SA6 were observed than those at Sal, while no significant differences were observed for other sites.At 10-20 cm depth, KO showed the highest TOC contents and SA10 also had higher TOC contents than the other sites.Although significant differences were found for soil TOC contents among all of the vegetation communities in the 20-50 cm depth, the discrepancies did not show any consistent trend.The soil POC contents had similar patterns at all sites and at each depth as those for soil TOC contents.In the 0-10 cm layer, Sal showed the largest soil DOC contents among all sites (p < 0.05).SA10 had more substantial soil DOC contents than MF (p < 0.05), while no significant differences were observed among the other sites.At 10-20 cm, MF had the lowest DOC contents and KO showed larger soil DOC levels than all of the other vegetation communities (p < 0.05), except for Sal.The soil DOC contents for Sal and SA1 at 20-30 cm were much larger than the other sites.The growth of vegetation significantly increased the MBC contents when compared with the bare mudflat (p < 0.05).SA1 had the largest soil MBC contents than the other sites at each depth (p < 0.05).The average total carbon pool for KO (0-50 cm) was 122.68 g C/m 2 (Figure 3), which was significantly higher (p < 0.01) than all of the other sites (48.08, 75.14, 76.37, 83.16, and 77.82 g C/m 2 for MF, Sal, SA1, SA6, and SA10, respectively).The bare MF had the lowest carbon pool when compared with all of the other sites (p < 0.01), while no significant differences were found among Sal and restored mangrove communities.The average total carbon pool for KO (0-50 cm) was 122.68 g C/m 2 (Figure 3), which was significantly higher (p < 0.01) than all of the other sites (48.08, 75.14, 76.37, 83.16, and 77.82 g C/m 2 for MF, Sal, SA1, SA6, and SA10, respectively).The bare MF had the lowest carbon pool when compared with all of the other sites (p < 0.01), while no significant differences were found among Sal and restored mangrove communities.The average total carbon pool for KO (0-50 cm) was 122.68 g C/m 2 (Figure 3), which was significantly higher (p < 0.01) than all of the other sites (48.08, 75.14, 76.37, 83.16, and 77.82 g C/m 2 for MF, Sal, SA1, SA6, and SA10, respectively).The bare MF had the lowest carbon pool when compared with all of the other sites (p < 0.01), while no significant differences were found among Sal and restored mangrove communities.

Relationship between Carbon Contents and Soil Physicochemical and Plant Community Age
The RDA results showed that the soils at 0-10 cm depth in different sites were well separated.MF, Sal, and young restored mangroves were similar, while KO was distinct from the other sites (Figure 4).Axis 1 and axis 2 accounted for 55.33% and 29.13% of the total variation, respectively.The first axis was positively correlated with mangrove age, while it was negatively correlated with bulk density and age of S. alterniflora.The second axis showed a positive relationship with the age of S. alterniflora and Eh values.
standard deviation of total carbon stock, letters indicate significant differences at p < 0.05).

Relationship between carbon contents and soil physicochemical and plant community age
The RDA results showed that the soils at 0-10 cm depth in different sites were well separated.MF, Sal, and young restored mangroves were similar, while KO was distinct from the other sites (Figure 4).Axis 1 and axis 2 accounted for 55.33% and 29.13% of the total variation, respectively.The first axis was positively correlated with mangrove age, while it was negatively correlated with bulk density and age of S. alterniflora.The second axis showed a positive relationship with the age of S. alterniflora and Eh values.

Effects of S. alterniflora and Restored Mangroves on Soil Organic Carbon Contents
The soil TOC contents under S. alterniflora ranged from 16.41 ± 1.21 to 25.42 ± 1.58 g/kg in this study, which was similar to the soil SOC contents of invasive S. alterniflora in the Zhangjiang mangrove ecosystem [10], but greater than the saltmarsh ecosystem in Jiangsu [17,31,38,39].This difference might relate to the invasion history [39] and geographic variation of the S. alterniflora biomass [14].When compared to the unvegetated mudflat, the TOC content at 0-20 cm was increased by 50%-107% in all layers of soil in the 20 years following the invasion of S. alterniflora.This was inconsistent with previous studies, demonstrating that S. alterniflora could enhance the SOC content by 1.69-5.3fold in the top 30 cm of soil [17] or by 73% in the 0-50 cm layers of soil [39] as compared with bare mudflat.The highest soil TOC content was found at 20-30 cm in S. alterniflora in this study, which was probably the result of the roots being mainly distributed in the top 30 cm [9].However, S. alterniflora only increased the TOC content within the surface 0-10 cm of soil in Zhangjiang mangrove reserve [10].Therefore, the effects of S. alterniflora on soil carbon accumulation related to above-and belowground biomass allocation and invasive duration, together with abiotic factors, such as tidal flooding, might determine the SOC distribution related to depth [39,40].

Effects of S. alterniflora and Restored Mangroves on Soil Organic Carbon Contents
The soil TOC contents under S. alterniflora ranged from 16.41 ± 1.21 to 25.42 ± 1.58 g/kg in this study, which was similar to the soil SOC contents of invasive S. alterniflora in the Zhangjiang mangrove ecosystem [10], but greater than the saltmarsh ecosystem in Jiangsu [17,31,38,39].This difference might relate to the invasion history [39] and geographic variation of the S. alterniflora biomass [14].When compared to the unvegetated mudflat, the TOC content at 0-20 cm was increased by 50%-107% in all layers of soil in the 20 years following the invasion of S. alterniflora.This was inconsistent with previous studies, demonstrating that S. alterniflora could enhance the SOC content by 1.69-5.3fold in the top 30 cm of soil [17] or by 73% in the 0-50 cm layers of soil [39] as compared with bare mudflat.The highest soil TOC content was found at 20-30 cm in S. alterniflora in this study, which was probably the result of the roots being mainly distributed in the top 30 cm [9].However, S. alterniflora only increased the TOC content within the surface 0-10 cm of soil in Zhangjiang mangrove reserve [10].Therefore, the effects of S. alterniflora on soil carbon accumulation related to above-and belowground biomass allocation and invasive duration, together with abiotic factors, such as tidal flooding, might determine the SOC distribution related to depth [39,40].
Mangrove structures are time and space dependent and carbon sequestration in mangrove soil can be determined by factors, such as climate and environmental conditions, vegetation community, and age of restored vegetation [8,26].The age of mangrove is important when quantifying the multi-decadal carbon balance, owing to the fact that the carbon accumulation is more influenced by stand age than other factors [7,26,41].Our RDA results also showed that the soil POC and TOC positively correlated with mangrove age (Figure 4).In the Shenzhen mangrove reserve, the younger mangroves, including exotic and native species (less than 20 years old), both showed significantly lower soil carbon content than the 70-year-old mature K. obovata ecosystem [7].Therefore, the time series of mangrove ages, especially the young mangrove forests with low carbon concentrations, should be considered to accurately assess the carbon accumulation ability of the mangrove ecosystem.
The soil TOC content of 10 year old S. apetala in this study was similar to the 2.2% global median organic carbon content that was estimated for mangrove soil [6].However, we also found that the mature K. obovata forest had 1.99-3.76times larger soil SOC content than the S. alterniflora community and restored S. apetala at a depth of 20 cm.The fast growing S. apetala can successfully inhibit the growth of S. alterniflora after physical removal [24,25], but it was suggested to have no superiority in soil carbon accumulation when compared with restored native mangrove species [7,27].In combination with the potential risk of invasion, the planting of exotic S. apetala is not recommended in the future.

Effects of Invasive S. alterniflora and Restored Mangroves on Soil Organic Carbon Fractions
The physical structure of soil can determine the organic matter sequestration and turnover under different vegetation [42][43][44].Soil aggregate structures are closely related to the amount of plant residuals and carbon contents [30] and they are influenced by biotic and environmental factors, such as vegetation type and land use changes [31,32,42].The RDA results in our study also suggested that the soil carbon contents were significantly influenced by the invasion/eradication of S. alterniflora and the mangrove restoration (Figure 4).Coarse fractions of soil organic matter (>53 µm), including macro-and micro-aggregates, are generally dominant among all of the fractions of soil.They explain most of the differences in organic carbon in soil among vegetation types in natural ecosystems [30].The plant residuals are first incorporated into the coarse fraction and are then transferred into the fine fraction [31,45].The fine fraction, including silt and clay aggregates, is prominent in long-term carbon sequestration due to its high proportion of old carbon, which had low accessibility to the microbiomass [11,31,43].
Soil POC contents accounted for more than 70% of soil TOC contents in this study and displayed similar variation with soil TOC content.Both of them were mainly derived from vegetation related sources.The relatively high proportion of POC showed that carbon derived from both mangroves and S. alterniflora is not easily transferred from the coarse fraction to the fine fraction due to the anaerobic conditions in coastal wetlands, and also perhaps to the refractory quality of the coarse fraction of organic carbon [17,31,44].Previous research also showed that the new carbon that was derived from invasive S. alterniflora was mainly distributed in macro-and micro-aggregates [16,46].We found that the mature mangrove forest reflected higher POC contents than the invasive S. alterniflora and short-term restored mangroves.Similar results were reported in that the coarse fractions in soil organic matter showed a larger rate of increase than the fine fraction in grassland invaded by woody plants [42].However, further separation of soil organic matter, especially the more stable and fertile humic substances [47], and long-term studies, are needed to determine the potential of soil aggregates for carbon sequestration during the chronosequence of the invasive S. alterniflora and restored mangroves [48].
Both the S. alterniflora invasion and the subsequent mangrove restoration increased the soil DOC and MBC contents relative to the mudflat.Positive relationships were observed among soil DOC, MBC, and TOC contents, indicating that the input of plant materials and the soil organic matter are the main sources for soil DOC and MBC [17,42].S. alterniflora had the highest soil DOC contents at 0-10 cm and 20-30 cm, which were mainly derived from the aboveground biomass and root material, respectively.These results showed that S. alterniflora could provide more available carbon for microbial activity and then potentially accelerate the organic carbon turnover rate.The highest MBC contents were found at the SA1 site, which might have resulted from the synergistic effects of invasion and mangrove establishment that supplied additional sources for the microbial community.The enhancement of invasive S. alterniflora on the MBC and DOC can cause the rapid decomposition of soil organic carbon, which is adverse to the soil carbon stability and storage in invaded coastal ecosystems [17,49].

Effects of S. alterniflora and Restored Mangroves on Soil Organic Carbon Stock
The invasion of S. alterniflora significantly increased the carbon stock at 0-50 cm soil depth as compared to the mudflat in this study.S. alterniflora significantly improved the total soil carbon density over mudflats and native marshes in 0-100 cm depths and the influences were more prominent with time following invasion in Jiangsu, China [39].However, spatial and temporal can affect the invasive consequence of S. alterniflora on carbon stock.For instance, the invasion of S. alterniflora did not yet alter the soil organic carbon pools at 0-60 cm depth in the Zhangjiang Estuary, China [10].
Mangroves generally have greater carbon storage and sequestration capability than S. alterniflora [40,50], but we found no significant differences in the carbon pool between S. alterniflora and mangrove forests, except for the mature K. obovata.Similar results were reported in Australian salt marshes encroached by mangroves: belowground carbon was not enhanced in 0-14-year mangroves, while the total carbon storage increased when the mangrove ages exceeded 30 years [51].The carbon stock of the mangroves in Qi'ao Island fell within the range that was previously reported for soil carbon stocks of mangrove ecosystems worldwide, ranging from 72 Mg C/ha in Saudi Arabia [52], to 126.14-657.61Mg C/ha in South China [53], to the highest 2063.6 Mg C/ha in Sulawesi [54].Although the soil carbon stocks were related to various environmental factors, such as degree of carbon inputs and export, biotic controls, and human pressures, the age of the ecosystem might be more significant in modulating mangrove carbon stocks [51].
In China, the fast growing exotic mangrove S. apetala has shown no advantages for carbon storage potential when compared to the native mangrove forests [7,27,54].In contrast, the native mangrove K. obovata not only controlled the growth of S. alterniflora, but it also had a larger soil carbon pool than the invasive S. alterniflora and gradually restored the food web complexity after 14 years of restoration [10,20,21].When considering the invasive potential and less prominent carbon stock ability of exotic S. apetala, native mangrove species are more appropriate in the future management in the control of S. alterniflora.
The carbon and nitrogen storage in coastal ecosystems may decrease if the mangrove vegetation is continually replaced by S. alterniflora [40], especially in areas that are subject to frequent anthropogenic disturbance [15].Furthermore, the spread of S. alterniflora can increase the emissions of greenhouse gases, which is detrimental to the "blue carbon" function of coastal wetlands [2,3,55,56].Therefore, a biological replacement method using mangrove species, especially native mangroves, to control S. alterniflora is recommended for the management of coastal wetlands [21,24,25].

Conclusions
The growth of S. alterniflora and mangroves can significantly increase the soil TOC, POC, and MBC content at 0-50 cm depths when compared to mudflats.S. alterniflora showed the highest soil DOC contents at 0-10 cm and 20-30 cm, indicating the impact of leaf litter and the developed root systems of S. alterniflora.The mixing of S. alterniflora and exotic S. apetala resulted in the highest soil MBC content in the one-year restored mangroves mixed with S. alterniflora.Both short-term invasion of S. alterniflora and mangrove restoration enhanced the total carbon pools of soil in comparison with mudflats.However, the S. alterniflora and young S. apetala forests had significantly lower soil TOC, POC contents, and total organic carbon pools than mature native mangroves.Biological replacement methods using mangrove species can enhance carbon storage in coastal wetland ecosystems.However, the native mangrove species are recommended.

Figure 1 .
Figure 1.Map of the study sites in Qi'ao Island in Guangdong, China.

Figure 1 .
Figure 1.Map of the study sites in Qi'ao Island in Guangdong, China.

Figure 2 .
Figure 2. Vertical variation of different carbon in the sediment at different sites.

Figure 2 .
Figure 2. Vertical variation of different carbon in the sediment at different sites.

ForestsFigure 2 .
Figure 2. Vertical variation of different carbon in the sediment at different sites.

Figure 3 .
Figure 3.Total soil organic carbon stock down to 50 cm at different sites (error bar represents standard deviation of total carbon stock, letters indicate significant differences at p < 0.05).

Figure 4 .
Figure 4. Redundancy Analysis (RDA) ordination showing the relationship between carbon contents of top 10 cm and environmental parameters of surface 0-10 cm soil (Sal means S. alterniflora age, Mangrove means mangrove age).

Figure 4 .
Figure 4. Redundancy Analysis (RDA) ordination showing the relationship between carbon contents of top 10 cm and environmental parameters of surface 0-10 cm soil (Sal means S. alterniflora age, Mangrove means mangrove age).

Table 1 .
Variation of soil physicochemical characteristics at different sites (n = 3, Mean ± SE.Different letters indicate significant differences at the same soil layer among sites).