Large quantities of organic carbon are stored in frozen soils (permafrost) within Arctic and sub-Arctic regions. A warming climate can induce environmental changes that accelerate the microbial breakdown of organic carbon and the release of the greenhouse gases carbon dioxide and methane. This feedback can accelerate climate change, but the magnitude and timing of greenhouse gas emission from these regions and their impact on climate change remain uncertain. Here we find that current evidence suggests a gradual and prolonged release of greenhouse gas emissions in a warming climate and present a research strategy with which to target poorly understood aspects of permafrost carbon dynamics.

Abstract. Soils and other unconsolidated deposits in the northern circumpolar permafrost region store large amounts of soil organic carbon (SOC). This SOC is potentially vulnerable to remobilization following soil warming and permafrost thaw, but SOC stock estimates were poorly constrained and quantitative error estimates were lacking. This study presents revised estimates of permafrost SOC stocks, including quantitative uncertainty estimates, in the 0–3 m depth range in soils as well as for sediments deeper than 3 m in deltaic deposits of major rivers and in the Yedoma region of Siberia and Alaska. Revised estimates are based on significantly larger databases compared to previous studies. Despite this there is evidence of significant remaining regional data gaps. Estimates remain particularly poorly constrained for soils in the High Arctic region and physiographic regions with thin sedimentary overburden (mountains, highlands and plateaus) as well as for deposits below 3 m depth in deltas and the Yedoma region. While some components of the revised SOC stocks are similar in magnitude to those previously reported for this region, there are substantial differences in other components, including the fraction of perennially frozen SOC. Upscaled based on regional soil maps, estimated permafrost region SOC stocks are 217 ± 12 and 472 ± 27 Pg for the 0–0.3 and 0–1 m soil depths, respectively (±95% confidence intervals). Storage of SOC in 0–3 m of soils is estimated to 1035 ± 150 Pg. Of this, 34 ± 16 Pg C is stored in poorly developed soils of the High Arctic. Based on generalized calculations, storage of SOC below 3 m of surface soils in deltaic alluvium of major Arctic rivers is estimated as 91 ± 52 Pg. In the Yedoma region, estimated SOC stocks below 3 m depth are 181 ± 54 Pg, of which 74 ± 20 Pg is stored in intact Yedoma (late Pleistocene ice- and organic-rich silty sediments) with the remainder in refrozen thermokarst deposits. Total estimated SOC storage for the permafrost region is ∼1300 Pg with an uncertainty range of ∼1100 to 1500 Pg. Of this, ∼500 Pg is in non-permafrost soils, seasonally thawed in the active layer or in deeper taliks, while ∼800 Pg is perennially frozen. This represents a substantial ∼300 Pg lowering of the estimated perennially frozen SOC stock compared to previous estimates.

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and Key Points: • Updated Community Land Model has more hydrological and ecological process fidelity and more comprehensive representation of land management. • The model is systematically evaluated using International Land Model Benchmarking system and shows marked improvement over prior versions. parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5. Plain Language Summary The Community Land Model (CLM) is the land component of the widely used Community Earth System Model (CESM). Here, we introduce model developments included in CLM version 5 (CLM5), the default land component for CESM2 which will be used for the Coupled Model Intercomparison Project (CMIP6). CLM5 includes many new and updated processes including (1) hydrology and snow features such as spatially explicit soil depth, canopy snow processes, a simple firn model, and a more mechanistic river model, (2) plant hydraulics and hydraulic redistribution, (3) revised nitrogen cycling with flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake, (4) expansion to six crop types (global) and time-evolving irrigated areas and fertilization rates, (5) improved urban building energy model, and (6) carbon isotopes. New optional features include a demographically structured dynamic vegetat...

International audiencePermafrost soils contain enormous amounts of organic carbon, which could act as a positive feedback to global climate change due to enhanced respiration rates with warming. We have used a terrestrial ecosystem model that includes permafrost carbon dynamics, inhibition of respiration in frozen soil layers, vertical mixing of soil carbon from surface to permafrost layers, and CH4 emissions from flooded areas, and which better matches new circumpolar inventories of soil carbon stocks, to explore the potential for carbon-climate feedbacks at high latitudes. Contrary to model results for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), when permafrost processes are included, terrestrial ecosystems north of 60°N could shift from being a sink to a source of CO2 by the end of the 21st century when forced by a Special Report on Emissions Scenarios (SRES) A2 climate change scenario. Between 1860 and 2100, the model response to combined CO2 fertilization and climate change changes from a sink of 68 Pg to a 27 + -7 Pg sink to 4 + -18 Pg source, depending on the processes and parameter values used. The integrated change in carbon due to climate change shifts from near zero, which is within the range of previous model estimates, to a climate-induced loss of carbon by ecosystems in the range of 25 + -3 to 85 + -16 Pg C, depending on processes included in the model, with a best estimate of a 62 + -7 Pg C loss. Methane emissions from high-latitude regions are calculated to increase from 34 Tg CH4/y to 41-70 Tg CH4/y, with increases due to CO2 fertilization, permafrost thaw, and warming-induced increased CH4 flux densities partially offset by a reduction in wetland extent

No abstract

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO 2 emissions from fossil fuels and industry (E FF ) are based on energy statistics and cement production data, while emissions from land-use change (E LUC ), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO 2 concentration is measured directly and its rate of growth (G ATM ) is computed from the annual changes in concentration. The mean ocean CO 2 sink (S OCEAN ) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S OCEAN is evaluated with data products based on surveys of ocean CO 2 measurements. The global residual terrestrial CO 2 sink (S LAND ) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO 2 , and land-cover change (some including nitrogen-carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ , reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014), E FF was 9.0 ± 0.5 GtC yr −1 , E LUC was 0.9 ± 0.5 GtC yr −1 , G ATM was 4.4 ± 0.1 GtC yr −1 , S OCEAN was 2.6 ± 0.5 GtC yr −1 , and S LAND was 3.0 ± 0.8 GtC yr −1 . For the year 2014 alone, E FF grew to 9.8 ± 0.5 GtC yr −1 , 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr −1 that took place during 2005-2014. Also, for 2014, E LUC was 1.1 ± 0.5 GtC yr −1 , G ATM was 3.9 ± 0.2 GtC yr −1 , S OCEAN was 2.9 ± 0.5 GtC yr −1 , and S LAND was 4.1 ± 0.9 GtC yr −1 . G ATM was lower in 2014 compared to the past decade (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014), reflecting a larger S LAND for that year. The global atmospheric CO 2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in E FF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0....

Abstract. Soils are a crucial component of the Earth system; they comprise a large portion of terrestrial carbon stocks, mediate the supply and demand of nutrients, and influence the overall response of terrestrial ecosystems to perturbations. In this paper, we develop a new soil biogeochemistry model for the Community Land Model, version 4 (CLM4). The new model includes a vertical dimension to carbon (C) and nitrogen (N) pools and transformations, a more realistic treatment of mineral N pools, flexible treatment of the dynamics of decomposing carbon, and a radiocarbon ( 14 C) tracer. We describe the model structure, compare it with site-level and global observations, and discuss the overall effect of the revised soil model on Community Land Model (CLM) carbon dynamics. Site-level comparisons to radiocarbon and bulk soil C observations support the idea that soil C turnover is reduced at depth beyond what is expected from environmental controls for temperature, moisture, and oxygen that are considered in the model. In better agreement with observations, the revised soil model predicts substantially more and older soil C, particularly at high latitudes, where it resolves a permafrost soil C pool. In addition, the 20th century-C dynamics of the model are more realistic than those of the baseline model, with more terrestrial C uptake over the 20th century due to reduced N downregulation and longer turnover times for decomposing C.

We analyze global climate model predictions of soil temperature (from the Coupled Model Intercomparison Project 5 (CMIP5) database) to assess the models' representation of currentclimate soil thermal dynamics, and their predictions of permafrost thaw during the 21 st century. We compare the models' predictions to observations of active layer thickness, air temperature, and soil temperature, and to theoretically-expected relationships between active layer thickness and air temperature annual mean and seasonal cycle amplitude. Models show a wide range of current permafrost areas, active layer statistics (cumulative distributions, correlations with mean annual air temperature and amplitude of seasonal air temperature cycle), and ability to accurately model the coupling between soil and air temperatures at high latitudes. Many of the between-model differences can be traced to differences in the coupling between either near-surface air and shallow soil temperatures, or between shallow and deeper (1m) soil temperatures, which in turn reflect differences in snow physics and soil hydrology. We compare the models to observational datasets to benchmark the permafrostrelevant physics of the models. The models show a wide range of predictions for permafrost loss: 2-66% for RCP2.6, 15-87% for RCP4.5, and 30-99% for RCP8.5. Normalizing the amount of permafrost loss by the amount of high-latitude warming in the RCP4.5 scenario, the models predict an absolute loss of 1.6 ± 0.7 million km 2 permafrost • C −1 high-latitude warming, or a fractional loss of 6-29 % • C −1 .1