Global mean sea-level rise (SLR), which during the last quarter century has occurred at an accelerating rate (1), averaging about +3 mm⋅y−1, threatens coastal communities and ecosystems worldwide. Adaptation measures accounting for the changing hazard, including building or raising permanent or movable structures such as surge barriers and sea walls, enhancing nature-based defenses such as wetlands, and selective retreat of populations and facilities from areas threatened by episodic flooding or permanent inundation, are being planned or implemented in several countries. Risk assessment for such adaptation efforts requires projections of future SLR, including careful characterization and evaluation of uncertainties (2) and regional projections that account for ocean dynamics, gravitational and rotational effects, and vertical land motion (3). During the nearly 40 y since the first modern, scientific assessments of SLR, understanding of the various causes of this rise has advanced substantially. Improvements during the past decade include closing the historic sea-level budget, attributing global mean SLR to human activities, confirming acceleration of SLR since the nineteenth century and during the satellite altimetry era, and developing analytical frameworks for estimating regional and local mean sea level and extreme water level changes. Nonetheless, long-term SLR projections remain acutely uncertain, in large part because of inadequate understanding of the potential future behaviors of the Greenland and Antarctic ice sheets and their responses to future global climate change. This limitation is especially troubling, given that the ice sheet influence on SLR has been increasing since the 1990s (4) and has overtaken mountain glaciers to become the largest barystatic (mass) contribution to SLR (5). In addition, for any given future climate scenario, the ice sheets constitute the component with the largest uncertainties by a substantial margin, especially beyond 2050 (6).
Advances since the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (7) include improved process understanding and representation in deterministic ice sheet models (8, 9), probabilistic projections calibrated against these models and the observational record (10), and new semiempirical models, based on the historical relationship between temperature and sea-level changes. Each of these approaches has limitations that stem from factors including poorly understood processes, poorly constrained boundary conditions, and a short (∼25 y) satellite observation record of ice sheets that does not capture the time scales of internal variability in the ice sheet climate system. As a consequence, it is unclear to what extent recent observed ice sheet changes (11) are a result of internal variability (ice sheet weather) or external forcing (ice sheet climate).
Where other methods are intractable for scientific or practical reasons, structured expert judgement (SEJ), using calibrated expert responses, provides a formal approach for estimating uncertain quantities according to current scientific understanding. It has been used in a wide range of applications, including natural and anthropogenic hazards such as earthquakes, volcanic eruptions, vector-borne disease spread, and nuclear waste security (12). That said, it should not be regarded as a substitute for fundamental research into driving processes, but instead as a source of complementary insights into the current state of knowledge and, in particular, the extent of the uncertainties (12). An SEJ study conducted in advance of the AR5 (13) (hereafter BA13) provided valuable insights into the uncertainties around ice sheet projections, as assessed at that time.
Since then, regional- and continental-scale, process-based modeling of ice sheets has advanced substantially (8, 9, 14–16), with the inclusion of new positive feedbacks that could potentially accelerate mass loss, and negative feedbacks that could potentially slow it. These include solid Earth and gravitational processes (17, 18), Antarctic marine ice cliff instability (19), and the influences of organic and inorganic impurities on the albedo of the Greenland Ice Sheet (20). The importance of these feedbacks is an area of continuing research. Therefore, alternative approaches must be exploited to assess future SLR and, critically, its associated uncertainties (21), to serve the more immediate needs of the science and policy communities.
Here, we report the findings of an SEJ exercise undertaken in 2018 via separate, 2-d workshops held in the United States and United Kingdom, involving 22 experts (hereafter SEJ2018). Details of how experts were selected are provided in SI Appendix, Note 1. The questions and format of the workshops were identical, so that the findings could be combined using an impartial weighting approach (Methods). SEJ (as opposed to other types of expert elicitation) weights each expert using objective estimates of their statistical accuracy and informativeness (22), determined using experts’ uncertainty evaluations over a set of seed questions from their field with ascertainable values (Methods). The approach is analogous to weighting climate models based on their skill in capturing a relevant property, such as the regional 20th century surface air temperature record (23). In SEJ, the synthetic expert (i.e., the performance weighted [PW] combination of all of the experts’ judgments) in general outperforms an equal weights (EW) combination in terms of statistical accuracy and informativeness, as illustrated in SI Appendix, Fig. S3. The approach is particularly effective at identifying those experts who are able to quantify their uncertainties with high statistical accuracy for specified problems rather than, for example, experts with restricted domains of knowledge or even high scientific reputation (12).
The participating experts quantified their uncertainties for three physical processes relevant to ice sheet mass balance: accumulation, discharge, and surface runoff. They did this for each of the Greenland, West Antarctic, and East Antarctic ice sheets (GrIS, WAIS, and EAIS, respectively), and for two schematic temperature change scenarios. The first temperature trajectory (denoted L) stabilized in 2100 at +2 °C above preindustrial global mean surface air temperature (defined as the average for 1850–1900), and the second (denoted H) stabilized at +5 °C (SI Appendix, Fig. S1). The experts generated values for four dates: 2050, 2100, 2200, and 2300. Experts also quantified the dependence between accumulation, runoff, and discharge within each of the three ice sheets, and between each ice sheet for discharge only, for the H scenario in 2100. We used temperature trajectories rather than emissions scenarios to isolate the experts’ judgements about the relationship between global mean surface air temperature change and ice sheet changes from judgements about climate sensitivity.
An important and unique element of SEJ2018 was the elicitation of intra- and inter-ice sheet dependencies (SI Appendix, Note 1.5). Two features of dependence were elicited: a central dependence and an upper tail dependence. The former captures the probability that one variable exceeds its median given that the other variable exceeds its median, whereas the latter captures the probability that one variable exceeds its 95th percentile given that the other exceeds its 95th percentile. It is well known that these two types of dependence are, in general, markedly different, a property that is not captured by the usual Gaussian dependence model. The latter always imposes tail independence, regardless of the degree of central dependence, and can produce large errors when applied inappropriately (24). For example, if GrIS discharge exceeds its 95th percentile, what is the probability that runoff will also exceed its 95th percentile? This probability may be substantially higher than the independent probability of 5%, and ignoring tail dependence may lead to underestimating the probability of high SLR contributions. On the basis of each expert’s responses, a joint distribution was constructed to capture the dependencies among runoff, accumulation, and discharge for GrIS, WAIS, and EAIS, with dependence structures chosen, per expert, to capture central and tail dependences (Methods and SI Appendix, Note 1.5). In BA13, heuristic dependency values were applied on the basis of simple assumptions about the response of processes to a common forcing.
To help interpret the findings, experts were also asked to provide qualitative and rank-order information on what they regard to be the leading processes that could influence ice dynamics and surface mass balance (snowfall minus ablation); henceforth, this is termed the descriptive rationale. Further details can be found in the SI Appendix. The combined sea-level contribution from all processes and ice sheets was determined assuming either independence or dependence. Here, we focus on the findings with dependence; we examine the effect of the elicited dependencies and the approach taken in SI Appendix, Note 1.5.
The ice sheet contributions were expressed as anomalies from the 2000–2010 mean states, which were predefined (SI Appendix, Table S7). The baseline sea-level contribution for this period was prescribed as 0.76 mm⋅y−1 (0.56, 0.20, and 0.00 mm⋅y−1 for GrIS, WAIS, and EAIS, respectively) and has been added to the elicited values discussed here. This is close to an observationally derived value of 0.79 mm⋅y−1 for the same period, which was published subsequently to the SEJ workshops (4).