Mechanisms of millennial-scale atmospheric CO2 change in numerical model simulations

Peter.Koehler [ at ]


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.

Item Type
Primary Division
Primary Topic
Publication Status
Eprint ID
DOI 10.1016/j.quascirev.2019.05.013

Cite as
Gottschalk, J. , Battaglia, G. , Fischer, H. , Frölicher, T. L. , Jaccard, S. L. , Jeltsch-Thömmes, A. , Joos, F. , Köhler, P. , Meissner, K. J. , Menviel, L. , Nehrbass-Ahles, C. , Schmitt, J. , Schmittner, A. , Skinner, L. C. and Stocker, T. F. (2019): Mechanisms of millennial-scale atmospheric CO2 change in numerical model simulations , Quaternary Science Reviews, 220 , pp. 30-74 . doi: 10.1016/j.quascirev.2019.05.013

[thumbnail of gottschalk2019qsr.pdf]

Download (11MB) | Preview

Add to AnyAdd to TwitterAdd to FacebookAdd to LinkedinAdd to PinterestAdd to Email


Geographical region

Research Platforms


Edit Item Edit Item