Dansgaard Oeschger Research Articles

Estimating biases during detection of leads and lags between climate elements across Dansgaard–Oeschger events

Abstract. Dansgaard–Oeschger (DO) events occurred throughout the last glacial period. Greenland ice cores show a rapid warming during each stadial to interstadial transition, alongside an abrupt loss of sea ice and a major reorganisation of the atmospheric circulation. Other records also indicate simultaneous abrupt changes to the oceanic circulation. Recently, an advanced Bayesian ramp-fitting method has been developed and used to investigate time lags between transitions in these different climate elements with a view to determining the relative order of these changes. Here, we critically review this method in both its original implementation and a new, extended implementation. Using ice core data, climate model output, and carefully synthesised data representing DO events, we demonstrate that both implementations of the method suffer from biases of up to 15 years. These biases mean that the method will tend to yield transition onsets that are too early. Further investigation of DO warming event records in climate models and ice core data reveals that the biases are on the same order of magnitude as potential timing differences between the abrupt transitions of different climate elements. Additionally, we find that higher-resolution records would not reduce these biases. We conclude that decadal-scale leads and lags between climate elements across DO events cannot be reliably detected, as we cannot exclude the possibility that they result solely from the biases we present here.

The Greenland spatial fingerprint of Dansgaard–Oeschger events in observations and models

Pleistocene Ice Ages display abrupt Dansgaard-Oeschger (DO) climate oscillations that provide prime examples of Earth System tipping points-abrupt transition that may result in irreversible change. Greenland ice cores provide key records of DO climate variability, but gas-calibrated estimates of the temperature change magnitudes have been limited to central and northwest Greenland. Here, we present ice-core δ15N-N2 records from south (Dye 3) and coastal east Greenland (Renland) to calibrate the local water isotope thermometer and provide a Greenland-wide spatial characterization of DO event magnitude. We combine these data with existing records of δ18O, deuterium excess, and accumulation rates to create a multiproxy "fingerprint" of the DO impact on Greenland. Isotope-enabled climate models have skill in simulating the observational multiproxy DO event impact, and we use a series of idealized simulations with such models to identify regions of the North Atlantic that are critical in explaining DO variability. Our experiments imply that wintertime sea ice variation in the subpolar gyre, rather than the commonly invoked Nordic Seas, is both a sufficient and a necessary condition to explain the observed DO impacts in Greenland, whatever the distal cause. Moisture-tagging experiments support the idea that Greenland DO isotope signals may be explained almost entirely via changes in the vapor source distribution and that site temperature is not a main control on δ18O during DO transitions, contrary to the traditional interpretation. Our results provide a comprehensive, multiproxy, data-model synthesis of abrupt DO climate variability in Greenland.

Variations in deep-sea methane seepage linked to millennial-scale changes in bottom water temperatures ~ 50–6 ka, NW Svalbard margin

During the last glaciation, the northern hemisphere experienced profound millennial-scale changes (termed Dansgaard-Oeschger (DO) events) in atmospheric and oceanic temperatures. In the North Atlantic, the fluctuations resulted in extremely unstable bottom water conditions with bottom water temperatures (BWT) varying up to > 5 °C. We have studied these changes in a core from 1,300 m water depth at Vestnesa Ridge, northwestern Svalbard margin to investigate a possible connection between BWT and seepage of methane from the seafloor covering the period ~ 50–6 ka. Beneath Vestnesa Ridge, gas hydrates containing vast amounts of methane are kept stable due to the high pressure and low temperatures. Release of gas is shown by numerous pockmarks on the seafloor. The pockmarks at 1,300 m water depth are presently inactive, but they bear witness of earlier activity. Our study shows that from ~ 50–6 ka, the core site experienced repeated increases in BWT and in the emissions of gas, both following the pattern of the DO events. This correspondence in time scale indicates that BWT was the primary forcing factor for the variability in methane release. However, the releases were delayed by up to > 1,000 years compared to the initial increase in BWT.

Characterisation of Dansgaard–Oeschger events in palaeoclimate time series using the matrix profile method

Abstract. Palaeoclimate time series, reflecting the state of Earth's climate in the distant past, occasionally display very large and rapid shifts showing abrupt climate variability. The identification and characterisation of these abrupt transitions in palaeoclimate records is of particular interest as this allows for understanding of millennial climate variability and the identification of potential tipping points in the context of current climate change. Methods that are able to characterise these events in an objective and automatic way, in a single time series, or across two proxy records are therefore of particular interest. In our study the matrix profile approach is used to describe Dansgaard–Oeschger (DO) events, abrupt warmings detected in the Greenland ice core, and Northern Hemisphere marine and continental records. The results indicate that canonical events DO-19 and DO-20, occurring at around 72 and 76 ka, are the most similar events over the past 110 000 years. These transitions are characterised by matching transitions corresponding to events DO-1, DO-8, and DO-12. They are abrupt, resulting in a rapid shift to warmer conditions, followed by a gradual return to cold conditions. The joint analysis of the δ18O and Ca2+ time series indicates that the transition corresponding to the DO-19 event is the most similar event across the two time series.

An overview of Upper Pleistocene coastal deposits on Mallorca island

This review paper delves into the geological history of Mallorca, focusing on the Middle to Late Pleistocene period and its impact on the landscape and climate of the island. The Quaternary, particularly, the Last Interglacial, represents a relatively short timeframe on the global chronological scale, yet it is notable for its climatic instability and complex alternation of glacial and interglacial phases, leading to fluctuations in sea levels. During the last interglacial, millennial-scale fluctuations, known as Dansgaard–Oeschger and Heinrich events, occurred. Aeolianites are thought to have formed during periods of sparse vegetation cover and abundant sediment availability from the sea bed or platform (Heinrich events), while colluvial deposits and paleosols are believed to represent periods of reduced sediment supply and enhanced pedogenesis (Dansgaard–Oeschger events). These different deposits reflect climatic shifts and sea-level changes, offering insights into the environmental history of the island. Fossil evidence, including Strombus bubonius (currently known as Thetystrombus lauts) or “Senegalese hosts,” reveals past warm climatic fluctuations, particularly during the Last Interglacial. Myotragus balearicus exemplifies evolutionary adaptation and isolation on the island, offering a unique perspective on its paleontological legacy. The integration of stratigraphic characteristics and literature research on Optically Stimulated Luminescence (OSL) dating techniques provides a comprehensive temporal framework, spanning from Marine Isotope Stage 6 (MIS 6) to Marine Isotope Stage 2 (MIS 2), enabling the precise dating of geological events and a stratigraphic correlation.

The Impact of CO2 and Climate State on Whether Dansgaard–Oeschger Type Oscillations Occur in Climate Models

AbstractGreenland ice core records feature Dansgaard–Oeschger (D‐O) events, which are abrupt warming episodes followed by gradual cooling during ice age climate. The three climate models used in this study (CCSM4, MPI‐ESM, and HadCM3) show spontaneous self‐sustained D‐O‐like oscillations (albeit with differences in amplitude, duration, and shape) in a remarkably similar, narrow window of carbon dioxide (CO2) concentration, roughly 185–230 ppm. This range matches atmospheric CO2 during Marine Isotopic Stage 3 (MIS 3: between 27.8 and 59.4 thousand of years BP, hereafter ka), a period when D‐O events were most frequent. Insights from the three climate models point to North Atlantic (NA) sea‐ice coverage as a key ingredient behind D‐O type oscillations, which acts as a “tipping element.” Other climate state properties such as Mean Atlantic Meridional Overturning Circulation strength, global mean temperature and salinity gradient in the Atlantic Ocean do not determine whether D‐O type behavior can occur in all three models.

Vegetation response in NW Mediterranean borderlands to the millennial-scale climate variability of the last glacial period

Deep-sea pollen records from the Western European margin show that during the Last Glacial period (115-27 ka), regional vegetation oscillated between steppe and open forest in response to the millennial scale climate variability, Dansgaard-Oeschger (D-O) cycles and Heinrich events (HE), and that the magnitude of the forest expansions during D-O warming events was modulated by orbital parameters. However, little is known about the vegetation response in the northwestern Mediterranean region. Here, we present a new well-chronologically constrained high-resolution marine pollen record from the Gulf of Lion documenting the NW Mediterranean vegetation response during Marine Isotope Stages (MIS) 4 and 3. The pollen record highlights the strongest forest expansions during D-O 17-16 and D-O 8 in NW Mediterranean borderlands, likely driven by minima in precession. In this region, Heinrich Stadials (HS) 5 and 4 are marked by steppe landscapes while temperate forest surprisingly characterized HS 6. We hypothesize that the combination of minima in precession and local atmospheric and marine processes allows the development of the temperate forest in NW Mediterranean while the expansion of open environments occurred in other areas of Western Europe. In addition, our results support the hypothesis of Picea survival in NW Mediterranean during MIS 4 to 3.

On the representation of multiplicative noise in modeling Dansgaard–Oeschger events

The interpretation of multiplicative noise in a stochastic differential equation in the context of data-driven inverse modeling is discussed. Application to the well-known paleoclimate phenomenon of Dansgaard–Oeschger events leads to qualitatively different ‘climate potentials’ in the case of the Itô or the Stratonovich interpretation of the stochastic integral. While a physical model is endowed with an interpretation from construction, whether implicitly or explicitly, inverse models derived from data do not afford such a luxury. In this case, a physical model must accompany the mathematical model equation in order to be able to choose a stochastic interpretation. This case study illustrates the differences between the two representations of stochastic noise and demonstrates the need for input from physical constraints when constructing conceptual stochastic models of the observed climate records.

Reconstruction of warm-season temperatures in central Europe during the past 60 000 years from lacustrine branched glycerol dialkyl glycerol tetraethers (brGDGTs)

Abstract. Millennial-scale climate variations during the last glacial period, such as Dansgaard–Oeschger (DO) cycles and Heinrich events, have been extensively studied using ice core and marine proxy records. However, there is a limited understanding of the magnitude of these temperature fluctuations in continental regions, and questions remain about the seasonal signal of these climate events. This study presents a 60 000-year-long temperature reconstruction based on branched glycerol dialkyl glycerol tetraethers (brGDGTs) extracted from lake sediments from the Eifel Volcanic Field, Germany. brGDGTs are bacterial membrane-spanning lipids that are known to have a strong relationship with temperature, making them suitable for temperature reconstructions. We test several temperature calibration models on modern samples taken from soils and multiple maar lakes. We find a negative bias in brGDGT-based temperature estimates associated with water depth and anoxic conditions that can be corrected for by accounting for a brGDGT isomer that is only produced in anoxic conditions. The corrected temperature reconstruction correlates with proxy and climate model estimates of temperature spanning the same time period, validating the calibration approach we selected. However, millennial-scale variability is significantly dampened in the brGDGT record, and in contrast to other Northern Hemisphere climate records, during several Heinrich stadials, temperatures actually increase. We demonstrate that these apparent discrepancies can be explained by the unique seasonal response of the brGDGT paleothermometer to temperatures of months above freezing (TMAF). Our data support the view that warm-season temperatures in Europe varied minimally during the last glacial period and that abrupt millennial-scale events were defined by colder, longer winters. Our continuous high-resolution temperature reconstruction provides important information about the magnitude of seasonal climate variability during the last glacial period that can be used to test climate models and inform studies of paleoecological change.

Statistical precursor signals for Dansgaard–Oeschger cooling transitions

Abstract. Given growing concerns about climate tipping points and their risks, it is important to investigate the capability of identifying robust precursor signals for the associated transitions. In general, the variance and short-lag autocorrelations of the fluctuations increase in a stochastically forced system approaching a critical or bifurcation-induced transition, making them theoretically suitable indicators to warn of such transitions. Paleoclimate records provide useful test beds if such a warning of a forthcoming transition could work in practice. The Dansgaard–Oeschger (DO) events are characterized by millennial-scale abrupt climate changes during the glacial period, manifesting most clearly as abrupt temperature shifts in the North Atlantic region. Some previous studies have found such statistical precursor signals for the DO warming transitions. On the other hand, statistical precursor signals for the abrupt DO cooling transitions have not been identified. Analyzing Greenland ice core records, we find robust and statistically significant precursor signals of DO cooling transitions in most of the interstadials longer than roughly 1500 years but not in the shorter interstadials. The origin of the statistical precursor signals is mainly related to so-called rebound events, humps in the temperature observed at the end of interstadial, some decades to centuries prior to the actual transition. We discuss several dynamical mechanisms that give rise to such rebound events and statistical precursor signals.

Distinct oxygenation modes of the Gulf of Oman over the past 43 000 years – a multi-proxy approach

Abstract. Changing climatic conditions can shape the strength and extent of the oxygen minimum zone (OMZ). The presence and variability of the OMZ in the Arabian Sea is of importance to the latter's ecosystem. The state of oxygenation has, for instance, an impact on the pelagic and benthic faunal community or the nitrogen and carbon cycles. It is important to understand the dynamics of the OMZ and related marine environmental conditions because of their climate feedbacks. In this study, we combined three independent proxies to reconstruct the oxygenation state of the water column and bottom water in the Gulf of Oman for the past 43 kyr approximately. This multi-proxy approach is done for the first time at the northeastern Oman margin located in the Gulf of Oman. We used bulk sedimentary nitrogen isotopes (δ15N) and the alkane ratio (lycopane +n-C35)/n-C31 and benthic foraminiferal faunal analysis to reconstruct the strength of the OMZ in the water column and bottom water oxygenation, respectively. Our results show that the Gulf of Oman experienced strong pronounced OMZ and bottom water deoxygenation during the Holocene. In contrast, during Marine Isotope Stage 2 (MIS 2), including the Last Glacial Maximum (LGM), the Gulf of Oman was very well ventilated, with a highly diverse benthic foraminiferal community. This may have been caused by stronger wind-induced mixing and better ventilation by oxygen-rich water masses. Our results also show moderate oxygenation during MIS 3, with deoxygenation events during most of the warmer Dansgaard–Oeschger (D–O) events. We propose two distinct oxygenation modes for the Gulf of Oman: (1) a stable period of either strongly pronounced water column OMZ and bottom water deoxygenation or well-oxygenated water column and bottom water conditions and (2) an unstable period of oscillating oxygenation states between moderately oxygenated (stadials) and deoxygenated (interstadials, D–O events) conditions. The unstable period may be triggered by an interstadial Atlantic meridional overturning circulation (AMOC) mode, which is required to initiate D–O events.

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.

Changes in the Red Sea overturning circulation during Marine Isotope Stage 3

Abstract. The oceanography of the Red Sea is controlled by the restricted exchange of water masses with the Indian Ocean and by high evaporation rates due to the arid climate of the surrounding land areas. In the northern Red Sea, the formation of oxygen-rich subsurface water ventilates the deeper parts of the basin, but little is known about the variability in this process in the past. The stable oxygen and carbon isotope records of epibenthic foraminifera from a sediment core of the central Red Sea and comparison with existing isotope records allow for the reconstruction of changes in the Red Sea overturning circulation (ROC) during Marine Isotope Stage 3. The isotope records imply millennial-scale variations in the ROC, in phase with the climate variability in the high northern latitudes. This suggests an immediate response of dense-water formation to the regional climate and hydrology of the northern Red Sea. Deep-water formation was intensified under the influence of cold and hyper-arid conditions during Heinrich stadials and was diminished during Dansgaard–Oeschger interstadials. While these changes are reflected in both stable oxygen and carbon isotope records, the latter data also exhibit changes in phase with the African–Indian monsoon system. The decoupling of the stable carbon and oxygen isotope records at the summer monsoon maximum centered around 55–60 ka may be associated with an increased inflow of nutrient-rich intermediate waters from the Arabian Sea to the central Red Sea. This process fueled local surface water productivity, resulting in enhanced remineralization of sinking organic matter and release of 12C at intermediate water depths.

Millennial-scale changes in the environment, ice conditions, and deep-water ventilation at the Detroit–Tenji Seamounts, Northwest Pacific, over the last 43 kyr

The Western Subarctic Gyre (WSG) in the North Pacific represents the termination of the Atlantic Meridional Overturning Circulation (AMOC). It has been strongly influenced by periodic glaciations of adjacent landmasses, orbital- and millennial-scale variability of its environment, and ventilation, which, to date, have been poorly studied. In this study, we investigated sediment core LV76-21 recovered from the Detroit Seamount using several productivity and lithological proxies, δ18O records of planktic and benthic foraminifera (δ18Opf and δ18Obf), AMS 14C data. We also used the published results from three well-dated cores from the Detroit and Tenji Seamounts. Comparison of the δ18Opf and δ18Obf records allowed us to establish iceberg discharge events 1 and 1’ at 39–36.8 and 42–40 ka, respectively, in the WSG over the last 43 kyr. This was accompanied by considerable lightning of δ18Opf relative to the coeval δ18Obf values. The lithological results have also shown variable sea-ice influence on the regional environment during Marine Isotope Stage 3 and particularly 2, with sea ice extensions at 28.5–25.5 ka and 17.6–16.6 ka, and its significant loss since the Bølling–Allerød warming. We used high resolution δ18Obf records from studied cores as a subtle proxy of contribution of North Atlantic Deep Water (NADW) to Lower Circumpolar Deep Water (LCDW), which originating from around Antarctica and have slightly different δ18Obf values, in order to reconstruction of its spatial variability in the context of regional deep-water ventilation. The contribution of NADW to LCDW in the ventilation of the western flank of Detroit Seamount significantly increased during the interglacial and long-lasting Dansgaard–Oeschger interstadials, which was synchronous with the AMOC intensification. Over the last 43 kyr, the deep water of its eastern flank has been constantly ventilated by the LCDW with low contribution of NADW.

Last glacial hydroclimate variability in the Yucatán Peninsula not just driven by ITCZ shifts

We reconstructed hydroclimate variability in the Yucatán Peninsula (YP) based on stalagmite oxygen and carbon isotope records from a well-studied cave system located in the northeastern YP, a region strongly influenced by Caribbean climate dynamics. The new stalagmite isotopic records span the time interval between 43 and 26.6 ka BP, extending a previously published record from the same cave system covering the interval between 26.5 and 23.2 ka BP. Stalagmite stable isotope records show dominant decadal and multidecadal variability, and weaker variability on millennial timescales. These records suggest significant precipitation declines in the broader Caribbean region during Heinrich events 4 and 3 of ice-rafted discharge into the North Atlantic, in agreement with the antiphase pattern of precipitation variability across the equator suggested by previous studies. On millennial timescales, the stalagmite isotope records do not show the distinctive saw-tooth pattern of climate variability observed in Greenland during Dansgaard–Oeschger (DO) events, but a pattern similar to North Atlantic sea surface temperature (SST) variability. We propose that shifts in the mean position of the Intertropical Convergence Zone (ITCZ), per se, are not the dominant driver of last glacial hydroclimate variability in the YP on millennial timescales but instead that North Atlantic SSTs played a dominant role. Our results support a negative climate feedback mechanism whereby large low latitude precipitation deficits resulting from AMOC slowdown would lead to elevated salinity in the Caribbean and ultimately help reactivate AMOC and Caribbean precipitation. However, because of the unique drivers of future climate in the region, predicted twenty-first century YP precipitation reductions are unlikely to be modulated by this negative feedback mechanism.

Global reorganization of atmospheric circulation during Dansgaard–Oeschger cycles

Ice core records from Greenland provide evidence for multiple abrupt cold-warm-cold events recurring at millennial time scales during the last glacial interval. Although climate variations resembling Dansgaard-Oeschger (DO) oscillations have been identified in climate archives across the globe, our understanding of the climate and ecosystem impacts of the Greenland warming events in lower latitudes remains incomplete. Here, we investigate the influence of DO-cold-to-warm transitions on the global atmospheric circulation pattern. We comprehensively analyze δ18O changes during DO transitions in a globally distributed dataset of speleothems and set those in context with simulations of a comprehensive high-resolution climate model featuring internal millennial-scale variations of similar magnitude. Across the globe, speleothem δ18O signals and model results indicate consistent large-scale changes in precipitation amount, moisture source, or seasonality of precipitation associated with the DO transitions, in agreement with northward shifts of the Hadley circulation. Furthermore, we identify a decreasing trend in the amplitude of DO transitions with increasing distances from the North Atlantic region. This provides quantitative observational evidence for previous suggestions of the North Atlantic region being the focal point for these archetypes of past abrupt climate changes.

Disparate energy sources for slow and fast Dansgaard–Oeschger cycles

Abstract. During the Late Pleistocene, Dansgaard–Oeschger (DO) cycles triggered warming events that were as abrupt as the present-day human-induced warming. However, in the absence of a periodic forcing operating on millennial timescales, the main energy sources of DO cycles remain debated. Here, we identify the energy sources of DO cycles by applying a bispectral analysis to the North Greenland Ice Core Project (NGRIP) oxygen isotope (δ18Oice) record; a 123 kyr long proxy record of air temperatures (Tair) over Greenland. For both modes of DO cyclicity – slow and fast – we detect disparate energy sources. Slow DO cycles, marked by multi-millennial periodicities in the 12.5 to 2.5 kyr bandwidth, receive energy from astronomical periodicities. Fast DO cycles, characterized by millennial periodicities in the 1.5 ± 0.5 kyr range, receive energy from centennial periodicities. We propose cryospheric and oceanic mechanisms that facilitate the transfer of energy from known sources to slow and fast DO cycles, respectively. Our findings stress the importance of understanding energy-transfer mechanisms across a broad range of timescales to explain the origins of climate cycles without primary periodic energy sources.

Aeolian process and climatic changes in loess records from the eastern Tibetan Plateau: Implications for paleoenvironmental dynamics since MIS 3

The Tibetan Plateau (TP) is extremely sensitive to climate change. However, climate variability during the last glaciation in the TP is not well-understood. The widely distributed loess deposits on the eastern TP serve as a significant archive for retrieving past climate changes and aeolian dust dynamics. Through detailed quartz optically stimulated luminescence and radiocarbon dating and the analysis of environmental proxies, we have established a reliable chronology of the Maerkang loess sequence in the eastern TP and reconstructed well-constrained histories of dust activities and environmental variation. Our results demonstrate that the accumulation of aeolian loess in this region occurred at least 60 thousand years ago. The dust mass accumulation rate (MAR) varied highly in the last glacial period, with high MARs in the Last Glacial Maximum and relatively low MARs in marine isotope stage (MIS) 3. Humid conditions occurred in MIS 3 (approximately 50–32 ka) due to the intensified Indian summer monsoon. The rapid climate variations of the Maerkang loess might be correlated to the Dansgaard–Oeschger cycles and the Heinrich events, which implies the significant impact of the North Atlantic climatic system on the TP. Our findings provide new evidence for a better understanding of the history of environmental variations since MIS 3 in the eastern TP.

The Problem with the Anthropocene: Kainos, Not Anthropos

The Problem with the Anthropocene: <i>Kainos</i>, Not <i>Anthropos</i>

Stable stadial and interstadial states of the last glacial's climate identified in a combined stable water isotope and dust record from Greenland

Abstract. During the last glacial interval, the Northern Hemisphere climate was punctuated by a series of abrupt changes between two characteristic climate regimes. The existence of stadial (cold) and interstadial (milder) periods is typically attributed to a hypothesised bistability in the glacial North Atlantic climate system, allowing for rapid transitions from the stadial to the interstadial state – the so-called Dansgaard–Oeschger (DO) events – and more gradual yet still fairly abrupt reverse shifts. The physical mechanisms driving these regime transitions remain debated. DO events are characterised by substantial warming over Greenland and a reorganisation of the Northern Hemisphere atmospheric circulation, which are evident from concomitant shifts in the δ18O ratios and dust concentration records from Greenland ice cores. Treating the combined δ18O and dust record obtained by the North Greenland Ice Core Project (NGRIP) as a realisation of a two-dimensional, time-homogeneous, and Markovian stochastic process, we present a reconstruction of its underlying deterministic drift based on the leading-order terms of the Kramers–Moyal equation. The analysis reveals two basins of attraction in the two-dimensional state space that can be identified with the stadial and interstadial regimes. The drift term of the dust exhibits a double-fold bifurcation structure, while – in contrast to prevailing assumptions – the δ18O component of the drift is clearly mono-stable. This suggests that the last glacial's Greenland temperatures should not be regarded as an intrinsically bistable climate variable. Instead, the two-regime nature of the δ18O record is apparently inherited from a coupling to another bistable climate process. In contrast, the bistability evidenced in the dust drift points to the presence of two stable circulation regimes of the last glacial's Northern Hemisphere atmosphere.