Search

Improving our confidence in future projections of sea-level rise requires models that can simulate ice-sheet evolution both in the future and in the geological past. A physically accurate treatment of large changes in ice-sheet geometry requires a proper treatment of processes near the margin, like grounding line dynamics, which in turn requires a high spatial resolution in that specific region, so that small-scale topographical features are resolved. This leads to a demand for computationally efficient models, where such a high resolution can be feasibly applied in simulations of 105-107 years in duration. Here, we present and evaluate a new ice-sheet model that solves the hybrid SIA-SSA approximation of the stress balance, including a heuristic rule for the groundingline flux. This is done on a dynamic adaptive mesh which is adapted to the modelled ice-sheet geometry during a simulation. Mesh resolution can be configured to be fine only at specified areas, such as the calving front or the grounding line, as well as specified point locations such as ice-core drill sites. This strongly reduces the number of grid points where the equations need to be solved, increasing the computational efficiency. A high resolution allows the model to resolve small geometrical features, such as outlet glaciers and sub-shelf pinning points, which can significantly affect large-scale ice-sheet dynamics. We show that the model reproduces the analytical solutions or model intercomparison benchmarks for a number of schematic ice-sheet configurations, indicating that the numerical approach is valid. Because of the unstructured triangular mesh, the number of vertices increases less rapidly with resolution than in a squaregrid model, greatly reducing the required computation time for high resolutions. A simulation of all four continental ice sheets during an entire 120 kyr glacial cycle, with a 4 km resolution near the grounding line, is expected to take 100- 200 wall clock hours on a 16-core system (1600-3200 core hours), implying that this model can be feasibly used for high-resolution palaeo-ice-sheet simulations., SCOPUS: ar.j, info:eu-repo/semantics/published

There is now much geological evidence that the Earth was fully glaciated during several periods in the geological past (about 700[Formula: see text]Myr ago) and attained a so-called Snowball Earth (SBE) state. Additional support for this idea has come from climate models of varying complexity that show transitions to SBE states and undergo hysteresis under changes in solar radiation. In this paper, we apply large-scale bifurcation analyses to a novel, fully-implicit Earth System Model of Intermediate Complexity (I-EMIC) to study SBE transitions. The I-EMIC contains a primitive equation ocean model, a model for atmospheric heat and moisture transport, a sea ice component and formulations for the adjustment of albedo over snow and ice. With the I-EMIC, high-dimensional branches of the SBE bifurcation diagram are obtained through parameter continuation. We are able to identify stable and unstable equilibria and uncover an intricate bifurcation structure associated with the ice-albedo feedback. Moreover, large-scale linear stability analyses are performed near major bifurcations, revealing the spatial nature of destabilizing perturbations.

The accurate quantification of uncertainties in regional sea-level projections is essential for guiding policy makers. As climate models do not currently simulate total sea level, these uncertainties must be quantified through summation of uncertainties in individual sea-level components. This summation depends on the correlation between the components, which has previously been prescribed or derived from each individual component's dependence on global mean surface temperature. In this study, we quantify, for the first time, regional correlations between sea-level components based on regional climate change projections. We compute regional sea-level projections consistent with climate projections from an ensemble of 14 Earth System Models. From the multi-model spread, we estimate the uncertainty in the regional climate's response to greenhouse forcing. To quantify the total uncertainty, we add the uncertainty in the response of sea-level components to this regional climate change. This approach reveals how regional climate processes impose correlations between sea-level components, affecting the total uncertainty. One example is an anti-correlation between North Atlantic sterodynamic change and Antarctic dynamic mass loss, suggesting a teleconnection established by the large-scale ocean circulation. We find that prescribed correlations, applied in the fifth assessment report of the Intergovernmental Panel on Climate Change, lead to a global overestimation in the uncertainty in regional sea-level projections on the order of 20%. Regionally, this overestimation exceeds 100%. We conclude that accurate uncertainty estimates of regional sea-level change must be based on projections of regional climate change and cannot be derived from global indicators such as global mean surface temperature.

Observations and models agree that the Greenland Ice Sheet (GrIS) surface mass balance (SMB) has decreased since the end of the 1990s due to an increase in meltwater runoff and that this trend will accelerate in the future. However, large uncertainties remain, partly due to different approaches for modelling the GrIS SMB, which have to weigh physical complexity or low computing time, different spatial and temporal resolutions, different forcing fields, and different ice sheet topographies and extents, which collectively make an inter-comparison difficult. Our GrIS SMB model intercomparison project (GrSMBMIP) aims to refine these uncertainties by intercomparing 13 models of four types which were forced with the same ERA-Interim reanalysis forcing fields, except for two global models. We interpolate all modelled SMB fields onto a common ice sheet mask at 1 km horizontal resolution for the period 1980-2012 and score the outputs against (1) SMB estimates from a combination of gravimetric remote sensing data from GRACE and measured ice discharge; (2) ice cores, snow pits and in situ SMB observations; and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting model deficiencies in an accurate representation of the GrIS ablation zone extent and processes related to surface melt and runoff. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of the same order as RCMs compared with observations and therefore remain useful tools for long-term simulations or coupling with ice sheet models. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present-day SMB relative to observations, suggesting that biases are not systematic among models and that this ensemble estimate can be used as a reference for current climate when carrying out future model developments. However, a higher density of in situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 m w.e. yr-1 due to large discrepancies in modelled snowfall accumulation.

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and inform on the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimated the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes and the forcings employed. This study presents results from 18 simulations from 15 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100, forced with different scenarios from the Coupled Model Intercomparison Project Phase 5 (CMIP5) representative of the spread in climate model results. The contribution of the Antarctic ice sheet in response to increased warming during this period varies between −7.8 and 30.0 cm of Sea Level Equivalent (SLE). The evolution of the West Antarctic Ice Sheet varies widely among models, with an overall mass loss up to 21.0 cm SLE in response to changes in oceanic conditions. East Antarctica mass change varies between −6.5 and 16.5 cm SLE, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional mass loss of 8 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 AOGCMs show an overall mass loss of 10 mm SLE compared to simulations done under present-day conditions, with limited mass gain in East Antarctica.

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015-2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between 7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between 6 :1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

Future sea-level change projections with process-based stand-alone ice sheet models are typically driven with surface mass balance (SMB) forcing derived from climate models. In this work we address the problems arising from a mismatch of the modelled ice sheet geometry with the geometry used by the climate model. We present a method for applying SMB forcing from climate models to a wide range of Greenland ice sheet models with varying and temporally evolving geometries. In order to achieve that, we translate a given SMB anomaly field as a function of absolute location to a function of surface elevation for 25 regional drainage basins, which can then be applied to different modelled ice sheet geometries. The key feature of the approach is the non-locality of this remapping process. The method reproduces the original forcing data closely when remapped to the original geometry. When remapped to different modelled geometries it produces a physically meaningful forcing with smooth and continuous SMB anomalies across basin divides. The method considerably reduces non-physical biases that would arise by applying the SMB anomaly derived for the climate model geometry directly to a large range of modelled ice sheet model geometries., SCOPUS: ar.j, info:eu-repo/semantics/published

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., info:eu-repo/semantics/published