Increase in Soil Carbon Pool Stability Rather Than Its Stock in Coastal Saline—Alkali Ditches following Reclamation Time

: Extensive drainage ditches are constructed to reduce soil salinity in reclaimed saline– alkali farmland, consequently forming plant growth hotspots and impacting soil carbon stocks therein. However, the investigation into changes in soil carbon stocks remains limited in these ditches. To address this, soil samples were collected from drainage ditches, which originated from the reclamation of saline–alkali farmland, at different reclamation years (the ﬁrst, seventh, and ﬁfteenth year). Moreover, fractions were separated from soil samples; a particle size separation method (particulate organic matter, POM; mineral–associated organic matter, MAOM) and a spatio–temporal substitution method were conducted to analyze the variations in soil carbon components and the underlying mechanisms. The results indicate that there were no signiﬁcant variations in the contents and stocks of soil organic carbon (SOC) and soil inorganic carbon (SIC) following the increase in reclamation time. However, in the POM fraction, the SOC content (SOC POM ) and stock signiﬁcantly decreased from 2.24 to 1.12 g kg − 1 and from 19.02 to 12.71 Mg ha − 1 , respectively. Conversely, in the MAOM fraction, the SOC content (SOC MAOM ) and stock signiﬁcantly increased from 0.65 to 1.70 g kg − 1 and from 5.30 to 12.27 Mg ha − 1 , respectively. The different changes in SOC POM and SOC MAOM , as well as the result of the structural equation model, showed a possible transformation process from SOC POM to SOC MAOM in the soil carbon pool under the driving force of reclamation time. The results in terms of the changes in soil carbon components demonstrate the stability rather than the stock of the soil carbon pool increase in coastal saline–alkali ditches following the excavation formation time. Although more long time series and direct evidence are needed, our ﬁndings further provide a case study for new knowledge about changes in the soil carbon pool within saline–alkali ditches and reveal the potential processes involved in the transformation of soil carbon components.


Introduction
The soil carbon pool is the largest carbon pool in terrestrial ecosystems, composed primarily of soil organic carbon (SOC) and soil inorganic carbon (SIC) [1,2].Changes in its stock have significant implications for global carbon cycling and climate change [3,4].The soil carbon pool is influenced by natural factors such as vegetation, soil quality, temperature, and precipitation [1,5].However, the significant impact of human activities, such as land management and land use changes, on the soil carbon pool should not be underestimated [4,6], particularly in the land use change after reclamation.For example, the irrational process of land reclamation disrupts soil aggregations [7], thereby affecting plant root growth and microbial activity [8], ultimately leading to a significant decrease in SOC and soil nutrient content [9].Therefore, investigations into the changes in soil carbon pools are crucial for formulating rational land management strategies, as well as for regulating the soil carbon stocks.
Divergence variations in soil carbon pools after land reclamation have been exhibited in the recent studies.It is characterized by significant increases [10,11], significant decreases [12], and decreases, followed by increases in SOC content [13].Furthermore, SIC content slightly increases with shorter reclamation times and decreases with longer reclamation times [11].Similar trends are observed in carbon stocks after reclamation [14,15].These inconsistencies are primarily due to the variations in soil types and vegetation types among the various experimental conditions, as well as the complex interactions between SOC and SIC.Generally, the stock of the soil carbon pool relies on a balance between the formation of SOC from plant residues and the mineralization of SOC [16].Plant residues and root exudates are crucial sources of SOC, thereby rendering the changes in the soil carbon pool easily influenced by vegetation conditions after reclamation, e.g., plant biomass and diversity after reclamation [17,18].Previous research has focused on driving factors impacting the soil carbon pool, such as reclamation time and land use type [13,19].However, investigating the changes in different soil carbon components and the mechanisms of SOC to SIC pool transformation remains limited, which hinders a deeper understanding of the processes and mechanisms of soil carbon stock changes therein.
Soil organic matter (SOM) is directly related to SOC and influences the stability of the changes in soil carbon pools [20,21].Due to the complexity of SOM, characterizing the different fractions of SOM has become a conventional technique for revealing its transformation processes [22][23][24].SOM is commonly categorized into two fractions: particulate organic matter (POM) and mineral-associated organic matter (MAOM), based on its size or density [25][26][27][28].POM primarily originates from plant residues [29,30], while MAOM is mainly derived from microbial metabolites [31,32].However, recent research suggests that soluble organic compounds also directly contribute to the formation of MAOM [16,33].Some studies suggest that the decomposition of POM generates soluble organic compounds which directly contribute to MAOM formation, thus linking POM and MAOM [28].Other studies propose that the formation of MAOM is independent of POM and relies on the input of soluble compounds from plant residues [34,35].Due to the complex relationship between POM and MAOM, previous investigation has attempted to quantify the association between POM and MAOM fractions with SOC or soil total carbon (TC) [27].However, few studies have examined the relationship between SIC and POM or MAOM fractions, as well as their potential transformations.Therefore, it is necessary to fully consider the changes in SOC and SIC of POM and MAOM fractions and investigate their transformation mechanisms.
The Yellow River Delta is formed by the sediment of the Yellow River.It is primarily comprised of coastal saline-alkali soil as a result of intense land-sea interactions [36].In recent years, the construction of numerous drainage ditches has been carried out to reduce soil salinity during land reclamation, with the unutilized land accounting for 9.10% or more of the total farmland area [37].The specific conditions of drainage ditches can impact soil carbon stocks.For example, waterlogging resulting from drainage alters the structure of soil microbial communities, inhibits microbial activity, and reduces the rate of SOM decomposition and mineralization [11,[38][39][40][41].Moreover, the accumulation of soil salinity in drainage ditches may restrict plant growth and microbial activities, thereby affecting changes in SOC stock and stability [42].Conversely, the unique structure of drainage ditches promotes material exchange and creates an improved saline-alkali environment, leading to the formation of plant growth hotspots.Considering the changes in the amount of carbon input by plants into the soil and the rate of SOM mineralization influenced by microbial activity, the presence of ditches can significantly influence the original soil carbon stocks and stability.The extent of salt-affected land in the Yellow River Delta region is approximately 240,000 ha [43].If all of them were reclaimed as saline-alkali farmland with drainage ditches, the changes in soil carbon stock should not be underestimated.
Therefore, the objective of this study is to investigate the processes and transformation mechanisms of changes in soil carbon stocks within saline-alkali ditches formed by farmland reclamation in the Yellow River Delta.We measured the carbon content and physicochemical properties of the drainage ditch soils in saline-alkali farmland after different reclamation times (the first, seventh, and fifteenth year).Furthermore, we attempted to explore the processes and mechanisms of changes in soil carbon components by separating POM and MAOM fractions.The present study provides a theoretical foundation for managing the soil carbon pool in recovered saline-alkali areas.carbon stocks and stability.The extent of salt-affected land in the Yellow River Delta region is approximately 240,000 ha [43].If all of them were reclaimed as saline-alkali farmland with drainage ditches, the changes in soil carbon stock should not be underestimated.Therefore, the objective of this study is to investigate the processes and transformation mechanisms of changes in soil carbon stocks within saline-alkali ditches formed by farmland reclamation in the Yellow River Delta.We measured the carbon content and physicochemical properties of the drainage ditch soils in saline-alkali farmland after different reclamation times (the first, seventh, and fifteenth year).Furthermore, we attempted to explore the processes and mechanisms of changes in soil carbon components by separating POM and MAOM fractions.The present study provides a theoretical foundation for managing the soil carbon pool in recovered saline-alkali areas.

Site Description
This research was conducted at the Yellow River Delta Research Center, Chinese Academy of Sciences, in Dongying, Shandong Province, China (37°40′N, 118°55′ E).The study area is situated in the alluvial deposition zone of the Yellow River Delta, characterized by saline-alkali soils.The area experiences a temperate continental monsoon climate, with an average annual temperature of 12.40 °C.The annual average precipitation ranges from 530-630 mm, with the majority occurring between June and September.Additionally, the annual evaporation is 1962 mm.Over recent years, extensive land development has occurred, resulting in the creation of numerous drainage ditches (Figure 1).The soil in this region has originated from the sediment of the Yellow River since 1855, and the main soil type is Calcaric Fluvisols, comprising 5.90% caly, 70.88% silt, and 23.22% sand.The vegetation in these drainage ditches primarily comprises reed (Phragmites australis (Cav.)Trin.ex Steud.),Suaeda salsa (Suaeda salsa (L.) Pall.), Imperata cylindrica (Imperata cylindrica (L.) P. Beauv.), Leymus chinensis (Leymus chinensis (Trin.), and Artemisia capillaris (Artemisia capillaris Thunb.).

Soil Sampling Collection
To investigate the changes in soil carbon stocks of saline-alkali ditches, three 1 km × 1 km sampling sites with different reclamation years, including the 1st, 7th, and 15th year, were selected in the saline-alkali farmland within a distance of 2 km from each other.Within each sampling site, four drainage ditches adjacent the farmlands were selected as sampling plots which had similar cross-section dimensions, ditch length, and vegetation coverage.Furthermore, within each sampling plot, three slope cross-sections were randomly selected.In each cross-section, five 1 m × 1 m subplots on the different positions, including the upper slope, middle slope, lower slope, slope crest, and furrow sole, were selected to collect the samples (Figure 1).For each subplot, we measured the height of the plants and assessed the vegetation coverage.Subsequently, all aboveground plants in the subplot were collected to evaluate the biomass.The vegetation coverage and growth conditions for sampling plots with the same reclamation year are presented in Table S1.Soil samples at depths of 0-20 cm, 20-40 cm, and 40-60 cm were collected using a soil auger.Due to the complex water and soil conditions in the furrow sole, soil samples were taken at depths of 0-20 cm and 20-40 cm.Soil samples obtained from the corresponding position and depth in the three slope sections of each plot were mixed to create a representative sample.In total, 168 composite samples were obtained.The mixed soil sample were sieved using a 2 mm mesh to remove plant debris and divided into two portions.The first portion, 500 g, was air-dried and used to measure soil parameters, including TC, SIC, SOC, POM, MAOM, soil total nitrogen (TN), soil total phosphorus (TP), soil available phosphorus (AP), soil pH, and soil electrical conductivity (EC).Moreover, to obtain the fractions of POM and MAOM, a particle size separation method was employed [25,44].Specifically, 10 g of the air-dried soil was dispersed in 30 mL of 5 g L −1 sodium hexametaphosphate via shaking for 15 h on a reciprocal shaker.The dispersed soil samples were passed through a 53 µm sieve and, after rinsing several times with water, the material that was retained on the sieve was defined as POM fraction (53-2000 µm) and the material passing through the sieve (<53 µm) was labeled as MAOM fraction.The divided soil was further utilized to determine the content of organic carbon (SOC POM ) and inorganic carbon (SIC POM ) in POM fraction soil, as well as organic carbon (SOC MAOM ) and inorganic carbon (SIC MAOM ) in MAOM fraction soil.The second portion, weighing 300 g, was stored at 4 • C and used to measure soil water content (SWC), soil dissolved organic carbon (DOC), soil soluble ammonium nitrogen (NH + 4 -N), and soil soluble nitrate nitrogen (NO − 3 -N).

Soil Sample Processing
TC and SIC contents were analyzed using an automated carbon-nitrogen analyzer (Primacs SNC100-IC-E, Skalar Analytical BV, Breda, The Netherlands).TC content was determined based on high-temperature combustion and non-dispersive infrared detection (NDIR), while SIC content was obtained via automatic acidification heating, gas purging, and non-dispersive infrared detection (NDIR).SOC content was calculated by subtracting SIC content from TC content.AP content was measured using a continuous flow analyzer (Auto Analyzer 3, Bran and Luebbe GmbH, Norderstedt, Germany).Soil pH and EC values were determined by a soil-water mixture with a ratio of 1:5.Soil pH was measured using a pH meter (PHS-3E, Leici Instruments, Shanghai, China) and EC was measured using an electrical conductivity meter (FE38-Standar, Mettler Toledo, Greifensee, Switzerland).TN content was determined using the Kjeldahl method, while TP content was determined using the H 2 SO 4 -HClO 4 digestion method.Both TN and TP were measured using the continuous flow analyzer.DOC content was determined using the method developed by Zhan and Zhou [45].NH + 4 -N and NO − 3 -N contents were determined using the KCl-indophenol blue colorimetric method and the dual-wavelength UV-Visible spectrophotometric method, respectively.

Calculation of Soil Carbon Stocks
The SOC and SIC stocks at different slope sampling positions are calculated according to Liu et al. [46] as follows: where M SOC is the SOC stocks (Mg ha −1 ); M SIC is the SIC stocks (Mg ha −1 ); D is soil depth (cm); B is the soil bulk density (g cm −3 ); C SOC is the SOC content (g kg −1 ); C SIC is the SIC content (g kg −1 ); K is the proportional coefficient.Specifically, the proportional coefficients for soil at the slope crest, upper slope, middle slope, lower slope, and furrow sole positions are 0.26, 0.26, 0.24, 0.17, and 0.07, respectively, calculated based on the proportion of each part to the whole ditch.
The soil TC stocks at different slope sampling positions are as follows: where M TC is the TC stocks (Mg ha −1 ).
The organic carbon and inorganic carbon stocks of soil in the POM and MAOM fractions are obtained as follows: M PSIC = M SIC × P ( 5) where M PSOC is the organic carbon stocks in the POM fraction (Mg ha −1 ); M PSIC is the inorganic carbon stocks in the POM fraction (Mg ha −1 ); M MSOC is the organic carbon stocks in the MAOM fraction (Mg ha −1 ); M MSIC is the inorganic carbon stocks in the MAOM fraction (Mg ha −1 ); and P is the ratio of the POM fraction soil in the bulk soil.

Data Analysis
After testing the normality (Shapiro-Wilk's test) and homogeneity of variances (Levene's test) for all measured variables, the differences in soil properties (TC, SOC, SIC, DOC, SOC POM , SOC MAOM , SIC POM , and SIC MAOM ; TN, NH + 4 -N, NO − 3 -N, TP, AP, pH, EC, and SWC) among different reclamation years, soil depths, and slope sampling positions were examined using a one−way analysis of variance (ANOVA) with the least significant difference (LSD) test.The differences in soil carbon stocks (M TC , M SOC , M SIC , M PSOC , M MSOC , M PSIC , and M MSIC ) were also tested.Additionally, a three-way ANOVA with the least significant difference (LSD) test was used to determine the effects of reclamation years, soil depths, slope sampling locations, and their interactions on the determined soil properties.Linear regression analysis was conducted to determine the overall and grouped relationships (based on reclamation years and soil depth) between organic carbon and inorganic carbon contents in different soil fractions (SOC vs. SIC; SOC POM vs. SIC POM ; and SOC MAOM vs. SIC MAOM ).Pearson correlation analysis was performed to assess the relationships between soil carbon contents (TC, SOC, SIC, DOC, SOC POM , SOC MAOM , SIC POM , and SIC MAOM ) and soil properties (TN, NH + 4 -N, NO - 3 -N, TP, AP, pH, EC, and SWC).All results were considered statistically significant at p < 0.05 and the statistical analyses were conducted using SPSS (version 22.0.0.0,SPSS, Chicago, IL, USA).
After conducting Detrended Correspondence Analysis (DCA) using Canoco 5 (Version 5.0, Microcomputer Power, Ithaca, NY, USA), Redundancy Analysis (RDA) was performed to further determine the relationship between soil carbon content and other determined soil properties and to calculate the contribution of soil properties to soil carbon content.Structural Equation Modeling (SEM) was constructed using AMOS (Version 28, SPSS, IBM Corp, Armonk, NY, USA) to investigate the direct and indirect effects of reclamation years on SOC POM , SOC MAOM , SIC POM , SIC MAOM , DOC, and SIC in drainage ditch soils.The goodness of fit of the structural equation model was evaluated based on the chi-square test (χ2-text), degrees of freedom (DF), p-value, goodness-of-fit index (GFI), adjusted goodness-of-fit index (AGFI), comparative fit index (CFI), standardized root mean square residual (SRMR), and root mean square error of approximation (RMSEA).

Changes in Physicochemical Properties of Saline-Alkali Ditch Soil after Reclamation
The changes in the soil properties in saline-alkali ditches after reclamation are presented in Table 1.The TC content exhibited a significant decrease from 13.90 to 13.76 g kg −1 with the increasing reclamation years (p < 0.05).No significant differences were observed in SOC and SIC contents among different reclamation years (p > 0.05).However, the SOC content significantly decreased from 3.30 g kg −1 over 0-20 cm to 2.82 g kg −1 over 40-60 cm (p > 0.05), while the DOC content exhibited a significant decrease from about 48 to 36 mg kg −1 with the increasing reclamation years (p < 0.05).On the other hand, the SOC MAOM content showed a significant increase from 0.65 to 1.70 g kg −1 with increasing reclamation years (p < 0.05), whereas both the SOC POM and SIC POM content experienced significant decreases from 2.24 to 1.12 g kg −1 and from 10.54 to 9.54 g kg −1 , respectively (p < 0.05).The TN and AP content significantly decreased from 143.07 to 112.60 mg kg −1 and from 26.39 to 20.30 mg kg −1 with increasing reclamation years, respectively (p < 0.05), while the NO − 3 -N content significantly increased from 19.46 to 20.13 mg kg −1 and from 0.48 to 0.55 g kg −1 , respectively (p < 0.05).Additionally, SWC significantly increased from 26.05 to 29.88% with the increasing reclamation time and from 25.07 to 26.45% with the increasing soil depth (p < 0.05).Soil pH showed a significant increase from 8.51 to 9.21 with increasing reclamation years (p < 0.05); however, a decrease with no significant differences was found in the NH + 4 -N content and EC (p > 0.05) The three−way ANOVA conducted on soil properties in saline−alkali ditches (Table 2) revealed that reclamation years had the most significant impact on soil properties.However, the contents of SIC, SIC POM , and NO − 3 -N did not display any significant responses to changes in reclamation years, soil depth, slope sampling position, or their interactions (p > 0.05).Furthermore, soil depth had a significant impact on DOC content and SWC (p < 0.01), whereas soil sampling location had a significant impact on the SOC MAOM content and SWC (p < 0.05).

Changes in Soil Carbon Stocks in Saline-Alkali Ditches after Reclamation
The soil carbon stocks in saline-alkali drainage ditches varied, as shown in Figure 2. M TC , M SOC , and M SIC (Figure 2a-c), along with M PSIC and M MSIC (Figure 2f,g), did not show any significant differences across different reclamation years (p > 0.05).However, M PSOC significantly decreased from 19.02 to 12.71 Mg ha −1 , while M MSOC significantly increased from 5.30 to 12.27 Mg ha −1 in the 15th year compared to the 1st year (p < 0.05).Furthermore, among the different soil sampling positions in the saline-alkali ditches, M MSOC exhibited no significant differences, while M PSIC showed no significant differences except in the 15th year (p > 0.05).Notably, the soil carbon stocks of other components varied significantly among the lower slope, furrow sole, and other soil sampling positions.Generally, the carbon stocks in the lower slope and furrow sole positions were significantly lower than those in the other soil sampling positions (p < 0.05).

Changes in Soil Carbon Stocks in Saline-Alkali Ditches after Reclamation
The soil carbon stocks in saline-alkali drainage ditches varied, as shown in Figure 2. MTC, MSOC, and MSIC (Figure 2a-c), along with MPSIC and MMSIC (Figure 2f,g), did not show any significant differences across different reclamation years (p > 0.05).However, MPSOC significantly decreased from 19.02 to 12.71 Mg ha −1 , while MMSOC significantly increased from 5.30 to 12.27 Mg ha −1 in the 15th year compared to the 1st year (p < 0.05).Furthermore, among the different soil sampling positions in the saline-alkali ditches, MMSOC exhibited no significant differences, while MPSIC showed no significant differences except in the 15th year (p > 0.05).Notably, the soil carbon stocks of other components varied significantly among the lower slope, furrow sole, and other soil sampling positions.Generally, the carbon stocks in the lower slope and furrow sole positions were significantly lower than those in the other soil sampling positions (p < 0.05).

Correlation between Different Components of SOC and SIC
The linear regression analysis results (Figure 3) demonstrated the positive correlation between SOC and SIC overall (Figure 3a; R 2 = 0.33), as well as between SOCMAOM and SIC-MAOM (Figure 3b; R 2 = 0.21), and between SOCPOM and SICPOM (Figure 3d; R 2 = 0.18, p < 0.05).Only the relationship between SOCPOM and SICPOM in the 15th year showed a significant negative correlation (Figure 3e; R 2 = 0.04, p < 0.05) among the different reclamation year groups.The correlations between the organic carbon and inorganic carbon contents of other components were significantly positive (Figure 3d-f; p < 0.05).Similarly, a significant positive correlation between the organic carbon and inorganic carbon contents was found in each component among the different soil depth groups (Figure 3g-i; p < 0.05).

Correlation between Different Components of SOC and SIC
The linear regression analysis results (Figure 3) demonstrated the positive correlation between SOC and SIC overall (Figure 3a; R 2 = 0.33), as well as between SOC MAOM and SIC MAOM (Figure 3b; R 2 = 0.21), and between SOC POM and SIC POM (Figure 3d; R 2 = 0.18, p < 0.05).Only the relationship between SOC POM and SIC POM in the 15th year showed a significant negative correlation (Figure 3e; R 2 = 0.04, p < 0.05) among the different reclamation year groups.The correlations between the organic carbon and inorganic carbon contents of other components were significantly positive (Figure 3d-f; p < 0.05).Similarly, a significant positive correlation between the organic carbon and inorganic carbon contents was found in each component among the different soil depth groups (Figure 3g-i; p < 0.05).

Direct and Indirect Effects of Reclamation Time on Different Soil Carbon Components
The SEM results showed that the reclamation years indirectly promoted the increase in SIC content by promoting the increase in SOC MAOM content, thereby indirectly increasing SIC MAOM content (Figure 4).The reclamation years and SOC POM content promoted SIC content by suppressing the increase in SIC POM and SIC MAOM contents.The reclamation years inhibited SOC MAOM and SIC MAOM contents by suppressing DOC content.Additionally, DOC content indirectly increased SIC content by promoting SOC MAOM and SIC MAOM contents.

Direct and Indirect Effects of Reclamation Time on Different Soil Carbon Components
The SEM results showed that the reclamation years indirectly promoted the increase in SIC content by promoting the increase in SOCMAOM content, thereby indirectly increasing SICMAOM content (Figure 4).The reclamation years and SOCPOM content promoted SIC content by suppressing the increase in SICPOM and SICMAOM contents.The reclamation years inhibited SOCMAOM and SICMAOM contents by suppressing DOC content.Additionally, DOC content indirectly increased SIC content by promoting SOCMAOM and SICMAOM contents.

The Primary Factors Affecting the Soil Carbon Content of Different Components
The results of the Pearson correlation analysis showed that TN, AP, SWC, pH, and EC were the primary factors influencing the measured soil carbon content (Figure 5).Moreover, the RDA results demonstrated that the first two axes accounted for 21.74% and 7.31% of the variation in soil carbon content, respectively (Figure S1).These findings suggest that these soil physicochemical parameters collectively explained 33.10% of the total variation (Table S2).Additionally, the contribution of each soil physicochemical parameter in explaining the changes in soil carbon content of different components were ranked as follows: TN > AP > EC > SWC > pH (Table 3).
Based on a comprehensive analysis of the Pearson correlation and RDA results (Figures 5 and S1), several significant positive correlations were observed between TN and the TC, SOC, SIC, and SIC POM contents (p < 0.05).These significant positive correlations were found between AP and the SOC, SIC, SOC MAOM , and SIC MAOM contents (p < 0.05).Conversely, pH exhibited significant negative correlations with TC, SOC, and SIC POM contents (p < 0.05), while SWC showed significant negative correlations with SOC and DOC contents (p < 0.05).Notably, SWC and pH demonstrated significant positive correlations with SOC MAOM content (p < 0.05), whereas EC exhibited a significant negative correlation with the SOC MAOM and SIC MAOM contents (p < 0.05).Based on a comprehensive analysis of the Pearson correlation and RDA results (Figures 5 and S1), several significant positive correlations were observed between TN and the TC, SOC, SIC, and SICPOM contents (p < 0.05).These significant positive correlations were found between AP and the SOC, SIC, SOCMAOM, and SICMAOM contents (p < 0.05).Conversely, pH exhibited significant negative correlations with TC, SOC, and SICPOM contents (p < 0.05), while SWC showed significant negative correlations with SOC and DOC con-

Dynamics of the Soil Carbon Pool in Saline−alkali Ditches after Reclamation
The results demonstrate that the increasing reclamation time did not have a significant impact on the SOC and SIC contents of the saline-alkali ditches (Table 1; Figure 2b,c).This unchanging trend in SOC and SIC contents with respect to reclamation time differs from previous research on reclaimed saline-alkali farmland [11,47,48].In farmland, SOC tends to increase due to fertilization and the accumulation of organic residues, while SIC decreases because of irrigation and precipitation [11].The reasons for these discrepant results are attributed to two main factors.Firstly, reclaimed saline-alkali ditches are subject to less influence from agricultural activities compared to farmland, which allows for a greater input of carbon from plants into the soil; however, this input is still limited in comparison to the effects of fertilization and the presence of organic residues.Secondly, the structure of the ditches leads to higher groundwater levels and moisture content compared to farmland, potentially inhibiting the microbial mineralization of soil carbon [49].
Additionally, the reclamation years did not affect the storage capacity of the soil carbon pool but influenced the compositions therein (Figure 2).These results indicate a unique carbon fixation process in coastal saline-alkali ditches.Specifically, the findings demonstrate that as the reclamation time increased, a notable decrease was found in the levels of SOC POM and SIC POM contents, while SOC MAOM content increased (Table 1).Similarly, the M PSOC showed a decreasing trend, whereas the M MSOC increased (Figure 2d,e).This result indicates that the changes in different soil carbon components are inherent factors for maintaining constant soil carbon stock, and it may be due to a conversion from SOC POM to SOC MAOM during the reclamation.In addition, the significant decrease in SOC POM and increase in SOC MAOM at the furrow sole position indicate greater stability of soil carbon in this position.This is because POM has a shorter residence time and faster decomposition rate in the soil compared to MAOM [50][51][52].
The reclamation years have a significant impact on the different fractions of soil carbon pools in saline-alkali ditches (Table 2).In addition, the response of different soil carbon components to the reclamation years, constituting an inherent factor, affect the changes in soil carbon pool (Table S2; Figure 4).In light of the correlation between SOC and SIC contents in the soil fractions of POM and MAOM and the SEM results driven by the reclamation years (Figures 3 and 4), the process of changes in soil carbon components of saline-alkali ditches after reclamation may be as follows (Figure 6).Firstly, after the reclamation, the saline−alkali soil conditions gradually improve with an increase in reclamation time (Table 1).This improvement leads to the colonization and growth of numerous plants (Table S1), resulting in the entrance of plant residues into the soil.These organic materials further decompose into SOC POM , SOC MAOM , and DOC through microbial metabolisms.Secondly, due to the different turnover rates of POM and MAOM, SOC POM rapidly accumulates between the 1st and 7th year (Table 1), meanwhile, it gradually recomposes into SIC POM and SOC MAOM through continuous biotic and abiotic processes.
Thirdly, between the 7th and 15th year, SOC POM contents/stocks decrease due to the continuous microbial transformation and consummation (Figure 1).Some of the carbon therein is converted into SOC MAOM (Table 1; Figure 3c,f,i) and is decomposed to CO 2 and/or carbonates, which may accelerate the formation of insoluble compounds, e.g., the SIC POM (Figures 3e and 4).Fourthly, SOC MAOM significantly increases with reclamation due to the inputs of plant residues and the transformation of SOC POM (Figure 4).The formation of SIC MAOM predominantly occurs through the transformation of SOC MAOM and DOC (Figure 4).These results together suggest that the stability of the soil carbon pool has been changed after the reclamation.
reclamation, the saline−alkali soil conditions gradually improve with an increase in reclamation time (Table 1).This improvement leads to the colonization and growth of numerous plants (Table S1), resulting in the entrance of plant residues into the soil.These organic materials further decompose into SOCPOM, SOCMAOM, and DOC through microbial metabolisms.Secondly, due to the different turnover rates of POM and MAOM, SOCPOM rapidly accumulates between the 1st and 7th year (Table 1), meanwhile, it gradually recomposes into SICPOM and SOCMAOM through continuous biotic and abiotic processes.Thirdly, between the 7th and 15th year, SOCPOM contents/stocks decrease due to the continuous microbial transformation and consummation (Figure 1).Some of the carbon therein is converted into SOCMAOM (Table1; Figure 3c,f,i) and is decomposed to CO2 and/or carbonates, which may accelerate the formation of insoluble compounds, e.g., the SICPOM (Figures 3e  and 4).Fourthly, SOCMAOM significantly increases with reclamation due to the inputs of plant residues and the transformation of SOCPOM (Figure 4).The formation of SICMAOM predominantly occurs through the transformation of SOCMAOM and DOC (Figure 4).These results together suggest that the stability of the soil carbon pool has been changed after the reclamation.With an increase in reclamation years of saline−alkali ditches, the stability of the soil carbon pool improves (Table 1; Figure 2).This stability is largely influenced by the proportion of stable carbon in the pool.Previous studies have shown that the stability of POM primarily occurs via the encapsulation or formation of microaggregates (53-250 µm) and the stabilization of microbial residues on silt and clay (<53 µm) surfaces [27,52,53].Therefore, the increase in MAOM contributes to the stability of the soil carbon pool, which is consistent with our own findings (Table 1; Figure 2d,e).Additionally, the SOC POM exhibited a negative correlation with the SIC POM in the 15th year, in contrast to other reclamation times.This suggests that the POM decreased as the reclamation time increased, leading to an improvement in the stability of the soil carbon pool.Although DOC is the main and more effective pathway for the formation of SOC MAOM [16,33], a significant decrease in DOC with reclamation years also indicates an improvement in the stability of the soil carbon pool (Table 1).However, the SEM results demonstrate that DOC has a promoting effect on the formation of SOC MAOM with no significance (p > 0.05), while it significantly promotes the formation of SIC MAOM (Figure 4).This could be attributed to changes in plant biomass and an increase in soil pH resulting from reclamation activities (Tables 1 and S1).
In recent years, several models have been proposed regarding the formation mechanisms of POM and MAOM [54,55].In these models, POM is expected to form from partially decomposed plant residues once they enter the soil.As a result, it is believed that POM does not transform as MAOM after decomposition [16,56].However, other studies suggest that during the process of POM decomposition, lower molecular weight soluble organic compounds are formed, such as DOC [33].Additionally, POM is mineralized by microbial metabolisms or bound to minerals to form MAOM [57,58].These different perceptions highlight the complexity of the POM to MAOM transformation process and suggest that it may not solely depend on the conversion of soluble organic compounds.Though the application of isotopic and chemical analysis is still necessary to confirm our research findings, our results indicate a possible transformation from SOC POM to SOC MAOM between POM and MAOM.Nonetheless, this suggests that in soil environments where DOC is easily removed, such as in frequently drained ditch environments, the transformation from POM to MAOM may be more likely to occur.Our results further indicate a unique carbon fixation process in coastal saline-alkali ditches, that is, although the storage has not increased following the reclamation years, the stability has increased.

The Primary Factors Affecting Soil Carbon Pools in Saline-alkali Ditches
For the saline-alkali ditches, TN and AP were identified as the predominant factors contributing to promoting the increase in different soil carbon components (Table 3; Figures 5 and Figure S1).This is mainly attributed to the fact that nitrogen and phosphorus are essential elements for plant growth [59]; plant growth may facilitate the increase in the soil carbon pool.However, no increase in soil carbon stocks in saline-alkali ditch is likely due to the opposite variations in TN and AP contents following reclamation (Table 1).In the early stage of saline-alkali farmland development, a large amount of fertilizer will be applied to improve soil nutrients.In view of the different migration processes of the elements of N and P, irrigation and salt drainage process can lead to the migration of N to the ditch, resulting in the phenomenon of total nitrogen increasing first and then decreasing following reclamation of the ditch [60].
Although EC, SWC, and pH inhibited the increase in different soil carbon components, their combined effect increased the soil carbon pool stability (Tables 1 and 3; Figures 5 and S1).Our results showed no significant decrease in EC, whereas significant decreases in pH and SWC were shown (Table 1).Moreover, the significant negative correlation between EC and SOC MAOM , and the significant positive correlation between SWC, pH, and SOC MAOM , suggest that EC, SWC, and pH all contribute to the accumulation of SOC MAOM (Figures 5 and S1).The changes in EC, SWC, and pH possibly inhibit microbial activity, making SOC MAOM easier to accumulate, which contributes to increasing the soil carbon pool stability.Additionally, SWC and pH were significantly negatively correlated with SOC and DOC contents, indicating that the increase in SWC and pH after reclamation inhibited the production of SOC and DOC (Table 1; Figures 5 and S1).It should be noted that the mechanisms behind these inhibitory effects are not the same.This is because excessive SWC creates water stress, which reduces plant net primary productivity and even leads to plant death [61][62][63].This, in turn, results in a reduction in the input of soil carbon.On the other hand, a higher soil pH (pH > 8) may increase the desorption and release of SOC and DOC from the soil, ultimately leading to a decrease in SOC and DOC contents [21,64].
In summary, our research results suggest that the main factors affecting the soil carbon pool in reclaimed saline-alkali ditches are TN, AP, EC, SWC, and pH.An increase in TN and AP contents raises carbon contents, while an increase in EC, SWC, and pH reduces carbon contents.However, the combined changes in these factors result in the overall stability of the soil carbon pool, with variations observed between POM and MAOM fractions in saline-alkali ditches.Although the continuous accumulation of SOC MAOM indicates an ongoing increase in soil carbon pool stability during reclamation (Figure 2), the build−up of SOC MAOM in the short−term reclamation process is still insufficient compared to the entire soil carbon pool.It requires the maintenance of high biomass in the saline-alkali ditches to increase the soil MAOM fraction.Compared to the saline−alkali farmland, saline-alkali ditches are hardly directly affected by agricultural activities after reclamation.These activities include physical soil disturbance caused by agricultural machinery and the application of fertilizers and pesticides, as well as irrigation and drainage.Instead, they are indirectly impacted, which causes the saline-alkali ditch soil carbon pool to be controlled by farmland management measures to some extent.Therefore, we suggest emptying as much water as possible from the drainage ditches following irrigation in the adjacent saline-alkali farmland.The drainage process would reduce the negative effects caused by the waterlogging and salinity in the ditches, which could accelerate soil nutrient utilization, promote plant biomass growth and microbial activities, and may increase the carbon stock of the saline-alkali ditches.
Although our data have shown the potential transformation of SOC POM to SOC MAOM in the soil carbon pool, more direct evidence is still required to confirm this process.To gain a more comprehensive and specific understanding of this process, it is necessary to utilize techniques such as isotopic and infrared spectroscopy, as well as different mathematical models to simulate soil carbon dynamics.In addition, to fully evaluate the soil carbon stocks in the saline-alkali ditches, more sufficient data relating to the entire soil profiles, equivalent soil masses, soil bulk density, and a sufficient number of samples are needed in the following studies [65,66].In conclusion, our results have enriched the academic understanding of the formation and transformation mechanisms of POM and MAOM.Our results further indicate a unique carbon fixation process in coastal saline-alkali ditches, that is, although the storage has not increased following the reclamation years, the stability has increased (Figure 2).As far as we know, the present study firstly provides a case study contributing to a comprehensive and reasonable evaluation of the impact of reclamation of saline-alkali farmland on the soil carbon pool, providing a theoretical basis for the management of the soil carbon pool in the reclaimed saline-alkali areas.

Conclusions
In this study, we investigated the dynamics of soil carbon pools in the saline-alkali ditches of the Yellow River Delta using the spatio-temporal substitution method.We also explored the changes in soil carbon components (POM and MAOM) and their stocks in drainage ditches of the reclaimed saline-alkali farmland.In the POM fraction, the SOC POM content and M PSOC significantly decreased from 2.24 to 1.12 g kg −1 and from 19.02 to 12.71 Mg ha −1 , respectively.Conversely, in the MAOM fraction, the SOC MAOM content and M MSOC significantly increased from 0.65 to 1.70 g kg −1 and from 5.30 to 12.27 Mg ha −1 , respectively.These results collectively indicate that the potential transformation of SOC POM to SOC MAOM and the mineralization process of SOC POM and SOC MAOM to SIC were the inherent processes that maintain the constant carbon stock of saline-alkali ditch soils.The changes in TN, AP, EC, SWC, and pH over the reclamation time were external factors.Notably, the accumulation of SOC MAOM after reclamation could enhance the stability of the saline-alkali ditch's soil carbon pool.The results undoubtedly demonstrate that although the storage has not increased following the reclamation years, the stability has increased.As far as we know, the present study firstly provides a case study contributing to a comprehensive and reasonable evaluation of the impact of reclamation of saline-alkali farmland on the soil carbon pool, providing a theoretical foundation for managing the soil carbon pool in recovered saline-alkali areas.
This research was conducted at the Yellow River Delta Research Center, Chinese Academy of Sciences, in Dongying, Shandong Province, China (37 • 40 N, 118 • 55 E).The study area is situated in the alluvial deposition zone of the Yellow River Delta, characterized by saline-alkali soils.The area experiences a temperate continental monsoon climate, with an average annual temperature of 12.40 • C. The annual average precipitation ranges from 530-630 mm, with the majority occurring between June and September.Additionally, the annual evaporation is 1962 mm.Over recent years, extensive land development has occurred, resulting in the creation of numerous drainage ditches (Figure 1).The soil in this region has originated from the sediment of the Yellow River since 1855, and the main soil type is Calcaric Fluvisols, comprising 5.90% caly, 70.88% silt, and 23.22% sand.The vegetation in these drainage ditches primarily comprises reed (Phragmites australis (Cav.)Trin.ex Steud.),Suaeda salsa (Suaeda salsa (L.) Pall.), Imperata cylindrica (Imperata cylindrica (L.) P. Beauv.), Leymus chinensis (Leymus chinensis (Trin.), and Artemisia capillaris (Artemisia capillaris Thunb.).Agronomy 2023, 13, x FOR PEER REVIEW 3 of 18

Figure 1 .
Figure 1.Distribution of sampling sites in the study area and the schematic diagram of the salinealkali ditch.

Figure 1 .
Figure 1.Distribution of sampling sites in the study area and the schematic diagram of the salinealkali ditch.

Figure 2 .
Figure 2. Effect of different reclamation years and slope sampling positions on (a) M TC , (b) M SOC , (c) M SIC , (d) M PSOC , (e) M MSOC , (f) M PSIC , and (g) M MSIC at average level per 20 cm in 0-60 cm soil depth.Different lowercase letters denote significant differences among different slope sampling positions in the same reclamation year (p < 0.05); different uppercase letters denote significant differences among different reclamation years (p < 0.05).(M TC , TC stock; M SOC , SOC stock; M SIC , SIC stock; M PSOC ).

Figure 3 .
Figure 3. Relationship between SOC and SIC (a,d,g), SOCPOM and SICPOM (b,e,h), and SOCMAOM and SICMAOM (c,f,i) under the grouping of reclamation years and depth.All fitting curves are at the significance level of p < 0.05.(SOC, soil organic carbon; SIC, soil inorganic carbon; SOCPOM, SOC in POM fraction; SOCMAOM, SOC in MAOM fraction; SICPOM, SIC in POM fraction; and SICMAOM, SIC in MAOM fraction).

Figure 3 .
Figure 3. Relationship between SOC and SIC (a,d,g), SOC POM and SIC POM (b,e,h), and SOC MAOM and SIC MAOM (c,f,i) under the grouping of reclamation years and depth.All fitting curves are at the significance level of p < 0.05.(SOC, soil organic carbon; SIC, soil inorganic carbon; SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; and SIC MAOM , SIC in MAOM fraction).Agronomy 2023, 13, x FOR PEER REVIEW 10 of 18

Figure 4 .
Figure 4. Direct and indirect effects of the reclamation years on SOCPOM, SOCMAOM, SICPOM, SICMAOM, DOC, and SIC in the saline-alkali ditches.The solid arrows represent a positive relationship and the dashed arrows represent a negative relationship, respectively.The width of the arrow is proportional to the strength of the relationship, and the adjacent numbers on arrows represent standardized path coefficients.R 2 represents the proportion of variance explained by the model.* p < 0.05, ** p < 0.01, and *** p < 0.001.(SIC, soil inorganic carbon; DOC, soil dissolved organic carbon; SOCPOM, SOC in POM fraction; SOCMAOM, SOC in MAOM fraction; SICPOM, SIC in POM fraction; and SICMAOM,

Figure 4 .
Figure 4. Direct and indirect effects of the reclamation years on SOC POM , SOC MAOM , SIC POM , SIC MAOM , DOC, and SIC in the saline-alkali ditches.The solid arrows represent a positive relationship

Figure 6 .
Figure 6.Conceptual framework for the conversion of soil organic carbon and inorganic carbon based on reclamation time and soil carbon components.The solid arrows represent a conversion

Figure 6 .
Figure 6.Conceptual framework for the conversion of soil organic carbon and inorganic carbon based on reclamation time and soil carbon components.The solid arrows represent a conversion relationship, and the dashed arrow represents no conversion relationship.The width of the arrow is proportional to the strength of the conversion relationship.(SIC, soil inorganic carbon; DOC, soil dissolved organic carbon; SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; SIC MAOM , SIC in MAOM fraction.).

Table 1 .
Soil properties of the saline-alkali ditches under the different reclamation years and soil depths.

Table 1 .
Cont.Values are mean ± standard error.Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.(TC, soil total carbon; SOC, soil organic carbon; SIC, soil inorganic carbon; DOC, soil dissolved organic carbon; SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; SIC MAOM , SIC in MAOM fraction; TN, soil total nitrogen; NH + 4 -N, soil soluble ammonium nitrogen; NO - 3 -N, soil soluble nitrate nitrogen; TP, soil total phosphorus; AP, soil available phosphorus; SWC, soil water content; pH, soil pH; and EC, soil electrical conductivity).

Table 2 .
The results (F values) of three-way ANOVA for soil properties in the saline-alkali ditches.

Table 3 .
Importance and significance test of soil properties' contribution rankings.