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.

Early to Late Pliocene Climate Change at Mid-Latitudinal North Atlantic Site U1313

Marine and Terrestrial Biomarker Records for the Last 35,000 Years at ODP Site 658C off NW Africa

Unraveling the mechanisms and implications of a stronger mid-Pliocene AMOC in PlioMIP2

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

Oxygenated Deep Waters Fed Early Atlantic Overturning Circulation Upon Antarctic Glaciation

Subsurface heat channel drove sea surface warming in the high‐latitude North Atlantic during the Mid‐Pleistocene Transition

Global Deepwater Circulation Between 2.4 and 1.7ma and Its Connection to the Onset of Northern Hemisphere Glaciation

Dynamics of the subpolar gyre and transition zone of the North Atlantic during the last glacial cycle

Records of melt water pulses and climate events from the northern South China Sea since late Last Glacial Maximum:Implication of terrigenous input

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

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

Transient Temperature Asymmetry Between Hemispheres in the Palaeogene Atlantic Ocean

Toward a high-resolution stable isotope stratigraphy of the last 3.4 million years: Sites 658 and 659 off Northwest Africa. In Ruddiman, W., Sarnthein, M. et al.

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

Late Miocene decoupling of oceanic warmth and atmospheric carbon dioxide forcing

Impacts of an active Pacific Meridional Overturning Circulation on the Pliocene climate and hydrological cycle

Re-organization of Pacific overturning circulation across the Miocene Climate Optimum

Neodymium Isotopic Evidence for Linked Changes in Southeast Atlantic and Southwest Pacific Circulation over the Last 200 Kyr

North Atlantic Temperature Change Across the Eocene‐Oligocene Transition From Clumped Isotopes

Climate- and gateway-driven cooling of Late Eocene to earliest Oligocene sea surface temperatures in the North Sea Basin

Multi‐Proxy Evidence for Atlantic Meridional Overturning Circulation (AMOC) Weakening During Deglaciations of the Past 150,000 Years