ALBERT

Description: Soil is currently thought to be a sink for carbon; however, the response of this sink to increasing levels of atmospheric carbon dioxide and climate change is uncertain. In this study, we analyzed soil organic carbon (SOC) changes from 11 Earth system models (ESMs) under the historical and high radiative forcing (RCP 8.5) scenarios between 1850 and 2100. We used a reduced complexity model based on temperature and moisture sensitivities to analyze the drivers of SOC losses. ESM estimates of SOC change over the 21st century (2090–2099 minus 1997–2006) ranged from a loss of 72 Pg C to a gain 253 Pg C with a multi-model mean gain of 63 Pg C. All ESMs showed cumulative increases in both NPP (15% to 59%) and decreases in SOC turnover times (15% to 28%) over the 21st century. Most of the model-to-model variation in SOC change was explained by initial SOC stocks combined with the relative changes in soil inputs and decomposition rates (R2 = 0.88, p

Description: Stocks of soil organic carbon represent a large component of the carbon cycle that may participate in climate change feedbacks, particularly on decadal and century scales. For Earth system models (ESMs), the ability to accurately represent the global distribution of existing soil carbon stocks is a prerequisite for predicting future carbon-climate feedbacks. We compared soil carbon predictions from 16 ESMs to empirical data from the Harmonized World Soil Database (HWSD) and Northern Circumpolar Soil Carbon Database (NCSCD). Model estimates of global soil carbon stocks ranged from 510 to 3050 Pg C, compared to an estimate of 890–1660 Pg C from the HWSD. Model predictions for the high latitudes fell between 60 and 800 Pg C, compared to 380–620 Pg C from the NCSCD and 290 Pg C from the HWSD. This 5.3-fold variation in global soil carbon across models compared to a 3.4-fold variation in net primary productivity (NPP) and a 3.8-fold variation in global soil carbon turnover times. The spatial distribution of soil carbon predicted by the ESMs was not well correlated with the HWSD (Pearson's correlations 〈 0.4, RMSE 9.4 to 22.8 kg C m−2), although model-data agreement generally improved at the biome scale. There was poor agreement between the HWSD and NCSCD datasets in northern latitudes (Pearson's correlation = 0.33), indicating uncertainty in empirical estimates of soil carbon. We found that a reduced complexity model dependent on NPP and soil temperature explained most of the spatial variation in soil carbon predicted by most ESMs (R2 values between 0.73 and 0.93). This result suggests that differences in soil carbon predictions between ESMs are driven primarily by differences in predicted NPP and the parameterization of soil carbon responses to NPP and temperature not by structural differences between the models. Future work should focus on accurately representing these driving variables and modifying model structure to include additional processes.

Description: Stocks of soil organic carbon represent a large component of the carbon cycle that may participate in climate change feedbacks, particularly on decadal and centennial timescales. For Earth system models (ESMs), the ability to accurately represent the global distribution of existing soil carbon stocks is a prerequisite for accurately predicting future carbon–climate feedbacks. We compared soil carbon simulations from 11 model centers to empirical data from the Harmonized World Soil Database (HWSD) and the Northern Circumpolar Soil Carbon Database (NCSCD). Model estimates of global soil carbon stocks ranged from 510 to 3040 Pg C, compared to an estimate of 1260 Pg C (with a 95% confidence interval of 890–1660 Pg C) from the HWSD. Model simulations for the high northern latitudes fell between 60 and 820 Pg C, compared to 500 Pg C (with a 95% confidence interval of 380–620 Pg C) for the NCSCD and 290 Pg C for the HWSD. Global soil carbon varied 5.9 fold across models in response to a 2.6-fold variation in global net primary productivity (NPP) and a 3.6-fold variation in global soil carbon turnover times. Model–data agreement was moderate at the biome level (R2 values ranged from 0.38 to 0.97 with a mean of 0.75); however, the spatial distribution of soil carbon simulated by the ESMs at the 1° scale was not well correlated with the HWSD (Pearson correlation coefficients less than 0.4 and root mean square errors from 9.4 to 20.8 kg C m−2). In northern latitudes where the two data sets overlapped, agreement between the HWSD and the NCSCD was poor (Pearson correlation coefficient 0.33), indicating uncertainty in empirical estimates of soil carbon. We found that a reduced complexity model dependent on NPP and soil temperature explained much of the 1° spatial variation in soil carbon within most ESMs (R2 values between 0.62 and 0.93 for 9 of 11 model centers). However, the same reduced complexity model only explained 10% of the spatial variation in HWSD soil carbon when driven by observations of NPP and temperature, implying that other drivers or processes may be more important in explaining observed soil carbon distributions. The reduced complexity model also showed that differences in simulated soil carbon across ESMs were driven by differences in simulated NPP and the parameterization of soil heterotrophic respiration (inter-model R2 = 0.93), not by structural differences between the models. Overall, our results suggest that despite fair global-scale agreement with observational data and moderate agreement at the biome scale, most ESMs cannot reproduce grid-scale variation in soil carbon and may be missing key processes. Future work should focus on improving the simulation of driving variables for soil carbon stocks and modifying model structures to include additional processes.

Description: Soil is currently thought to be a sink for carbon; however, the response of this sink to increasing levels of atmospheric carbon dioxide and climate change is uncertain. In this study, we analyzed soil organic carbon (SOC) changes from 11 Earth system models (ESMs) contributing simulations to the Coupled Model Intercomparison Project Phase 5 (CMIP5). We used a reduced complexity model based on temperature and moisture sensitivities to analyze the drivers of SOC change for the historical and high radiative forcing (RCP 8.5) scenarios between 1850 and 2100. ESM estimates of SOC changed over the 21st century (2090–2099 minus 1997–2006) ranging from a loss of 72 Pg C to a gain of 253 Pg C with a multi-model mean gain of 65 Pg C. Many ESMs simulated large changes in high-latitude SOC that ranged from losses of 37 Pg C to gains of 146 Pg C with a multi-model mean gain of 39 Pg C across tundra and boreal biomes. All ESMs showed cumulative increases in global NPP (11 to 59%) and decreases in SOC turnover times (15 to 28%) over the 21st century. Most of the model-to-model variation in SOC change was explained by initial SOC stocks combined with the relative changes in soil inputs and decomposition rates (R2 = 0.89, p 〈 0.01). Between models, increases in decomposition rate were well explained by a combination of initial decomposition rate, ESM-specific Q10-factors, and changes in soil temperature (R2 = 0.80, p 〈 0.01). All SOC changes depended on sustained increases in NPP with global change (primarily driven by increasing CO2). Many ESMs simulated large accumulations of SOC in high-latitude biomes that are not consistent with empirical studies. Most ESMs poorly represented permafrost dynamics and omitted potential constraints on SOC storage, such as priming effects, nutrient availability, mineral surface stabilization, and aggregate formation. Future models that represent these constraints are likely to estimate smaller increases in SOC storage over the 21st century.