Rapid northern hemisphere ice sheet melting during the penultimate deglaciation

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

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

Past 200 kyr hydroclimate variability in the western Mediterranean and its connection to the African Humid Periods

The Iberian Peninsula is located at the intersection between the subtropical and temperate climate zones and the paleoclimate records from this region are key to elucidate the varying humidity and changing dominance of atmospheric circulation patterns in the Mediterranean-North African region in the past. Here we present a quantitative hydroclimate reconstruction for the last ca. 200 kyr from southern Iberian Peninsula based on pollen data from the Padul lake sediment record. We use the newly developed Scale-normalized Significant Zero crossing (SnSiZer) method to detect not only the statistically significant precipitation changes but also to estimate the relative magnitude of these oscillations in our reconstruction. We identify six statistically significant main humid phases, termed West Mediterranean Humid Periods (WMHP 1–6). These humid periods correlate with other West/Central Mediterranean paleohydrological records, suggesting that similar climatic factors affected different areas of the Mediterranean. In addition, the WMPHs are roughly coeval with the African Humid Periods (AHPs) during high seasonality, suggesting the same North Atlantic ocean-atmospheric dynamics and orbital forcing as main drivers of both areas. In contrast, during low seasonality periods, the West Mediterranean still appears to be affected by the westerlies and the local Mediterranean rainfall systems with moderate-to-high precipitation, whereas West Africa was characterized by droughts.

The Potential of Tufa as a Tool for Paleoenvironmental Research—A Study of Tufa from the Zrmanja River Canyon, Croatia

Tufa is a fresh-water surface calcium carbonate deposit precipitated at or near ambient temperature, and commonly contains the remains of macro- and microphytes. Many Holocene tufas are found along the Zrmanja River, Dalmatian karst, Croatia. In this work we present radiocarbon dating results of older tufa that was found for the first time at the Zrmanja River near the Village of Sanaderi. Tufa outcrops were observed at different levels, between the river bed and up to 26 m above its present level. Radiocarbon dating of the carbonate fraction revealed ages from modern, at the river bed, up to 40 kBP ~20 m above its present level. These ages fit well with the hypothesis that the Zrmanja River had a previous surface connection with the Krka River, and changed its flow direction toward the Novigrad Sea approximately 40 kBP (Marine Isotope Stage 3). Radiocarbon AMS dating of tufa organic residue yielded a maximum conventional age of 17 kBP for the highest outcrop position indicating probable penetration of younger organic material to hollow tufa structures, as confirmed by radiocarbon analyses of humin extracted from the samples. Stable carbon isotope composition (δ13C) of the carbonate fraction of (−10.4 ± 0.6)‰ and (−9.7 ± 0.8)‰ for the Holocene and the older samples, respectively, indicate the autochthonous origin of the carbonate. The δ13C values of (−30.5 ± 0.3)‰ and (−29.6 ± 0.6)‰ for organic residue, having ages <500 BP and >5000 BP, respectively, suggest a unique carbon source for photosynthesis, mainly atmospheric CO2, with an indication of the Suess effect in δ13C during last centuries. The oxygen isotopic composition (δ18O) agrees well with deposition of tufa samples in two stages, the Holocene (−8.02 ± 0.72‰) and “old” (mainly MIS 3 and the beginning of MIS 2) (−6.89 ± 0.34‰), suggesting a ~4 °C lower temperature in MIS 3 compared to the current one.

Speleothem Records of the Hydroclimate Variability throughout the Last Glacial Cycle from Manita peć Cave (Velebit Mountain, Croatia)

We present stable carbon (δ13C) and oxygen (δ18O) isotope records from two partially coeval speleothems from Manita peć Cave, Croatia. The cave is located close to the Adriatic coast (3.7 km) at an elevation of 570 m a.s.l. The site experienced competing Mediterranean and continental climate influences throughout the last glacial cycle and was situated close to the ice limit during the glacial phases. U-Th dating constrains the growth history from Marine Isotope Stage (MIS) 5 to MIS 3 and the transition from MIS 2 to MIS 1. 14C dating was used to estimate the age of the youngest part of one stalagmite found to be rich in detrital thorium and thus undatable by U-Th. On a millennial scale, δ18O variations partly mimic the Dansgaard–Oeschger interstadials recorded in Greenland ice cores (Greenland Interstadials, GI) from GI 22 to GI 13. We interpret our δ18O record as a proxy for variations in precipitation amount and/or moisture sources, and the δ13C record is interpreted as a proxy for changes in soil bioproductivity. The latter indicates a generally reduced vegetation cover towards MIS 3–MIS 4, with shifts of ~8‰ and approaching values close to those of the host rock. However, even during the coldest phases, when a periglacial setting and enhanced aridity sustained long-residence-time groundwater, carbonic-acid dissolution remains the driving force of the karstification processes. Speleothem morphology follows changes in environmental conditions and complements regional results of submerged speleothems findings. Specifically, narrow sections of light porous spelaean calcite precipitated during the glacial/stadial sea-level lowstands, while the warmer and wetter conditions were marked with compact calcite and hiatuses in submerged speleothems due to sea-level highstands. Presumably, the transformation of this littoral site to a continental one with somewhat higher amounts of orographic precipitation was a site-specific effect that masked regional environmental changes.

Magnesium in subaqueous speleothems as a potential palaeotemperature proxy

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

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

Millennial climate oscillations controlled the structure and evolution of Termination II

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

Stable Isotope Hydrology of Cave Groundwater and Its Relevance for Speleothem-Based Paleoenvironmental Reconstruction in Croatia

Speleothems deposited from cave drip waters retain, in their calcite lattice, isotopic records of past environmental changes. Among other proxies, δ18O is recognized as very useful for this purpose, but its accurate interpretation depends on understanding the relationship between precipitation and drip water δ18O, a relationship controlled by climatic settings. We analyzed water isotope data of 17 caves from different latitudes and altitudes in relatively small but diverse Croatian karst regions in order to distinguish the dominant influences. Drip water δ18O in colder caves generally shows a greater resemblance to the amount-weighted mean of precipitation δ18O compared to warmer sites, where evaporation plays an important role. However, during glacial periods, today’s ‘warm’ sites were cold, changing the cave characteristics and precipitation δ18O transmission patterns. Superimposed on these settings, each cave has site-specific features, such as morphology (descending or ascending passages), altitude and infiltration elevation, (micro) location (rain shadow or seaward orientation), aquifer architecture (responsible for the drip water homogenization) and cave atmosphere (governing equilibrium or kinetic fractionation). This necessitates an individual approach and thorough monitoring for best comprehension.

Speleothem record attests to stable environmental conditions during Neanderthal-modern human turnover in southern Italy

The causes of Neanderthal-modern human (MH) turnover are ambiguous. While potential biocultural interactions between the two groups are still little known, it is clear that Neanderthals in southern Europe disappeared about 42 thousand years ago (ka) after cohabitation for ~3,000 years with MH. Among a plethora of hypotheses on Neanderthal extinction, rapid climate changes during the Middle to Upper Palaeolithic transition (MUPT) are regarded as a primary factor. Here we show evidence for stable climatic and environmental conditions during the MUPT in a region (Apulia) where Neanderthals and MH coexisted. We base our findings on a rare glacial stalagmite deposited between ~106 and ~27 ka, providing the first continuous western Mediterranean speleothem palaeoclimate archive for this period. The uninterrupted growth of the stalagmite attests to the constant availability of rainfall and vegetated soils, while its δ¹³C-δ¹⁸O palaeoclimate proxies demonstrate that Apulia was not affected by dramatic climate oscillations during the MUPT. Our results imply that, because climate did not play a key role in the disappearance of Neanderthals in this area, Neanderthal-MH turnover must be approached from a perspective that takes into account climatic and environmental conditions favourable for both species.

Chaos and periodicities in a climatic time series of the Iberian Margin

We analyze the time series of the temperature of the sedimentary core MD01-2443 originating from the Iberian Margin with a duration of 420 kyr. The series has been tested for unit-root and a long term trend is estimated. We identify four significant periodicities together with a low climatic activity every 100 kyr, and these were associated with internal and external forcings. Also, we identify a high-frequency fast component that acts on top of a nonlinear, irreversible slow-changing dynamics. We find the presence of chaos in the climate of the Iberian Margin by means of a neural network asymptotic test on the largest Lyapunov exponent. The analysis suggests that the chaotic dynamics is associated with the fast high-frequency component. We also carry out a statistical analysis of the dimensionality of the attractor. Our results confirm the possibility that periodic behavior and chaos may coexist on different time scales. This could lead to different degrees of predictability in the climate system according to the characteristic time scales and/or phase-space locations.

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

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

Last Interglacial Climate in Northern Sweden—Insights from a Speleothem Record

Continental records with absolute dates of the timing and progression of climatic conditions during the Last Interglacial (LIG) from northern Europe are rare. Speleothems from northern Europe have a large potential as archives for LIG environmental conditions since they were formed in sheltered environments and may be preserved beneath ice sheets. Here, we present δ13C and δ18O values from speleothem Kf-21, from Korallgrottan in Jämtland (northwest Sweden). Kf-21 is dated with five MC-ICPMS U-Th dates with errors smaller than ~1 ka. Kf-21 started forming at ~130.2 ka and the main growth phase with relatively constant growth rates lasted from 127.3 ka to 124.4 ka, after which calcite formation ceased. Both δ13C and δ18O show rapid shifts but also trends, with a range of values within their Holocene counterparts from Korallgrottan. Our results indicate an early onset of the LIG in northern Europe with ice-free conditions at ~130 ka. Higher growth rates combined with more negative δ18O values between ~127.3 and 126.8 ka, interpreted here as warmer and more humid conditions, as well as indications of a millennial-scale cold spell centered at 126.2 ka, resemble findings from speleothem records from other parts of Europe, highlighting that these were regional scale climatic patterns.

800-kyr land temperature variations modulated by vegetation changes on Chinese Loess Plateau

The complicity of long-term land surface temperature (LST) changes has been under investigated and less understood, hindering our understanding of the history and mechanism of terrestrial climate change. Here, we report the longest (800 thousand years) LSTs based on distributions of soil fossil bacterial glycerol dialkyl glycerol tetraethers preserved in well-dated loess-paleosol sequences at the center of the Chinese Loess Plateau. We have found a previously-unrecognized increasing early and prolonged warming pattern toward the northwestern plateau at the onset of the past seven deglaciations, corresponding to the decrease in vegetation coverage, suggesting underlying surface vegetation or lack of has played an important role in regulating LSTs, superimposed on the fundamental global glacial–interglacial changes. Our results support that LSTs in semi-humid and semi-arid regions with little vegetation will be more sensitive to the anticipated global temperature rise, while improving vegetation coverage would reduce LSTs and thus ecological impacts.

Modern observation indicates that vegetation cover could modulate land surface temperatures substantially. Here the authors demonstrate that such vegetation feedbacks could be clearly identified in the Chinese Loess Plateau land surface temperature records during past cool periods when vegetation cover was reduced.

Speleothems from the Middle East: An Example of Water Limited Environments in the SISAL Database

The Middle East (ME) spans the transition between a temperate Mediterranean climate in the Levant to hyper-arid sub-tropical deserts in the southern part of the Arabian Peninsula (AP), with the complex alpine topography in the northeast feeding the Euphrates and Tigris rivers which support life in the Southeastern Fertile Crescent (FC). Climate projections predict severe drying in several parts of the ME in response to global warming, making it important to understand the controls of hydro-climate perturbations in the region. Here we discuss 23 ME speleothem stable oxygen isotope (δ18Occ) records from 16 sites from the SISAL_v1 database (Speleothem Isotope Synthesis and Analysis database), which provide a record of past hydro-climatic variability. Sub-millennial changes in ME δ18Occ values primarily indicate changes in past precipitation amounts the result of the main synoptic pattern in the region, specifically Mediterranean cyclones. This pattern is superimposed on change in vapor source δ18O composition. The coherency (or lack thereof) between regional records is reviewed from Pleistocene to present, covering the Last Glacial Maximum (~22 ka), prominent events during deglaciation, and the transition into the Holocene. The available δ18Occ time-series are investigated by binning and normalizing at 25-year and 200-year time windows over the Holocene. Important climatic oscillations in the Holocene are discussed, such as the 8.2 ka, 4.2 ka and 0.7 ka (the Little Ice Age) Before Present events. Common trends in the normalized anomalies are tested against different climate archives. Finally, recommendations for future speleothem-based research in the region are given along with comments on the utility and completeness of the SISAL database.

East Asian hydroclimate modulated by the position of the westerlies during Termination I

Speleothem oxygen isotope records have revolutionized our understanding of the paleo East Asian monsoon, yet there is fundamental disagreement on what they represent in terms of the hydroclimate changes. We report a multiproxy speleothem record of monsoon evolution during the last deglaciation from the middle Yangtze region, which indicates a wetter central eastern China during North Atlantic cooling episodes, despite the oxygen isotopic record suggesting a weaker monsoon. We show that this apparent contradiction can be resolved if the changes are interpreted as a lengthening of the Meiyu rains and shortened post-Meiyu stage, in accordance with a recent hypothesis. Model simulations support this interpretation and further reveal the role of the westerlies in communicating the North Atlantic influence to the East Asian climate.

Enhanced climate instability in the North Atlantic and southern Europe during the Last Interglacial

Considerable ambiguity remains over the extent and nature of millennial/centennial-scale climate instability during the Last Interglacial (LIG). Here we analyse marine and terrestrial proxies from a deep-sea sediment sequence on the Portuguese Margin and combine results with an intensively dated Italian speleothem record and climate-model experiments. The strongest expression of climate variability occurred during the transitions into and out of the LIG. Our records also document a series of multi-centennial intra-interglacial arid events in southern Europe, coherent with cold water-mass expansions in the North Atlantic. The spatial and temporal fingerprints of these changes indicate a reorganization of ocean surface circulation, consistent with low-intensity disruptions of the Atlantic meridional overturning circulation (AMOC). The amplitude of this LIG variability is greater than that observed in Holocene records. Episodic Greenland ice melt and runoff as a result of excess warmth may have contributed to AMOC weakening and increased climate instability throughout the LIG.

It is important to establish a baseline for natural climate variability under relatively warm conditions. Here we show that the Last Interglacial in the North Atlantic and southern Europe was characterized by enhanced climate instability relative to the pre-industrial Holocene.

Evidence of thermophilisation and elevation-dependent warming during the Last Interglacial in the Italian Alps

Thermophilisation is the response of plants communities in mountainous areas to increasing temperatures, causing an upward migration of warm-adapted (thermophilic) species and consequently, the timberline. This greening, associated with warming, causes enhanced evapotranspiration that leads to intensification of the hydrological cycle, which is recorded by hydroclimate-sensitive archives, such as stalagmites and flowstones formed in caves. Understanding how hydroclimate manifests at high altitudes is important for predicting future water resources of many regions of Europe that rely on glaciers and snow accumulation. Using proxy data from three coeval speleothems (stalagmites and flowstone) from the Italian Alps, we reconstructed both the ecosystem and hydrological setting during the Last Interglacial (LIG); a warm period that may provide an analogue to a near-future climate scenario. Our speleothem proxy data, including calcite fabrics and the stable isotopes of calcite and fluid inclusions, indicate a +4.3 ± 1.6 °C temperature anomaly at ~2000 m a.s.l. for the peak LIG, with respect to present-day values (1961–1990). This anomaly is significantly higher than any low-altitude reconstructions for the LIG in Europe, implying elevation-dependent warming during the LIG. The enhanced warming at high altitudes must be accounted for when considering future climate adaption strategies in sensitive mountainous regions.

Reconciliation of the Devils Hole climate record with orbital forcing

The driving force behind Quaternary glacial-interglacial cycles and much associated climate change is widely considered to be orbital forcing. However, previous versions of the iconic Devils Hole (Nevada) subaqueous calcite record exhibit shifts to interglacial values ~10,000 years before orbitally forced ice age terminations, and interglacial durations ~10,000 years longer than other estimates. Our measurements from Devils Hole 2 replicate virtually all aspects of the past 204,000 years of earlier records, except for the timing during terminations, and they lower the age of the record near Termination II by ~8000 years, removing both ~10,000-year anomalies. The shift to interglacial values now broadly coincides with the rise in boreal summer insolation, the marine termination, and the rise in atmospheric CO2, which is consistent with mechanisms ultimately tied to orbital forcing.

Bipolar seesaw control on last interglacial sea level

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

North Atlantic storm track changes during the Last Glacial Maximum recorded by Alpine speleothems

The European Alps are an effective barrier for meridional moisture transport and are thus uniquely placed to record shifts in the North Atlantic storm track pattern associated with the waxing and waning of Late-Pleistocene Northern Hemisphere ice sheets. The lack of well-dated terrestrial proxy records spanning this time period, however, renders the reconstruction of past atmospheric patterns difficult. Here we present a precisely dated, continuous terrestrial record of meteoric precipitation in Europe between 30 and 14.7 ka. In contrast to present-day conditions, our speleothem data provide strong evidence for preferential advection of moisture from the South across the Alps supporting a southward shift of the storm track during the local Last Glacial Maximum (that is, 26.5–23.5 ka). Moreover, our age control indicates that this circulation pattern preceded the Northern Hemisphere precession maximum by ~3 ka, suggesting that obliquity may have played a considerable role in the Alpine ice aggradation.

Insights into Late-Pleistocene Northern Hemisphere storm track variability are hampered by a lack of well-dated proxy records. Here, the authors present a precisely dated record of meteoric precipitation between 30 and 14.7 ka, and show that obliquity may have played a vital role in Alpine glacier advance.

Luminescence dating of glaciofluvial deposits linked to the penultimate glaciation in the Eastern Alps

During the penultimate glaciation vast areas of the Alps were glaciated, with piedmont glaciers protruding into the foreland. In the easternmost part of the northward draining valleys of the Alps, the glaciers did not reach the foreland, but formed valley glaciers confined by the mountainous terrain. This also applies to the Ybbs valley, where samples for luminescence dating out of glaciofluvial gravel accumulations were taken at three locations along the present day river course. In a highly dynamic depositional environment, such as a glacier-fed river system, incomplete resetting of the luminescence signal is possible, in particular when transport distances are short. In such cases, quartz usually is the preferred mineral over feldspar, especially if dose rates are low and may theoretically allow obtaining quartz ages beyond 150 ka. Because previous research has shown, and as corroborated within this study, quartz from the research area exhibits analytical problems in the high age range. Therefore luminescence properties of coarse grain (100-200 μm) quartz and in addition K-rich feldspar were investigated with the aim to reconstruct the chronology of the glacial processes within the Ybbs catchment area. Issues of incomplete bleaching were pIRIR225 encountered and addressed by comparing quartz OSL, fading corrected K feldspar IR50 and pIRIR225 to identify reliable ages. Depositional ages based on quartz OSL and feldspar pIRIR225 signals revealed deposition of ice marginal kame terraces and glaciofluvial foreland terraces during late to middle MIS 6. In combination with results from previous studies, we could reconstruct the valley evolution during the Riss glaciation. Newly gained luminescence ages of the deglaciation in the easternmost Alps coincide with OSL dated deglaciation events in the Western Alps, indicating that climatic change along the north side of the Alps happened simultaneously.