Contribution To Sea Level Rise Research Articles

Coupled ice–ocean interactions during future retreat of West Antarctic ice streams in the Amundsen Sea sector

Abstract. The Amundsen Sea sector has some of the fastest-thinning ice shelves in Antarctica, caused by high, ocean-driven basal melt rates, which can lead to increased ice streamflow, causing increased sea level rise (SLR) contributions. In this study, we present the results of a new synchronously coupled ice-sheet–ocean model of the Amundsen Sea sector. We use the Wavelet-based, Adaptive-grid, Vertically Integrated ice sheet model (WAVI) to solve for ice velocities and the Massachusetts Institute of Technology general circulation model (MITgcm) to solve for ice thickness and three-dimensional ocean properties, allowing for full mass conservation in the coupled ice–ocean system. The coupled model is initialised in the present day and run forward under idealised warm and cold ocean conditions with a fixed ice front. We find that Thwaites Glacier dominates the future SLR from the Amundsen Sea sector, with a SLR that evolves approximately quadratically over time. The future evolution of Thwaites Glacier depends on the lifespan of small pinning points that form during the retreat. The rate of melting around these pinning points provides the link between future ocean conditions and the SLR from this sector and will be difficult to capture without a coupled ice–ocean model. Grounding-line retreat leads to a progressively larger Thwaites Ice Shelf cavity, leading to a positive trend in total melting, resulting from the increased ice basal surface area. Despite these important sensitivities, Thwaites Glacier retreats even in a scenario with zero ocean-driven melting. This demonstrates that a tipping point may have been passed in these simulations and some SLR from this sector is now committed.

Responses of the Pine Island and Thwaites glaciers to melt and sliding parameterizations

Abstract. The Pine Island and Thwaites glaciers are the two largest contributors to sea level rise from Antarctica. Here we examine the influence of basal friction and ice shelf basal melt in determining projected losses. We examine both Weertman and Coulomb friction laws with explicit weakening as the ice thins to flotation, which many friction laws include implicitly via the effective pressure. We find relatively small differences with the choice of friction law (Weertman or Coulomb) but find losses to be highly sensitive to the rate at which the basal traction is reduced as the area upstream of the grounding line thins. Consistent with earlier work on Pine Island Glacier, we find sea level contributions from both glaciers to vary linearly with the melt volume averaged over time and space, with little influence from the spatial or temporal distribution of melt. Based on recent estimates of melt from other studies, our simulations suggest that the combined melt-driven and sea level rise contribution from both glaciers may not exceed 10 cm by 2200, although the uncertainty in model parameters allows for larger increases. We do not include other factors, such as ice shelf breakup, that might increase loss, or factors such as increased accumulation and isostatic uplift that may mitigate loss.

Coupling MAR (Modèle Atmosphérique Régional) with PISM (Parallel Ice Sheet Model) mitigates the positive melt–elevation feedback

Abstract. The Greenland Ice Sheet is a key contributor to sea level rise. By melting, the ice sheet thins, inducing higher surface melt due to lower surface elevations, accelerating the melt coming from global warming. This process is called the melt–elevation feedback and can be considered by using two types of models: either (1) atmospheric models, which can represent the surface mass balance (SMB), or SMB estimates resulting from simpler models such as positive degree day models or (2) ice sheet models representing the surface elevation evolution. The latter ones do not represent the surface mass balance explicitly as well as polar-oriented climate models. A new coupling between the MAR (Modèle Atmosphérique Régional) regional climate model and the PISM (Parallel Ice Sheet Model) ice sheet model is presented here following the CESM2 (Community Earth System Model; SSP5-8.5, Shared Socioeconomic Pathway) scenario until 2100 at the MAR lateral boundaries. The coupling is extended to 2200 with a stabilised climate (+7 ∘C compared to 1961–1990) by randomly sampling the last 10 years of CESM2 to force MAR and reaches a sea level rise contribution of 64 cm. The fully coupled simulation is compared to a one-way experiment where surface topography remains fixed in MAR. However, the surface mass balance is corrected for the melt–elevation feedback when interpolated on the PISM grid by using surface mass balance vertical gradients as a function of local elevation variations (offline correction). This method is often used to represent the melt–elevation feedback and prevents a coupling which is too expensive in computation time. In the fully coupled MAR simulation, the ice sheet morphology evolution (changing slope and reducing the orographic barrier) induces changes in local atmospheric patterns. More specifically, wind regimes are modified, as well as temperature lapse rates, influencing the melt rate through modification of sensible heat fluxes at the ice sheet margins. We highlight mitigation of the melt lapse rate on the margins by modifying the surface morphology. The lapse rates considered by the offline correction are no longer valid at the ice sheet margins. If used (one-way simulation), this correction implies an overestimation of the sea level rise contribution of 2.5 %. The mitigation of the melt lapse rate on the margins can only be corrected by using a full coupling between an ice sheet model and an atmospheric model.

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.

Rapid and synchronous response of outlet glaciers to ocean warming on the Barents Sea coast, Novaya Zemlya

Abstract The Arctic is a hotspot for climate warming, making it crucial to quantify the sea level rise contribution from its ice masses. Novaya Zemlya's ice caps are the largest glacier complex in Europe and are a major contributor to contemporary sea level rise. Here we show that Novaya Zemlya outlet glaciers on the Barents Sea coast respond rapidly and consistently to oceanic forcing at annual timescales, likely due to their exposure to Atlantic Water variability. Glaciers on the Kara Sea show more variable response, likely reflecting their reduced exposure to Atlantic Water. Data demonstrate that the pause in glacier retreat previously observed on Novaya Zemlya between 2013 and 2015 has not persisted and that these changes correspond to ocean temperature variability on the Barents Sea coast. We document a marked shift to warmer air and ocean temperatures, and reduced sea ice concentrations from 2005 onwards. Although we identify ocean warming as the primary trigger for glacier retreat, we suggest that multi-year thinning, driven by the shift towards warmer air temperatures since 2005, pre-conditioned Novaya Zemlya's glaciers to retreat. Despite commonality in the timing of outlet glacier retreat, the magnitude is highly variable during rapid retreat phases, which we attribute to glacier-specific factors.

A framework for time-dependent ice sheet uncertainty quantification, applied to three West Antarctic ice streams

Abstract. Ice sheet models are the main tool to generate forecasts of ice sheet mass loss, a significant contributor to sea level rise; thus, knowing the likelihood of such projections is of critical societal importance. However, to capture the complete range of possible projections of mass loss, ice sheet models need efficient methods to quantify the forecast uncertainty. Uncertainties originate from the model structure, from the climate and ocean forcing used to run the model, and from model calibration. Here we quantify the latter, applying an error propagation framework to a realistic setting in West Antarctica. As in many other ice sheet modelling studies we use a control method to calibrate grid-scale flow parameters (parameters describing the basal drag and ice stiffness) with remotely sensed observations. Yet our framework augments the control method with a Hessian-based Bayesian approach that estimates the posterior covariance of the inverted parameters. This enables us to quantify the impact of the calibration uncertainty on forecasts of sea level rise contribution or volume above flotation (VAF) due to the choice of different regularization strengths (prior strengths), sliding laws, and velocity inputs. We find that by choosing different satellite ice velocity products our model leads to different estimates of VAF after 40 years. We use this difference in model output to quantify the variance that projections of VAF are expected to have after 40 years and identify prior strengths that can reproduce that variability. We demonstrate that if we use prior strengths suggested by L-curve analysis, as is typically done in ice sheet calibration studies, our uncertainty quantification is not able to reproduce that same variability. The regularization suggested by the L curves is too strong, and thus propagating the observational error through to VAF uncertainties under this choice of prior leads to errors that are smaller than those suggested by our two-member “sample” of observed velocity fields.

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.

The acceleration of sea-level rise along the coast of the Netherlands started in the 1960s

Abstract. The global acceleration of sea-level rise (SLR) during the 20th century is now established. On the local scale, this is harder to establish as several drivers of SLR play a role, which can mask the acceleration. Here, we study the rate of SLR along the coast of the Netherlands from the average of six tide gauge records covering the period 1890–2021. To isolate the effects of the wind field variations and the nodal tide from the local sea-level trend, we use four generalised additive models (GAMs) which include different predictive variables. From the sea-level trend estimates, we obtain the continuous evolution of the rate of SLR and its uncertainty over the observational period. The standard error in the estimation of the rate of SLR is reduced when we account for nodal-tide effects and is reduced further when we also account for the wind effects, meaning these provide better estimates of the rate of SLR. A part of the long-term SLR is due to wind forcing related to a strengthening and northward shift of the jet stream, but this SLR contribution decelerated over the observational period. Additionally, we detect wind-forced sea-level variability on multidecadal timescales with an amplitude of around 1 cm. Using a coherence analysis, we identify both the North Atlantic Oscillation and the Atlantic Multidecadal Variability as its drivers. Crucially, accounting for the nodal-tide and wind effects changes the estimated rate of SLR, unmasking an SLR acceleration that started in the 1960s. Our best-fitting GAM, which accounts for nodal and wind effects, yields a rate of SLR of about 1.72.21.3 mm yr−1 in 1900–1919 and 1.51.91.2 mm yr−1 in 1940–1959 compared to 2.93.52.4 mm yr−1 in 2000–2019 (where the lower and upper bounds denote the 5th and 95th percentiles). If we discount the nodal tide, wind and fluctuation effects and assume a constant rate of SLR, then the probability (p value) of finding a rate difference between 1940–1959 and 2000–2019 of at least our estimate is smaller than 1 %. Consistent with global observations and the expectations based on the physics of global warming, our results show unequivocally that SLR along the Dutch coast has accelerated since the 1960s.

Seasonal variability in Antarctic ice shelf velocities forced by sea surface height variations

Abstract. Antarctica's ice shelves resist the flow of grounded ice towards the ocean through “buttressing” arising from their contact with ice rises, rumples, and lateral margins. Ice shelf thinning and retreat reduce buttressing, leading to increased delivery of mass to the ocean that adds to global sea level. Ice shelf response to large annual cycles in atmospheric and oceanic processes provides opportunities to study the dynamics of both ice shelves and the buttressed grounded ice. Here, we explore whether seasonal variability of sea surface height (SSH) can explain observed seasonal variability of ice velocity. We investigate this hypothesis using several time series of ice velocity from the Ross Ice Shelf (RIS), satellite-based estimates of SSH seaward of the RIS front, ocean models of SSH under and near RIS, and a viscous ice sheet model. The observed annual changes in RIS velocity are of the order of 1–10 m a−1 (roughly 1 % of mean flow). The ice sheet model, forced by the observed and modelled range of SSH of about 10 cm, reproduces the observed velocity changes when sufficiently large basal drag changes near the grounding line are parameterised. The model response is dominated by grounding line migration but with a significant contribution from SSH-induced tilt of the ice shelf. We expect that climate-driven changes in the seasonal cycles of winds and upper-ocean summer warming will modify the seasonal response of ice shelves to SSH and that nonlinear responses of the ice sheet will affect the longer trend in ice sheet response and its potential sea-level rise contribution.

Future sea-level projections with a coupled atmosphere-ocean-ice-sheet model

Climate-forced, offline ice-sheet model simulations have been used extensively in assessing how much ice-sheets can contribute to future global sea-level rise. Typically, these model projections do not account for the two-way interactions between ice-sheets and climate. To quantify the impact of ice-ocean-atmosphere feedbacks, here we conduct greenhouse warming simulations with a coupled global climate-ice-sheet model of intermediate complexity. Following the Shared Socioeconomic Pathway (SSP) 1-1.9, 2-4.5, 5-8.5 emission scenarios, the model simulations ice-sheet contributions to global sea-level rise by 2150 of 0.2 ± 0.01, 0.5 ± 0.01 and 1.4 ± 0.1 m, respectively. Antarctic ocean-ice-sheet-ice-shelf interactions enhance future subsurface basal melting, while freshwater-induced atmospheric cooling reduces surface melting and iceberg calving. The combined effect is likely to decelerate global sea-level rise contributions from Antarctica relative to the uncoupled climate-forced ice-sheet model configuration. Our results demonstrate that estimates of future sea-level rise fundamentally depend on the complex interactions between ice-sheets, icebergs, ocean and the atmosphere.

The collapse of the Cordilleran–Laurentide ice saddle and early opening of the Mackenzie Valley, Northwest Territories, Canada, constrained by 10Be exposure dating

Abstract. Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ∼15.8 ka (17.1–14.6 ka), ∼1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ∼14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ∼11.2 m between 16 and 13 ka.

Surface melt on the Shackleton Ice Shelf, East Antarctica (2003–2021)

Abstract. Melt on the surface of Antarctic ice shelves can potentially lead to their disintegration, accelerating the flow of grounded ice to the ocean and raising global sea levels. However, the current understanding of the processes driving surface melt is incomplete, increasing uncertainty in predictions of ice shelf stability and thus of Antarctica's contribution to sea-level rise. Previous studies of surface melt in Antarctica have usually focused on either a process-level understanding of melt through energy-balance investigations or used metrics such as the annual number of melt days to quantify spatiotemporal variability in satellite observations of surface melt. Here, we help bridge the gap between work at these two scales. Using daily passive microwave observations from the AMSR-E and AMSR-2 sensors and the machine learning approach of a self-organising map, we identify nine representative spatial distributions (“patterns”) of surface melt on the Shackleton Ice Shelf in East Antarctica from 2002/03–2020/21. Combined with output from the RACMO2.3p3 regional climate model and surface topography from the REMA digital elevation model, our results point to a significant role for surface air temperatures in controlling the interannual variability in summer melt and also reveal the influence of localised controls on melt. In particular, prolonged melt along the grounding line shows the importance of katabatic winds and surface albedo. Our approach highlights the necessity of understanding both local and large-scale controls on surface melt and demonstrates that self-organising maps can be used to investigate the variability in surface melt on Antarctic ice shelves.

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.

Relative sea level rise contributions to flood construction levels in British Columbia

This paper proposes an updated approach to incorporating relative sea level rise in the estimation of flood construction levels in British Columbia, which are required in coastal planning and infrastructure development. Initially, current guidelines for British Columbia are outlined, and recent information on projected relative sea level rise is summarized. Then, recognizing the uncertainty in sea level rise projections, a probabilistic model of estimating the contribution of relative sea level rise to the flood construction level is described, relying on an available theoretical formulation that is corroborated by the application of a Monte Carlo simulation method. An updated approach to incorporating relative sea level rise in the estimation of flood construction levels is thereby recommended, relying on available data portals for relative sea level rise. Finally, the estimation of the probability of overtopping is described, and the implications of this for making informed decisions on flood construction levels with respect to risk assessments are indicated.

Quantifying the Impact of Bedrock Topography Uncertainty in Pine Island Glacier Projections for This Century

AbstractThe predicted Antarctic contribution to global‐mean sea‐level rise is one of the most uncertain among all major sources. Partly this is because of instability mechanisms of the ice flow over deep basins. Errors in bedrock topography can substantially impact the projected resilience of glaciers against such instabilities. Here we analyze the Pine Island Glacier topography to derive a statistical model representation. Our model allows for inhomogeneous and spatially dependent uncertainties and avoids unnecessary smoothing from spatial averaging or interpolation. A set of topography realizations is generated representing our best estimate of the topographic uncertainty in ice sheet model simulations. The bedrock uncertainty alone creates a 5%–25% uncertainty in the predicted sea level rise contribution at year 2100, depending on friction law and climate forcing. Pine Island Glacier simulations on this new set are consistent with simulations on the BedMachine reference topography but diverge from Bedmap2 simulations.

Brief communication: Increased glacier mass loss in the Russian High Arctic (2010–2017)

Abstract. Glaciers in the Russian High Arctic have been subject to extensive atmospheric warming due to global climate change, yet their contribution to sea level rise has been relatively small over the past decades. Here we show surface elevation change measurements and geodetic mass balances of 93 % of all glacierized areas of Novaya Zemlya, Severnaya Zemlya, and Franz Josef Land using interferometric synthetic aperture radar measurements taken between 2010 and 2017. We calculate an overall mass loss rate of -22±6 Gt a−1, corresponding to a sea level rise contribution of 0.06±0.02 mm a−1. Compared to measurements prior to 2010, mass loss of glaciers on the Russian archipelagos has doubled in recent years.

Projected increases in surface melt and ice loss for the Northern and Southern Patagonian Icefields

The Northern Patagonian Icefield (NPI) and the Southern Patagonian Icefield (SPI) have increased their ice mass loss in recent decades. In view of the impacts of glacier shrinkage in Patagonia, an assessment of the potential future surface mass balance (SMB) of the icefields is critical. We seek to provide this assessment by modelling the SMB between 1976 and 2050 for both icefields, using regional climate model data (RegCM4.6) and a range of emission scenarios. For the NPI, reductions between 1.5 m w.e. (RCP2.6) and 1.9 m w.e. (RCP8.5) were estimated in the mean SMB during the period 2005–2050 compared to the historical period (1976–2005). For the SPI, the estimated reductions were between 1.1 m w.e. (RCP2.6) and 1.5 m w.e. (RCP8.5). Recently frontal ablation estimates suggest that mean SMB in the SPI is positively biased by 1.5 m w.e., probably due to accumulation overestimation. If it is assumed that frontal ablation rates of the recent past will continue, ice loss and sea-level rise contribution will increase. The trend towards lower SMB is mostly explained by an increase in surface melt. Positive ice loss feedbacks linked to increasing in meltwater availability are expected for calving glaciers.

The Role of Beach Morphology and Mid-Century Climate Change Effects on Wave Runup and Storm Impact on the Northern Yucatan Coast

Wave runup is a relevant parameter to determine the storm impact on barrier islands. Here, the role of the beach morphology on wave runup and storm impact was investigated at four coastal communities located on the northern Yucatan coast. Current wave conditions based on regional wind simulations, topo-bathymetric transects measured at each location, and a nonlinear wave transformation model were employed to reconstruct multi-year runup time series. Dune morphology features and extreme water levels (excluding storm surge contributions) were further employed to determine the storm impact at each site for different return periods. Despite the similar offshore conditions along the coast, extreme water levels (i.e., runup and setup) showed intersite differences that were mainly ascribed to subaerial and submerged morphological features. Numerical results showed that the average surf zone beach slope, sandbars, berm, and dune elevation played an important role in controlling extreme water levels and storm impact at the study sites under the present climate. Moreover, in order to assess the potential effect of climate change on coastal flooding, we analyzed wave runup and storm impact in the best-preserved site by considering wave conditions and sea level rise (SLR) projections under the RCP 8.5 scenario. Modelling results suggest no significant increase in the storm impact regime between the present and future conditions in the study area unless SLR is considered. It was found that to accurately estimate SLR contribution, it should be incorporated into mean sea level prior to performing numerical wave runup simulations, rather than simply adding it to the resulting wave-induced water levels.

Greater Greenland Ice Sheet contribution to global sea level rise in CMIP6

Future climate projections show a marked increase in Greenland Ice Sheet (GrIS) runoff during the 21st century, a direct consequence of the Polar Amplification signal. Regional climate models (RCMs) are a widely used tool to downscale ensembles of projections from global climate models (GCMs) to assess the impact of global warming on GrIS melt and sea level rise contribution. Initial results of the CMIP6 GCM model intercomparison project have revealed a greater 21st century temperature rise than in CMIP5 models. However, so far very little is known about the subsequent impacts on the future GrIS surface melt and therefore sea level rise contribution. Here, we show that the total GrIS sea level rise contribution from surface mass loss in our high-resolution (15 km) regional climate projections is 17.8 ± 7.8 cm in SSP585, 7.9 cm more than in our RCP8.5 simulations using CMIP5 input. We identify a +1.3 °C greater Arctic Amplification and associated cloud and sea ice feedbacks in the CMIP6 SSP585 scenario as the main drivers. Additionally, an assessment of the GrIS sea level contribution across all emission scenarios highlights, that the GrIS mass loss in CMIP6 is equivalent to a CMIP5 scenario with twice the global radiative forcing.

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.