Biome Q10 and Dryness

Temperature sensitivity of soil respiration (Q10) is a critical parameter in carbon cycle models with important implications for climate-carbon feedbacks in the 21st century. The common assumption of a constant Q10, usually with a value of 2.0, was shown to be invalid by a previous model-data fusion study that reported biome-specific values of this parameter. We extend the previous analysis by demonstrating that these biome-level values of Q10 also are a function of dryness (R 2 = 0.54). When tundra and cultivated lands are excluded, the correlation is much stronger (R 2 = 0.92). Therefore dryness is the primary driver for variability in respiration-temperature sensitivity in forest and grassland ecosystems. This finding has important implications for the response of the terrestrial carbon cycle to climate change, as it implies that the increasing dryness would potentially accelerate the respiration temperature sensitivity feedback.


Introduction
Globally, soil respiration releases CO 2 annually at a rate that is over an order of magnitude larger than anthropogenic releases [1]. Although soil heterotrophic respiration is currently balanced or slightly exceeded by terrestrial net primary productivity (NPP), relatively small changes in this large flux could have large impacts on the global net carbon balance. The most important climate driver of soil respiration is temperature, and increasing temperature is likely to induce a positive feedback between climate and the carbon cycle. Uncertainty about the strength of this feedback is a primary source of uncertainty for predicted behaviour of terrestrial carbon sinks in the latter half of this century [2,3]. Additionally, this sensitivity of respiration to temperature was shown to vary as functions of temperature, substrate, soil moisture and/or biome [4,5]. Despite this, many global carbon cycle models assume constant temperature sensitivity, or they assume a sensitivity that depends on only a limited subset of these factors.
Here we extend the analysis of Zhou et al. [4], who used an inversion approach to assimilate worldwide soil respiration measurements and measured soil organic carbon into a widely used carbon cycle model. The authors of that study concluded that Q 10 is a function of biome, and that the assumption of a constant Q 10 results in an underestimation of the respiration-temperature feedback intensity by 25%. Upon further analysis of this unique dataset, we find that the respiration-temperature sensitivity is also a strong function of dryness at the biome level. This finding has important implications for the behaviour of the carbon cycle in a changing climate: increasing dryness, which is likely in a warming climate [6], may increase the respiration-temperature sensitivity, accelerating decomposition and providing a stronger feedback to the climate system.

Methods
This analysis focuses on Q 10 , the parameter controlling the temperature-dependence of soil respiration in the following way: where   R T and ref R are soil respiration at measured temperature (T) and reference temperature ( ref T ), respectively, and 10 Q is a factor by which respiration is multiplied when temperature is increased by 10˚C. At Q 10 = 1 respiration would be independent of temperature, while larger Q 10 values indicates a stronger temperature dependence. In many ecosystem models, Q 10 is treated as a constant, the most common value being 2. However, considerable variation (1.3 to ~10) in Q 10 values have been reported by numerous investigations [7][8][9]. Since van't Hoff introduced Q 10 in 1898 it has been debated weather Q 10 is a universal constant and what controls Q 10 [1][2][3][4][5][6][7][8][9]. Since soil emissions of CO 2 are expected to have a positive feedback on global warming, modelling of climate change, its consequences and control strategies requires a clearer understanding of the Q 10 of soils.
We examined climate control of Q 10 at the level of biomes globally. Spatially resolved Q 10 values were estimated at a resolution of 1˚ by 1˚ using a model-data fusion technique to assimilate worldwide soil respiration measurements and measured soil organic carbon [4]. The model used in that analysis was the Carnegie-Ames-Stanford Approach (CASA) model [10,11], which includes a CENTURY-based soil carbon module that simulates soil organic carbon processes using two carbon storage pools, in addition to litter and microbial pools. Following Zhou et al. [4], we averaged the estimated Q 10 values for each biome except Desert and Shrub & Bare Ground, both of which are subject to prolonged periods of desiccation. Average temperature, precipitation and net radiation were estimated for each 1˚ by 1˚ grid cell and these climate variables were used to estimate biome dryness [12]: where n R (MJ·m −2 ·yr −1 ) is an annual sum of net radiation, L (2.5 MJ·kg −1 ) is a latent heat coefficient and P (mm·yr −1 ) is the total annual precipitation. We then linked the biome-level Q 10 values to biome-level climate and dryness data ( Table 1).

Results and Discussion
We found that biome-level Q 10 is significantly correlated to dryness (R 2 = 0.54, Figure 1) or to precipitation (R 2 = 0.45, Figure omitted) and is independent of net radiation and temperature. We expect that the correlation of Q 10 with dryness is stronger than with precipitation because the soil moisture content is determined not only by precipitation (input), but also by energy available for evaporation (output). In our analysis, both tundra and cultivated soils appear to be outliers. When these systems are excluded, the correlation of biome-level Q 10 with dryness is much higher (R 2 = 0.92).
It is not surprising that the Q 10 value for agricultural soil is an outlier (Figure 1) because cultivation (including tillage, fertilization, irrigation and drainage) accelerates soil respiration in ways that are not adequately captured by CASA [7]. The high Q 10 value for tundra may be a consequence of the non-linearity of respiration with respect to temperature, particularly as T approaches 0˚C. In systems with permafrost, which are poorly represented by CASA, the high Q 10 value is also likely related to depth of the active layer. As the soil column warms, the depth of unfrozen soil increases and exposes more soil organic matter to decomposition [13]. In actuality this increases the base respiration rate, but the strong correlation of this effect with temperature causes it to be interpreted by the model-data fusion technique as a higher Q 10 value. Thus both cultivated soil and tundra Q 10 values are as high as grasslands even though the dryness of these systems is similar to deciduous forests.
Our results demonstrate that the temperature sensitivity of soil respiration-aggregated at the biome level, is controlled by dryness, and that soils in different biomes respond differently to dryness. As can be seen in Figure  1, the temperature sensitivity of forest soil respiration to dryness is much less than that of grasslands, which may also related to substrate quality (woody and non-woody components). This implies that conversion of forest to pasture and agriculture would increase temperature sensitivity of soil respiration, accelerating CO 2 emissions Table 1
and global warming [14,15]. Furthermore, we reiterate the findings of Zhou et al. [4] that models of soil responses to climate change, at least at the biome level, should not assume Q 10 = 2, but need to accommodate the sensitivity of soil respiration in different soil types to dryness. Thus, the feedback to terrestrial ecosystems from respiration temperature sensitivity may accelerate more than previously predicted if climate change causes an increase in dryness. Uncertainties may result from the methods of Zhou et al. [4], in which the International Geosphere-Biosphere Programme (IGBP) soil organic carbon (SOC) dataset was used to optimize the Q 10 parameter in the CASA model. A significant part of the total SOC is relatively recalcitrant and resistant to decomposition, while soil respiration is mostly produced from the newly shed plant litter, and surface SOC. This may lead to an underestimate of the optimized Q 10 because the decomposition of the recalcitrant SOC is slower than assumed by the model [16]. However, at the global scale, this underestimation of the optimized Q 10 is mitigated because Zhou et al. [4] constrained the optimization process such that the global mean optimized Q 10 is equal to the global mean value of Q 10 (1.72) of soil respiration measurements in major ecosystems of the world as reported by Raich et al. [1]. While local biases in Q 10 may result if the modeled SOC quality is not correct, comparisons against site-level soil respiration data in Zhou et al. [4] show consistent improvement when optimized values are used compared to a globally constant value. Theoretically, the biome Q 10 reflects apparent temperature sensitivity that is controlled by environmental constraints [16]. Our results demonstrate that dryness is the most important control on Q 10 among environmental constraints at biome level.