The ice sheets covering Antarctica and Greenland present the greatest uncertainty in, and largest potential contribution to, future sea level rise. The uncertainty arises from a paucity of suitable observations covering the full range of ice sheet behaviors, incomplete understanding of the influences of diverse processes, and limitations in defining key boundary conditions for the numerical models. To investigate the impact of these uncertainties on ice sheet projections we undertook a structured expert judgement study. Here, we interrogate the findings of that study to identify the dominant drivers of uncertainty in projections and their relative importance as a function of ice sheet and time. We find that for the 21st century, Greenland surface melting, in particular the role of surface albedo effects, and West Antarctic ice dynamics, specifically the role of ice shelf buttressing, dominate the uncertainty. The importance of these effects holds under both a high-end 5°C global warming scenario and another that limits global warming to 2°C. During the 22nd century the dominant drivers of uncertainty shift. Under the 5°C scenario, East Antarctic ice dynamics dominate the uncertainty in projections, driven by the possible role of ice flow instabilities. These dynamic effects only become dominant, however, for a temperature scenario above the Paris Agreement 2°C target and beyond 2100. Our findings identify key processes and factors that need to be addressed in future modeling and observational studies in order to reduce uncertainties in ice sheet projections.

The conventional sea level budget (SLB) equates changes in sea surface height with the sum of ocean mass and steric change, where solid-Earth movements are included as corrections but limited to the impact of glacial isostatic adjustment. However, changes in ocean mass load also deform the ocean bottom elastically. Until the early 2000s, ocean mass change was relatively small, translating into negligible elastic ocean bottom deformation (OBD), hence neglected in the SLB equation. However, recently ocean mass has increased rapidly; hence, OBD is no longer negligible and likely of similar magnitude to the deep steric sea level contribution. Here, we use a mass-volume framework, which allows the ocean bottom to respond to mass load, to derive a SLB equation that includes OBD. We discuss the theoretical appearance of OBD in the SLB equation and its implications for the global SLB.

Despite considerable advances in process understanding, numerical modeling, and the observational record of ice sheet contributions to global mean sea-level rise (SLR) since the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change, severe limitations remain in the predictive capability of ice sheet models. As a consequence, the potential contributions of ice sheets remain the largest source of uncertainty in projecting future SLR. Here, we report the findings of a structured expert judgement study, using unique techniques for modeling correlations between inter- and intra-ice sheet processes and their tail dependences. We find that since the AR5, expert uncertainty has grown, in particular because of uncertain ice dynamic effects. For a +2 °C temperature scenario consistent with the Paris Agreement, we obtain a median estimate of a 26 cm SLR contribution by 2100, with a 95th percentile value of 81 cm. For a +5 °C temperature scenario more consistent with unchecked emissions growth, the corresponding values are 51 and 178 cm, respectively. Inclusion of thermal expansion and glacier contributions results in a global total SLR estimate that exceeds 2 m at the 95th percentile. Our findings support the use of scenarios of 21st century global total SLR exceeding 2 m for planning purposes. Beyond 2100, uncertainty and projected SLR increase rapidly. The 95th percentile ice sheet contribution by 2200, for the +5 °C scenario, is 7.5 m as a result of instabilities coming into play in both West and East Antarctica. Introducing process correlations and tail dependences increases estimates by roughly 15%.

We clarify key aspects of the evaluation, by de Vries and van de Wal (2015), of our expert elicitation paper on the contributions of ice sheet melting to sea level rise due to future global temperature rise scenarios (Bamber and Aspinall 2013), and extend the conversation with further analysis of their proposed approach for combining expert uncertainty judgments.

The glacial history of Antarctica during the most recent Milankovitch cycles is poorly constrained relative to the Northern Hemisphere. As a consequence, the contribution of mass changes in the Antarctic ice sheet to global sea-level change and the prediction of its future evolution remain uncertain. The process of Glacial Isostatic Adjustment (GIA) represents the ongoing response of the solid Earth to the Late-Pleistocene deglaciation and, therefore, provides information about Antarctic glacial history. Moreover, insufficient knowledge of GIA hampers the determination of present-day changes in the Antarctic mass balance through satellite gravity measurements. Previous studies have laid the theoretical foundation for distinguishing between signals of ongoing GIA and contemporary ice mass change through the combination of satellite gravimetry and satellite altimetry. This distinction is made possible by the fact the GIA-induced changes (involving relatively dense rock) will produce a different combination of topography and gravity change than those produced by variations in ice or firn thickness (due to the lower density of these materials); however, no conclusive results have been produced to date. Here we show that, by combining laser altimetry and gravity data from the ICESat and GRACE satellite missions over the period March 2003-March 2008, the GIA contribution can indeed be isolated. The inferred GIA signal over the Antarctic continent, which represents the first result derived from direct observations by satellite techniques, strongly supports Late-Pleistocene ice models derived from glacio-geologic studies. The GIA impact on GRACE-derived estimates of mass balance is found to be 100 ± 67 Gt/yr.