The 1991 review “Ocean circulation and climate change” by Klaus Hasselmann (Tellus B 43:82–103) closed with five unresolved issues spanning coupled ocean–atmosphere modeling: alternative Younger Dryas forcings (freshwater routing, winds) untested in coupled simulations; ~century-long warming integrations too short to separate transient ocean heat uptake from equilibrium climate change, including possible unstable transitions between circulation equilibria; unvalidated acceleration methods such as burst integration for long ocean time scales; the absence of an operational combined wind–wave data-assimilation system; and the inability of ocean carbon-cycle models to partition paleoclimatic CO2 change uniquely among competing mechanisms. Subsequent work has substantially addressed the first four (coupled deglacial transient simulations; multi-centennial integrations and standardized transient/equilibrium diagnostics; validated acceleration and burst-coupling methods; operational wave and coupled ocean–atmosphere data assimilation), while the fifth survives essentially unchanged: the quantitative mechanistic partitioning of the roughly 80–100 ppmv glacial–interglacial atmospheric CO2 change among ocean carbon-cycle mechanisms remains unresolved, as current assessments confirm. The task below addresses that surviving open core.

Background

The ocean is the slow, heat- and carbon-storing component of the climate system on time scales from years to millennia, and its circulation is controlled by exchanges of buoyancy (heat and freshwater) with the atmosphere, particularly where deep water forms at high latitudes. An ocean general circulation model (OGCM) solves the fluid equations for this circulation globally; when an OGCM is coupled to an atmospheric general circulation model, the resulting coupled model historically suffered from climate drift (an unrealistic slow evolution of the simulated state), often suppressed by flux correction (artificial adjustments of the surface heat and freshwater fluxes). Hasselmann's 1991 review described experiments with the Hamburg Large-Scale Geostrophic (LSG) OGCM showing that modest changes in high-latitude surface buoyancy forcing shift the simulated global circulation between regimes dominated by North Atlantic and by Antarctic deep-water formation — the sensitivity that motivates this task.

The review listed five unresolved issues. First, its Younger Dryas experiments — the Younger Dryas (YD) being the abrupt cold interval roughly 12,900–11,700 years ago that interrupted the last deglaciation, generally attributed to a weakening of the Atlantic meridional overturning circulation (AMOC, the global-scale overturning driven by deep-water formation in the North Atlantic and around Antarctica) — were ocean-only simulations with prescribed meltwater pulses (measured in Sverdrups; 1 Sv = one million cubic meters per second), so alternative causes such as the spatial distribution of freshwater input or wind changes remained untested in coupled models. Second, coupled warming integrations under the IPCC business-as-usual scenario covered only about 100 years, too short to separate the transient delay of warming by ocean heat uptake from changes in the asymptotic equilibrium climate state, including possible unstable transitions to a different circulation equilibrium. Third, extending coupled simulations over the ocean's full dynamic time-scale range required acceleration methods — techniques that advance slow ocean tracers with distorted time steps to save computer time — of which burst integration was still unvalidated. Fourth, a fully operational analysis system combining atmospheric modeling, wave modeling and satellite wind and wave observations was yet to be implemented and validated. Fifth, the Hamburg three-dimensional ocean carbon-cycle model, validated against the GEOSECS chemical-tracer survey of the world ocean, reproduced the paleoclimatic carbon record only when several model parameters were tuned jointly, leaving the mechanistic partition of past atmospheric CO2 changes unresolved. This fifth thread is the surviving open core and is stated in the problem below.

The carbon-cycle background needed for that core: atmospheric CO2 is controlled on glacial time scales mainly by the ocean, which contains about sixty times more carbon than the atmosphere and land biota combined. The biological (soft-tissue) pump transfers carbon downward through the settling of organic matter produced by plankton; the carbonate pump is the analogous transfer by calcium-carbonate shells, which also alters alkalinity (the ocean's capacity to hold dissolved inorganic carbon, restored over millennia by dissolution of deep-sea sediments, a process called carbonate compensation); the solubility pump is the temperature-driven dissolution of CO2 in cold, sinking waters. Nutrients carried into the ocean interior without being used by biology are called preformed, while those regenerated in the interior by remineralization are called regenerated; the preformed fraction measures how completely the biological pump harvests the ocean's nutrient supply. Ventilation is the rate at which surface waters are exchanged with the deep ocean, which sets how quickly deep carbon re-equilibrates with the atmosphere. Iron fertilization refers to the limiting role of wind-borne, dust-delivered iron in the nutrient-rich but iron-poor Southern Ocean. Ice-core records show that glacial-period atmospheric CO2 was roughly 80–100 ppmv lower than the preindustrial value of about 280 ppmv (parts per million by volume).

In the model experiments reviewed in 1991 (using the transport scheme of Maier-Reimer and Hasselmann, 1987), sensitivity tests in which individual parameters were varied separately showed that changes in deep-ocean ventilation produced the largest single modeled contribution to atmospheric CO2, but no single parameter change reproduced the suite of paleoclimatic carbon-cycle indices; only a tuned, time-dependent combination of several parameter changes could. That is the paper's own statement of the partitioning problem.

Three decades of later work have refined but not closed it. The candidate mechanisms are now standard: Southern Ocean stratification and ventilation change, sea-ice-limited air–sea gas exchange, iron-fertilization-driven nutrient utilization, solubility change from temperature and salinity, whole-ocean alkalinity change with carbonate compensation, and the boundary terms of terrestrial biosphere carbon storage and weathering fluxes (Sigman and Boyle, 2000; Kohfeld and Ridgwell, 2009; Fischer et al., 2010; Sigman et al., 2010). Mechanism-specific quantifications exist — the deglacial iron-fertilization contribution bounded by Bauska et al. (2016), buoyancy-forced ventilation changes quantified by Watson et al. (2015), preformed-nutrient proxy evidence for pump efficiency by Vollmer et al. (2022) — and the mechanisms interact nonlinearly, so their combined CO2 effect is not the sum of isolated contributions (Hain et al., 2010). The most recent assessments state that identifying which mechanisms dominated remains unresolved (Sigman and Hain, 2024).

Problem Statement

Question. Determine the quantitative mechanistic partition of the glacial–interglacial change in atmospheric CO2 — the roughly 80–100 ppmv lower concentration during the Last Glacial Maximum (about 20,000 years ago) relative to the preindustrial level, together with its time evolution over the last glacial cycle and deglaciation — among the competing mechanisms: (i) changes in deep-ocean ventilation and stratification, especially of the Southern Ocean and other high-latitude deep-water formation regions; (ii) changes in the efficiency of the biological pump, including iron-fertilization-driven nutrient utilization; (iii) changes in air–sea gas exchange, for example through sea-ice coverage; (iv) solubility changes from ocean temperature and salinity; (v) whole-ocean alkalinity change with carbonate compensation; and (vi) the boundary contributions of terrestrial biosphere carbon storage and weathering fluxes.

The quantitative partition assigns each mechanism an amplitude contribution with uncertainty under an explicitly declared decomposition convention; it reproduces the observed ice-core CO2 time series and the co-varying carbon-cycle indices — atmospheric and oceanic carbon-isotope ratios (δ¹³C), nutrient- and oxygen-based proxies of preformed versus regenerated carbon, radiocarbon-based ventilation ages, and carbonate-system proxies — within their observational uncertainties; and it excludes the principal alternative partitions rather than exhibiting one tuned admissible combination of parameter changes.

References

  1. Klaus Hasselmann (1991). Ocean circulation and climate change. Tellus B 43(4), 82–103. DOI: 10.1034/j.1600-0889.1991.t01-2-00008.x (also published in Tellus A 43(4); open-access version DOI 10.3402/tellusa.v43i4.11939). https://doi.org/10.1034/j.1600-0889.1991.t01-2-00008.x
  2. E. Maier-Reimer and K. Hasselmann (1987). Transport and storage of CO2 in the ocean — an inorganic ocean-circulation carbon cycle model. Climate Dynamics 2(2), 63–90. DOI: 10.1007/BF01054491. https://doi.org/10.1007/BF01054491
  3. Uwe Mikolajewicz, Thomas J. Crowley, Andreas Schiller, Reinhard Voss (1997). Modelling teleconnections between the North Atlantic and North Pacific during the Younger Dryas. Nature 387, 384–387. DOI: 10.1038/387384a0. https://doi.org/10.1038/387384a0
  4. Peter U. Clark, Shawn J. Marshall, Garry K. C. Clarke, Steven W. Hostetler, Joseph M. Licciardi, James T. Teller (2001). Freshwater forcing of abrupt climate change during the last glaciation. Science 293(5528), 283–287. DOI: 10.1126/science.1062517. https://doi.org/10.1126/science.1062517
  5. Z. Liu, B. L. Otto-Bliesner, F. He, E. C. Brady, R. Tomas, P. U. Clark, A. E. Carlson, J. Lynch-Stieglitz, W. Curry, E. Brook, D. Erickson, R. Jacob, J. Kutzbach, J. Cheng (2009). Transient simulation of last deglaciation with a new mechanism for Bølling-Allerød warming. Science 325(5938), 310–314. DOI: 10.1126/science.1171041. https://doi.org/10.1126/science.1171041
  6. S. Manabe and R. J. Stouffer (1988). Two stable equilibria of a coupled ocean-atmosphere model. Journal of Climate 1(9), 841–866. DOI: 10.1175/1520-0442(1988)001<0841:TSEOAC>2.0.CO;2. https://doi.org/10.1175/1520-0442(1988)001%3C0841%3ATSEOAC%3E2.0.CO%3B2
  7. S. Manabe and R. J. Stouffer (1994). Multiple-century response of a coupled ocean-atmosphere model to an increase of atmospheric carbon dioxide. Journal of Climate 7(1), 5–23. DOI: 10.1175/1520-0442(1994)007<0005:MCROAC>2.0.CO;2. https://doi.org/10.1175/1520-0442(1994)007%3C0005%3AMCROAC%3E2.0.CO%3B2
  8. R. J. Stouffer and S. Manabe (1999). Response of a coupled ocean–atmosphere model to increasing atmospheric carbon dioxide: sensitivity to the rate of increase. Journal of Climate 12(8), 2224–2237. DOI: 10.1175/1520-0442(1999)012<2224:ROACOA>2.0.CO;2. https://doi.org/10.1175/1520-0442(1999)012%3C2224%3AROACOA%3E2.0.CO%3B2
  9. J. M. Gregory, W. J. Ingram, M. A. Palmer, G. S. Jones, P. A. Stott, R. B. Thorpe, J. A. Lowe, T. C. Johns, K. D. Williams (2004). A new method for diagnosing radiative forcing and climate sensitivity. Geophysical Research Letters 31(3), L03205. DOI: 10.1029/2003GL018747. https://doi.org/10.1029/2003GL018747
  10. Stefan Rahmstorf (1995). Bifurcations of the Atlantic thermohaline circulation in response to changes in the hydrological cycle. Nature 378, 145–149. DOI: 10.1038/378145a0. https://doi.org/10.1038/378145a0
  11. E. Hawkins, R. S. Smith, L. C. Allison, J. M. Gregory, T. J. Woollings, H. Pohlmann, B. de Cuevas (2011). Bistability of the Atlantic overturning circulation in a global climate model and links to ocean freshwater transport. Geophysical Research Letters 38(10), L10605. DOI: 10.1029/2011GL047208. https://doi.org/10.1029/2011GL047208
  12. Peter Ditlevsen and Susanne Ditlevsen (2023). Warning of a forthcoming collapse of the Atlantic meridional overturning circulation. Nature Communications 14, 4254. DOI: 10.1038/s41467-023-39810-w. https://doi.org/10.1038/s41467-023-39810-w
  13. Gokhan Danabasoglu (2004). A comparison of global ocean general circulation model solutions obtained with synchronous and accelerated integration methods. Ocean Modelling 7(3–4), 323–341. DOI: 10.1016/j.ocemod.2003.10.001. https://doi.org/10.1016/j.ocemod.2003.10.001
  14. Stephan J. Lorenz and Gerrit Lohmann (2004). Acceleration technique for Milankovitch type forcing in a coupled atmosphere-ocean circulation model: method and application for the Holocene. Climate Dynamics 23(7–8), 727–743. DOI: 10.1007/s00382-004-0469-y. https://doi.org/10.1007/s00382-004-0469-y
  15. Marco van Hulten, Christoph Heinze, Jörg Schwinger, Jerry Tjiputra (2022). Long-term climate simulation in NorESM: burst-coupling the sediment in the BLOM/iHAMOCC ocean module. arXiv:2211.14087. https://arxiv.org/abs/2211.14087
  16. Lars-Anders Breivik and Magnar Reistad (1994). Assimilation of ERS-1 altimeter wave heights in an operational numerical wave model. Weather and Forecasting 9(3), 440–451. DOI: 10.1175/1520-0434(1994)009<0440:AOAWHI>2.0.CO;2. https://doi.org/10.1175/1520-0434(1994)009%3C0440%3AAOAWHI%3E2.0.CO%3B2
  17. Peter Janssen (2000). ECMWF wave modeling and satellite altimeter wave data. In: David Halpern (ed.), Satellites, Oceanography and Society, Elsevier Oceanography Series, pp. 35–56. DOI: 10.1016/S0422-9894(00)80004-5. ISBN 9780444505019. https://doi.org/10.1016/S0422-9894(00)80004-5
  18. Philip A. Browne, Patricia de Rosnay, Hao Zuo, Andrew Bennett, Andrew Dawson (2019). Weakly coupled ocean-atmosphere data assimilation in the ECMWF NWP system. Remote Sensing 11(3), 234. DOI: 10.3390/rs11030234. https://doi.org/10.3390/rs11030234
  19. Daniel M. Sigman and Edward A. Boyle (2000). Glacial/interglacial variations in atmospheric carbon dioxide. Nature 407, 859–869. DOI: 10.1038/35038000. https://doi.org/10.1038/35038000
  20. Karen E. Kohfeld and Andy Ridgwell (2009). Glacial-interglacial variability in atmospheric CO2. In: Surface Ocean–Lower Atmosphere Processes, Geophysical Monograph Series, pp. 251–286. DOI: 10.1029/2008GM000845. https://doi.org/10.1029/2008GM000845
  21. Hubertus Fischer, Jochen Schmitt, Dieter Lüthi, Thomas F. Stocker, Tobias Tschumi, Payal Parekh, Fortunat Joos, Peter Köhler, Christoph Völker, Rainer Gersonde, Carlo Barbante, Martine Le Floch, Dominique Raynaud, Eric Wolff (2010). The role of Southern Ocean processes in orbital and millennial CO2 variations – A synthesis. Quaternary Science Reviews 29(1–2), 193–205. DOI: 10.1016/j.quascirev.2009.06.007. https://doi.org/10.1016/j.quascirev.2009.06.007
  22. Mathis P. Hain, Daniel M. Sigman, Gerald H. Haug (2010). Carbon dioxide effects of Antarctic stratification, North Atlantic Intermediate Water formation, and subantarctic nutrient drawdown during the last ice age: diagnosis and synthesis in a geochemical box model. Global Biogeochemical Cycles 24(4). DOI: 10.1029/2010GB003790. https://doi.org/10.1029/2010GB003790
  23. Daniel M. Sigman, Mathis P. Hain, Gerald H. Haug (2010). The polar ocean and glacial cycles in atmospheric CO2 concentration. Nature 466, 47–55. DOI: 10.1038/nature09149. https://doi.org/10.1038/nature09149
  24. Andrew J. Watson, Geoffrey K. Vallis, Maxim Nikurashin (2015). Southern Ocean buoyancy forcing of ocean ventilation and glacial atmospheric CO2. Nature Geoscience 8, 861–864. DOI: 10.1038/ngeo2538. https://doi.org/10.1038/ngeo2538
  25. Thomas K. Bauska, Daniel Baggenstos, Edward J. Brook, Alan C. Mix, Shaun A. Marcott, Vasilii V. Petrenko, Hinrich Schaefer, Jeffrey P. Severinghaus, James E. Lee (2016). Carbon isotopes characterize rapid changes in atmospheric carbon dioxide during the last deglaciation. Proceedings of the National Academy of Sciences 113(13), 3465–3470. DOI: 10.1073/pnas.1513868113. https://doi.org/10.1073/pnas.1513868113
  26. Julia Gottschalk, Luke C. Skinner, Jörg Lippold, Hendrik Vogel, Norbert Frank, Samuel L. Jaccard, Claire Waelbroeck (2016). Biological and physical controls in the Southern Ocean on past millennial-scale atmospheric CO2 changes. Nature Communications 7, 11539. DOI: 10.1038/ncomms11539. https://doi.org/10.1038/ncomms11539
  27. T. D. Vollmer, T. Ito, J. Lynch-Stieglitz (2022). Proxy-based preformed phosphate estimates point to increased biological pump efficiency as primary cause of Last Glacial Maximum CO2 drawdown. Paleoceanography and Paleoclimatology 37(11), e2021PA004339. DOI: 10.1029/2021PA004339. https://doi.org/10.1029/2021PA004339
  28. Daniel M. Sigman and Mathis P. Hain (2024). Ocean oxygen, preformed nutrients, and the cause of the lower carbon dioxide concentration in the atmosphere of the Last Glacial Maximum. Paleoceanography and Paleoclimatology 39(1), e2023PA004775. DOI: 10.1029/2023PA004775. https://doi.org/10.1029/2023PA004775
  29. Frederik Schenk, Minna Väliranta, Francesco Muschitiello, Lev Tarasov, Maija Heikkilä, Svante Björck, Jenny Brandefelt, Arne V. Johansson, Jens-Ove Näslund, Barbara Wohlfarth (2018). Warm summers during the Younger Dryas cold reversal. Nature Communications 9, 1634. DOI: 10.1038/s41467-018-04071-5. https://doi.org/10.1038/s41467-018-04071-5
  30. Defang You, Ruediger Stein, Gerrit Lohmann, Ahmadreza Masoum, Kirsten Fahl, Junjie Wu, Rebecca Jackson (2026). Surface freshening in the subpolar North Atlantic sustaining the weakened AMOC during the late Younger Dryas. Science Advances 12(1), eadv6220. DOI: 10.1126/sciadv.adv6220. https://doi.org/10.1126/sciadv.adv6220