Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

Description: The remaining carbon budget represents the total amount of CO2 that can still be emitted in the future while limiting global warming to a given temperature target. Remaining carbon budget estimates range widely, however, and this uncertainty can be used to either trivialize the most ambitious mitigation targets by characterizing them as impossible, or to argue that there is ample time to allow for a gradual transition to a low-carbon economy. Neither of these extremes is consistent with our best understanding of the policy implications of remaining carbon budgets. Understanding the scientific and socio-economic uncertainties affecting the size of the remaining carbon budgets, as well as the methodological choices and assumptions that underlie their calculation, is essential before applying them as a policy tool. Here we provide recommendations on how to calculate remaining carbon budgets in a traceable and transparent way, and discuss their uncertainties and implications for both international and national climate policies.

Description: The amount of additional future temperature change following a complete cessation of CO2 emissions is a measure of the unrealized warming to which we are committed due to CO2 already emitted to the atmosphere. This "zero emissions commitment" (ZEC) is also an important quantity when estimating the remaining carbon budget - a limit on the total amount of CO2 emissions consistent with limiting global mean temperature at a particular level. In the recent IPCC Special Report on Global Warming of 1.5 degrees C, the carbon budget framework used to calculate the remaining carbon budget for 1.5 degrees C included the assumption that the ZEC due to CO2 emissions is negligible and close to zero. Previous research has shown significant uncertainty even in the sign of the ZEC. To close this knowledge gap, we propose the Zero Emissions Commitment Model Intercomparison Project (ZECMIP), which will quantify the amount of unrealized temperature change that occurs after CO2 emissions cease and investigate the geophysical drivers behind this climate response. Quantitative information on ZEC is a key gap in our knowledge, and one that will not be addressed by currently planned CMIP6 simulations, yet it is crucial for verifying whether carbon budgets need to be adjusted to account for any unrealized temperature change resulting from past CO2 emissions. We request only one top-priority simulation from comprehensive general circulation Earth system models (ESMs) and Earth system models of intermediate complexity (EMICs) - a branch from the 1% CO2 run with CO2 emissions set to zero at the point of 1000 PgC of total CO2 emissions in the simulation - with the possibility for additional simulations, if resources allow. ZECMIP is part of CMIP6, under joint sponsorship by C4MIP and CDR-MIP, with associated experiment names to enable data submissions to the Earth System Grid Federation. All data will be published and made freely available.

Description: High-latitude soils store vast amounts of perennially frozen and therefore inert organic matter. With rising global temperatures and consequent permafrost degradation, a part of this carbon store will become available for microbial decay and eventual release to the atmosphere. We have developed a simplified, two-dimensional multi-pool model to estimate the strength and timing of future carbon dioxide (CO2) and methane (CH4) fluxes from newly thawed permafrost carbon (i.e. carbon thawed when temperatures rise above pre-industrial levels). We have especially simulated carbon release from deep deposits in Yedoma regions by describing abrupt thaw under thermokarst lakes. The computational efficiency of our model allowed us to run large, multi-centennial ensembles under various scenarios of future warming to express uncertainty inherent to simulations of the permafrost-carbon feedback. Under moderate warming of the representative concentration pathway (RCP) 2.6 scenario, cumulated CO2 fluxes from newly thawed permafrost carbon amount to 20 to 58 petagrammes of carbon (Pg-C) (68% range) by the year 2100 and reach 40 to 98 Pg-C in 2300. The much larger permafrost degradation under strong warming (RCP8.5) results in cumulated CO2 release of 42–141 and 157–313 Pg-C (68% ranges) in the years 2100 and 2300, respectively. Our estimates do only consider fluxes from newly thawed permafrost but not from soils already part of the seasonally thawed active layer under preindustrial climate. Our simulated methane fluxes contribute a few percent to total permafrost carbon release yet they can cause up to 40% of total permafrost-affected radiative forcing in the 21st century (upper 68% range). We infer largest methane emission rates of about 50 Tg-CH4 year–1 around the mid of the 21st century when simulated thermokarst lake extent is at its maximum and when abrupt thaw under thermokarst lakes is accounted for. CH4 release from newly thawed carbon in wetland-affected deposits is only discernible in the 22nd and 23rd century because of the absence of abrupt thaw processes. We further show that release from organic matter stored in deep deposits of Yedoma regions does crucially affect our simulated circumpolar methane fluxes. The additional warming through the release from newly thawed permafrost carbon proved only slightly dependent on the pathway of anthropogenic emission and amounts about 0.03–0.14 °C (68% ranges) by end of the century. The warming increased further in the 22nd and 23rd century and was most pronounced under the RCP6.0 scenario with adding 0.16–0.39°C (68% range) to simulated global mean surface air temperatures in the year 2300.

Description: High-latitude soils store vast amounts of perennially frozen and therefore inert organic matter. With rising global temperatures and consequent permafrost degradation, a part of this carbon stock will become available for microbial decay and eventual release to the atmosphere. We have developed a simplified, two-dimensional multi-pool model to estimate the strength and timing of future carbon dioxide (CO2) and methane (CH4) fluxes from newly thawed permafrost carbon (i.e. carbon thawed when temperatures rise above pre-industrial levels). We have especially simulated carbon release from deep deposits in Yedoma regions by describing abrupt thaw under newly formed thermokarst lakes. The computational efficiency of our model allowed us to run large, multi-centennial ensembles under various scenarios of future warming to express uncertainty inherent to simulations of the permafrost carbon feedback. Under moderate warming of the representative concentration pathway (RCP) 2.6 scenario, cumulated CO2 fluxes from newly thawed permafrost carbon amount to 20 to 58 petagrams of carbon (Pg-C) (68% range) by the year 2100 and reach 40 to 98 Pg-C in 2300. The much larger permafrost degradation under strong warming (RCP8.5) results in cumulated CO2 release of 42 to 141 Pg-C and 157 to 313 Pg-C (68% ranges) in the years 2100 and 2300, respectively. Our estimates only consider fluxes from newly thawed permafrost, not from soils already part of the seasonally thawed active layer under pre-industrial climate. Our simulated CH4 fluxes contribute a few percent to total permafrost carbon release yet they can cause up to 40% of total permafrost-affected radiative forcing in the 21st century (upper 68% range). We infer largest CH4 emission rates of about 50 Tg-CH4 per year around the middle of the 21st century when simulated thermokarst lake extent is at its maximum and when abrupt thaw under thermokarst lakes is taken into account. CH4 release from newly thawed carbon in wetland-affected deposits is only discernible in the 22nd and 23rd century because of the absence of abrupt thaw processes. We further show that release from organic matter stored in deep deposits of Yedoma regions crucially affects our simulated circumpolar CH4 fluxes. The additional warming through the release from newly thawed permafrost carbon proved only slightly dependent on the pathway of anthropogenic emission and amounts to about 0.03–0.14 °C (68% ranges) by end of the century. The warming increased further in the 22nd and 23rd century and was most pronounced under the RCP6.0 scenario, adding 0.16 to 0.39 °C (68% range) to simulated global mean surface air temperatures in the year 2300.

Description: With rising global temperatures and consequent permafrost degradation a part of old carbon stored in high latitude soils will become available for microbial decay and eventual release to the atmosphere. To estimate the strength and timing of future carbon dioxide and methane fluxes from newly thawed permafrost carbon, we have developed a simplified, two-dimensional multi-pool model. As large amounts of soil organic matter are stored in depths below three meters, we have also simulated carbon release from deep deposits in Yedoma regions. For this purpose we have modelled abrupt thaw under thermokarst lakes which can unlock large amounts of soil carbon buried deep in the ground. The computational efficiency of our 2-D model allowed us to run large, multi-centennial ensembles of differing scenarios of future warming to express uncertainty inherent to simulations of the permafrost-carbon feedback. Our model simulations, which are constrained by multiple lines of recent observations, suggest cumulated CO2 fluxes from newly thawed permafrost until the year 2100 of 20-58 Pg-C under moderate warming (RCP2.6), and of 42–141Pg-C under strong warming (RCP8.5). Under intense thermokarst activity, our simulated methane fluxes proved substantial and caused up to 40 % of total permafrost-affected radiative forcing in the 21st century. By quantifying CH4 contributions from different pools and depth levels, we discuss the role of thermokarst dynamics in affecting future Arctic carbon release. The additional global warming through the release from newly thawed permafrost carbon proved only slightly dependent on the pathway of anthropogenic emission in our simulations and reached about 0.1 °C by end of the century. The long-term, permafrost-affected global warming increased further in the 22nd and 23rd century, reaching a maximum of about 0.4°C in the year 2300.

Description: The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.