Description: One of the most dramatic signs of ongoing global change is the mass loss of the Greenland Ice Sheet and the resulting rise in sea level, whereby most of the recent ice sheet mass loss can be attributed to an increase in meltwater runoff. The retreat and thinning of Greenland glaciers has been caused by rising air and ocean temperatures over the past decades. Despite the global scale impact of the changing ice sheet balance, estimates of glacial runoff in Greenland rarely extend past several decades, thus limiting our understanding of long-term glacial response to temperature. Here we present a 42-year long annually resolved red coralline algal Mg/Ca proxy temperature record from a southwestern Greenland fjord, with temperature ranging from 1.5 to 4 °C (standard error = 1.06 °C). This temperature time series in turn tracks the general trend of glacial runoff from four West Greenland glaciers discharging freshwater into the fjord (all p 〈 0.001). The algal time series further exhibits significant correlations to Irminger Sea temperature patterns, which are transmitted to western Greenland fjords via the West Greenland Current. The 42-year long record demonstrates the potential of annual increment forming coralline algae, which are known to live up to 650 years and which are abundant along the Greenland coastline, for reconstructing time series of sea surface temperature.

Description: Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continent’s interior1,2,3, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves4,5,6,7,8,9,10,11. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85% of Antarctica’s coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4±61 Gt yr−1. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59% in 10 years to reach 132±60 Gt yr−1 in 2006. In the Peninsula, losses increased by 140% to reach 60±46 Gt yr−1 in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance.

Description: The Greenland Ice Sheet (GrIS) mass loss has been accelerating at a rate of about 20 ± 10 Gt/yr2 since the end of the 1990's, with around 60 % of this mass loss directly attributed to enhanced surface meltwater runoff. However, in the climate and glaciology communities, different approaches exist on how to model the different surface mass balance (SMB) components using: (1) complex physically-based climate models which are computationally expensive; (2) intermediate complexity energy balance models; (3) simple and fast positive degree day models which base their inferences on statistical principles and are computationally highly efficient. Additionally, many of these models compute the SMB components based on different spatial and temporal resolutions, with different forcing fields as well as different ice sheet topographies and extents, making inter-comparison difficult. In the GrIS SMB model intercomparison project (GrSMBMIP) we address these issues by forcing each model with the same data (i.e., the ERA-Interim reanalysis) except for two global models for which this forcing is limited to the oceanic conditions, and at the same time by interpolating all modelled results onto a common ice sheet mask at 1 km horizontal resolution for the common period 1980–2012. The SMB outputs from 13 models are then compared over the GrIS to (1) SMB estimates using a combination of gravimetric remote sensing data from GRACE and measured ice discharge, (2) ice cores, snow pits, in-situ SMB observations, and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Our results reveal that the mean GrIS SMB of all 13 models has been positive between 1980 and 2012 with an average of 340 ± Gt/yr, but has decreased at an average rate of −7.3 Gt/yr2 (with a significance of 96 %), mainly driven by an increase of 8.0 Gt/yr2 (with a significance of 98 %) in meltwater runoff. Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting the need for accurate representation of the GrIS ablation zone extent and processes driving the surface melt. In addition, a higher density of in-situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 mWE/yr due to large discrepancies in modelled snowfall accumulation. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of same order than RCMs with observations and remain then useful tools for long-term simulations. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present day SMB relative to observations, suggesting that biases are not systematic among models.

Description: The land ice contribution to global mean sea level rise has not yet been predicted1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models2,3,4,5,6,7,8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

Description: Iceberg calving has been assumed to be the dominant cause of mass loss for the Antarctic ice sheet, with previous estimates of the calving flux exceeding 2,000 gigatonnes per year. More recently, the importance of melting by the ocean has been demonstrated close to the grounding line and near the calving front. So far, however, no study has reliably quantified the calving flux and the basal mass balance (the balance between accretion and ablation at the ice-sheet base) for the whole of Antarctica. The distribution of fresh water in the Southern Ocean and its partitioning between the liquid and solid phases is therefore poorly constrained. Here we estimate the mass balance components for all ice shelves in Antarctica, using satellite measurements of calving flux and grounding-line flux, modelled ice-shelf snow accumulation rates and a regional scaling that accounts for unsurveyed areas. We obtain a total calving flux of 1,321 ± 144 gigatonnes per year and a total basal mass balance of -1,454 ± 174 gigatonnes per year. This means that about half of the ice-sheet surface mass gain is lost through oceanic erosion before reaching the ice front, and the calving flux is about 34 per cent less than previous estimates derived from iceberg tracking. In addition, the fraction of mass loss due to basal processes varies from about 10 to 90 per cent between ice shelves. We find a significant positive correlation between basal mass loss and surface elevation change for ice shelves experiencing surface lowering and enhanced discharge. We suggest that basal mass loss is a valuable metric for predicting future ice-shelf vulnerability to oceanic forcing.

Description: A 21-year record is presented of surface mass balance measurements along the K-transect. The series covers the period 1990-2011. Data are available at 8 sites along a transect over an altitude range of 390 - 1850 m at approximately 67° N in West Greenland. The surface mass balance gradient is on average 3.8 x 10**-3 m w.e./m, and the mean equilibrium line altitude is 1553 m a.s.l. Only the lower 3 sites within 10 km of the margin experience a significant increasing trend in the ablation over the entire period.

Description: The firn density, temperature and liquid water content of the Greenland ice sheet have been modelled with the IMAU-FDM firn model. IMAU-FDM is forced at the surface with the latest output of the regional climate model RACMO2.3p2. The data is on a horizontal grid of 11x11 km and covers 1960-2016 with a 10-day temporal resolution. Here, time series of the firn air content (vertically integrated difference between firn and ice density (= 917 kg m-3)) and 10-m firn temperature are provided. All other IMAU-FDM output is available from the authors without conditions.