Dansgaard-Oeschger Events Research Articles

Early-warning signals for Dansgaard-Oeschger events in a high-resolution ice core record

The Dansgaard–Oeschger (DO) events, as observed in oxygen isotope ratios from the North Greenland Ice Core Project (NGRIP) record, are an outstanding example of past abrupt climate transitions. Their physical cause remains debated, and previous research indicated that they are not preceded by classical early-warning signals (EWS). Subsequent research hypothesized that the DO events are caused by bifurcations of physical mechanisms operating at decadal timescales, and proposed to search for EWS in the high-frequency fluctuation levels. Here, a time series with 5-year resolution is obtained from the raw NGRIP record, and significant numbers of EWS in terms of variance and autocorrelation increases are revealed in the decadal-scale variability. Wavelet analysis indicates that the EWS are most pronounced in the 10–50-year periodicity band, confirming the above hypothesis. The DO events are hence neither directly noise-induced nor purely externally forced, which provides valuable constraints regarding potential physical causes.

A spatiotemporal reconstruction of sea-surface temperatures in the North Atlantic during Dansgaard–Oeschger events 5–8

Abstract. Here, we establish a spatiotemporal evolution of the sea-surface temperatures in the North Atlantic over Dansgaard–Oeschger (DO) events 5–8 (approximately 30–40 kyr) using the proxy surrogate reconstruction method. Proxy data suggest a large variability in North Atlantic sea-surface temperatures during the DO events of the last glacial period. However, proxy data availability is limited and cannot provide a full spatial picture of the oceanic changes. Therefore, we combine fully coupled, general circulation model simulations with planktic foraminifera based sea-surface temperature reconstructions to obtain a broader spatial picture of the ocean state during DO events 5–8. The resulting spatial sea-surface temperature patterns agree over a number of different general circulation models and simulations. We find that sea-surface temperature variability over the DO events is characterized by colder conditions in the subpolar North Atlantic during stadials than during interstadials, and the variability is linked to changes in the Atlantic Meridional Overturning circulation and in the sea-ice cover. Forced simulations are needed to capture the strength of the temperature variability and to reconstruct the variability in other climatic records not directly linked to the sea-surface temperature reconstructions. This is the first time the proxy surrogate reconstruction method has been applied to oceanic variability during MIS3. Our results remain robust, even when age uncertainties of proxy data, the number of available temperature reconstructions, and different climate models are considered. However, we also highlight shortcomings of the methodology that should be addressed in future implementations.

Using Holo-Hilbert spectral analysis to quantify the modulation of Dansgaard-Oeschger events by obliquity

Astronomical forcing (obliquity and precession) has been thought to modulate Dansgaard-Oeschger (DO) events, yet the detailed quantification of such modulations has not been examined. In this study, we apply the novel Holo-Hilbert Spectral Analysis (HHSA) to five polar ice core records, quantifying astronomical forcing's time-varying amplitude modulation of DO events and identifying the preferred obliquity phases for large amplitude modulations. The unique advantages of HHSA over the widely used windowed Fourier spectral analysis for quantifying astronomical forcing's nonlinear modulations of DO events is first demonstrated with a synthetic data that closely resembles DO events recorded in Greenland ice cores (NGRIP, GRIP, and GISP2 cores on GICC05 modelext timescale). The analysis of paleoclimatic proxies show that statistically significantly more frequent DO events, with larger amplitude modulation in the Greenland region, tend to occur in the decreasing phase of obliquity, especially from its mean value to its minimum value. In the eastern Antarctic, although statistically significantly more DO events tend to occur in the decreasing obliquity phase in general, the preferred phase of obliquity for large amplitude modulation on DO events is a segment of the increasing phase near the maximum obliquity, implying that the physical mechanisms of DO events may be different for the two polar regions. Additionally, by using cross-spectrum and magnitude-squared analyses, Greenland DO mode at a timescale of about 1400 years leads the Antarctic DO mode at the same timescale by about 1000 years.

Novel automated inversion algorithm for temperature reconstruction using gas isotopes from ice cores

Abstract. Greenland past temperature history can be reconstructed by forcing the output of a firn-densification and heat-diffusion model to fit multiple gas-isotope data (δ15N or δ40Ar or δ15Nexcess) extracted from ancient air in Greenland ice cores using published accumulation-rate (Acc) datasets. We present here a novel methodology to solve this inverse problem, by designing a fully automated algorithm. To demonstrate the performance of this novel approach, we begin by intentionally constructing synthetic temperature histories and associated δ15N datasets, mimicking real Holocene data that we use as “true values” (targets) to be compared to the output of the algorithm. This allows us to quantify uncertainties originating from the algorithm itself. The presented approach is completely automated and therefore minimizes the “subjective” impact of manual parameter tuning, leading to reproducible temperature estimates. In contrast to many other ice-core-based temperature reconstruction methods, the presented approach is completely independent from ice-core stable-water isotopes, providing the opportunity to validate water-isotope-based reconstructions or reconstructions where water isotopes are used together with δ15N or δ40Ar. We solve the inverse problem T(δ15N, Acc) by using a combination of a Monte Carlo based iterative approach and the analysis of remaining mismatches between modelled and target data, based on cubic-spline filtering of random numbers and the laboratory-determined temperature sensitivity for nitrogen isotopes. Additionally, the presented reconstruction approach was tested by fitting measured δ40Ar and δ15Nexcess data, which led as well to a robust agreement between modelled and measured data. The obtained final mismatches follow a symmetric standard-distribution function. For the study on synthetic data, 95 % of the mismatches compared to the synthetic target data are in an envelope between 3.0 to 6.3 permeg for δ15N and 0.23 to 0.51 K for temperature (2σ, respectively). In addition to Holocene temperature reconstructions, the fitting approach can also be used for glacial temperature reconstructions. This is shown by fitting of the North Greenland Ice Core Project (NGRIP) δ15N data for two Dansgaard–Oeschger events using the presented approach, leading to results comparable to other studies.

Abrupt or gradual? Change point analysis of the late Pleistocene–Holocene climate record from Chew Bahir, southern Ethiopia

AbstractWe used a change point analysis on a late Pleistocene–Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from ~15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from ~4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than ~100 yr) at ~13,260 yr and lasted until ~11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between ~45–20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.

Temperature and seawater isotopic controls on two stalagmite records since 83 ka from maritime Japan

Millennial-scale interstadial Dansgaard-Oeschger (D-O) cycles and Heinrich (H) stadial events are pronounced paleoclimatic features during the last glacial period, which were first demonstrated in the North Atlantic region. These stadial and interstadial events are expressed in marine and terrestrial high-resolution records elsewhere in the world, but the magnitude and mode of the regional climate changes are still poorly quantified. Here we present new replicated stalagmite δ18O profiles from two caves in central Japan, which extend back to 83.4 ka. The records clearly display the H7 to H3 events, but not D-O cycles. An important feature of the two Japanese stalagmites is the small difference (∼2.9‰) in δ18O values between the mid-Holocene and the Last Glacial Maximum (LGM). Long-term trends of the stalagmite δ18O values at the more maritime site generally follow that of the δ18O record of seawater, which is responsible for ∼1.1‰ of the ∼2.9‰ difference between mid-Holocene and LGM. The remaining 1.8‰ in the difference can be accounted for by +9 °C of warming between the LGM and mid-Holocene and -3 °C cooling at H events, which are comparable with the previous estimates of land paleo-temperature in the Japanese Islands. The attenuated isotopic signal associated with D-O interstadials indicates that the warming in the Atlantic did not significantly transfer to the maritime Japan. These unique features of the isotopic records of the Japanese stalagmites are due to the geographic position at the vicinity of the moisture source, Kuroshio warm current.

Random and externally controlled occurrences of Dansgaard–Oeschger events

Abstract. Dansgaard–Oeschger (DO) events constitute the most pronounced mode of centennial to millennial climate variability of the last glacial period. Since their discovery, many decades of research have been devoted to understand the origin and nature of these rapid climate shifts. In recent years, a number of studies have appeared that report emergence of DO-type variability in fully coupled general circulation models via different mechanisms. These mechanisms result in the occurrence of DO events at varying degrees of regularity, ranging from periodic to random. When examining the full sequence of DO events as captured in the North Greenland Ice Core Project (NGRIP) ice core record, one can observe high irregularity in the timing of individual events at any stage within the last glacial period. In addition to the prevailing irregularity, certain properties of the DO event sequence, such as the average event frequency or the relative distribution of cold versus warm periods, appear to be changing throughout the glacial. By using statistical hypothesis tests on simple event models, we investigate whether the observed event sequence may have been generated by stationary random processes or rather was strongly modulated by external factors. We find that the sequence of DO warming events is consistent with a stationary random process, whereas dividing the event sequence into warming and cooling events leads to inconsistency with two independent event processes. As we include external forcing, we find a particularly good fit to the observed DO sequence in a model where the average residence time in warm periods are controlled by global ice volume and cold periods by boreal summer insolation.

Climate dynamics during the penultimate glacial period recorded in a speleothem from Kanaan Cave, Lebanon (central Levant)

AbstractLittle is known about terrestrial climate dynamics in the Levant during the penultimate interglacial-glacial period. To decipher the palaeoclimatic history of the Marine Oxygen Isotope Stage (MIS) 6 glacial period, a well-dated stalagmite (~194 to ~154 ka) from Kanaan Cave on the Mediterranean coast in Lebanon was analyzed for its petrography, growth history, and stable isotope geochemistry. A resolved climate record has been recovered from this precisely U–Th dated speleothem, spanning the late MIS 7 and early MIS 6 at low resolution and the mid–MIS 6 at higher resolution. The stalagmite grew discontinuously from ~194 to ~163 ka. More consistent growth and higher growth rates between ~163 and ~154 ka are most probably linked to increased water recharge and thus more humid conditions. More distinct layering in the upper part of the speleothem suggests strong seasonality from ~163 ka to ~154 ka. Short-term oxygen and carbon isotope excursions were found between ~155 and ~163 ka. The inferred Kanaan Cave humid intervals during the mid–MIS 6 follow variations of pollen records in the Mediterranean basins and correlate well with the synthetic Greenland record and East Asian summer monsoon interstadial periods, indicating short warm/wet periods similar to the Dansgaard-Oeschger events during MIS 4–3 in the eastern Mediterranean region.

The Antarctic Centennial Oscillation: A Natural Paleoclimate Cycle in the Southern Hemisphere That Influences Global Temperature

We report a previously-unexplored natural temperature cycle recorded in ice cores from Antarctica—the Antarctic Centennial Oscillation (ACO)—that has oscillated for at least the last 226 millennia. Here we document the properties of the ACO and provide an initial assessment of its role in global climate. We analyzed open-source databases of stable isotopes of oxygen and hydrogen as proxies for paleo-temperatures. We find that centennial-scale spectral peaks from temperature-proxy records at Vostok over the last 10,000 years occur at the same frequencies (±2.4%) in three other paleoclimate records from drill sites distributed widely across the East Antarctic Plateau (EAP), and >98% of individual ACOs evaluated at Vostok match 1:1 with homologous cycles at the other three EAP drill sites and conversely. Identified ACOs summate with millennial periodicity to form the Antarctic Isotope Maxima (AIMs) known to precede Dansgaard-Oeschger (D-O) oscillations recorded in Greenland ice cores. Homologous ACOs recorded at the four EAP drill sites during the last glacial maximum appeared first at lower elevations nearest the ocean and centuries later on the high EAP, with latencies that exceed dating uncertainty >30-fold. ACO homologs at different drill sites became synchronous, however, during the warmer Holocene. Comparative spectral analysis suggests that the millennial-scale AIM cycle declined in period from 1500 to 800 years over the last 70 millennia. Similarly, over the last 226 millennia ACO repetition period (mean 352 years) declined by half while amplitude (mean 0.67 °C) approximately doubled. The period and amplitude of ACOs oscillate in phase with glacial cycles and related surface insolation associated with planetary orbital forces. We conclude that the ACO: encompasses at least the EAP; is the proximate source of D-O oscillations in the Northern Hemisphere; therefore affects global temperature; propagates with increased velocity as temperature increases; doubled in intensity over geologic time; is modulated by global temperature variations associated with planetary orbital cycles; and is the probable paleoclimate precursor of the contemporary Antarctic Oscillation (AAO). Properties of the ACO/AAO are capable of explaining the current global warming signal.

A high-resolution paleosecular variation record from Black Sea sediments indicating fast directional changes associated with low field intensities during marine isotope stage (MIS) 4

A total of nine sediment cores recovered from the Archangelsky Ridge in the SE Black Sea were systematically subjected to intense paleo- and mineral magnetic analyses. Besides 16 accelerator mass spectrometry (AMS) 14C ages available for another core from this area, dating was accomplished by correlation of short-term warming events during the last glacial monitored by high-resolution X-ray fluorescence (XRF) scanning as maxima in both Ca/Ti and K/Ti ratios in Black Sea sediments to the so-called ‘Dansgaard–Oeschger events’ recognized from Greenland ice cores. Thus, several hiatuses could be identified in the various cores during the last glacial/interglacial cycle. Finally, core sections documenting marine isotope stage (MIS) 4 at high resolution back to 69 ka were selected for detailed analyses. At 64.5 ka, according to obtained results from Black Sea sediments, the second deepest minimum in relative paleointensity during the past 69 ka occurred, with the Laschamp geomagnetic excursion at 41 ka being associated with the lowest field intensities. The field minimum during MIS 4 is associated with large declination swings beginning about 3 ka before the minimum. While a swing to 50°E is associated with steep inclinations (50–60°) according to the coring site at 42°N, the subsequent declination swing to 30°W is associated with shallow inclinations of down to 40°. Nevertheless, these large deviations from the direction of a geocentric axial dipole field (I=61°, D=0°) still can not yet be termed as ‘excursional’, since latitudes of corresponding virtual geomagnetic poles (VGP) only reach down to 51.5°N (120°E) and 61.5°N (75°W), respectively. However, these VGP positions at opposite sides of the globe are linked with VGP drift rates of up to 0.2° per year in between. These extreme secular variations might be the mid-latitude expression of a geomagnetic excursion with partly reversed inclinations found at several sites much further North in Arctic marine sediments between 69°N and 81°N. Thus, the pronounced intensity minimum at 64.5 ka and described directional variations might be the effect of a weak geomagnetic field with a multi-polar geometry in the middle of MIS 4.

Analytical constraints on layered gas trapping and smoothing of atmospheric variability in ice under low-accumulation conditions

Abstract. We investigate for the first time the loss and alteration of past atmospheric information from air trapping mechanisms under low-accumulation conditions through continuous CH4 (and CO) measurements. Methane concentration changes were measured over the Dansgaard–Oeschger event 17 (DO-17, ∼ 60 000 yr BP) in the Antarctic Vostok 4G-2 ice core. Measurements were performed using continuous-flow analysis combined with laser spectroscopy. The results highlight many anomalous layers at the centimeter scale that are unevenly distributed along the ice core. The anomalous methane mixing ratios differ from those in the immediate surrounding layers by up to 50 ppbv. This phenomenon can be theoretically reproduced by a simple layered trapping model, creating very localized gas age scale inversions. We propose a method for cleaning the record of anomalous values that aims at minimizing the bias in the overall signal. Once the layered-trapping-induced anomalies are removed from the record, DO-17 appears to be smoother than its equivalent record from the high-accumulation WAIS Divide ice core. This is expected due to the slower sinking and densification speeds of firn layers at lower accumulation. However, the degree of smoothing appears surprisingly similar between modern and DO-17 conditions at Vostok. This suggests that glacial records of trace gases from low-accumulation sites in the East Antarctic plateau can provide a better time resolution of past atmospheric composition changes than previously expected. We also developed a numerical method to extract the gas age distributions in ice layers after the removal of the anomalous layers based on comparison with a weakly smoothed record. It is particularly adapted for the conditions of the East Antarctic plateau, as it helps to characterize smoothing for a large range of very low-temperature and low-accumulation conditions.

Abrupt climatic events recorded by the Ili loess during the last glaciation in Central Asia: Evidence from grain-size and minerals

The loess record of Central Asia provides an important archive of regional climate and environmental changes. In contrast to the widely investigated loess deposits in the Chinese Loess Plateau, Central Asian loess–paleosol sequences remain poorly understood. Here, we present an aeolian loess section in the southern Ili Basin. Based on granularity and mineralogical analyses, we reconstruct climatic changes during the last glaciation. The results indicated that most of the abrupt climatic events (such as Dansgaard-Oeschger events and Heinrich events) were imprinted in this loess section, although their amplitudes and ages showed some differences. Compared with the millennial oscillations recoded in loess and stalagmites in East Asia, the arid Central Asia responded more sensitively to the warming events than to the cooling events. The shifting trajectory of westerlies across Central Asia played an important role in dust deposition during the stadials. The North Atlantic climatic signals may have been transmitted from Central Asia to the East Asian monsoon regions via the westerlies.

Contrasting shell growth strategies in two Mediterranean bivalves revealed by oxygen-isotope ratio geochemistry: The case of Pecten jacobaeus and Glycymeris pilosa

High-resolution stable-isotope ratio data (δ18O, δ13C) were used to study growth strategies of two bivalve species, Pecten jacobaeus (calcitic shell) and Glycymeris pilosa (aragonitic shell) from the North Adriatic Sea. The principal objectives of this study were to identify the period of the year when the growth line is formed in the shell of two target species, to identify the main growing season of these two species, to identify the environmental drivers of shell growth, and to evaluate the potential applicability of δ18O and δ13C values for the reconstruction of environmental variability. Samples were collected from the North Adriatic Sea by commercial bean trawl (P. jacobaeus, December 2013 and January 2014, N=4) and SCUBA diver (Glycymeris pilosa, March 2016, N=3). Samples for the oxygen (δ18O) and carbon (δ13C) isotope composition of the calcium carbonate were collected by drilling the outer shell layer across several annual cycles. Temporal and spatial temperature and salinity values inside the investigated area were simulated using the 3D numerical ocean model - ROMS. The δ18O cycles corresponded to the number of seasonal growth marks observed on the external shell surface of both target species, thereby confirming the annual periodicity of these growth patterns. In February 2012, extreme cooling of the water column accompanied by dense water formation occurred in the Adriatic Sea - an event recorded by P. jacobaeus shells. This study indicates that P. jacobaeus and G. pilosa have contrasting shell growth strategies. Pecten jacobaeus grows during winter and slows shell growth during the warmest part of the year, and thereby may be an interesting archive for winter conditions. Due to its longevity and continuous growth during the warmest part of the year, G. pilosa is a promising archive for the reconstruction of summer seawater temperatures.

The southern Norwegian Sea during the last 45 ka: hydrographical reorganizations under changing ice‐sheet dynamics

ABSTRACTThe last glacial millennial Heinrich and Dansgaard–Oeschger events have been the subject of many studies in recent decades. However, significant uncertainties remain, notably regarding associated oceanic mechanisms. Here we investigate the evolution of hydrographical conditions in the 3–41 ka cal BP interval of core MD99‐2285 located in the southern Norwegian Sea. Our results reveal an unusual scheme with: (i) warm Greenland interstadials marked by a relatively cold and homogeneous upper water column seasonally ice‐covered and active deep water formation, (ii) cold Greenland and Heinrich stadials marked by a warm, low‐salinity and nearly sea ice‐free surface layer with iceberg calving, separated by a strong halocline from a colder and saltier subsurface layer, as well as strongly reduced deep convection, and (iii) close to modern Holocene oceanic conditions, but showing significant differences from interstadial conditions. Hence, our work strongly suggests different deep‐water convection modes in the Nordic Seas, with, notably, interstadial deep water formation mainly forced by brine release. According to our results, the transition from the glacial mode of deep water formation to the modern open ocean convection mode might have occurred at the end of deglaciation/beginning of the Holocene, when both modes probably occurred.

Synchronous precipitation reduction in the American Tropics associated with Heinrich 2

During the last ice age temperature in the North Atlantic oscillated in cycles known as Dansgaard-Oeschger (D-O) events. The magnitude of Caribbean hydroclimate change associated with D-O variability and particularly with stadial intervals, remains poorly constrained by paleoclimate records. We present a 3.3 thousand-year long stalagmite δ18O record from the Yucatan Peninsula (YP) that spans the interval between 26.5 and 23.2 thousand years before present. We estimate quantitative precipitation variability and the high resolution and dating accuracy of this record allow us to investigate how rainfall in the region responds to D-O events. Quantitative precipitation estimates are based on observed regional amount effect variability, last glacial paleotemperature records, and estimates of the last glacial oxygen isotopic composition of precipitation based on global circulation models (GCMs). The new precipitation record suggests significant low latitude hydrological responses to internal modes of climate variability and supports a role of Caribbean hydroclimate in helping Atlantic Meridional Overturning Circulation recovery during D-O events. Significant in-phase precipitation reduction across the equator in the tropical Americas associated with Heinrich event 2 is suggested by available speleothem oxygen isotope records.

A spatio-temporal reconstruction of sea-surface temperatures in the North Atlantic during Dansgaard-Oeschger events 5–8

Here we establish a spatio-temporal evolution of the sea-surface temperatures in the North Atlantic over Dansgaard Oeschger (DO) events 5–8 (c. 30–40 ka) using the proxy surrogate reconstruction method. Proxy data suggest a large variability in North Atlantic sea-surface temperatures during the DO-events of the last glacial period. However, proxy data availability is limited and cannot provide a full spatial picture of the oceanic changes. Therefore, we combine fully coupled, general circulation model simulations with planktic foraminifera based sea-surface temperature reconstructions to obtain a broader spatial picture of the ocean state during DO-events 5–8. The resulting spatial sea-surface temperature patterns agree over a number of different general circulation models and simulations. We find that sea-surface temperature variability over the DO-events is characterized by colder conditions in the subpolar North Atlantic during stadials than during interstadials, and the variability is linked to changes in the Atlantic Meridional Overturning circulation, and in the sea-ice cover. Forced simulations are needed to capture the strength of the temperature variability and to reconstruct the variability in other climatic records not directly linked to the sea-surface temperature reconstructions. Our results are robust to uncertainties in the age models of the proxy data, the number of available temperature reconstructions, and over a range of climate models.

Glacial to Holocene changes in sea surface temperature and seawater δ18O in the northern Indian Ocean

Sea surface temperature (SST) and seawater δ18Osw records were generated from sediment cores located in the southern Bay of Bengal (SBOB) and the northeastern Arabian Sea (NEAS) to understand glacial to Holocene changes in the hydrography of these regions. This was accomplished through the use of paired δ18O and Mg/Ca measurements in planktic foraminifera (Globigerinoides ruber) from two sediment cores; SK157-14 in the SBOB (lat. 5°11′N; long. 90°05′E; water-depth 3306m) and SK148-21 in the NEAS (lat. 21°29′N, long. 67°01′E; water-depth 1900m). The results suggest significant changes in SST and δ18Osw since last glacial period. The glacial SSTs in the SBOB and the NEAS were lower by ~2–2.5°C relative to the Holocene. The deglacial period in both cores is characterized by significant variations in SSTs and seawater δ18Osw. The Dansgaard-Oeschger (DO) cycles and Heinrich events are strongly expressed in SST and δ18Osw records of the SBOB core. These records support a strong control of atmospheric-oceanic changes in the northern high latitude on thermal state of the SBOB. The SST time series in both these regions indicates a sudden and abrupt increase in temperature at the end of last glacial maximum. The SST time series in the NEAS core SK 148-21 reveals an intensification of the northeast monsoon during the last glacial maximum.

Regional seesaw between the North Atlantic and Nordic Seas during the last glacial abrupt climate events

Abstract. Dansgaard–Oeschger oscillations constitute one of the most enigmatic features of the last glacial cycle. Their cold atmospheric phases have been commonly associated with cold sea-surface temperatures and expansion of sea ice in the North Atlantic and adjacent seas. Here, based on dinocyst analyses from the 48–30 ka interval of four sediment cores from the northern Northeast Atlantic and southern Norwegian Sea, we provide direct and quantitative evidence of a regional paradoxical seesaw pattern: cold Greenland and North Atlantic phases coincide with warmer sea-surface conditions and shorter seasonal sea-ice cover durations in the Norwegian Sea as compared to warm phases. Combined with additional palaeorecords and multi-model hosing simulations, our results suggest that during cold Greenland phases, reduced Atlantic meridional overturning circulation and cold North Atlantic sea-surface conditions were accompanied by the subsurface propagation of warm Atlantic waters that re-emerged in the Nordic Seas and provided moisture towards Greenland summit.

Characteristics of the Heinrich 1 abrupt climate event inferred from a speleothem record from Xianyun Cave, Fujian Province

The climate system of the Earth experienced a series of abrupt changes such as Heinrich events and Dansgaard-Oeschger cycles during the last glaciation. Among them, the Heinrich 1 event (H1 event) is the most prominent cold event. It was recorded in Greenland ice cores, deep sea sediments, loess, stalagmites and lacustrine sediments and was proposed to be a global climate event. During H1 event (~18–15 ka), the ice raft extended southward and large amount of cold and fresh water was integrated into the North Atlantic Ocean. It led to the reduction of the Atlantic Meridional Overturning Circulation (AMOC), cooling of the northern hemisphere, and the weakening of the East Asian summer monsoon (EASM) as the results of the southward shift of the Inter-tropical Convergence Zone (ITCZ). Pausata et al. (2011) used a numerical climate model with an embedded oxygen-isotope model to assess what caused the shifts in the oxygen-isotope signature of precipitation during a climate perturbation designed to mimic the Heinrich event. The simulations reproduced the enrichment of δ 18O over Northern India and East Asia evident in speleothem records during Heinrich events. They then concluded that changes in the δ 18O of Chinese stalagmites associated with Heinrich events reflect changes in the intensity of Indian rather than East Asian monsoon precipitations. The simulation result also shows no significant change in δ 18O in precipitation over southeastern China, where its climate is influenced by EASM only. However, there is a lack of well-dated stalagmite record covering the H1 event from the southeastern China to test this hypothesis. Here we present a high resolution stalagmite isotope record (XYⅣ-3) from Xianyun Cave in western Fujian Province to reveal the detailed precipitation history of the EASM from 16.8 to 15.2 ka BP, based on eight high-precision 230Th ages and 256 oxygen isotope data, yielding an average resolution of 7 a. This high resolution stalagmite δ 18O record constrained detailed changes of the EASM during H1 event. It is characterized by a two-step positive shift process beginning at ~16.2 ka BP with the amplitude of 1.8‰, which indicates a significant weakened monsoon period. This result does not support the model simulation which was proposed by Pausata et al. (2011) that the δ 18O of Chinese stalagmites reflected changes of the Indian Summer Monsoon rainfall. Indeed our record suggests that changes in δ 18O of Chinese stalagmites reflect changes of EASM, especially during Henrich events. Comparing with stalagmite records from Hulu Cave and Qingtian Cave, stalagmites records from all 3 caves captured the H1 abrupt climate changes within dating errors. However, changes in Xianyun Cave are more prolonged (~260 a) than the abrupt (18–19 a) changes in Qingtian and Hulu caves. We conclude that the prolonged process of the Xianyun Cave is influenced, to some extent, by ocean-atmosphere coupling processes over the tropical western Pacific Ocean. The hydroclimate of the tropical west Pacific may cause the different responses. Due to the influence of sea-surface temperature changes in the tropical western Pacific warm pool and ENSO, the response of δ 18O in precipitation in Fujian Province evident by our stalagmites to North Atlantic Ocean change would last longer than other EASM regions.

Abrupt North Atlantic circulation changes in response to gradual CO2 forcing in a glacial climate state

Glacial climate is marked by abrupt, millennial-scale climate changes known as Dansgaard–Oeschger cycles. The most pronounced stadial coolings, Heinrich events, are associated with massive iceberg discharges to the North Atlantic. These events have been linked to variations in the strength of the Atlantic meridional overturning circulation. However, the factors that lead to abrupt transitions between strong and weak circulation regimes remain unclear. Here we show that, in a fully coupled atmosphere–ocean model, gradual changes in atmospheric CO2 concentrations can trigger abrupt climate changes, associated with a regime of bi-stability of the Atlantic meridional overturning circulation under intermediate glacial conditions. We find that changes in atmospheric CO2 concentrations alter the transport of atmospheric moisture across Central America, which modulates the freshwater budget of the North Atlantic and hence deep-water formation. In our simulations, a change in atmospheric CO2 levels of about 15 ppmv—comparable to variations during Dansgaard–Oeschger cycles containing Heinrich events—is sufficient to cause transitions between a weak stadial and a strong interstadial circulation mode. Because changes in the Atlantic meridional overturning circulation are thought to alter atmospheric CO2 levels, we infer that atmospheric CO2 may serve as a negative feedback to transitions between strong and weak circulation modes.