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.

The Middle Pleistocene Transition by frequency locking and slow ramping of internal period

The bifurcation structure and noise induced transitions in the Pleistocene glacial cycles

Modelling Coupled Oscillations of Volcanic CO2 Emissions and Glacial Cycles

Glacial abrupt climate change as a multi-scale phenomenon resulting from monostable excitable dynamics

Interpretation and Preliminary Simulation of the 40 Ka Periodicity of the Quaternary Glaciation

Bayesian Data Analysis for Revealing Causes of the Middle Pleistocene Transition

Distinguishing Between Insolation‐Driven and Phase‐Locked 100‐Kyr Ice Age Scenarios Using Example Models

Direct astronomical influence on abrupt climate variability

Revisiting the mid-Pleistocene transition ocean circulation crisis

Oscillators and relaxation phenomena in Pleistocene climate theory

Nonlinear resonance in an overdamped Duffing oscillator as a model of paleoclimate oscillations

Obliquity and precession as pacemakers of Pleistocene deglaciations

Strengthening Atlantic Inflow Across the Mid‐Pleistocene Transition

Abrupt climate and weather changes across time scales

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

A Punctuated Equilibrium Analysis of the Climate Evolution of Cenozoic exhibits a Hierarchy of Abrupt Transitions

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

Coupling between Grand cycles and Events in Earth’s climate during the past 115 million years

Deep-sea hiatus record reveals orbital pacing by 2.4 Myr eccentricity grand cycles

A carbon cycle coupled climate model of Neoproterozoic glaciation: Explicit carbon cycle with stochastic perturbations