Antarctica Ice Core Research Articles

Two likely stratospheric volcanic eruptions in the 1450s C.E. found in a bipolar, subannually dated 800 year ice core record

AbstractAn 800 year volcanic record is constructed from high‐resolution chemical analysis of recently obtained West Antarctica and central Greenland ice cores. The high accuracy and precision of the ice core chronologies are a result of dating by annual ice layer counting. Nineteen bipolar volcanic signals in this record represent large, explosive eruptions in the tropics with probable climatic impact. One of the two bipolar volcanic signals dated at 1453 and 1459 is probably left by the eruption of the submarine volcano Kuwae in the tropical Pacific, one of the largest volcanic eruptions in the last millennium. The discovery of the two signals in the 1450s casts doubt on the eruption year of 1452 or 1453 for Kuwae based on previous ice core records. The volcanic sulfate deposition patterns in this bipolar record suggest that the later signal is likely from the Kuwae eruption in 1458, although a firm attribution is not possible. Sulfur isotope composition in the volcanic sulfate in the central Greenland cores indicates that both eruptions in the 1450s injected sulfur gases into the stratosphere with probable impact on the global climate. These results are in agreement with tree ring records showing two short cold episodes during that decade. The bipolar volcanic record supports the hypothesis that unusually active volcanism in the thirteenth century contributed to the onset of the Little Ice Age and another active period in the mid fifteenth century may have helped to sustain the Little Ice Age.

Mineral dust variability in central West Antarctica associated with ozone depletion

Abstract. We present here data of mineral dust variability retrieved from an ice core of the central West Antarctic, spanning the last five decades. Main evidence provided by the geochemical analysis is that northerly air mass incursions to the coring site, tracked by insoluble dust microparticles, have declined over the past 50 yr. This result contrasts with dust records from ice cores reported to the coastal West Antarctic that show increases since mid-20th century. We attribute this difference to regional climatic changes due to the ozone depletion and its implications to westerly winds. We found that the diameters of insoluble microparticles in the central West Antarctica ice core are significantly correlated with cyclone depth (energy) and wind intensity around Antarctica.

Factors controlling the nitrate in the DT-401 ice core in eastern Antarctica

Nitrate, an oxidized product of NO x preserved in the polar ice cores, has often been used to estimate past changes of the atmospheric nitrogen cycle. A 102.65 m ice core drilled at DT-401 (79°01′S, 77°00′E) in 1999 provides an opportunity to investigate the possible influencing factors for the nitrate budget in the eastern Antarctica. We studied the relationship between the δ 18O (representing the temperature), accumulation rate, volcanic deposition and the astronomical factors (such as the solar activities, supernovae, etc.), and the nitrate variation along the whole duration (2680-year) of the ice core. Prominent impacts of the accumulation rate acting on the nitrate flux rather than the concentration were detected. However, no significant correlation was found between the δ 18O and the nitrate deposition characteristics (concentration and the flux variations). Volcanic deposition can significantly affect the deposition of nitrate with a decreasing trend accompanied by the nss-SO4 −2 (volcanic signal) peak values. Impacts of the solar activities on the nitrate deposition can be detected at this site, and three prominent periodicities (16.6, 24.0 and 102.0 yr) were found for the nitrate concentration variations. Six climatic events (Dalton Minimum, Maunder Minimum, Sporer Minimum, Wolf Minimum, Oort Minimum and Medieval Maximum) during the past 1150 years were observed with lower nitrate values for the foregoing five events and higher value for the last one.

Microbial analyses of ancient ice core sections from greenland and antarctica.

Ice deposited in Greenland and Antarctica entraps viable and nonviable microbes, as well as biomolecules, that become temporal atmospheric records. Five sections (estimated to be 500, 10,500, 57,000, 105,000 and 157,000 years before present, ybp) from the GISP2D (Greenland) ice core, three sections (500, 30,000 and 70,000 ybp) from the Byrd ice core, and four sections from the Vostok 5G (Antarctica) ice core (10,500, 57,000, 105,000 and 105,000 ybp) were studied by scanning electron microscopy, cultivation and rRNA gene sequencing. Bacterial and fungal isolates were recovered from 10 of the 12 sections. The highest numbers of isolates were found in ice core sections that were deposited during times of low atmospheric CO2, low global temperatures and low levels of atmospheric dust. Two of the sections (GISP2D at 10,500 and 157,000 ybp) also were examined using metagenomic/metatranscriptomic methods. These results indicated that sequences from microbes common to arid and saline soils were deposited in the ice during a time of low temperature, low atmospheric CO2 and high dust levels. Members of Firmicutes and Cyanobacteria were the most prevalent bacteria, while Rhodotorula species were the most common eukaryotic representatives. Isolates of Bacillus, Rhodotorula, Alternaria and members of the Davidiellaceae were isolated from both Greenland and Antarctica sections of the same age, although the sequences differed between the two polar regions.

Trace elements and Pb isotope records in Dome C (East Antarctica) ice over the past 800,000 years

Trace elements (V, Cr, Mn, Fe, Co, Cu, Zn, As, Rb, Sr, Mo, Cd, Sb, Ba, Tl, Pb, Bi, Th and U) and Pb isotopic compositions from the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core have been determined using inductively coupled plasma sector field mass spectrometry (ICP-SFMS) and thermal ionization mass spectrometry (TIMS), covering the period from ~533 kyr BP to ~800 kyr BP, respectively. Our data have enabled us to extend the previous EDC records of trace elements and Pb isotopes from the Holocene back to the Marine Isotopic Stage 20.2, ~800 kyr BP. We here discuss the EDC records of Ba, Rb, Mo, Sb, Cd, Tl, Bi and Pb isotopes. Crustal elements such as Ba and Rb show well defined variations in concentrations in relation to climatic conditions with lower values during the interglacial periods and much higher values during the coldest periods of the last eight climatic cycles. Volcanogenic Cd, Tl and Bi show a less pronounced relationship between concentrations and climatic conditions. The isotopic signatures of Pb suggest that changes in the provenance of dust reaching the East Antarctic Plateau from Potential Source Areas occurred during the interglacial periods before the MBE. Our data suggest that the main factors affecting deposition fluxes and sources of natural trace elements over Antarctica are most likely linked to a progressive coupling of the climates of Antarctica and lower latitudes over the past 800 kyr.

High-resolution variations in size, number and arrangement of air bubbles in the EPICA DML (Antarctica) ice core

AbstractWe investigated the large-scale (10–1000 m) and small-scale (mm–cm) variations in size, number and arrangement of air bubbles in the EPICA Dronning Maud Land (EDML) (Antarctica) ice core, down to the end of the bubble/hydrate transition (BHT) zone. On the large scale, the bubble number density shows a general correlation with the palaeo-temperature proxy, δ18O, and the dust concentration, which means that in Holocene ice there are fewer bubbles than in glacial ice. Small-scale variations in bubble number and size were identified and compared. Above the BHT zone there exists a strong anticorrelation between bubble number density and mean bubble size. In glacial ice, layers of high number density and small bubble size are linked with layers with high impurity content, identified as cloudy bands. Therefore, we regard impurities as a controlling factor for the formation and distribution of bubbles in glacial ice. The anticorrelation inverts in the middle of the BHT zone. In the lower part of the BHT zone, bubble-free layers exist that are also associated with cloudy bands. The high contrast in bubble number density in glacial ice, induced by the impurities, indicates a much more pronounced layering in glacial firn than in modern firn.

High-resolution interpolar difference of atmospheric methane around the Last Glacial Maximum

Abstract. Reconstructions of past atmospheric methane concentrations are available from ice cores from both Greenland and Antarctica. The difference observed between the two polar methane concentration levels represents a valuable constraint on the geographical location of the methane sources. Here we present new high-resolution methane records from the North Greenland Ice Core Project (NGRIP) and the European Project for Ice Coring in Antarctica (EPICA) Dronning Maud Land (EDML) ice cores covering Termination 1, the Last Glacial Maximum, and parts of the last glacial back to 32 000 years before present. Due to the high resolution of the records, the synchronisation between the ice cores from NGRIP and EDML is considerably improved, and the interpolar concentration difference of methane is determined with unprecedented precision and temporal resolution. Relative to the mean methane concentration, we find a rather stable positive relative interpolar difference throughout the record with its minimum value of 3.7 ± 0.7 % between 21 900–21 200 years before present, which is higher than previously estimated in this interval close to the Last Glacial Maximum. This implies that Northern Hemisphere boreal wetland sources were never completely shut off during the peak glacial, as suggested from previous bipolar methane concentration records. Starting at 21 000 years before present, i.e. several millennia prior to the transition into the Holocene, the relative interpolar difference becomes even more positive and stays at a fairly stable level of 6.5 ± 0.8 % during Termination 1. We thus find that the boreal and tropical methane sources increased by approximately the same factor during Termination 1. We hypothesise that latitudinal shifts in the Intertropical Convergence Zone (ITCZ) and the monsoon system contribute, either by dislocation of the methane source regions or, in case of the ITCZ, also by changing the relative atmospheric volumes of the Northern and Southern Hemispheres, to the subtle variations in the relative interpolar concentration difference on glacial/interglacial as well as on millennial time scales.

Sulphate–climate coupling over the past 300,000 years in inland Antarctica

Sulphate aerosols, particularly micrometre-sized particles of sulphate salt and sulphate-adhered dust, can act as cloud condensation nuclei, leading to increased solar scattering that cools Earth's climate. Evidence for such a coupling may lie in the sulphate record from polar ice cores, but previous analyses of melted ice-core samples have provided only sulphate ion concentrations, which may be due to sulphuric acid. Here we present profiles of sulphate salt and sulphate-adhered dust fluxes over the past 300,000 years from the Dome Fuji ice core in inland Antarctica. Our results show a nearly constant flux of sulphate-adhered dust through glacial and interglacial periods despite the large increases in total dust flux during glacial maxima. The sulphate salt flux, however, correlates inversely with temperature, suggesting a climatic coupling between particulate sulphur and temperature. For example, the total sulphate salt flux during the Last Glacial Maximum averages 5.78 mg m(-2) yr(-1), which is almost twice the Holocene value. Although it is based on a modern analogue with considerable uncertainties when applied to the ice-core record, this analysis indicates that the glacial-to-interglacial decrease in sulphate would lessen the aerosol indirect effects on cloud lifetime and albedo, leading to an Antarctic warming of 0.1 to 5 kelvin.

Ranges of moisture-source temperature estimated from Antarctic ice cores stable isotope records over glacial–interglacial cycles

Abstract. A single isotope ratio (δD or δ18O) of water is widely used as an air-temperature proxy in Antarctic ice core records. These isotope ratios, however, do not solely depend on air-temperature but also on the extent of distillation of heavy isotopes out of atmospheric water vapor from an oceanic moisture source to a precipitation site. The temperature changes at the oceanic moisture source (Δ Tsource) and at the precipitation site (Δ Tsite) can be retrieved by using deuterium-excess (d) data. A new d record from Dome Fuji, Antarctica spanning the past 360 000 yr is presented and compared with records from Vostok and EPICA Dome C ice cores. In previous studies, to retrieve Δ Tsource and Δ Tsite information, different linear regression equations were proposed using theoretical isotope distillation models. A major source of uncertainty lies in the coefficient of regression, βsite which is related to the sensitivity of d to Δ Tsite. We show that different ranges of temperature and selections of isotopic model outputs may increase the value of βsite by more than a factor of two. To explore the impacts of this coefficient on reconstructed temperatures, we apply for the first time the exact same methodology to the isotope records from the three Antarctica ice cores. We show that uncertainties in the βsite coefficient strongly affect (i) the glacial–interglacial magnitude of Δ Tsource; (ii) the imprint of obliquity in Δ Tsource and in the site-source temperature gradient. By contrast, we highlight the robustness of Δ Tsite reconstruction using water isotopes records.

9,400 years of cosmic radiation and solar activity from ice cores and tree rings

Understanding the temporal variation of cosmic radiation and solar activity during the Holocene is essential for studies of the solar-terrestrial relationship. Cosmic-ray produced radionuclides, such as (10)Be and (14)C which are stored in polar ice cores and tree rings, offer the unique opportunity to reconstruct the history of cosmic radiation and solar activity over many millennia. Although records from different archives basically agree, they also show some deviations during certain periods. So far most reconstructions were based on only one single radionuclide record, which makes detection and correction of these deviations impossible. Here we combine different (10)Be ice core records from Greenland and Antarctica with the global (14)C tree ring record using principal component analysis. This approach is only possible due to a new high-resolution (10)Be record from Dronning Maud Land obtained within the European Project for Ice Coring in Antarctica in Antarctica. The new cosmic radiation record enables us to derive total solar irradiance, which is then used as a proxy of solar activity to identify the solar imprint in an Asian climate record. Though generally the agreement between solar forcing and Asian climate is good, there are also periods without any coherence, pointing to other forcings like volcanoes and greenhouse gases and their corresponding feedbacks. The newly derived records have the potential to improve our understanding of the solar dynamics and to quantify the solar influence on climate.

Ice and Bedrock Characteristics Underneath Dome C (Antarctica) From Radio Echo Sounding Data Analysis

The radio echo sounding (RES) system is one of the most widely used active remote sensing techniques for polar ice sheet exploration, including bedrock morphology studies and subglacial lake investigations. Recently, bedrock characterization has been improved through the analysis of radar echo strength. In this paper, the analysis of the RES signal amplitude has been used to collect information about the controversial problem of electromagnetic ice absorption to highlight areas of high reflectivity variation, ascribable to wet ice-bedrock interfaces. A method to distinguish a wet or dry bedrock-ice interface using a model to describe the internal ice absorption is proposed and discussed. Moreover, the comparison between the ice absorption rates from RES measurements and from European Project for Ice Coring in Antarctica Dome C (Antarctica) ice core conductivity data, the signal amplitude contributions of internal ice layers, and different kinds of rock interface is evaluated. Encouraged by the results, the data analysis led to obtaining a bedrock reflectivity variation map of the Dome C area. This map outlined a wide dispersion of wet/dry rock interfaces in the studied area, indicating the possibility of flowing water along both sides of the Concordia Trench.

Modeled methanesulfonic acid (MSA) deposition in Antarctica and its relationship to sea ice

[1] Methanesulfonic acid (MSA) has previously been measured in ice cores in Antarctica as a proxy for sea ice extent and Southern Hemisphere circulation. In a series of chemical transport model (GEOS-Chem) sensitivity experiments, we identify mechanisms that control the MSA concentrations recorded in ice cores. Sea ice is linked to MSA via dimethylsulfide (DMS), which is produced biologically in the surface ocean and known to be particularly concentrated in the sea ice zone. Given existing ocean surface DMS concentration data sets, the model does not demonstrate a strong relationship between sea ice and MSA deposition in Antarctica. The variability of DMS emissions associated with sea ice extent is small (11–30%) due to the small interannual variability of sea ice extent. Wind plays a role in the variability in DMS emissions, but its contribution relative to that of sea ice is strongly dependent on the assumed DMS concentrations in the sea ice zone. Atmospheric sulfur emitted as DMS from the sea ice undergoes net transport northward. Our model runs suggest that DMS emissions from the sea ice zone may account for 26–62% of MSA deposition at the Antarctic coast and 36–95% in inland Antarctica. Though our results are sensitive to model assumptions, it is clear that an improved understanding of both DMS concentrations and emissions from the sea ice zone are required to better assess the impact of sea ice variability on MSA deposition to Antarctica.

Dynamic boundary-monsoon intensity hypothesis: evidence from the deglacial Amazon River discharge record

Glacioeustatic- and temperature-corrected planktonic foraminiferal oxygen isotope (∆δ 18O) records from ODP Site 942 on the Amazon Fan provide a means of monitoring past changes in the outflow of the Amazon River. This study focuses on the last deglaciation and reveals that during this period there were significant variations in the outflow, which implies large changes in moisture availability in the Amazon Basin. Aridity in the Amazon Basin seems to occur between 20.5 ka (calendar) to 17.0 ka and 13.6 ka to 11 ka. The second arid period correlates with the start of the Antarctic Cold Reversal and aridity continues throughout the Younger Dryas period. We find that the large-scale trends in Amazon River outflow are dissimilar to high-latitude variability in either hemisphere. Instead high-resolution variations correlate with the δ 18O difference between Greenland and Antarctica ice core temperature records. This suggests a link between Hemispheric temperature gradients and moisture availability over the Amazon. Based on our results and previously published work we present a new testable ‘dynamic boundary-monsoon intensity hypothesis’, which suggests that tropical moisture is not a simple belt that moves north or south. Rather, the northern and southern boundaries of the South American Summer Monsoon (SASM) are independently dynamic and driven by temperature gradients within their individual hemispheres. The intensity of rainfall within the SASM, however, is driven by precessionally modulated insolation and the resultant convection strength. Combining these two influences produces the dynamic heterogenic changes in the moisture availability observed over tropical South America since the Last Glacial Maximum.

Microstructural evolution in the fine-grained region of the Siple Dome (Antarctica) ice core

AbstractAn in-depth analysis of seven samples from the Siple Dome (Antarctica) ice core, using optical microscopy and electron backscatter diffraction, illustrates rotational recrystallization or polygonization in the fine-grained region of the core between 700 and 800 m. Between 640 and 700 m, the microstructure is characterized by a bimodal grain-size distribution and a broken girdle fabric with evidence of polygonization. From 727 to 770 m, mean grain size decreases and a single-maximum fabric is found, and, by 790 m, mean grain size has again increased and a multiple-maxima fabric manifests itself. We compare grain-size distribution, c- and a-axis orientation, and misorientation between adjacent grains. We found that misorientations between adjacent grains in the 727–770 m region were predominantly low-angle and typically around a common a-axis, suggesting polygonization. This conclusion is supported by radar evidence of a physical disturbance at 757 m, which may be correlated with higher than usual strain in the 700–800 m range. Below 770 m, larger less regular misorientations and textural evidence show that migration recrystallization is the primary recrystallization mechanism.

Deep ice cores: the need for going back in time

After the success of EPICA, the European Project for Ice coring in Antarctica which, at Dome C (East Antarctica) has provided access to climate and environmental records covering the last 800 ka (thousands of years), the ice core community is now engaged in the challenge to obtain older records. Obtaining a 1.5 million year (Ma) record of climate and greenhouse gases is one of the priorities of the International Partnerships in Ice Core Sciences (IPICS) largely motivated by the yet unexplained shift in climate cyclicity around the Mid-Pleistocene transition (MPT) around 1.2 Ma ago. From EPICA results and recently published articles, we further examine how such a 1.5 Ma old ice core will help, and is indeed indispensable, to depict and understand this transition.

Seasonal variability in the input of lead, barium and indium to Law Dome, Antarctica

Lead (Pb) isotopic compositions and concentrations, and barium (Ba) and indium (In) concentrations have been determined at monthly resolution in five Law Dome (coastal Eastern Antarctica) ice core sections dated from ∼1757 AD to ∼1898 AD. ‘Natural’ background Pb concentrations in ∼1757 AD average ∼0.2 pg g −1 and can be attributed to mineral dust and volcanic emissions, with 206Pb/ 207Pb ratios reaching up to 1.266 ± 0.002. From ∼1887 AD to ∼1898 AD, Pb concentrations reached ∼5 pg g −1 and 206Pb/ 207Pb ratios decreased to 1.058 ± 0.001 as a result of additional inputs of Pb from anthropogenic sources. Seasonal variability in the late 1880s has been investigated by decoupling volcanic Pb from the total measured Pb concentrations, revealing spring and autumn maxima, and consistent winter minima, in anthropogenic Pb and mineral dust (Ba) concentrations. We link this variability to the annual cycle in the position and strength of the Antarctic Circumpolar Trough and, the Southern Ocean westerly winds to the north of the trough region. During the autumn and spring seasons, these systems increase in strength, transporting more impurity laden air from the Southern Hemisphere continental regions to Eastern Antarctica and Law Dome. As this Pb is isotopically identical to that emitted from south-eastern Australia (Broken Hill, Port Pirie) this implies a relatively direct air trajectory pathway from southern Australia to Law Dome (Eastern Antarctica).

CO 2 and O 2/N 2 variations in and just below the bubble–clathrate transformation zone of Antarctic ice cores

CO 2 measurements on the EPICA (European Project for Ice Coring in Antarctica) DML ice core in depth levels just below the bubble ice–clathrate ice transformation zone (1230–2240 m depth) were performed. In the youngest part (1230–1600 m), they reveal variations of up to 25 ppmv around the mean atmospheric concentration within centimetres, corresponding to a snow deposition interval of a few years. Similar results are found at corresponding depth regions of the Dome C and the Talos Dome ice cores. Since we can exclude all hitherto known processes altering the concentration of CO 2 in ice cores, we present a hypothesis about spatial fractionation of air components related to episodically increasing clathrate formation followed by diffusion processes from bubbles to clathrates. This hypothesis is supported by optical line-scan observations and by O 2/N 2 measurements at the same depth where strong CO 2 variations are detected. Below the clathrate formation zone, this small-scale fractionation process is slowly smoothed out, most likely by diffusion, regaining the initial mean atmospheric concentration. Although this process compromises the representativeness of a single CO 2 measurement on small ice samples in the clathrate formation zone of an ice core, it does not affect the mean atmospheric CO 2 concentration if CO 2 values are averaged over a sufficiently long depth scale (> 10 cm in case of the EPICA DML ice core).

A 2680 year volcanic record from the DT‐401 East Antarctic ice core

Volcanic signals recorded in the Antarctic and Greenland ice cores can provide useful information on past explosive volcanism and its impact. In this study, we carried out a continuous sulfate analysis of a 102.65 m East Antarctic ice core (DT‐401, dated as 2682 years) and identified 36 extensive volcanic eruption signals using Cole‐Dai's method, which gives an average of 1.4 eruptions per century, consistent with the results from the Plateau Remote (PR‐B) ice core. When the record is divided into three parts, the latest millennium (1999–1000 A.D.), the middle millennium (999–1 A.D.), and the earliest 682 years (0 A.D. to 682 B.C.), it is found that there were more volcanic eruptions that occurred during the latest millennium (19 eruptions) than during the middle millennium (10 eruptions) of the record and that the intensities of the eruptions in the latest millennium are markedly larger than those in the middle one. There were only seven events recorded in the earliest 682 years, but their intensities were greater, and nearly half of the eruptions had a similar intensity to Tambora's (1815 A.D.), which differs from the PR‐B record. It is also found that volcanism and its average accumulation rate were lower during the “Little Ice Age” than during the “Medieval Warm Period.” Comparison of volcanic records between DT‐401 and other Antarctica ice cores (PR‐B, Dome C, DT‐263, and Byrd) show that in the East Antarctica area with its lower accumulation rates, postdepositional effects may play an important role in the deposition of the sulfate.

Constraint of the CO2 rise by new atmospheric carbon isotopic measurements during the last deglaciation

The causes of the ∼80 ppmv increase of atmospheric carbon dioxide (CO2) during the last glacial‐interglacial climatic transition remain debated. We analyzed the parallel evolution of CO2 and its stable carbon isotopic ratio (δ13CO2) in the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core to bring additional constraints. Agreeing well but largely improving the Taylor Dome ice core record of lower resolution, our δ13CO2 record is characterized by a W shape, with two negative δ13CO2 excursions of 0.5‰ during Heinrich 1 and Younger Dryas events, bracketing a positive δ13CO2 peak during the Bølling/Allerød warm period. The comparison with marine records and the outputs of two C cycle box models suggest that changes in Southern Ocean ventilation drove most of the CO2 increase, with additional contributions from marine productivity changes on the initial CO2 rise and δ13CO2 decline and from rapid vegetation buildup during the CO2 plateau of the Bølling/Allerød.

The precession phase of the boreal summer monsoon as viewed from the eastern Mediterranean (ODP Site 968)

The astronomical timescale of the Eastern Mediterranean Plio–Pleistocene builds on tuning of sapropel layers to Northern Hemisphere summer insolation maxima. A 3000-year precession lag has become instrumental in the tuning procedure as radiocarbon dating revealed that the midpoint of the youngest sapropel, S1, in the early Holocene occurred approximately 3000 years after the insolation maximum. The origin of the time lag remains elusive, however, because sapropels are generally linked to maximum African monsoon intensities and transient climate modeling results indicate an in-phase behavior of the African monsoon relative to precession forcing. Here we present new high-resolution records of bulk sediment geochemistry and benthic foraminiferal oxygen isotopes from ODP Site 968 in the Eastern Mediterranean. We show that the 3000-year precession time lag of the sapropel midpoints is consistent with (1) the global marine isotope chronology, (2) maximum (monsoonal) precipitation conditions in the Mediterranean region and China derived from radiometrically dated speleothem records, and (3) maximum atmospheric methane concentrations in Antarctica ice cores. We show that the time lag relates to the occurrence of precession-paced North Atlantic cold events, which systematically delayed the onset of strong boreal summer monsoon intensity. Our findings may also explain a non-stationary behavior of the African monsoon over the past 3 million years due to more frequent and intensive cold events in the Late Pleistocene.