The enhanced vegetation growth by climate warming plays a pivotal role in amplifying the seasonal cycle of atmospheric CO 2 at northern lands (>50° N) since 1960s. However, the correlation between vegetation growth, temperature and seasonal amplitude of atmospheric CO 2 concentration have become elusive with the slowed increasing trend of vegetation growth and weakened temperature control on CO 2 uptake since late 1990s. Here, based on in situ atmospheric CO 2 concentration records from the Barrow observatory site, we found a slowdown in the increasing trend of the atmospheric CO 2 amplitude from 1990s to mid-2000s. This phenomenon was associated with the paused decrease in the minimum CO 2 concentration ([CO 2 ] min ), which was significantly correlated with the slowdown of vegetation greening and growing-season length extension. We then showed that both the vegetation greenness and growing-season length were positively correlated with spring but not autumn temperature over the northern lands. Furthermore, such asymmetric dependences of vegetation growth upon spring and autumn temperature cannot be captured by the state-of-art terrestrial biosphere models. These findings indicate that the responses of vegetation growth to spring and autumn warming are asymmetric, and highlight the need of improving autumn phenology in the models for predicting seasonal cycle of atmospheric CO 2 concentration.

A realistic simulation of snow cover and its thermal properties are important for accurate modelling of permafrost. We analyse simulated relationships between air and near-surface (20 cm) soil temperatures in the Northern Hemisphere permafrost region during winter, with a particular focus on snow insulation effects in nine land surface models, and compare them with observations from 268 Russian stations. There are large cross-model differences in the simulated differences between near-surface soil and air temperatures (ΔT; 3 to 14 °C), in the sensitivity of soil-to-air temperature (0.13 to 0.96 °C °C-1), and in the relationship between ΔT and snow depth. The observed relationship between ΔT and snow depth can be used as a metric to evaluate the effects of each model's representation of snow insulation, hence guide improvements to the model's conceptual structure and process parameterisations. Models with better performance apply multilayer snow schemes and consider complex snow processes. Some models show poor performance in representing snow insulation due to underestimation of snow depth and/or overestimation of snow conductivity. Generally, models identified as most acceptable with respect to snow insulation simulate reasonable areas of near-surface permafrost (13.19 to 15.77 million km2). However, there is not a simple relationship between the sophistication of the snow insulation in the acceptable models and the simulated area of Northern Hemisphere near-surface permafrost, because several other factors, such as soil depth used in the models, the treatment of soil organic matter content, hydrology and vegetation cover, also affect the simulated permafrost distribution.

We conducted a model-based assessment of changes in permafrost area and carbon storage for simulations driven by RCP4.5 and RCP8.5 projections between 2010 and 2299 for the northern permafrost region. All models simulating carbon represented soil with depth, a critical structural feature needed to represent the permafrost carbon-climate feedback, but that is not a universal feature of all climate models. Between 2010 and 2299, simulations indicated losses of permafrost between 3 and 5 million km² for the RCP4.5 climate and between 6 and 16 million km² for the RCP8.5 climate. For the RCP4.5 projection, cumulative change in soil carbon varied between 66-Pg C (10¹⁵-g carbon) loss to 70-Pg C gain. For the RCP8.5 projection, losses in soil carbon varied between 74 and 652 Pg C (mean loss, 341 Pg C). For the RCP4.5 projection, gains in vegetation carbon were largely responsible for the overall projected net gains in ecosystem carbon by 2299 (8- to 244-Pg C gains). In contrast, for the RCP8.5 projection, gains in vegetation carbon were not great enough to compensate for the losses of carbon projected by four of the five models; changes in ecosystem carbon ranged from a 641-Pg C loss to a 167-Pg C gain (mean, 208-Pg C loss). The models indicate that substantial net losses of ecosystem carbon would not occur until after 2100. This assessment suggests that effective mitigation efforts during the remainder of this century could attenuate the negative consequences of the permafrost carbon-climate feedback.

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m−2 yr−1), most models produced higher NPP (309 ± 12 g C m−2 yr−1) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m−2 yr−1), which mainly resulted from differences in simulated maximum monthly GPP (GPPmax). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C (Vcmax_25), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPPmax as well as their sensitivity to climate change.

A significant portion of the large amount of carbon (C) currently stored in soils of the permafrost region in the Northern Hemisphere has the potential to be emitted as the greenhouse gases CO2 and CH4 under a warmer climate. In this study we evaluated the variability in the sensitivity of permafrost and C in recent decades among land surface model simulations over the permafrost region between 1960 and 2009. The 15 model simulations all predict a loss of near-surface permafrost (within 3 m) area over the region, but there are large differences in the magnitude of the simulated rates of loss among the models (0.2 to 58.8 × 103 km2 yr−1). Sensitivity simulations indicated that changes in air temperature largely explained changes in permafrost area, although interactions among changes in other environmental variables also played a role. All of the models indicate that both vegetation and soil C storage together have increased by 156 to 954 Tg C yr−1 between 1960 and 2009 over the permafrost region even though model analyses indicate that warming alone would decrease soil C storage. Increases in gross primary production (GPP) largely explain the simulated increases in vegetation and soil C. The sensitivity of GPP to increases in atmospheric CO2 was the dominant cause of increases in GPP across the models, but comparison of simulated GPP trends across the 1982–2009 period with that of a global GPP data set indicates that all of the models overestimate the trend in GPP. Disturbance also appears to be an important factor affecting C storage, as models that consider disturbance had lower increases in C storage than models that did not consider disturbance. To improve the modeling of C in the permafrost region, there is the need for the modeling community to standardize structural representation of permafrost and carbon dynamics among models that are used to evaluate the permafrost C feedback and for the modeling and observational communities to jointly develop data sets and methodologies to more effectively benchmark models.