Search

The total air content (TAC) of polar ice cores has long been considered a potential proxy for past ice sheet elevation. Recent work, however, has shown that a variety of other factors also influence this parameter. In this paper we present a high-resolution TAC record from the South Pole ice core (SPC14) covering the last 54 000 years and discuss the implications of the data for interpreting TAC from ice cores. The SPC14 TAC record shows multiple features of interest, including (1) long-term orbital-scale variability, (2) millennial-scale variability in the Holocene and last glacial period, and (3) a period of stability from 35 to 25 ka. The longer, orbital-scale variations in TAC are highly correlated with integrated summer insolation (ISI), corroborating the potential of TAC to provide an independent dating tool via orbital tuning. Large millennial-scale variability in TAC during the last glacial period is positively correlated with past accumulation rate reconstructions as well as δ15 N-N 2 , a firn thickness proxy. These TAC variations are too large to be controlled by direct effects of temperature and too rapid to be tied to elevation changes. We propose that grain size metamorphism near the firn surface explains these changes. We note, however, that at sites with different climate histories than the South Pole, TAC variations may be dominated by other processes. Our observations of millennial-scale variations in TAC show a different relationship with accumulation rate than observed at sites in Greenland. [ABSTRACT FROM AUTHOR]

The North Sea Last Interglacial sea level is sensitive to the fingerprint of mass loss from polar ice sheets. However, the signal is complicated by the influence of glacial isostatic adjustment driven by Penultimate Glacial Period ice-sheet changes, and yet these ice-sheet geometries remain significantly uncertain. Here, we produce new reconstructions of the Eurasian ice sheet during the Penultimate Glacial Maximum (PGM) by employing large ensemble experiments from a simple ice-sheet model that depends solely on basal shear stress, ice extent, and topography. To explore the range of uncertainty in possible ice geometries, we use a parameterised shear-stress map as input that has been developed to incorporate bedrock characteristics and the influence of ice-sheet basal processes. We perform Bayesian uncertainty quantification, utilising Gaussian process emulation, to calibrate against global ice-sheet reconstructions of the Last Deglaciation and rule out combinations of input parameters that produce unrealistic ice sheets. The refined parameter space is then applied to the PGM to create an ensemble of constrained 3D Eurasian ice-sheet geometries. Our reconstructed PGM Eurasian ice-sheet volume is 48±8 m sea-level equivalent (SLE). We find that the Barents–Kara Sea region displays both the largest mean volume and volume uncertainty of 24±8 m SLE while the British–Irish sector volume of 1.7±0.2 m SLE is the smallest. Our new workflow may be applied to other locations and periods where ice-sheet histories have limited empirical data. [ABSTRACT FROM AUTHOR]

We investigate the sensitivity of a distributed glacier surface mass and energy balance model using a variance-based analysis, for two distinct periods of the last glacial cycle: the present day (PD) and the Last Glacial Maximum (LGM). The results can be summarized in three major findings: the sensitivity towards individual model parameters and parameterizations is as variable in space as it is in time. The model is most sensitive to uncertainty related to atmospheric emissivity and the down-welling longwave radiation. While the turbulent latent heat flux has a sizable contribution to the surface mass balance uncertainty in central Greenland today, it dominates over the entire ice sheet during the cold climate of the LGM, in spite of its low impact on the overall surface mass balance of the Greenland ice sheet in the modern climate. We conclude that quantifying the model sensitivity is very helpful for tuning free model parameters because it clarifies the relative importance of individual parameters and highlights interactions between them that need to be considered. [ABSTRACT FROM AUTHOR]

Outburst floods from subglacial lakes beneath the Antarctic Ice Sheet modulate ice-flow velocities over periods of months to years. Although subglacial lake drainage events have been observed from satellite-altimetric data, little is known about their role in the long-term evolution of ice-sheet basal hydrology. Here, we systematically map and model past water flow through an extensive area containing over 1000 subglacial channels and 19 former lake basins exposed on over 19 000 km 2 of seafloor by the retreat of Pine Island and Thwaites glaciers, West Antarctica. At 507 m wide and 43 m deep on average, the channels offshore of present-day Pine Island and Thwaites glaciers are approximately twice as deep, 3 times as wide, and cover an area over 400 times larger than the terrestrial meltwater channels comprising the Labyrinth in the Antarctic Dry Valleys. The channels incised into bedrock offshore of contemporary Pine Island and Thwaites glaciers would have been capable of accommodating discharges of up to 8.8×106 m 3 s -1. We suggest that the channels were formed by episodic discharges from subglacial lakes trapped during ice-sheet advance and retreat over multiple glacial periods. Our results document the widespread influence of episodic subglacial drainage events during past glacial periods, in particular beneath large ice streams similar to those that continue to dominate contemporary ice-sheet discharge. [ABSTRACT FROM AUTHOR]

The European Alps, the cradle of pioneering glacial studies, are one of the regions where geological markers of past glaciations are most abundant and well-studied. Such conditions make the region ideal for testing numerical glacier models based on simplified ice flow physics against field-based reconstructions and vice versa. Here, we use the Parallel Ice Sheet Model (PISM) to model the entire last glacial cycle (120-0 ka) in the Alps, using horizontal resolutions of 2 and 1 km. Climate forcing is derived using two sources: present-day climate data from WorldClim and the ERA-Interim reanalysis; time-dependent temperature offsets from multiple palaeo-climate proxies. Among the latter, only the European Project for Ice Coring in Antarctica (EPICA) ice core record yields glaciation during marine oxygen isotope stages 4 (69-62 ka) and 2 (34-18 ka). This is spatially and temporally consistent with the geological reconstructions, while the other records used result in excessive early glacial cycle ice cover and a late Last Glacial Maximum. Despite the low variability of this Antarcticbased climate forcing, our simulation depicts a highly dynamic ice sheet, showing that Alpine glaciers may have advanced many times over the foreland during the last glacial cycle. Ice flow patterns during peak glaciation are largely governed by subglacial topography but include occasional transfluences through the mountain passes. Modelled maximum ice surface is on average 861m higher than observed trimline elevations in the upper Rhône Valley, yet our simulation predicts little erosion at high elevation due to cold-based ice. Finally, despite the uniform climate forcing, differences in glacier catchment hypsometry produce a timetransgressive Last Glacial Maximum advance, with some glaciers reaching their modelled maximum extent as early as 27 ka and others as late as 21 ka. [ABSTRACT FROM AUTHOR]

Systematically monitoring Greenland's outlet glaciers is central to understanding the timescales over which their flow and sea level contributions evolve. In this study we use data from the new Sentinel-1a/b satellite constellation to generate 187 velocity maps, covering four key outlet glaciers in Greenland: Jakobshavn Isbræ, Petermann Glacier, Nioghalvfjerdsfjorden, and Zachariæ Isstrøm. These data provide a new high temporal resolution record (6-day averaged solutions) of each glacier's evolution since 2014, and resolve recent seasonal speedup periods and inter-annual changes in Greenland outlet glacier speed with an estimated certainty of 10 %. We find that since 2012, Jakobshavn Isbræ has been decelerating, and now flows approximately 1250myr-1 (10 %), slower than 5 years previously, thus reversing an increasing trend in ice velocity that has persisted during the last decade. Despite this, we show that seasonal variability in ice velocity remains significant: up to 750myr-1 (14 %) at a distance of 12 km inland of the terminus. We also use our new dataset to estimate the duration of speedup periods (80-95 days) and to demonstrate a strong relationship between ice front position and ice flow at Jakobshavn Isbræ, with increases in speed of ~1800myr-1 in response to 1 km of retreat. Elsewhere, we record significant seasonal changes in flow of up to 25% (2015) and 18% (2016) at Petermann Glacier and Zachariæ Isstrøm, respectively. This study provides a first demonstration of the capacity of a new era of operational radar satellites to provide frequent and timely monitoring of ice sheet flow, and to better resolve the timescales over which glacier dynamics evolve. [ABSTRACT FROM AUTHOR]

After more than a century of geological research, the Cordilleran ice sheet of North America remains among the least understood in terms of its former extent, volume, and dynamics. Because of the mountainous topography on which the ice sheet formed, geological studies have often had only local or regional relevance and shown such a complexity that ice-sheet-wide spatial reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (62.2-56.9 ka) and 2 (23.2-16.9 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.7 ka). [ABSTRACT FROM AUTHOR]

We present a new seafloor map for the northern Antarctic Peninsula (AP), including swath multibeam data sets from five national programs. Our map allows for the examination and interpretation of Last Glacial Maximum (LGM) paleo-ice-flow paths developed on the seafloor from the preservation of mega-scale glacial lineations, drumlinized features, and selective linear erosion. We combine this with terrestrial observations of flow direction to place constraints on ice divides and ice domes on the AP continental shelf during the LGM time interval. The results show a flow bifurcation as ice exits the Larsen B embayment. Flow emanating off the Seal Nunataks (including Robertson Island) is directed toward the southeast, then eastward as the flow transits toward the Robertson Trough. A second, stronger "streaming flow" is directed toward the southeast, then southward as ice overflowed the tip of the Jason Peninsula to reach the southern perimeter of the embayment. Our reconstruction also refines the extent of at least five other distinct paleo-ice-stream systems that, in turn, serve to delineate seven broad regions where contemporaneous ice domes must have been centered on the continental shelf at LGM. Our reconstruction is more detailed than other recent compilations because we followed specific ice-flow indicators and have kept tributary flow paths parallel. [ABSTRACT FROM AUTHOR]

The article offers information on the ice melting of the Northern Greenland during the interglaciation. It states that the stimulation on the Greenland ice sheet (GrIS) during the Eemian interglaciation with the three-dimensional ice sheet model SICOPOLIS was performed. It mentions that the northeastern part of the GrIS is unstable and retreats significantly, despite moderate melt rates.