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.

Persistent orbital influence on millennial climate variability through the Pleistocene

Millennial climate variability and organic matter accumulation under icehouse conditions

Direct astronomical influence on abrupt climate variability

Millennial-scale East Asian monsoon variability of the last glacial deduced from annually laminated sediments from Lake Sihailongwan, N.E. China

Obliquity and precession as pacemakers of Pleistocene deglaciations

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

Orbital- and millennial-scale Asian winter monsoon variability across the Pliocene–Pleistocene glacial intensification

Millennial-Scale Monsoon Variability Modulated by Low-Latitude Insolation During the Last Glaciation

Mineral dust coupled with climate‐carbon cycle on orbital timescales over the past 4 Ma

Solar Imprints on Asian Inland Moisture Fluctuations over the Last Millennium

Dominant eccentricity cycles in paleoenvironmental variabilities recorded by pelagic sediments in the western Pacific during 15–11 Ma

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

Astronomical influence on Miocene continental environments in Tibet, east-central Asia

Co-evolution of orbital-to-millennial East Asian monsoon and westerlies variability over the past 150 ka

Millennial resolution late Miocene northern China precipitation record spanning astronomical analogue interval to the future

A review of orbital-scale monsoon variability and dynamics in East Asia during the Quaternary

Orbital and Millennial‐Scale Cycles Through the Hirnantian (Late Ordovician) in Southern China

North Atlantic influence on the glacial amplitude of East Asian millennial-scale monsoon variability

High-sedimentation-rate loess records: A new window into understanding orbital- and millennial-scale monsoon variability

Multi-centennial climate variability: Observational, theoretical and modeling studies

Orbital Modulation of an Intensified Hydrological Cycle During the Paleocene-Eocene Thermal Maximum