Abstract:Climate during the last glacial period was marked by abrupt instability on millennial timescales that included large swings of temperature in and around Greenland (Daansgard–Oeschger events) and smaller, more gradual changes in Antarctica (AIM events). Less is known about the existence and nature of similar variability during older glacial periods, especially during the early Pleistocene when glacial cycles were dominantly occurring at 41 kyr intervals compared to the much longer and deeper glaciations of the more recent period. Here, we report a continuous millennially resolved record of stable isotopes of planktic and benthic foraminifera at IODP Site U1385 (the "Shackleton Site") from the southwestern Iberian margin for the last 1.5 million years, which includes the Middle Pleistocene Transition (MPT). Our results demonstrate that millennial climate variability (MCV) was a persistent feature of glacial climate, both before and after the MPT. Prior to 1.2 Ma in the early Pleistocene, the amplitude of MCV was modulated by the 41 kyr obliquity cycle and increased when axial tilt dropped below 23.5 ∘ and benthic δ 18 O exceeded ∼3.8 ‰ (corrected to Uvigerina ), indicating a threshold response to orbital forcing. Afterwards, MCV became focused mainly on the transitions into and out of glacial states (i.e. inceptions and terminations) and during times of intermediate ice volume. After 1.2 Ma, obliquity continued to play a role in modulating the amplitude of MCV, especially during times of glacial inceptions, which are always associated with declining obliquity. A non-linear role for obliquity is also indicated by the appearance of multiples (82, 123 kyr) and combination tones (28 kyr) of the 41 kyr cycle. Near the end of the MPT ( ∼0.65 Ma), obliquity modulation of MCV amplitude wanes as quasi-periodic 100 kyr and precession power increase, coinciding with the growth of oversized ice sheets on North America and the appearance of Heinrich layers in North Atlantic sediments. Whereas the planktic δ 18 O of Site U1385 shows a strong resemblance to Greenland temperature and atmospheric methane (i.e. Northern Hemisphere climate), millennial changes in benthic δ 18 O closely follow the temperature history of Antarctica for the past 800 kyr. The phasing of millennial planktic and benthic δ 18 O variation is similar to that observed for MIS 3 throughout much of the record, which has been suggested to mimic the signature of the bipolar seesaw – i.e. an interhemispheric asymmetry between the timing of cooling in Antarctica and warming in Greenland. The Iberian margin isotopic record suggests that bipolar asymmetry was a robust feature of interhemispheric glacial climate variations for at least the past 1.5 Ma despite changing glacial boundary conditions. A strong correlation exists between millennial increases in planktic δ 18 O (cooling) and decreases in benthic δ 13 C , indicating that millennial variations in North Atlantic surface temperature are mirrored by changes in deep-water circulation and remineralization of carbon in the abyssal ocean. We find strong evidence that climate variability on millennial and orbital scales is coupled across different timescales and interacts in both directions, which may be important for linking internal climate dynamics and external astronomical forcing.

Centennial-scale variations in the carbon cycle enhanced by high obliquity

Synchronous Centennial Abrupt Events in the Ocean and Atmosphere During the Last Deglaciation

Centennial-scale changes in the global carbon cycle during the last deglaciation

Southern Ocean drives multidecadal atmospheric CO2 rise during Heinrich Stadials

Rapid shifts in circulation and biogeochemistry of the Southern Ocean during deglacial carbon cycle events.

Exploring Cyclic Changes of the Ocean Carbon Reservoir

Greenhouse gases in the Earth system: a palaeoclimate perspective

CO2 Storage and Release in the Deep Southern Ocean on Millennial to Centennial Timescales

Carbon isotope budget indicates biological disequilibrium dominated ocean carbon storage at the Last Glacial Maximum

Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary

Astronomical Constraints on Global Carbon-Cycle Perturbation During Oceanic Anoxic Event 2 (OAE2)

Abrupt climate and weather changes across time scales

Carbon Cycle and Circulation Change in the North Pacific Ocean at the Initiation of Northern Hemisphere Glaciation Constrained by Boron‐Based Proxies in Diatoms

Changes in Deep Pacific Circulation and Carbon Storage During the Pliocene-Pleistocene Transition

Modelling Coupled Oscillations of Volcanic CO2 Emissions and Glacial Cycles

A 1.5-million-year record of orbital and millennial climate variability in the North Atlantic

Orbital‐Scale Global Ocean Sea Surface Temperatures Coupling With Cryosphere‐Carbon Cycle Changes Over the Past 4 Million Years

Covariation of Deep Antarctic Pacific Oxygenation and Atmospheric CO2 During the Last 770 Kyr

The nature of deep overturning and reconfigurations of the silicon cycle across the last deglaciation

The grandest of them all: the Lomagundi–Jatuli Event and Earth's oxygenation

Direct astronomical influence on abrupt climate variability