Articles | Volume 15, issue 12
The Cryosphere, 15, 5705–5715, 2021
The Cryosphere, 15, 5705–5715, 2021

Brief communication 17 Dec 2021

Brief communication | 17 Dec 2021

Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level

Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level
Andy Aschwanden1,, Timothy C. Bartholomaus2,, Douglas J. Brinkerhoff3,, and Martin Truffer1, Andy Aschwanden et al.
  • 1Geophysical Institute, University of Alaska Fairbanks, Fairbanks, Alaska, USA
  • 2Department of Earth and Spatial Sciences, University of Idaho, Moscow, Idaho, USA
  • 3Department of Computer Science, University of Montana, Missoula, MT, USA
  • These authors contributed equally to this work.

Correspondence: Andy Aschwanden (


Accurately projecting mass loss from ice sheets is of critical societal importance. However, despite recent improvements in ice sheet models, our analysis of a recent effort to project ice sheet contribution to future sea level suggests that few models reproduce historical mass loss accurately and that they appear much too confident in the spread of predicted outcomes. The inability of models to reproduce historical observations raises concerns about the models' skill at projecting mass loss. Here we suggest that uncertainties in the future sea level contribution from Greenland and Antarctica may well be significantly higher than reported in that study. We propose a roadmap to enable a more realistic accounting of uncertainties associated with such forecasts and a formal process by which observations of mass change should be used to refine projections of mass change. Finally, we note that tremendous government investment and planning affecting tens to hundreds of millions of people is founded on the work of just a few tens of scientists. To achieve the goal of credible projections of ice sheet contribution to sea level, we strongly believe that investment in research must be commensurate with the scale of the challenge.

1 Sea level rise predictions from ice sheet loss

Global sea level rose during the 20th century more than 3 times faster than at any time during the last 2000 years (Kopp et al.2016). Over the last several decades, mass loss from the Greenland ice sheet has been the fastest-growing contributor to this rise (Chen et al.2017; Rietbroek et al.2016), currently tracking the upper-end estimates of the Intergovernmental Panel on Climate Change's (IPCC) fifth assessment report (AR5; IPCC2013). Sea level rise driven by global warming is expected to continue over the coming century, potentially flooding 14–322 million people per year in 2100 and reducing annual global gross domestic production by as much as 9 % (Hinkel et al.2014). To guide planning for and mitigation of anticipated damages, the IPCC published a suite of sea level rise projections for the remainder of the 21st century: its sixth assessment report (AR6; Masson-Delmotte et al.2021). Effective planning for coming sea level rise necessitates that these estimates be credible but also that they be accompanied by a defensible assessment of uncertainty (Moon et al.2020).

Ice sheet models have emerged as the de facto standard for generating estimates of ice sheet contribution to sea level rise and are the basis for AR6's estimates for the next century, particularly when coupled to models that simulate relevant forcing from the atmosphere and ocean. However, defensible estimates of the uncertainty in estimates produced by these models remain one of the most challenging goals of scientific inquiry. Here we take ice sheet models to comprise computer code used to predict change in ice sheet mass, including ice dynamics, models of surface mass balance, and models of the ice–ocean boundary.

AR6's estimate of mass change for both Greenland and Antarctica over the next century is based primarily on the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) (Nowicki et al.2016, 2020; Payne et al.2021; Seroussi et al.2020; Goelzer et al.2020). Through generous collaboration and leadership, 21 groups from around the world contributed 37 different models of Greenland and Antarctic ice sheet change through a set of core and optional experiments and corresponding historical simulations. These simulations had been performed before the latest socio-economic scenarios and climate models were available. Edwards et al. (2021) statistically reproduced the response of these models to climate change and then estimated glacier and ice sheet response to the most recent scenarios and climate model outputs. These Edwards et al. (2021) projections comprise the AR6 estimate for Greenland. For Antarctica, the AR6 estimate is created by averaging an alternative model intercomparison (Levermann et al.2020) with these emulated ISMIP6 projections.

We believe that the results of ISMIP6 represent the state of the art in terms of understanding ice sheet model variability and the breadth of behavior that ice sheet models encompass, and the lessons that it provides will be far-reaching both on their own and also with respect to planning additional collaborative efforts to assess uncertainty in sea level prediction. Nonetheless, it is our perspective that the IPCC AR6 estimates of sea level contribution from ice sheets that are based on ISMIP6 may not be accurate and that the accompanying uncertainties do not reflect the true breadth of uncertainties associated with ice sheet change. Our skepticism is based on the premise that accurate predictions of the cryosphere's contribution to sea level require that models

  1. fully characterize uncertainties in model structure, parameters, initial conditions, and boundary conditions;

  2. yield simulations that fit observations within observational uncertainty.

If the first point is not satisfied, then predictive uncertainties are likely to be underestimated. If the second condition is not satisfied, then the distribution of model predictions is likely to be biased relative to reality. For the purpose of credibly projecting mass change, we assert that the models must accurately reproduce observed mass change. Although such validation is insufficient on its own to instill confidence in projections, it is a necessary condition for projection. These points are not just of academic importance; they can lead to a false sense of security when planning coastal infrastructure and preparing for future sea level rise, with potentially dire consequences.

Our concern is illustrated by comparing the ISMIP6 simulations of mass loss from the Greenland and Antarctic ice sheets between 2000 and 2025 (Goelzer et al.2020) with observations of mass loss (The IMBIE Team2019) (see “Methods” section for details). This 25-year period begins with 15 years of ISMIP6 historical simulations, during which modeling groups were free to select climate forcings necessary to bring their modeled ice sheets to the start of the projection period in 2015. To visualize how the historical simulations impact the projections, we also show the first 10 years of the projection period, during which surface mass balance and temperature anomalies were imposed uniformly on all ice sheet models (Nowicki et al.2020).

For Greenland, a clear picture emerges (Fig. 1), where ISMIP6 simulations systematically underestimate recent relative cumulative mass loss when compared to observations. Indeed, the 95th percentile of the ISMIP6 experiments follows the observed record of mass loss, which implies that either the model uncertainty is underestimated, the ensemble is a biased predictor of cumulative mass change, or both. For Antarctica, the picture is substantially different, with the spread in simulations much larger than observational uncertainty. In this case, many simulations that are not consistent with observations yield a predictive variance that is also not reliable, albeit in the direction of too little specificity.

Figure 1Observed and simulated historical mass changes from the Antarctic ice sheet (AIS) and Greenland ice sheet (GIS) between 2000 and 2020 in gigatons (Gt) and centimeters of sea level equivalent (cm SLE). A consensus estimate of observed mass changes (The IMBIE team2018; The IMBIE Team2019) is plotted in blue along with their respective uncertainties (shaded). The ensembles of ISMIP6 (Goelzer et al.2020; Seroussi et al.2020) historical simulations and projections are plotted with dark-gray lines, and the 5th to 95th percentile mass loss rates are shown as a 90 % credibility interval with light-gray shading. Due to the large variance in ISMIP6 historical simulations for Antarctica, the uncertainties in IMBIE are not visible in the plot.


In an effort to assess why these patterns appear and to guide future efforts, we recast the problem of ice sheet simulation in terms of the characterization of joint probability distributions and assess how our two conditions above relate to this viewpoint. We then sketch a path forward for robustly characterizing the potential ice sheet contribution to sea level over the coming century. We note that our commentary is not intended to be a comprehensive review of uncertainty quantification in ice sheet modeling but rather to serve as a first attempt to define a consistent language around which community efforts can be discussed. We also note that the commentary should not be seen as a criticism of the ISMIP6 effort, which we regard as important and successful and in which some of us have actively participated, but rather as pointing towards important work still to be done.

2 Quantifying uncertainties

For the practical problem of predicting the ice sheets' contribution to sea level, we find it useful to adopt a probabilistic framework. In that framework, we seek to establish a credibility bound (say 90 %) between which sea level contribution will fall with that pre-supposed probability. Such an interval can readily be constructed from a probability density function (PDF) for the cryosphere's contribution to global sea level by computing quantiles, and thus this is the function that experiments aiming to quantify sea level contribution must correctly characterize. We write this predictive distribution as

(1) P ( Δ z | F ) ,

where Δz is sea level contribution, and represents climate forcing scenarios (i.e., greenhouse gas emissions) expressed as, for example, Representative Concentration Pathways (RCP) or Shared Socioeconomic Pathways, which should also be characterized by their own PDF. In this short communication we will not address the issue of uncertainty in the specification of a forcing scenario (The IMBIE Team2019) but concentrate instead on the uncertainties arising solely from ice sheet models.

While interpretation of Pz|ℱ) is straightforward, its accurate construction is a grand scientific challenge. The standard approach involves running computer programs that approximately solve mathematical equations describing our best understanding of ice sheet physics. In the best case, all facets of a physical system are known (including initial and boundary conditions), the equations describing those systems are complete and deterministic, and the mechanism of solution is perfect. In this idealized situation there is no predictive uncertainty in sea level contribution, and Pz|ℱ) can be characterized with a single model run. In practice, several types of uncertainties complicate the issue and introduce bias and variance in the predictions. In the following, we discuss these different categories of uncertainty as they pertain to the problem of sea level contribution.

2.1 Model uncertainty

The equations used to describe the physical processes in models are invariably an idealization of reality. Indeed, all models are subject to some degree of model error; some physical processes are represented incompletely, while others are omitted altogether. For example, the impact of subglacial hydrology on basal motion (e.g., Bueler and van Pelt2015) and the effect of ice mélange (Amundson et al.2010; Joughin et al.2020) and iceberg calving (Amaral et al.2020) on terminus position remain poorly or not represented in numerical ice sheet models, leading to potentially large model uncertainty (sometimes called structural uncertainty). Some hypothesized processes that are often not included in ice sheet models may lead to critical dynamic instabilities that could deeply affect model evolution (DeConto et al.2021; Sadai et al.2020). On the other hand, the omission of frictional stresses from wind over an ice sheet surface yields a model that is incorrect, yet the resulting error is negligible. Unfortunately, assessing the non-negligible drivers of model inadequacy is a long and arduous process. The different choices that modelers make in this regard leads to an implicitly defined probability distribution P(ℳ), where a particular model is a random (although very likely biased) sample from that distribution. Such model error affects the distribution over sea level contribution as

(2) P ( Δ z | F ) = P ( Δ z | F , M ) P ( M ) d M .

Monte Carlo approximation of this integral is exceptionally challenging because drawing a single sample from P(ℳ) requires the development of a new and (ostensibly) independent ice sheet model, an effort which can take years. However, because many ice sheet models have been developed in parallel, it is now possible to approximate it through model intercomparison (e.g., Bindschadler et al.2013; Nowicki et al.2016; Seroussi et al.2019; Levermann et al.2020) or through structured expert judgment (Bamber et al.2019). We note that an implicit assumption of using models for prediction is that the true data generating process should be contained in P(), which is a questionable assumption indeed.

2.2 Initial state uncertainty

Decade- to century-scale forecasts of ice sheet behavior are sensitive to the initial state, similar to numerical weather forecasts (Vaughan and Arthern2007; Aschwanden et al.2013; Aðalgeirsdóttir et al.2014). Unfortunately, observations alone are insufficient to define an initial state since not all aspects of an ice sheet state are observable to begin with, necessitating the use of data assimilation to combine sparse observational data with models of varying complexity.

Similar to model uncertainty, initial state uncertainty affects the distribution over sea level contribution as

(3) P ( Δ z | F ) = P ( Δ z | F , I ) P ( I ) d I ,

where is an initial state.

Details vary from model to model but generally include initial conditions for the conservation of mass (ice thickness and extent), momentum (basal stress distribution), and energy (temperature or enthalpy).

2.3 Parametric uncertainty

Due to computational and conceptual constraints, there are limits to the level of detail at which processes can be simulated in ice sheet models predicting sea level contribution. For example, the fast and small-scale fracture processes that occur at a marine ice sheet's calving front are more complex than can be reasonably captured in a large-scale model with practical time steps. This gives rise to parameters k=k1,,kN, where N is the number of parameters, specific to a given model, which may include different parameterizations than others. These parameters are explicit numerical values that act as the bridge between un-simulated small-scale processes and their integrated effects at a practical computational scale. Unfortunately, accurate numerical values for such parameters generally do not exist. This lack of knowledge induces parametric uncertainty; for example, different values of thermal conductivity within firn might lead to different predictions of sea level contribution. The predictive distribution under parametric uncertainty is

(4) P ( Δ z | F ) = P ( Δ z | F , k ) P ( k ) d k ,

where P(k) is the probability distribution over a given model's parameter values, which we assume to be independent of scenario. Multiple works have approximately evaluated the equation above for either Greenland or Antarctica using Monte Carlo simulation, which is computationally challenging but conceptually simple (NIAS et al.2016; Schlegel et al.2018; Aschwanden et al.2019; Hill et al.2021): sample a large number of parameter values from P(k), and compute sea level contribution for each sample.

2.4 Aleatoric uncertainty

Ice sheet models additionally have aleatoric uncertainty; i.e., they are subject to irreducibly random processes, most notably the chaotic dynamics present in both atmospheric and oceanic forcings. The predictive distribution under this kind of uncertainty can be decomposed as

(5) P ( Δ z | F ) = P ( Δ z | f ) P ( f | F ) d f ,

where f represents a specific realization of a random forcing, and P(f|ℱ) is its probability distribution under scenario . Due to the relatively slow response time of the cryosphere to such forcings, aleatoric uncertainty often contributes little variance to predictions in sea level contribution over practical timescales of decades to centuries. However, in circumstances where these forcings may interact with a critical glaciological instability like the marine ice sheet instability (Mercer1978), aleatoric uncertainty has the tendency of producing “fat tails”, effectively biasing ice sheet evolution towards more extreme mass loss scenarios (Robel et al.2019). While only a few studies have characterized the distribution over ice sheet responses to aleatoric uncertainty (e.g., Hoffman et al.2019), and its influence is not precisely known, Monte Carlo simulation can be used to understand the effects of this kind of uncertainty when multiple realizations of forcings are available.

3 Assessing the ISMIP6 ensemble through the probabilistic lens

The response of an ice sheet to a given forcing may be estimated with earth system models directly. At present, however, earth system models with built-in interactive ice sheets remain in their infancy (Vizcaino2014) and are not yet able to adequately resolve ice sheet processes such as grounding line migration at the necessary resolution, requiring intermediate steps. A common approach, pursued by Goelzer et al. (2020), involves general circulation models to calculate how the global climate responds to a given forcing , regional climate models to downscale the global climate response to the ice sheet scale, and process models and parameterizations (e.g., surface energy balance models, calving models, or frontal ablation models) to interface with ice sheet models. To make the daunting task of estimating ice sheet response to different forcings a tractable community effort, a certain degree of standardization, streamlining, and simplification was necessary. The ISMIP6 steering committee and its working groups prepared data sets that could be used by individual modeling groups, including but not limited to preparing oceanic (Slater et al.2019, 2020) and atmospheric (Nowicki et al.2020) boundary conditions.

Here we consider uncertainty within the ISMIP6 experimental protocols (Goelzer et al.2020; Seroussi et al.2020) through the probabilistic framework outlined above.

3.1 Incomplete consideration of uncertainty

ISMIP6 integrates over the model uncertainties, including models of ice sheet dynamics, surface mass balance, and ice front position. It does not integrate over uncertainty in parameters. We note that the ISMIP6 protocol allowed modelers to submit as many model setups as they deemed appropriate; in practice, however, each group contributed one to three setups. While it is difficult to gauge the magnitude of the resulting underestimation in predictive variance, Aschwanden et al. (2019) suggest that the parametric uncertainty for the Greenland ice sheet (inter-quartile range) at 2100 could be up to 0.3 and 12.9 cm SLE for RCP 2.6 and 8.5, respectively, which is larger than the model uncertainty suggested by the ISMIP6 experiments (0.8 and 3.4 cm SLE, respectively). If one takes the distributions of Aschwanden et al. (2019) over model parameters as representing reasonable a priori estimates of uncertainty, then the variance in ISMIP6's predictive distribution may be substantially underestimated. Similarly, aleatoric uncertainty is not considered, which has the potential to underestimate mass loss, particularly when dynamic instabilities are likely to play a large role in ice sheet evolution (Robel et al.2019). Uncertainty also emerges from model initial conditions. A strength of the ISMIP6 protocol was the independence of different modeling groups to select their model initialization protocol.

Taken together, neglecting these additional uncertainties leads to an underestimation of variability in ensemble projections. As a secondary consideration, when comparing model predictions to observations as in Fig. 1, this has the effect of ascribing misfit between modeled predictions to model uncertainty when one of these alternative sources of uncertainty may just as likely be the culprit.

3.2 A biased sample over models

The implicit hypothesis made when accounting for model error using an ensemble approach is that each model is an independent sample from P(ℳ), where the mode of P(ℳ) is the true data generating process (i.e., reality). However, the models included in the ensemble are not likely to be independent: they share many critical features like numerical methods, parameterizations, and a joint omission of potentially important physical processes that have not yet been discovered. We emphatically note that this is not a methodological criticism: it is a challenge that exists generally in science, with analogous situations arising in climate modeling (Qian et al.2016). We note also that such biases may also arise from incorrectly specified prior distributions over parameters and forcings. Nonetheless, the challenge remains real, as does its potential effect on the credibility of and uncertainty in sea level rise projections. As shown in Fig. 1, ensemble predictions of mass loss are biased relative to present observations. While the accurate reproduction of observed mass change was not a goal of Goelzer et al. (2020), the credible projection of future mass change was a stated goal. However, there is no reason to believe that a prediction that is biased now does not remain biased in its future predictions.

4 A path forward

While we do not consider AR6's use of ISMIP6 (Goelzer et al.2020; Seroussi et al.2020) and its downstream analysis (Edwards et al.2021) appropriate for use as the consensus estimate of the ice sheets' contribution to sea level over the next century, it remains a powerful blueprint for the collaborative efforts that the ice sheet modeling community is able to achieve. Building upon the multi-model ensemble approach of ISMIP6, below we offer suggestions on how to more completely account for uncertainties in intercomparison projects.

4.1 Accounting for all sources of uncertainty

While modeling efforts have captured aleatoric, parametric, initial state, and model uncertainties independently, an effective projection of the ice sheets' contribution to sea level must incorporate all of these sources simultaneously by approximately computing

(6) P ( Δ z | F ) = P ( Δ z | f , k , M ) × P ( f | F ) P ( k | M ) P ( I | M ) P ( M ) × d k d f d I d M .

To do this, we envision a multi-model ensemble similar to the effort of ISMIP6, but with each model contributing an ensemble of simulations using random parameter values drawn from consensus estimates of the uncertainties associated with parametrically defined physics (cf. Aschwanden et al.2019; Bulthuis et al.2019) and with random realizations of climate and ocean forcings (for example, an explicit ensemble of fields characterizing the probability distribution of surface and ocean temperatures) developed in collaboration with their respective modeling communities (cf. Robel et al.2019). We anticipate that such an effort will yield a distribution of sea level projections that is much broader, and thus less certain, than that presented in recent sea level rise projections (IPCC2019). However, we feel that only through modeling what may be considered “unlikely” projections will our community accurately quantify the a priori variance in the projections of numerical ice sheet models.

4.2 Conditioning simulations on observations

While accounting for all sources of uncertainty produces a prior distribution over model projections that appropriately acknowledges the current limits of our scientific understanding, it does not ameliorate the problem of inherent biases in the sampled forcings, parameters, and models. Scientists can add specificity and value to the projected distribution by taking advantage of additional information, such as the observations illustrated in Fig. 1. To address both of these problems simultaneously, we advocate for conditioning ensemble predictions on relevant observations (Aschwanden et al.2013). One way of doing this is through Bayes' theorem (often called Bayesian calibration), which states that

(7) P ( Δ z | F , O ) = P ( O | Δ z , F ) P ( Δ z | F ) P ( O | Δ z , F ) P ( Δ z | F ) d Δ z ,

where 𝒪 is a set of observations; P(O|Δz,F) is the likelihood that some simulation associated with sea level contribution prediction Δz agrees with observations; and P(Δz|F,O) is the posterior predictive distribution of sea level, which can be thought of as the prior ensemble (Eq. 6) filtered by relevant data.

All ice sheet models already perform this calibration for certain subsets of available observations, e.g., by calibration of basal traction or other parameters to yield observed surface velocity or ice geometry within observational uncertainty. This is a necessary step for models and in many cases part of model specification. However, for the purposes of projecting ice mass change, we argue that the most salient observations on which to condition the prior distribution are measurements of mass change itself (Aschwanden et al.2013; Aðalgeirsdóttir et al.2014). It is worth noting that a few works have already performed this Bayesian calibration on observations similar to the mass change observations of Fig. 1, particularly over ensembles meant to capture parametric uncertainty (Nias et al.2019; Gilford et al.2020). These studies should be used as a model for future efforts.

Conditioning on observations also requires carefully accounting for the complicated relationship between the timescales of variability in model physics, forcings, and observational uncertainty; the appropriate timescale over which simulations need to show agreement with observations is not (yet) known. The further back in time, the more spatially and temporally sparse observations become, and the larger their associated uncertainties are. Nonetheless, reliable observations of mass change are now available on the decadal timescale (see Fig. 1), reducing the likelihood of mistakenly fitting models to short-term fluctuations in weather and ocean conditions. Fortunately, the record of detailed accurate observations is growing continually, soon spanning a climatology (30 years).

By accounting for the broad range of potential a priori uncertainties in model projections and then ascribing predictive weight only to those models that demonstrate skill at reproducing observations, the path towards realistic distributions of sea level contribution over the next century is within reach. Without a large, but realistic, spread of model outcomes it might well be possible that an insufficient number of models remain after fitting to observations.

4.3 Complementary efforts

Projections made with numerical “high-fidelity” models are computationally expensive, and creating ensemble simulations of sufficient size is limited by the availability of computational resources. Training surrogate models (“emulators”) with the output of the high-fidelity models can help better characterize sea level contribution probability distribution functions.

It is worth noting that recent efforts have used ISMIP6 as a basis for further analysis, in particular by training a surrogate model on the ISMIP6 and GlacierMIP output that effectively acts as an interpolant (Edwards et al.2021). While this interpolant is an effective tool for querying the predictive distribution of sea level contribution as a function of time and climate scenario as quantified by ISMIP6, it inherits the same challenges as its antecedent, namely a lack of accounting for all uncertainty types and a mechanism for bias correction.

Modern machine learning methods show promise to complement established numerical research tools in earth system science in general and ice sheet modeling in particular (Reichstein et al.2019; Edwards et al.2019; Brinkerhoff et al.2021; Gilford et al.2020; Edwards et al.2021; Jouvet et al.2021, e.g., ). If numerical and statistical models are paired carefully and skillfully with structured expert judgment (Bamber et al.2019), credible projections of ice sheet contribution to sea level are within reach.

Figure 2Observed and simulated historical mass changes from the Antarctic ice sheet (AIS) and Greenland ice sheet (GIS) between 2000 and 2020 in gigatons (Gt) and centimeters of sea level equivalent (cm SLE). A consensus estimate of observed mass changes (The IMBIE team2018; The IMBIE Team2019) is plotted in blue along with their respective uncertainties (shaded). The ensembles of ISMIP6 (Goelzer et al.2020; Seroussi et al.2020) historical simulations and projections are plotted with dark-gray lines, and the 5th to 95th percentile mass loss rates are shown as a 90 % credibility interval with light-gray shading. Here the control simulation was removed, in alignment with Seroussi et al. (2020) and Goelzer et al. (2020).


5 Meeting the challenge

The potential economic impact of rising sea level has been estimated at over USD 1 trillion (Diaz and Moore2017), and major world economies including the US consider investing trillions of US dollars to prepare for and avert further climate change (Blumer2020). Contrasting these staggering numbers, the current funding for research related to sea level rise remains miserly, although exact numbers are not readily available and vary from country to country. During an ISMIP6 planning effort in September 2018, participating modeling groups were polled as to how many simulations they could execute in support of projecting ice sheet contributions to sea level rise. Several groups, none of whom were receiving funding to support these simulations, estimated that they could run 5–10 simulations scheduled amongst their existing commitments. Most likely as a consequence of inadequate funding, no group submitted more than three distinct model setups despite ISMIP6 encouraging participants to explore parametric uncertainties. Projecting future sea level is an effort too severely under-resourced to meet its mission, and yet millions of lives and trillions of US dollars depend on an accurate, reliable answer.

In order to assess potential impacts of sea level rise, we urgently need to be able to deliberately quantify and then systematically reduce uncertainties. Ice sheet modeling, like climate modeling before it, developed from efforts to address basic science questions. However, despite major advances in the capabilities of ice sheet models and expanding appreciation for the importance of their projections, the funding model of modest grants to address basic science and accomplish incremental model development along the way is unchanged. Internationally, governments directly support development, maintenance, and operation of the earth system models that serve as the foundation for CMIP6 (Eyring et al.2016), and this financial support has contributed to a suite of models that now convincingly reproduce observed climate variability (Jones et al.2013). It is time to similarly bring ice sheet modeling, both standalone and embedded in earth system models, to an operational level and support it with the funding the problem deserves.

The ambitious characterization of uncertainties and ensemble conditioning we propose requires a massive international and inter-agency effort in both model development and improved observational capabilities. We call for professional support for the largely computational sea level projection effort. These resources, in the form of dedicated developers and high-performance computing time, will free up scientists to continue basic science, while the global community receives the applied science (i.e., reliable sea level projections) it needs.

The past 2 decades have shown that ice sheets react to climate far more rapidly than previously thought (Rignot and Kanagaratnam2006; Joughin et al.2014). The study of glaciers and ice sheets has moved from a fringe scientific exercise to a central question of major global economic significance. In response to COVID-19, USD 18 billion flowed from the US government to fund vaccine development (Tozzi et al.2020). Appropriate resourcing is possible. While the emergent threat of sea level rise is less abrupt than that from COVID-19, a similarly serious effort is required to reduce uncertainties in sea level projections.

Data availability

We downloaded the scalar time series produced by ISMIP6 for the Antarctic and Greenland ice sheet from Zenodo with the digital object identifiers (Goelzer2020) and (Seroussi2020), respectively. The data are split into a historical period (pre-2015) and the projection period (2015–2100). For the projections we used the version where the control simulation was not removed; however a version where the control simulation was removed is shown below (Fig. 2). Removal of the control simulation is intended to account for unforced model drift and mass loss committed as a result of non-equilibrium ice sheet conditions at the start of the simulations. Committed sea level rise is estimated to add an additional 6 mm to simulated sea level rise by 2100 (Price et al.2011; Goelzer et al.2020) and thus has little impact on the low bias of recent, simulated mass loss. For observations, we used a multi-method consensus estimate (The IMBIE team2018; The IMBIE Team2019). To analyze the data, we created a Jupyter notebook (, Aschwanden et al.2021).

Author contributions

AA, TCB, DJB, and MT jointly conceived the study and contributed equally to the writing. AA performed the data analysis with input from the other authors.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to thank editor Olaf Eisen, Alexander Robel, Nicolas Jourdain, and the ISMIP6 steering committee for their comments, which helped to improve this paper.

Review statement

This paper was edited by Olaf Eisen and reviewed by Alexander Robel and Nicolas Jourdain.


Amaral, T., Bartholomaus, T. C., and Enderlin, E. M.: Evaluation of Iceberg Calving Models Against Observations From Greenland Outlet Glaciers, J. Geophys. Res.-Earth, 125, 1–29,, 2020. a

Amundson, J. M., Fahnestock, M., Truffer, M., Brown, J., Lüthi, M. P., and Motyka, R. J.: Ice mélange dynamics and implications for terminus stability, Jakobshavn Isbræ, Greenland, J. Geophys. Res., 115, F01005,, 2010. a

Aschwanden, A., Aðalgeirsdóttir, G., and Khroulev, C.: Hindcasting to measure ice sheet model sensitivity to initial states, The Cryosphere, 7, 1083–1093,, 2013. a, b, c

Aschwanden, A., Fahnestock, M. A., Truffer, M., Brinkerhoff, D. J., Hock, R., Khroulev, C., Mottram, R., and Khan, S. A.: Contribution of the Greenland Ice Sheet to sea level over the next millennium, Science Advances, 5, eaav9396,, 2019. a, b, c, d

Aschwanden, A., Brinkerhoff, D. C., Bartholomaus, T. C., and Truffer, M.: Scripts for “Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level”, v0.1, Zenodo [code],, 2021. a

Aðalgeirsdóttir, G., Aschwanden, A., Khroulev, C., Boberg, F., Mottram, R., and Lucas-Picher, P.: Role of model initialization for projections of 21st-century Greenland ice sheet mass loss, J. Glaciol., 60, 782–794,, 2014. a, b

Bamber, J. L., Oppenheimer, M., Kopp, R. E., Aspinall, W. P., and Cooke, R. M.: Ice sheet contributions to future sea-level rise from structured expert judgment, P. Natl. Acad. Sci. USA, 166, 11195–11200,, 2019. a, b

Bindschadler, R. A., Nowicki, S., Abe-Ouchi, A., Aschwanden, A., Choi, H., Fastook, J., Granzow, G., Greve, R., Gutowski, G., Herzfeld, U., Jackson, C., Johnson, J., Khroulev, C., Levermann, A., Lipscomb, W. H., Martin, M. A., Morlighem, M., Parizek, B. R., Pollard, D., Price, S. F., Ren, D., Saito, F., Sato, T., Seddik, H., Seroussi, H., Takahashi, K., Walker, R., and Wang, W. L.: Ice-sheet model sensitivities to environmental forcing and their use in projecting future sea level (the SeaRISE project), J. Glaciol., 59, 195–224,, 2013. a

Blumer, B.: To Cut Emissions to Zero, U.S. Needs to Make Big Changes in Next 10 Years, available at: (last access: 27 April 2021), 2020. a

Brinkerhoff, D., Aschwanden, A., and Fahnestock, M.: Constraining subglacial processes from surface velocity observations using surrogate-based Bayesian inference, J. Glaciol., 67, 385–403,, 2021. a

Bueler, E. and van Pelt, W.: Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model version 0.6, Geosci. Model Dev., 8, 1613–1635,, 2015. a

Bulthuis, K., Arnst, M., Sun, S., and Pattyn, F.: Uncertainty quantification of the multi-centennial response of the Antarctic ice sheet to climate change, The Cryosphere, 13, 1349–1380,, 2019. a

Chen, X., Zhang, X., Church, J. A., Watson, C. S., King, M. A., Monselesan, D., Legresy, B., and Harig, C.: The increasing rate of global mean sea-level rise during 1993–2014, Nat. Clim. Change, 7, 492–495,, 2017. a

DeConto, R. M., Pollard, D., Alley, R. B., Velicogna, I., Gasson, E., Gomez, N., Sadai, S., Condron, A., Gilford, D. M., Ashe, E. L., Kopp, R. E., Li, D., and Dutton, A.: The Paris Climate Agreement and future sea-level rise from Antarctica, Nature, 593, 83–89,, 2021. a

Diaz, D. and Moore, F.: Quantifying the economic risks of climate change, Nat. Clim. Change, 7, 774–782,, 2017. a

Edwards, T. L., Brandon, M. A., Durand, G., Edwards, N. R., Golledge, N. R., Holden, P. B., Nias, I. J., Payne, A. J., Ritz, C., and Wernecke, A.: Revisiting Antarctic ice loss due to marine ice-cliff instability, Nature, 566, 58–64,, 2019. a

Edwards, T. L., Nowicki, S., Marzeion, B., Hock, R., Goelzer, H., Seroussi, H., Jourdain, N. C., Slater, D. A., Turner, F. E., Smith, C. J., McKenna, C. M., Simon, E., Abe-Ouchi, A., Gregory, J. M., Larour, E., Lipscomb, W. H., Payne, A. J., Shepherd, A., Agosta, C., Alexander, P., Albrecht, T., Anderson, B., Asay-Davis, X., Aschwanden, A., Barthel, A., Bliss, A., Calov, R., Chambers, C., Champollion, N., Choi, Y., Cullather, R., Cuzzone, J., Dumas, C., Felikson, D., Fettweis, X., Fujita, K., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huss, M., Huybrechts, P., Immerzeel, W., Kleiner, T., Kraaijenbrink, P., Le clec’h, S., Lee, V., Leguy, G. R., Little, C. M., Lowry, D. P., Malles, J.-H., Martin, D. F., Maussion, F., Morlighem, M., O’Neill, J. F., Nias, I., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Radić, V., Reese, R., Rounce, D. R., Rückamp, M., Sakai, A., Shafer, C., Schlegel, N.-J., Shannon, S., Smith, R. S., Straneo, F., Sun, S., Tarasov, L., Trusel, L. D., Van Breedam, J., van de Wal, R., van den Broeke, M., Winkelmann, R., Zekollari, H., Zhao, C., Zhang, T., and Zwinger, T.: Projected land ice contributions to twenty-first-century sea level rise, Nature, 593, 74–82,, 2021. a, b, c, d, e

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958,, 2016. a

Gilford, D. M., Ashe, E. L., DeConto, R. M., Kopp, R. E., Pollard, D., and Rovere, A.: Could the Last Interglacial Constrain Projections of Future Antarctic Ice Mass Loss and Sea-Level Rise?, J. Geophys. Res.-Earth, 125, e2019JF005418,, 2020. a, b

Goelzer, H.: Dataset for “The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6”, v7.0.1, Zenodo [data set],, 2020. a

Goelzer, H., Nowicki, S., Payne, A., Larour, E., Seroussi, H., Lipscomb, W. H., Gregory, J., Abe-Ouchi, A., Shepherd, A., Simon, E., Agosta, C., Alexander, P., Aschwanden, A., Barthel, A., Calov, R., Chambers, C., Choi, Y., Cuzzone, J., Dumas, C., Edwards, T., Felikson, D., Fettweis, X., Golledge, N. R., Greve, R., Humbert, A., Huybrechts, P., Le clec'h, S., Lee, V., Leguy, G., Little, C., Lowry, D. P., Morlighem, M., Nias, I., Quiquet, A., Rückamp, M., Schlegel, N.-J., Slater, D. A., Smith, R. S., Straneo, F., Tarasov, L., van de Wal, R., and van den Broeke, M.: The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6, The Cryosphere, 14, 3071–3096,, 2020. a, b, c, d, e, f, g, h, i, j

Hill, E. A., Rosier, S. H. R., Gudmundsson, G. H., and Collins, M.: Quantifying the potential future contribution to global mean sea level from the Filchner–Ronne basin, Antarctica, The Cryosphere, 15, 4675–4702,, 2021. a

Hinkel, J., Lincke, D., Vafeidis, A. T., Perrette, M., Nicholls, R. J., Tol, R. S. J., Marzeion, B., Fettweis, X., Ionescu, C., and Levermann, A.: Coastal flood damage and adaptation costs under 21st century sea-level rise, P. Natl. Acad. Sci. USA, 111, 3292–3297,, 2014. a

Hoffman, M. J., Asay-Davis, X., Price, S. F., Fyke, J., and Perego, M.: Effect of Subshelf Melt Variability on Sea Level Rise Contribution From Thwaites Glacier, Antarctica, J. Geophys. Res.-Earth, 124, 2798–2822,, 2019. a

IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2013. a

IPCC: Special Report: The Ocean and Cryosphere in a Changing Climate, in press, available at: (last access: 30 November 2021), 2019. a

Jones, G. S., Stott, P. A., and Christidis, N.: Attribution of observed historical near-surface temperature variations to anthropogenic and natural causes using CMIP5 simulations, J. Geophys. Res.-Atmos., 118, 4001–4024,, 2013. a

Joughin, I., Smith, B. E., and Medley, B.: Marine Ice Sheet Collapse Potentially Underway for the Thwaites Glacier Basin, West Antarctica, Science, 344, 735–738,, 2014. a

Joughin, I., Shean, D. E., Smith, B. E., and Floricioiu, D.: A decade of variability on Jakobshavn Isbræ: ocean temperatures pace speed through influence on mélange rigidity, The Cryosphere, 14, 211–227,, 2020. a

Jouvet, G., Cordonnier, G., Kim, B., Lüthi, M., Vieli, A., and Aschwande, A.: Deep learning speeds up ice flow modelling by several orders of magnitude, J. Glaciol., in press, 2021. a

Kopp, R. E., Kemp, A. C., Bittermann, K., Horton, B. P., Donnelly, J. P., Gehrels, W. R., Hay, C. C., Mitrovica, J. X., Morrow, E. D., and Rahmstorf, S.: Temperature-driven global sea-level variability in the Common Era, P. Natl. Acad. Sci. USA, 113, E1434–E1441,, 2016. a

Levermann, A., Winkelmann, R., Albrecht, T., Goelzer, H., Golledge, N. R., Greve, R., Huybrechts, P., Jordan, J., Leguy, G., Martin, D., Morlighem, M., Pattyn, F., Pollard, D., Quiquet, A., Rodehacke, C., Seroussi, H., Sutter, J., Zhang, T., Van Breedam, J., Calov, R., DeConto, R., Dumas, C., Garbe, J., Gudmundsson, G. H., Hoffman, M. J., Humbert, A., Kleiner, T., Lipscomb, W. H., Meinshausen, M., Ng, E., Nowicki, S. M. J., Perego, M., Price, S. F., Saito, F., Schlegel, N.-J., Sun, S., and van de Wal, R. S. W.: Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2), Earth Syst. Dynam., 11, 35–76,, 2020. a, b

Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B. (Eds.): Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge Univerity Press, Cambridge, 2021. a

Mercer, J. H.: West Antarctic ice sheet and CO2 greenhouse effect: a threat of disaster, Nature, 271, 321–325,, 1978. a

Moon, T., Scambos, T., Abdalati, W., Ahlstrøm, A. P., Bindschadler, R., Gambill, J., Heimbach, P., Hock, R., Langley, K., Miller, I., and Truffer, M.: Ending a Sea of Confusion: Insights and Opportunities in Sea-Level Change Communication, Environment: Science and Policy for Sustainable Development, 62, 4–15,, 2020. a

NIAS, I. J., CORNFORD, S. L., and PAYNE, A. J.: Contrasting the modelled sensitivity of the Amundsen Sea Embayment ice streams, J. Glaciol., 62, 552–562,, 2016. a

Nias, I. J., Cornford, S. L., Edwards, T. L., Gourmelen, N., and Payne, A. J.: Assessing Uncertainty in the Dynamical Ice Response to Ocean Warming in the Amundsen Sea Embayment, West Antarctica, Geophys. Res. Lett., 46, 11253–11260,, 2019. a

Nowicki, S., Goelzer, H., Seroussi, H., Payne, A. J., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Alexander, P., Asay-Davis, X. S., Barthel, A., Bracegirdle, T. J., Cullather, R., Felikson, D., Fettweis, X., Gregory, J. M., Hattermann, T., Jourdain, N. C., Kuipers Munneke, P., Larour, E., Little, C. M., Morlighem, M., Nias, I., Shepherd, A., Simon, E., Slater, D., Smith, R. S., Straneo, F., Trusel, L. D., van den Broeke, M. R., and van de Wal, R.: Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models, The Cryosphere, 14, 2331–2368,, 2020. a, b, c

Nowicki, S. M. J., Payne, A., Larour, E., Seroussi, H., Goelzer, H., Lipscomb, W., Gregory, J., Abe-Ouchi, A., and Shepherd, A.: Ice Sheet Model Intercomparison Project (ISMIP6) contribution to CMIP6, Geosci. Model Dev., 9, 4521–4545,, 2016. a, b

Payne, A. J., Nowicki, S., Abe‐Ouchi, A., Agosta, C., Alexander, P., Albrecht, T., Asay‐Davis, X., Aschwanden, A., Barthel, A., Bracegirdle, T. J., Calov, R., Chambers, C., Choi, Y., Cullather, R., Cuzzone, J., Dumas, C., Edwards, T. L., Felikson, D., Fettweis, X., Galton‐Fenzi, B. K., Goelzer, H., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Munneke, P. K., Larour, E., Le clec'h, S., Lee, V., Leguy, G., Lipscomb, W. H., Little, C. M., Lowry, D. P., Morlighem, M., Nias, I., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Rückamp, M., Schlegel, N., Seroussi, H., Shepherd, A., Simon, E., Slater, D., Smith, R. S., Straneo, F., Sun, S., Tarasov, L., Trusel, L. D., Van Breedam, J., Wal, R., Broeke, M., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: Future Sea Level Change Under Coupled Model Intercomparison Project Phase 5 and Phase 6 Scenarios From the Greenland and Antarctic Ice Sheets, Geophys. Res. Lett., 48, 1–13,, 2021. a

Price, S. F., Payne, A. J., Howat, I. M., and Smith, B. E.: Committed sea-level rise for the next century from Greenland ice sheet dynamics during the past decade, P. Natl. Acad. Sci. USA, 108, 8978–8983,, 2011. a

Qian, Y., Jackson, C., Giorgi, F., Booth, B., Duan, Q., Forest, C., Higdon, D., Hou, Z. J., and Huerta, G.: Uncertainty Quantification in Climate Modeling and Projection, B. Am. Meteorol. Soc., 97, 821–824,, 2016. a

Reichstein, M., Camps-Valls, G., Stevens, B., Jung, M., Denzler, J., Carvalhais, N., and Prabhat: Deep learning and process understanding for data-driven Earth system science, Nature, 566, 195–204,, 2019. a

Rietbroek, R., Brunnabend, S.-E., Kusche, J., Schröter, J., and Dahle, C.: Revisiting the contemporary sea-level budget on global and regional scales., P. Natl. Acad. Sci. USA, 113, 1504–1509,, 2016. a

Rignot, E. and Kanagaratnam, P.: Changes in the Velocity Structure of the Greenland Ice Sheet, Science, 311, 986–990,, 2006. a

Robel, A. A., Seroussi, H., and Roe, G. H.: Marine ice sheet instability amplifies and skews uncertainty in projections of future sea-level rise, P. Natl. Acad. Sci. USA, 116, 14887–14892,, 2019. a, b, c

Sadai, S., Condron, A., DeConto, R., and Pollard, D.: Future climate response to Antarctic Ice Sheet melt caused by anthropogenic warming, Science Advances, 6, 1–9,, 2020. a

Schlegel, N.-J., Seroussi, H., Schodlok, M. P., Larour, E. Y., Boening, C., Limonadi, D., Watkins, M. M., Morlighem, M., and van den Broeke, M. R.: Exploration of Antarctic Ice Sheet 100-year contribution to sea level rise and associated model uncertainties using the ISSM framework, The Cryosphere, 12, 3511–3534,, 2018. a

Seroussi, H.: Scripts for “ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century”, v1, Zenodo [data set],, 2020. a

Seroussi, H., Nowicki, S., Simon, E., Abe-Ouchi, A., Albrecht, T., Brondex, J., Cornford, S., Dumas, C., Gillet-Chaulet, F., Goelzer, H., Golledge, N. R., Gregory, J. M., Greve, R., Hoffman, M. J., Humbert, A., Huybrechts, P., Kleiner, T., Larour, E., Leguy, G., Lipscomb, W. H., Lowry, D., Mengel, M., Morlighem, M., Pattyn, F., Payne, A. J., Pollard, D., Price, S. F., Quiquet, A., Reerink, T. J., Reese, R., Rodehacke, C. B., Schlegel, N.-J., Shepherd, A., Sun, S., Sutter, J., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., and Zhang, T.: initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6, The Cryosphere, 13, 1441–1471,, 2019. a

Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070,, 2020. a, b, c, d, e, f

Slater, D. A., Straneo, F., Felikson, D., Little, C. M., Goelzer, H., Fettweis, X., and Holte, J.: Estimating Greenland tidewater glacier retreat driven by submarine melting, The Cryosphere, 13, 2489–2509,, 2019. a

Slater, D. A., Felikson, D., Straneo, F., Goelzer, H., Little, C. M., Morlighem, M., Fettweis, X., and Nowicki, S.: Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution, The Cryosphere, 14, 985–1008,, 2020. a

The IMBIE team: Mass balance of the Antarctic Ice Sheet from 1992 to 2017, Nature, 558, 219–222,, 2018. a, b, c

The IMBIE Team: Mass balance of the Greenland Ice Sheet from 1992 to 2018, Nature, 579, 233–239,, 2019. a, b, c, d, e

Tozzi, J., Griffin, R., and Stein, S.: Trump Administration Dips Into Protective Gear, CDC Funds to Fund Vaccine Push, available at: (last access: 27 April 2021), 2020. a

Vaughan, D. G. and Arthern, R. J.: Why Is It Hard to Predict the Future of Ice Sheets?, Science, 315, 1503–1504,, 2007.  a

Vizcaino, M.: Ice sheets as interactive components of Earth System Models: progress and challenges, Wiley Interdisciplinary Reviews: Climate Change, 1–12,, 2014. a

Short summary
Estimating how much ice loss from Greenland and Antarctica will contribute to sea level rise is of critical societal importance. However, our analysis shows that recent efforts are not trustworthy because the models fail at reproducing contemporary ice melt. Here we present a roadmap towards making more credible estimates of ice sheet melt.