Enhanced Soil Carbon Stability through Alterations in Components of Particulate and Mineral-Associated Organic Matter in Reclaimed Saline–Alkali Drainage Ditches

: Soil carbon content and stability are primarily influenced by the stabilization of particulate organic matter (POM) and mineral-associated organic matter (MAOM). Despite extensive research on the stabilization processes of POM and MAOM carbon components under various land-use types, the investigation into stabilization processes of soil carbon remains limited in saline–alkali soils. Therefore, we collected soil samples from different positions of saline–alkali drainage ditches at four reclamation times (the first, seventh, fifteenth, and thirtieth year) to determine their carbon content and physicochemical properties. Moreover, POM and MAOM fractions were separated from soil samples, and Fourier transform infrared spectra (FTIR) were used to investigate changes in their chemical composition. The results showed that with increasing reclamation time, the soil total carbon and soil organic carbon (SOC) contents significantly increased from 14 to 15 and 2.9 to 5.5 g kg − 1 , respectively. In contrast, soil inorganic carbon content significantly decreased from 11 to 9.6 g kg − 1 . Notably, the changes in soil carbon components following the increasing reclamation time were primarily observed in the furrow sole at a depth of 20–40 cm. While the SOC content of the POM fraction (SOC POM ) decreased significantly, the SOC content of the MAOM fraction (SOC MAOM ) increased significantly. These alterations were largely dominated by drainage processes after reclamation instead of a possible conversion from SOC POM to SOC MAOM . FTIR results revealed that MAOM was greatly influenced by the reclamation time more than POM was, but the change in both POM and MAOM contributed to an increase in soil carbon stability. Our findings will deepen the comprehension of soil carbon stabilization processes in saline–alkali drainage ditches after reclamation and offer a research framework to investigate the stability processes of soil carbon components via alterations in POM and MAOM fractions.


Introduction
Soil organic carbon (SOC) is a crucial component of the terrestrial soil carbon pool, accounting for approximately two-thirds of the carbon stocks in terrestrial ecosystems, surpassing the combined carbon stock in the global atmosphere and vegetation [1][2][3].Maintaining and increasing SOC content is crucial for enhancing ecosystem productivity, mitigating climate change, and other related aspects [4,5].However, the dynamics and stability of SOC are changed due to continuous alterations in climate, vegetation conditions, and land-use changes [6][7][8].Land-use changes, particularly the conversion of fallow Agronomy 2024, 14, 869 2 of 17 land to farmland, are likely to lead to SOC losses [9,10].Despite previous extensive studies suggesting that SOC sequestration is a key measure to address global climate change [11][12][13], the current understanding of SOC stability processes under land-use changes remains limited.
Considering the regulatory role of soil organic matter (SOM) in the formation, sequestration, and stabilization of SOC, solely understanding the stability processes of SOC from a SOC perspective may have limitations [14][15][16][17].SOM can be categorized into two main components, particulate organic matter (POM) and mineral-associated organic matter (MAOM), which exhibit distinct turnover rates and stabilization processes [16][17][18][19].Despite the MAOM fraction being acknowledged as an important pool for long-term C stabilization in the soil [20][21][22], recent research has emphasized the stability processes with various carbon components.Islam et al. [23] proposed a stability framework for SOM, suggesting that root exudates, litter, and microbial products collectively contribute to the formation of POM and MAOM through different pathways.On the other hand, Liao et al. [17] argued that the formation and stabilization of soil carbon follow the sequence of "litter→POM→MAOM".The stability of SOM depends not only on the properties of POM and MAOM but also on the interactions of physical and chemical processes in the soil conditions [24].The inconsistency in these stabilization processes is largely attributed to variations in environmental factors and is impacted by changes in the SOC stability through POM and MAOM components.Consequently, it is crucial to examine how POM and MAOM components respond to environmental factors in the stabilization processes.
In the context of global change, the impact of land-use change on POM and MAOM components has been confirmed [16].Early research indicates that in undisturbed soils, the agricultural planting process following land reclamation can lead to significant differences between POM and MAOM, making POM more susceptible to loss than MAOM [25][26][27].However, when considering land-use changes in saline-alkali soils, specifically the salinealkali drainage ditches formed after saline-alkali land reclamation, alterations in soil carbon components are primarily influenced by vegetation status and drainage processes.On the one hand, excessive soil salinity in saline-alkali areas restrains plant growth, reducing the carbon input from plants into the soil [28,29]; on the other hand, drainage activities often lead to frequent wet-dry cycles in drainage ditch soils, facilitating the breakdown of soil aggregates, impacting microbial activity, and exposing and decomposing, physically protected SOM [30,31].The collective changes in these environmental factors are likely to influence the alterations and stability of soil carbon components.While current studies concentrate on investigating the stability processes and mechanisms of POM and MAOM components in various contexts [3,18,32,33], the research primarily focuses on conventional land uses like farmland, forests, and shrubs, with minimal attention given to the carbon stabilization processes in saline-alkali soils.
A comprehensive understanding of the stable process of soil carbon components is needed for examining the changes in their chemical composition.Infrared spectroscopy is commonly utilized for the qualitative analysis of functional groups in SOM and the quantitative analysis of SOC [34,35].In the physical fractionation of SOM, Fourier transform infrared (FTIR) spectroscopy serves as a supplementary method to gain insights into the chemical composition of soil carbon components [36].Furthermore, it facilitates the qualitative analysis of the functional group composition of SOM without changing the soil chemical composition [37].The continuous advancements in infrared spectroscopy technology have seen the rise of diffuse reflectance (DRIFT) and attenuated total reflection (ATR) FITR spectroscopy as alternative options to FTIR spectroscopy [38].While DRIFT and ATR FITR spectroscopy allow spectral analysis without the need for traditional KBr pellet preparation, they exhibit minimal changes in the spectrum compared to transmission FTIR spectra under the same analysis conditions [39,40].At present, FTIR spectroscopy is a more established method for determining the relative composition of specific functional groups in SOM, maintaining comparability with DRIFT and ATR FITR spectroscopy [38,41].Furthermore, to enhance the precision of FTIR spectroscopy, the utilization of a signif-icant amount of spectroscopy data is imperative to minimize systematic errors during spectral scanning.
Therefore, the objective of this study is to reveal the stabilization processes of soil carbon components in reclaimed saline-alkali ditches and their responses to environmental factors.To address this, the agricultural farmland in the Yellow River Delta region was chosen as the study area.We measured soil carbon contents and physicochemical properties of the saline-alkali drainage ditches at four reclamation times (the first, seventh, fifteenth, and thirtieth years).In addition, to gain further insight into the stable process of soil carbon components, the POM and MAOM fractions were isolated and analyzed for changes in their chemical composition using FTIR spectroscopy.The findings of this study will contribute to the development of a theoretical foundation for effectively enhancing the carbon content and stability of saline-alkali soils.

Site Description
The study area is situated in Dongying City, Shandong Province, China, near the research center base of the Yellow River Delta, Chinese Academy of Sciences (37 • 40 ′ N, 118 • 55 ′ E).With an average annual temperature of 12.60 • C and an average annual rainfall of around 465 mm, it has a moderate continental monsoon climate.There is an annual variation in evaporation of 1900-2400 mm.The soil in this study area comes from the alluvial deposits of the Yellow River.The soil type is classed as sandy loam by the United States Department of Agriculture, with sand particles (0.05-2.00 mm) accounting for 23.22%, silt particles (0.002-0.05 mm) accounting for 70.88%, and clay particles (<0.002 mm) accounting for 5.90%.Numerous drainage ditches have been formed as a result of the reclamation of saline-alkali farmland in recent decades.Phragmites australis (Cav.)Trin.ex Steud., Suaeda salsa (L.) Pall., Imperata cylindrica (L.) P. Beauv., Leymus chinensis (Trin.)Tzvel., Artemisia capillaris Thunb., and Tamarix chinensis Lour are the main plant species found in these drainage ditches.

Soil Sampling Processing
Four reclamation years (the first, seventh, fifteenth, and thirtieth years) within a 3 km radius of one another in the study area were used to choose sample locations in order to examine differences in the soil carbon components of saline-alkali drainage ditches (Figure 1).Four nearby drainage ditches with comparable cross-sectional dimensions, channel lengths, and plant covers were selected as sampling plots at each sampling location.Three randomly selected 1 m × 1 m subplots were chosen at various positions within each sampling plot, including the furrow sole, middle slope, lower slope, upper slope, and slope crest.For every subplot, plant biomass was measured by harvesting the above-ground biomass after plant height and vegetation covering were noted.The types and growth status of native vegetation in drainage ditches with different reclamation years are provided in Table 1.Within each subplot, soil samples were taken with a soil auger at three different depths: 0-20 cm, 20-40 cm, and 40-60 cm.Soil samples were taken from this location at depths of 0-20 cm and 20-40 cm due to the complicated soil and water conditions at the furrow sole.Soil samples from the same position and depth inside the three slope sections were mixed to generate a composite sample that represented each drainage ditch after plant litter and roots were removed using a 2 mm sieve.There were 224 composite samples gathered in all.

Fraction Separation, Soil Carbon Content, and Physicochemical Properties Analyses
According to Cotrufo et al. [19], POM and MAOM fractions were separated by the physical fractionation method.In this method, 10 g of air-dried soil was dispersed in a 30 mL solution containing 5 g L −1 of (NaPO3)6.The soil suspension was then shaken on a reciprocating shaker for 15 h to distribute soil particles.Subsequently, the distributed soil mixture was filtered using a 53 µm sieve and repeatedly washed several times with deionized water.The material retained on the sieve was identified as the POM fraction (53-2000 µm), while the material that passed through the sieve was identified as the MAOM fraction (<53 µm).
The soil total carbon (TC) and soil inorganic carbon (SIC) contents were measured using an automated carbon nitrogen analyzer (Primacs SNC100-IC-E, Skalar Analytical BV, Breda, The Netherlands).The TC content was determined through high-temperature combustion and nondispersive infrared detection (NDIR), while the SIC content was determined through automatic acidification, heating, gas stripping, and NDIR detection.The SOC content was calculated by subtracting the SIC from the TC content.Moreover, the analysis also included the determination of organic carbon (SOCPOM and SOCMAOM) and inorganic carbon (SICPOM and SICMAOM) contents in the POM and MAOM fractions using the same analytical methodology.Using a continuous flow analyzer (Auto Analyzer 3, Bran and Luebbe GmbH, Norderstedt, Germany), the accessible phosphorus (AP) level of the soil was determined.Deionized water was used to create a 1:5 soil-water mixture.A pH meter (PHS-3E, Leici Instruments, Shanghai, China) and a conductivity meter (FE38-Standard, Mettler Toledo, Greifensee, Switzerland) were used to measure the soil's pH and electrical conductivity (EC), respectively.Soil total nitrogen (TN) content was determined using the Kjeldahl method, while the soil total phosphorus (TP) content was determined using the H2SO4-HClO4 digestion method.Continuous flow analyzers were used to measure both TN and TP.Using a modified technique created by Zhan and Zhou, the dissolved organic carbon (DOC) content of the soil was ascertained [42].The KCIindophenol blue colorimetric method and the dual-wavelength UV spectrophotometric approach were used to quantify the amounts of soil soluble ammonium nitrogen (NH 4 + -N) and soil soluble nitrate nitrogen (NO 3 -N).A UV-visible spectrophotometer was used for all three measurements.

Chemical Composition Analysis of Carbon Components
A FTIR spectrometer (Nicolet iS10, ThermoFisher, Madison, WI, USA) was employed to analyze the infrared spectra of soil samples from the POM and MAOM fractions.Before conducting infrared spectrum scanning, soil samples from distinct sections with the same reclamation time were individually mixed (56 samples each) to minimize errors during

Fraction Separation, Soil Carbon Content, and Physicochemical Properties Analyses
According to Cotrufo et al. [19], POM and MAOM fractions were separated by the physical fractionation method.In this method, 10 g of air-dried soil was dispersed in a 30 mL solution containing 5 g L −1 of (NaPO 3 ) 6 .The soil suspension was then shaken on a reciprocating shaker for 15 h to distribute soil particles.Subsequently, the distributed soil mixture was filtered using a 53 µm sieve and repeatedly washed several times with deionized water.The material retained on the sieve was identified as the POM fraction (53-2000 µm), while the material that passed through the sieve was identified as the MAOM fraction (<53 µm).
The soil total carbon (TC) and soil inorganic carbon (SIC) contents were measured using an automated carbon nitrogen analyzer (Primacs SNC100-IC-E, Skalar Analytical BV, Breda, The Netherlands).The TC content was determined through high-temperature combustion and nondispersive infrared detection (NDIR), while the SIC content was determined through automatic acidification, heating, gas stripping, and NDIR detection.The SOC content was calculated by subtracting the SIC from the TC content.Moreover, the analysis also included the determination of organic carbon (SOC POM and SOC MAOM ) and inorganic carbon (SIC POM and SIC MAOM ) contents in the POM and MAOM fractions using the same analytical methodology.Using a continuous flow analyzer (Auto Analyzer 3, Bran and Luebbe GmbH, Norderstedt, Germany), the accessible phosphorus (AP) level of the soil was determined.Deionized water was used to create a 1:5 soil-water mixture.A pH meter (PHS-3E, Leici Instruments, Shanghai, China) and a conductivity meter (FE38-Standard, Mettler Toledo, Greifensee, Switzerland) were used to measure the soil's pH and electrical conductivity (EC), respectively.Soil total nitrogen (TN) content was determined using the Kjeldahl method, while the soil total phosphorus (TP) content was determined using the H 2 SO 4 -HClO 4 digestion method.Continuous flow analyzers were used to measure both TN and TP.Using a modified technique created by Zhan and Zhou, the dissolved organic carbon (DOC) content of the soil was ascertained [42].The KCIindophenol blue colorimetric method and the dual-wavelength UV spectrophotometric approach were used to quantify the amounts of soil soluble ammonium nitrogen (NH + 4 -N) and soil soluble nitrate nitrogen (NO − 3 -N).A UV-visible spectrophotometer was used for all three measurements.

Chemical Composition Analysis of Carbon Components
A FTIR spectrometer (Nicolet iS10, ThermoFisher, Madison, WI, USA) was employed to analyze the infrared spectra of soil samples from the POM and MAOM fractions.Before conducting infrared spectrum scanning, soil samples from distinct sections with the same reclamation time were individually mixed (56 samples each) to minimize errors during the scanning process.The samples were then combined with KBr in a 1:100 ratio, ground, and compressed into pellet molds using an agate mortar.Subsequently, the Fourier transform infrared spectrometer was used to measure the infrared transmittance spectra of the pellets with a spectroscopically pure KBr pellet serving as the background reference.The spectral range was set between 4000 and 400 cm −1 , with a resolution of 4 cm −1 and a scanning frequency of 32 times.The analysis of the infrared spectra was performed using OMNIC software (v.8.2).
In this study, we chose Liu et al.'s [43] research method based on previous research and review [44,45] due to the similarity of soil deposition and vegetation restoration conditions in our study area.To examine alterations in the chemical composition of soil carbon fractions during the reclamation process, we distinguished four specific functional groups in the infrared spectra of POM and MAOM soil fractions.These functional groups include alcohol C-OH/C-NH bonds (peak height at 3424 cm −1 , range 3300-3700 cm −1 ), aliphatic C-H bonds (peak height at 2925 and 2850 cm −1 , range 2800-2950 cm −1 ), aromatic C=C/C=N bonds (peak height at 1638 cm −1 , range 1350-1700 cm −1 ), and polysaccharide C-O bonds (peak height at 1030 cm −1 , range 1000-1100 cm −1 ).Peak areas were calculated after normalizing all FTIR spectra to represent the relative proportions of different functional groups.The extended multiplicative signal correction method, as detailed by Margenot et al. [46], was applied to adjust the baseline of the spectra and eliminate biases and noise.
Additionally, this study employed the methodology outlined by Pärnpuu et al. [47] to assess the decomposition extent of various carbon components and their reactions to the water environment in saline-alkali drainage ditches.The evaluation included parameters such as soil hydrophobicity (W), hydrophobicity index (HI), and decomposition degree index (DDI).W was computed by calculating the ratio of the sum of peak heights at 2921 and 2852 cm −1 to the peak height at 1633 cm −1 in the infrared spectra [48,49].HI was determined by evaluating the ratio of the sum of peak heights at 2921 and 2852 cm −1 to the SOC content in the infrared spectra [50].DDI was calculated by dividing the peak height at 1633 cm −1 by the sum of peak heights at 2921 and 2852 cm −1 in the infrared spectra [51][52][53].

Data Analysis
After checking the normality (Shapiro-Wilk's test) and homogeneity of variances (Levene's test) for all measured indicators, one-way ANOVA with least significant difference (LSD) was used to examine the differences in soil properties (TC, SOC, SIC, DOC, TN, NH + 4 -N, NO − 3 -N, TP, AP, pH, EC, SWC) of different reclamation years, soil carbon component contents (SOC POM , SOC MAOM , SIC POM , SIC MAOM ), relative proportions of different functional groups (alcohol C-OH/C-NH, aliphatic C-H, aromatic C=C/C=N, polysaccharide C-O), and W, HI, and DDI among different reclamation years, soil depths, and sampling positions.A significance level of p < 0.05 was considered statistically significant for all values, and all statistical analyses were conducted using SPSS (v22.0.0.0) software.
In this study, a multiple regression model was employed to assess the effects and relative importance of vegetation status (dry weight and coverage), soil nitrogen (TN, NH + 4 -N, and NO − 3 -N), soil phosphorus (TP and AP), and soil and water conditions (SWC, EC, and pH) on different soil carbon components.To standardize all environmental factors and response variables, a model was constructed using the MuMIn package in the R language (version 4.3.2),incorporating all possible combinations of initial influencing factors.The models were subsequently ranked based on the Akaike information criterion (AIC) derived from maximum likelihood fitting, and the fit best model with ∆AIC < 2 was selected.Finally, model averaging was employed to estimate the parameters and relevant p-values.The relative effects of each influencing factor were calculated and compared with the effects of all parameter estimates.The analyses were implemented using the MuMIn, rdacca.hp,ggplot2, and cowplot packages in R language.

Changes in Vegetation Condition and Soil Physicochemical Properties after Reclamation
Changes in vegetation status after reclamation are presented in Table 1.Plant height, dry weight, and vegetation coverage had significant variations during the different reclamation years.As reclamation time increased, plant height, dry weight, and coverage significantly increased from 73 to 150 cm, 0.46 to 0.76 kg, and 54% to 97%, respectively (p < 0.05).Additionally, the species of main vegetation in the saline-alkali ditches increased gradually.
Changes in soil properties following reclamation are presented in Table 2.As reclamation time increased, TN and TP contents significantly increased from 140 to 340 mg kg −1 and 0.48 to 0.53 g kg −1 , respectively (p < 0.05).NO − 3 -N content increased from 19 mg kg −1 to 20 mg kg −1 in the 1st year compared to the 15th year of reclamation, subsequently decreasing to 8.4 mg kg −1 in the 30th year (p < 0.05).AP content decreased significantly from 26 to 12 mg kg −1 (p < 0.05), while NH + 4 -N content exhibited no significant variations across different reclamation years (p > 0.05).SWC and pH significantly increased from 26% to 30% and 8.5 to 9.1, respectively, whereas EC decreased significantly from 3600 to 290 µs cm −1 (p < 0.05).

Changes in Soil Carbon Components after Reclamation
Overall, significant changes in soil carbon fractions following reclamation can be observed (Table 2, Figure 2).With increasing reclamation time, TC and SOC contents significantly increased from 14 to 15 g kg −1 and 2.9 to 5.5 g kg −1 , respectively, while SIC content decreased significantly from 11 to 9.6 g kg −1 (p < 0.05).DOC content decreased significantly from 48 to 47 mg kg −1 in the 1st year compared to the 7th year of reclamation, then significantly increased to 78 mg kg −1 in the 30th year (p < 0.05).Furthermore, SOC POM , SOC MAOM , and SIC MAOM contents significantly increased from 2.2 to 3.6 g kg −1 , 0.7 to 1.9 g kg −1 , and 1.2 to 6.3 g kg −1 , respectively, while SIC POM content decreased from 9.9 g kg −1 to 3.5 g kg −1 (p < 0.05).Specifically, the changes in soil carbon components after reclamation mainly occurred in the soil layer at a depth of 20-40 cm (Figure 3).From the 7th year to the 15th year, there was a significant decrease in SOC POM content and a significant increase in SOC MAOM content in furrow sole.From the 7th year to the 30th year, there was a significant decrease in SOC MAOM content (p < 0.05).

Changes in Soil Carbon Components after Reclamation
Overall, significant changes in soil carbon fractions following reclamation can be observed (Table 2, Figure 2).With increasing reclamation time, TC and SOC contents significantly increased from 14 to 15 g kg −1 and 2.9 to 5.5 g kg −1 , respectively, while SIC content decreased significantly from 11 to 9.6 g kg −1 (p < 0.05).DOC content decreased significantly from 48 to 47 mg kg −1 in the 1st year compared to the 7th year of reclamation, then significantly increased to 78 mg kg −1 in the 30th year (p < 0.05).Furthermore, SOCPOM, SOCMAOM, and SICMAOM contents significantly increased from 2.2 to 3.6 g kg −1 , 0.7 to 1.9 g kg −1 , and 1.2 to 6.3 g kg −1 , respectively, while SICPOM content decreased from 9.9 g kg −1 to 3.5 g kg −1 (p < 0.05).Specifically, the changes in soil carbon components after reclamation mainly occurred in the soil layer at a depth of 20-40 cm (Figure 3).From the 7th year to the 15th year, there was a significant decrease in SOCPOM content and a significant increase in SOCMAOM content in furrow sole.From the 7th year to the 30th year, there was a significant decrease in SOCMAOM content (p < 0.05).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).

Changes in Soil Carbon Components after Reclamation
Overall, significant changes in soil carbon fractions following reclamation can be observed (Table 2, Figure 2).With increasing reclamation time, TC and SOC contents significantly increased from 14 to 15 g kg −1 and 2.9 to 5.5 g kg −1 , respectively, while SIC content decreased significantly from 11 to 9.6 g kg −1 (p < 0.05).DOC content decreased significantly from 48 to 47 mg kg −1 in the 1st year compared to the 7th year of reclamation, then significantly increased to 78 mg kg −1 in the 30th year (p < 0.05).Furthermore, SOCPOM, SOCMAOM, and SICMAOM contents significantly increased from 2.2 to 3.6 g kg −1 , 0.7 to 1.9 g kg −1 , and 1.2 to 6.3 g kg −1 , respectively, while SICPOM content decreased from 9.9 g kg −1 to 3.5 g kg −1 (p < 0.05).Specifically, the changes in soil carbon components after reclamation mainly occurred in the soil layer at a depth of 20-40 cm (Figure 3).From the 7th year to the 15th year, there was a significant decrease in SOCPOM content and a significant increase in SOCMAOM content in furrow sole.From the 7th year to the 30th year, there was a significant decrease in SOCMAOM content (p < 0.05).

Chemical Composition in POM and MAOM Fractions
Four specific functional group peak areas of POM and MAOM in the soil showed significant variation among different reclamation years (p < 0.05), while they did not show significant differences among different soil depths and sampling locations (p > 0.05), as indicated by FTIR results (Tables 3, S1 and S2, Figure 4).The relative proportion of the polysaccharides C-O was the highest, while the relative proportion of aromatic C=C/C=N was the lowest.In the POM fraction, the relative proportion of alcohol C-OH/C-NH was greater than alkanes C-H, while in the MAOM fraction, the relative proportion of C-H was greater than C-OH/C-NH.With increasing reclamation time, there was no significant change observed in the C-O groups (p > 0.05).The C=C/C=N groups significantly decreased in the POM fraction but showed an increase in the MAOM fraction (p < 0.05).The C-OH/C-NH groups did not show a significant difference in the POM fraction (p > 0.05).However, C-OH/C-NH groups decreased from the 1st to the 15th year and significantly increased from the 15th to the 30th year of reclamation in the MAOM fraction (p < 0.05).Moreover, the C-H groups showed a significant increase (p < 0.05).The changes in the W, HI, and DDI of POM and MAOM fractions varied significantly with increasing reclamation years (Table 4).Conversely, no significant differences were observed among different soil depths and sampling locations (Tables S3 and S4).In the POM fraction, W decreased significantly in the 15th year of reclamation and increased significantly in the 30th year (p < 0.05); HI decreased significantly in the 30th year of reclamation (p < 0.05), while DDI remained unchanged (p > 0.05).In the MAOM fraction, W decreased significantly in the 15th year of reclamation and increased significantly in the 30th year (p < 0.05); DDI increased significantly in the 15th year of reclamation and decreased significantly in the 30th year (p < 0.05), with no significant changes observed for HI (p > 0.05).

The Main Factors Influencing the Changes in Soil Carbon Components
The multiple regression models for different carbon components are presented in Figure 5.In the SOCPOM model, a significant positive correlation was observed between SOCPOM and TN, while significant negative correlations were found with SWC, pH, and coverage (p < 0.05).Vegetation status, soil nitrogen, soil phosphorus, and the soil and water conditions explained 8.2%, 40%, 0.86%, and 1.3% of the variance in SOCPOM, respectively.In the case of SOCMAOM, significant positive correlations were identified between SOCMAOM and SWC, pH, and coverage, while a significant negative correlation was

The Main Factors Influencing the Changes in Soil Carbon Components
The multiple regression models for different carbon components are presented in Figure 5.In the SOC POM model, a significant positive correlation was observed between SOC POM and TN, while significant negative correlations were found with SWC, pH, and coverage (p < 0.05).Vegetation status, soil nitrogen, soil phosphorus, and the soil and water conditions explained 8.2%, 40%, 0.86%, and 1.3% of the variance in SOC POM , respectively.In the case of SOC MAOM , significant positive correlations were identified between SOC MAOM and SWC, pH, and coverage, while a significant negative correlation was observed with NH + 4 -N (p < 0.05).Vegetation status, soil nitrogen, and the soil and water conditions explained 13%, 7.6%, and 9.4% of the variance in SOC MAOM , respectively.For the SIC POM model, significant positive correlations were found with TN, SWC, NH + 4 -N, and AP, while a significant negative correlation was observed with pH (p < 0.05).Vegetation status, soil nitrogen, soil phosphorus, and the soil and water conditions contributed 6.3%, 43%, 11%, and 2.5% to the variance in SIC POM , respectively.Regarding the SIC MAOM model, it exhibited a significant positive correlation with coverage and a significant negative correlation with NH + 4 -N (p < 0.05).Vegetation status and soil nitrogen accounted for 6.5% and of the variance in SIC MAOM , respectively.
observed with NH 4 + -N (p < 0.05).Vegetation status, soil nitrogen, and the soil and water conditions explained 13%, 7.6%, and 9.4% of the variance in SOCMAOM, respectively.For the SICPOM model, significant positive correlations were found with TN, SWC, NH 4 + -N, and AP, while a significant negative correlation was observed with pH (p < 0.05).Vegetation status, soil nitrogen, soil phosphorus, and the soil and water conditions contributed 6.3%, 43%, 11%, and 2.5% to the variance in SICPOM, respectively.Regarding the SICMAOM model, it exhibited a significant positive correlation with coverage and a significant negative correlation with NH 4 + -N (p < 0.05).Vegetation status and soil nitrogen accounted for 6.5% and 50% of the variance in SICMAOM, respectively.

Discussion
As a consequence of the reclamation of saline-alkali drainage ditches, alterations in soil carbon content and components have occurred, alongside an improvement in the stability of the soil carbon pool.This study specifically reveals that TC and SOC contents in the saline-alkali drainage ditches rose as the reclamation time advanced, whereas SIC content decreased (Table 2).This contrasts with previous studies on reclaimed saline-alkali farmland, which reported an increase in both SOC and SIC contents [54].The changes in vegetation status and SWC of drainage ditches largely account for this discrepancy.On

Discussion
As a consequence of the reclamation of saline-alkali drainage ditches, alterations in soil carbon content and components have occurred, alongside an improvement in the stability of the soil carbon pool.This study specifically reveals that TC and SOC contents in the saline-alkali drainage ditches rose as the reclamation time advanced, whereas SIC content decreased (Table 2).This contrasts with previous studies on reclaimed saline-alkali farmland, which reported an increase in both SOC and SIC contents [54].The changes in vegetation status and SWC of drainage ditches largely account for this discrepancy.On the one hand, drainage ditches experience less disturbance compared to farmland and exhibit better vegetation status (Table 1), leading to increased carbon input from plants to the soil.On the other hand, the drainage process in saline-alkali ditches results in higher SIC loss due to increased carbonate dissolution caused by higher SWC (Table 2).This leads to more SIC being dissolved into the water body or groundwater of the ditches [11,55,56].However, the nuanced changes in soil carbon components more accurately illustrate the variations in SOC and SIC contents (Figure 2).As reclamation time increased, the SOC POM , SOC MAOM , and SIC MAOM contents increased, while the SIC POM content decreased, showing a more significant change compared to SIC MAOM .This suggests an overall increase in SOC content and a decrease in SIC content (Table 2).Moreover, the increasing SOC MAOM content also suggests an enhancement in stable components of saline-alkali drainage ditch soils (Figure 2d,e), ultimately enhancing the stability of the soil carbon pool [57][58][59].
The extended turnover time and enhanced stability of MAOM in comparison to POM establish MAOM as a crucial component in soil carbon sequestration [16,32,60].In this study, the proportion of unstable functional groups (polysaccharide C-O and alcohol C-OH/C-NH) was higher in the POM fraction than in the MAOM fraction (Figure 4).This finding is consistent with previous research and implies that MAOM possesses greater stability than POM [43,61].Changes in the peak area of functional groups in both POM and MAOM fractions further illustrate the enhancing stability of carbon components with the increasing reclamation time (Table 3).Specifically, as reclamation time increased, the C=C/C=N groups decreased and C-H groups increased in the POM fraction, while both C=C/C=N and C-H groups increased in the MAOM fraction (Table 3).This trend implies that increased reclamation time enhances the stability of both POM and MAOM fractions.Moreover, the changes in POM and MAOM components likely occurred through changes in C=C/C=N and C-H groups, as indicated by the decrease in SOC POM content and the increase in SOC MAOM content in the 15th year of reclamation (Figure 2b,c).Notably, the higher C-H groups in POM are associated with plant contributions [32,62,63], whereas the comparable distribution of C-H groups in both MAOM and POM fractions indicates that saline-alkali ditch plants contribute similarly to both organic matter fractions (Table 3).
By quantifying the degree of decomposition in the POM and MAOM fractions and their response to the water environment in saline-alkali drainage ditches, the results further demonstrate an increase in the stability of soil carbon components (Table 4).Specifically, the concentration of W in both the POM and MAOM fractions initially decreased from the 7th to the 15th year of reclamation and then increased from the 15th to the 30th year.This variation is primarily influenced by fluctuations in SWC.In anaerobic conditions, oxygen-containing functional groups serve as a source of oxygen, leading to an increase in C-H groups and a decrease in C-O groups, consequently elevating the concentration of W [64][65][66].This increase in W reflects enhanced aggregate stability [67].The inverse relationship between SWC and W confirms this association and implies that the aggregate stability of both the POM and MAOM fractions initially diminishes and then amplifies over time.Moreover, in the POM fraction, DDI remains constant; meanwhile, in the MAOM fraction, it decreases from the 7th to the 15th year of reclamation and subsequently increases from the 15th to the 30th year, consistently exceeding the values of the POM fraction (Table 4).This discrepancy indicates differences in the degree of decomposition between the POM and MAOM fractions.Despite changes occurring in the MAOM fraction, its stability consistently surpasses that of the POM fraction.The stable nature of the POM fraction further emphasizes the pronounced influence of the reclamation process on MAOM.
Although FTIR spectroscopy has enhanced our understanding of the stability changes in soil carbon components through chemical composition analysis, there are still limitations in the current FTIR spectroscopy analysis methods for the physically separated POM and MAOM fractions from soil.The obtained infrared spectroscopic data from these fractions may still be subject to uncertainties.One challenge is the potential peak shifts in the infrared spectra of POM and MAOM fractions due to the complex particle composition of soil samples, which can lead to misinterpretation of peaks as SOM or soil moisture.Additionally, KBr, a commonly used solid dispersant in infrared spectroscopic analysis, may not be completely inert and could react with metal halides present in the soil.However, FTIR spectroscopy continues to be a valuable method for qualitatively analyzing the composition of SOM functional groups.Our study offers additional insights into the stabilization process of soil carbon fractions, such as POM and MAOM, by analyzing a large number of FTIR spectra.To enhance the accuracy and reliability of infrared spectroscopic data, future research should focus on refining the resolution of using KBr pellets for infrared spectroscopic analysis and augmenting the spectral information with mineral fraction identification tables specific to soil samples.
The stabilization of soil carbon components was mainly affected by vegetation status and soil and water conditions in saline-alkali ditches after reclamation.Multiple regression analysis revealed that vegetation status and soil nitrogen play key roles in influencing these changes (Figure 5).This can be attributed to the fact that TN promotes SOC POM and SIC POM contents.Additionally, coverage facilitates SOC MAOM and SIC MAOM .Consequently, an increase in soil nitrogen not only improves vegetation growth status but also enhances the carbon input into the soil by plants [68].Improving vegetation growth status and increasing soil nitrogen levels thus prove to be effective approaches for augmenting carbon content in ditch soils.The inhibition of SIC POM by NH + 4 -N and the promotion of SOC MAOM and SIC MAOM in drained ditches indicate a likely association between NH + 4 -N and the generation of SIC, possibly promoting the formation and stabilization of MAOM.This relationship arises from the long-term fertilization processes in farmland, in which nutrient elements inevitably migrate into the soil and water via drainage.While this process enhances vegetation growth in ditches, the soil acidification resulting from nitrogen fertilizer utilization can lead to the loss of SIC [69][70][71][72].The distinct responses of SOC POM and SOC MAOM contents to variations in SWC and pH suggest that the combined influence of SWC and pH promotes SOC MAOM while inhibiting SOC POM .Increased SWC and pH reduce microbial activity, thereby enhancing the stability of SOC MAOM [73][74][75].Conversely, elevated SWC and pH levels inhibit plant growth, as excessive water stress affects normal plant development and reduces carbon input into the soil by plants [76,77].Additionally, a high pH soil condition (pH > 8) can result in SOC desorption and greater losses of relatively labile SOC POM [78,79].
Early models proposed that MAOM is formed from POM decomposition [1,80].However, recent research suggests separate formation and stabilization pathways for POM and MAOM, with no consistent evidence of POM converting into MAOM [81][82][83].Notably, a key observation is that the decrease in POM is not parallel with the increase in MAOM [83].In our study, the decrease in SOC POM and the increase in SOC MAOM contents occurred in the 15th year of reclamation, primarily at the 20-40 cm depth of the furrow sole (Figures 2 and 3g).However, this does not necessarily indicate a conversion from POM to MAOM, as it is largely influenced by the structural characteristics of the saline-alkali ditches.The plant soluble inputs accumulate at the furrow sole during drainage processes in saline-alkali ditches, forming MAOM when they interact with soil mineral surfaces at depths of 20-40 cm [23,[84][85][86].Conversely, plant structural inputs are unlikely to concentrate at the furrow sole during drainage processes; instead, more POM is formed relative to the furrow sole [23,81,86].
Enhancing soil carbon content and stability in saline-alkali ditches primarily relies on long-term drainage processes.To effectively reduce soil salinity, numerous drainage ditches have been constructed, leading to the influx of nutrients from farmland soils into the soil and water within ditches during extended drainage.Subsequent enhancements in vegetation status and soil and water conditions have led to an amount of carbon in the ditch soil.Despite the loss in SIC and reduced SOC stability caused by soil acidification from fertilizer application, the overall impact on the carbon content and stability in saline-alkaline ditch soil is relatively minor.Moreover, the potential of ditch structures to promote soil carbon sequestration may be underestimated.With the improvement in the soil and water conditions in the reclaimed ditches, vegetation restoration has led to the gradual formation and stabilization of aggregates, providing physical protection for soil carbon and promoting soil carbon sequestration [87][88][89].Specifically, in the initial drainage phase, the ditch structure promotes the accumulation of SOC MAOM at the furrow sole (Figures 3g and 6).The higher SWC and decreased pH at the furrow sole decrease microbial activity (Figure 3h), consequently increasing the stability of soil carbon.Conversely, in the last drainage phase, the decreased drainage frequency and volume result in alternate wetting and drying at the furrow sole soils (Figure 6).These alternating changes in the soil environment impact microbial activity, leading to the exposure and decomposition of SOC MAOM [30,31].Thus, to effectively enhance the soil carbon sequestration capacity of saline-alkali ditches, it is recommended to maintain water levels in the ditches beyond the 15th year of reclamation and to prevent mechanical damage to the ditch structure caused by agricultural management practices.
soil and water conditions in the reclaimed ditches, vegetation restoration has led to the gradual formation and stabilization of aggregates, providing physical protection for soil carbon and promoting soil carbon sequestration [87][88][89].Specifically, in the initial drainage phase, the ditch structure promotes the accumulation of SOCMAOM at the furrow sole (Figures 3g and 6).The higher SWC and decreased pH at the furrow sole decrease microbial activity (Figure 3h), consequently increasing the stability of soil carbon.Conversely, in the last drainage phase, the decreased drainage frequency and volume result in alternate wetting and drying at the furrow sole soils (Figure 6).These alternating changes in the soil environment impact microbial activity, leading to the exposure and decomposition of SOCMAOM [30,31].Thus, to effectively enhance the soil carbon sequestration capacity of saline-alkali ditches, it is recommended to maintain water levels in the ditches beyond the 15th year of reclamation and to prevent mechanical damage to the ditch structure caused by agricultural management practices.Following the reclaimed saline-alkali ditches, this study observed a general rise in the soil carbon pool, with the soil organic carbon pool making a more significant contribution compared to the soil inorganic carbon pool.Therefore, when concentrating on the response of soil carbon pools to environmental factors, it is crucial to focus on the distinct roles played by SOC and SIC in the soil carbon stabilization process.Moreover, the enhancement of vegetation conditions and soil and water conditions in saline-alkali ditches has altered the soil carbon components after reclamation, thereby improving soil carbon stability, largely due to long-term drainage processes.Although our findings shed light on these shifts in carbon components, future research will need to employ isotope tracing methods to reveal the underlying stabilization processes behind these alterations.To further elucidate the stabilization mechanisms of soil carbon components, future studies are needed to investigate the changes in soil carbon components under diverse environmental factors.Emphasizing field experiments over laboratory studies will aid in holistically considering the processes and pathways involved in soil carbon component stabilization.This approach will offer invaluable practical insights for the informed development of soil carbon pool management strategies.Following the reclaimed saline-alkali ditches, this study observed a general rise in the soil carbon pool, with the soil organic carbon pool making a more significant contribution compared to the soil inorganic carbon pool.Therefore, when concentrating on the response of soil carbon pools to environmental factors, it is crucial to focus on the distinct roles played by SOC and SIC in the soil carbon stabilization process.Moreover, the enhancement of vegetation conditions and soil and water conditions in saline-alkali ditches has altered the soil carbon components after reclamation, thereby improving soil carbon stability, largely due to long-term drainage processes.Although our findings shed light on these shifts in carbon components, future research will need to employ isotope tracing methods to reveal the underlying stabilization processes behind these alterations.To further elucidate the stabilization mechanisms of soil carbon components, future studies are needed to investigate the changes in soil carbon components under diverse environmental factors.Emphasizing field experiments over laboratory studies will aid in holistically considering the processes and pathways involved in soil carbon component stabilization.This approach will offer invaluable practical insights for the informed development of soil carbon pool management strategies.

Conclusions
In summary, this study presents a case study that investigates the processes underlying the stability of soil carbon components based on POM and MAOM fractions.The findings demonstrate that the reclamation process of saline-alkali drainage ditches changes the soil carbon content and components, ultimately enhancing soil carbon stability.As the reclamation time increased, there was a significant increase in TC and SOC contents, along with a significant decrease in SIC content.Within the POM fraction, there was a notable increase in SOC content and a decrease in SIC content.In the MAOM fraction, both SOC MAOM and SIC MAOM contents showed a significant increase.In addition, a decrease in SOC POM con-tent and an increase in SOC MAOM content with increasing reclamation time at the 20-40 cm depth of the furrow sole was largely dominated by drainage processes.FTIR results revealed that MAOM was greatly influenced by the reclamation time more than POM was, but the changes in both POM and MAOM contributed to an increase in soil carbon stability.Enhancing soil carbon content and stability in saline-alkali soils can be effectively achieved by utilizing the theoretical foundation provided by these findings, which deepen our understanding of soil carbon dynamics and the influence of environmental factors.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy14040869/s1, Figure S1: Typical FTIR spectra for POM and MAOM fractions with different reclamation years, soil depths, and soil sampling positions; Table S1: The specific peak areas of four functional groups in POM and MAOM fractions by FTIR with different soil depths; Table S2: The specific peak areas of four functional groups in POM and MAOM fractions by FTIR with different soil sampling positions; Table S3: Soil organic matter water repellence (W), hydrophobicity index (HI), and degree of decomposition (DDI) of POM and MAOM fractions with different soil depths; Table S4: Soil organic matter water repellence (W), hydrophobicity index (HI), and degree of decomposition (DDI) of POM and MAOM fractions with different soil sampling positions.
Author Contributions: Writing-review and editing, X.L., Z.L. and J.L.; methodology, X.L., Z.L. and Y.G.; validation, Z.L. and J.L.; data curation, X.L. and Y.G.; conceptualization, Z.L. and J.L.; project administration, Z.L.All authors have read and agreed to the published version of the manuscript.

Figure 1 .
Figure 1.Location of the study area and the distribution of sampling sites with different reclamation years.

Figure 1 .
Figure 1.Location of the study area and the distribution of sampling sites with different reclamation years.

Figure 2 .
Figure 2. Distribution of the soil carbon components (a), SOC POM (b), SOC MAOM (c), SIC POM (d), and SIC MAOM (e).Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.(SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; SIC MAOM , SIC in MAOM fraction.)

Figure 3 .
Figure 3. Soil carbon components distribution in soil sampling positions with different reclamation years and soil depths.Different lowercase letters of the same variable indicate significant differences between the different slope sampling positions at the level of p < 0.05.(SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; SIC MAOM , SIC in MAOM fraction.)

Figure 4 .
Figure 4. Relative proportions of the four specific functional groups in the POM and MAOM fractions with different reclamation years (a), soil depths (b), and soil sampling positions (c).

Figure 4 .
Figure 4. Relative proportions of the four specific functional groups in the POM and MAOM fractions with different reclamation years (a), soil depths (b), and soil sampling positions (c).

Figure 5 .
Figure 5. Relative effects of multiple environmental factors on SOCPOM (a,b), SOCMAOM (c,d), SICPOM (e,f), and SICMAOM (g,h).The averaged parameter estimates (standardized regression coefficients) of the model environmental factors are shown with their associated 95% confidence intervals (a,c,e,g) along with the relative importance of each environmental factor, expressed as the percentage of explained variance (b,d,f,h) for SOCPOM, SOCMAOM, SICPOM, and SICMAOM.The relative effect of the environmental factors is calculated as the ratio between the parameter estimate of the environmental factor and the sum of all parameter estimates, and it is expressed as a percentage.The vegetation status includes dry weight and coverage; soil nitrogen includes TN, NH 4 + -N, and SOCPOM, SOC in

Figure 5 .
Figure 5. Relative effects of multiple environmental factors on SOC POM (a,b), SOC MAOM (c,d), SIC POM (e,f), and SIC MAOM (g,h).The averaged parameter estimates (standardized regression coefficients) of the model environmental factors are shown with their associated 95% confidence intervals (a,c,e,g) along with the relative importance of each environmental factor, expressed as the percentage of explained variance (b,d,f,h) for SOC POM , SOC MAOM , SIC POM , and SIC MAOM .The relative effect of the environmental factors is calculated as the ratio between the parameter estimate of the environmental factor and the sum of all parameter estimates, and it is expressed as a percentage.The vegetation status includes dry weight and coverage; soil nitrogen includes TN, NH + 4 -N, and SOC POM , SOC in POM fraction; SOC MAOM , SOC in MAOM fraction; SIC POM , SIC in POM fraction; SIC MAOM , SIC in MAOM fraction; NO − 3 -N; soil phosphorus includes TP and AP; soil and water condition include SWC, pH, and EC.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.* p < 0.05.

Figure 6 .
Figure 6.Soil carbon stabilization process in saline-alkali ditches at different drainage stages after reclamation.

Figure 6 .
Figure 6.Soil carbon stabilization process in saline-alkali ditches at different drainage stages after reclamation.

Funding:
This research was funded by the National Natural Science Foundation of China (No. 32101596), the Key R&D Program of Shandong Province, China (2022TZXD0043), the Strategic Priority Research Program of the Chinese Academy of Sciences (no. XDA26050202), and the Shandong Provincial Natural Science Foundation (no.ZR2021QC055).

Table 1 .
The types and growth status of native vegetation in drainage ditches with different reclamation years.
Values are mean ± standard error.Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.

Table 2 .
Soil properties of the saline-alkali ditches with different reclamation years.
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; 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 3 .
The specific peak areas of four functional groups in POM and MAOM fractions by FTIR.
Values are mean ± standard error.Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.
Values are mean ± standard error.Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.
Values are mean ± standard error.Different lowercase letters of the same variable indicate significant differences at the level of p < 0.05.