The permafrost zone is expected to be a substantial carbon source to the atmosphere, yet large-scale models currently only simulate gradual changes in seasonally thawed soil. Abrupt thaw will probably occur in <20% of the permafrost zone but could affect half of permafrost carbon through collapsing ground, rapid erosion and landslides. Here, we synthesize the best available information and develop inventory models to simulate abrupt thaw impacts on permafrost carbon balance. Emissions across 2.5 million km(2) of abrupt thaw could provide a similar climate feedback as gradual thaw emissions from the entire 18 million km(2) permafrost region under the warming projection of Representative Concentration Pathway 8.5. While models forecast that gradual thaw may lead to net ecosystem carbon uptake under projections of Representative Concentration Pathway 4.5, abrupt thaw emissions are likely to offset this potential carbon sink. Active hillslope erosional features will occupy 3% of abrupt thaw terrain by 2300 but emit one-third of abrupt thaw carbon losses. Thaw lakes and wetlands are methane hot spots but their carbon release is partially offset by slowly regrowing vegetation. After considering abrupt thaw stabilization, lake drainage and soil carbon uptake by vegetation regrowth, we conclude that models considering only gradual permafrost thaw are substantially underestimating carbon emissions from thawing permafrost.

Northern peatlands have accumulated large stocks of organic carbon (C) and nitrogen (N), but their spatial distribution and vulnerability to climate warming remain uncertain. Here, we used machine-learning techniques with extensive peat core data (n > 7,000) to create observation-based maps of northern peatland C and N stocks, and to assess their response to warming and permafrost thaw. We estimate that northern peatlands cover 3.7 ± 0.5 million km² and store 415 ± 150 Pg C and 10 ± 7 Pg N. Nearly half of the peatland area and peat C stocks are permafrost affected. Using modeled global warming stabilization scenarios (from 1.5 to 6 °C warming), we project that the current sink of atmospheric C (0.10 ± 0.02 Pg C⋅y⁻¹) in northern peatlands will shift to a C source as 0.8 to 1.9 million km² of permafrost-affected peatlands thaw. The projected thaw would cause peatland greenhouse gas emissions equal to ∼1% of anthropogenic radiative forcing in this century. The main forcing is from methane emissions (0.7 to 3 Pg cumulative CH₄-C) with smaller carbon dioxide forcing (1 to 2 Pg CO₂-C) and minor nitrous oxide losses. We project that initial CO₂-C losses reverse after ∼200 y, as warming strengthens peatland C-sinks. We project substantial, but highly uncertain, additional losses of peat into fluvial systems of 10 to 30 Pg C and 0.4 to 0.9 Pg N. The combined gaseous and fluvial peatland C loss estimated here adds 30 to 50% onto previous estimates of permafrost-thaw C losses, with southern permafrost regions being the most vulnerable.

Tundra soils account for 50% of global stocks of soil organic carbon (SOC), and it is expected that the amplified climate warming in high latitude could cause loss of this SOC through decomposition. Decomposed SOC could become hydrologically accessible, which increase downstream dissolved organic carbon (DOC) export and subsequent carbon release to the atmosphere, constituting a positive feedback to climate warming. However, DOC export is often neglected in ecosystem models. In this paper, we incorporate processes related to DOC production, mineralization, diffusion, sorption-desorption, and leaching into a customized arctic version of the dynamic ecosystem model LPJ-GUESS in order to mechanistically model catchment DOC export, and to link this flux to other ecosystem processes. The extended LPJ-GUESS is compared to observed DOC export at Stordalen catchment in northern Sweden. Vegetation communities include flood-tolerant graminoids (Eriophorum) and Sphagnum moss, birch forest and dwarf shrub communities. The processes, sorption-desorption and microbial decomposition (DOC production and mineralization) are found to contribute most to the variance in DOC export based on a detailed variance-based Sobol sensitivity analysis (SA) at grid cell-level. Catchment-level SA shows that the highest mean DOC exports come from the Eriophorum peatland (fen). A comparison with observations shows that the model captures the seasonality of DOC fluxes. Two catchment simulations, one without water lateral routing and one without peatland processes, were compared with the catchment simulations with all processes. The comparison showed that the current implementation of catchment lateral flow and peatland processes in LPJ-GUESS are essential to capture catchment-level DOC dynamics and indicate the model is at an appropriate level of complexity to represent the main mechanism of DOC dynamics in soils. The extended model provides a new tool to investigate potential interactions among climate change, vegetation dynamics, soil hydrology and DOC dynamics at both stand-alone to catchment scales.