Distinct temperature sensitivity of soil carbon decomposition in forest organic layer and mineral soil.

The roles of substrate availability and quality in determining temperature sensitivity (Q10) of soil carbon (C) decomposition are still unclear, which limits our ability to predict how soil C storage and cycling would respond to climate change. Here we determined Q10 in surface organic layer and subsurface mineral soil along an elevation gradient in a temperate forest ecosystem. Q10 was calculated by comparing the times required to respire a given amount of soil C at 15 and 25°C in a 350-day incubation. Results indicated that Q10 of the organic layer was 0.22-0.71 (absolute difference) higher than Q10 of the mineral soil. Q10 in both the organic layer (2.5-3.4) and the mineral soil (2.1-2.8) increased with decreasing substrate quality during the incubation. This enhancement of Q10 over incubation time in both layers suggested that Q10 of more labile C was lower than that of more recalcitrant C, consistent with the Arrhenius kinetics. No clear trend of Q10 was found along the elevation gradient. Because the soil organic C pool of the organic layer in temperate forests is large, its higher temperature sensitivity highlights its importance in C cycling under global warming.

Soil is the largest reservoir of terrestrial carbon (C) compared to vegetation and atmosphere1 and has diverse soil organic C (SOC) compounds with different stability2. Concerns about the clin class="Disease">mate–C cycle feedback under current global climate change have catalyzed efforts on SOC decomposition and its temperature sensitivity (Q10)23456, which represents change in C decomposition rate to a 10°C temperature increase. However, as yet there is no consensus on the relationship of SOC decomposition with temperature owing to various stability of different SOC components257. As a result, doubts have been cast on the assumed similar sensitivities of discrete C pools to temperature in soil C models8.

Most uncertainty in predicting responses of C storage to changing climate comes from temperature sensitivity of recalcitrant soil C. Compared to more labile C, recalcitrant C has been found to have similar49 or lower310 Q10 of C decomposition. However, substantial researches in recent years supported the Arrhenius kinetic theory27 and indicated that recalcitrant C with higher activation energy has greater Q10 than labile C111213141516. Because soil C quality generally decreases with increasing depth through mineral soil profile due to possible differences in physical and/or chemical protection from C decomposition and microbial community composition17, the increase in Q10 with increasing C recalcitrance, also known as the “C quality-dependent hypothesis” of Q10, has thus been tested at different depths of mineral soils1416. In boreal and temperate forests with organic forest floor, a large amount of SOC is stored in the organic layer besides the mineral soil. The organic layer and the mineral soil have different rates of C input, accumulation, and turnover, resulting in different substrate availability and quality for C decomposition1819. Substrate availability, besides substrate quality, could affect Q10 based on the Michaelis–Menten kinetics22021. Consequently, the temperature sensitivity of C decomposition are likely to differ between the two soil layers1322.

It has been suggested that Q10 of C decomposition may be latitude-dependent, with greater Q10 at higher latitude2324. The altitudinal pattern of C decomposition with temperature change may be comparable with that caused by latitudinal gradient1422. Therefore, elevation gradients in mountains are well recognized as “natural experiments” to explore responses of SOC to temperature changes25. Some studies suggested Q10 values of both labile and recalcitrant C increased along the elevation gradient in Mountain Wuyi of China1416. On the contrary, other studies found no altitudinal trend for Q10 of C decomposition2226. These inconsistencies may be explained by further understanding of the interactions of intrinsic (e.g. substrate availability) and extrinsic controlling factors (e.g. soil texture) of temperature sensitivity of C decomposition.

Here we collected soils from the northern slope of Mountain Changbai in n class="Chemical">Northeast China27 (Table 1) to investigate variations of Q10 of decomposition along the elevation gradient. Soil samples from both the organic and the mineral layers were incubated for 350 days and Q10 of decomposition was determined by the approach derived from Conant et al.11 (see Methods). The specific questions we aim to address in this study include: (1) Does Q10 of C decomposition of both the mineral layer and the organic layer support the “C quality-dependent hypothesis”? (2) Whether and how does Q10 of the organic layer differ from that of the mineral soil along the elevation gradient?

Table 1

Sites description

Elevation (m)	Position	MAT (°C)	MAP (mm)	Tree category	Vegetation type	Dominant species	Soil type	The depth of organic layer§
795	42°23′N, 128°05′E	2.0	702	Mixed coniferous broad-leaved forest	Deciduous broadleaf and evengreen needleleaf	Pinus koraiensis Tilia amuresis Fraxinus mandshurica	Dark brown forest soil	4.9 ± 0.2a
1102	42°11′N, 128°09′E	0.4	783	Broad-leaved Korean pine forest	Deciduous broadleaf and evengreen needleleaf	Pinus koraiensis Picea jezoensis	Brown coniferous forest soil	3.7 ± 0.2b
1248	42°08′N, 128°07′E	−0.3	824	Korean pine-spruce-fir forest	Evengreen needleleaf	Abies nephrolepis Picea jezoensis Larix gmellini	Brown coniferous forest soil	3.6 ± 0.3b
1510	42°05′N, 128°04′E	−1.7	904	Dark coniferous spruce-fir forest	Evengreen needleleaf	Picea jezoensis Abies nephrolepis	Brown coniferous forest soil	3.3 ± 0.1b
1690	42°04′N, 128°03′E	−2.6	964	Ermans birch-spruce-fir forest	Evengreen needleleaf	Picea jezoensis Abies nephrolepis Betula ermanii	Brown coniferous forest soil	3.6 ± 0.6b
1900	42°03′N, 128°04′E	−3.7	1038	Ermans birch forest	Deciduous broadleaf	Betula ermanii	Soddy forest soil	3.9 ± 0.4b

§Data are mean ± SE (n = 4). Lowercase letters next to data indicate significant differences (p < 0.05) in the depths of the organic layer among the 6 elevations determined by t test.

Results

Soil characteristics

Soil organic C (SOC) and total N (Tn class="Chemical">N) in the organic layer were 1.9–6.8 and 1.8–6.1 times higher (both p < 0.05) than those in the mineral soil across the 6 elevation sites (Fig. 1a and b, Table 2). Correspondingly, C:N ratio (p < 0.05) in the organic layer was higher compared to the mineral soil (Fig. 1c). In the organic layers, SOC ranged from 90.3 g kg−1 at 1900 m to 165.2 g kg−1 at 1510 m (Fig. 1a). However, TN in the organic layer was not influenced by elevation whereas its values varied from 6.5 g kg−1 at 1248 m to 8.9 g kg−1 at 795 m (Fig. 1b). SOC and TN in the mineral soil increased from 17.9 to 46.6 g kg−1 and from 1.1 to 3.7 g kg−1 along the elevation gradients, respectively. Soil pH and texture were also different among the 6 sites (all p < 0.05) and between the two soil layers (all p < 0.05 except for silt content) (Fig. 1d–g, Table 2). Soil organic layer had higher clay content and lower sand content than the mineral soil across the 6 sites (Fig. 1e and g).

Figure 1

Soil properties and Q10 in the organic layer and the mineral soil along the elevation gradient (mean ± SE, n = 4).

*: p < 0.05, which indicates significant differences between the two soil layers.

Table 2

Results (F-values) of two-way ANOVA for soil properties and Q10. *: p < 0.05, ns: p > 0.05

					Soil texture
Source of variance	SOC	TN	C:N ratio	pH	Sand	Silt	Clay	Q10–1st	Q10–2nd	Q10–3rd	Q10–4th	Q10–5th	Q10–6th
Soil layer (SL)	201.5*	247.4*	18.6*	10.6*	34.5*	0.3ns	244.2*	13.5*	2.6ns	15.8*	20.3*	17.2*	8.2*
Elevation (E)	3.6*	2.9*	71.7*	16.4*	8.3*	8.6*	7.1*	5.8*	1.3ns	2.4ns	3.1*	5.1*	4.8*
SL × E	3.6*	3.1*	1.4ns	9.8*	4.8*	3.1*	21.4*	1.2ns	2.6*	4.4*	9.1*	6.5*	5.1*

SOC decomposition

Soil respiration rates decreased by 343.8% for the organic layer and by 575.8% for the mineral layer after the 350-day incubation period at 15°C (both p < 0.05, Fig. 2a, Table 3). Both respiration rates and cumulative proportion of soil C respired incubated at 25°C were higher than those at 15°C during the whole incubation period (both p < 0.05, Fig. 2). Means of soil respiration rates over the whole incubation period was 7.6 times greater (p < 0.05) in the organic layer than in the mineral soil across the incubation temperatures and the elevation gradients (Fig. 2a). Consequently, the cumulative proportion of C respired after 350 days of incubation in the organic layer (22.1% at 25°C and 11.3% at 15°C, respectively) were 1.6 times higher (both p < 0.05) than those in the mineral soil (13.4% at 25°C and 7.0% at 15°C, respectively) (Fig. 2b). In addition, both of the two measured variables were different among sites with different elevations (both p < 0.05, Table 3). Soil C decomposition rates incubated at 15°C in the organic layer and the mineral soil, for example, ranged from 30.0 µg CO2-C g soil−1 day−1 at 1690 m to 86.9 µg n class="Chemical">CO2-C g soil−1 day−1 at 1510 m and from 3.2 µg CO2-C g soil−1 day−1 at 795 m to 12.3 µg CO2-C g soil−1 day−1 at 1900 m, respectively.

Figure 2

Soil respiration rate (a) and cumulative proportion of soil C respired (b) in the organic layer and the mineral soil across the 6 elevations at different incubation temperature during the incubation (mean ± SE, n = 6).

Table 3

Results (F-values) of repeated-measures ANOVA for soil respiration rate and cumulative proportion of soil C respired among the 22 sampling times during the incubation.*: p < 0.05, ns: p > 0.05

Source of variance	Sampling time	Incubation temperature (T)	Soil layer (SL)	Elevation (E)	T × SL	T × E	SL × E	T × SL × E
Soil respiration rate	271.4*	159.2*	848.4*	28.7*	100.1*	2.5*	20.4*	1.7ns
Cumulative proportion of soil C respired	1899.9*	219.8*	141.8*	11.9*	23.9*	1.5ns	4.2**	0.6ns

Temperature sensitivity

The temperature sensitivity (Q10) of SOC decomposition increased with increasing proportion of C respired over the incubation period in both the organic layer and the mineral soil across the elevation gradient (Fig. 3). The Q10 values in the two soil layers were estimated from 1 to 6% of the proportion of respired C to initial SOC by a step of 1% and named as Q10–1st, Q10–2nd, Q10–3rd, Q10–4th, Q10–5th, or Q10–6th, respectively. Q10–3rd, for example, was determined by dividing the time taken to respire the third 1% of initial C after 2% of initial C was respired at 15°C by that at 25°C. When more than 5% and 4% of the initial SOC had been decomposed in the organic layer and the mineral soil, respectively, Q10 of the 1% SOC (Q10–5th for the organic layer and Q10–4th for the mineral layer) began to be different from Q10–1st significantly (Fig. 3). Q10 values of the organic layer (2.5–3.4) were always higher (all p < 0.05 except for Q10–2nd) compared to those of the mineral soil (2.1–2.8) during the incubation (Fig. 3, Table 2). There was no altitudinal trend of Q10 (Q10-1st–Q10-6th) in either soil layer. For example, the highest Q10–1st value of the organic layer and the mineral soil occurred at 1248 m (3.3) and 1102 m (2.4), respectively (Fig. 1h). In addition, elevation also interacted with soil layer to influence Q10 (Table 2).

Figure 3

Changes in Q10 associated with proportion of soil C respired in the organic layer and the mineral soil along the elevations during the incubation (mean ± SE, n = 6).

Bars with different upper- and lower-case letters are significantly different (p < 0.05) in the organic layer and the mineral soil, respectively.

Discussion

The alteration of soil C quality can arise through the progressive depletion of labile C and the increasing contribution of recalcitrant C to SOC decomposition during the long-term incubation1114 compared to the relatively short incubation duration4. In this 350-day incubation study, the curves of changes in SOC decomposition rates (Fig. 2a) revealed that respiration declined sharply in the first 3 weeks and then started to settle and reduce slowly in both the organic layer and the mineral soil. This observed pattern was consistent with previous long-term laboratory incubation studies11141516. The enhancement of Q10 over incubation time in both the organic layer (from 2.5 in Q10–1st to 3.4 in Q10–6th) and the mineral soil (from 2.1 in Q10–1st to 2.8 in Q10–6th) (Fig. 3) suggested that the temperature sensitivity of more labile C was lower than that of more recalcitrant C. This result could be well explained by the Arrhenius equation (k = A exp(–E/RT)), where k is the reaction rate of C decomposition, A is the frequency or pre-exponential factor (the theoretical reaction rate at E = 0), E is the activation energy of C components (J mol−1), R is the gas constant (8.314 J mol−1 K−1), and T is the absolute temperature (in Kelvin). It indicated that higher E was required for decomposition of more recalcitrant C due to less reaction of stabilized substrates27. Q10, defined as the factor by which k increases with increasing temperature by 10°C, should increase with increasing E and C recalcitrance. Therefore, our results support the “C quality-dependent hypothesis”, which has been proven in temperate grasn class="Disease">sland1115, subtropical forest1416, boreal forest13, and cropland12. Our findings in the temperate forest, together with these previous reports point to the critical role of the relatively recalcitrant C fraction with higher amounts and higher temperature sensitivity in affecting C balance and cycling under further climate warming.

We found constantly higher Q10 values in the surface organic layer than the subsurface mineral soil during the incubation (Fig. 3), which was distinct from previous findings of increasing Q10 with soil depths in the mineral soil profiles141617. In addition to Arrhenius kinetics, Q10 of decomposition may be affected by Michaelis–Menten kinetics simultaneously. We believe this finding may be explained by the Michaelis–Menten kinetics (k = Vmax × [S]/(Km + [S]), where k is the reaction rate, Vmax is the maximum rate of enzyn class="Disease">matic activity at a given temperature, [S] is the substrate availability (substrate concentration at active site of enzyme), and Km is the Michaelis–Menten constant, representing the affinity of enzymes for the substrates expressed as the substrate concentration at which the reaction rate equals Vmax/2. Both Vmax and Km are temperature dependent and their temperature sensitivities can neutralize each other, called “cancelling effect”. When the substrate is least limiting ([S] ≫ Km), this “cancelling effect” is eliminated, causing higher Q10. Therefore, the higher Q10 values in the organic layer than in the mineral soil could be due to the reduction in the “cancelling effect” between Vmax and Km by the significantly higher C availability in the organic layer1320.

The Q10 values observed in both the organic layer and the mineral soil during the incubation (Fig. 3) were larger than 2.0, which has been typically used as Q10 value of different SOC pools in current soil C models3. Therefore, care should be taken when applying this generally used Q10 value to different ecosystems. More importantly, various Q10 values in different soil C pools should be considered to be incorporated into soil C models for better simulation of C decomposition in these forests. In particular, soil organic layer and mineral layer should have different Q10 values when modeling temperate forest with thick organic layers. In addition, Q10 values only reflect the relative changes in decomposition rates as a function of temperature. The absolute amounts of C that can be released from soils are the most important values to be predicted under climate change. Therefore, both the absolute and relative rates in response to temperature change should be considered in combination228. In this study, the high amount of SOC and Q10 of C decomposition in the organic layer (Fig. 1a and 3) could have additive effects on C release, suggesting the extreme importance of the organic layer.

We did not find a clear trend of Q10 along the elevation gradient, which is in accordance with findings in temperate and boreal forests2226, whereas increase of Q10 with increasing elevation was observed in subtropical forests1416. These contradictory results could be explained by the different controls of Q10 under different conditions. Soil properties have been postulated to play roles in regulating Q10 due to their influences on decomposability of soil C15. For example, in mineral soil of subtropical forests, Xu et al.14 and Wang et al.16 found that both microbial biomass C and soil C:N ratio increased with increasing elevation, which could result in the positive relationship of Q10 with elevation. While in our study, SOC, Tn class="Chemical">N, and soil C:N ratio did not show any altitudinal patterns in either the organic layer or the mineral soil (Fig. 1a and b), which may be the reason why no clear trend of Q10 was observed along the elevation gradient. For the organic soils in Niklińska and Klimek26, we re-calculated soil C:N ratios based on the soil C and N concentrations in their paper and found no effect of elevation on C:N ratio. Therefore, how elevation may affect of Q10 decomposition actually depends on how elevation influences quality and quantity of soil organic matter and/or activity and abundance of decomposers.

In conclusion, our results suggested the organic layer had higher Q10 than the mineral soil, which is consistent with the Michaelis-Menten kinetics. This highlights the necessity of separating the organic layer from the mineral soil to explore the responses of soil C decomposition to temperature. When the two layers were considered separately, the Arrhenius kinetic theory was supported by our results which found decomposition Q10 increased with decreasing soil C lability. Because the SOC pool of the organic layer in temperate forests is large, its higher temperature sensitivity highlights its importance in C cycling under global warming. The different Q10 values observed in the two soil layers should be incorporated into soil C models for better prediction of responses of C decomposition to temperature change in temperate ecosystems.

Methods

Site description and soil sampling

The study was conducted on the northern slope of Mountain Changbai in n class="Chemical">Northeast China (41°42′–42°25′N, 127°42′–128°17′E). Mean annual temperature declined from 2.0 to −3.7°C and mean annual precipitation increased from 702 to 1038 mm with increasing altitude from 795 to 1900 m. Forest types were especially distinct due to climatic changes along the elevation gradient. A 20 × 20 m2 plot was set at each of the six major vegetation types with different elevation. The detailed site descriptions were shown in Table 1.

Soil samples of the organic layer and the upper 10 cm of the mineral soil were collected at 4 points randomly within each plot in September 2011. The organic layer was distinguished from the mineral soil by its morphology (including soil color, texture, and consistency). Soils were black or dark brown in color and less friable in the organic layer compared to the mineral soil. In addition, the depths of the organic layer were 3.3–4.9 cm across the 6 elevations (Table 1). Each fresh soil sample was sieved with a 2-mm sieve to remove rocks and root materials, thoroughly hand-mixed, and divided into two subsamples. One subsample was stored at 4°C until incubation began. The other subsample was air-dried to measure soil properties. SOC, total n class="Chemical">nitrogen (TN), and C:N ratio were determined on oven-dried (60°C) soil subsample by an elemental analyzer (VARIO EL III, Germany). Soil pH was measured with a Calomel electrode at 1:5 soil-to-water ratio. Soil texture was measured by the pipette method29. Soil water content in both soil layers and field water capacity of mineral soil samples were also determined.

Soil incubation

Fresh soil samples of the organic layer (50 g with original soil moisture) and the mineral layer (150 g with 60% of field capacity) were incubated in 1 L Mason jars at 15 and 25°C, respectively. Empty jars without soil samples were used as controls. The soil samples were pre-incubated for five days before measurements of SOC decomposition by the alkali absorption method14. In detail, soil respired CO2 was trapped by 30 ml n class="Chemical">NaOH (1M) solution in an open vial placed in each sealed jar. Samples were taken 22 times (2, 5, 9, 14, 21, 28, 35, 49, 63, 77, 91, 105, 126, 147, 168, 189, 210, 238, 266, 294, 322, 350 days after the incubation) with different intervals. The amount of respired CO2 was determined by titration of NaOH with 1 M HCl to pH = 8.3 in the presence of BaCl2. Samples were flushed with compressed air to allow replenishment of O2 and were remoistened to maintain moisture after each interval. The cumulative proportion of soil C respired at a given sampling time was calculated as the summed amounts of C respired before the given sampling time divided by SOC.

Temperature sensitivity (Q10) of SOC decomposition during the incubation was calculated using the method developed by Conant et al.11: where t15 and t25 are the time required to respire a given amount of soil C during the incubation at 15 and 25°C, respectively.

Statistical analysis

Shapiro-Wilk and Levene's tests were used to analyze the data for normality and homogeneity of variance, and if necessary, data were transformed by the square root or logarithmic transformation. Two-way An class="Chemical">NOVAs were used to analyze the effects of soil layer and elevation on the Q10 values with a given proportion of cumulative C respired during the incubation and soil characteristics (SOC, TN, C:N ratio, pH value, and texture). The differences between the organic and the mineral layers at each site were also determined by paired t tests. The cumulative proportion of C respired and soil respiration rate during the incubation were analyzed by repeated-measures ANOVA and tested for sampling time, soil layer, elevation, and their interactions. All statistical analyses were performed by SAS 9.2 software (SAS Institute Inc., Cary, NC, USA).

Author Contributions

E.B. and W.X. designed the experiment. W.X. conducted the measurements, data analyses and wrote the manuscript. W.L., P.J. and H.W. assisted with the experiments. E.B. and W.X. reviewed the manuscript.