The origin and bifurcation structure of abrupt millennial-scale climate transitions under steady external solar forcing and in the absence of atmospheric synoptic variability is studied by means of a global coupled model of intermediate complexity. We show that the origin of Dansgaard-Oeschger type oscillations in the model is caused by the weaker northward oceanic heat transport in the Atlantic basin. This is in agreement with previous studies realized with much simpler models, based on highly idealized geometries and simplified physics. The existence of abrupt millennial-scale climate transitions during glacial times can therefore be interpreted as a consequence of the weakening of the negative temperature-advection feedback. This is confirmed through a series of numerical experiments designed to explore the sensitivity of the bifurcation structure of the Atlantic meridional overturning circulation to increased atmospheric CO2 levels under glacial boundary conditions. Contrasting with the cold, stadial, phases of millennial oscillations, we also show the emergence of strong interdecadal variability in the North Atlantic sector during warm interstadials. The instability driving these interdecadal-interstadial oscillations is shown to be identical to that found in ocean-only models forced by fixed surface buoyancy fluxes, that is, a large-scale baroclinic instability developing in the vicinity of the western boundary current in the North Atlantic. Comparisons with modern observations further suggest a physical mechanism similar to that driving the 30–40 years time scale associated with the Atlantic multidecadal oscillation.
