Deposition Flux, Stocks of C, N, P, S, and Their Ecological Stoichiometry in Coastal Wetlands With Three Plant Covers

The depositional flux of coastal wetlands and the deposition rate of biogenic elements greatly affect the carbon sink storage. Ecological stoichiometry is an important ecological indicator, which can simply and intuitively indicate the biogeochemical cycle process of the region. This study investigated the soil deposition flux, stocks, and ecological stoichiometric ratios of C, N, P, and S under different water and salt conditions based on 137Cs dating technology in the Yellow River Delta (YRD) of China. The results showed that the deposition fluxes were 0.38 cm/year for PV wetlands, 1.08 cm/yr for PA wetlands, and 1.06 cm/yr for SS wetlands. Similarly, PA wetlands showed higher deposition fluxes of C, N, and S compared with SS and PV wetlands. PA wetlands had higher stocks of C (5.86 kg/m2), N (0.36 kg/m2) and S (0.36 kg/m2) in the top 1-m soil layer compared with PV and SS wetlands. However, the highest deposition rate of P (9.82 g/yr/m2) was observed in SS wetlands among the three wetlands. Three accumulative hotspots of C, N, and S in soil profiles of PA and SS wetlands were observed at soil depths of 0–10, 40–60, and 90–100 cm, whereas one accumulative hotspot of P was at the soil depth of 10–12 cm in SS wetlands and 80–82 cm in PA wetlands. PV wetlands showed higher accumulations of C, P, and S in the top 10 cm soil layer and N at the soil depth of 90–100 cm. The higher top concentration factors in these three wetlands indicated that the dominant input of plant residues was the main reason. The ratios of C/N and C/N/P of each sampling site were higher in the surface soils and decreased with depth. The ratios of C/P and N/P were larger in the surface layer (0–20 cm), the middle layer (40–60 cm), and the deep layer (90–100 cm). The ratios of N/P and C/N/P were relatively lower, indicating that these studied wetlands were N-limited ecosystems. The results implied that the coastal wetlands in the YRD have huge storage potential of biogenic elements as blue carbon ecosystems.


INTRODUCTION
Coastal wetlands, sea-land interlaced belts and buffer zones are important areas for the strong exchange of materials and energy between the ocean and land. They play a vital role in regulating global environmental changes and serving as sources and sinks of carbon (C), nitrogen (N), phosphorus (P), and sulfur (S). C, N, P, and S are important constituent elements in the soil. Their contents and stocks can directly affect the soil quality, nutrient cycling, and ecological functions of the coastal wetland ecosystems Feng et al., 2017;Lu et al., 2018).
The regulation of biological elements with different soil profiles in the coastal wetland is not obvious because coastal wetlands are subject to the combined effects of erosion and sedimentation under the two-way regulation of tidal current and runoff. Taking soil organic carbon (SOC) as an example, the SOC pool is composed of the active pool and the inert pool, which leads to different stability in the soil profile due to high spatial heterogeneity. In addition, soil animals and microorganisms could decompose SOC and nitrogen into CO 2 , CH 4 or N 2 O, which would exacerbate the heterogeneity of soil carbon or nitrogen pools (Wu et al., 2013. Many studies have shown that SOC is mainly distributed in the upper soil layer of the soil profile and it tends to gradually decrease from the surface to the bottom (Wang et al., 2010). However, the current research on the SOC storage of coastal wetlands is mainly concentrated on the surface soil (0-30 cm) of the Yellow River Delta (YRD) region, and there are few studies on the contribution rate of the deep SOC storage to the carbon pool (Luo et al., 2020). Therefore, it is not enough to accurately reflect the soil carbon pool of coastal wetlands by studying the biogenic elements in surface soils. For the study of the distribution form of soil biogenic elements and ecological stoichiometry in the soil, a deeper soil profile should be selected.
The ecological stoichiometric ratio provides important insights for studying the energy and material cycles in estuarine ecosystems (Coynel et al., 2016). In addition, each ecological stoichiometric ratio is of great significance to the biogeochemical cycle of C, N, P, and S in the soil. For example, the soil C/N ratio is a sensitive indicator used to measure soil quality (Elser et al., 2003); the C/P ratio and C/S ratio are considered to be symbols of the mineralization ability of phosphorus and sulfur in the soil. Previous studies have focused on the stoichiometric ratios of soil and the effects of various environmental factors, such as vegetation types and rainfall, on the ecological stoichiometric ratios of soil at different soil depths in forests, karst, and grasslands (Yang et al., 2014;Fan et al., 2015;Wang et al., 2018). Some researchers have presented the impacts of the ecological stoichiometric ratios on various ecological functions in the ecosystems and the interactions between microbes and the ecological stoichiometric ratios in the abovementioned ecosystems. Although the stoichiometric ratios in coastal wetlands are more important due to the impacts of historical sedimentation and strong exchange between marine and terrestrial materials of coastal wetlands, only few studies have been conducted on the relations between sedimentation and ecological stoichiometry in coastal wetlands.
The Yellow River is one of the rivers with the highest sediment content in the world, and the relationship between its sedimentation flux and the stocks of soil biogenic elements is also a problem that needs to be revealed. As the most fragile and sensitive new delta in the world, the YRD has extremely precious ecological value. In recent years, governance along the Yellow River has also received great attention from the Chinese government. Therefore, it is necessary to study the soil deposition rate of the YRD and enhance the sequestration capacity of these elements and soil quality management of the coastal wetlands. The primary objectives of this study were to: (1) analyze the depth of the distributions of the contents and stocks of C, N, P, and S along soil profiles with different water and salt conditions in the YRD; (2) investigate the deposition rates of C, N, P, and S in the soil using the 137 Cs dating technique in coastal wetlands; and (3) identify the changes in ecological stoichiometric ratios of C, N, P, and S with depth along the soil profiles.

Study Area
The study area is located in the YRD (N 37 • 35 -38 • 12 , E 118 • 33 -119 • 20 ), Dongying city, Shandong Province, China. It has a warm temperate semihumid monsoon climate with four distinct seasons, rain and heat at the same time, small regional climate differences, and a frost-free period of 196 days (Cui et al., 2008). The annual average sunshine hours are 2,590-2,830 h, and the  annual average temperature is 12.1 • C. The annual precipitation is 551.6 mm, 70% of which is concentrated in the period from June to August. The annual evaporation is 1,962 mm. The spring evaporation is strong, and the evaporation accounts for approximately 51.7% of the year. The drought index is as high as 3.56 . The main vegetation types are herbs, such as Pteris violata, Phragmites australis, Triarrhena sacchariflora, Suaeda salsa, Myriophyllum spicatum, and Limoninum sinense; shrubs, such as Tamarix chinensis; and trees, such as Salix matsudana (Cheng et al., 2021). There are more than 40 families, 110 genera, and 160 species. Phragmites australis, Pteris violate, and Suaeda salsa are the main dominant species in this area.
Three sampling sites were selected along one sampling belt perpendicular to the riverbed in the coastal wetlands on the north bank of the Yellow River (Figure 1), including Pteris violata wetlands (PV), Phragmites australis wetlands (PA), and Suaeda salsa wetlands (SS). Among them, PA wetlands can be affected by underground seawater and freshwater when the water and sediment regulation project of the YRD is implemented in July. Comparatively, SS wetlands are dominantly affected by tidal water.

Sample Collection and Analysis
The soil samples were collected by excavating the soil profiles at a depth of 0-100 cm at intervals of 2 cm. After removing visible plant residues and stones, the soil samples were airdried for 2 or 3 weeks. One replicate was sampled at each point for analyzing relevant indicators. One part of the air-dried samples were ground through sieving using a 20-mesh for the determination of 137 Cs labeling. Another part of the sample was further ground to pass through a 100-mesh sieve for determining soil properties. SOC was measured by the potassium dichromate dilution thermal-colorimetric method. Total nitrogen (N) was measured on an elemental analyzer (CHNOS elemental analyzer, Vario EL, Germany). Total phosphorus (P) and total sulfur (S) were determined by inductively coupled plasma atomic emission spectroscopy (ICP/AES) after digestion in an HClO 4 -HNO 3 -HF mixture in Teflon tubes. Quality assurance and quality control were assessed using duplicates, method blanks, and standard reference materials (GBW07401) from the Chinese Academy of Measurement Sciences with each batch of samples (1 blank and 1 standard for every 10 samples). The recoveries of samples spiked with standards ranging from 95 to 106%. 137 Cs value in the soil was determined using a high-purity germanium gamma spectrometer in the laboratory of the Geography Department, Beijing Normal University.

Stocks and Topsoil Concentration Factors
The stocks of C, N, P, and S were calculated as follows: where CS, NS, PS, and SS are the stocks of C, N, P, and S per unit area (kg/m 2 ), respectively; C i , N i , P i , and S i are the average contents of C, N, P and S in the soil layer i (i = 1, 2, 3, 4, and 5); and T is the thickness of soil layer i (cm). The topsoil concentration factor is calculated by the stocks in the top 10 cm and the stocks in the top 100 cm.

Statistical Analysis and Graphing
One-way ANOVA was used to identify the significant difference in C, N, P, and S contents and stocks between different wetlands. The difference was considered to be significant if P < 0.05. Data processing was performed using Excel 2019 software package software, and graphs were created by Origin 2017 software package.

Contents and Stocks of C, N, P and S
The highest average contents of C, N, and S in the PA wetland were 3.91, 0.24, and 0.24 g/kg, respectively, and the PV wetland had the highest average contents of P (0.65 g/kg) (Figure 2). SOC content in PA wetlands was significantly higher than that in PV wetlands (P < 0.05). Three layers with high concentrations of C, N, and S in soil profiles were observed in PA and SS wetlands at soil depths of 0-10 cm, 40-60 cm, and 90-100 cm. In PV wetlands, higher concentrations of C (5.29 g/kg) and S (0.36 g/kg) appeared in surface soils (0-10 cm) and N (0.48 g/kg) at soil depths of 90-100 cm, respectively. However, P contents showed different depth distributions in three wetlands, with an accumulative peak at the soil depth of 4-6 cm (3.70 g/kg) in PV wetlands, 10-12 cm (2.58 g/kg) in SS wetlands, and 80-82 cm (3.36 g/kg) in PA wetlands. SOC stocks in 1-m depth reached approximately 6 kg/m 2 in PA wetlands, slightly higher than that in the PV (3.55 kg/m 2 ) and SS (5 kg/m 2 ) wetlands (Figure 3). In contrast, low soil N stocks (0.36 kg/m 2 , PA; 0.22 kg/m 2 , SS; and 0.32 kg/m 2 , PV) were observed in the YRD. Unlike C and N, three wetlands showed similar soil P stocks (0.92 kg/m 2 ) in 1-m depth. Similar to N, PA wetlands exhibited the largest soil S stocks (0.36 kg/m 2 ), followed by SS wetlands (0.31 kg/m 2 ), while the lowest soil S stocks (0.23 kg/m 2 ) in PV wetlands were observed. SOC stocks in three accumulative hotspots in PA wetlands were 97.79 g/m 2 in top 10 cm soils, 148.25 g/m 2 for the 40-60 cm soil depth, and 66.32 g/m 2 for the 90-100 cm soil depth, also 64.67, 113.86, and 59.25 g/m 2 in SS wetlands, respectively. Soil TN stocks in three accumulative hotspots were 4.09 g/m 2 (0-10 cm), 12.14 g/m 2 (40-60 cm), and 5.36 g/m 2 (90-100 cm) in PA wetlands, with TS stocks of 3.87, 11.50, and 3.65 g/m 2 , respectively. In SS wetlands, TN and TS stocks were 2.05 and 5.03 g/m 2 in the top 10 cm, 5.88 and 2.81 g/m 2 for the 40-60 cm soil depth, and 3.52 and 1.99 g/m 2 for the 90-100 cm soil depth. Soil P stocks in the accumulative hotspot were 9.85 g/m 2 for the 10-20 cm soil depth in SS wetland and 11.25 g/m 2 for the 80-90 cm soil depth in the PA wetland. In contrast, the accumulative hotspot of C, P, and S stocks appeared in top 10 cm soils of PV wetlands (with the levels of 74.06, 11.67, and 3.83 g/m 2 , respectively), whereas the accumulative hotspot of N appeared in the 90-100 cm soil depth (6.96 g/m 2 ).

Soil Concentration Factors and Deposition Rates
The top concentration factors of C, N, P, and S in each wetland were all more than 0.10, which indicated that the accumulation of these elements was mainly affected by plant secretions (Table 1). Among them, the topsoil concentration factor of C fell within the scope of 0.16-0.22, following the order PV > PA > SS. As for soil N, their topsoil concentration factors in the three wetlands were relatively low. The followed order of topsoil concentration factors of both P and S was PV > SS > PA. The topsoil concentration factors of C, N, P, and S in the PA wetlands were higher than that of PV and SS wetlands. After 137 Cs isotope dating, the depth of the soil layer represented in the year 1963 was identified at the 20cm soil depth in PV wetland, 56-cm soil depths in PA, and 58-cm soil depths in SS wetlands. As shown in Figure 4, the deposition rate of sediments in PV wetland was 0.38 cm/yr, 1.08 cm/yr in PA wetland, and 1.06 cm/yr in SS wetland.
Deposition Fluxes of C, N, P and S Different deposition rates of C, N, P, and S in three wetlands were also observed. PA wetlands exhibited the highest deposition rate of C (71.17 g/yr/m 2 ), N (4.35 g/yr/m 2 ), and S (4.34 g/yr/m 2 ), while the lowest deposition rate appeared in PV wetlands ( Figure 5). In contrast, the highest deposition rate of P followed the order SS wetlands (9.82 g/yr/m 2 ) > PA wetlands (9.70 g/yr/m 2 ) > PV wetlands (4.08 g/yr/m 2 ). The contents of four elements in three wetlands showed a similar pattern with depth along soil profiles (Figure 2). Except for N, the contents of C, P, and S in the top 15-cm soils were relatively high, then declined significantly at approximately 20-cm soil depth and kept stable. However, at soil depth of approximately 40 cm, the contents of C, N, and S in the soil increased significantly again, but the change of P content was not obvious at this depth. At approximately 60-cm soil depth, the contents of these four elements again fell to a relatively stable level, to the soil layer of approximately 90-100 cm.

Ecological Stoichiometric Ratios
The soil C/N ratios in the three wetlands remained relatively stable (Figure 4). PV wetlands showed a higher range of the C/N ratios than PA and SS wetlands. This indicated that the C/N ratios at higher ranges in low-salinity wetlands were lower than that in middle-salinity and high-salinity wetlands. The highest C/P ratios were observed in PA wetlands, followed by those in SS wetlands, while the lowest value appeared in PV wetlands. Similarly, both C/N and C/P ratios were higher in surface soils of the three wetlands, and higher ratios also appeared at soil depths of 40-60 and 90-100 cm in PA and SS wetlands. However, this was not the case in PV wetlands. PA wetlands showed higher C/S ratios in upper soils than deeper soils, while a fluctuating change was observed along soil profiles in PV and SS wetlands. The C/N/P of SS wetlands was significantly higher than that of PA wetlands; the C/N/P of PA wetlands was significantly higher than that of PV wetlands (P < 0.05). N/P ratios decreased along soil profiles in PV wetlands, and a fluctuation was observed in PA and SS wetlands with higher ratios at the soil depths of 40-60 and 90-100 cm.

Depth Distributions of the Contents of C, N, P and S Along Soil Profiles
Higher contents of C, N, P, and S in surface soils (0-10 cm) of the three wetlands were generally observed, which was associated with exogenous inputs of plant litters, sedimentation by freshwater, and seawater input (Saintilan et al., 2013). The contents of C and N in the surface soils were higher than that in the deep soils in coastal wetlands . The Frontiers in Ecology and Evolution | www.frontiersin.org reason might be that the surface soils are more easily affected by environmental changes (Salome et al., 2010). In addition, higher contents of C, N, and S in the three wetlands were also observed at soil depths of 40-60 and 90-100 cm, which was associated with historical input and sedimentation. Over the past years, historical New wetlands on the east coast of the U.S.
1.58 ± 0.11 Autumn Krull and Craft, 2009 Mature wetlands on the east coast of the U.S.
5.83 ± 0.25 Autumn Krull and Craft, 2009 plant death and sedimentation due to flooding or Yellow River runoff might contribute to the elemental accumulation in deep soils. According to the annual Yellow River Yearbook, the Yellow River was diverted in 1964, which may have caused a large number of plants to be buried underground at this stage and converted into biogenic elements, such as C, N, P, and S. Another reason might be related to the high content of silt and clay in these layers since soil C content has a linear relationship with the composition content of clay + silt in the YRD (Zhao et al., 2020). Additionally, microbial activity and leaching could also cause C, N, and S in the surface soil to migrate downward (Wang et al., 2010).

Soil Deposition and Stocks of C, N, P and S in Soil Profiles
Compared with PA and SS wetlands, PV wetlands showed lower soil deposition at the rate of 0.38 cm/yr, which could be ascribed to the strong erosion by river water in PV wetland and higher sediment deposition in PA (1.08 cm/yr) and SS (1.06 cm/yr) wetlands affected by the flow-sediment regulation in July and sediment input by tidal flow . Our results fell within the range by DeLaune et al. (2003) that the deposition rate of Louisiana estuarine wetlands affected by freshwater input was 0.10-1.11 cm/yr. Ding et al. (2016) also reported that the deposition rate of tidal flats in the YRD was 0.58 cm/yr, which was lower than that of PA and SS wetlands and higher than that of PV wetland observed in this study. This indicated that tidal flats underwent much stronger hydraulic erosion by tidal flow than coastal salt marshes (Ding et al., 2016). Moreover, river flow might have a much stronger erosion than tidal flow in the study area.
In general, the stocks of C, N, P, and S in PA wetlands were higher than those of PV and SS wetlands. This might be associated with less erosion by runoff or tidal flow in PA wetlands. Moreover, higher plant biomass of Phragmites australis and high litter inputs in PA wetlands than PV and SS wetlands could explain higher levels of these elements (Lu et al., 2018). Additionally, higher salinity in SS wetland may inhibit the accumulation of SOC (Zhao et al., 2017). The SOC stocks in this study were 1.47-1.89 kg/m 2 at the top 30 cm soils, which was within the range of the previous results (1.17-2.14 kg/m 2 ) (Yu et al., 2013). Compared with other Chinese coastal wetlands, the SOC storage in the YRD was lower than that of the Chongming Dongtan coastal wetland (2.32 kg/m 2 ) (Jiang et al., 2015), but higher than that in the Shanghai tidal flats (1.38 kg/m 2 ) (Shi et al., 2010; Table 2). The possible explanation was that the YRD is a newly formed wetland, vegetation types and different runoff or tidal flow conditions would also lead to different SOC storage in different regions. The organic carbon storage in the YRD is close to that of the new wetlands, but much lower than the mature wetlands on the eastern coast of the United States (Krull and Craft, 2009). This was because the vegetation of mature wetlands had gradually evolved into more advanced large plants, and more SOC was imported into the soil. This also indicated that wetland protection and restoration will help to improve the YRD's carbon storage.
The SOC content of the three accumulative hotspots of PA wetlands accounted for 16.69% (0-10 cm), 25.30% (40-60 cm), 11.32% (90-100 cm) of the total soil profile; and 13.97% (0-10 cm), 24.60% (40-60 cm), 12.97% (90-100 cm) in the SS wetlands. The proportion of SOC in the three accumulative hotspots of PA and SS wetlands exceeded 40% (PA: 53.31%, SS: 51.36%). SOC stocks of PV wetlands in the top 10 cm soils accounted for 19.1% of the soil profiles. The SOC storage in the 0-10 cm hot zone showed PV > PA > SS, which may be related to the difference of soil salinity. In the deep hot zone (90-100 cm), it was just the opposite (PV < PA < SS). The results showed that the increase in salinity might accelerate the downward migration of SOC.
The proportion of N in the three accumulative hotspots of PA and SS wetlands exceeded 40% (PA: 59.6%, SS: 52.8%). The proportion of S in the three accumulative hotspots of PA and SS wetlands exceeded 40% (PA: 52.8%, SS: 52.0%). The distribution of P is relatively uniform in the entire section, and only a significant jump occurs in some soil layers. Similar to the Chongming Dongtan wetlands (Jiang et al., 2015), the nitrogen content of the YRD also showed a high peak at a depth of about 10 cm. In addition, according to the C/N comparison with other wetlands, the YRD is a nitrogen-restricted area (Cheng et al., 2021). Compared with PV wetlands, the stocks of C and S in both PA and SS wetlands increased significantly at soil depths below 40 cm. The possible reason was higher historical inputs and underground seawater inputs. Additionally, the root input of these wetland plants is mainly concentrated at a depth above 40 cm. The assimilation and utilization of nutrients by the root system make the soil at 0-40 cm less neutralized by C and S .
Top concentration factors for three wetlands showed higher stocks of C, N, P, and S in surface soils (0-10 cm), indicating that plant-soil cycling is the dominant factor influencing the biogeochemical processes of nutrients in the YRD. This was in agreement with the results of Bai et al. (2016). The topsoil concentration factor of C, P, and S in PV wetlands was higher than that of PA and SS wetlands. This may be associated with the effects of freshwater and seawater in PV and SS wetlands, respectively. Compared with SS wetland impacted by tidal water, the salinity of SS wetland was higher than that of PV due to the impacts of tidal flooding in SS wetlands with worse plant growth. Appropriate water conditions and lower salinity would have a positive effect on the topsoil concentration factor since water and salinity conditions can greatly influence vegetation covers and plant growth. Therefore, SS wetlands had a lower surface enrichment factor than PV wetlands. Additionally, higher salinity in SS wetlands might enhance the leaching effect of soil nutrients .

Changes in Ecological Stoichiometric Ratios Along Soil Profiles
The C/N ratio of the soil has been proven to be an important indicator for the ability of microorganisms to decompose soil organic matter (Canfora et al., 2017). If C/N < 9 it represents the stage of accelerating the decomposition of soil organic matter; C/N ratio between 9 and 11 indicates a dynamic equilibrium stage; and C/N > 11 can indicate the process of complete humification (Thomsen et al., 2008). Parolari and Porporato (2016) also proposed that when C/N > 10, the mineralization of SOC began to be restricted. Generally, the C/N ratios in this study were within the range of 10-50, except for a few low-value points. This indicated that soil organic matter is in the process of complete humification in the YRD, and the mineralization of organic carbon in this area had been inhibited to a certain extent. Cleveland and Liptzin (2007) observed that the soil C/N ratio was relatively consistent across various ecosystem types, although the soil was characterized by high biological diversity, structural complexity, and spatial heterogeneity (Yang et al., 2010). Our results in the current study showed that the C/N ratio in the surface soils was higher than that in deeper soils in PV and SS wetlands (P < 0.05). Although there was no statistically significant value, the soil C/N ratios tended to decrease with soil depth in PV and SS wetlands. However, the soil C:N ratios were relatively unstable in the PA wetland.
The C/P ratio is usually an indicator for the mineralization ability of soil organic phosphorus, and C/S represents the effect of the microbial biomass of the site on the availability of soil S (Heinze et al., 2010). When the ratios of C/P and C/S in a certain area are high, the limits of P and S during the decomposition of soil organic matter are not conducive for the growth of plants in the area. In contrast, low C/P and C/S ratios will contribute to the release of nutrients during the decomposition of organic matter and increase the available P and S levels in the soil. In addition, Reddy and DeLaune (2008) presented that the initial S mineralization would occur when the C/S ratio is less than 200, and the initial S immobilization will occur when the C/S ratio is in access of 400. In this study, the low ratios of C/P (2-7) and C/S (10-35) indicated that the mineralization of P and S, instead of immobilization, is the main process in the YRD (Figure 4), resulting in lower accumulation in this region. Compared with other coastal wetlands in China, the soil nutrients in the YRD are relatively low (Zhang et al., 2012(Zhang et al., , 2013, indicating that the exchange of carbon and nutrients between the YRD and the external environment might be very active (Cao et al., 2015). Soil N/P and C/N/P also could be used as indicators of nutrient restriction types. In the current study, The N/P ratios in three wetlands were less than 1.2, which was similar to the results of Gao et al. (2012). Similarly, the C/N/P ratios were between 0 and 0.1. Therefore, the YRD is a nitrogen-limited area, and the accumulated SOC content in this area was relatively less than N and P.

CONCLUSION
The deposition and ecological stoichiometric of C, N, P, and S were investigated in coastal wetlands along a sampling belt of the YRD in this study. The results showed three accumulative hotspots of C, N, and S in soil profiles in PA and SS wetlands (0-10, 40-60, and 90-100 cm), and one accumulative hotspot of P in SS (10-20 cm) and PA (80-90 cm) wetlands. In PV wetlands, one accumulative hotspot of C, P, and S appeared in the top 10 cm, whereas N at the soil depth of 90-100 cm. Generally, due to the erosion of rivers, higher soil deposition rates and deposition fluxes of C, N, P, and S were observed in PA and SS wetlands compared with PV wetlands. Topsoil concentration factors showed that the area was mainly affected by the input of plant residues, while the leaching effect gradually strengthens with increasing salinities. The ecological stoichiometry further verified that coastal wetlands in the YRD were nitrogen-limited ecosystems and in the stage of accelerating the decomposition of soil organic matter. Improving the degree of salinization in this area, increasing the vegetation coverage of YRD, and reducing the erosion effect of rivers on soil will not only help to improve the ecological stability of this area but also increase the blue carbon storage of YRD and make better use of YRD's carbon sink function.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
SD and CW contributed to the writing of the manuscript. SD, QZ, and JJ involved in the field work and data collection. JB and QZ contributed to concept of study. JB, CW, YG, JJ, GZ, and CY provided important guidance on methods and writing. All authors contributed critically to drafts and gave approval for publication.

FUNDING
This study was financially supported by the Joint Funds of the National Natural Science Foundation of China (No. U2006215) and the National Natural Science Foundation of China (No. 42107490).