Future Sea Level Contribution Research Articles

Cloud- and ice-albedo feedbacks drive greater Greenland Ice Sheet sensitivity to warming in CMIP6 than in CMIP5

Abstract. The Greenland Ice Sheet (GrIS) has been losing mass since the 1990s as a direct consequence of rising temperatures and has been projected to continue to lose mass at an accelerating pace throughout the 21st century, making it one of the largest contributors to future sea-level rise. The latest Coupled Model Intercomparison Project Phase 6 (CMIP6) models produce a greater Arctic amplification signal and therefore also a notably larger mass loss from the GrIS when compared to the older CMIP5 projections, despite similar forcing levels from greenhouse gas emissions. However, it is also argued that the strength of regional factors, such as melt–albedo feedbacks and cloud-related feedbacks, will partly impact future melt and sea-level rise contribution, yet little is known about the role of these regional factors in producing differences in GrIS surface melt projections between CMIP6 and CMIP5. In this study, we use high-resolution (15 km) regional climate model simulations over the GrIS performed using the Modèle Atmosphérique Régional (MAR) to physically downscale six CMIP5 Representative Concentration Pathway (RCP) 8.5 and five CMIP6 Shared Socioeconomic Pathway (SSP) 5-8.5 extreme high-emission-scenario simulations. Here, we show a greater annual mass loss from the GrIS at the end of the 21st century but also for a given temperature increase over the GrIS, when comparing CMIP6 to CMIP5. We find a greater sensitivity of Greenland surface mass loss in CMIP6 centred around summer and autumn, yet the difference in mass loss is the largest during autumn with a reduction of 27.7 ± 9.5 Gt per season for a regional warming of +6.7 ∘C and 24.6 Gt per season more mass loss than in CMIP5 RCP8.5 simulations for the same warming. Assessment of the surface energy budget and cloud-related feedbacks suggests a reduction in high clouds during summer and autumn – despite enhanced cloud optical depth during autumn – to be the main driver of the additional energy reaching the surface, subsequently leading to enhanced surface melt and mass loss in CMIP6 compared to CMIP5. Our analysis highlights that Greenland is losing more mass in CMIP6 due to two factors: (1) a (known) greater sensitivity to greenhouse gas emissions and therefore warmer temperatures and (2) previously unnotified cloud-related surface energy budget changes that enhance the GrIS sensitivity to warming.

The evolution of future Antarctic surface melt using PISM-dEBM-simple

Abstract. It is virtually certain that Antarctica's contribution to sea-level rise will increase with future warming, although competing mass balance processes hamper accurate quantification of the exact magnitudes. Today, ocean-induced melting underneath the floating ice shelves dominates mass losses, but melting at the surface will gain importance as global warming continues. Meltwater at the ice surface has crucial implications for the ice sheet's stability, as it increases the risk of hydrofracturing and ice-shelf collapse that could cause enhanced glacier outflow into the ocean. Simultaneously, positive feedbacks between ice and atmosphere can accelerate mass losses and increase the ice sheet's sensitivity to warming. However, due to long response times, it may take hundreds to thousands of years until the ice sheet fully adjusts to the environmental changes. Therefore, ice-sheet model simulations must be computationally fast and capture the relevant feedbacks, including the ones at the ice–atmosphere interface. Here we use the novel surface melt module dEBM-simple (a slightly modified version of the “simple” diurnal Energy Balance Model) coupled to the Parallel Ice Sheet Model (PISM, together referred to as PISM-dEBM-simple) to estimate the impact of 21st-century atmospheric warming on Antarctic surface melt and ice dynamics. As an enhancement compared to the widely adopted positive degree-day (PDD) scheme, dEBM-simple includes an implicit diurnal cycle and computes melt not only from the temperature, but also from the influence of solar radiation and changes in ice albedo, thus accounting for the melt–albedo feedback. We calibrate PISM-dEBM-simple to reproduce historical and present-day Antarctic surface melt rates given by the regional atmospheric climate model RACMO2.3p2 and use the calibrated model to assess the range of possible future surface melt trajectories under Shared Socioeconomic Pathway SSP5-8.5 warming projections until the year 2100. To investigate the committed impacts of the enhanced surface melting on the ice-sheet dynamics, we extend the simulations under fixed climatological conditions until the ice sheet has reached a state close to equilibrium with its environment. Our findings reveal a substantial surface-melt-induced speed-up in ice flow associated with large-scale elevation reductions in sensitive ice-sheet regions, underscoring the critical role of self-reinforcing ice-sheet–atmosphere feedbacks in future mass losses and sea-level contribution from the Antarctic Ice Sheet on centennial to millennial timescales.

A thicker Antarctic ice stream during the mid-Pliocene warm period

Ice streams regulate most ice mass loss in Antarctica. Determining ice stream response to warmer conditions during the Pliocene could provide insights into their future behaviour, but this is hindered by a poor representation of subglacial topography in ice-sheet models. We address this limitation using a high-resolution model for Dronning Maud Land (East Antarctica). We show that contrary to dynamic thinning of the region’s ice streams following ice-shelf collapse, the largest ice stream, Jutulstraumen, thickens by 700 m despite lying on a retrograde bed slope. We attribute this counterintuitive thickening to a shallower Pliocene subglacial topography and inherent high lateral stresses at its flux gate. These conditions constrict ice drainage and, combined with increased snowfall, allow ice accumulation upstream. Similar stress balances and increased precipitation projections occur across 27% of present-day East Antarctica, and understanding how lateral stresses regulate ice-stream discharge is necessary for accurately assessing Antarctica’s future sea-level rise contribution.

Effects of extreme melt events on ice flow and sea level rise of the Greenland Ice Sheet

Abstract. Over the past decade, Greenland has experienced several extreme melt events, the most pronounced ones in the years 2010, 2012 and 2019. With progressing climate change, such extreme melt events can be expected to occur more frequently and potentially become more severe and persistent. So far, however, projections of ice loss and sea level change from Greenland typically rely on scenarios which only take gradual changes in the climate into account. Using the Parallel Ice Sheet Model (PISM), we investigate the effect of extreme melt events on the overall mass balance of the Greenland Ice Sheet and the changes in ice flow, invoked by the altered surface topography. As a first constraint, this study estimates the overall effect of extreme melt events on the cumulative mass loss of the Greenland Ice Sheet. We find that the sea level contribution from Greenland might increase by 2 to 45 cm (0.2 % to 14 %) by the year 2300 if extreme events occur more frequently in the future under a Representative Concentration Pathway 8.5 (RCP8.5) scenario, and the ice sheet area might be reduced by an additional 6000 to 26 000 km2 by 2300 in comparison to future warming scenarios without extremes. In conclusion, projecting the future sea level contribution from the Greenland Ice Sheet requires consideration of the changes in both the frequency and intensity of extreme events. It is crucial to individually address these extremes at a monthly resolution as temperature forcing with the same excess temperature but evenly distributed over longer timescales (e.g., seasonal) leads to less sea level rise than for the simulations of the resolved extremes. Extremes lead to additional mass loss and thinning. This, in turn, reduces the driving stress and surface velocities, ultimately dampening the ice loss attributed to ice flow and discharge. Overall, we find that the surface elevation feedback largely amplifies melting for scenarios with and without extremes, with additional mass loss attributed to this feedback having the greatest impact on projected sea level.

Damage detection on antarctic ice shelves using the normalised radon transform

Areas of structural damage mechanically weaken Antarctic ice shelves. This potentially preconditions ice shelves for disintegration and enhanced grounding line retreat. The development of damage and its feedback on marine ice sheet dynamics has been identified as key to future ice shelf stability and sea level contributions from Antarctica. However, it is one of the least understood processes that impact ice shelf instability since quantifying damage efficiently and accurately is a challenging task. Challenges relate to the complex surface of Antarctica, variations in viewing-illumination geometry, snow or cloud cover and variable signal-to-noise levels in satellite imagery. Therefore, automated damage assessment approaches require careful pre- and post-processing, lacking the option to be applied to wider spatiotemporal domains. Simultaneously, studies that use manual mapping are usually limited due to the effort required for extensive mapping, which either results in a limited spatial domain or the use of low resolution data.This study proposes the NormalisEd Radon transform Damage detection (NeRD) method to detect damage features and their orientations from multi-source satellite imagery. NeRD performs robust, high resolution, large-scale damage assessments. NeRD is applied to the ice shelves in the Amundsen Sea Embayment (ASE) and validated with both manually labelled and existing fracture maps. Validation shows that NeRD detects damage with high recall and provides an accurate physical representation of multi-scale damage features and their orientation. Sensitivity analyses indicate NeRD is robust to different resolution parameter settings. NeRD consistently detects damage for different data sources ranging from optical Landsat 7/8 and Sentinel-2 optical to Synthetic Aperture Radar Sentinel-1 data. Therefore, NeRD paves the way for synergistic multi-source damage detection that overcomes remaining limitations from individual sources. Results show varying damage patterns on the ice shelves in the ASE area in austral summer 2020–2021, with most damage located on the Pine Island, Crosson and Thwaites ice shelves. We show a damage increase on the Pine Island ice shelf between 2013–2019, and display advection and rotation of crevasses. The detected damage orientation can provide insight in the type of crevasse opening mode and the development of damage over time. The damage maps produced with NeRD can help evaluate ice sheet models or machine learning approaches, improving our understanding of damage evolution.

Stabilizing effect of bedrock uplift on retreat of Thwaites Glacier, Antarctica, at centennial timescales

Viscoelastic rebound of the solid Earth upon the removal of ice loads has the potential to inhibit marine ice sheet instability, thereby forestalling ice-sheet retreat and global mean sea-level rise. The timescale over which the solid Earth - ice sheet system responds to changes in ice thickness and bedrock topography places a strong control on the spatiotemporal influence of this negative feedback mechanism. In this study, we assess the impact of solid-earth rheological structure on model projections of the retreat of Thwaites Glacier, West Antarctica, and the concomitant sea-level rise by coupling the dynamic ice sheet model MALI to a regional glacial isostatic adjustment (GIA) model. We test the sensitivity of model projections of ice-sheet retreat and associated sea-level rise across a range of four different solid-earth rheologies, forced by standard ISMIP6 ocean and atmospheric datasets for the RCP8.5 climate scenario. These model parameters are applied to 500-year, coupled ice-sheet - GIA simulations. For the mantle viscosity best supported by observations, the negative GIA feedback leads to a reduction in mass loss that remains above 20% after about a hundred years. Mass-loss reduction peaks at 50% around 2300, which is when a control simulation without GIA experiences its maximum rate of retreat. For a weaker solid-earth rheology that is unlikely but compatible with observational uncertainty, mass loss reduction remains above 50% after 2150. At 2100, mass loss reduction is 10% for the best-fit rheology and 25% for the weakest rheology. At the same time, we estimate that water expulsion from the rebounding solid Earth beneath the ocean near Thwaites Glacier may increase sea-level rise by up to 20% at five centuries. Additionally, the reduction in ice-sheet retreat caused by GIA is substantially reduced under stronger climate forcings, suggesting that the stabilizing feedback of GIA will also be an indirect function of emissions scenario. We hypothesize that feedbacks between the solid Earth - ice sheet system are controlled by a competition between the spatial extent and timescale of bedrock uplift relative to the rate of grounded ice retreat away from the region of most rapid unloading. Although uncertainty in solid-earth rheology leads to large uncertainty in future sea-level rise contribution from Thwaites Glacier, under all plausible parameters the GIA effects are too large to be ignored for future projections of Thwaites Glacier of more than a century.

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

Abstract. 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.

Two decades of dynamic change and progressive destabilization on the Thwaites Eastern Ice Shelf

Abstract. The Thwaites Eastern Ice Shelf (TEIS) buttresses the eastern grounded portion of Thwaites Glacier through contact with a pinning point at its seaward limit. Loss of this ice shelf will promote further acceleration of Thwaites Glacier. Understanding the dynamic controls and structural integrity of the TEIS is therefore important to estimating Thwaites' future sea-level contribution. We present a ∼ 20-year record of change on the TEIS that reveals the dynamic controls governing the ice shelf's past behaviour and ongoing evolution. We derived ice velocities from MODIS and Sentinel-1 image data using feature tracking and speckle tracking, respectively, and we combined these records with ITS_LIVE and GOLIVE velocity products from Landsat-7 and Landsat-8. In addition, we estimated surface lowering and basal melt rates using the Reference Elevation Model of Antarctica (REMA) DEM in comparison to ICESat and ICESat-2 altimetry. Early in the record, TEIS flow dynamics were strongly controlled by the neighbouring Thwaites Western Ice Tongue (TWIT). Flow patterns on the TEIS changed following the disintegration of the TWIT around 2008, with a new divergence in ice flow developing around the pinning point at its seaward limit. Simultaneously, the TEIS developed new rifting that extends from the shear zone upstream of the ice rise and increased strain concentration within this shear zone. As these horizontal changes occurred, sustained thinning driven by basal melt reduced ice thickness, particularly near the grounding line and in the shear zone area upstream of the pinning point. This evidence of weakening at a rapid pace suggests that the TEIS is likely to fully destabilize in the next few decades, leading to further acceleration of Thwaites Glacier.

The triggers of the disaggregation of Voyeykov Ice Shelf (2007), Wilkes Land, East Antarctica, and its subsequent evolution

AbstractThe weakening and/or removal of floating ice shelves in Antarctica can induce inland ice flow acceleration. Numerical modelling suggests these processes will play an important role in Antarctica's future sea-level contribution, but our understanding of the mechanisms that lead to ice tongue/shelf collapse is incomplete and largely based on observations from the Antarctic Peninsula and West Antarctica. Here, we use remote sensing of structural glaciology and ice velocity from 2001 to 2020 and analyse potential ocean-climate forcings to identify mechanisms that triggered the rapid disintegration of ~2445 km2of ice mélange and part of the Voyeykov Ice Shelf in Wilkes Land, East Antarctica between 27 March and 28 May 2007. Results show disaggregation was pre-conditioned by weakening of the ice tongue's structural integrity and was triggered by mélange removal driven by a regional atmospheric circulation anomaly and a less extensive latent-heat polynya. Disaggregation did not induce inland ice flow acceleration, but our observations highlight an important mechanism through which floating termini can be removed, whereby the break-out of mélange and multiyear landfast sea ice triggers disaggregation of a structurally-weak ice shelf. These observations highlight the need for numerical ice-sheet models to account for interactions between sea-ice, mélange and ice shelves.

Future sea level contribution from Antarctica inferred from CMIP5 model forcing and its dependence on precipitation ansatz

Abstract. Various observational estimates indicate growing mass loss at Antarctica's margins as well as heavier precipitation across the continent. Simulated future projections reveal that heavier precipitation, falling on Antarctica, may counteract amplified iceberg discharge and increased basal melting of floating ice shelves driven by a warming ocean. Here, we test how the ansatz (implementation in a mathematical framework) of the precipitation boundary condition shapes Antarctica's sea level contribution in an ensemble of ice sheet simulations. We test two precipitation conditions: we either apply the precipitation anomalies from CMIP5 models directly or scale the precipitation by the air temperature anomalies from the CMIP5 models. In the scaling approach, it is common to use a relative precipitation increment per degree warming as an invariant scaling constant. We use future climate projections from nine CMIP5 models, ranging from strong mitigation efforts to business-as-usual scenarios, to perform simulations from 1850 to 5000. We take advantage of individual climate projections by exploiting their full temporal and spatial structure. The CMIP5 projections beyond 2100 are prolonged with reiterated forcing that includes decadal variability; hence, our study may underestimate ice loss after 2100. In contrast to various former studies that apply an evolving temporal forcing that is spatially averaged across the entire Antarctic Ice Sheet, our simulations consider the spatial structure in the forcing stemming from various climate patterns. This fundamental difference reproduces regions of decreasing precipitation despite general warming. Regardless of the boundary and forcing conditions applied, our ensemble study suggests that some areas, such as the glaciers from the West Antarctic Ice Sheet draining into the Amundsen Sea, will lose ice in the future. In general, the simulated ice sheet thickness grows along the coast, where incoming storms deliver topographically controlled precipitation. In this region, the ice thickness differences are largest between the applied precipitation methods. On average, Antarctica shrinks for all future scenarios if the air temperature anomalies scale the precipitation. In contrast, Antarctica gains mass in our simulations if we apply the simulated precipitation anomalies directly. The analysis reveals that the mean scaling inferred from climate models is larger than the commonly used values deduced from ice cores; moreover, it varies spatially: the highest scaling is across the East Antarctic Ice Sheet, and the lowest scaling is around the Siple Coast, east of the Ross Ice Shelf. The discrepancies in response to both precipitation ansatzes illustrate the principal uncertainty in projections of Antarctica's sea level contribution.

The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6

Abstract. The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean.

Damage accelerates ice shelf instability and mass loss in Amundsen Sea Embayment

Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are among the fastest changing outlet glaciers in West Antarctica with large consequences for global sea level. Yet, assessing how much and how fast both glaciers will weaken if these changes continue remains a major uncertainty as many of the processes that control their ice shelf weakening and grounding line retreat are not well understood. Here, we combine multisource satellite imagery with modeling to uncover the rapid development of damage areas in the shear zones of Pine Island and Thwaites ice shelves. These damage areas consist of highly crevassed areas and open fractures and are first signs that the shear zones of both ice shelves have structurally weakened over the past decade. Idealized model results reveal moreover that the damage initiates a feedback process where initial ice shelf weakening triggers the development of damage in their shear zones, which results in further speedup, shearing, and weakening, hence promoting additional damage development. This damage feedback potentially preconditions these ice shelves for disintegration and enhances grounding line retreat. The results of this study suggest that damage feedback processes are key to future ice shelf stability, grounding line retreat, and sea level contributions from Antarctica. Moreover, they underline the need for incorporating these feedback processes, which are currently not accounted for in most ice sheet models, to improve sea level rise projections.

Intermittent structural weakening and acceleration of the Thwaites Glacier Tongue between 2000 and 2018

AbstractEvolving conditions at the terminus of Thwaites Glacier will be important in determining the rate of its future sea-level contribution over the coming decades. Here, we use remote-sensing observations to investigate recent changes (2000–2018) in the structure and velocity of Thwaites Glacier and its floating tongue. We show that the main trunk of Thwaites Glacier has accelerated by 38% over this period, while its previously intact floating tongue has transitioned to a weaker mélange of fractured icebergs bounded by sea ice. However, the rate of structural weakening and acceleration was not uniform across the observational period and we identify two periods of rapid acceleration and structural weakening (2006–2012; 2016–2018), separated by a period of deceleration and re-advance of the structurally-intact shear margin boundary (2012–2015). The timing of these accelerations/decelerations strongly suggests a link to variable ocean forcing. The weakened tongue now has some dependency on landfast sea ice for structural integrity and is vulnerable to changes in landfast ice persistency. Future reductions in landfast sea ice could manifest from changes in climate and/or the imminent removal of the B-22A iceberg from the Thwaites embayment. Such changes could have important implications for the integrity of the ice tongue and future glacier discharge.

Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution

Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.

Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) – Part 2: Parameter ensemble analysis

Abstract. The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation–age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model–data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016). This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4±4.1 m (or 6.5±2.0×106km3), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.

Glacier algae accelerate melt rates on the south-western Greenland Ice Sheet

Abstract. Melting of the Greenland Ice Sheet (GrIS) is the largest single contributor to eustatic sea level and is amplified by the growth of pigmented algae on the ice surface, which increases solar radiation absorption. This biological albedo-reducing effect and its impact upon sea level rise has not previously been quantified. Here, we combine field spectroscopy with a radiative-transfer model, supervised classification of unmanned aerial vehicle (UAV) and satellite remote-sensing data, and runoff modelling to calculate biologically driven ice surface ablation. We demonstrate that algal growth led to an additional 4.4–6.0 Gt of runoff from bare ice in the south-western sector of the GrIS in summer 2017, representing 10 %–13 % of the total. In localized patches with high biomass accumulation, algae accelerated melting by up to 26.15±3.77 % (standard error, SE). The year 2017 was a high-albedo year, so we also extended our analysis to the particularly low-albedo 2016 melt season. The runoff from the south-western bare-ice zone attributed to algae was much higher in 2016 at 8.8–12.2 Gt, although the proportion of the total runoff contributed by algae was similar at 9 %–13 %. Across a 10 000 km2 area around our field site, algae covered similar proportions of the exposed bare ice zone in both years (57.99 % in 2016 and 58.89 % in 2017), but more of the algal ice was classed as “high biomass” in 2016 (8.35 %) than 2017 (2.54 %). This interannual comparison demonstrates a positive feedback where more widespread, higher-biomass algal blooms are expected to form in high-melt years where the winter snowpack retreats further and earlier, providing a larger area for bloom development and also enhancing the provision of nutrients and liquid water liberated from melting ice. Our analysis confirms the importance of this biological albedo feedback and that its omission from predictive models leads to the systematic underestimation of Greenland's future sea level contribution, especially because both the bare-ice zones available for algal colonization and the length of the biological growth season are set to expand in the future.

Cloud microphysics and circulation anomalies control differences in future Greenland melt

Recently, the Greenland Ice Sheet (GrIS) has become the main source of barystatic sea-level rise1,2. The increase in the GrIS melt is linked to anticyclonic circulation anomalies, a reduction in cloud cover and enhanced warm-air advection3–7. The Climate Model Intercomparison Project fifth phase (CMIP5) General Circulation Models (GCMs) do not capture recent circulation dynamics; therefore, regional climate models (RCMs) driven by GCMs still show significant uncertainties in future GrIS sea-level contribution, even within one emission scenario5,8–10. Here, we use the RCM Modele Atmospherique Regional to show that the modelled cloud water phase is the main source of disagreement among future GrIS melt projections. We show that, in the current climate, anticyclonic circulation results in more melting than under a neutral-circulation regime. However, we find that the GrIS longwave cloud radiative effect is extremely sensitive to the modelled cloud liquid-water path, which explains melt anomalies of +378 Gt yr–1 (+1.04 mm yr–1 global sea level equivalent) in a +2 °C-warmer climate with a neutral-circulation regime (equivalent to 21% more melt than under anticyclonic circulation). The discrepancies between modelled cloud properties within a high-emission scenario introduce larger uncertainties in projected melt volumes than the difference in melt between low- and high-emission scenarios11. Greenland Ice Sheet melt is contributing to sea-level rise; however, uncertainties exist about its future contributions. A regional climate model shows that clouds are the primary cause of this uncertainty, with melt varying significantly depending on the cloud water phase and atmospheric circulation.

Exploration of Antarctic Ice Sheet 100-year contribution to sea level rise and associated model uncertainties using the ISSM framework

Abstract. Estimating the future evolution of the Antarctic Ice Sheet (AIS) is critical for improving future sea level rise (SLR) projections. Numerical ice sheet models are invaluable tools for bounding Antarctic vulnerability; yet, few continental-scale projections of century-scale AIS SLR contribution exist, and those that do vary by up to an order of magnitude. This is partly because model projections of future sea level are inherently uncertain and depend largely on the model's boundary conditions and climate forcing, which themselves are unknown due to the uncertainty in the projections of future anthropogenic emissions and subsequent climate response. Here, we aim to improve the understanding of how uncertainties in model forcing and boundary conditions affect ice sheet model simulations. With use of sampling techniques embedded within the Ice Sheet System Model (ISSM) framework, we assess how uncertainties in snow accumulation, ocean-induced melting, ice viscosity, basal friction, bedrock elevation, and the presence of ice shelves impact continental-scale 100-year model simulations of AIS future sea level contribution. Overall, we find that AIS sea level contribution is strongly affected by grounding line retreat, which is driven by the magnitude of ice shelf basal melt rates and by variations in bedrock topography. In addition, we find that over 1.2 m of AIS global mean sea level contribution over the next century is achievable, but not likely, as it is tenable only in response to unrealistically large melt rates and continental ice shelf collapse. Regionally, we find that under our most extreme 100-year warming experiment generalized for the entire ice sheet, the Amundsen Sea sector is the most significant source of model uncertainty (1032 mm 6σ spread) and the region with the largest potential for future sea level contribution (297 mm). In contrast, under a more plausible forcing informed regionally by literature and model sensitivity studies, the Ronne basin has a greater potential for local increases in ice shelf basal melt rates. As a result, under this more likely realization, where warm waters reach the continental shelf under the Ronne ice shelf, it is the Ronne basin, particularly the Evans and Rutford ice streams, that are the greatest contributors to potential SLR (161 mm) and to simulation uncertainty (420 mm 6σ spread).

Simulation of the future sea level contribution of Greenland with a new glacial system model

Abstract. We introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. The aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961–1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961–1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation–surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm until the year 2300, including our simulated additional sea level rise due to elevation–surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.

Projections of Future Sea Level Contributions from the Greenland and Antarctic Ice Sheets: Challenges Beyond Dynamical Ice Sheet Modeling

As Earth's climate warms, rising sea levels are becoming a great concern. Providing substantiated and well-informed guidance on the amount and rate of future sea level rise is important, but remains challenging. Observations of the Greenland and Antarctic Ice Sheets-Earth's largest freshwater reservoirs-reveal their ongoing, rapid, and complex changes in response to an evolving climate. The improved understanding of ice sheet behavior resulting from observations is driving the development of ice sheet models, and is allowing better simulation of past, present, and future ice sheet evolution. However, insight into future changes requires a better understanding of how ice sheets interact with other components of the Earth system and associated feedbacks. As climate models start to include dynamical ice sheet components, our understanding of such interactions and the feedback mechanisms will advance. These new developments in ice sheet modeling and their implementation in climate models are timely, as observations indicate an accelerating contribution of ice sheets to sea level rise.