A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

On the reduced sensitivity of the Atlantic overturning to Greenland ice sheet melting in projections: a multi-model assessment

Large uncertainties exist concerning the impact of Greenland ice sheet melting on the Atlantic meridional overturning circulation (AMOC) in the future, partly due to different sensitivity of the AMOC to freshwater input in the North Atlantic among climate models. Here we analyse five projections from different coupled ocean–atmosphere models with an additional 0.1 Sv (1 Sv = 10 6 m3/s) of freshwater released around Greenland between 2050 and 2089. We find on average a further weakening of the AMOC at 26°N of 1.1 ± 0.6 Sv representing a 27 ± 14% supplementary weakening in 2080–2089, as compared to the weakening relative to 2006–2015 due to the effect of the external forcing only. This weakening is lower than what has been found with the same ensemble of models in an identical experimen - tal set-up but under recent historical climate conditions. This lower sensitivity in a warmer world is explained by two main factors. First, a tendency of decoupling is detected between the surface and the deep ocean caused by an increased thermal stratification in the North Atlantic under the effect of global warming. This induces a shoaling of ocean deep ventilation through convection hence ventilating only intermediate levels. The second important effect concerns the so-called Canary Current freshwater leakage; a process by which additionally released fresh water in the North Atlantic leaks along the Canary Current and escapes the convection zones towards the subtropical area. This leakage is increasing in a warming climate, which is a consequence of decreasing gyres asymmetry due to changes in Ekman rumping. We suggest that these modifications are related with the northward shift of the jet stream in a warmer world. For these two reasons the AMOC is less susceptible to freshwater perturbations (near the deep water formation sides) in the North Atlantic as compared to the recent historical climate conditions. Finally, we propose a bilinear model that accounts for the two former processes to give a conceptual explanation about the decreasing AMOC sensitivity due to freshwater input. Within the limit of this bilinear model, we find that 62 ± 8% of the reduction in sensitivity is related with the changes in gyre asymmetry and freshwater leakage and 38 ± 8% is due to the reduction in deep ocean ventilation associated with the increased stratification in the North Atlantic.

Elevation lines of modelled ice sheet extension during last glacial maximum

A numerical ice-sheet model was used to reconstruct the Late Weichselian glaciation of the Eurasian High Arctic, between Franz Josef Land and Severnaya Zemlya. An ice sheet was developed over the entire Eurasian High Arctic so that ice flow from the central Barents and Kara seas toward the northern Russian Arctic could be accounted for. An inverse approach to modeling was utilized, where ice-sheet results were forced to be compatible with geological information indicating ice-free conditions over the Taymyr Peninsula during the Late Weichselian. The model indicates complete glaciation of the Barents and Kara seas and predicts a "maximum-sized" ice sheet for the Late Weichselian Russian High Arctic. In this scenario, full-glacial conditions are characterized by a 1500-m-thick ice mass over the Barents Sea, from which ice flowed to the north and west within several bathymetric troughs as large ice streams. In contrast to this reconstruction, a "minimum" model of glaciation involves restricted glaciation in the Kara Sea, where the ice thickness is only 300 m in the south and which is free of ice in the north across Severnaya Zemlya. Our maximum reconstruction is compatible with geological information that indicates complete glaciation of the Barents Sea. However, geological data from Severnaya Zemlya suggest our minimum model is more relevant further east. This, in turn, implies a strong paleoclimatic gradient to colder and drier conditions eastward across the Eurasian Arctic during the Late Weichselian.

(Appendix A) Geochemical composition of snow and eolian mineral dust from Berkner Island ice sheet, Antarctica

The Sr and Nd isotopic composition of dust extracted from recent snow layers at the top of Berkner Island ice sheet (located within the Filchner-Ronne Ice Shelf at the southern end of the Weddell Sea) enables us, for the first time, to document dust provenance in Antarctica outside the East Antarctic Plateau (EAP) where all previous studies based on isotopic fingerprinting were carried out. Berkner dust displays an overall crust-like isotopic signature, characterized by more radiogenic 87Sr/86Sr and much less radiogenic 143Nd/144Nd compared to dust deposited on the EAP during glacial periods. Differences with EAP interglacial dust are not as marked but still significant, indicating that present-day Berkner dust provenance is distinct, at least to some extent, from that of the dust reaching the EAP. The fourteen snow-pit sub-seasonal samples that were obtained span a two-year period (2002-2003) and their dust Sr and Nd isotopic composition reveals that multiple sources are at play over a yearly time period. Southern South America, Patagonia in particular, likely accounts for part of the observed spring/summer dust deposition maxima, when isotopic composition is shifted towards 'younger' isotopic signatures. In the spring, possible additional inputs from Australian sources would also be supported by the data. Most of the year, however, the measured isotopic signatures would be best explained by a sustained background supply from putative local sources in East Antarctica, which carry old-crust-like isotopic fingerprints. Whether the restricted East Antarctic ice-free areas produce sufficient eolian material has yet to be substantiated however. The fact that large (> 5 µm) particles represent a significant fraction of the samples throughout the entire time-series supports scenarios that involve contributions from proximal sources, either in Patagonia and/or Antarctica (possibly including snow-free areas in the Antarctic Peninsula and other areas as well). This also indicates that additional dust transport, which does not reach the EAP, must occur at low-tropospheric levels to this coastal sector of Antarctica.

(Table 1) GRACE-derived mass balance of the Greenland ice sheet 2002-2010

We re-evaluate the Greenland mass balance for the recent period using low-pass Independent Component Analysis (ICA) post-processing of the Level-2 GRACE data (2002-2010) from different official providers (UTCSR, JPL, GFZ) and confirm the present important ice mass loss in the range of -70 and -90 Gt/y of this ice sheet, due to negative contributions of the glaciers on the east coast. We highlight the high interannual variability of mass variations of the Greenland Ice Sheet (GrIS), especially the recent deceleration of ice loss in 2009-2010, once seasonal cycles are robustly removed by Seasonal Trend Loess (STL) decomposition. Interannual variability leads to varying trend estimates depending on the considered time span. Correction of post-glacial rebound effects on ice mass trend estimates represents no more than 8 Gt/y over the whole ice sheet. We also investigate possible climatic causes that can explain these ice mass interannual variations, as strong correlations between GRACE-based mass balance and atmosphere/ocean parallels are established: (1) changes in snow accumulation, and (2) the influence of inputs of warm ocean water that periodically accelerate the calving of glaciers in coastal regions and, feed-back effects of coastal water cooling by fresh currents from glaciers melting. These results suggest that the Greenland mass balance is driven by coastal sea surface temperature at time scales shorter than accumulation.

Seismic investigations of the Greenland ice sheet at EGIG points

The compressional and the shear wave velocities in the Greenland ice sheet are derived from seismic records of the EGIG 1959. Further the variation of velocities in the firn and the dependance of Poisson's ratio from depth are determined. At Station Centrale, two P-waves are recorded from underground layers. Their velocities show that the ice basement consists of crystalline rocks. The P-wave velocities derived from reflections agree well with those obtained by refraction shooting. From this agreement results that the ice is ± homogenous and ± isotropic for Pwaves. The elastic constants for isotropic ice are calculated. Finally the temperature dependence of the velocities is discussed.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

(Table 1) Basin area, average late-summer speedup, number of sinks and runoff of West Greenland Ice Sheet catchments

We use interferometric synthetic aperture radar observations recorded in a land-terminating sector of western Greenland to characterise the ice sheet surface hydrology and to quantify spatial variations in the seasonality of ice sheet flow. Our data reveal a non-uniform pattern of late-summer ice speedup that, in places, extends over 100 km inland. We show that the degree of late-summer speedup is positively correlated with modelled runoff within the 10 glacier catchments of our survey, and that the pattern of late-summer speedup follows that of water routed at the ice sheet surface. In late-summer, ice within the largest catchment flows on average 48% faster than during winter, whereas changes in smaller catchments are less pronounced. Our observations show that the routing of seasonal runoff at the ice sheet surface plays an important role in shaping the magnitude and extent of seasonal ice sheet speedup.

Simulated temperature changes as a result of Antarctic ice sheet growth and CO2 reduction

During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere-ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.