Deglacial rapid sea level rises caused by ice-sheet saddle collapses

Enhanced subglacial discharge from Antarctica during meltwater pulse 1A

Subglacial discharge from the Antarctic Ice Sheet (AIS) likely played a crucial role in the loss of the ice sheet and the subsequent rise in sea level during the last deglaciation. However, no direct proxy is currently available to document subglacial discharge from the AIS, which leaves significant gaps in our understanding of the complex interactions between subglacial discharge and ice-sheet stability. Here we present deep-sea coral 234U/238U records from the Drake Passage in the Southern Ocean to track subglacial discharge from the AIS. Our findings reveal distinctively higher seawater 234U/238U values from 15,400 to 14,000 years ago, corresponding to the period of the highest iceberg-rafted debris flux and the occurrence of the meltwater pulse 1A event. This correlation suggests a causal link between enhanced subglacial discharge, synchronous retreat of the AIS, and the rapid rise in sea levels. The enhanced subglacial discharge and subsequent AIS retreat appear to have been preconditioned by a stronger and warmer Circumpolar Deep Water, thus underscoring the critical role of oceanic heat in driving major ice-sheet retreat.

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.

The timing, duration and magnitude of the 8.2 ka event in global speleothem records

Abrupt events are a feature of many palaeoclimate records during the Holocene. The best example is the 8.2 ka event, which was triggered by a release of meltwater into the Labrador Sea and resulted in a weakening of poleward heat transport in the North Atlantic. We use an objective method to identify rapid climate events in globally distributed speleothem oxygen isotope records during the Holocene. We show that the 8.2 ka event can be identified in >70% of the speleothem records and is the most coherent signal of abrupt climate change during the last 12,000 years. The isotopic changes during the event are regionally homogenous: positive oxygen isotope anomalies are observed across Asia and negative anomalies are seen across Europe, the Mediterranean, South America and southern Africa. The magnitude of the isotopic excursions in Europe and Asia are statistically indistinguishable. There is no significant difference in the duration and timing of the 8.2 ka event between regions, or between the speleothem records and Greenland ice core records. Our study supports a rapid and global climate response to the 8.2 ka freshwater pulse into the North Atlantic, likely transmitted globally via atmospheric teleconnections.

Helicobacter pylori's historical journey through Siberia and the Americas

The peopling of Siberia and the Americas is intriguing for archaeologists, linguists, and human geneticists, but despite significant recent developments, many details remain controversial. Here, we provide insights based on genetic diversity within Helicobacter pylori, a bacterium that infects 50% of all humans. H. pylori strains were collected from across eastern Eurasia and the Americas. Sequence analyses indicated that Siberia contains both anciently diverged and recently admixed bacteria, supporting both human persistence over the last glacial maximum and more recent human recolonization. We inferred a single migration across the Bering land bridge, accompanied by a dramatic reduction in effective population size, followed by bidirectional Holocene gene flow between Asia and the Americas.

The gastric bacterium Helicobacter pylori shares a coevolutionary history with humans that predates the out-of-Africa diaspora, and the geographical specificities of H. pylori populations reflect multiple well-known human migrations. We extensively sampled H. pylori from 16 ethnically diverse human populations across Siberia to help resolve whether ancient northern Eurasian populations persisted at high latitudes through the last glacial maximum and the relationships between present-day Siberians and Native Americans. A total of 556 strains were cultivated and genotyped by multilocus sequence typing, and 54 representative draft genomes were sequenced. The genetic diversity across Eurasia and the Americas was structured into three populations: hpAsia2, hpEastAsia, and hpNorthAsia. hpNorthAsia is closely related to the subpopulation hspIndigenousAmericas from Native Americans. Siberian bacteria were structured into five other subpopulations, two of which evolved through a divergence from hpAsia2 and hpNorthAsia, while three originated though Holocene admixture. The presence of both anciently diverged and recently admixed strains across Siberia support both Pleistocene persistence and Holocene recolonization. We also show that hspIndigenousAmericas is endemic in human populations across northern Eurasia. The evolutionary history of hspIndigenousAmericas was reconstructed using approximate Bayesian computation, which showed that it colonized the New World in a single migration event associated with a severe demographic bottleneck followed by low levels of recent admixture across the Bering Strait.

A reconciled solution of Meltwater Pulse 1A sources using sea-level fingerprinting

The most rapid global sea-level rise event of the last deglaciation, Meltwater Pulse 1A (MWP-1A), occurred ∼14,650 years ago. Considerable uncertainty regarding the sources of meltwater limits understanding of the relationship between MWP-1A and the concurrent fast-changing climate. Here we present a data-driven inversion approach, using a glacio-isostatic adjustment model to invert for the sources of MWP-1A via sea-level constraints from six geographically distributed sites. The results suggest contributions from Antarctica, 1.3 m (0-5.9 m; 95% probability), Scandinavia, 4.6 m (3.2-6.4 m) and North America, 12.0 m (5.6-15.4 m), giving a global mean sea-level rise of 17.9 m (15.7-20.2 m) in 500 years. Only a North American dominant scenario successfully predicts the observed sea-level change across our six sites and an Antarctic dominant scenario is firmly refuted by Scottish isolation basin records. Our sea-level based results therefore reconcile with field-based ice-sheet reconstructions.

Magnitude of the 8.2 ka event freshwater forcing based on stable isotope modelling and comparison to future Greenland melting

The northern hemisphere experienced an abrupt cold event ~ 8200 years ago (the 8.2 ka event) that was triggered by the release of meltwater into the Labrador Sea, and resulting in a weakening of the poleward oceanic heat transport. Although this event has been considered a possible analogue for future ocean circulation changes due to the projected Greenland Ice Sheet (GIS) melting, large uncertainties in the amount and rate of freshwater released during the 8.2 ka event make such a comparison difficult. In this study, we compare sea surface temperatures and oxygen isotope ratios from 28 isotope-enabled model simulations with 35 paleoproxy records to constrain the meltwater released during the 8.2 ka event. Our results suggest that a combination of 5.3 m of meltwater in sea level rise equivalent (SLR) released over a thousand years, with a short intensification over ~ 130 years (an additional 2.2 m of equivalent SLR) due to routing of the Canadian river discharge, best reproduces the proxy anomalies. Our estimate is of the same order of magnitude as projected future GIS melting rates under the high emission scenario RCP8.5.

Sea level fingerprinting of the Bering Strait flooding history detects the source of the Younger Dryas climate event

During the Last Glacial Maximum, expansive continental ice sheets lowered globally averaged sea level ~130 m, exposing a land bridge at the Bering Strait. During the subsequent deglaciation, sea level rose rapidly and ultimately flooded the Bering Strait, linking the Arctic and Pacific Oceans. Observational records of the Bering Strait flooding have suggested two apparently contradictory scenarios for the timing of the reconnection. We reconcile these enigmatic datasets using gravitationally self-consistent sea-level simulations that vary the timing and geometry of ice retreat between the Laurentide and Cordilleran Ice Sheets to the southwest of the Bering Strait to fit observations of a two-phased flooding history. Assuming the datasets are robust, we demonstrate that their reconciliation requires a substantial melting of the Cordilleran and western Laurentide Ice Sheet from 13,000 to 11,500 years ago. This timing provides a freshwater source for the widely debated Younger Dryas cold episode (12,900 to 11,700 years ago).

Labrador Sea freshening at 8.5 ka BP caused by Hudson Bay Ice Saddle collapse

A significant reduction in the Atlantic Meridional Overturning Circulation and rapid northern Hemisphere cooling 8200 years ago have been linked to the final melting of the Laurentide Ice Sheet. Although many studies associated this cold event with the drainage of Lake Agassiz-Ojibway, recent model simulations have shown that the Hudson Bay Ice Saddle collapse would have had much larger effects on the Atlantic Meridional Overturning Circulation than the lake outburst itself. Based on a combination of Mg/Ca and oxygen isotope ratios of benthic foraminifera, this study presents the first direct evidence of a major Labrador shelfwater freshening at 8.5 ka BP, which we associate with the Hudson Bay Ice Saddle collapse. The freshening is preceded by a subsurface warming of the western Labrador Sea, which we link to the strengthening of the West Greenland Current that could concurrently have accelerated the ice saddle collapse in Hudson Bay.

The exact freshwater scenario that caused the 8.2 ka cold event is still debated. This study presents new evidence for a severe Labrador shelfwater freshening 8500 years ago that was caused by the Hudson Bay Ice Saddle collapse following a warming in subsurface waters.

The silicon cycle impacted by past ice sheets

Globally averaged riverine silicon (Si) concentrations and isotope composition (δ³⁰Si) may be affected by the expansion and retreat of large ice sheets during glacial−interglacial cycles. Here we provide evidence of this based on the δ³⁰Si composition of meltwater runoff from a Greenland Ice Sheet catchment. Glacier runoff has the lightest δ³⁰Si measured in running waters (−0.25 ± 0.12‰), significantly lower than nonglacial rivers (1.25 ± 0.68‰), such that the overall decline in glacial runoff since the Last Glacial Maximum (LGM) may explain 0.06–0.17‰ of the observed ocean δ³⁰Si rise (0.5–1.0‰). A marine sediment core proximal to Iceland provides further evidence for transient, low-δ³⁰Si meltwater pulses during glacial termination. Diatom Si uptake during the LGM was likely similar to present day due to an expanded Si inventory, which raises the possibility of a feedback between ice sheet expansion, enhanced Si export to the ocean and reduced CO₂ concentration in the atmosphere, because of the importance of diatoms in the biological carbon pump.

The role ice sheets play in the silica cycle over glacial−interglacial timescales remains unclear. Here, based on the measurement of silica isotopes in Greenland meltwater and a nearby marine sediment core, the authors suggest expanding ice sheets considerably increased isotopically light silica in the oceans.

Glacial isostatic uplift of the European Alps

Following the last glacial maximum (LGM), the demise of continental ice sheets induced crustal rebound in tectonically stable regions of North America and Scandinavia that is still ongoing. Unlike the ice sheets, the Alpine ice cap developed in an orogen where the measured uplift is potentially attributed to tectonic shortening, lithospheric delamination and unloading due to deglaciation and erosion. Here we show that ∼90% of the geodetically measured rock uplift in the Alps can be explained by the Earth’s viscoelastic response to LGM deglaciation. We modelled rock uplift by reconstructing the Alpine ice cap, while accounting for postglacial erosion, sediment deposition and spatial variations in lithospheric rigidity. Clusters of excessive uplift in the Rhône Valley and in the Eastern Alps delineate regions potentially affected by mantle processes, crustal heterogeneity and active tectonics. Our study shows that even small LGM ice caps can dominate present-day rock uplift in tectonically active regions.

For half a century, the cause for recent uplift of the European Alps has been debated. Here, the authors show that ∼90% of the geodetically measured rock uplift in the Alps can be explained by the Earth’s viscoelastic response to ice melting after the Last Glacial Maximum.

Ocean mixing and ice-sheet control of seawater 234U/238U during the last deglaciation

Seawater ²³⁴U/²³⁸U provides global-scale information about continental weathering and is vital for marine uranium-series geochronology. Existing evidence supports an increase in ²³⁴U/²³⁸U since the last glacial period, but the timing and amplitude of its variability has been poorly constrained. Here we report two seawater ²³⁴U/²³⁸U records based on well-preserved deep-sea corals from the low-latitude Atlantic and Pacific Oceans. The Atlantic ²³⁴U/²³⁸U started to increase before major sea-level rise and overshot the modern value by 3 per mil during the early deglaciation. Deglacial ²³⁴U/²³⁸U in the Pacific converged with that in the Atlantic after the abrupt resumption of Atlantic meridional overturning. We suggest that ocean mixing and early deglacial release of excess ²³⁴U from enhanced subglacial melting of the Northern Hemisphere ice sheets have driven the observed ²³⁴U/²³⁸U evolution.

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.

Ice stream activity scaled to ice sheet volume during Laurentide Ice Sheet deglaciation

The contribution of the Greenland and West Antarctic ice sheets to sea level has increased in recent decades, largely owing to the thinning and retreat of outlet glaciers and ice streams. This dynamic loss is a serious concern, with some modelling studies suggesting that the collapse of a major ice sheet could be imminent or potentially underway in West Antarctica, but others predicting a more limited response. A major problem is that observations used to initialize and calibrate models typically span only a few decades, and, at the ice-sheet scale, it is unclear how the entire drainage network of ice streams evolves over longer timescales. This represents one of the largest sources of uncertainty when predicting the contributions of ice sheets to sea-level rise. A key question is whether ice streams might increase and sustain rates of mass loss over centuries or millennia, beyond those expected for a given ocean-climate forcing. Here we reconstruct the activity of 117 ice streams that operated at various times during deglaciation of the Laurentide Ice Sheet (from about 22,000 to 7,000 years ago) and show that as they activated and deactivated in different locations, their overall number decreased, they occupied a progressively smaller percentage of the ice sheet perimeter and their total discharge decreased. The underlying geology and topography clearly influenced ice stream activity, but--at the ice-sheet scale--their drainage network adjusted and was linked to changes in ice sheet volume. It is unclear whether these findings can be directly translated to modern ice sheets. However, contrary to the view that sees ice streams as unstable entities that can accelerate ice-sheet deglaciation, we conclude that ice streams exerted progressively less influence on ice sheet mass balance during the retreat of the Laurentide Ice Sheet.

SEA-LEVEL RISE. Sea-level rise due to polar ice-sheet mass loss during past warm periods

Interdisciplinary studies of geologic archives have ushered in a new era of deciphering magnitudes, rates, and sources of sea-level rise from polar ice-sheet loss during past warm periods. Accounting for glacial isostatic processes helps to reconcile spatial variability in peak sea level during marine isotope stages 5e and 11, when the global mean reached 6 to 9 meters and 6 to 13 meters higher than present, respectively. Dynamic topography introduces large uncertainties on longer time scales, precluding robust sea-level estimates for intervals such as the Pliocene. Present climate is warming to a level associated with significant polar ice-sheet loss in the past. Here, we outline advances and challenges involved in constraining ice-sheet sensitivity to climate change with use of paleo-sea level records.

Gradual demise of a thin southern Laurentide ice sheet recorded by Mississippi drainage

At the Last Glacial Maximum (LGM), about 21,000 years before present, land-based ice sheets held enough water to reduce global mean sea level by 130 metres. Yet after decades of study, major uncertainties remain as to the distribution of that ice. Here we test four reconstructions of North American deglacial ice-sheet history by quantitatively connecting them to high-resolution oxygen isotope (δ(18)O) records from the Gulf of Mexico using a water mixing model. For each reconstruction, we route meltwater and seasonal runoff through the time-evolving Mississippi drainage basin, which co-evolves with ice geometry and changing topography as ice loads deform the solid Earth and produce spatially variable sea level in a process known as glacial isostatic adjustment. The δ(18)O records show that the Mississippi-drained southern Laurentide ice sheet contributed only 5.4 ± 2.1 metres to global sea level rise, of which 0.66 ± 0.07 metres were released during the meltwater pulse 1A event 14,650-14,310 years before present, far less water than previously thought. In contrast, the three reconstructions based on glacial isostatic adjustment overpredict the δ(18)O-based post-LGM meltwater volume by a factor of 1.6 to 3.6. The fourth reconstruction, which is based on ice physics, has a low enough Mississippi-routed meltwater discharge to be consistent with δ(18)O constraints, but also contains the largest LGM North American ice volume. This suggests that modelling based on ice physics may be the best way of matching isotopic records while also sequestering enough water in the North American ice sheets to match the observed LGM sea level fall.

Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume

The growth and reduction of Northern Hemisphere ice sheets over the past million years is dominated by an approximately 100,000-year periodicity and a sawtooth pattern (gradual growth and fast termination). Milankovitch theory proposes that summer insolation at high northern latitudes drives the glacial cycles, and statistical tests have demonstrated that the glacial cycles are indeed linked to eccentricity, obliquity and precession cycles. Yet insolation alone cannot explain the strong 100,000-year cycle, suggesting that internal climatic feedbacks may also be at work. Earlier conceptual models, for example, showed that glacial terminations are associated with the build-up of Northern Hemisphere 'excess ice', but the physical mechanisms underpinning the 100,000-year cycle remain unclear. Here we show, using comprehensive climate and ice-sheet models, that insolation and internal feedbacks between the climate, the ice sheets and the lithosphere-asthenosphere system explain the 100,000-year periodicity. The responses of equilibrium states of ice sheets to summer insolation show hysteresis, with the shape and position of the hysteresis loop playing a key part in determining the periodicities of glacial cycles. The hysteresis loop of the North American ice sheet is such that after inception of the ice sheet, its mass balance remains mostly positive through several precession cycles, whose amplitudes decrease towards an eccentricity minimum. The larger the ice sheet grows and extends towards lower latitudes, the smaller is the insolation required to make the mass balance negative. Therefore, once a large ice sheet is established, a moderate increase in insolation is sufficient to trigger a negative mass balance, leading to an almost complete retreat of the ice sheet within several thousand years. This fast retreat is governed mainly by rapid ablation due to the lowered surface elevation resulting from delayed isostatic rebound, which is the lithosphere-asthenosphere response. Carbon dioxide is involved, but is not determinative, in the evolution of the 100,000-year glacial cycles.