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.

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

Strengthening Atlantic Inflow Across the Mid‐Pleistocene Transition

Transient Temperature Asymmetry Between Hemispheres in the Palaeogene Atlantic Ocean

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

Revisiting the mid-Pleistocene transition ocean circulation crisis

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

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

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

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

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

Mid-Pleistocene Climate Transition Triggered by Antarctic Ice Sheet Growth

Arctic and Antarctic forcing of ocean interior warming during the last deglaciation

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

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

Late Holocene thermohaline perturbation of the N-Atlantic Subpolar Gyre linked to exceptional Greenland Ice Sheet melting between 4.4 and 4.0 ka BP

Reconstruction of the late-Holocene changes in the Sub-Arctic Front position at the Reykjanes Ridge, north Atlantic

Asynchronous warming and δ18O evolution of deep Atlantic water masses during the last deglaciation

The North Atlantic Subpolar Ocean Dynamics during the Past 21,000 Years

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

Worldwide consequences of a mid-Holocene cold event in the Nordic Seas

North Atlantic surface ocean warming and salinization in response to middle Eocene greenhouse warming