1 Introduction
Forest ecosystems play a key role in global carbon storage by responsible for a total amount of 861 ± 66 Pg C, of which 44% is stored as soil organic carbon (SOC) (
Pan et al., 2011). Therefore, small changes in forest SOC pool may lead to a substantial fluctuation in atmospheric carbon dioxide (CO
2) concentration, and subsequently affect global climate. Soil microbes use SOC as a source of energy and convert a fraction of carbon (C) into atmosphere as CO
2, a process known as soil organic carbon decomposition or microbial respiration (
Rm) (
Karhu et al., 2014). Global warming may promote
Rm and reduce the ability of forest soil to act as a carbon sink, leading to a positive C cycle-climate feedback (
Bradford et al., 2016). As a key parameter in assessing the response of
Rm to global warming, temperature sensitivity (or
Q10 hereafter) of
Rm represents a factor by which
Rm increases in response to a temperature increase of 10°C (
Walker et al., 2018). Therefore, a better understanding of the controls on
Rm and
Q10 is undoubtedly crucial for improving the representative of C cycle-climate feedbacks in Earth system models under global warming (
Niu et al., 2021). Increasing plant species diversity (PSD) benefits SOC accumulation (
Chen et al., 2023;
Qian et al., 2023), and has been identified as an effective strategy for maintaining or increasing forest SOC pool under global changes (
Ontl et al., 2020). Though increasing PSD has been found to stimulate
Rm, and vice versa (
Chen and Chen, 2019), whether and how PSD is related to
Q10 remain unknown.
Both
Rm and
Q10 are affected by a variety of biotic and abiotic factors, including microbial community abundance and composition, C quantity and quality, SOC protection via mineral association or aggregation (
Qin et al., 2019), but their main controls may be divergent (
Chen et al., 2022). For instance,
Wang et al. (2016) observed that microbial properties strongly affected
Rm while
Q10 was controlled by substrate quality.
Rm is an important process during soil microbial catabolism, so that factors that affect microbial activity would subsequently impact
Rm. Soil temperature and water content directly or indirectly affect
Rm through impacting microbial activities or changing substrate availability (
Ma et al., 2019). Soil pH (
Wong et al., 2008) and substrate availability (
He et al., 2023) also exert significant effects on
Rm by influencing microbial abundance, activity or community composition.
Q10 may be affected by multiple factors including substrate quality, mineral protection of SOC and microbial properties (
Dash et al., 2019). According to the C quality-temperature hypothesis, substrate quality is considered as one of the major drivers in regulating
Q10 (
Liu et al., 2021). The decomposition of low-quality C substrates requires high activation energy, so that is more sensitive to changes in temperature (
Lefèvre et al., 2014). Although many findings are consistent with C quality-temperature hypothesis, cumulative evidence shows that mineral protection would attenuate
Q10 by limiting substrate availability (
Gentsch et al., 2018). On one hand, minerals promote the formation of soil aggregates, which weaken
Q10 through limiting the accessibility of substrates to microorganisms or extracellular enzymes (
Qin et al., 2019). On the other hand, SOC can be absorbed to mineral surface by covalent bonding or electrostatic bonding, so that its decomposition and in turn
Q10 would decline (
Porras et al., 2017). Microbial physiologic characteristics, community structure and composition also affect
Q10 (
Walker et al., 2018). For example, the
Q10 values increased with a higher proportion of K-selected microbial communities (
Li et al., 2021), while microbial diversity was found to negatively impact
Q10 (
Qin et al., 2021).
PSD significantly affects a range of soil biotic and abiotic factors (
Luo et al., 2019) (See Supplementary Materials, Fig. S1). Many studies have shown that SOC and
Rm would be enhanced by increasing PSD, since high PSD promotes detritus inputs via stimulating aboveground and underground plant productivity (
Lange et al., 2015;
Chen et al., 2019a). Meanwhile, PSD directly or indirectly affects soil microbial properties and
Rm by altering the chemical properties of plant detritus inputs (
Cotrufo et al., 2015). Since the detritus composition and quality differ among plant species, they would subsequently change the level and stability of SOC and microbial features (
Jewell et al., 2017) (Fig. S1). More nutrients are needed to satisfy plant nutrient demand owing to stimulated plant productivity under higher PSD, so that plant would invest more resources to enhance root biomass to acquire more nutrients including calcium (Ca) and magnesium (Mg) from deeper soil with the nutrients eventually returning to top soil layer through litter inputs (
Furey and Tilman, 2021). Meanwhile, greater root biomass would be accompanied by more organic acid release, so that increase soil calcium and other minerals via promoting bedrock weathering (
Beugnon et al., 2021). Minerals such as iron (Fe), aluminum (Al), Ca or Mg can absorb SOC and form mineral-organic associations, which may contribute to SOC stabilization and in turn restrain
Rm and
Q10 (
McGrath et al., 2022) (Fig. S1). The alteration in soil minerals by PSD may affect soil moisture and pH that are crucial factors affecting
Rm and
Q10 (
Wu et al., 2020) (Fig. S1). However, though quite a few studies have determined the effects of PSD on
Rm (
Chen and Chen, 2019), relatively few studies have explored the effect of PSD on
Q10 (
Luan et al., 2018).
Luan et al. (2018) reported that increasing PSD suppressed
Q10 in temperate forests, but the underlying mechanisms were not explored. Considering that increasing PSD is an effective way to enhance SOC accumulation during vegetation restoration or afforestation, knowledge about how increasing PSD alters the temperature sensitivity of SOC decomposition is urgently needed under the background of climate change.
Natural PSD gradients in mature forests have often been used to explore the impacts of increasing PSD on ecosystem function and process (
Qian et al., 2023;
Zhu et al., 2023). In the current study, 43 plots covering a natural gradient of PSD were selected in a subtropical karst forest in south-west China (
Qian et al., 2023;
Zhu et al., 2023) to explore the effects of increasing PSD on
Rm and
Q10. The major objectives were to address: 1) how
Rm and
Q10 respond to increasing PSD; 2) what mechanisms underline the effects of PSD on
Rm and
Q10.
2 Materials and methods
2.1 Experimental design and sample collection
This study was carried out in Mulun National Nature Reserve (25°06′09″−25°12′25″ N, 107°53′29″−108°05′45″ E), southwest China. The region is located in the subtropical monsoon climate zone with a mean annual temperature of 19°C and mean annual precipitation of 1389 mm (
Zhu et al., 2023). The selected forest is a subtropical karst evergreen and deciduous broad-leafed forest with the underlain bedrock being a mixture of limestone and dolomite. The soil is calcareous with soil types being Cambisols, Luvisols, or Leptosols depending on topographic positions (
Zhu et al., 2023). Forty-five plots with an area of 20 m × 20 m each were established in the valley and lower slope with the elevation ranging from 446 m to 521 m (See Supplementary Materials, Fig. S2) and soil type being Luvisols (
Zhu et al., 2023). The interval between any two adjacent plots is 40 m. All the plants with a diameter at the breast height greater than 1 cm in each plot were measured, labeled and recorded. Shannon-Weiner index was calculated to represent PSD, which ranged from 0.15 to 3.53. Detailed information of the dominant plant species in each plot is presented in Supplementary Data 1.
Field sampling was carried out between July 28 and August 1, 2020. In each plot, 16 mineral soil samples were randomly collected with a soil corer (5 cm in inner diameter) to a depth of 10 cm after removing the surface organic matter. The soil samples in each plot were thoroughly mixed into one composite sample and sieved to 2 mm with stones and plant residues removed. Each composite sample was divided into three parts. One part was stored at 4°C for the determination of Rm, Q10, dissolved organic C (DOC). The second part was stored at −40°C for measuring microbial properties. The third part was air-dried to analyze physicochemical variables, including soil pH, texture, SOC, exchangeable Ca, exchangeable Mg exchangeable Fe, exchangeable Al. It should be noted that only 43 plots were included due to soil missing in the current study.
2.2 Determination of soil microbial respiration and its temperature sensitivity
Fresh soil samples of 12 g each (dry weight equivalent) were weighted into 125 mL glass jars with two replicates per sample. The incubation was conducted in a biochemical incubator (LRH-150F, Shanghai Daohan Industry Co., Ltd, Shanghai, China) with the temperature being adjusted manually. Each jar was covered with Parafilm
® to prevent evaporation after the soil moisture was adjusted to 60% water holding capacity (WHC), and the soils were pre-incubated at 25°C under dark conditions for 7 days. A scheme of changing temperatures was adopted to determine soil microbial respiration and its temperature sensitivity (
Liu et al., 2017;
Liu et al., 2019;
Liu et al., 2021). The incubation temperature was changed from 15°C to 35°C and then down to 15°C with a changing step of 5°C, i.e., 15°C – 20°C – 25°C – 30°C – 35°C – 30°C – 25°C – 20°C – 15°C. When reaching a new temperature, an equilibration period of 3 h was applied to allow soil microbes to adapt to the new temperature (
Ding et al., 2016). Soil moisture was regularly checked and adjusted if necessary, during the incubation period. During the equilibration period of 3 h, each jar was covered with Parafilm
® to let air exchange between the headspace and the ambient atmosphere but reduce evaporation. Upon the end of equilibration period, each jar was sealed with an air-tight lid with a rubber septum for the collection of headspace gas. Twenty ml of laboratory air was immediately injected to the headspace of each jar and mixed well using a 60-mL syringe, followed by 20 mL headspace air being collected and stored in a 12-mL vacuumed vial (839W, Labco, England). A second headspace gas sample was collected from each jar after 18 h, 12 h, 8 h, 3 h, and 1 h following the collection of the first sample for the incubation temperature of 15°C, 20°C, 25°C, 30°C, and 35°C, respectively. The air-tight lid was replaced by the Parafilm
® after the second headspace gas collection, followed by a new equilibration period under another incubation temperature. The experiment only lasted for one round of incubation with changing temperature (~5 days). The sealed incubation period was determined based on a preliminary experiment to avoid the CO
2 concentration in the jar from exceeding 0.3% (
Jewell et al., 2017). The lips were replaced by Parafilm
® again after each Rm determination. CO
2 concentrations were analyzed using an Agilent 7890A gas chromatograph equipped with a thermal conductivity detector (Agilent Technologies, Santa Clara, CA, USA). The soil microbial respiration (
Rm, μg C·g
−1 soil·h
−1) was calculated using Eq. (1):
where
M is the molar mass of C (g·mol
−1); 22.4 is the molar mass of gas in the standard conditions (273 K, 1013 hPa);
T0 and
P0 are the air temperature (K) and pressure (hPa) under standard conditions, respectively;
T and
P are the determined air temperature (K) and pressure (hPa), respectively; ΔC and Δ
t are the changes of CO
2 concentrations (ppm) and time (h), respectively, between the two sampling time points for each
Rm determination;
V is the headspace volume of jar (L); and
m is the mass of dry soil (g). Since soil microbial respiration has been found to increase up to a threshold of ~25°C across all non-desert biomes (
Carey et al., 2016), microbial respiration under 25°C was used to assess the impact of increasing PSD on
Rm. Similarly, the respiration under 20°C or 25°C from incubation with changing temperatures has been used to assess
Rm variation or response to environmental change (
Wang et al., 2016).
The exponential function was used to fit the correlation of Rm with temperature (R2 > 0.95), and Q10 was calculated using Eqs. (2) and (3):
where T is temperature (°C); and B and k are model parameters.
2.3 Measurements of auxiliary variables
Litter, root and soil physicochemical properties were determined using methods described in companion studies (
Duan et al., 2023;
Zhu et al., 2023) (See Supplementary Materials, Table S1). To investigate the relationship between
Rm or
Q10 and mineral protection, we calculated the clay to SOC content ratio (Clay/SOC) to evaluate how clay particles contribute to SOC protection (
Fang et al., 2021;
Duan et al., 2023).
The functional group composition of SOC was determined using a Fourier transform infrared spectrometer (Nicolet is50, Thermo Fisher, USA). Samples were first dried under 60°C until completely dry, then ground to a fine power. Spectra were recorded by averaging 40 scans at 4 cm
−1 resolution over the range of 4000 to 400 cm
−1. Specifically, polysaccharides (Poly), phenolic and aliphatic (Ali), and aromatics and aromatic or aliphatic carboxylates (Aro) were represented by the absorption bands at 950–1170 cm
−1 (1033 cm
−1), 1420 cm
−1 and 1630 cm
−1, respectively. The relative peak ratios of 1630/1033 (Aro/Poly) and 1420/1033 (Ali/Poly) were used as recalcitrance indices of SOC (
Wang et al., 2016;
Qin et al., 2021), with greater ratios implying higher recalcitrance.
Microbial abundance and community composition were assessed by phospholipid fatty acids (PLFAs). Details of PLFAs assay are presented in a companion study (
Li et al., 2023). The F/B ratio represent fungal-to-bacterial biomass ratio.
2.4 Statistical analysis
The normality of the data was checked using Shapiro-Wilk test, and logarithmical transformation was applied if necessary. Correlation analyses were conducted to examine the relationships of Rm (the respiration at 25°C) or Q10 with the measured biotic and abiotic variables. Random forest analysis was then conducted to identify the major variables affecting Rm and Q10. To clarify the direct and indirect drivers of Rm or Q10, the partial least squares path model (PLS-PM) was utilized. Prior to PLS-PM, all selected variables were standardized. To reduce collinearity among the predictor variables, principal component analysis (PCA) was conducted for substrate properties, mineral protection and microbial traits with the first component (PC1) was used for PLS-PM (See Supplementary Materials, Table S2). All the possible models based on the priori model (Fig. S1) were tested with the non-significant pathways removed. The statistical analyses were implemented in R software 4.2.2.
3 Results
Soil microbial respiration ranged from 0.12−0.59 μg C·g−1 soil·h−1, with an average of 0.30 ± 0.11 μg C·g−1 soil·h−1 (Fig. 1(a)). Q10 ranged from 3.30 to 7.16, with a mean value of 4.87 ± 0.97 (Fig. 1(b)). Rm was significantly and positively correlated with PSD (P < 0.01) (Fig. 1(c)), but no significant correlation was found between Q10 and PSD (Fig. 1(d)).
Pearson correlation analysis demonstrated that Rm was significantly and positively correlated with plant properties (litter C/P, litter N/P and fine root biomass), soil water content (SWC), substrate properties (SOC, DOC, Aro/Poly) and microbial biomass (total PLFAs, bacterial PLFAs and fungal PLFAs). In contrast, Rm exhibited negative correlations with Clay/SOC and F/B (Fig. 2(a)). Q10 was positively correlated with pH, Ali/Poly, (Ca + Mg)exe and microbial biomass (total PLFAs, bacterial PLFAs, and fungal PLFAs), but negatively correlated with (Fe + Al)o (Fig. 2(b)). Random forest analysis revealed that Clay/SOC, F/B, SWC, SOC and total PLFAs were identified as the strong explanatory factors for Rm, with the total explanation being 48.59% (Fig. 3(a)). For Q10, pH, F/B, fine root biomass, SOC were selected as the strong explanatory factors, with the total explanation being 25.02% (Fig. 3(b)).
According to the partial least squares path model, increasing PSD stimulated Rm via promoting SWC and total PLFAs but suppressing F/B, owing to increased substrate level caused by improved fine root biomass but lowered Clay/SOC (Fig. 4(a)). The total effects of PSD, fine root biomass, Clay/SOC, substrate properties, SWC and microbial traits were 0.43, 0.15, −0.48, 0.64, 0.46, and 0.38, respectively (Fig. 4(b)). Increasing PSD had neutral effect on Q10 because of the contrasting effects of soil biotic and abiotic variables. On one hand, increasing PSD suppressed Q10 via decreasing F/B owing to increased subtrate level caused by promoted fine root biomass and mineral protection (Fig. 4(c)). On the other hand, increasing PSD stimulated Q10 by enhancing soil pH through elevating (Ca + Mg)exe (Fig. 4(c)). The total effects of PSD, fine root biomass, mineral protection, substrate properties, pH and F/B were 0.15, −0.04, 0.28, −0.10, 0.58 and −0.23, respectively (Fig. 4(d)). The partial least squares path models explained 58% and 37% of the variance in Rm and Q10, respectively (Figs. 4(b) and 4(d)).
4 Discussion
4.1 Increasing plant species diversity stimulates soil microbial respiration
The current study shows that increasing PSD significantly stimulated
Rm. Consistently,
Rm was found to increase as PSD increased (
Chen et al., 2019a), but decrease when PSD lost (
Chen and Chen, 2019).
Rm is usually influenced by soil microbiome, substrate quantity and quality, mineral protection and environmental factors (
Sáez-Sandino et al., 2023). These are applicable to the present study. The stimulated input of plant residues or root exudates under higher PSD leads to greater levels of substrate quantity (
Chen et al., 2019b;
Fang et al., 2024). Though plant residue inputs were not determined in the current study, they were most likely stimulated as suggested by the increased fine root biomass under higher PSD. Higher substrate availability would benefit microbial growth and activity, and in turn stimulate microbial biomass and respiration (
Prommer et al., 2020). This is corroborated by the current study. Additionally, the soils under higher PSD would have greater water holding capacity, which is positively correlated with soil organic matter (
Li et al., 2007), and hence would have higher SWC. Relatively high SWC is necessary for maintaining microbial activity and
Rm (
Schimel, 2018). Consistently, both SWC and
Rm were stimulated as PSD increased in the current study.
Soil microbiomes represent the ultimate actors in the decomposition of SOC with higher microbial biomass usually accompanied by great
Rm under similar conditions (
Sáez-Sandino et al., 2023). Additionally, carbon use strategies differ between fungi and bacteria, with fungi having a lower respiration quotient and a higher carbon use efficiency than bacteria (
Soares and Rousk, 2019;
Zhang et al., 2020). For example, Zhang et al. (
2020) observed a negative correlation between F/B ratio and
Rm in forest ecosystems. Some studies show that F/B ratio increases with PSD (
Beugnon et al., 2021;
Chen et al., 2019b). In contrast, F/B ratio decreased as PSD increased in the current study, likely owing to the increase of labile organic C (e.g., DOC), which favors the growth of bacteria (
Huang et al., 2021). Therefore, the increase of microbial biomass and decrease of F/B can partly explain the stimulated
Rm under higher PSD in the current study.
4.2 Increasing plant species diversity does not alter temperature sensitivity of microbial respiration
The current study demonstrates that increasing PSD had no significant effect on
Q10 due to the contrasting effects of soil biotic and abiotic variables. On one hand, increasing PSD stimulated
Q10 by enhancing soil pH through increasing exchangeable Ca and Mg. Two possible mechanisms exist for the increase of exchangeable Ca and Mg under higher PSD. First, increasing PSD stimulates the root investment to obtain nutrients including Ca and Mg from deep soil, and the nutrients would be returned to the top soil layer via detritus inputs (
Furey and Tilman, 2021). In the current study, fine root biomass was significantly stimulated under higher PSD as reported by a companion study (
Zhu et al., 2023). Second, the increase of fine root biomass usually stimulates release of low-molecular-weight organic acids, which subsequently participate in the weathering of primary minerals or bedrocks (a mixture of limestone and dolomite in the current study), so that increase the release of dissolve Ca and Mg (
Qin et al., 2019;
Beugnon et al., 2021). High concentrations of exchangeable Ca may contribute to the stabilization of SOC through inner or outer sphere interactions and benefit SOC accumulation (
Rowley et al., 2018). As exchangeable Ca and Mg accumulate, the poorly crystalline Fe and Al minerals would be transformed to more crystalline phases and their reactivity accordingly decreases (
Goldberg, 2014). As basic cations, the increase of soil exchangeable Ca and Mg would increase soil pH (
Leiva Soto et al., 2023). The increased pH subsequently stimulated
Q10 in the current study, consistent with some previous researches (
Ding et al., 2016;
Liu et al., 2017). Although the direct linkage between soil pH and
Q10 has not be well explored, it is well known that microbial activity and community structure are regulated by soil pH (
Rousk et al., 2010;
Meyer et al., 2018). For example, some studies show that soil pH is the dominant and direct control on
Q10 probably owing to its effect on microbial activity and community structure (
Min et al., 2014;
Meyer et al., 2018). Nevertheless, both positive (
Min et al., 2014;
Ding et al., 2016;
Liu et al., 2017;
Komarova et al., 2022) and negative (
Min et al., 2014;
Meyer et al., 2018) effects of soil pH on
Q10 have been observed likely caused by the differential effects of soil pH on microbial activity and structure under different soil conditions. The mechanism underlying the positive effect of soil pH on
Q10 deserve further investigation.
On the other hand, we find that PSD suppressed
Q10 via decreasing F/B ratio and increasing SOC quantity and recalcitrance. According to the C quality-temperature hypothesis,
Q10 should increase with substrate recalcitrance (
Jewell et al., 2017;
Wang et al., 2018). With regard to the current study, SOC recalcitrance may have increased under higher PSD according to the increased ratio of aliphatics/polysaccharides (
Wang et al., 2016). Consistent with the C quality-temperature hypothesis, significantly positive correlation was observed between the ratio of aliphatics/polysaccharides and
Q10 (Fig. 2(b)). Similarly,
Wang et al. (2018) found that
Q10 increased with indices of SOC recalcitrance in coniferous forests. Although the recalcitrance increased with PSD, labile C fraction as indexed by DOC also increased with PSD in the current study (See Supplementary Materials, Fig. S3). The increase in labile C availability with PSD could shift microbial community composition toward bacterial dominance, as bacteria preferentially utilize labile substrates, thereby reducing the F/B ratio (
Huang et al., 2021). Given that fungi are the primary decomposers of recalcitrant C (
Fontaine et al., 2011;
Ma et al., 2024), a lower F/B ratio may reduce the microbial capacity to utilize recalcitrant C pools. This shift in decomposition pathways could subsequently lower
Q10 of soil organic matter decomposition, as recalcitrant C typically exhibits higher
Q10 values due to its greater activation energy requirements for breakdown (
Li et al., 2021;
Pei et al., 2024). In summary, our study explored the effects of PSD on
Rm and its temperature sensitivity in a subtropical karst forest with calcareous soil. Increasing PSD stimulated
Rm by increasing substrate quantity, microbial abundance and SWC but decreasing F/B. Nevertheless, increasing PSD had no significant effect on
Q10 due to the contrasting effects from soil biotic and abiotic variables, i.e., stimulating
Q10 by increasing soil pH, but suppressing
Q10 via decreasing F/B. These findings suggest that increasing PSD would stimulate microbial respiration along with soil organic carbon accumulation, but not alter the carbon cycle-climate feedback when implementing ecological restoration or afforestation. Additionally, our study, for the first time, provides the mechanism underlying the response of temperature sensitivity of microbial respiration to PSD, which should be considered when assessing the impacts of PSD on soil organic carbon dynamics under climate change.