Abstract
Low moisture conditions result in substantially more soil inorganic carbon (SIC) than soil organic carbon (SOC) in drylands. However, whether and how changes in moisture affect the temperature response of SIC in drylands are poorly understood. Here, we report that the temperature sensitivity of SIC dissolution increases but that of SOC decomposition decreases with increasing natural aridity from 30 dryland sites along a 4,500 km aridity gradient in northern China. To directly test the effects of moisture changes alone, a soil moisture control experiment also revealed opposite moisture effects on the temperature sensitivities of SIC and SOC. Moreover, we found that the temperature sensitivity of SIC was primarily regulated by pH and base cations, whereas that of SOC was mainly regulated by physicochemical protection along the aridity gradient. Given the overall increases in aridity in a warming world, our findings highlight that drought may exacerbate dryland soil carbon loss from SIC under warming.
Similar content being viewed by others
Introduction
Drylands are regions where the aridity index is below 0.65 (ref. 1), and they cover approximately 41% of the terrestrial surface2. Global warming is predicted to increase the aridity of terrestrial ecosystems worldwide3, resulting in an 11–23% increase in the global area of drylands within the 21st century4. However, these ecosystems are considered fragile5 and vulnerable to aridity changes6,7. The presence of soil inorganic carbon (SIC) is primarily regulated by parent material8, and the low moisture conditions in drylands favor a higher ratio of SIC to soil organic carbon (SOC) given that plant inputs are low, with approximately 2–10 times more SIC storage than SOC in these ecosystems9. Global SIC is estimated to be approximately 1237 Pg C to a depth of 2 m (ref. 10), with as much as 95% of the SIC stored in drylands11. Historically, SIC is generally considered very stable12, and thus, previous studies on the temperature response of soil carbon (C) have focused entirely on SOC13,14. Similar to the biochemical reaction of SOC, increasing evidence has indicated that the SIC process is temperature-dependent15, and SIC, therefore may contribute substantially to warming-induced soil C losses8,16,17. A recent synthesis of 28 studies has shown that SIC-derived CO2 contributed 27% of total CO2 emissions from calcareous soils16. However, it is not yet clear how temperature changes affect this vast SIC pool in drylands, where more frequent and extreme heat waves are predicted during the 21st century7.
The general projected trend of global drylands is toward drying under climate change18,19. Although it is well known that drying decreases soil CO2 emissions by inhibiting both SOC decomposition and SIC dissolution16, how changes in moisture affect the temperature response of these two processes remains unexplored. Previous studies on soil moisture effects on the temperature sensitivity (Q10) of total CO2 emissions (Q10_total) have revealed inconsistent results; positive20, negative21, and no effect22 have all been reported. The causes of this controversy may partly stem from the confounding moisture effects on the Q10 of SOC-derived CO2 emissions (Q10_SOC) and that of SIC-derived CO2 emissions (Q10_SIC)16. A decrease in soil moisture may suppress Q10_SOC by decreasing substrate availability23,24. Soil moisture can also exert crucial controls on Q10_SOC by regulating physicochemical protection25 and microbial communities26. In comparison, the release of C during SIC dissolution is mainly a chemical process, which is described by the following equations:
Soil moisture changes can directly drive the CaCO3–CO2–\({{{\mbox{HCO}}}}_{3}^{-}\) equilibrium equations to promote or inhibit CaCO3 dissolution17,27, indicating that moisture effects on SIC dissolution are not linear or stagnant and that reprecipitation processes may also occur. For example, a decrease in soil moisture would inhibit CaCO3 dissolution16, and the observed net SIC-derived CO2 emissions can be especially low under low moisture conditions, as CO2 is also consumed during carbonate dissolution28. Soil moisture can also indirectly affect SIC dissolution by mediating soil pH and/or base cations (e.g., Ca2+ and Mg2+)16 that can shift the reactions represented in Eqs. (1) and (2). Accordingly, a drop in soil H+ owing to the enhancement of soil pH resulting from soil moisture decrease29 will lead to the reactions represented in Eqs. (1) and (2) to proceed to the left, leading to low SIC-derived CO2 emission rates. Given that a low CO2 emission rate is more sensitive to environmental changes (e.g., temperature), moisture decrease and/or pH or base cation increase may enhance the temperature response of SIC dissolution.
Moreover, as SOC decomposition affects the soil CO2 concentration30, factors (e.g., physicochemical protection25 and substrate24) that affect SOC decomposition may also mediate SIC dissolution and its temperature response. An increase in SOC decomposition will lead to more active CO2 sources in soil for HCO3− production, which will restrict CaCO3 dissolution. Consequently, soil physicochemical protection and substrate can regulate SIC processes by mediating SOC decomposition. Although mineral protection of Ca bridges and/or Fe oxides has been shown to largely inhibit SOC decomposition and its temperature sensitivity14,31, its effects on Q10_SIC remain unknown. These direct and indirect moisture effects are not independent but coexist temporally and spatially16, while the main drivers and their differences in regulating Q10_SOC and Q10_SIC are poorly understood. Until now, no attempt has been made to examine moisture effects on Q10_SOC and Q10_SIC across a large moisture gradient in drylands. This knowledge gap urgently needs to be filled since future climate change is predicted to largely affect dryland moisture regimes18 and, consequently, the climate–C cycle feedbacks in these water-limited ecosystems.
Here, we hypothesized that (i) Q10_SOC decreased, but Q10_SIC increased with decreasing moisture content, and (ii) Q10_SOC was mainly regulated by physicochemical protection, while Q10_SIC was primarily regulated by chemical properties (e.g., pH and cation exchange capacity (CEC)). To test these hypotheses, we conducted two experiments (Fig. 1): a natural aridity gradient and a moisture control treatment. In the first experiment, soil moisture regime differences were evaluated by sampling soils from 30 sites across a wide aridity index (the ratio of precipitation to potential evapotranspiration, ranging from 0.04 to 0.59) along an approximately 4500 km east–west transect in the drylands of northern China; in this experiment, Q10_SOC and Q10_SIC were determined with field moisture conditions. Considerable differences in Q10 and its controls were expected to exist throughout the soil profile32,33 owing to the large differences in soil biotic and abiotic factors34,35. To test whether moisture effects on Q10_SOC and Q10_SIC persist among different soil depths, soils from the topsoil (0–10 cm) and subsoil (35–50 cm) were collected at each site. In the second experiment, to directly test the effects of only moisture changes on Q10_SOC and Q10_SIC, we conducted a moisture control experiment by incubating soils under different moisture conditions of 20%, 40%, and 60% water holding capacity (WHC). To determine the main drivers and their differences associated with variations in Q10_SOC and Q10_SIC that were determined under field moisture conditions along the aridity gradient, we analyzed various potential factors related to climate (mean annual temperature (MAT) and aridity index), physical (SOC stored in particulate organic matter (OC-POM) and mineral-associated organic matter (OC-MAOM) fractions, and SOC associated with Ca bridges (OC-Ca) and Fe oxides (OC-Fe)), chemical (pH, CEC, Ca2+, and Mg2+) and substrate (quantity, quality and availability) properties.
Results
Moisture effects on SOC- and SIC-derived CO2 emissions and their Q 10 values
Isotopic measurement of 13C content is considered an effective approach for distinguishing different CO2 sources16; SOC is less enriched with heavier 13C than SIC, and thus, the δ13C value of SOC-derived CO2 differs substantially from that of SIC-derived CO215,17. In this study, natural isotope technology was adopted to separate total soil CO2 emissions into SOC- and SIC-derived sources. Across the aridity gradient in drylands, the average contribution of SIC-derived CO2 to total CO2 emissions was 7.2% and 11.1% in the topsoil and subsoil, respectively, at 20 °C (Supplementary Fig. 1). SOC- and SIC-derived CO2 emissions differed as a function of soil moisture, showing that they both decreased with increasing aridity (Supplementary Fig. 2). Similar results were observed from the moisture control experiment, with lower emission rates under lower moisture contents (Supplementary Fig. 3).
To reveal moisture effects on the temperature response of SOC and SIC, we first evaluated changes in Q10_SOC and Q10_SIC along the aridity gradient; to do this, soils were incubated under field moisture conditions. Opposing patterns of Q10_SOC and Q10_SIC were observed in response to aridity changes, showing that in both soil layers, Q10_SOC decreased significantly with drying (that is, decreasing aridity index), whereas Q10_SIC increased significantly with drying (P < 0.001, Fig. 2). This was further verified by the moisture control experiment; given the predicted overall aridity increase in drylands in a warmer world, soils were incubated under different moisture contents of 20%, 40%, and 60% WHC. The moisture control experiment also showed that Q10_SOC was lower under lower experimental moisture conditions, but the opposite was true for Q10_SIC (P < 0.01, Fig. 3).
Factors regulating the Q 10_SOC and Q 10_SIC along the aridity gradient
We next explored factors that regulate the variations in the Q10 of SOC- and SIC-derived CO2 emissions along the aridity gradient. Potential influencing factors were determined, comprising groups of factors related to climate (i.e., MAT and aridity index), physical (i.e., SOC stored in the POM and MAOM fractions and SOC associated with Ca oxides and Fe bridges), chemical (i.e., pH, CEC, Ca2+ and Mg2+), substrate (i.e., substrate quantity of SOC and SIC contents and substrate availability of C availability index (CAI), and substrate quality of SOC decomposability (DSOC)). Correlation analysis revealed that Q10_SOC was linked to climate, soil physical, and substrate at both soil depths (Fig. 4). Specifically, Q10_SOC was positively correlated with the aridity index, OC-POM, and CAI but negatively correlated with OC-MAOM, OC-Fe, OC-Ca, and DSOC (Fig. 4). Q10_SIC was linked to climate and soil chemical properties at both soil depths, showing that Q10_SIC was negatively correlated with the aridity index but positively correlated with pH, Ca2+, and Mg2+ (Fig. 4). A similar phenomenon was observed when Q10_SOC and Q10_SIC were determined under common moisture conditions (Supplementary Fig. 4).
A structural equation modeling (SEM) was further constructed to assess the direct and indirect effects of these factors (i.e., climate, physical, chemical, and substrate properties) on the Q10 of SOC- and SIC-derived CO2 (Fig. 5). The results showed that among all the groups of factors tested, climate, soil physical and substrate had direct effects on Q10_SOC (Fig. 5a), with greater total effects of climate and soil physical variables than other factors (Fig. 5c). For Q10_SIC, climate and soil chemical properties had significant direct effects (P < 0.05, Fig. 5b), and they had greater total effects than other factors (Fig. 5d). Similar results were observed in the subsoil (Supplementary Fig. 5).
Discussion
On the basis of large-scale sampling and isotope approaches, over the short-term time scale, we observed substantial SIC contributions to soil total CO2 emissions. In long-term monitoring studies at national or even global scales, some recent studies have also observed SIC losses36,37. Given that SIC accumulation usually takes substantially more time than SOC38, SIC losses are thus more impactful for the soil C reservoirs than that of SOC in these water-limited ecosystems. Moreover, we observed that absolute SOC- and SIC-derived CO2 emissions both decreased, but the contribution of SIC to total CO2 emissions increased with decreasing moisture along the aridity gradient; this might be because the microbial process of SOC decomposition is more sensitive to moisture changes than the chemical process of SIC dissolution17,39. Our results thus indicate that the dissolution of SIC is more important than previously thought in regulating atmospheric CO2 concentrations8, and if future climate change accelerates aridity in drylands18, the contribution of SIC-derived CO2 to total CO2 emissions may become even more substantial.
Temperature sensitivity represents a key parameter in many previously developed biogeochemical models that have simulated C emissions40, and small inaccuracies in this parameter can result in large errors41. However, total CO2 emissions were typically measured in the majority of previous studies on Q10, and the calculated value was considered the Q10 of SOC decomposition13,14,32,42. For acidic soils, measurements of total CO2 emissions are sufficient for evaluating the Q10 of SOC decomposition33; for calcareous soils, however, this would bias our understanding, as Q10_total was higher than Q10_SOC. Although the absolute rate of SIC-derived CO2 is low compared to that of SOC-derived CO2 emissions, the high-temperature sensitivity of SIC and the vast SIC stock in drylands11 can also harbor great potential for regulating climate–C cycle feedbacks in drylands. Our results provide further evidence that moisture has opposite effects on the temperature response of SOC decomposition and SIC dissolution, which is a crucial step forward in gaining an understanding of climate–C cycle feedbacks in drylands. The general projected global trend for drylands predicts drying18, and our findings of positive Q10_SOC–moisture but negative Q10_SIC–moisture relationships suggest that drought may exacerbate warming-induced soil C loss from inorganic C in drylands.
Total SOC- and SIC-derived CO2 emissions are roughly estimated at 20.4 Pg C year−1 in drylands (assuming that soil respiration in drylands accounts for 38.6% of global soil respiration43, 60% of which is from the heterotrophic component44), and SIC-derived CO2 contributes to approximately 27.0% of total CO2 emissions16. Using Q10_SOC to represent Q10_SIC would underestimate warming-induced SIC-derived CO2 emissions (3.2 Pg C year–1) by approximately 25.6% compared to that estimated (4.3 Pg C year−1) using the higher Q10_SIC under 4 °C of warming. Moreover, considering the different responses of Q10 to moisture changes (Q10_SOC decreases by 0.47 and Q10_SIC increases by 0.39 per 0.1 decrease in the aridity index; Fig. 2), the net increase in soil C due to the lower Q10_SOC (1.5 Pg C year−1) would be offset by approximately 26.7% due to the higher Q10_SIC (0.4 Pg C year−1) under 0.1 decreases in the aridity index. Consequently, although these are rough estimates, they highlight the importance of separately representing Q10_SOC and Q10_SIC and their different responses to moisture changes to improve projections of climate–C cycle feedback in drylands.
This study has further identified differential mechanisms of Q10_SOC and Q10_SIC along the aridity gradient. Soil physicochemical protection primarily regulated the temperature sensitivity of SOC decomposition. The effect of soil physicochemical protection on Q10_SOC could be due to the constraints of mineral protection on SOC availability and/or enzyme activity. Specifically, MAOM fractions can restrict oxygen diffusion and lead to the compartmentalization of organic C substrates from enzymes, and these processes can be enhanced by Ca bridges and/or Fe oxides31. In addition, Ca bridges and/or Fe oxides can constrain substrate availability by forming inner- and outer-sphere cation bridging between the negatively charged phyllosilicates and SOC45. Either of these processes may suppress the temperature response of SOC decomposition24,34. Consistent with this speculation, Q10_SOC was positively correlated with OC-POM but negatively associated with OC-MAOM, OC-Fe, and OC-Ca along the aridity gradient (Fig. 4). In addition to physicochemical protection, aridity-induced changes in substrate also exerted roles in regulating Q10_SOC. Consistent with previous studies on moist soils46,47, we observed that low substrate quality was associated with high Q10_SOC, as indicated by the negative relationship of Q10_SOC with DSOC (Fig. 4), suggesting that the C quality-temperature hypothesis13 is also applicable in these water-limited ecosystems.
However, Q10_SIC was primarily regulated by aridity-induced changes in the soil chemical properties of soil pH and base cations along the aridity gradient, showing that higher pH and CEC enhance the temperature response of SIC. This is because a higher pH and/or base cation (e.g., Ca2+ and Mg2+) concentration may enhance the reverse reactions represented in Eqs. (1) and (2) toward the absorption of CO2 into the soil solution. Given that a low CO2 emission rate might be more sensitive to temperature changes, high pH and/or base cations can thus enhance Q10_SIC. Soil pH largely determines the stability of SIC48, and we observed a positive correlation between pH and Q10_SIC but not Q10_SOC; this demonstrates that the usually observed positive pH–Q10 relationship that measures total CO2 emissions33,49 may result from pH-induced changes to the temperature responses of SIC dissolution but not to microbial processes of SOC decomposition. The generally predicted drying may increase soil pH in drylands and thus further enhance the temperature response of SIC. Although an overall increase in aridity is predicted in drylands in a warmer world18, some dryland regions are becoming increasingly prone to flooding50, leading to increases in pedogenic carbonate accumulation51; however, flooding may also result in losses of both SOC and SIC through soil erosion in these areas52.
In conclusion, our findings provide evidence for opposing moisture effects on the temperature response of SOC and SIC in drylands, suggesting that drying will further enhance the temperature response of SIC but weaken that of SOC (Fig. 6). This may partly explain the substantial loss of SIC pools over the past few decades on the continental scale53 under conditions of drying. Additionally, we identified differential mechanisms regulating the temperature responses of SOC and SIC (Fig. 6). As drought is expected to enhance soil alkalinity54, warming-induced SIC losses may be gradually enhanced by ongoing drought. In contrast, global changes in widespread nitrogen (N) deposition and/or acid rain are expected to promote soil acidification55, possibly weakening warming-induced SIC losses. Although some studies have shown that soil C in drylands is resistant to N deposition56,57, a recent study concluded that global N fertilization results in releases of 7.5 × 1012 g C year−1 from carbonates37, and this value may be even underestimated8. Nevertheless, our finding of the positive pH–Q10_SIC relationship suggests that SIC losses attributed to N deposition or acid rain may be weakened in a drying world. Therefore, to gain a better understanding of climate–C cycle feedbacks at an ecosystem level in drylands, future work should assess the potential effects of multiple global change factors on soil physicochemical (e.g., physical protection and pH) and biological (e.g., microbial community composition and functions) conditions and consequently their linkages with soil organic and inorganic C cycling.
Methods
Study area and field sampling
Soils were collected from 30 sites along an approximately 4500 km east–west transect with a longitudinal gradient of 81.02–123.53°E in northern China (Fig. 1). The climate is predominantly arid and semiarid continental; this transect covers the majority of the drylands in China and is the main reservoir of SIC in China38. The MAT and MAP ranged from −1.2 to 10.0 °C and from 46 to 486 mm, respectively, resulting in a wide range of aridity, with aridity indices (the ratio of precipitation to potential evapotranspiration) ranging from 0.04 to 0.59 at these sites; no significant relationship between MAT and aridity index was observed along the aridity gradient (P = 0.128, Supplementary Fig. 6). Soil properties for these sites were also greatly different; for example, soil pH and SIC content were significantly and negatively correlated with aridity index (P < 0.05, Supplementary Fig. 7). Therefore, this transect provided an ideal natural platform for studying moisture effects on the temperature sensitivity of SOC and SIC in drylands.
For soil sampling, three random locations (>20 m apart from each other) were chosen at each site. Because of the differences in environmental constraints, soil physical and chemical properties, and substrate and microbial properties among soil depths34,35, there were likely to be considerable differences in Q10 and its controlling factors throughout the soil profile32. Thus, soils from different depths of the topsoil (0–10 cm) and subsoil (35–50 cm) were collected. Soil samples were then passed through a 2-mm sieve, and soils from the three random locations were gently mixed to produce a homogeneous composite sample for each depth at each site. The composite sample was divided into three subsamples: one part was air-dried and processed for measurements of soil physical, chemical, and some substrate properties; another part was stored at 4 °C for soil incubations; the third part was stored at −20 °C for some microbial property analyses.
Temperature sensitivity assessments
The Q10 of SOC- and SIC-derived CO2 emissions for the 60 soils (30 sites × 2 soil depths) was determined using a laboratory incubation experiment with a short‐term dynamic temperature ramping method58, which could minimize the different depletion of soil C pools59,60 and microbial adaptation61 to different temperatures relative to separate soil incubations at different but constant temperatures. To test moisture effects on the Q10_SOC and Q10_SIC, we used two experiments: 1) the soils were incubated under field moisture conditions, and the relationship of Q10 with the aridity index was tested; and 2) the soils were incubated under some common moisture content conditions of 20%, 40%, and 60% WHC, and the difference in Q10 among different moisture conditions was tested. For the first experiment, the soils were quickly sieved to 2 mm at 4 °C after they were transported to the laboratory to minimize moisture changes and were then incubated to estimate Q10. For the second experiment, soils sieved to 2 mm were adjusted to different moisture conditions of 20%, 40%, and 60% WHC by adding deionized water (with a pH of approximately 7.3) and were then incubated. For soil incubation for both experiments, 50 g of dry‐weight fresh soil, with four experimental replicates, was maintained under field moisture conditions for the first experiment or adjusted to different moisture conditions (20%, 40%, and 60% WHC) for the second experiment and incubated in 250 ml jars. After a 2-week preincubation period at 20 °C to minimize disturbances from soil packing and rewatering, the jars with soils were incubated at 5–30 °C with a stepwise increase of 5 °C to perform the dynamic temperature ramping incubation32. After the soils were adjusted to the new target temperature and equilibrated for 3 h, the soil containers were sealed and flushed with CO2-free air following previous studies39,62; 15 ml of headspace gas was removed for the initial CO2 analysis, and 15 ml of CO2-free air was immediately injected into the jars to allow them to equalize to atmospheric pressure. After incubation for 4–72 h (depending on the target temperature), another 15 ml gas sample was collected. The CO2 concentration and δ13C value of these gas samples were analyzed using a gas isotope analyzer (G2201-20i, Picarro, USA).
SOC- and SIC-derived CO2 was determined using a two-end-member mixing model:
where fSIC is the contribution of SIC-derived CO2 to total CO2 emissions and δ13CCO2_total, δ13CCO2_SOC and δ13CCO2_SIC are the δ13C values for total CO2 release, SOC-derived CO2 and SIC-derived CO2, respectively. We assumed that the δ13C value was the same for SIC and SIC-derived CO2 and for SOC and SOC-derived CO215,16,17. The δ13C value of CO2 was then corrected because of differences in the fractionation of δ13C at different temperatures63. Moreover, the δ13C value of CO2 released was 4.4‰ lower than that of CO2 in the soil because of the fractionation induced by molecular diffusion64.
A prior test showed that SOC- and SIC-derived CO2 emissions increased exponentially with increases in incubation temperature (the fitting coefficient R2 > 0.95 for SOC-derived CO2 emissions and R2 > 0.85 for SIC-derived CO2 emissions). The temperature sensitivity of SOC- and SIC-derived CO2 emissions was then calculated as follows:
where R is the rate of SOC-derived and SIC-derived CO2 (μg C g−1 soil h−1), T is the incubation temperature (°C), and B and k are model fitting parameters.
Although the soil incubation experiment allowed us to reveal a general pattern and the mechanisms of moisture effects on Q10_SIC and Q10_SOC across large scales, we were aware that there were some possible influences, such as CO2-free air flushing and soil sieving. To test the possible effects of CO2-free air flushing on Q10_SOC and Q10_SIC estimation, we conducted a supplementary experiment using six representative soils (see Supplementary Text). The results showed that CO2-free air flushing significantly underestimated SOC- and SIC-derived CO2 emissions (P < 0.05, Supplementary Fig. 8); this might be because some of the SOC- and SIC-derived CO2 might go into the pore space of soil, leading to an underestimation of SOC- and SIC-derived CO2 emissions. However, CO2-free air flushing had no significant effects on Q10_SIC and Q10_SOC (P > 0.05, Supplementary Fig. 8); this might be because Q10 is a ratio for the temperature response of CO2 emissions, resulting in limited effects of CO2-free air flushing on Q10 values. Moreover, we conducted another supplementary experiment using intact and sieved soils to test the effects of soil sieving on Q10_SOC and Q10_SIC (see Supplementary Text). The results showed that sieving did not exert significant effects on SOC- and SIC-derived CO2 emissions and their Q10 values (P > 0.05, Supplementary Fig. 9). This might be because the soils in drylands are usually sandy, resulting in limited effects from sieving. A recent study has also shown that soil sieving had no substantial effects on the Q10 of total CO2 emissions65.
Climate variables
MAT and MAP data were obtained from the WorldClim (https://www.worldclim.org/). Aridity indices were obtained from the Global Aridity Index and Potential Evapotranspiration Climate database (https://cgiarcsi.community/).
Soil property analyses
To explore the factors regulating the temperature response of SIC- and SOC-derived CO2 along the aridity gradient, soil physical, chemical, and substrate properties were determined for the 60 soil samples (30 sites × 2 soil depths) collected across the natural aridity gradient.
Soil physical properties
The physical properties of POM and MAOM and the SOC associated with Ca bridges (OC-Ca) and Fe oxides (OC-Fe) were determined. A fractionation technique was adopted to estimate the SOC stored in the POM and MAOM fractions. Air-dried soil was separated into light and high-density fractions with sodium polytungstate solution (1.60 g cm−3)66; the high-density fractions were then wet sieved to collect POM (>53 μm) and MAOM (<53 μm)67. Moreover, to determine OC-Ca and OC-Fe contents68, the high-density fractions were extracted using 0.5 M Na2SO4 to release OC-Ca; the remaining residues were then extracted with citrate–bicarbonate–dithionite and sodium chloride for the treatment and control groups, respectively, and the differences in SOC content between the two groups were treated as the OC-Fe measurements. The SOC contents in these fractionations were ultimately determined using an elemental analyzer (Multi EA 4000, Analytik Jena, Germany) after inorganic C was removed with 1 M HCl.
Soil chemical properties
The chemical properties of pH, CEC, Ca2+, and Mg2+ were determined. Soil pH was determined using a pH electrode (Seven Excellence S479-uMix, Mettler-Toledo, Switzerland) in a 1:2.5 soil:water suspension. Ca2+ and Mg2+ contents were measured by using inductively coupled plasma‒optical emission spectrometry69. CEC was determined by using a microplate reader (Synergy 2, BioTek, USA) following extraction using [Co(NH3)6]Cl3 (ref. 70).
Substrate properties
The substrate quantity included the SOC and SIC contents. The SOC content was analyzed using an elemental analyzer (Multi EA 4000, Analytik Jena, Germany) after inorganic C was removed with 1 M HCl. SIC was determined by a pressure calcimeter method71. Briefly, 0.5 g of soil was mixed with 2 mL 6 M HCl and reacted in a closed reaction vessel. Two hours later, the pressure was determined using the pressure transducer and voltage meter, and then the carbonate concentration was calculated using a calibration curve, which was obtained in the same way using known quantities of CaCO3. The SIC content was finally determined by multiplying by a coefficient of 0.12, which is the mass proportion of C in calcium carbonate. The substrate availability was indicated by the C availability index, which was defined as the ratio of the basal respiration to the substrate-induced respiration24. A 60 g L−1 glucose solution was added for the substrate-induced respiration, and deionized water was added in the same manner for the basal respiration rate; respiration rates at 20 °C were estimated for the added glucose and ambient-substrate treatments. The substrate quality was indicated by SOC decomposability (DSOC), that is the SOC decomposition rate per unit of SOC content per hour. DSOC was calculated by the ratio of B (the parameter from Eq. (4)) to SOC content.
Statistical analyses
A paired-sample t-test was applied to examine the differences in Q10_SOC or Q10_SIC among different experimental moisture treatments (20%, 40%, and 60% WHC). Correlation analysis was conducted to test the correlations of Q10_SOC and Q10_SIC with each variable tested. Additionally, SEM was conducted to partition the direct and indirect effects of climate, physical, chemical, and substrate properties on Q10_SOC and Q10_SIC. Because the variables within each of these groups were closely correlated, principal component analyses were conducted to create a multivariate functional index prior to SEM analyses14,72. The first component was used for the combined group properties in the SEM analysis. The maximum likelihood estimation method was used to fit the data in the SEM analysis. The selection of the final model was based on the p-value, χ2 test, root-mean-squared error of approximation, and goodness-of-fit index73. The SEM was conducted using AMOS 21.0 software (Amos Development Corporation, Chicago, IL).
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
Supplementary Information is available online. The Q10 value and soil properties data are available at https://doi.org/10.5281/zenodo.10370941.
References
UNCCD. United Nations Convention to Combat Desertification—Global Land Outlook (UNCCD, 2017).
Reynolds, J. F. et al. Global desertification: building a science for dryland development. Science 316, 847–851 (2007).
Fu, Q. & Feng, S. Responses of terrestrial aridity to global warming. J. Geophys. Res. Atmos. 119, 7863–7875 (2014).
Huang, J., Yu, H., Guan, X., Wang, G. & Guo, R. Accelerated dryland expansion under climate change. Nat. Clim. Change 6, 166–171 (2016).
Li, C. et al. Drivers and impacts of changes in China’s drylands. Nat. Rev. Earth Environ. 2, 858–873 (2021).
Thornton, P. K., Ericksen, P. J., Herrero, M. & Challinor, A. J. Climate variability and vulnerability to climate change: a review. Glob. Change Biol. 20, 3313–3328 (2014).
Middleton, N. & Sternberg, T. Climate hazards in drylands: a review. Earth Sci. Rev. 126, 48–57 (2013).
Zamanian, K. & Kuzyakov, Y. Contribution of soil inorganic carbon to atmospheric CO2: more important than previously thought. Glob. Change Biol. 25, E1–E3 (2019).
Ferdush, J. & Paul, V. A review on the possible factors influencing soil inorganic carbon under elevated CO2. CATENA 204, 105434 (2021).
Plaza, C. et al. Soil resources and element stocks in drylands to face global issues. Sci. Rep. 8, 13788 (2018).
Lal, R. Carbon cycling in global drylands. Curr. Clim. Change Rep. 5, 221–232 (2019).
Raza, S. et al. Inorganic carbon losses by soil acidification jeopardize global efforts on carbon sequestration and climate change mitigation. J. Clean. Prod. 315, 128036 (2021).
Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173 (2006).
Qin, S. et al. Temperature sensitivity of permafrost carbon release mediated by mineral and microbial properties. Sci. Adv. 7, eabe3596 (2021).
Chevallier, T., Cournac, L., Hamdi, S., Gallali, T. & Bernoux, M. Temperature dependence of CO2 emissions rates and isotopic signature from a calcareous soil. J. Arid Environ. 135, 132–139 (2016).
Sun, Z., Meng, F. & Zhu, B. Influencing factors and partitioning methods of carbonate contribution to CO2 emissions from calcareous soils. Soil Ecol. Lett. 5, 6–20 (2022).
Inglima, I. et al. Precipitation pulses enhance respiration of Mediterranean ecosystems: the balance between organic and inorganic components of increased soil CO2 efflux. Glob. Change Biol. 15, 1289–1301 (2009).
Dai, A. Increasing drought under global warming in observations and models. Nat. Clim. Change 3, 52–58 (2013).
Huang, J. et al. Global semi-arid climate change over last 60 years. Clim. Dyn. 46, 1131–1150 (2016).
Craine, J. M. & Gelderman, T. M. Soil moisture controls on temperature sensitivity of soil organic carbon decomposition for a mesic grassland. Soil Biol. Biochem. 43, 455–457 (2011).
Meyer, N., Welp, G. & Amelung, W. The temperature sensitivity (Q10) of soil respiration: controlling factors and spatial prediction at regional scale based on environmental soil classes. Glob. Biogeochem. Cycles 32, 306–323 (2018).
Bradley-Cook, J. I., Petrenko, C. L., Friedland, A. J. & Virginia, R. A. Temperature sensitivity of mineral soil carbon decomposition in shrub and graminoid tundra, west Greenland. Clim. Change Responses 3, 1–15 (2016).
Davidson, E. A., Janssens, I. A. & Luo, Y. On the variability of respiration in terrestrial ecosystems: moving beyond Q10. Glob. Change Biol. 12, 154–164 (2006).
Gershenson, A., Bader, N. E. & Cheng, W. Effects of substrate availability on the temperature sensitivity of soil organic matter decomposition. Glob. Change Biol. 15, 176–183 (2009).
Huang, W., Ye, C., Hockaday, W. C. & Hall, S. J. Trade‐offs in soil carbon protection mechanisms under aerobic and anaerobic conditions. Glob. Change Biol. 26, 3726–3737 (2020).
Wang, G. et al. Soil moisture drives microbial controls on carbon decomposition in two subtropical forests. Soil Biol. Biochem. 130, 185–194 (2019).
Lardner, T., George, S. & Tibbett, M. Interacting controls on innate sources of CO2 efflux from a calcareous arid zone soil under experimental acidification and wetting. J. Arid Environ. 122, 117–123 (2015).
Liu, Z., Dreybrodt, W. & Liu, H. Atmospheric CO2 sink: silicate weathering or carbonate weathering? Appl. Geochem. 26, S292–S294 (2011).
Wang, C. et al. Aridity threshold in controlling ecosystem nitrogen cycling in arid and semi-arid grasslands. Nat. Commun. 5, 4799 (2014).
Fierer, N., Allen, A. S., Schimel, J. P. & Holden, P. A. Controls on microbial CO2 production: a comparison of surface and subsurface soil horizons. Glob. Change Biol. 9, 1322–1332 (2003).
Singh, M. et al. Stabilization of soil organic carbon as influenced by clay mineralogy. Adv. Agron. 148, 33–84 (2018).
Li, J. et al. Rising temperature may trigger deep soil carbon loss across forest ecosystems. Adv. Sci. 7, 2001242 (2020).
Li, J. et al. Depth dependence of soil carbon temperature sensitivity across Tibetan permafrost regions. Soil Biol. Biochem. 126, 82–90 (2018).
Qin, S. et al. Temperature sensitivity of SOM decomposition governed by aggregate protection and microbial communities. Sci. Adv. 5, eaau1218 (2019).
Fierer, N., Schimel, J. P. & Holden, P. A. Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 35, 167–176 (2003).
Zamanian, K., Zhou, J. & Kuzyakov, Y. Soil carbonates: the unaccounted, irrecoverable carbon source. Geoderma 384, 114817 (2021).
Zamanian, K., Zarebanadkouki, M. & Kuzyakov, Y. Nitrogen fertilization raises CO2 efflux from inorganic carbon: a global assessment. Glob. Change Biol. 24, 2810–2817 (2018).
Mi, N. et al. Soil inorganic carbon storage pattern in China. Glob. Change Biol. 14, 2380–2387 (2008).
Schindlbacher, A. et al. Contribution of carbonate weathering to the CO2 efflux from temperate forest soils. Biogeochemistry 124, 273–290 (2015).
Todd-Brown, K. E. et al. Causes of variation in soil carbon simulations from CMIP5 Earth system models and comparison with observations. Biogeosciences 10, 1717–1736 (2013).
Zhou, T., Shi, P., Hui, D. & Luo, Y. Global pattern of temperature sensitivity of soil heterotrophic respiration (Q10) and its implications for carbon‐climate feedback. J. Geophys. Res. Biogeosci. 114, G02016 (2009).
Li, J. et al. Biogeographic variation in temperature sensitivity of decomposition in forest soils. Glob. Change Biol. 26, 1873–1885 (2020).
Warner, D. L., Bond-Lamberty, B., Jian, J., Stell, E. & Vargas, R. Spatial predictions and associated uncertainty of annual soil respiration at the global scale. Glob. Biogeochem. Cycles 33, 1733–1745 (2019).
Throop, H. L., Seely, M. K. & Marufu, V. J., Summer Drylands Program, P. Multiple scales of spatial heterogeneity control soil respiration responses to precipitation across a dryland rainfall gradient. Plant Soil 453, 423–443 (2020).
Rowley, M. C., Grand, S. & Verrecchia, É. P. Calcium-mediated stabilisation of soil organic carbon. Biogeochemistry 137, 27–49 (2018).
Li, H. et al. Temperature sensitivity of SOM decomposition is linked with a K‐selected microbial community. Glob. Change Biol. 27, 2763–2779 (2021).
Craine, J. M., Fierer, N. & McLauchlan, K. K. Widespread coupling between the rate and temperature sensitivity of organic matter decay. Nat. Geosci. 3, 854–857 (2010).
Raza, S. et al. Dramatic loss of inorganic carbon by nitrogen‐induced soil acidification in Chinese croplands. Glob. Change Biol. 26, 3738–3751 (2020).
Li, J. et al. Carbon quality mediates the temperature sensitivity of soil organic carbon decomposition in managed ecosystems. Agric. Ecosyst. Environ. 250, 44–50 (2017).
Farquharson, F. A. K., Meigh, J. R. & Sutcliffe, J. V. Regional flood frequency analysis in arid and semi-arid areas. J. Hydrol. 138, 487–501 (1992).
Ortiz, A. C. et al. Dryland irrigation increases accumulation rates of pedogenic carbonate and releases soil abiotic CO2. Sci. Rep. 12, 464 (2022).
Xu, T. et al. Grassland degradation with saline-alkaline reduces more soil inorganic carbon than soil organic carbon storage. Ecol. Indic. 131, 108194 (2021).
Song, X. et al. Significant loss of soil inorganic carbon at the continental scale. Natl Sci. Rev. 9, nwab120 (2022).
Jiao, F., Shi, X., Han, F. & Yuan, Z. Increasing aridity, temperature and soil pH induce soil C-N-P imbalance in grasslands. Sci. Rep. 6, 19601 (2016).
Tian, D. & Niu, S. A global analysis of soil acidification caused by nitrogen addition. Environ. Res. Lett. 10, 024019 (2015).
Osborne, B. B. et al. Biogeochemical and ecosystem properties in three adjacent semi-arid grasslands are resistant to nitrogen deposition but sensitive to edaphic variability. J. Ecol. 110, 1615–1631 (2022).
Lafuente, A. et al. Simulated nitrogen deposition influences soil greenhouse gas fluxes in a Mediterranean dryland. Sci. Total Environ. 737, 139610 (2020).
Li, J. et al. Key microorganisms mediate soil carbon-climate feedbacks in forest ecosystems. Sci. Bull. 66, 2036–2044 (2021).
Fang, C., Smith, P., Moncrieff, J. B. & Smith, J. U. Similar response of labile and resistant soil organic matter pools to changes in temperature. Nature 433, 57–59 (2005).
Karhu, K. et al. Temperature sensitivity of soil respiration rates enhanced by microbial community response. Nature 513, 81–84 (2014).
Li, J., Pei, J., Fang, C., Li, B. & Nie, M. Thermal adaptation of microbial respiration persists throughout long-term soil carbon decomposition. Ecol. Lett. 26, 1803–1814 (2023).
Cardinael, R. et al. Organic carbon decomposition rates with depth and contribution of inorganic carbon to CO2 emissions under a Mediterranean agroforestry system. Eur. J. Soil Sci. 71, 909–923 (2020).
Cerling, T. E. The stable isotopic composition of modern soil carbonate and its relationship to climate. Earth Planet. Sci. Lett. 71, 229–240 (1984).
Cerling, T. E., Solomon, D. K., Quade, J. & Bowman, J. R. On the isotopic composition of carbon in soil carbon dioxide. Geochim. Cosmochim. Acta 55, 3403–3405 (1991).
Meyer, N., Welp, G. & Amelung, W. Effect of sieving and sample storage on soil respiration and its temperature sensitivity (Q10) in mineral soils from Germany. Biol. Fert. Soils 55, 825–832 (2019).
Jia, Y. et al. Plant and microbial pathways driving plant diversity effects on soil carbon accumulation in subtropical forest. Soil Biol. Biochem. 161, 108375 (2021).
Lavallee, J. M., Soong, J. L. & Cotrufo, M. F. Conceptualizing soil organic matter into particulate and mineral‐associated forms to address global change in the 21st century. Glob. Change Biol. 26, 261–273 (2020).
Ye, C. et al. Reconciling multiple impacts of nitrogen enrichment on soil carbon: plant, microbial and geochemical controls. Ecol. Lett. 21, 1162–1173 (2018).
Choi, J. C., Lee, M., Chun, Y., Kim, J. & Oh, S. Chemical composition and source signature of spring aerosol in Seoul, Korea. J. Geophys. Res. Atmos. 106, 18067–18074 (2001).
Ciesielski, H. & Sterckeman, T. Determination of cation exchange capacity and exchangeable cations in soils by means of cobalt hexamine trichloride. Effects of experimental conditions. Agronomie 17, 1–7 (1997).
Sherrod, L. A., Dunn, G., Peterson, G. A. & Kolberg, R. L. Inorganic carbon analysis by modified pressure-calcimeter method. Soil Sci. Soc. Am. J. 66, 299–305 (2002).
Chen, L. et al. Regulation of priming effect by soil organic matter stability over a broad geographic scale. Nat. Commun. 10, 5112 (2019).
Schermelleh-Engel, K., Moosbrugger, H. & Müller, H. Evaluating the fit of structural equation models: tests of significance and descriptive goodness-of-fit measures. Methods Psychol. Res. 8, 23–74 (2003).
Acknowledgements
We thank Xinxin Xu for the assistance in soil sampling and laboratory analyses and Hao Liu for the help with the drawing of the distribution map of soil sampling sites. This work was supported by the Science & Technology Fundamental Resources Investigation Program (2023FY100100), the National Natural Science Foundation of China (92251305, 32101377, and 32101336), the Program of Shanghai Academic/Technology Research Leader (21XD1420700), the Science and Technology Plan Project of Shanghai (23DZ1202700), and the Natural Science Foundation of Shanghai (23ZR1404400).
Author information
Authors and Affiliations
Contributions
J.L. and M.N. designed the research. J.L. conducted the overall experiment with the assistance of J.P.; J.L. performed the overall analysis. J.L. wrote the first draft; J.L., M.N., B.L., and C.F. contributed substantially to reviewing and editing the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks David R. Elliott, Kristen Ball, and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, J., Pei, J., Fang, C. et al. Drought may exacerbate dryland soil inorganic carbon loss under warming climate conditions. Nat Commun 15, 617 (2024). https://doi.org/10.1038/s41467-024-44895-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-44895-y