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. Significance Over many millennia, northern peatlands have accumulated large amounts of carbon and nitrogen, thus cooling the global climate. Over shorter timescales, peatland disturbances can trigger losses of peat and release of greenhouses gases. Despite their importance to the global climate, peatlands remain poorly mapped, and the vulnerability of permafrost peatlands to warming is uncertain. This study compiles over 7,000 field observations to present a data-driven map of northern peatlands and their carbon and nitrogen stocks. We use these maps to model the impact of permafrost thaw on peatlands and find that warming will likely shift the greenhouse gas balance of northern peatlands. At present, peatlands cool the climate, but anthropogenic warming can shift them into a net source of warming.

Scenarios that limit global warming to below 2 °C by 2100 assume significant land-use change to support large-scale carbon dioxide (CO₂) removal from the atmosphere by afforestation/reforestation, avoided deforestation, and Biomass Energy with Carbon Capture and Storage (BECCS). The more ambitious mitigation scenarios require even greater land area for mitigation and/or earlier adoption of CO₂ removal strategies. Here we show that additional land-use change to meet a 1.5 °C climate change target could result in net losses of carbon from the land. The effectiveness of BECCS strongly depends on several assumptions related to the choice of biomass, the fate of initial above ground biomass, and the fossil-fuel emissions offset in the energy system. Depending on these factors, carbon removed from the atmosphere through BECCS could easily be offset by losses due to land-use change. If BECCS involves replacing high-carbon content ecosystems with crops, then forest-based mitigation could be more efficient for atmospheric CO₂ removal than BECCS.

This paper describes the rationale and the protocol of the first component of the third simulation round of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a, http://www.isimip.org, last access: 2 November 2023) and the associated set of climate-related and direct human forcing data (CRF and DHF, respectively). The observation-based climate-related forcings for the first time include high-resolution observational climate forcings derived by orographic downscaling, monthly to hourly coastal water levels, and wind fields associated with historical tropical cyclones. The DHFs include land use patterns, population densities, information about water and agricultural management, and fishing intensities. The ISIMIP3a impact model simulations driven by these observation-based climate-related and direct human forcings are designed to test to what degree the impact models can explain observed changes in natural and human systems. In a second set of ISIMIP3a experiments the participating impact models are forced by the same DHFs but a counterfactual set of atmospheric forcings and coastal water levels where observed trends have been removed. These experiments are designed to allow for the attribution of observed changes in natural, human, and managed systems to climate change, rising CH₄ and CO₂ concentrations, and sea level rise according to the definition of the Working Group II contribution to the IPCC AR6.

Physical and chemical stabilisation mechanisms are now known to play a critical role in controlling carbon (C) storage in mineral soils, leading to suggestions that climate warming-induced C losses may be lower than previously predicted. By analysing > 9,000 soil profiles, here we show that, overall, C storage declines strongly with mean annual temperature. However, the reduction in C storage with temperature was more than three times greater in coarse-textured soils, with limited capacities for stabilising organic matter, than in fine-textured soils with greater stabilisation capacities. This pattern was observed independently in cool and warm regions, and after accounting for potentially confounding factors (plant productivity, precipitation, aridity, cation exchange capacity, and pH). The results could not, however, be represented by an established Earth system model (ESM). We conclude that warming will promote substantial soil C losses, but ESMs may not be predicting these losses accurately or which stocks are most vulnerable.

Understanding future changes in the terrestrial carbon cycle is important for reliable projections of climate change and impacts on ecosystems. It is well known that nitrogen (N) could limit plants' response to increased atmospheric carbon dioxide and it is therefore important to include a representation of the N cycle in Earth system models. Here we present the implementation of the terrestrial nitrogen cycle in the Joint UK Land Environment Simulator (JULES) – the land surface scheme of the UK Earth System Model (UKESM). Two configurations are discussed – the first one (JULES-CN) has a bulk soil biogeochemical model and the second one is a development configuration that resolves the soil biogeochemistry with depth (JULES-CN_layer). In JULES the nitrogen (N) cycle is based on the existing carbon (C) cycle and represents all the key terrestrial N processes in a parsimonious way. Biological N fixation is dependent on net primary productivity, and N deposition is specified as an external input. Nitrogen leaves the vegetation and soil system via leaching and a bulk gas loss term. Nutrient limitation reduces carbon-use efficiency (CUE – ratio of net to gross primary productivity) and can slow soil decomposition. We show that ecosystem level N limitation of net primary productivity (quantified in the model by the ratio of the potential amount of C that can be allocated to growth and spreading of the vegetation compared with the actual amount achieved in its natural state) falls at the lower end of the observational estimates in forests (approximately 1.0 in the model compared with 1.01 to 1.38 in the observations). The model shows more N limitation in the tropical savanna and tundra biomes, consistent with the available observations. Simulated C and N pools and fluxes are comparable to the limited available observations and model-derived estimates. The introduction of an N cycle improves the representation of interannual variability of global net ecosystem exchange, which was more pronounced in the C-cycle-only versions of JULES (JULES-C) than shown in estimates from the Global Carbon Project. It also reduces the present-day CUE from a global mean value of 0.45 for JULES-C to 0.41 for JULES-CN and 0.40 for JULES-CN_layer, all of which fall within the observational range. The N cycle also alters the response of the C fluxes over the 20th century and limits the CO₂ fertilisation effect, such that the simulated current-day land C sink is reduced by about 0.5 Pg C yr⁻¹ compared to the version with no N limitation. JULES-CN_layer additionally improves the representation of soil biogeochemistry, including turnover times in the northern high latitudes. The inclusion of a prognostic land N scheme marks a step forward in functionality and realism for the JULES and UKESM models.