Linking soil organic carbon mineralization to soil physicochemical properties and bacterial alpha diversity at different depths following land use changes

Anthropogenic land use changes (LUCs) impart intensifying impacts on soil organic carbon (SOC) turnover, leading to uncertainty concerning SOC mineralization patterns and determining whether soils act as “source” or “sink” in the global carbon budget. Therefore, understanding the SOC mineralization characteristics of different LUC patterns and their potential influencing factors is crucial. An indoor incubation experiment was conducted to study the SOC mineralization patterns and their relevance to soil physicochemical properties, soil enzyme activity, SOC fractions, and bacterial alpha diversity. The soils were collected from two layers of five typical LUC patterns in Yellow Sea Forest Park, including four that were converted from wheat–corn rotation systems [a gingko plantation (G), a metasequoia plantation (M), a gingko–wheat–corn agroforestry system (GW), and a gingko–metasequoia system (GM)] and a traditional wheat–corn system (W). LUCs had significant and diverse impacts on the SOC content and SOC fraction contents and on soil enzyme activity. The cumulative SOC mineralization was significantly higher in the M systen than in the W and GW systems at 0–20 cm depth and higher in the G system than in the GW system at 20–40 cm depth after 60-day incubation. The mineralization ratio was highest in the W system and lowest in the GW system. The soil pH and bulk density had a significant negative correlation with the cumulative SOC mineralization, while the soil bacterial Shannon index had a significant positive correlation with cumulative SOC mineralization. Multiple stepwise linear regression analysis showed that the SOC mineralization potential was dominantly explained by the bacterial Shannon index and operational taxonomic units (OTUs). The GW system had lower potentially mineralizable SOC and higher SOC stability. Additionally, the incubation time and cumulative SOC mineralization were well fitted by the first-order kinetic equation. LUCs significantly changed SOC mineralization characteristics and the results highlighted the important roles of the bacterial community in soil carbon cycling, which contributes to the fundamental understanding of SOC turnover regulation.


Introduction
Global average temperatures are expected to increase by approximately 2-6 °C above the temperatures measured before the Industrial Revolution by 2100, and the trend of temperature increase on a global scale is parallel to the trend of changes in atmospheric carbon dioxide (CO 2 ) (Al-Ghussain 2019; Letcher 2019).Soil organic carbon (SOC) pool is the largest organic carbon sink in the terrestrial biosphere; its extent depends on the dynamic balance between plant carbon input and microbial carbon outputs, and small changes in SOC will have a significant impact on atmospheric CO 2 concentration (Schmidt et al. 2011;Scharlemann et al. 2014;Bradford et al. 2016).Thus, the Intergovernmental Panel on Climate Change (IPCC) has indicated that both CO 2 removal and critical reductions in greenhouse gas (GHG) emissions from anthropogenic land use changes (LUCs) should be adopted to limit global warming to within 1.5 °C by 2100 (IPCC 2018;Anderson et al. 2019;Beillouin et al. 2022).Particularly, global warming will also accelerate the decomposition of SOC, thus creating a strong positive feedback with climate change (Wang et al. 2022).The soil carbon sink capacity has received extensive attention, and thus increasing SOC storage and reducing GHG emissions from soil are effective ways to address the climate change issues (Smith et al. 2016).The SOC stock exhibits significant spatial heterogeneity but is mainly affected by local-scale land use patterns and has been confirmed to be vulnerable to LUCs (Don et al. 2011).
Due to the possible transformation of soil carbon sinks and sources, a two-pronged approach, referring to broadening soil carbon input and reducing soil carbon expenditure, has become the basic strategy for carbon sequestration (Garcia-Pausas et al. 2008;Caddeo et al. 2019).LUCs driven by anthropogenic activities have induced the release of 1.6 Gt C year −1 into the atmosphere in the past decade (Friedlingstein et al. 2020).However, Poeplau and Don (2013) obtained an estimate of 21 ± 13 Mg•ha −1 of SOC accumulation after croplands have been changed to forests.Similarly, afforestation plays a significant role in carbon sequestration, as 39% of terrestrial ecosystem C is stored in forest ecosystems, 21.4% is stored as stand biomass and 78.6% is stored in the soil (Huang et al. 2012;Rahman et al. 2017;Zhang et al. 2021b).In general, afforested agricultural lands increase the SOC stock by 29%; this increase is relatively stable and further steadily increases with increasing tree age (Don et al. 2011).These studies have demonstrated the high C sequestration potential of afforestation, according to the aforementioned strategy, to further tap the C sequestration potential, and thoroughly elucidating the SOC output patterns is also of great significance.Namely, soil respiration is also key process affecting the global C balance and soil nutrient cycle, which is primarily driven by the microbial decomposition of SOC, and increases or decreases in the soil respiration rate can determine whether soil serves as a "source" or "sink" (Chen et al. 2018).There is a positive correlation between soil respiration and bacterial richness in general, and thus, changes in bacterial communities following LUCs definitely alter SOC mineralization patterns, mainly because microbes are the drivers of SOC mineralization and the microbial activities directly affect the SOC turnover rate (Tardy et al. 2015;Barnett et al. 2020).
SOC mineralization reflects the SOC decomposition process and is closely related to nutrient release, CO 2 emission, and soil property maintenance (Li et al. 2010).The specific mechanisms underlying SOC mineralization are not yet fully understood, but some key questions regarding SOC mineralization patterns have been answered.Temperature and moisture are known to be fundamental factors affecting carbon mineralization (Huang and Hall 2017;Schlüter et al. 2022).Intermediate moisture facilitates oxygen supply and substrate diffusion, and appropriate temperatures increase microbial activity, thus affecting SOC mineralization (Wang et al. 2016;Huang and Hall 2017).Other environmental factors (e.g., soil pH, particle size distribution, nutrients, enzyme activity, and iron oxides) and anthropogenic factors (e.g., land use patterns, tillage measures, exogenous organic matter additions, and fertilization) also critically impact SOC mineralization (Autret et al. 2020;Jin et al. 2020;Wei et al. 2020;Jeewani et al. 2021;Patel et al. 2021;Zhang et al. 2021a).The abovementioned influencing factors and their interactions make it incredibly challenging to comprehensively understand SOC mineralization (Rong et al. 2021).LUCs can lead to changes in the physical, chemical and biological conditions of the soil, which will promote or inhibit the release of SOC, possibly due to the altered soil microbial activity and/or community structure, and in turn alters the SOC decomposition rate.Therefore, it is important to investigate the pathways and mechanisms of LUC affecting SOC mineralization to gain a deeper understanding of the carbon sequestration and emission reduction role of afforestation and to help achieve carbon neutrality.
Laboratory incubation methods are widely used to study SOC mineralization, although incubation at constant temperature and humidity levels cannot reflect the actual SOC mineralization patterns in the field (Moinet and Millard 2020).However, these methods can still provide a better understanding of the SOC mineralization patterns under different land use patterns from the same origin.In this study, we selected four land use patterns converted from a wheat-corn rotation system and one legacy wheat-corn rotation system in the Yellow Sea Forest Park as the study area.We hypothesized that LUCs can affect SOC mineralization by regulating soil physicochemical properties and bacterial communities.Specifically, the main hypotheses were as follows: (1) LUCs would alter the SOC mineralization rate, i.e., lower SOC mineralization rates would be observed in the converted land use patterns than in the wheat-corn rotation system, and these differences would be controlled by certain factors among these land use patterns; and (2) the cumulative C mineralization would be controlled by different factors at different soil depths.To test these hypotheses, a 60-day incubation experiment was conducted.The results will deepen our understanding of the patterns and mechanism of SOC mineralization after LUCs in this area or similar areas.

Study site description
The study site is located in Yellow Sea Forest Park in the city of Dongtai, Jiangsu Province (32°51′-32°52′N, 120°48′-120°50′E, 5 m above sea level), and has a subtropical monsoon climate.The region is characterized by an average annual temperature of 15.6 °C, total annual precipitation of 1061.2 mm, and a sunshine duration of 2130.5 h.The soils at the study site are classified as fluvisols according to the World Reference Base (IUSS Working Group WRB 2015).The vegetation in the park mainly consists of deciduous broad-leaved forests and evergreen broad-leaved forests, with a forest coverage rate of 66.6% and a total area of 4156.93 ha.
This study included one wheat (Triticum aestivum)corn (Zea mays) rotation system with a > 30-year history, abbreviated herein as W, and four forest land types that were wheat-corn rotation systems prior to afforestation.More details are listed in Table 1.Specifically, the ginkgo (Ginkgo biloba) trees at the G, GW, and GM sites were planted in 2002, and the metasequoia (Metasequoia glyptostroboides) trees at the M and GM sites were planted in 2010.The harvest of most crop residues (mainly wheat and corn straw residues) at the W and GW sites occurs every year, leaving only small amounts of residue.

Field sampling
Three 10 m × 10 m plots were randomly selected within each land use pattern, and the stand edge areas were always avoided to prevent any interactions among adjacent land use patterns and eliminate any stand edge effects.Sampling took place in mid-May 2018, when the wheat in the W and GW systems was at the beginning of the grain-filling period.According to our previous study, the SOC content decreased greatly at a depth of 20 cm (Guo et al. 2018a).Thus, stratified sampling was conducted at 0-20 cm and 20-40 cm soil depths in the abovementioned plots following an "S"-type pattern.Each sample was mixed with 10 soil cores that were randomly collected from each soil depth with a cutting ring.The sample was immediately sieved through a 2-mm screen, and fine roots and litter debris were removed.A total of 30 samples (3 replicates × 2 depths × 5 land use patterns) were collected; a small portion of the samples were placed into a dry-ice box, and the rest were placed into an ice box.After transporting them to the lab, the samples in the dry-ice box were stored in a -81 °C ultralow-temperature refrigerator for 16S rRNA sequencing.The samples in the ice box were divided into two subsamples.One subsample was then stored in a 4 °C refrigerator for SOC fraction and soil enzyme activity determination, and the other subsample was air-dried indoors for further incubation and soil analysis.

Laboratory incubation
The incubation experiment performed herein included five land-use patterns and two soil depths with three replicates.Specifically, 15 g of sieved, air-dried soil was evenly placed in a 250-mL jar, deionized water was added to maintain 60% of the field water capacity, and three empty jars were used as controls.After preculturing for 7 days (with daily hydration during this period), a 60-day incubation experiment was conducted using an incubator  2, 3, 4, 5, 6, 7, 9, 11, 13, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 46, 50, 54, 58, and 60 days) during the incubation process.Then, excess 1 M BaCl 2 was added to precipitate CO 2− 3 , and the residual NaOH in the solution was titrated with 1 M HCl in the presence of excess BaCl 2 .Cumulative mineralization was expressed in mg CO 2 -C kg −1 dry soil (Li et al. 2022a).After each CO 2 release measurement, the jar was kept open for 2 h to ensure sufficient O 2 availability for soil microbial respiration.The SOC mineralization ratio was calculated as the cumulative CO 2 emissions divided by the SOC content (g CO 2 -C g −1 SOC), and the metabolic quotient (qCO 2 , mg CO 2 -C g −1 MBC h −1 ) was calculated according to Xiao et al. (2017): where ΣC i is the cumulative CO 2 emission (mg CO 2 -C kg −1 ), MBC is the soil MBC content (g kg −1 ), and T is the incubation time.

Soil analysis
Classical approaches were employed to determine the pH, bulk density, SOC, total nitrogen (TN), total phosphorus (TP), and total potassium (TK) contents of the soil samples.Briefly, the soil pH was determined using the potentiometric method (soil:water ratio of 1:2.5, w/v); the ring knife method was used to determine the bulk density; SOC was determined using the potassium dichromate oxidation external heating method; TN was determined using the semimicro Kjeldahl method; TP was determined through sodium hydroxide melting molybdenum antimony colorimetry; and TK was determined using the sodium hydroxide melting flame photometer method.The results of these indices and the soil bacterial community composition and diversity have been published in Guo et al. (2021).The analysis methods for microbial biomass carbon (MBC), dissolved organic carbon (DOC), and readily oxidizable organic carbon (ROC) were previously described in Guo et al. (2018a).
The soil urease, alkaline phosphatase, and catalase activities were determined based on the methods of Guan (1986).Briefly, the soil alkaline phosphatase activity was determined by the sodium benzene phosphate colorimetric method, soil urease activity was determined by the sodium diphenyl phosphate colorimetric method, and catalase activity was determined by the potassium permanganate titration method.

Kinetic modeling
A first-order kinetic model was used to calculate the key parameters of SOC mineralization, and the equation can be expressed as follows (Li et al. 2022a): where t is the incubation time, C t is the accumulated amount of SOC mineralized at time t, k is the mineralization rate constant, and C 0 is the potentially mineralizable SOC.In addition, the proportion of potentially mineralizable SOC to total SOC can be calculated with the following formula: C 0 /SOC.

Statistical analysis
SPSS software (version 20.0, Chicago, USA) was used to calculate the mean and standard deviations of the triplicate samples.Two-way analysis of variance (ANOVA) was used to analyze the effects of different land-use patterns and soil depths on the resulting SOC fractions, enzyme activity, and mineralization-relevant indicators.Before conducting ANOVA analyses, the distribution normality and variance uniformity were tested.Pearson's correlation analysis was performed to analyze the correlations between the cumulative SOC mineralization and soil SOC fractions, physicochemical properties, enzyme activity, and bacterial α-diversity indices.The first-order kinetics fitting curve was derived using Origin software (version 2018, Northampton, USA).

SOC and fractions and extracellular soil enzyme activity
Compared with the W system, LUCs increased the MBC by 15.41-77.02%and by 87.50-156.31% at the 0-20 cm and 20-40 cm soil depths, respectively.Only the GW system at the 0-20 cm soil depth and the M system at the 20-40 cm soil depth exhibited higher ROC contents than those in the W system, with an increase by 25.86% and 43.72%, respectively.The DOC contents in the G system were highest at both the 0-20 cm and 20-40 cm soil depths, reaching 261.16 and 283.72 mg•kg −1 , respectively.The two-way ANOVA results showed that LUCs had significant influences on the MBC, ROC, and DOC, while the soil depth significantly influenced the MBC and ROC contents (Fig. 1, P < 0.05).The lowest urease activity was found in the W system at the 0-20 cm soil depth, only 14.14-29.20% of other systems.The GW system had 114.66% and 10.15% higher alkaline phosphatase activity than the W system at the 0-20 cm and 20-40 cm soil depths, respectively.The two-way ANOVA results showed that LUCs exerted a significant influence on the alkaline phosphatase and catalase activity and that soil depth significantly influenced the urease, alkaline phosphatase and catalase activity (Fig. 2, P < 0.05).

Cumulative SOC mineralization
During this 60-day indoor incubation experiment, the cumulative SOC mineralization significantly varied among the five land-use patterns (Fig. 3).At the 0-20 cm depth, the cumulative mineralization followed the order of M > G > GM > W > GW, and the cumulative amount was highest in the M system (118.80mg CO 2 -C•kg −1 dry soil), which was 1.17 times that of W system and 1.20 times that of GW system.The cumulative SOC mineralization amount of W and GW systems was similar and significantly lower than those in other systems (P < 0.05).At the 20-40 cm depth, the cumulative SOC mineralization followed the order of G > W > GM > M > GW, and the cumulative amount under the G system was 78.08 mg CO 2 -C•kg −1 dry soil, which was 1.29 times that of the GW system.The cumulative SOC mineralization decreased with increasing soil depth, with the M system showing the greatest decrease of 55.26 mg of CO 2 -C•kg −1 dry soil and the W system showing the smallest decrease of 31.13 mg CO 2 -C•kg −1 dry soil.The Pearson correlation analysis results showed that the cumulative SOC mineralization was significantly positively correlated with ROC, SOC (R 2 = 0.764, P = 0.01; R 2 = 0.643, P = 0.045), MBC, and C/N (R 2 = 0.467, P = 0.174; Fig. 4).

SOC mineralization rate
All land use patterns showed rapid increases and great fluctuations in SOC mineralization rate during the initial incubation phase; maximum rates were observed at 15-18 d (for the 0-20 cm depth) and at 9 d (for the 20-40 cm depth), and after that the efflux gradually decreased and then stabilized (Fig. 5).At the 0-20 cm depth, the maximum rate measured under the M system was 13.47 mg CO 2 -C•kg −1 •d −1 , which was 94.86% and 121.18% higher than those obtained under the W and G systems, respectively.At the 20-40 cm depth, the maximum rate measured under the M system was 8.79 mg CO 2 -C•kg −1 •d −1 , which was 7.67% and 36.69%higher than those obtained under the W and G systems, respectively.The average SOC mineralization rates showed significant differences among the different land-use patterns and between the two soil depths (P < 0.05), as well as the SOC mineralization ratio and the microbial metabolic quotient (Fig. 6).The mineralization ratio was highest in the W system and lowest in the GW system, and the W system had the highest microbial metabolic quotient, reaching 41.42 mg CO 2 -C•g −1 MBC•h −1 .

SOC mineralization fitting parameters
The first-order kinetic equation of incubation time and cumulative SOC mineralization fitted well, with fitting coefficients (R 2 ) ranging from 0.931 to 0.986 (Table 2).

Relationships between cumulative SOC mineralization and soil properties, enzyme activity, and soil bacterial α-diversity
Relationships between the cumulative SOC mineralization and the soil properties, enzyme activity, and soil bacterial α-diversity were analyzed (Table 3).The soil pH and bulk density were significantly negatively correlated with the cumulative SOC mineralization (P < 0.05).The soil bacterial Shannon index was significantly positively correlated with cumulative SOC mineralization (P < 0.05).
Multiple stepwise linear regressions were used to further analyze the relationships between the cumulative SOC mineralization and soil properties, enzyme activity, and soil bacterial α-diversity indices (Table 4).The SOC mineralization potential was well explained by the Shannon index (68.9%,Model 1) and jointly explained by the Shannon index and OTUs (83.5%, Model 2).

Discussion
Overall, reducing the CO 2 expenditure (a lower SOC mineralization level) in the soil is always effective strategy for combating climate change (Montagnini and Nair 2004;Hontoria et al. 2016).LUCs always lead to changes in the physical, chemical and biological conditions of the soil, which may promote or inhibit the SOC release.Thus, exploring the SOC mineralization rules and their variations among land use patterns is critical for accurately predicting SOC dynamics and the C sequestration potential following LUCs.Generally, the results obtained in our study were partly contrary to our first hypothesis but they verified the second hypothesis.Exactly, different LUC patterns had diverse impacts on SOC mineralization The G, M, and GM systems had significantly higher cumulative SOC mineralization values than the other two systems in the 0-20 cm soil layer.GW had the highest SOC content and lowest cumulative SOC mineralization, indicating the small proportion of substrate available for microbial decomposition (Zhang et al. 2021a).Another explanation is that the wheat straws generally have higher C/N ratio and more recalcitrant C components, which hindered the litter decomposition and C mineralization (Guo et al. 2018b).The cumulative SOC mineralization at the 0-20 cm and 20-40 cm soil depths after 60-day incubation ranged from 99.08-118.80mg CO 2 -C•kg −1 dry soil and from 60.79-78.08mg CO 2 -C•kg −1 dry soil, respectively.The mineralization ratio was relatively low, likely because the incubated soils represented original soils without any addition of organic substrates, such as straw or biochar (Liu et al. 2019;Roy et al. 2022).Afforestation leads to changes of aboveground vegetation, root distribution and density, litter quality and quantity, and then affects the microbial community composition, thus resulting in SOC and fractions variations, all of these will affect the SOC mineralization patterns, especially in the topsoil layer (Tian et al. 2016;Chavarria et al. 2018;Guo et al. 2021).Likewise, Chia et al. (2017) reported that a greater labile SOC pool and larger soil microbial populations enhanced SOC mineralization, especially in the first 10 years of afforestation (Chen et al. 2016).The gradually increasing proportion of recalcitrant SOC (aromatic-like substances) maintains pace with the forest stand development, which may lead to a reduction in the SOC mineralization rate (Huang et al. 2019).Thus, the greater cumulative SOC mineralization in the M system might be ascribed to the higher substrate availability and the resultant increases of microbial biomass and activity in this system, which is conducive to microbial decomposition.In our study, there were significant differences in the cumulative SOC mineralization between the two soil depths, with lower values in the subsoil layer (Neculman et al. 2021).Liang et al. (2018) suggested that the lower cumulative SOC mineralization in subsoils may be due to delayed mineralization rather than the provisioning of greater protection through mineral binding or a lack of metabolic potential.In our study, cumulative SOC mineralization was found to be positively correlated with ROC, SOC, MBC, and C/N, indicating that labile organic C pools are an important factor determining the soil CO 2 emission (Liu et al. 2019;Lu et al. 2020;Zhang et al. 2021a).Despite the decreasing proportion of labile C pools following afforestation, the total SOC amount increased, leading to a lower mineralization ratio but increased cumulative SOC mineralization.Moreover, a poor correlation between the cumulative SOC mineralization and DOC was also observed in Entisols by Ma et al. (2019); their finding was partly explained by the relatively high proportion of recalcitrant compounds in the DOC.The microbial biomass and substrate availability are always the dominant factors of SOC mineralization following afforestation.Conversely, the increased microbial growth and extracellular enzyme activities will facilitate the degradation of SOC.As demonstrated by Barnett et al. (2022), LUCs have profound impacts on the soil C cycle path through the microbial food web, and these differences always linked to C mineralization performance.Microbial biomass influences C cycling by altering resource pathways between microbes with diverse life history strategies, in which low-biomass soils favor fastgrowing organisms that exhibit repeated cycles, while high-biomass soils favor microbes that adapt to efficient resource access and less dynamic growth cycles.These provide a framework for improving our understanding of microbially mediated terrestrial carbon cycling.
In line with previous studies, marked peaks and fluctuations in the mineralization rate were also observed in our incubation experiments, regardless of the soil properties or land use patterns (Zhou et al. 2013;Fanin and Bertrand 2016).The results showed that the maximum mineralization rate at the 0-20 cm soil depth was higher than that at the 20-40 cm soil depth, and great differences were observed among the different land-use patterns.Notably, greater litter and root inputs and higher microbial activities occurred in the topsoil layer, thus affecting soil particle aggregation and SOC quality and quantity and resulting in relatively high cumulative SOC mineralization in the topsoil layer (Chavarria et al. 2018;Nadal-Romero et al. 2016).The diverse SOC mineralization rates were directly or indirectly caused by long-term changes in biotic and abiotic factors and their interactions following LUCs.The relatively higher SOC mineralization rates in the afforested systems, especially at the subsoil depths, could be partially explained by the greater root biomass and deeper root distribution in these land use types (Kuzyakov and Xu 2013).Previous studies have shown that the main drivers of SOC changes in different tillage systems are the quantity and quality of organic inputs rather than any change in soil tillage practices, which, in fact, may not affect the specific SOC mineralization rates (Autret et al. 2020).The degradation of labile C contributed to the rapid increase in the SOC mineralization rate observed over the first 9 or 15 days at the 20-40 cm or 0-20 cm soil depths, and the subsequent reductions in the active C pool proportion led to subsequent declines in the SOC mineralization rate (Munda et al. 2018).Hence, soils with low available SOC contents reached the maximal SOC mineralization rates rapidly, while soils with high available SOC contents require a longer period to consume these substrates, which would result in a lag before peak C mineralization is achieved.If the process was further divided into two stages, the initial period (days 0-33) mainly relied on the colonization ability of microbial communities and active C pool mineralization, while the second phase was mainly constrained by the labile C availability in the soil (Luo et al. 2017;Tang et al. 2018).It can also be interpreted as the mineralization patterns consist of the first phase of readily available SOC decomposition and the second phase of secondary transformations of microbial products (Barnett et al. 2022).Growing evidence has shown that soil physicochemical properties and microbial communities strongly affect the SOC mineralization rate (You et al. 2014;Fabian et al. 2017).The qCO 2 level has been proven to be a significant indicator of the microbial assimilation efficiency; a relatively high qCO 2 indicates a low utilization efficiency of carbon sources.Here, the qCO 2 levels at both soil depths in the W system were higher than those in the other systems, demonstrating that more carbon sources are used by microorganisms for respiration  rather than sustaining microbial development in this system, thus limiting the development of microbial communities and inhibiting soil quality improvements and SOC sequestration (Huang et al. 2019).Accordingly, SOC mineralization was directly promoted by increased microbial C assimilation, and the relative abundances of genes associated with labile C decomposition may have also played important roles (Wang et al. 2023).
Additionally, the first-order kinetic equation employed herein distinctly showed the amount of C lost at our study site.Furthermore, the relatively low C 0 /SOC values reflected the low mineralization rate and high SOC sequestration capacity to a certain extent (Sarkar et al. 2021).The decreased C 0 /SOC in the plantations and agroforestry systems after LUCs indicated that the afforestation promoted SOC transformation but reduced soil C loss due to mineralization.In our study, the cumulative SOC mineralization was found to be significantly correlated to the soil pH, bulk density, and bacterial Shannon index.When relatively high pH values were measured in the alkaline range, they always induced a reduction in microbial activity and then limited the decomposition of plant-derived macromolecules (Yang et al. 2022).Furthermore, in agreement with our results, Hu et al. (2020) indicated that the soil pH, texture and microbial communities greatly contribute to the persistence of SOM, especially the portion controlled by microorganisms.Thus, the significant effects of soil physicochemical properties, soil substrates, and soil bacterial communities and their interactions on the cumulative SOC mineralization and mineralization rate emphasize the complex regulatory mechanisms associated with biotic and abiotic factors on the accumulation and turnover of soil C (You et al. 2016;Tang et al. 2018).The results of the multiple stepwise linear regression analysis further showed that the SOC mineralization potential could be well explained by the Shannon index (68.9%,Model 1) and jointly explained by the Shannon index and OTUs (83.5%, Model 2); these findings corroborated those of Li et al. (2022b), who reported that the SOC mineralization is mainly controlled by the bacterial community composition but may also be regulated by the organic C supply in the subsoil.The bacterial community composition and networks affect enzyme activity involved in the degradation of C compounds and SOC mineralization (Wang et al. 2021).Regardless of the analyzed system or soil layer, the available organic C is exhausted by the late incubation stage, leading to changes in the abundances of keystone species and thereby increasing the utilization of refractory organic matters (Blagodatskaya and Kuzyakov 2008).To obtain empirical evidence of the relationships between the microbial community and SOC mineralization, sequencing and functional microbial diversity analyses should be further studied, especially the relationships between abundance of copiotrophic/oligotrophic taxa and the mineralization patterns, and the microbial variations during the incubation (Wu et al. 2023).Furthermore, due to the important role of fungi in recalcitrant  SOC decomposition, further studies should be conducted to verify their assumed role (Dong et al. 2023).

Conclusion
Our results confirmed that both LUCs and soil depths significantly affected the SOC mineralization rate, while different LUC patterns had diverse impacts on SOC mineralization, and the cumulative mineralization in the metasequoia plantation and ginkgo plantation was significantly higher than those in the wheat-corn and gingko-wheat-corn systems.The greater cumulative SOC mineralization in the topsoil was probably due to the greater organic C supply and the diverse bacterial community structure.Through the kinetic modeling, we also found that the gingko-wheat-corn system had lower potentially mineralizable organic C and higher SOC stability, which indicated its greater C sequestration potential and the specific mechanism that needs to be thoroughly studied in situ.Moreover, SOC mineralization potential was well explained by the bacterial Shannon index and the OTUs, revealing the response of SOC mineralization to bacterial communities.Overall, the results showed that LUCs significantly changed the SOC mineralization characteristics and highlighted the important roles of microbial processes driving soil carbon cycle, which contributes to the fundamental understanding of SOC turnover regulation.Our findings provide a scientific basis for understanding the C budget in this area, but the mechanisms of such CO 2 emission and the growth strategies of microorganisms warrant further investigation.

Fig. 1
Fig. 1 Contents of soil organic carbon and fractions of different land-use patterns at 0-20 cm and 20-40 cm soil depths.G, W, M, GW and GM represent ginkgo plantations, wheat-corn systems, metasequoia plantations, ginkgo-wheat-corn systems and ginkgo-metasequoia systems, respectively (values are the means ± standard deviations, n = 3).Different lowercase letters indicate significant differences among different land use patterns at the same soil depth (P < 0.05) The potentially mineralizable SOC (C 0 ) values ranged from 62.70 to 148.04 mg•kg −1 .The M and G systems had the highest C 0 values at the 0-20 cm and 20-40 cm soil depths (148.04 mg•kg −1 and 118.46 mg•kg −1 , respectively), while the GW system had the lowest C 0 value at both soil depths.The topsoil layer always had higher C 0 values and lower mineralization rate constants (k) than the 20-40 cm depth.The dynamics of cumulative SOC mineralization following incubation time at the 0-20 cm and 20-40 cm soil depths under different land-use patterns and the corresponding kinetic simulations are shown in Fig. 7.

Fig. 2
Fig. 2 Soil urease (a), alkaline phosphatase (b), and catalase (c) activity of different land use patterns at 0-20 cm and 20-40 cm soil depths (values are the means ± standard deviations, n = 3).Different lowercase letters indicate significant differences among land use patterns at the P = 0.05 level

Fig. 3
Fig. 3 Cumulative SOC mineralization under different land-use patterns during 60-day incubation at 25 °C (values are the means ± standard deviations, n = 3).a Cumulative SOC mineralization at 0-20 cm soil depth; b cumulative SOC mineralization at 20-40 cm soil depth

Fig. 4
Fig. 4 Pearson correlations between the cumulative SOC mineralization and microbial biomass carbon (a), and readily oxidizable organic carbon (b), and dissolved organic carbon (c), and soil organic carbon (d), and ratio of soil organic carbon and total nitrogen (e)

Fig. 7
Fig. 7 Changes in SOC cumulative mineralization with incubation time and dynamic simulation based on the first-order kinetic equation at 0-20 cm and 20-40 cm soil depths under different land-use patterns

Table 1
Basic characteristics of the study sites maintained at a constant temperature of 25°C.A small beaker containing 15 mL 1 M NaOH was used to absorb the released CO 2 .The beakers were removed to determine the absorbed CO 2 at different intervals (at exactly 1,

Table 2
Kinetic parameters of organic carbon mineralization in the topsoil and subsoil under different land use patternsThe k is the mineralization rate constant; C 0 is the potentially mineralizable organic carbon content in the soil; and C 0 /SOC is the proportion of potentially mineralizable organic carbon to total organic carbon

Table 3
Pearson's correlation analysis results between the cumulative SOC mineralization and soil properties, enzyme activity, and soil bacterial α-diversity levelsThe bold values reflecte the significant correlation between the cumulative SOC mineralization and specific soil property (*0.01 < P < 0.05, **P < 0.01) BD indicates bulk density; TN indicates total nitrogen; TP indicates total phosphorus; TK indicates total potassium; the soil bacterial α-diversity, operational taxonomic units (OTUs), Shannon's diversity index, ACE, and Chao values are also listed

Table 4
Multiple stepwise regression analysis results