The Impact of Sea Embankment Reclamation on Greenhouse Gas GHG Fluxes and Stocks in Invasive Spartina alterniflora and Native Phragmites australis Wetland Marshes of East China

The introduction of embankment seawalls to limit the expansion of the exotic C4 perennial grass Spartina alteniflora Loisel in eastern China’s coastal wetlands has more than doubled in the past decades. Previous research focused on the impact of sea embankment reclamation on the soil organic carbon (C) and nitrogen (N) stocks in salt marshes, whereas no study attempted to assess the impact of sea embankment reclamation on greenhouse gas (GHG) fluxes in such marshes. Here we examined the impact of sea embankment reclamation on GHG stocks and fluxes of an invasive Spartina alterniflora and native Phragmites australis dominated salt marsh in the Dongtai wetlands of China’s Jiangsu province. Sea embankment reclamation significantly decreased soil total organic C by 54.0% and total organic N by 73.2%, decreasing plant biomass, soil moisture, and soil salinity in both plants’ marsh. It increased CO2 emissions by 38.2% and 13.5%, and reduced CH4 emissions by 34.5% and 37.1%, respectively, in the Spartina alterniflora and Phragmites australis marshes. The coastal embankment wall also significantly increased N2O emission by 48.9% in the Phragmites australis salt marsh and reduced emissions by 17.2% in the Spartina alterniflora marsh. The fluxes of methane CH4 and carbon dioxide CO2 were similar in both restored and unrestored sections, whereas the fluxes of nitrous oxide N2O were substantially different owing to increased nitrate as a result of N-loading. Our findings show that sea embankment reclamation significantly alters coastal marsh potential to sequester C and N, particularly in native Phragmites australis salt marshes. As a result, sea embankment reclamation essentially weakens native and invasive saltmarshes’ C and N sinks, potentially depleting C and N sinks in coastal China’s wetlands. Stakeholders and policymakers can utilize this scientific evidence to strike a balance between seawall reclamation and invasive plant expansion in coastal wetlands.


Introduction
Global climate warming due to increased greenhouse gas (GHG) emissions from human activities is one of the most urgent challenges in ecology and climate change study [1]. Two of the most severe threats to world biodiversity are habitat alteration and the introduction and establishment of invasive species, especially in wetlands. Natural wetlands have been lost at a rate of 54-57% since pre-historic times, although it might have been as high as whereas tidal inundation enhances the net capture of CO 2 in a Phragmites salt marsh [41]. Plant invasions may alter the size and quality of the SOC pool, posing a threat to ecosystem function and the global carbon cycle [42]. Therefore, the time scale must be considered in assessing the relative constancy of C and N storage in coastal wetland soils [43]. Coastal wetlands and boreal peatlands are substantial C sinks (tons C/ha) according to global estimates of wetland pools. However, the lack of worldwide estimates per wetland type renders comparison to global C pool projections of saltmarshes, peatlands, and mangroves arduous [44]. The seemingly contradictory findings on GHG source-sink dynamics among coastal wetland ecosystems due to site-specific environmental circumstances [45] makes the quantification of the magnitudes and dynamics of coastal wetland GHG fluxes/stocks at regional and/or global scales essential [45,46].
Invasive alien plant species and coastal embankments have been studied in relation to the ecosystem dynamics of GHG fluxes and stocks in coastal wetlands. The diversity of plant species has a direct impact on GHG flux, as demonstrated by China's coastal wetlands, which pose a significant threat to global change [47]. The average CO 2 , CH 4 , and N 2 O fluxes (ecosystem respiration) through coastal wetlands of China are determined to be 388.76 ± 42.28 mg m −2 h −1 , 2.20 ± 0.31 mg m −2 h −1 , and 16.44 ± 2.96 mg m −2 h −1 , respectively [47]. Contrary to other reports, the restored wetlands were found to be sources of CO 2 release rather than sinks [25,27]. In a Typha plant-dominated mesocosm, CH 4 releases were linked considerably with soil carbon, nitrogen, and surface biomass, increasing methane emissions up to three times 45.9 ± 16.7 mg CH 4 -C m −2 h −1 under high water elevations (Lawrence et al., 2017 [48]), whereas in an estuarine wetland in Ohio, USA, it was 219.4 g m −2 yr −1 (Rey-Sanchez et al., 2018 [34]). These were significantly higher than average yearly 18.7, 3.26, and 5.26 g m −2 yr −1 CH 4 flow from S. alterniflora, P. australis, and C. malaccensis predominant flora in a tropical coastal estuarine wetland in China, suggesting that S. alterniflora's upper stalk, increasing CH 4 output, and extra efficient plant CH 4 delivery all contribute to the plant's higher CH 4 emission (Tong et al., 2012 [29]). Spartina soil had annual average CH 4 and N 2 O fluxes of 13.77 and 1.14 mol m −2 h −1 , respectively, which exceeds those of Kandelia, Sonneratia, and Cyperus, suggesting that the invasion of exotic wetland plants could transform coastal soils into a significant source of greenhouse gases while also increasing soil carbon deposition (Chen et al., 2015 [26]). Given the fact that CH 4 fluxes from mangrove, salt marsh, and seagrass environments vary greatly, they all function as net methane emitters., ranging from 67.33 to 72,867.83 mol CH 4 m −2 day −1 for mangrove; 92.60 to 94,129.68 mol CH 4 m −2 day −1 for marshland; and 1.25-401.50 mol CH 4 m −2 day −1 for seagrass [49]. Total soil N 2 O emissions from S. alterniflora communities in China are estimated to be approximately 0.06 Tg N 2 O y −1 , accounting for about 0.60% of total N 2 O emissions (9.6-10.8 Tg N 2 O y −1 ) (Gao et al., 2019 [50]).
Vast sea embankments have been built in S. alterniflora marsh to prevent its expansion and restore native P. australis marshes, which sea embankments have also confined as they grow inland [9]. Seawall-reclaimed wetland areas in the Yangtze River estuary increased by an average of 24,000 ha yr −1 , resulting in GHG flux rates of CO 2 and CH 4 that were 83.7 × 10 3 -164 × 10 3 mg m −2 yr −1 and 268-81.20 × 10 3 mg m −2 yr −1 , respectively [51]. Sea embankment reclamation significantly altered SOC and TN contents by 0.062% and 1.46%, respectively, compared to the values in the natural wetland of 0.053% and 1.30% (Zhou and Bi 2020 [52]); they decreased soil total organic C (233 to 495 g m −2 ) by 50.60% and total organic N (22 to 46 g m −2 ) by 49.99% (Yang et al., 2016 [53]), and they reduced soil TOC by 57% and TON by 59% in S. alterniflora salt marsh (Yang et al., 2017 [38]). A constructed seawall along the Coastal Line of Lake Nakaumi, Japan, influenced greenhouse gas emissions, with a high efflux of 245 mg CH 4 m −2 h −1 during the day and 725 mg CO 2 m −2 h −1 during the night (Hirota et al., 2007 [54]). According to a worldwide meta-analysis using a database of 209 sites to assess the impact of reclamation and land use cover changes of diverse wetland types, natural wetlands were found to be net sinks of atmospheric CO 2 and net sources of CH 4 and N 2 O [55].
Reclamation by coastal embankments may alter GHG dynamics, according to our hypothesis, thus affecting plant development and soil physiochemical characteristics in invasive S. alteniflora and native P. australis salt marshes. The aim of this research was to evaluate the following hypothesis and questions: (1) The construction of a sea embankment in local habitats will dramatically increase greenhouse gas fluxes between the soil and the atmosphere in invasive S. alteniflora and native P. australis marshes; (2) whether soil carbon stocks vary between S. alteniflora and P. australis marshes in their response to sea embankment reclamation; and (3) whether there are differences in GHG fluxes and stocks between seawall restored S. alteniflora and P. australis marshes and adjoining unrestored S. alteniflora and P. australis marshes.

Study Area
The research site is in Dongtai, Jiangsu Province, China, and is close to China's Yellow Sea Coast (Figures 1 and 2). With a shoreline of 95 km, Dongtai is situated at 32 • 37 -33 • 0 N. This is a transitional zone with a monsoon climate and an average yearly temperature of 15 • C, mean temperature range -0. 3 [56]. Mudflats in the low tide zone and salt marshes in the middle and high tidal zones make up Dongtai's estuary wetlands [57]. The offshore sediments of Dongtai are the largest in China and Jiangsu's most extensive land reserve. Their land area grew from 1268 km 2 in 1992 to 1300 km 2 in 1977. They enclose a sea body that is more than 200 km long (N-S) and 90 km wide (E-W). Even though actual evaluation planting on the offshore dunes failed, the 1996 breakthrough in seed dispersal and germination led to spontaneous spread throughout the 1 km wide coastal channel and was considered promising [58].
Since its introduction, the total area of S. alteniflora spread has increased dramatically. The plot, which was originally 3 km from the Dongtai-Dafeng border, had expanded into Dafeng tideland by 1992. S. alteniflora marsh with elevations ranging from 3.5 to 1.8 m has advanced about 1 km seaward, based on the 1.8 m elevation as the seaward marsh border. Its lush growth was determined to be best at heights of 2.0-2.5 m. The marsh was 1.1 km wide and covered 70% of the area [58]. Due to increased elevation, its landward portion was eventually replaced by Suaeda salsa. To combat the risk of Spartina alterniflora invasion, the Dongtai coastal administration began in 2011 erecting a 6163 km long seawall out of a total length of 15,339 km (Figures 3 and 4). The study area's soils are made up of marine and alluvial deposit materials, consisting of silt, sand, clay, and gravel, as well as a significant amount of organic matter. Its soil and natural plants were in fair shape at the time of enclosing. Phragmites australis, Suaeda salsa, Spartina alterniflora, and Imperata cylindrica var. covered approximately 80% of the soil [58].  Since its introduction, the total area of S. alteniflora spread has increased dramatically. The plot, which was originally 3 km from the Dongtai-Dafeng border, had expanded into Dafeng tideland by 1992. S. alteniflora marsh with elevations ranging from 3.5 to 1.8 m has advanced about 1 km seaward, based on the 1.8 m elevation as the seaward marsh border.

GHG Sampling and Flux Measurement
From January to September 2020, closed static chambers were used to monitor CO 2 , CH 4 , and N 2 O fluxes eight different times to compare P. australis and S. alteniflora GHG emissions. Gas samples were taken from both sites simultaneously using the enclosed chamber method, and all sampling was completed 2 h before the least tidal cycle (G. C. Chen et al., 2010 [59]). Stainless steel collars of diameter 20 cm were driven 5 cm deep into the soil for at least 36 h before taking GHG readings. On the cylinder's side surface, a metal exhaust pipe of 1 mm diameter was installed (approximately 8 cm above the soil) to enable normalization of pressure disturbances between the chamber and surrounding environment. To produce a gas tight grip around the exhaust pipe, silicone adhesive was used. In January/February, three observations were taken on vegetated sediments and two observations taken on natively unvegetated sediments for both restored and unrestored marshes. Measurements were collected three times a month on both vegetated and unvegetated sediments in populations of both plant species from March to September. We assessed respiration without photosynthesis and saw changes in N 2 O and CH 4 fluxes attributable to plant-catalyzed gas transfer by vegetation cover in unvegetated soils that were free of vegetation cover. Unvegetated soils were incidentally identical to vegetated strata and were as near as 1 m apart. Because the study marsh's plant density was low, we could insert gas flux flanges in between plants and take unvegetated observations within plant stands. With a diameter of 3 cm, the gas sampling port was located on top of the chamber. The hole was carefully covered with a plastic cork and closed with silicone glue once the chamber was lowered into the soil. Gas samples were obtained by piercing the plastic with the connected injectable needle into the chamber with a 10 mL glass syringe. An electric fan was installed in the chamber to ensure that the inside air was well mixed. We kept the air temperature within the chambers under control during the summer observations by enclosing the chamber top with cotton covers.
Gas samples of 25 mL were taken in duplicates immediately after securing the chamber for four more time periods, ranging from 30 to 60 min, by incorporating a 60 mL plastic syringe with a faucet through a plastic tube placed in the cylinder's side. The 25 mL gas samples were instantly introduced into the 12 mL Exetainer vials (Labco, CN) following the use of ultra-pure helium (UHP He) as a flush and then drained. To evaluate the soil-atmosphere fluxes of three greenhouse gases, we used the gas chromatography (GC) approach to assess concentrations of GHG within 30 days. Based on standards, the injector, column, and detector had temperatures of 100 • C, 70 • C, and 320 • C, accordingly sensitive to low levels of N 2 O (Y. Chen et al., 2015 [26]). CH 4 and CO 2 concentrations were measured using the same GC system with a thermal conductivity detector (TCD). With nitrogen as the carrier gas, at a 20 mL min −1 flow rate, temperatures of 80 • C, 50 • C, and 80 • C were maintained for the injector, column, and detector, respectively (G. C. Chen et al., 2011 [60]). By correlating the peak area of gas samples to that of standard curves, the concentration of greenhouse gases was calculated. CO 2 , CH 4 , and N 2 O had relative standard deviations of 2.8%, 2.7%, and 2.4%, respectively, in replicate standard measurements (Y. Chen et al., 2015 [26]). Each greenhouse gas's soil-atmosphere flux (F, mol m −2 h −1 ) was calculated using the following formula (Chen et al., 2015) [26], Chen et al., 2016 [61]): where V (m 3 ) is the volume of air remaining inside the closed chamber after lowering the open end into the soil; ∆G c (h −1 ) is the variation in gas concentration inside the chamber; A (m 2 ) is the surface area of soil enclosed by the chamber; and ρ (m 3 mol −1 ) is the volume of gas at air pressure. N 2 O and CH 4 GHG (CO 2 -eq) fluxes (Fe, mg CO 2 m −2 h −1 ) were calculated using the following formula (Chen et al., 2015 [26]): where F GHG is the greenhouse gas flux, M GHG is the greenhouse gas molecular weight, and GWP (GHG) is the GHG global warming potential (with global warming potential of 23 and 296 CH 4 for and N 2 O, respectively, on a time horizon of 100 years) (IPCC 2014 [62]).

Biomass and Soil Sampling
Following the completion of gas sampling, biomass and soil sampling were conducted. Three 50 cm by 50 cm quadrats were established in each transect of the coastline with restored and unrestored S. alteniflora and P. australis marshes to collect plant leaves, stems, and litter, and three soil segments (10 cm length by 10 cm breadth by 60 cm depth) were unearthed to extract roots. Four 40 m by 40 m transects were deployed in October 2020 in the seawall restored S. alteniflora marsh (RSA), P. australis salt marsh (RPA), unrestored S. alteniflora salt marsh (USA), and P. australis salt marsh (UPA) ( Figure 1). According to assessments of Thematic Mapper satellite photos, the restored and unrestored S. alteniflora salt marshes in the survey area infringed on mudflats for a significant number of years (Chung et al., 2004 [58]). The density and height of the plants were measured on-site. Separate samples of aboveground biomass (AGB), belowground biomass (BGB), and litterfall biomass (LB) were taken. The AGB included the overall weight of leaves, branches, and stems. The sections for the BGB quantification were 1 m × 1 m in size, with the depth of each section varying depending on the rooting depth.Phragmites australis roots are mainly found at ≥60 cm depth, and that of other species shallower in depth (Lu et al., 2019 [36]). The soil in each plot was scooped and gently cleaned through a fine screen to gather all the roots. The entire litter weight within the plot was included in the LB. Employing a balance (BSA124S-CW Sartorius), the fresh weight of the AGB, BGB, and LB parts from each quadrat were measured independently. For further analysis, each component was sub-sampled to roughly 100 g fresh weight. To establish a consistent weight, samples were oven-dried for at least 72 h at 60 • C, and their dry weight quantified to ascertain the water content. Each transect was divided into three 2 m by 2 m plots, with three spots chosen at random for sampling soil at depths of 0-10 cm, 10-20 cm, and 20-30 cm within every plot, and all soil samples were well mixed to obtain a blend for each section.

Laboratory Analyses
After filtering each soil sample through a 100-mesh sieve and properly flushing with water, the residual roots were amassed. For biomass measurement, a total of nine 1 m by 1 m plots were randomly chosen and garnered for each of the two main species. Fresh samples were properly cleaned and oven-dried for at least 72 h at 60 • C to quantify leaf, stem, litter, root, and total biomass. Forceps were used to extract organic materials from soil samples, which were then sieved at 2 mm and oven-dried to a consistent weight at 50 • C to evaluate the soil water content. The dry weight of the total sample span was divided by the volume of the core section to determine bulk density for each sample. An elemental analyzer (elementar UNICUBE) was used to determine the carbon and nitrogen content. Based on the above-measured carbon, nitrogen, and biomass values, the carbon pool of each portion was determined. A glass membrane electrode was used to determine soil pH in combination with a 1:2.5 soil-to-water ratio. A conductivity meter was used to determine the soil salinity combined with a 1:5 soil-to-water ratio (W. Yang, Li, et al., 2016 [9]).

Data Analysis
Stocks of organic and inorganic carbon in the sediments were examined following a prior evaluation [63]. The accrued carbon stocks in all layers were used to calculate the pools of inorganic and organic carbon (Mg C ha −1 ) in 1 m soil strata, and the carbon pool in each stratum was calculated as the product of carbon metadata (% dry weight) and dry bulk density (g dry soil mL −1 ). To calculate the quantity of CO 2 net sequestered, estimates of both the soil inorganic and organic carbon stocks were utilized. In soil formation, the net total CO 2 flux from the air to the soil can be computed using the following formula (Lu et al., 2019 [36]): where C OC is the organic carbon density in the sediment (mol C mL soil −1 ); C IC is the inorganic carbon density in the sediment (mol C mL soil −1 ); Ψ is the predicted CO 2 and CaCO 3 gas exchange reaction ratio; and CO 2 netseq refer to the moles of CO 2 net sequestered in mL of soil. The value of Ψ = 0.6 for shallow coastal sites was used in this estimation [64,65]. All data were reported as the average of three replicates' averages and standard deviations. Prior to testing with analysis of variance (ANOVA), the data for each variable were examined for homogeneity and then log-transformed if necessary. Using ANOVA and Tukey HSD multiple comparisons, the differences in plant characteristics and carbon content of the AGB, BGB, and LB between the two plant marshes were investigated. The carbon fluxes and pools (AGB, BGB, and LB) were then estimated using the biomass and carbon content data. SPSS version 16.0 was used for all statistical analyses (SPSS Inc., Chicago, IL, USA). Obtained carbon storage and flux data were then compared to the carbon fluxes and stocks published in Web of Science relevant papers.

Plant Biological Characteristics and Biomass Carbon Pool
S. alteniflora was taller (F = 22.6, p = 0.001) and had a higher plant density (F = 12.5, p < 0.001) than P. australis. Plant root (0-60 cm soil depth) and stem biomass were considerably lower in the restored P. australis salt marsh compared to the unrestored P. australis salt marsh, whereas AGB and total biomass were significantly lower in the restored S. alteniflora salt marsh than in the unrestored S. alteniflora salt marsh (Figure 5a,b). Conversely, litter biomass was much higher in the restored P. australis salt marsh compared to the unrestored. However, there were no statistically significant variations in leaf biomass and total biomass between the unrestored and restored P. australis salt marshes ( Figure 5b). Overall, the unrestored S. alteniflora marsh had the most AGB, BGB, and LB. Both AGB and LB carbon content were identical for both plants (both p < 0.001), whereas the BGB carbon content was significantly higher for UPA (31.16 ± 0.8%) and RPA (29.34 ± 0.8%) than for USA (9.94 ± 1.7%) or RSA (7.46 ± 1.8%). This suggests that P. australis roots play a significant role in carbon sequestration, regardless of the presence of the sea embankment. Overall, the highest AGB carbon content was found in unrestored S. alteniflora, the highest BGB carbon content was found in unrestored P. australis, and the highest LB carbon content was found in unrestored S. alteniflora (p < 0.001). The data show that the sea embankment affects the carbon sequestration potential of both P. australis and S. alteniflora plants, notwithstanding the slight variations in the restored and unrestored biomass carbon content marshes.

Soil Carbon C and Nitrogen N Pools
Both water content and bulk density showed substantial correlations between ha tats and core depths. The P. australis habitats had much lower mean water content th the S. alteniflora marsh. In both restored and unrestored marshes, water content increas steadily with increasing core depth (Table 1), whereas bulk density increased significan with core depth for the first 40 cm in all marshes. Reclamation by sea embankments h a significant impact on SOC, SIC, and total N ( Table 2). SOC, SIC, and total N concent tions were positively correlated across species in unrestored and restored salt marsh and all had significant positive correlations with soil moisture and salinity (Table  Nonetheless, there were substantial negative linkages observed between varying orga carbon and nitrogen components in the soil, as well as soil bulk density and pH (Table  Accumulation of SOC, SIC, and total N were much lower in the restored S. alteniflora s marsh than in the unrestored. In the unrestored and restored S. alteniflora salt marsh the SOC stock varied from 397 to 994 g m −2 and 398 to 633 g m −2 , respectively, whereas the P. australis marshes, it varied from 438 to 831 g m −2 and 247 to 496 g m −2 , respectiv

Soil Carbon C and Nitrogen N Pools
Both water content and bulk density showed substantial correlations between habitats and core depths. The P. australis habitats had much lower mean water content than the S. alteniflora marsh. In both restored and unrestored marshes, water content increased steadily with increasing core depth (Table 1), whereas bulk density increased significantly with core depth for the first 40 cm in all marshes. Reclamation by sea embankments had a significant impact on SOC, SIC, and total N ( Table 2). SOC, SIC, and total N concentrations were positively correlated across species in unrestored and restored salt marshes, and all had significant positive correlations with soil moisture and salinity (Table 3). Nonetheless, there were substantial negative linkages observed between varying organic carbon and nitrogen components in the soil, as well as soil bulk density and pH ( Table 3). Accumulation of SOC, SIC, and total N were much lower in the restored S. alteniflora salt marsh than in the unrestored. In the unrestored and restored S. alteniflora salt marshes, the SOC stock varied from 397 to 994 g m −2 and 398 to 633 g m −2 , respectively, whereas in the P. australis marshes, it varied from 438 to 831 g m −2 and 247 to 496 g m −2 , respectively (Figure 5e,f). In the S. alteniflora marsh, SIC stock varied from 1108 to 1786 g m −2 USA and 485 to 914 g m −2 RSA, respectively, whereas in the P. australis marshes it varied from 1007 to 1650 g m −2 UPA and 309 to 622 g m −2 RPA, respectively. Nitrogen N content varied from 33 to 91 g m −2 in the USA, 24 to 59 g m −2 in the RSA, whereas for P. australis, it varied from 45 to 88 g m −2 UPA and 36 to 74 g m −2 RPA, respectively (Figure 5e,f). In the restored S. alterniflora salt marsh, SOC, SIC, and N stocks in the 0-60 cm depth declined by 20.5%, 52.9%, and 41.2%, respectively, compared to the unrestored marsh (Figure 5g,h), whereas it declined by 50.8%, 58.9%, and 11.4%, respectively, in the P. australis marshes.
With the exception of SOC in the S. alteniflora 20-60 cm depth and N stocks in the P. australis 20-40 cm soil-depth range, significant variations in soil (0-60 cm depth) SIC and TC stocks were identified in both the unrestored and restored salt marshes of both plants. Significantly low TOC ratios in the USA and RSA marshes and low TN ratios in the UPA and RPA marshes show that the SOC content of the S. alteniflora salt marsh and the N content of the P. australis wetlands were not significantly influenced by reclamation seawalls. Overall, our findings revealed that habitat and core depth significantly impacted inorganic carbon, organic carbon, and total nitrogen ( Table 1). For organic carbon and total nitrogen, there was a strong link between biome and core depth (p < 0.05), and the reverse was true for inorganic carbon content (p > 0.05). Table 1. Soil physical properties and carbon and nitrogen content (mean ± SE, n = 54) in the unreclaimed and seawall-reclaimed S. alterniflora and P. australis salt marshes. Lower-case superscript letters ( a,b ) indicate statistical significance at p < 0.05 between the unreclaimed and sea embankment-reclaimed salt marshes in the same community and at the same soil depth.

Greenhouse Gas (GHG) Fluxes
Both vegetation type and reclamation forms significantly impacted greenhouse gas fluxes between the soil and the atmosphere (Table 4). CO 2 fluxes for both species were similar, whereas there were differences between embankment restored and unrestored segments and over the span of the year. The atmospheric flux of all three GHGs was positively correlated with soil temperature, water content, biomass, TOC, TIC, and TN but negatively correlated with soil pH (Table 5). Slightly above a quarter (27%) of the fluxes recorded were nil or beyond observation limits. CO 2 fluxes were either zero or positive in January to May of both transect sections, at which photosynthesis is limited by composition, indicating net CO 2 release (Figure 6a,b).    In the months of January to May, CO 2 fluxes from embankment restored parts were also either zero or positive, identical to unrestored sections, but immediately AGB became significantly available for photosynthesis to counteract respiration rates, and CO 2 fluxes turned negative (Figure 6a,b). CO 2 fluxes from restored sections were much higher in April, July, and September than fluxes from the restored sections in June, July, and September. In July, CO 2 uptake in restored segments for both plants were significantly higher (13.98 m mol m −2 h −1 RSA and 11.34 m mol m −2 h −1 RPA) than in any other months. According to Pearson correlation analysis, the most critical determinants of CO 2 fluxes were live AGB and soil temperature, which combined accounted for over half of the variance in CO 2 flux data ( Table 5).  (Figure 6e,f). Although seasonal fluctuation is still evident, the importance of sediment temperature in the correlation analysis shows that other seasonal variations would be expected. Summer months showed a modest net positive emission rate of 105 ± 48 n mol m −2 h −1 , whereas cooler months had a lower net absorption rate of −64 ± 47 n mol m −2 h −1 (mean ± standard error) in the unrestored marshes. The results of the correlation analysis indicate that sediment temperature and live BGB influence N 2 O flux rate (Table 5). Overall, these indicators predicted about half of the variation in N 2 O fluxes (49%).

Spatial Variations and Determinants of Soil Carbon and Nitrogen Stocks
The robust relationships observed in both restored and natural wetlands between SOC and plant community cover, overall vegetation cover, BGB, soil C:N, and soil salinity support the theory that all have a significant impact on SOC [52,66,67]. Findings show that in the 0-60 cm soil depth of S. alteniflora salt marsh, sea embankment reclamation reduced SIC by 33.6% and total carbon TC by 30.5%, whereas in the P. australis marsh, reductions were 15.4% and 23.5%, respectively (Figure 5g,h). Decrease in dry weight carbon content was slightly significant, 25% in the BGB of the S. alteniflora marsh compared to the other biomass categories (Figure 5g,h). Plant-residuals input and SOM degradation are well-known factors influencing sequestration of soil organic carbon and nitrogen [9,52]. Preceding research found that embankment walls can limit S. alteniflora development and inhibit tidal flooding flows and salt marsh inundation caused by precipitation or irrigation. As such, S. alteniflora is hampered from completing its regular life cycle in freshwater habitats, whereas mild salinity promotes its growth by accelerating flowering [68,69]. Our study corroborates these findings, as shown in Table 1 and Figure 5a-d. In the S. alteniflora salt marsh, reclamation by seawalls greatly degraded soil salinity (0-60 cm depth) from 1.52% to 0.40%, LB from 1012 g m −2 to 569.6 g m −2 , and root biomass from 1078 g m −2 to 551.3 g m −2 . In comparison to native P. australis, Scirpus mariqueter, S. salsa, and cyperus marshes of eastern China, S. alteniflora invasion often promotes soil organic carbon and nitrogen sequestration by enhancing LB and BGB feeds and decreasing SOM breakdown [9,26,69,70].
The lesser litter production and root biomass resulted in a considerable decline in soil SOC and TN in the embankment-restored marsh (Table 1). There are two possible explanations for these alterations. First, the recovered wetland had decreased production recovery and yield rates compared to native intertidal wetland because it is a relatively dry habitat [71]. Furthermore, because of the increased rate of organic carbon sequestration in wetland soils, native marsh can influence prolonged carbon capture [72]. Second, the restored S. alteniflora salt marsh diminished LB and BGB soil inputs, resulting in soil organic carbon and nitrogen stocks degradation (Figure 5e-h).
Additionally, because of the hypoxic soil conditions that support lengthy deposition of SOM, marshes with high water content or water tables are favorable for SOM accumulation [9,73]. Recurrent tidal flooding of salt marshes can be hampered by sea embankments, resulting in significant declines in soil water levels [23]. The soil moisture in the restored S. alteniflora salt marsh dropped dramatically, from 46.67% to 27.752% (Table 1), potentially increasing SOM breakdown and speeding up the deficiency of soil organic carbon and nitrogen. This could speed up the degradation of SOM and the loss of soil organic carbon and nitrogen [9]. The average SOC and TN were lower in the embankment-restored marsh transects than in the restored marsh (Table 1), which could elucidate why the seawallrestored marsh has more spatial variability of plant species. The mean values of soil SOC and TN may be affected by habitat degradation in an embankment-restored wetland. These findings also emphasized the necessity of characterizing regional patterns of soil characteristics and uneven vegetation.

Environmental Drivers of Greenhouse Gas (GHG) Fluxes
Although China's wetlands constitute a net carbon sink, they also release CO 2 and CH 4 from the soil, making them a substantial source of atmospheric carbon. However, the volume of emissions differed greatly depending on the type of wetland. The tropical location, expanse of mangrove vegetation, periodic shoreline flooding, and moderate weather are all factors that contribute to the highest carbon emissions in coastal wetlands. Higher elevations and latitudes decreased soil CO 2 and CH 4 emissions in marsh wetlands and may be responsible for the low annual temperature and a delayed growing season for grassland plants [74]. Furthermore, the shortage of aquatic plants and low CO 2 and CH 4 emissions in river and lake wetlands could be linked to their lengthy flooding regimes [75]. In addition to the impact of photosynthesis on fluxes of CO 2 , we observed significant fluctuations in N 2 O fluxes across the restored sections of the S. alterniflora and P. australis marshes (Figure 6e,f). Though unrestored sections had more live belowground biomass, the matching of biogeochemical properties was attributable to the plants' closeness to the unvegetated plots. Because the existence or absence of plants had no effect on observations, CH 4 or CO 2 fluxes were attributed to biomass characteristics. As a result, we infer that wetland plants influenced GHG emissions by modifying the soil instead of plant-induced diffusion [24]. Live AGB was negatively linked to fluxes of CO 2 (i.e., the higher the biomass, the more negative the flux) and strongly linked to soil temperature in the correlation analysis. Although live AGB was the most critical indicator of GHGs (p < 0.001), sediment temperature was another critical predictor, particularly in unrestored areas, according to correlation analysis.
The highest CO 2 emissions obtained (587.4 mg m −2 h −1 ) in the unrestored S. alteniflora marshes occurred in later summer month (Figure 6a), which coincided with the hottest weather and driest sediments (Table 4). This value far exceeds the average CO 2 emissions (388.76 ± 42.28 mg m −2 h −1 ) across coastal wetlands in China [47], indicating the seawall embankment played a significant role in carbon emissions. However, [24] observed that invasive P. australis marsh emits up to 660 mg m −2 h −1 CO 2 , whereas the unvegetated marsh could emit up to 880 mg m −2 h −1 of CO 2 . Thus, net CO 2 uptake in embankment restored areas was highest during the time of intense biomass production in June and July but fell significantly in August and September (Figure 6a,b). This occurrence is influenced by a peak in ecosystem respiration and delayed photosynthetic uptake senescence.
The considerable spatial-temporal changes in CH 4 flux were explained by a positive correlation between CH 4 flux and temperature and water content in both reclamation types (Table 5). Furthermore, the positive correlation between CH 4 emissions and increased biomass is explained by plant-mediated transport and root exudation-stimulated methanogens [76]. The CH 4 fluxes in this study ranged from −3.09 to 200.32 µ mol m −2 h −1 , which is within the typical range of 2200 ± 310 µ mol m −2 h −1 of methane CH 4 emissions across China's coastal wetlands. Our results are similarly consistent with data from other salt marshes such as 4.8-70.5 µ mol m −2 h −1 [29]; 97.5-313.75 µ mol m −2 h −1 [77] for mangrove forests; 2.00 to 21.70 µ mol m −2 h −1 [26]; 11.9-5168.6 µ mol m −2 h −1 in Hong Kong [59]; and 5.73-201.88 µ mol m −2 h −1 in East India [78]. By discharging exudates and producing plant litterfall, soil CH 4 emission represents a source of methanogenic substrate [29] and can be controlled by vegetation in three ways, as exemplified in [79]. The marginally significant difference observed in the regression between biomass and CH 4 emissions for the embankment restored and unrestored sections explain how plant factors such as root exudation are most likely responsible for some disparities in CH 4 flux rates across marshes. Because they promote favorable circumstances for nitrogen cycling activities, salt marshes are essential nitrogen sinks in coastal zones. Most marshes, however, are minor N 2 O sinks because denitrifiers can use N 2 O when nitrate is scarce [80]. Because of the increased availability of energetically more suitable nitrate due to N-loading, marshes might suddenly shift from net sinks to net sources of N 2 O [81]. Such was observed in this study, as shown in Figure 6e,f. Between January and July, N 2 O emissions were consistent between the restored S. alterniflora and P. australis marshes. However, an abrupt shift was observed in the P. australis marsh, as N 2 O emission rose from 16.21 n mol m −2 h −1 in July to 534.16 n mol m −2 h −1 in August, compared to −88.44 n mol m −2 h −1 to 80.27 n mol m −2 h −1 in the restored S. alterniflora marsh. These findings corroborate the evidence that S. alteniflora relatively lowers N 2 O emission compared to P. australis, and thus can be regarded as a species for tidal zone stabilization [82]. The peak N 2 O emission (534.16 n mol m −2 h −1 ) obtained in this study exceeds the average N 2 O fluxes (373.6 n mol m −2 h −1 ) estimated across coastal wetlands in China [47]. The contribution of P. australis rhizosphere oxygenation to N 2 O flux rates has also been substantiated in a freshwater wetland investigation [83]. As a result, it is reasonable to infer that embankment reclamation significantly increases N 2 O emissions in a P. australis salt marsh as observed in this study. As a result of the increased CO 2 content in the chamber air, the CO 2 gradient between the soil and the chamber is reduced, lowering the detected fluxes. Lower fluxes could be attributed to a cooler soil temperature inside the static chamber. During the experiment, soil temperature was monitored hourly outside the chambers, which was approximately 4 • C higher than the temperature within the static chambers [84]. Such diurnal variations in air temperature affects GHG estimates. The marshy environment influenced chamber positioning relative to plant tissue, which influenced the interpretation of flux data, especially in the case of CO 2 , where perhaps microbial respiration and also root and shoot respiration and photosynthesis must be properly balanced.
Owing to the anoxic and oxic conditions influenced by the above noted factors, marshes are most certainly the most significant source of soil carbon emissions in China's wetlands [74]. China's wetlands emit less CO 2 than the United States, but more than Canada, Europe, and Russia, as well as the global mean. Conversely, soil CH 4 in China's wetlands (0.29 Mg C ha −1 year −1 ) is considerably less than those in America, Canada and the global average, but higher than that of Europe [74]. These variances could be explained by changes in climate and soil conditions between China and these other nations, resulting in disparities in the distribution patterns of each wetland category.

Invasion/Reclamation Context
In China, the pace of seawall construction and expansion for the primary goal of coastal land reclamation [5] and control of Spartina alterniflora invasion [85] has more than twice expanded in the past decades [86]. This resulted in the loss of at least 10,520 km 2 of habitats in the Yellow Sea [7] and triggered other ecological impacts of climate change, GHG emissions, marine pollution, etc. Our study is the first to concurrently examine the effects of a seawall embankment on GHG fluxes and stocks in S. alteniflora and P. australis dominated saltmarsh. Compared to a study that examined the patterns and drivers of greenhouse gases and flux in china's coastal wetlands [47], the average emissions of carbon dioxide (CO 2 ), methane (CH 4 ), and nitrous oxide (N 2 O) were estimated to be 388.76 mg m −2 h −1 , 2.20 mgm −2 h −1 , and 16.44 µgm −2 h −1 respectively. Although [47] did not differentiate nor compare between reclaimed and unreclaimed wetlands, their estimates were relatively less than the averages found in the reclaimed marshes of this study (404.76 mg m −2 h −1 , 3.06 mgm −2 h −1 , and 23.38 µgm −2 h −1 , respectively). This essentially justifies our hypothesis that seawall embankment alters GHG emissions in reclaimed wetlands, because no study assessed these effects in the area. Studies by [9,52,85] examined only the impact of seawall embankment reclamation on soil organic carbon and nitrogen pools. In all these scenarios, the seawall reclamation essentially altered soil C and N accumulation, salinity, soil moisture, total biomass, and microbial biomass, especially in the invasive plant community, as observed in this study. S. alteniflora invasion has been found to cause a shift in methanogen community composition and CH 4 production capability [87], increasing CH 4 emissions in China [29]. Hence, the erection of coastal embankments, which is becoming a more prevalent method of restricting the spread of Spartina alterniflora, provides an opportunity to reduce CH 4 emissions, as observed. However, our findings on GHG fluxes show that embankment reclamation significantly increased CO 2 and reduced CH 4 emissions in both reclamation types and increased N 2 O emissions in the restored P. australis marsh in the warmer months. As for N 2 O emissions, the potential of embankment reclamation to increase emissions in a native marsh, particularly in the warmer months, is attributable to the significance of rhizosphere oxygenation in P. australis on N 2 O flux rates [83]. Hence, it is imperative to note that these findings may not be indicative of other parts of the world because the functions of invasive and native species are opposite in those areas; hence further investigation is needed.

Management Implications
In an invasive plant-dominated salt marsh, seawall construction can have negative implications such as decreased soil organic carbon and nitrogen buildup, which can lead to weak soil nitrogen nitrification/denitrification, among other effects [9,52,85]. In contrast, salt marshes, when preserved in front of seawalls, levees, and other "grey" coastal defenses, can offer protection strata that significantly mitigate floods and storms [88][89][90]. Paddy fields and seawalls can help to mitigate the economic damage inflicted by severe storms [91], restoring wetlands beyond the seawalls, which provides a healthier water habitat [92], whereas open dikes protect native salt marshes, resulting in almost no loss of ecosystem functions and services [6]. The eradication of S. alteniflora resulted in reduced CO 2 , CH 4 , and N 2 O emissions in areas it earlier occupied, whereas CH 4 fluxes, in contrast, increased dramatically in places where the native species P. australis was maintained [93]. Similarly, reclamation of a S. alteniflora dominated salt marsh by embankment seawalls significantly reduced soil total organic carbon, total organic nitrogen, soil salinity, litter, moisture, and root biomass effects [9,52,85]. In this study, sea embankment enhanced CO 2 emissions and decreased CH4 emissions in the Spartina alterniflora and Phragmites australis marshes, respectively. In addition, the coastal embankment wall increased N 2 O emissions in the Phragmites australis salt marsh while decreasing emissions in the Spartina alterniflora marsh. Methane CH 4 and carbon dioxide CO 2 fluxes were similar in both restored and unrestored portions, whereas nitrous oxide N 2 O fluxes were significantly different owing to enhanced nitrate due to N-loading. In both plants' marsh, sea embankment reclamation diminished soil total organic C and total organic N, as well as plant biomass, soil moisture, and soil salinity. The three gases studied have their fluxes and stocks varied in both reclamation types and were linked to soil temperature, reclamation type, anoxic conditions, and living biomass. Seasonal variations in CO 2 flows were observed, as expected. As predicted in salt marshes, CH 4 production rates were moderate, but they varied substantially, indicating that CH 4 generation is highly spatially variable. The GHG stock soil measurements were obtained during the autumn period of the season, which could limit our conclusions. As such, further studies should consider the spring, summer, and winter seasons, such that spatiotemporal variations of GHG stocks in seawall reclaimed wetlands relative to alien plants invasions will yield more concrete conclusions.
Our findings corroborate the evidence that coastal reclamation by seawalls could curb S. alteniflora invasion while also endangering biodiversity, risks in terms of climate change mitigation, and long-term sustainability of coastal wetlands. Stakeholders in charge of ecosystem eradication initiatives and other types of restoration must therefore define their objectives and criteria for success, as evaluating the implications of coastal reclamation management options in the context of global climate change mitigation becomes more challenging in today's threatened ecosystems. At both the national and municipal levels, thorough wetland control mechanisms are required to accomplish a goal of "zero net loss" and long-term sustainability of coastal wetlands. Offsetting anthropogenic disruption and restoring damaged wetlands can enable sustained C budgets and spur global efforts to combat global warming and climate change. These findings help move us closer to the UN's goal for a steady transition to a carbon-free future and eventual GHG pollution eradication. This scientific information offers stakeholders and policymakers leverage to balance seawall reclamation and invasive plant spread in coastal wetlands. There are opportunities for further research on balancing the tradeoffs between the synergistic impacts of reclamation by invasive plant species and reclamation seawalls.