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: Carbon cycle models used to calculate the marine reservoir age of the non-polar surface ocean (called Marine20) out of IntCal20, the compilation of atmospheric C, have so far neglected a key aspect of the millennial-scale variability connected with the thermal bipolar seesaw: changes in the strength of the Atlantic meridional overturning circulation (AMOC) related to Dansgaard/Oeschger and Heinrich events. Here we implement such AMOC changes in the carbon cycle box model BICYCLE-SE to investigate how model performance over the last 55 kyr is affected, in particular with respect to available 14C and CO2 data. Constraints from deep ocean 14C data suggest that the AMOC in the model during Heinrich stadial 1 needs to be highly reduced or even completely shutdown. Ocean circulation and sea ice coverage combined are the processes that almost completely explain the simulated changes in deep ocean 14C age, and these are also responsible for a glacial drawdown of ∼60 ppm of atmospheric CO2. We find that the implementation of abrupt reductions in AMOC during Greenland stadials in the model setup that was previously used for the calculation of Marine20 leads to differences of less than ±100 14C yrs. The representation of AMOC changes therefore appears to be of minor importance for deriving non-polar mean ocean radiocarbon calibration products such as Marine20, where atmospheric carbon cycle variables are forced by reconstructions. However, simulated atmospheric CO2 exhibits minima during AMOC reductions in Heinrich stadials, in disagreement with ice core data. This mismatch supports previous suggestions that millennial-scale changes in CO2 were probably not driven directly by the AMOC, but rather by biological and physical processes in the Southern Ocean and by contributions from variable land carbon storage.

Description: The concentration of radiocarbon (14C) differs between ocean and atmosphere. Radiocarbon determinations from samples which obtained their 14C in the marine environment therefore need a marine-specific calibration curve and cannot be calibrated directly against the atmospheric-based IntCal20 curve. This paper presents Marine20, an update to the internationally agreed marine radiocarbon age calibration curve that provides a non-polar global-average marine record of radiocarbon from 0–55 cal kBP and serves as a baseline for regional oceanic variation. Marine20 is intended for calibration of marine radiocarbon samples from non-polar regions; it is not suitable for calibration in polar regions where variability in sea ice extent, ocean upwelling and air-sea gas exchange may have caused larger changes to concentrations of marine radiocarbon. The Marine20 curve is based upon 500 simulations with an ocean/atmosphere/biosphere box-model of the global carbon cycle that has been forced by posterior realizations of our Northern Hemispheric atmospheric IntCal20 14C curve and reconstructed changes in CO2 obtained from ice core data. These forcings enable us to incorporate carbon cycle dynamics and temporal changes in the atmospheric 14C level. The box-model simulations of the global-average marine radiocarbon reservoir age are similar to those of a more complex three-dimensional ocean general circulation model. However, simplicity and speed of the box model allow us to use a Monte Carlo approach to rigorously propagate the uncertainty in both the historic concentration of atmospheric 14C and other key parameters of the carbon cycle through to our final Marine20 calibration curve. This robust propagation of uncertainty is fundamental to providing reliable precision for the radiocarbon age calibration of marine based samples. We make a first step towards deconvolving the contributions of different processes to the total uncertainty; discuss the main differences of Marine20 from the previous age calibration curve Marine13; and identify the limitations of our approach together with key areas for further work. The updated values for ΔR, the regional marine radiocarbon reservoir age corrections required to calibrate against Marine20, can be found at the data base http://calib.org/marine/.

Description: Highlights • Strong ocean stratification in the Fram Strait during the late glacial (33–26 ka). • Breakup of ocean stratification during the LGM (26–20 ka) due to enhanced upwelling. • No extreme aging of 〉6000 yr in the Arctic Mediterranean as previously suggested. Abstract The present-day ocean ventilation in the Arctic Mediterranean (Nordic Seas and Arctic Ocean), via transformation of northward inflowing warm Atlantic surface water into cold deep water, affects regional climate, atmospheric circulation and carbon storage in the deep ocean. Here we study the glacial evolution of the Arctic Mediterranean circulation and its influence on glacial climate using radiocarbon reservoir-age reconstructions on deep-sea cores from the Fram Strait that cover the late glacial period (33,000–20,000 yr ago; 33–20 ka). Our results show high Benthic-Planktic 14C age differences of ∼1500 14C years 33–26.5 ka suggesting significant water column stratification between ∼100–2600 m water depth, and reduction and/or shoaling of deep-water formation. This phase was followed by break-up of the stratification during the Last Glacial Maximum (LGM; 26–20 ka), with Benthic-Planktic 14C age differences of ∼250 14C years, likely due to enhanced upwelling. These ocean circulation changes potentially contributed to the final intensification phase of glaciation via positive cryosphere-atmosphere-ocean circulation-carbon cycle feedbacks. Our data also do not support ‘extreme aging’ of 〉6000 14C years in the deep Arctic Mediterranean, and appear to rule out the proposed outflow of very old Arctic Ocean water to the Nordic Seas during the LGM and to the subpolar North Atlantic Ocean during the deglacial period.

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.