800,000 years of abrupt climate variability

Millennial-scale climate variability of the Asian summer monsoon over the last 690,000 years: insights from cave records

The Asian summer monsoon (ASM) is an important component of the global climate system. Cave oxygen isotope (δ18O) records from the region have well characterized the ASM millennial-scale climate variability (MCV) over the last 640 ka, with especially detailed insights for the most recent 60 ka, but little is known about ASM variability beyond the U-Th dating limit of ∼640 ka. Furthermore, questions remain regarding the climatic significance of millennial-scale ASM variability recorded among various climate archives, particularly in the context of the "orbital-scale paradox". Here, we present high-resolution and U-Pb dated cave δ18O records from two Chinese caves spanning 690-600 ka BP (before present, where present = 1950 CE). These records reveal coupling between millennial-scale ASM weakening, North Atlantic cooling, and Antarctic warming, essentially mirroring the pattern observed during the last 640 ka, despite a potential change in the Atlantic meridional overturning circulation (AMOC) caused by a switch in the freshwater route to the North Atlantic at ∼640 ka BP. Comparisons of MCV amplitudes among different ASM proxies show remarkable disparities, suggesting that each proxy reflects different aspects of the ASM. In our records, declines in summer insolation repeatedly triggered millennial-scale weak ASM events near the mid-precession band, associated with AMOC weakening, rather than only at interglacial terminations. Additionally, our analysis highlights the critical roles of atmospheric CO2 and global ice volume conditions in shaping the ASM variability during ∼688-635 ka BP.

Obliquity and precession forcing of the amplitude of millennial-scale East Asian monsoon variability during the late Miocene

Asian monsoon systems influence billions of people and understanding past monsoon variability and dynamics is instructive in predicting future trajectories. Recent studies have demonstrated that insolation, plus interrelated changes in ice-sheet extent and volume, particularly in the Northern Hemisphere (NH), have controlled the magnitude of millennial-scale East Asian summer monsoon (EASM) variability during the Quaternary. However, it is unclear how the EASM varied at both orbital and millennial timescales, and whether orbital-scale variations impacted millennial-scale variations, during intervals lacking large permanent NH ice sheets at times when CO2 levels were close to the present-day value of approximately 400 ppm. Here, we present high-resolution (∼1-kyr) dry-wet variation records from late Miocene eolian sediments in the Jianzha Basin of the northeastern Tibetan Plateau, spanning an interval when CO2 levels were persistently close to the modern level and the NH lacked large permanent ice sheets. Our results reveal orbital-scale forcing interwoven with millennial cycles. In contrast to dominant precession and eccentricity forcing of the EASM at orbital scales, the amplitude variations of millennial-scale EASM cycles exhibit strong obliquity (and its modulating cycle) forcing and weak precession forcing. This pattern is different from the pattern observed in either the Quaternary or early Miocene, which we attribute to the effects of different boundary conditions.

790,000 years of millennial-scale Cape Horn Current variability and interhemispheric linkages

Millennial-scale variations in the strength and position of the Antarctic Circumpolar Current exert considerable influence on the global meridional overturning circulation and the ocean carbon cycle. The mechanistic understanding of these variations is still incomplete, partly due to the scarcity of sediment records covering multiple glacial-interglacial cycles with millennial-scale resolution. Here, we present high-resolution current strength and sea surface temperature records covering the past 790,000 years from the Cape Horn Current as part of the subantarctic Antarctic Circumpolar Current system, flowing along the Chilean margin. Both temperature and current velocity data document persistent millennial-scale climate variability throughout the last eight glacial periods with stronger current flow and warmer sea surface temperatures coinciding with Antarctic warm intervals. These Southern Hemisphere changes are linked to North Atlantic millennial-scale climate fluctuations, plausibly involving changes in the Atlantic thermohaline circulation. The variations in the Antarctic Circumpolar Current system are associated with atmospheric CO2 changes, suggesting a mechanistic link through the Southern Ocean carbon cycle.

Distinct roles for precession, obliquity, and eccentricity in Pleistocene 100-kyr glacial cycles

Identifying the specific roles of precession, obliquity, and eccentricity in glacial-interglacial transitions is hindered by imprecise age control. We circumvent this problem by focusing on the morphology of deglaciation and inception, which we show depends strongly on the relative phasing of precession versus obliquity. We demonstrate that although both parameters are important, precession has more influence on deglacial onset, whereas obliquity is more important for the attainment of peak interglacial conditions and glacial inception. We find that the set of precession peaks (minima) responsible for terminations since 0.9 million years ago is a subset of those peaks that begin (i.e., the precession parameter starts decreasing) while obliquity is increasing. Specifically, termination occurs with the first of these candidate peaks to occur after each eccentricity minimum. Thus, the gross morphology of 100-thousand-year (100-kyr) glacial cycles appears largely deterministic.

Asynchronicity of deglacial permafrost thawing controlled by millennial-scale climate variability

Permafrost is a potentially important source of deglacial carbon release alongside deep-sea carbon outgassing. However, limited proxies have restricted our understanding in circumarctic regions and the last deglaciation. Tibetan Plateau (TP), the Earth's largest low-latitude and alpine permafrost region, remains underexplored. Using speleothem growth phases, we reconstruct TP permafrost thawing history over the last 500,000 years, standardizing chronology to investigate Northern Hemisphere permafrost thawing patterns. We find that, unlike circumarctic permafrost, TP permafrost generally initiates thawing at the onset of deglaciations, coinciding with Weak Monsoon Intervals and sluggish Atlantic Meridional Overturning Circulation (AMOC) during Terminal Stadials. Modeling elaborates that the associated Asian monsoon weakening induces anomalous TP warming through local cloud-precipitation-soil moisture feedback. This, combined with high-latitude cooling, results in asynchronous boreal permafrost thawing. During the last deglaciation, however, anomalous AMOC variability delayed TP and advanced circumarctic permafrost thawing. Our results indicate that permafrost carbon release, influenced by millennial-scale AMOC variability, may have been a non-trivial contributor to deglacial CO2 rise.

Nonlinear comparative analysis of Greenland and Antarctica ice cores data

We analyze the temperature time series of the EPICA Dome C ice cores in Antarctica and of the Greenland project, Summit, with durations of 800 000 and 248 000 years, respectively, with a recent mathematical tool defined through the Fourier phases of the series, known as the J-index. This data driven index can differentiate between purely random dynamics and dynamics with a deterministic component. It is sensitive to nonlinear components and robust to the presence of noise. Our J-index data analysis shows that both Greenland and Antarctica climatic fluctuations possess deterministic traits and suggests the presence of an underlying nonlinear dynamics. Furthermore, in both regions, it reveals the simultaneous occurrence of an important global event known as the "Pelukian transgression." For Antarctica, it also detects the marine isotopic stage 11. Additionally, our calculation of the time series Hurst exponents and our detrended fluctuation analysis show the presence of long-range persistent correlations for Antarctica and anti-persistent correlations for Greenland. For the latter case, our fractal dimension determinations are indicative of a more complex climatic dynamics in Greenland with respect to Antarctica. Our results are encouraging for further development of climate variability deterministic models for these regions.

Dansgaard-Oeschger cycles of the penultimate and last glacial period recorded in stalagmites from Türkiye

The last glacial period is characterized by abrupt climate oscillations, also known as Dansgaard-Oeschger (D-O) cycles. However, D-O cycles remain poorly documented in climate proxy records covering the penultimate glacial period. Here we present highly resolved and precisely dated speleothem time series from Sofular Cave in northern Türkiye to provide clear evidence for D-O cycles during Marine Isotope Stage (MIS) 6 as well as MIS 2-4. D-O cycles are most clearly expressed in the Sofular carbon isotope time series, which correlate inversely with regional sea surface temperature (SST) records from the Black Sea. The pacing of D-O cycles is almost twice as long during MIS 6 compared to MIS 2-4, and could be related to a weaker Atlantic Meridional Overturning Circulation (AMOC) and a different mean climate during MIS 6 compared to MIS 2-4, leading most likely to a higher threshold for the occurrence of D-O cycles.

A reliable benchmark of the last 640,000 years millennial climate variability

How often have past climates undergone abrupt transitions? While our understanding of millennial variability during the past 130,000 years is well established, with precise dates available, such information on previous climate cycles is limited. To address this question, we identified 196 abrupt transitions in the δ18O record of the well-dated Chinese composite speleothem for the last 640,000 years. These results correspond to abrupt changes in the strength of the East Asian Monsoon, which align with the Greenland stadials and interstadials observed in the North Atlantic region during the last 130,000 years before present. These precise dates of past abrupt climate changes constitute a reliable and necessary benchmark for Earth System models used to study future climate scenarios.

Reversed Holocene temperature-moisture relationship in the Horn of Africa

Anthropogenic climate change is predicted to severely impact the global hydrological cycle1, particularly in tropical regions where agriculture-based economies depend on monsoon rainfall2. In the Horn of Africa, more frequent drought conditions in recent decades3,4 contrast with climate models projecting precipitation to increase with rising temperature5. Here we use organic geochemical climate-proxy data from the sediment record of Lake Chala (Kenya and Tanzania) to probe the stability of the link between hydroclimate and temperature over approximately the past 75,000 years, hence encompassing a sufficiently wide range of temperatures to test the 'dry gets drier, wet gets wetter' paradigm6 of anthropogenic climate change in the time domain. We show that the positive relationship between effective moisture and temperature in easternmost Africa during the cooler last glacial period shifted to negative around the onset of the Holocene 11,700 years ago, when the atmospheric carbon dioxide concentration exceeded 250 parts per million and mean annual temperature approached modern-day values. Thus, at that time, the budget between monsoonal precipitation and continental evaporation7 crossed a tipping point such that the positive influence of temperature on evaporation became greater than its positive influence on precipitation. Our results imply that under continued anthropogenic warming, the Horn of Africa will probably experience further drying, and they highlight the need for improved simulation of both dynamic and thermodynamic processes in the tropical hydrological cycle.

The PaleoJump database for abrupt transitions in past climates

Tipping points (TPs) in Earth's climate system have been the subject of increasing interest and concern in recent years, given the risk that anthropogenic forcing could cause abrupt, potentially irreversible, climate transitions. Paleoclimate records are essential for identifying past TPs and for gaining a thorough understanding of the underlying nonlinearities and bifurcation mechanisms. However, the quality, resolution, and reliability of these records can vary, making it important to carefully select the ones that provide the most accurate representation of past climates. Moreover, as paleoclimate time series vary in their origin, time spans, and periodicities, an objective, automated methodology is crucial for identifying and comparing TPs. To address these challenges, we introduce the open-source PaleoJump database, which contains a collection of carefully selected, high-resolution records originating in ice cores, marine sediments, speleothems, terrestrial records, and lake sediments. These records describe climate variability on centennial, millennial and longer time scales and cover all the continents and ocean basins. We provide an overview of their spatial distribution and discuss the gaps in coverage. Our statistical methodology includes an augmented Kolmogorov-Smirnov test and Recurrence Quantification Analysis; it is applied here, for illustration purposes, to selected records in which abrupt transitions are automatically detected and the presence of potential tipping elements is investigated. These transitions are shown in the PaleoJump database along with other essential information about the records, including location, temporal scale and resolution, as well as temporal plots. This open-source database represents, therefore, a valuable resource for researchers investigating TPs in past climates.

Decadal-centennial-scale solar-linked climate variations and millennial-scale internal oscillations during the Early Cretaceous

Understanding climate variability and stability under extremely warm 'greenhouse' conditions in the past is essential for future climate predictions. However, information on millennial-scale (and shorter) climate variability during such periods is scarce, owing to a lack of suitable high-resolution, deep-time archives. Here we present a continuous record of decadal- to orbital-scale continental climate variability from annually laminated lacustrine deposits formed during the late Early Cretaceous (123-120 Ma: late Barremian-early Aptian) in southeastern Mongolia. Inter-annual changes in lake algal productivity for a 1091-year interval reveal a pronounced solar influence on decadal- to centennial-scale climatic variations (including the ~ 11-year Schwabe cycle). Decadally-resolved Ca/Ti ratios (proxy for evaporation/precipitation changes) for a ~ 355-kyr long interval further indicate millennial-scale (~ 1000-2000-yr) extreme drought events in inner-continental areas of mid-latitude palaeo-Asia during the Cretaceous. Millennial-scale oscillations in Ca/Ti ratio show distinct amplitude modulation (AM) induced by the precession, obliquity and short eccentricity cycles. Similar millennial-scale AM by Milankovitch cycle band was also previously observed in the abrupt climatic oscillations (known as Dansgaard-Oeschger events) in the 'intermediate glacial' state of the late Pleistocene, and in their potential analogues in the Jurassic 'greenhouse'. Our findings indicate that external solar activity forcing was effective on decadal-centennial timescales, whilst the millennial-scale variations were likely amplified by internal process such as changes in deep-water formation strength, even during the Cretaceous 'greenhouse' period.

Milankovitch theory and monsoon

The widely accepted “Milankovitch theory” explains insolation-induced waxing and waning of the ice sheets and their effect on the global climate on orbital timescales. In the past half century, however, the theory has often come under scrutiny, especially regarding its “100-ka problem.” Another drawback, but the one that has received less attention, is the “monsoon problem,” which pertains to the exclusion of monsoon dynamics in classic Milankovitch theory even though the monsoon prevails over the vast low-latitude (∼30° N to ∼30° S) region that covers half of the Earth’s surface and receives the bulk of solar radiation. In this review, we discuss the major issues with the current form of Milankovitch theory and the progress made at the research forefront. We suggest shifting the emphasis from the ultimate outcomes of the ice volume to the causal relationship between changes in northern high-latitude insolation and ice age termination events (or ice sheet melting rate) to help reconcile the classic “100-ka problem.” We discuss the discrepancies associated with the characterization of monsoon dynamics, particularly the so-called “sea-land precession-phase paradox” and the “Chinese 100-ka problem.” We suggest that many of these discrepancies are superficial and can be resolved by applying a holistic “monsoon system science” approach. Finally, we propose blending the conventional Kutzbach orbital monsoon hypothesis, which calls for summer insolation forcing of monsoons, with Milankovitch theory to formulate a combined “Milankovitch-Kutzbach hypothesis” that can potentially explain the dual nature of orbital hydrodynamics of the ice sheet and monsoon systems, as well as their interplays and respective relationships with the northern high-latitude insolation and inter-tropical insolation differential.

•

Orbital-scale climate variations of Earth are dictated by ice sheet and monsoon

•

Views of “monsoon system science” reinforce the Kutzbach monsoon hypothesis

•

A unified Milankovitch-Kutzbach hypothesis better explains the orbital dual nature

Changes in North Atlantic deep-water oxygenation across the Middle Pleistocene Transition

The oxygen concentrations of oceanic deep-water and atmospheric carbon dioxide (pCO₂) are intrinsically linked through organic carbon remineralization and storage as dissolved inorganic carbon in the deep sea. We present a high-resolution reconstruction of relative changes in oxygen concentration in the deep North Atlantic for the past 1.5 million years using the carbon isotope gradient between epifaunal and infaunal benthic foraminifera species as a proxy for paleo-oxygen. We report a significant (>40 micromole per kilogram) reduction in glacial Atlantic deep-water oxygenation at ~960 thousand to 900 thousand years ago that coincided with increased continental ice volume and a major change in ocean thermohaline circulation. Paleo-oxygen results support a scenario of decreasing deep-water oxygen concentrations, increased respired carbon storage, and a reduction in glacial pCO₂ across the Middle Pleistocene Transition.

Rapid northern hemisphere ice sheet melting during the penultimate deglaciation

The rate and consequences of future high latitude ice sheet retreat remain a major concern given ongoing anthropogenic warming. Here, new precisely dated stalagmite data from NW Iberia provide the first direct, high-resolution records of periods of rapid melting of Northern Hemisphere ice sheets during the penultimate deglaciation. These records reveal the penultimate deglaciation initiated with rapid century-scale meltwater pulses which subsequently trigger abrupt coolings of air temperature in NW Iberia consistent with freshwater-induced AMOC slowdowns. The first of these AMOC slowdowns, 600-year duration, was shorter than Heinrich 1 of the last deglaciation. Although similar insolation forcing initiated the last two deglaciations, the more rapid and sustained rate of freshening in the eastern North Atlantic penultimate deglaciation likely reflects a larger volume of ice stored in the marine-based Eurasian Ice sheet during the penultimate glacial in contrast to the land-based ice sheet on North America as during the last glacial.

Stalagmites from NW Iberia record the rapid demise of large ice sheets during the penultimate deglaciation, and reveal decadal-scale feedbacks between warming and ice melting.

700,000 years of tropical Andean glaciation

Our understanding of the climatic teleconnections that drove ice-age cycles has been limited by a paucity of well-dated tropical records of glaciation that span several glacial–interglacial intervals. Glacial deposits offer discrete snapshots of glacier extent but cannot provide the continuous records required for detailed interhemispheric comparisons. By contrast, lakes located within glaciated catchments can provide continuous archives of upstream glacial activity, but few such records extend beyond the last glacial cycle. Here a piston core from Lake Junín in the uppermost Amazon basin provides the first, to our knowledge, continuous, independently dated archive of tropical glaciation spanning 700,000 years. We find that tropical glaciers tracked changes in global ice volume and followed a clear approximately 100,000-year periodicity. An enhancement in the extent of tropical Andean glaciers relative to global ice volume occurred between 200,000 and 400,000 years ago, during sustained intervals of regionally elevated hydrologic balance that modified the regular approximately 23,000-year pacing of monsoon-driven precipitation. Millennial-scale variations in the extent of tropical Andean glaciers during the last glacial cycle were driven by variations in regional monsoon strength that were linked to temperature perturbations in Greenland ice cores¹; these interhemispheric connections may have existed during previous glacial cycles.

Analysis of a continuous and independently dated record of glaciation in the tropical Andes spanning 700,000 years shows that Andean glaciation follows patterns of global ice volume change, with a periodicity of approximately 100,000 years.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

An In-Depth Analysis of Physical Blue and Green Water Scarcity in Agriculture in Terms of Causes and Events and Perceived Amenability to Economic Interpretation

An analytical review of physical blue and green water scarcity in terms of agricultural use, and its amenability to economic interpretation, is presented, employing more than 600 references. The main definitions and classifications involved and information about reserves and resources are critically analyzed, blue and green water scarcity are examined along with their interchange, while their causal connection with climate in general is analyzed along with the particular instances of Europe, Africa, Asia and the WANA region. The role of teleconnections and evaporation/moisture import-export is examined as forms of action at a distance. The human intervention scarcity driver is examined extensively in terms of land use land cover change (LULCC), as well as population increase. The discussion deals with following critical problems: green and blue water availability, inadequate accessibility, blue water loss, unevenly distributed precipitation, climate uncertainty and country level over global level precedence. The conclusion singles out, among others, problems emerging from the inter-relationship of physical variables and the difficulty to translate them into economic instrumental variables, as well as the lack of imbedding uncertainty in the underlying physical theory due to the fact that country level measurements are not methodically assumed to be the basic building block of regional and global water scarcity.

Millennial scale feedbacks determine the shape and rapidity of glacial termination

Within the Late Pleistocene, terminations describe the major transitions marking the end of glacial cycles. While it is established that abrupt shifts in the ocean/atmosphere system are a ubiquitous component of deglaciation, significant uncertainties remain concerning their specific role and the likelihood that terminations may be interrupted by large-amplitude abrupt oscillations. In this perspective we address these uncertainties in the light of recent developments in the understanding of glacial terminations as the ultimate interaction between millennial and orbital timescale variability. Innovations in numerical climate simulation and new geologic records allow us to highlight new avenues of research and identify key remaining uncertainties such as sea-level variability.

Speleothem record of mild and wet mid-Pleistocene climate in northeast Greenland

The five interglacials before the Mid-Brunhes Event (MBE) [c.430 thousand years (ka) ago] are generally considered to be globally cooler than those post-MBE. Inhomogeneities exist regionally, however, which suggest that the Arctic was warmer than present during Marine Isotope Stage (MIS) 15a. Using the first speleothem record for the High Arctic, we investigate the climatic response of northeast Greenland between c.588 and c.549 ka ago. Our results indicate an enhanced warmth of at least +3.5°C relative to the present, leading to permafrost thaw and increased precipitation. We find that δ¹⁸O of precipitation was at least 3‰ higher than today and recognize two local cooling events (c.571 and c.594 ka ago) thought to be caused by freshwater forcing. Our results are important for improving understanding of the regional climatic response leading up to the MBE and specifically provide insights into the climatic response of a warmer Arctic.

A late Pleistocene dataset of Agulhas Current variability

The interocean transfer of thermocline water between the Indian and the Atlantic Oceans known as 'Agulhas leakage' is of global significance as it influences the Atlantic Meridional Overturning Circulation (AMOC) on different time scales. Variability in the Agulhas Current regime is key in shaping hydroclimate on the adjacent coastal areas of the African continent today as well as during past climates. However, the lack of long, continuous records from the proximal Agulhas Current region dating beyond the last glacial cycle prevents elucidation of its role in regional and wider global climate changes. This is the first continuous record of hydrographic variability (SST; δ18Osw) from the Agulhas Current core region spanning the past 270,000 years. The data set is analytical sound and provides a solid age model. As such, it can be used by paleoclimate scientists, archaeologists, and climate modelers to evaluate, for example, linkages between the Agulhas Current system and AMOC dynamics, as well as connections between ocean heat transport and Southern African climate change in the past and its impact on human evolution.

Regional impacts of climate change and its relevance to human evolution

The traditional concept of long and gradual, glacial-interglacial climate changes during the Quaternary has been challenged since the 1980s. High temporal resolution analysis of marine, terrestrial and ice geological archives has identified rapid, millennial- to centennial-scale, and large-amplitude climatic cycles throughout the last few million years. These changes were global but have had contrasting regional impacts on the terrestrial and marine ecosystems, with in some cases strong changes in the high latitudes of both hemispheres but muted changes elsewhere. Such a regionalization has produced environmental barriers and corridors that have probably triggered niche contractions/expansions of hominin populations living in Eurasia and Africa. This article reviews the long- and short-timescale ecosystem changes that have punctuated the last few million years, paying particular attention to the environments of the last 650,000 years, which have witnessed key events in the evolution of our lineage in Africa and Eurasia. This review highlights, for the first time, a contemporaneity between the split between Denisovan and Neanderthals, at ~650-400 ka, and the strong Eurasian ice-sheet expansion down to the Black Sea. This ice expansion could form an ice barrier between Europe and Asia that may have triggered the genetic drift between these two populations.

Magnesium in subaqueous speleothems as a potential palaeotemperature proxy

Few palaeoclimate archives beyond the polar regions preserve continuous and datable palaeotemperature proxy time series over multiple glacial-interglacial cycles. This hampers efforts to develop a more coherent picture of global patterns of past temperatures. Here we show that Mg concentrations in a subaqueous speleothem from an Italian cave track regional sea-surface temperatures over the last 350,000 years. The Mg shows higher values during warm climate intervals and converse patterns during cold climate stages. In contrast to previous studies, this implicates temperature, not rainfall, as the principal driver of Mg variability. The depositional setting of the speleothem gives rise to Mg partition coefficients that are more temperature dependent than other calcites, enabling the effect of temperature change on Mg partitioning to greatly exceed the effects of changes in source-water Mg/Ca. Subaqueous speleothems from similar deep-cave environments should be capable of providing palaeotemperature information over multiple glacial-interglacial cycles.

Few palaeoclimate archives beyond the polar regions preserve continuous and datable paleotemperature proxy time series over multiple glacial-interglacial cycles. Here, the authors show that Mg concentrations in a subaqueous speleothem from an Italian cave track regional sea-surface temperatures over the last 350,000 years.

Millennial climate oscillations controlled the structure and evolution of Termination II

The controls that affect the structure and timing of terminations are still poorly understood. We studied a tufa deposit from the Iberian Peninsula that covers Termination II (T-II) and whose chronology was synchronized to speleothem records. We used the same chronology to synchronize ocean sediments from the North Atlantic to correlate major climate events in a common timescale. We identify two stages within T-II. The first stage started with the increase of boreal summer integrated solar insolation, and during this stage three millennial climate oscillations were recorded. These oscillations resulted from complex ocean-atmosphere interactions in the Nordic seas, caused by the progressive decay of Northern Hemisphere ice-sheets. The second stage commenced after a glacial outburst that caused the collapse of the Thermohaline Circulation, a massive Heinrich event, and the onset of the Bipolar Seesaw Mechanism (BSM) that eventually permitted the completion of T-II. The pace of the millennial oscillations during the first stage of T-II controlled the onset of the second stage, when the termination became a non-reversible and global phenomenon that accelerated the deglaciation. During the last the two terminations, the BSM was triggered by different detailed climate interactions, which suggests the occurrence of different modes of terminations.

Abrupt CO2 release to the atmosphere under glacial and early interglacial climate conditions

Pulse-like carbon dioxide release to the atmosphere on centennial time scales has only been identified for the most recent glacial and deglacial periods and is thought to be absent during warmer climate conditions. Here, we present a high-resolution carbon dioxide record from 330,000 to 450,000 years before present, revealing pronounced carbon dioxide jumps (CDJ) under cold and warm climate conditions. CDJ come in two varieties that we attribute to invigoration or weakening of the Atlantic meridional overturning circulation (AMOC) and associated northward and southward shifts of the intertropical convergence zone, respectively. We find that CDJ are pervasive features of the carbon cycle that can occur during interglacial climate conditions if land ice masses are sufficiently extended to be able to disturb the AMOC by freshwater input.

Smoothed millennial-scale palaeoclimatic reference data as unconventional comparison targets: Application to European loess records

Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of δ18O data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last ~15 ka and between ~50-30 ka. Between ~30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.

Persistent influence of obliquity on ice age terminations since the Middle Pleistocene transition

Radiometric dating of glacial terminations over the past 640,000 years suggests pacing by Earth's climatic precession, with each glacial-interglacial period spanning four or five cycles of ~20,000 years. However, the lack of firm age estimates for older Pleistocene terminations confounds attempts to test the persistence of precession forcing. We combine an Italian speleothem record anchored by a uranium-lead chronology with North Atlantic ocean data to show that the first two deglaciations of the so-called 100,000-year world are separated by two obliquity cycles, with each termination starting at the same high phase of obliquity, but at opposing phases of precession. An assessment of 11 radiometrically dated terminations spanning the past million years suggests that obliquity exerted a persistent influence on not only their initiation but also their duration.

Millennial-Scale Climate Variability and Dinoflagellate-Cyst-Based Seasonality Changes Over the Last ~150 kyrs at "Shackleton Site" U1385

During the last glacial period, climate conditions in the North Atlantic region were determined by the alternation of relatively warm interstadials and relatively cool stadials, with superimposed rapid warming (Dansgaard-Oeschger) and cooling (Heinrich) events. So far little is known about the impact of these rapid climate shifts on the seasonal variations in sea surface temperature (SST) within the North Atlantic region. Here, we present a high-resolution seasonal SST record for the past 152 kyrs derived from Integrated Ocean Drilling Program "Shackleton" Site U1385, offshore Portugal. Assemblage counts of dinoflagellates cysts (dinocysts) in combination with a modern analog technique (MAT), and regression analyses were used for the reconstructions. We compare our records with previously published SST records from the same location obtained from the application of MAT on planktonic foraminifera. Our dinocyst-based reconstructions confirm the impression of the Greenland stadials and interstadials offshore the Portuguese margin and indicate increased seasonal contrast of temperature during the cold periods of the glacial cycle (average 9.0 °C, maximum 12.2 °C) with respect to present day (5.1 °C), due to strong winter cooling by up to 8.3 °C. Our seasonal temperature reconstructions are in line with previously published data, which showed increased seasonality due to strong winter cooling during the Younger Dryas and the Last Glacial Maximum over the European continent and North Atlantic region. In addition, we show that over longer time scales, increased seasonal contrasts of temperature remained characteristic of the colder phases of the glacial cycle.

Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison

Earlier large-scale Greenland ice sheet sea-level projections (e.g., those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of the initMIP-Greenland intercomparison exercise is to compare, evaluate and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of 1) the initial present-day state of the ice sheet and 2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly), and should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap, but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

Millennial and orbital scale variability of the South American Monsoon during the penultimate glacial period

The presence of large, rapid climate oscillations is the most prominent feature of the Earth’s last glacial period. These oscillations are observed throughout the Northern Hemisphere and into the Southern Hemisphere tropics. Whether similar oscillations are typical of prior glacial periods, however, has not been well established. Here, we present results of a study of the South American Summer Monsoon system that covers nearly the entire penultimate glacial period, from 195 to 135 ky BP. We use a well-dated, high-resolution (~50 y) time series of oxygen isotopes to show that the precession of the earth’s orbit is the primary control on monsoon intensity. After removing the precession signal we observe millennial oscillations that are very similar in amplitude and structure to the Dansgaard/Oeschger cycles of the last interglacial and that match well a synthetic reconstruction of millennial variability. Time series analyses shows that the most prominent of the observed cycles occur at considerably longer frequency (~3500 y) that the Dansgaard/Oeschger cycles from Marine Isotope Stages 2–4.

Climate-driven ecological stability as a globally shared cause of Late Quaternary megafaunal extinctions: the Plaids and Stripes Hypothesis

Controversy persists about why so many large-bodied mammal species went extinct around the end of the last ice age. Resolving this is important for understanding extinction processes in general, for assessing the ecological roles of humans, and for conserving remaining megafaunal species, many of which are endangered today. Here we explore an integrative hypothesis that asserts that an underlying cause of Late Quaternary megafaunal extinctions was a fundamental shift in the spatio-temporal fabric of ecosystems worldwide. This shift was triggered by the loss of the millennial-scale climate fluctuations that were characteristic of the ice age but ceased approximately 11700 years ago on most continents. Under ice-age conditions, which prevailed for much of the preceding 2.6 Ma, these radical and rapid climate changes prevented many ecosystems from fully equilibrating with their contemporary climates. Instead of today's 'striped' world in which species' ranges have equilibrated with gradients of temperature, moisture, and seasonality, the ice-age world was a disequilibrial 'plaid' in which species' ranges shifted rapidly and repeatedly over time and space, rarely catching up with contemporary climate. In the transient ecosystems that resulted, certain physiological, anatomical, and ecological attributes shared by megafaunal species pre-adapted them for success. These traits included greater metabolic and locomotory efficiency, increased resistance to starvation, longer life spans, greater sensory ranges, and the ability to be nomadic or migratory. When the plaid world of the ice age ended, many of the advantages of being large were either lost or became disadvantages. For instance in a striped world, the low population densities and slow reproductive rates associated with large body size reduced the resiliency of megafaunal species to population bottlenecks. As the ice age ended, the downsides of being large in striped environments lowered the extinction thresholds of megafauna worldwide, which then increased the vulnerability of individual species to a variety of proximate threats they had previously tolerated, such as human predation, competition with other species, and habitat loss. For many megafaunal species, the plaid-to-stripes transition may have been near the base of a hierarchy of extinction causes whose relative importances varied geographically, temporally, and taxonomically.

In and out of glacial extremes by way of dust-climate feedbacks

Mineral dust aerosols cool Earth directly by scattering incoming solar radiation and indirectly by affecting clouds and biogeochemical cycles. Recent Earth history has featured quasi-100,000-y, glacial-interglacial climate cycles with lower/higher temperatures and greenhouse gas concentrations during glacials/interglacials. Global average, glacial maxima dust levels were more than 3 times higher than during interglacials, thereby contributing to glacial cooling. However, the timing, strength, and overall role of dust-climate feedbacks over these cycles remain unclear. Here we use dust deposition data and temperature reconstructions from ice sheet, ocean sediment, and land archives to construct dust-climate relationships. Although absolute dust deposition rates vary greatly among these archives, they all exhibit striking, nonlinear increases toward coldest glacial conditions. From these relationships and reconstructed temperature time series, we diagnose glacial-interglacial time series of dust radiative forcing and iron fertilization of ocean biota, and use these time series to force Earth system model simulations. The results of these simulations show that dust-climate feedbacks, perhaps set off by orbital forcing, push the system in and out of extreme cold conditions such as glacial maxima. Without these dust effects, glacial temperature and atmospheric CO₂ concentrations would have been much more stable at higher, intermediate glacial levels. The structure of residual anomalies over the glacial-interglacial climate cycles after subtraction of dust effects provides constraints for the strength and timing of other processes governing these cycles.

Pollen from the Deep-Sea: A Breakthrough in the Mystery of the Ice Ages

Pollen from deep-sea sedimentary sequences provides an integrated regional reconstruction of vegetation and climate (temperature, precipitation, and seasonality) on the adjacent continent. More importantly, the direct correlation of pollen, marine and ice indicators allows comparison of the atmospheric climatic changes that have affected the continent with the response of the Earth's other reservoirs, i.e., the oceans and cryosphere, without any chronological uncertainty. The study of long continuous pollen records from the European margin has revealed a changing and complex interplay between European climate, North Atlantic sea surface temperatures (SSTs), ice growth and decay, and high- and low-latitude forcing at orbital and millennial timescales. These records have shown that the amplitude of the last five terrestrial interglacials was similar above 40°N, while below 40°N their magnitude differed due to precession-modulated changes in seasonality and, particularly, winter precipitation. These records also showed that vegetation response was in dynamic equilibrium with rapid climate changes such as the Dangaard-Oeschger (D-O) cycles and Heinrich events, similar in magnitude and velocity to the ongoing global warming. However, the magnitude of the millennial-scale warming events of the last glacial period was regionally-specific. Precession seems to have imprinted regions below 40°N while obliquity, which controls average annual temperature, probably mediated the impact of D-O warming events above 40°N. A decoupling between high- and low-latitude climate was also observed within last glacial warm (Greenland interstadials) and cold phases (Greenland stadials). The synchronous response of western European vegetation/climate and eastern North Atlantic SSTs to D-O cycles was not a pervasive feature throughout the Quaternary. During periods of ice growth such as MIS 5a/4, MIS 11c/b and MIS 19c/b, repeated millennial-scale cold-air/warm-sea decoupling events occurred on the European margin superimposed to a long-term air-sea decoupling trend. Strong air-sea thermal contrasts promoted the production of water vapor that was then transported northward by the westerlies and fed ice sheets. This interaction between long-term and shorter time-scale climatic variability may have amplified insolation decreases and thus explain the Ice Ages. This hypothesis should be tested by the integration of stochastic processes in Earth models of intermediate complexity.

Enhanced Arctic Amplification Began at the Mid-Brunhes Event ~400,000 years ago

Arctic Ocean temperatures influence ecosystems, sea ice, species diversity, biogeochemical cycling, seafloor methane stability, deep-sea circulation, and CO₂ cycling. Today’s Arctic Ocean and surrounding regions are undergoing climatic changes often attributed to “Arctic amplification” – that is, amplified warming in Arctic regions due to sea-ice loss and other processes, relative to global mean temperature. However, the long-term evolution of Arctic amplification is poorly constrained due to lack of continuous sediment proxy records of Arctic Ocean temperature, sea ice cover and circulation. Here we present reconstructions of Arctic Ocean intermediate depth water (AIW) temperatures and sea-ice cover spanning the last ~ 1.5 million years (Ma) of orbitally-paced glacial/interglacial cycles (GIC). Using Mg/Ca paleothermometry of the ostracode Krithe and sea-ice planktic and benthic indicator species, we suggest that the Mid-Brunhes Event (MBE), a major climate transition ~ 400–350 ka, involved fundamental changes in AIW temperature and sea-ice variability. Enhanced Arctic amplification at the MBE suggests a major climate threshold was reached at ~ 400 ka involving Atlantic Meridional Overturning Circulation (AMOC), inflowing warm Atlantic Layer water, ice sheet, sea-ice and ice-shelf feedbacks, and sensitivity to higher post-MBE interglacial CO₂ concentrations.

Regional and global sea-surface temperatures during the last interglaciation

The last interglaciation (LIG, 129 to 116 thousand years ago) was the most recent time in Earth's history when global mean sea level was substantially higher than it is at present. However, reconstructions of LIG global temperature remain uncertain, with estimates ranging from no significant difference to nearly 2°C warmer than present-day temperatures. Here we use a network of sea-surface temperature (SST) records to reconstruct spatiotemporal variability in regional and global SSTs during the LIG. Our results indicate that peak LIG global mean annual SSTs were 0.5 ± 0.3°C warmer than the climatological mean from 1870 to 1889 and indistinguishable from the 1995 to 2014 mean. LIG warming in the extratropical latitudes occurred in response to boreal insolation and the bipolar seesaw, whereas tropical SSTs were slightly cooler than the 1870 to 1889 mean in response to reduced mean annual insolation.

Pyrite sulfur isotopes reveal glacial-interglacial environmental changes

The sulfur biogeochemical cycle plays a key role in regulating Earth's surface redox through diverse abiotic and biological reactions that have distinctive stable isotopic fractionations. As such, variations in the sulfur isotopic composition (δ³⁴S) of sedimentary sulfate and sulfide phases over Earth history can be used to infer substantive changes to the Earth's surface environment, including the rise of atmospheric oxygen. Such inferences assume that individual δ³⁴S records reflect temporal changes in the global sulfur cycle; this assumption may be well grounded for sulfate-bearing minerals but is less well established for pyrite-based records. Here, we investigate alternative controls on the sedimentary sulfur isotopic composition of marine pyrite by examining a 300-m drill core of Mediterranean sediments deposited over the past 500,000 y and spanning the last five glacial-interglacial periods. Because this interval is far shorter than the residence time of marine sulfate, any change in the sulfur isotopic record preserved in pyrite (δ³⁴S_pyr) necessarily corresponds to local environmental changes. The stratigraphic variations (>76‰) in the isotopic data reported here are among the largest ever observed in pyrite, and are in phase with glacial-interglacial sea level and temperature changes. In this case, the dominant control appears to be glacial-interglacial variations in sedimentation rates. These results suggest that there exist important but previously overlooked depositional controls on sedimentary sulfur isotope records, especially associated with intervals of substantial sea level change. This work provides an important perspective on the origin of variability in such records and suggests meaningful paleoenvironmental information can be derived from pyrite δ³⁴S records.

Timing and nature of AMOC recovery across Termination 2 and magnitude of deglacial CO2 change

Large amplitude variations in atmospheric CO₂ were associated with glacial terminations of the Late Pleistocene. Here we provide multiple lines of evidence suggesting that the ∼20 p.p.m.v. overshoot in CO₂ at the end of Termination 2 (T2) ∼129 ka was associated with an abrupt (≤400 year) deepening of Atlantic Meridional Overturning Circulation (AMOC). In contrast to Termination 1 (T1), which was interrupted by the Bølling-Allerød (B-A), AMOC recovery did not occur until the very end of T2, and was characterized by pronounced formation of deep waters in the NW Atlantic. Considering the variable influences of ocean circulation change on atmospheric CO₂, we suggest that the net change in CO₂ across the last 2 terminations was approximately equal if the transient effects of deglacial oscillations in ocean circulation are taken into account.

Differences in the sequence and timing of ocean circulation changes across glacial terminations could affect the magnitude of deglacial atmospheric CO₂ rise. Here, the authors argue that late ocean circulation recovery during the penultimate deglaciation (T2) led to a larger rise in CO₂ compared with T1.

State dependence of climatic instability over the past 720,000 years from Antarctic ice cores and climate modeling

Climatic variabilities on millennial and longer time scales with a bipolar seesaw pattern have been documented in paleoclimatic records, but their frequencies, relationships with mean climatic state, and mechanisms remain unclear. Understanding the processes and sensitivities that underlie these changes will underpin better understanding of the climate system and projections of its future change. We investigate the long-term characteristics of climatic variability using a new ice-core record from Dome Fuji, East Antarctica, combined with an existing long record from the Dome C ice core. Antarctic warming events over the past 720,000 years are most frequent when the Antarctic temperature is slightly below average on orbital time scales, equivalent to an intermediate climate during glacial periods, whereas interglacial and fully glaciated climates are unfavourable for a millennial-scale bipolar seesaw. Numerical experiments using a fully coupled atmosphere-ocean general circulation model with freshwater hosing in the northern North Atlantic showed that climate becomes most unstable in intermediate glacial conditions associated with large changes in sea ice and the Atlantic Meridional Overturning Circulation. Model sensitivity experiments suggest that the prerequisite for the most frequent climate instability with bipolar seesaw pattern during the late Pleistocene era is associated with reduced atmospheric CO2 concentration via global cooling and sea ice formation in the North Atlantic, in addition to extended Northern Hemisphere ice sheets.

The Atlantic Meridional Overturning Circulation and Abrupt Climate Change

Abrupt changes in climate have occurred in many locations around the globe over the last glacial cycle, with pronounced temperature swings on timescales of decades or less in the North Atlantic. The global pattern of these changes suggests that they reflect variability in the Atlantic meridional overturning circulation (AMOC). This review examines the evidence from ocean sediments for ocean circulation change over these abrupt events. The evidence for changes in the strength and structure of the AMOC associated with the Younger Dryas and many of the Heinrich events is strong. Although it has been difficult to directly document changes in the AMOC over the relatively short Dansgaard-Oeschger events, there is recent evidence supporting AMOC changes over most of these oscillations as well. The lack of direct evidence for circulation changes over the shortest events leaves open the possibility of other driving mechanisms for millennial-scale climate variability.

Strong coupling of Asian Monsoon and Antarctic climates on sub-orbital timescales

There is increasing evidence that millennial-scale climate variability played an active role on orbital-scale climate changes, but the mechanism for this remains unclear. A ²³⁰Th-dated stalagmite δ¹⁸O record between 88 and 22 thousand years (ka) ago from Yongxing Cave in central China characterizes changes in Asian monsoon (AM) strength. After removing the 65°N insolation signal from our record, the δ¹⁸O residue is strongly anti-phased with Antarctic temperature variability on sub-orbital timescales during the Marine Isotope Stage (MIS) 3. Furthermore, once the ice volume signal from Antarctic ice core records were removed and extrapolated back to the last two glacial-interglacial cycles, we observe a linear relationship for both short- and long-duration events between Asian and Antarctic climate changes. This provides the robust evidence of a link between northern and southern hemisphere climates that operates through changes in atmospheric circulation. We find that the weakest monsoon closely associated with the warmest Antarctic event always occurred during the Terminations. This finding, along with similar shifts in the opal flux record, suggests that millennial-scale events play a key role in driving the deglaciation through positive feedbacks associated with enhanced upwelling and increasing CO₂.

The Asian monsoon over the past 640,000 years and ice age terminations

Oxygen isotope records from Chinese caves characterize changes in both the Asian monsoon and global climate. Here, using our new speleothem data, we extend the Chinese record to cover the full uranium/thorium dating range, that is, the past 640,000 years. The record's length and temporal precision allow us to test the idea that insolation changes caused by the Earth's precession drove the terminations of each of the last seven ice ages as well as the millennia-long intervals of reduced monsoon rainfall associated with each of the terminations. On the basis of our record's timing, the terminations are separated by four or five precession cycles, supporting the idea that the '100,000-year' ice age cycle is an average of discrete numbers of precession cycles. Furthermore, the suborbital component of monsoon rainfall variability exhibits power in both the precession and obliquity bands, and is nearly in anti-phase with summer boreal insolation. These observations indicate that insolation, in part, sets the pace of the occurrence of millennial-scale events, including those associated with terminations and 'unfinished terminations'.

Metrological challenges for measurements of key climatological observables: Oceanic salinity and pH, and atmospheric humidity. Part 1: Overview

Water in its three ambient phases plays the central thermodynamic role in the terrestrial climate system. Clouds control Earth's radiation balance, atmospheric water vapour is the strongest "greenhouse" gas, and non-equilibrium relative humidity at the air-sea interface drives evaporation and latent heat export from the ocean. On climatic time scales, melting ice caps and regional deviations of the hydrological cycle result in changes of seawater salinity, which in turn may modify the global circulation of the oceans and their ability to store heat and to buffer anthropogenically produced carbon dioxide. In this paper, together with three companion articles, we examine the climatologically relevant quantities ocean salinity, seawater pH and atmospheric relative humidity, noting fundamental deficiencies in the definitions of those key observables, and their lack of secure foundation on the International System of Units, the SI. The metrological histories of those three quantities are reviewed, problems with their current definitions and measurement practices are analysed, and options for future improvements are discussed in conjunction with the recent seawater standard TEOS-10. It is concluded that the International Bureau of Weights and Measures, BIPM, in cooperation with the International Association for the Properties of Water and Steam, IAPWS, along with other international organisations and institutions, can make significant contributions by developing and recommending state-of-the-art solutions for these long standing metrological problems in climatology.

Eastern South African hydroclimate over the past 270,000 years

Processes that control the hydrological balance in eastern South Africa on orbital to millennial timescales remain poorly understood because proxy records documenting its variability at high resolution are scarce. In this work, we present a detailed 270,000 year-long record of terrestrial climate variability in the KwaZulu-Natal province based on elemental ratios of Fe/K from the southwest Indian Ocean, derived from X-ray fluorescence core scanning. Eastern South African climate variability on these time scales reflects both the long-term effect of regional insolation changes driven by orbital precession and the effects associated with high-latitude abrupt climate forcing over the past two glacial-interglacial cycles, including millennial-scale events not previously identified. Rapid changes towards more humid conditions in eastern South Africa as the Northern Hemisphere entered phases of extreme cooling were potentially driven by a combination of warming in the Agulhas Current and shifts of the subtropical anticyclones. These climate oscillations appear coherent with other Southern Hemisphere records but are anti-phased with respect to the East Asian Monsoon. Numerical modelling results reveal that higher precipitation in the KwaZulu-Natal province during precession maxima is driven by a combination of increased local evaporation and elevated moisture transport into eastern South Africa from the coast of Mozambique.

Icebergs not the trigger for North Atlantic cold events

Abrupt climate change is a ubiquitous feature of the Late Pleistocene epoch. In particular, the sequence of Dansgaard-Oeschger events (repeated transitions between warm interstadial and cold stadial conditions), as recorded by ice cores in Greenland, are thought to be linked to changes in the mode of overturning circulation in the Atlantic Ocean. Moreover, the observed correspondence between North Atlantic cold events and increased iceberg calving and dispersal from ice sheets surrounding the North Atlantic has inspired many ocean and climate modelling studies that make use of freshwater forcing scenarios to simulate abrupt change across the North Atlantic region and beyond. On the other hand, previous studies identified an apparent lag between North Atlantic cooling events and the appearance of ice-rafted debris over the last glacial cycle, leading to the hypothesis that iceberg discharge may be a consequence of stadial conditions rather than the cause. Here we further establish this relationship and demonstrate a systematic delay between pronounced surface cooling and the arrival of ice-rafted debris at a site southwest of Iceland over the past four glacial cycles, implying that in general icebergs arrived too late to have triggered cooling. Instead we suggest that--on the basis of our comparisons of ice-rafted debris and polar planktonic foraminifera--abrupt transitions to stadial conditions should be considered as a nonlinear response to more gradual cooling across the North Atlantic. Although the freshwater derived from melting icebergs may provide a positive feedback for enhancing and or prolonging stadial conditions, it does not trigger northern stadial events.

Bipolar seesaw control on last interglacial sea level

Our current understanding of ocean-atmosphere-cryosphere interactions at ice-age terminations relies largely on assessments of the most recent (last) glacial-interglacial transition, Termination I (T-I). But the extent to which T-I is representative of previous terminations remains unclear. Testing the consistency of termination processes requires comparison of time series of critical climate parameters with detailed absolute and relative age control. However, such age control has been lacking for even the penultimate glacial termination (T-II), which culminated in a sea-level highstand during the last interglacial period that was several metres above present. Here we show that Heinrich Stadial 11 (HS11), a prominent North Atlantic cold episode, occurred between 135 ± 1 and 130 ± 2 thousand years ago and was linked with rapid sea-level rise during T-II. Our conclusions are based on new and existing data for T-II and the last interglacial that we collate onto a single, radiometrically constrained chronology. The HS11 cold episode punctuated T-II and coincided directly with a major deglacial meltwater pulse, which predominantly entered the North Atlantic Ocean and accounted for about 70 per cent of the glacial-interglacial sea-level rise. We conclude that, possibly in response to stronger insolation and CO2 forcing earlier in T-II, the relationship between climate and ice-volume changes differed fundamentally from that of T-I. In T-I, the major sea-level rise clearly post-dates Heinrich Stadial 1. We also find that HS11 coincided with sustained Antarctic warming, probably through a bipolar seesaw temperature response, and propose that this heat gain at high southern latitudes promoted Antarctic ice-sheet melting that fuelled the last interglacial sea-level peak.

Permafrost thawing as a possible source of abrupt carbon release at the onset of the Bølling/Allerød

One of the most abrupt and yet unexplained past rises in atmospheric CO₂ (>10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, ~14,600 years ago. Here we use a U/Th-dated record of atmospheric Δ¹⁴C from Tahiti corals to provide an independent and precise age control for this CO₂ rise. We also use model simulations to show that the release of old (nearly ¹⁴C-free) carbon can explain these changes in CO₂ and Δ¹⁴C. The Δ¹⁴C record provides an independent constraint on the amount of carbon released (~125 Pg C). We suggest, in line with observations of atmospheric CH₄ and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

Ice core records show evidence for an abrupt, and thus far unexplained, increase in atmospheric CO₂ levels ~14,600 years ago. Here, the authors combine ice core data, a precisely dated decline in atmospheric ¹⁴C and numerical simulations, and propose thawing permafrost as a possible source of this event.

Oceanographic dynamics and the end of the last interglacial in the subpolar North Atlantic

The last interglacial interval was terminated by the inception of a long, progressive glaciation that is attributed to astronomically influenced changes in the seasonal distribution of sunlight over the earth. However, the feedbacks, internal dynamics, and global teleconnections associated with declining northern summer insolation remain incompletely understood. Here we show that a crucial early step in glacial inception involves the weakening of the subpolar gyre (SPG) circulation of the North Atlantic Ocean. Detailed new records of microfossil foraminifera abundance and stable isotope ratios in deep sea sediments from Ocean Drilling Program site 984 south of Iceland reveal repeated, progressive cold water-mass expansions into subpolar latitudes during the last peak interglacial interval, marine isotope substage 5e. These movements are expressed as a sequence of progressively extensive southward advances and subsequent retreats of a hydrographic boundary that may have been analogous to the modern Arctic front, and associated with rapid changes in the strength of the SPG. This persistent millennial-scale oceanographic oscillation accompanied a long-term cooling trend at a time of slowly declining northern summer insolation, providing an early link in the propagation of those insolation changes globally, and resulting in a rapid transition from extensive regional warmth to the dramatic instability of the subsequent ∼ 100 ka.