Younger Dryas cooling and the Greenland climate response to CO2

Contextualizing wild cereal harvesting at Middle Palaeolithic Ghar-e Boof in the southern Zagros

A stratigraphic sequence from Ghar-e Boof, a cave site in Iran, covering a period of c. 80,000-30,000 BP and containing more than 20,000 seed and chaff remains, allows a detailed study of the use of annual seed species of Palaeolithic hunter-gatherer groups and its evolution under the influence of changing environmental conditions. Taxonomic changes in the archaeobotanical assemblage and the stable carbon isotope data of pistachio support a considerable change in environmental conditions over the sequence from MIS 5a to MIS 3. The exceptional dominance of wild ancestors of modern crop species, including glume wheat and large-seeded legumes from Middle Palaeolithic layers AH VI (OSL ranges 72-81 ka BP), coincides broadly with the transition from MIS 5a to MIS 4. With the beginning of MIS 4 these taxa are strongly reduced, corresponding with a strong decrease in global CO2 concentrations and in the Δ13C values of Pistacia khinjuk/atlantica from the site. Wild glume wheat completely disappears after Middle Palaeolithic AH Vb and never reappears at the site. We hypothesize that the Middle Palaeolithic niche that allowed the harvesting and consumption of wild cereals and legumes ended with a destabilization of the vegetation in early MIS 4.

Donor-Acceptor Activation of Carbon Dioxide

The activation and functionalization of carbon dioxide entails great interest related to its abundance, low toxicity and associated environmental problems. However, the inertness of CO2 has posed a challenge towards its efficient conversion to added-value products. In this review we discuss one of the strategies that have been widely used to capture and activate carbon dioxide, namely the use of donor-acceptor interactions by partnering a Lewis acidic and a Lewis basic fragment. This type of CO2 activation resembles that found in metalloenzymes, whose outstanding performance in catalytically transforming carbon dioxide encourages further bioinspired research. We have divided this review into three general sections based on the nature of the active sites: metal-free examples (mainly formed by frustrated Lewis pairs), main group-transition metal combinations, and transition metal heterobimetallic complexes. Overall, we discuss one hundred compounds that cooperatively activate carbon dioxide by donor-acceptor interactions, revealing a wide range of structural motifs.

Possible obliquity-forced warmth in southern Asia during the last glacial stage

Orbital-scale global climatic changes during the late Quaternary are dominated by high-latitude influenced ~100,000-year global ice-age cycles and monsoon influenced ~23,000-year low-latitude hydroclimate variations. However, the shortage of highly-resolved land temperature records remains a limiting factor for achieving a comprehensive understanding of long-term low-latitude terrestrial climatic changes. Here, we report paired mean annual air temperature (MAAT) and monsoon intensity proxy records over the past 88,000 years from Lake Tengchongqinghai in southwestern China. While summer monsoon intensity follows the ~23,000-year precession beat found also in previous studies, we identify previously unrecognized warm periods at 88,000-71,000 and 45,000-22,000 years ago, with 2-3 °C amplitudes that are close to our recorded full glacial-interglacial range. Using advanced transient climate simulations and comparing with forcing factors, we find that these warm periods in our MAAT record probably depends on local annual mean insolation, which is controlled by Earth's ~41,000-year obliquity cycles and is anti-phased to annual mean insolation at high latitudes. The coincidence of our identified warm periods and intervals of high-frequent dated archaeological evidence highlights the importance of temperature on anatomically modern humans in Asia during the last glacial stage.

Abrupt Heinrich Stadial 1 cooling missing in Greenland oxygen isotopes

Abrupt climate changes during the last deglaciation have been well preserved in proxy records across the globe. However, one long-standing puzzle is the apparent absence of the onset of the Heinrich Stadial 1 (HS1) cold event around 18 ka in Greenland ice core oxygen isotope δ¹⁸ O records, inconsistent with other proxies. Here, combining proxy records with an isotope-enabled transient deglacial simulation, we propose that a substantial HS1 cooling onset did indeed occur over the Arctic in winter. However, this cooling signal in the depleted oxygen isotopic composition is completely compensated by the enrichment because of the loss of winter precipitation in response to sea ice expansion associated with AMOC slowdown during extreme glacial climate. In contrast, the Arctic summer warmed during HS1 and YD because of increased insolation and greenhouse gases, consistent with snowline reconstructions. Our work suggests that Greenland δ¹⁸ O may substantially underestimate temperature variability during cold glacial conditions.

Vegetation feedback causes delayed ecosystem response to East Asian Summer Monsoon Rainfall during the Holocene

One long-standing issue in the paleoclimate records is whether East Asian Summer Monsoon peaked in the early Holocene or mid-Holocene. Here, combining a set of transient earth system model simulations with proxy records, we propose that, over northern China, monsoon rainfall peaked in the early Holocene, while soil moisture and tree cover peaked in the mid-Holocene. The delayed ecosystem (soil moisture and tree cover) response to rainfall is caused by the vegetation response to winter warming and the subsequent feedback with soil moisture. Our study provides a mechanism for reconciling different evolution behaviors of monsoon proxy records; it sheds light on the driving mechanism of the monsoon evolution and monsoon-ecosystem feedback over northern China, with implications to climate changes in other high climate sensitivity regions over the globe.

StableClim, continuous projections of climate stability from 21000 BP to 2100 CE at multiple spatial scales

Paleoclimatic data are used in eco-evolutionary models to improve knowledge of biogeographical processes that drive patterns of biodiversity through time, opening windows into past climate–biodiversity dynamics. Applying these models to harmonised simulations of past and future climatic change can strengthen forecasts of biodiversity change. StableClim provides continuous estimates of climate stability from 21,000 years ago to 2100 C.E. for ocean and terrestrial realms at spatial scales that include biogeographic regions and climate zones. Climate stability is quantified using annual trends and variabilities in air temperature and precipitation, and associated signal-to-noise ratios. Thresholds of natural variability in trends in regional- and global-mean temperature allow periods in Earth’s history when climatic conditions were warming and cooling rapidly (or slowly) to be identified and climate stability to be estimated locally (grid-cell) during these periods of accelerated change. Model simulations are validated against independent paleoclimate and observational data. Projections of climatic stability, accessed through StableClim, will improve understanding of the roles of climate in shaping past, present-day and future patterns of biodiversity.

Measurement(s)	climate change • climate • temperature of air • volume of hydrological precipitation
Technology Type(s)	computational modeling technique • digital curation
Factor Type(s)	timing of temperature and precipitation estimates
Sample Characteristic - Environment	climate system
Sample Characteristic - Location	Earth (planet)

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12831935

Oceanic forcing of penultimate deglacial and last interglacial sea-level rise

Sea-level histories during the two most recent deglacial-interglacial intervals show substantial differences^1-3 despite both periods undergoing similar changes in global mean temperature^4,5 and forcing from greenhouse gases⁶. Although the last interglaciation (LIG) experienced stronger boreal summer insolation forcing than the present interglaciation⁷, understanding why LIG global mean sea level may have been six to nine metres higher than today has proven particularly challenging². Extensive areas of polar ice sheets were grounded below sea level during both glacial and interglacial periods, with grounding lines and fringing ice shelves extending onto continental shelves⁸. This suggests that oceanic forcing by subsurface warming may also have contributed to ice-sheet loss^9-12 analogous to ongoing changes in the Antarctic^13,14 and Greenland¹⁵ ice sheets. Such forcing would have been especially effective during glacial periods, when the Atlantic Meridional Overturning Circulation (AMOC) experienced large variations on millennial timescales¹⁶, with a reduction of the AMOC causing subsurface warming throughout much of the Atlantic basin^9,12,17. Here we show that greater subsurface warming induced by the longer period of reduced AMOC during the penultimate deglaciation can explain the more-rapid sea-level rise compared with the last deglaciation. This greater forcing also contributed to excess loss from the Greenland and Antarctic ice sheets during the LIG, causing global mean sea level to rise at least four metres above modern levels. When accounting for the combined influences of penultimate and LIG deglaciation on glacial isostatic adjustment, this excess loss of polar ice during the LIG can explain much of the relative sea level recorded by fossil coral reefs and speleothems at intermediate- and far-field sites.

Detecting past changes in vegetation resilience in the context of a changing climate

Anthropogenic climate change is continuously altering ecological responses to disturbance and must be accounted for when examining ecological resilience. One way to measure resilience in ecological datasets is by considering the amount and duration of change from a baseline created by perturbations, such as disturbances like wildfire. Recovery occurs when ecological conditions return to equilibrium, meaning that no subsequent changes can be attributed to the effects of the disturbance, but climate change often causes the recovered state to differ from the previous baseline. The palaeoecological record provides an opportunity to examine these expectations because palaeoclimates changed continuously; few periods existed when environmental conditions were stationary. Here we demonstrate a framework for examining resilience in palaeoecological records against the backdrop of a non-stationary climate by considering resilience as two components of (i) resistance (magnitude of change) and (ii) recovery (time required to return) to predicted equilibrium values. Measuring these components of resilience in palaeoecological records requires high-resolution fossil (e.g. pollen) records, local palaeoclimate reconstructions, a model to predict ecological change in response to climate change, and disturbance records measured at the same spatial scale as the ecological (e.g. vegetation history) record. Resistance following disturbance is measured as the deviation of the fossil record from the ecological state predicted by the palaeoclimate records, and recovery time is measured as the time required for the fossil record to return to predicted values. We show that some cases may involve nearly persistent equilibrium despite large climate changes, but that others can involve a shift to a new state without any complete recovery.

North-south dipole in winter hydroclimate in the western United States during the last deglaciation

During the termination of the last glacial period the western U.S. experienced exceptionally wet conditions, driven by changes in location and strength of the mid-latitude winter storm track. The distribution of modern winter precipitation is frequently characterized by a north-south wet/dry dipole pattern, controlled by interaction of the storm track with ocean-atmosphere conditions over the Pacific and Atlantic Oceans. Here we show that a dipole pattern of similar geographic extent persisted and switched sign during millennial-scale abrupt climate changes of the last deglaciation, based on a new lake level reconstruction for pluvial Lake Chewaucan (northwestern U.S.), and a compilation of regional paleoclimate records. This suggests the dipole pattern is robust, and one mode may be favored for centuries, thereby creating persistent contrasting wet/dry conditions across the western U.S. The TraCE-21k climate model simulation shows an equatorward enhancement of winter storm track activity in the northeastern Pacific, favoring wet conditions in southwestern U.S. during the second half of Heinrich Stadial 1 (16.1–14.6 ka) and consistent with paleoclimate evidence. During the Bølling/Allerød (14.6–12.8 ka), the northeastern Pacific storm track contracted poleward, consistent with wetter conditions concentrated poleward toward the northwest U.S.

Impact of abrupt sea ice loss on Greenland water isotopes during the last glacial period

The Dansgaard–Oeschger events contained in Greenland ice cores constitute the archetypal record of abrupt climate change. An accurate understanding of these events hinges on interpretation of Greenland records of oxygen and nitrogen isotopes. We present here the important results from a suite of modeled Dansgaard–Oeschger events. These simulations show that the change in oxygen isotope per degree of warming becomes smaller during larger events. Abrupt reductions in sea ice also emerge as a strong control on ice core oxygen isotopes because of the influence on both the moisture source and the regional temperature increase. This work confirms the significance of sea ice for past abrupt warming events.

Greenland ice cores provide excellent evidence of past abrupt climate changes. However, there is no universally accepted theory of how and why these Dansgaard–Oeschger (DO) events occur. Several mechanisms have been proposed to explain DO events, including sea ice, ice shelf buildup, ice sheets, atmospheric circulation, and meltwater changes. DO event temperature reconstructions depend on the stable water isotope (δ18O) and nitrogen isotope measurements from Greenland ice cores: interpretation of these measurements holds the key to understanding the nature of DO events. Here, we demonstrate the primary importance of sea ice as a control on Greenland ice core δ18O: 95% of the variability in δ18O in southern Greenland is explained by DO event sea ice changes. Our suite of DO events, simulated using a general circulation model, accurately captures the amplitude of δ18O enrichment during the abrupt DO event onsets. Simulated geographical variability is broadly consistent with available ice core evidence. We find an hitherto unknown sensitivity of the δ18O paleothermometer to the magnitude of DO event temperature increase: the change in δ18O per Kelvin temperature increase reduces with DO event amplitude. We show that this effect is controlled by precipitation seasonality.

Instability of the Northeast Greenland Ice Stream over the last 45,000 years

The sensitivity of the Northeast Greenland Ice Stream (NEGIS) to prolonged warm periods is largely unknown and geological records documenting such long-term changes are needed to place current observations in perspective. Here we use cosmogenic surface exposure and radiocarbon ages to determine the magnitude of NEGIS margin fluctuations over the last 45 kyr (thousand years). We find that the NEGIS experienced slow early Holocene ice-margin retreat of 30–40 m a⁻¹, likely as a result of the buttressing effect of sea-ice or shelf-ice. The NEGIS was ~20–70 km behind its present ice-extent ~41–26 ka and ~7.8–1.2 ka; both periods of high orbital precession index and/or summer temperatures within the projected warming for the end of this century. We show that the NEGIS was smaller than present for approximately half of the last ~45 kyr and is susceptible to subtle changes in climate, which has implications for future stability of this ice stream.

The outlet glaciers that comprise the Northeast Greenland Ice Stream (NEGIS) have experienced accelerated retreat in recent years, yet their longterm stability remains unclear. Here, via cosmogenic surface exposure and radiocarbon ages, the authors investigate the stability of the NEGIS for the past 45 kyr.

Warm summers during the Younger Dryas cold reversal

The Younger Dryas (YD) cold reversal interrupts the warming climate of the deglaciation with global climatic impacts. The sudden cooling is typically linked to an abrupt slowdown of the Atlantic Meridional Overturning Circulation (AMOC) in response to meltwater discharges from ice sheets. However, inconsistencies regarding the YD-response of European summer temperatures have cast doubt whether the concept provides a sufficient explanation. Here we present results from a high-resolution global climate simulation together with a new July temperature compilation based on plant indicator species and show that European summers remain warm during the YD. Our climate simulation provides robust physical evidence that atmospheric blocking of cold westerly winds over Fennoscandia is a key mechanism counteracting the cooling impact of an AMOC-slowdown during summer. Despite the persistence of short warm summers, the YD is dominated by a shift to a continental climate with extreme winter to spring cooling and short growing seasons.

Mechanisms causing the Younger Dryas cold reversal have been questioned by inconsistencies between proxy and modelling results. Here, the authors show that the concept of a strong North Atlantic Ocean cooling event as major driver is consistent with warm European summers caused by intensified atmospheric blocking.

Abrupt Bølling warming and ice saddle collapse contributions to the Meltwater Pulse 1a rapid sea level rise

Elucidating the source(s) of Meltwater Pulse 1a, the largest rapid sea level rise caused by ice melt (14-18 m in less than 340 years, 14,600 years ago), is important for understanding mechanisms of rapid ice melt and the links with abrupt climate change. Here we quantify how much and by what mechanisms the North American ice sheet could have contributed to Meltwater Pulse 1a, by driving an ice sheet model with two transient climate simulations of the last 21,000 years. Ice sheet perturbed physics ensembles were run to account for model uncertainties, constraining ice extent and volume with reconstructions of 21,000 years ago to present. We determine that the North American ice sheet produced 3-4 m global mean sea level rise in 340 years due to the abrupt Bølling warming, but this response is amplified to 5-6 m when it triggers the ice sheet saddle collapse.

Correlation and anti-correlation of the East Asian summer and winter monsoons during the last 21,000 years

Understanding the past significant changes of the East Asia Summer Monsoon (EASM) and Winter Monsoon (EAWM) is critical for improving the projections of future climate over East Asia. One key issue that has remained outstanding from the paleo-climatic records is whether the evolution of the EASM and EAWM are correlated. Here, using a set of long-term transient simulations of the climate evolution of the last 21,000 years, we show that the EASM and EAWM are positively correlated on the orbital timescale in response to the precessional forcing, but are anti-correlated on millennial timescales in response to North Atlantic melt water forcing. The relation between EASM and EAWM can differ dramatically for different timescales because of the different response mechanisms, highlighting the complex dynamics of the East Asian monsoon system and the challenges for future projection.

Future projection of changes in the East Asia Summer and Winter Monsoon are hindered by a lack of understanding of past variability. Here, using longterm transient simulations, the authors show that the monsoons respond in phase to precessional forcing, yet out of phase millennial-scale North Atlantic forcing.

Coherent changes of southeastern equatorial and northern African rainfall during the last deglaciation

During the last deglaciation, wetter conditions developed abruptly ~14,700 years ago in southeastern equatorial and northern Africa and continued into the Holocene. Explaining the abrupt onset and hemispheric coherence of this early African Humid Period is challenging due to opposing seasonal insolation patterns. In this work, we use a transient simulation with a climate model that provides a mechanistic understanding of deglacial tropical African precipitation changes. Our results show that meltwater-induced reduction in the Atlantic meridional overturning circulation (AMOC) during the early deglaciation suppressed precipitation in both regions. Once the AMOC reestablished, wetter conditions developed north of the equator in response to high summer insolation and increasing greenhouse gas (GHG) concentrations, whereas wetter conditions south of the equator were a response primarily to the GHG increase.

Greenland temperature response to climate forcing during the last deglaciation

Greenland ice core water isotopic composition (δ(18)O) provides detailed evidence for abrupt climate changes but is by itself insufficient for quantitative reconstruction of past temperatures and their spatial patterns. We investigate Greenland temperature evolution during the last deglaciation using independent reconstructions from three ice cores and simulations with a coupled ocean-atmosphere climate model. Contrary to the traditional δ(18)O interpretation, the Younger Dryas period was 4.5° ± 2°C warmer than the Oldest Dryas, due to increased carbon dioxide forcing and summer insolation. The magnitude of abrupt temperature changes is larger in central Greenland (9° to 14°C) than in the northwest (5° to 9°C), fingerprinting a North Atlantic origin. Simulated changes in temperature seasonality closely track changes in the Atlantic overturning strength and support the hypothesis that abrupt climate change is mostly a winter phenomenon.