Description: Phytoplankton forms the base of the marine food web by transforming CO2 into organic carbon via photosynthesis. Some of the organic carbon is then transferred through the food web and exported into the deep ocean, a process known as the biological carbon pump. Despite the importance of phytoplankton for marine ecosystems and the global carbon cycle, projections of phytoplankton biomass in response to climate change differ strongly across Earth system models, illustrating uncertainty in our understanding of the underlying processes. Differences are especially large in the Southern Ocean, a region that is notoriously difficult to represent in models. Here, we argue that water column-integrated phytoplankton biomass in the Southern Ocean is projected to largely remain unchanged under climate change by the CMIP6 multi-model ensemble because of a shifting balance of bottom-up and top-down processes driven by a shoaling mixed layer depth. A shallower mixed layer is projected to improve growth conditions and consequently weaken bottom-up control. In addition to enhanced phytoplankton growth, the shoaling of the mixed layer also compresses phytoplankton closer to the surface and promotes zooplankton grazing efficiency, thus intensifying top-down control. Overall, our results suggest that while changes in bottom-up conditions stimulate enhanced growth, intensified top-down control opposes an increase in phytoplankton and becomes increasingly important for phytoplankton response under climate change in the Southern Ocean.

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: Numerical models are important tools for understanding the processes and feedbacks in the Earth system, including those involving changes in atmospheric CO2 (CO2,atm) concentrations. Here, we compile 55 published model studies (consisting of 778 individual simulations) that assess the impact of six forcing mechanisms on millennial-scale CO2,atm variations: changes in freshwater supply to the North Atlantic and Southern Ocean, the strength and position of the southern-hemisphere westerlies, Antarctic sea ice extent, and aeolian dust fluxes. We generally find agreement on the direction of simulated CO2,atm change across simulations, but the amplitude of change is inconsistent, primarily due to the different complexities of the model representation of Earth system processes. When freshwater is added to the North Atlantic, a reduced Atlantic Meridional Overturning Circulation (AMOC) is generally accompanied by an increase in Southern Ocean- and Pacific overturning, reduced Antarctic sea ice extent, spatially varying export production, and changes in carbon storage in the Atlantic (rising), in other ocean basins (generally decreasing) and on land (more varied). Positive or negative CO2,atm changes are simulated during AMOC minima due to a spatially and temporally varying dominance of individual terrestrial and oceanic drivers (and compensating effects between them) across the different models. In contrast, AMOC recoveries are often accompanied by rising CO2,atm levels, which are mostly driven by ocean carbon release (albeit from different regions). The magnitude of simulated CO2,atm rise broadly scales with the duration of the AMOC perturbation (i.e., the stadial length). When freshwater is added to the Southern Ocean, reduced deep-ocean ventilation drives a CO2,atm drop via reduced carbon release from the Southern Ocean. Although the impacts of shifted southern-hemisphere westerlies are inconsistent across model simulations, their intensification raises CO2,atm via enhanced Southern Ocean Ekman pumping. Increased supply of aeolian dust to the ocean, and thus iron fertilisation of marine productivity, consistently lowers modelled CO2,atm concentrations via more efficient nutrient utilisation. The magnitude of CO2,atm change in response to dust flux variations, however, largely depends on the complexity of models' marine ecosystem and iron cycle. This especially applies to simulations forced by Antarctic sea ice changes, in which the direction of simulated CO2,atm change varies greatly across model hierarchies. Our compilation highlights that no single (forcing) mechanism can explain observed past millennial-scale CO2,atm variability, and identifies important future needs in coupled carbon cycle-climate modelling to better understand the mechanisms governing CO2,atm changes in the past.

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 Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

Description: Numerical models are important tools for understanding the processes and feedbacks in the Earth system, including those involving changes in atmospheric CO2 (CO2,atm) concentrations. Here, we compile 55 published model studies (consisting of 778 individual simulations) that assess the impact of six forcing mechanisms on millennial-scale CO2,atm variations: changes in freshwater supply to the North Atlantic and Southern Ocean, the strength and position of the southern-hemisphere westerlies, Antarctic sea ice extent, and aeolian dust fluxes. We generally find agreement on the direction of simulated CO2,atm change across simulations, but the amplitude of change is inconsistent, primarily due to the different complexities of the model representation of Earth system processes. When freshwater is added to the North Atlantic, a reduced Atlantic Meridional Overturning Circulation (AMOC) is generally accompanied by an increase in Southern Ocean- and Pacific overturning, reduced Antarctic sea ice extent, spatially varying export production, and changes in carbon storage in the Atlantic (rising), in other ocean basins (generally decreasing) and on land (more varied). Positive or negative CO2,atm changes are simulated during AMOC minima due to a spatially and temporally varying dominance of individual terrestrial and oceanic drivers (and compensating effects between them) across the different models. In contrast, AMOC recoveries are often accompanied by rising CO2,atm levels, which are mostly driven by ocean carbon release (albeit from different regions). The magnitude of simulated CO2,atm rise broadly scales with the duration of the AMOC perturbation (i.e., the stadial length). When freshwater is added to the Southern Ocean, reduced deep-ocean ventilation drives a CO2,atm drop via reduced carbon release from the Southern Ocean. Although the impacts of shifted southern-hemisphere westerlies are inconsistent across model simulations, their intensification raises CO2,atm via enhanced Southern Ocean Ekman pumping. Increased supply of aeolian dust to the ocean, and thus iron fertilisation of marine productivity, consistently lowers modelled CO2,atm concentrations via more efficient nutrient utilisation. The magni- tude of CO2,atm change in response to dust flux variations, however, largely depends on the complexity of models' marine ecosystem and iron cycle. This especially applies to simulations forced by Antarctic sea ice changes, in which the direction of simulated CO2,atm change varies greatly across model hierarchies. Our compilation highlights that no single (forcing) mechanism can explain observed past millennial-scale CO2,atm variability, and identifies important future needs in coupled carbon cycle-climate modelling to better understand the mechanisms governing CO2,atm changes in the past.

Description: The methods to quantify equilibrium climate sensitivity are still debated. We collect millennial-length simulations of coupled climate models and show that the global mean equilibrium warming is higher than those obtained using extrapolation methods from shorter simulations. Specifically, 27 simulations with 15 climate models forced with a range of CO2 concentrations show a median 17% larger equilibrium warming than estimated from the first 150 years of the simulations. The spatial patterns of radiative feedbacks change continuously, in most regions reducing their tendency to stabilizing the climate. In the equatorial Pacific, however, feedbacks become more stabilizing with time. The global feedback evolution is initially dominated by the tropics, with eventual substantial contributions from the mid-latitudes. Time-dependent feedbacks underscore the need of a measure of climate sensitivity that accounts for the degree of equilibration, so that models, observations, and paleo proxies can be adequately compared and aggregated to estimate future warming.