Search Results
Description and validation of the ice-sheet model Yelmo (version 1.0)
2020, Robinson, Alexander, Alvarez-Solas, Jorge, Montoya, Marisa, Goelzer, Heiko, Greve, Ralf, Ritz, Catherine
We describe the physics and features of the ice-sheet model Yelmo, an open-source project intended for collaborative development. Yelmo is a thermomechanical model, solving for the coupled velocity and temperature solutions of an ice sheet simultaneously. The ice dynamics are currently treated via a “hybrid” approach combining the shallow-ice and shallow-shelf/shelfy-stream approximations, which makes Yelmo an apt choice for studying a wide variety of problems. Yelmo's main innovations lie in its flexible and user-friendly infrastructure, which promotes portability and facilitates long-term development. In particular, all physics subroutines have been designed to be self-contained, so that they can be easily ported from Yelmo to other models, or easily replaced by improved or alternate methods in the future. Furthermore, hard-coded model choices are eschewed, replaced instead with convenient parameter options that allow the model to be adapted easily to different contexts. We show results for different ice-sheet benchmark tests, and we illustrate Yelmo's performance for the Antarctic ice sheet.
The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model
2020, Reese, Ronja, Levermann, Anders, Albrecht, Torsten, Seroussi, Hélène, Winkelmann, Ricarda
Mass loss from the Antarctic Ice Sheet constitutes the largest uncertainty in projections of future sea level rise. Ocean-driven melting underneath the floating ice shelves and subsequent acceleration of the inland ice streams are the major reasons for currently observed mass loss from Antarctica and are expected to become more important in the future. Here we show that for projections of future mass loss from the Antarctic Ice Sheet, it is essential (1) to better constrain the sensitivity of sub-shelf melt rates to ocean warming and (2) to include the historic trajectory of the ice sheet. In particular, we find that while the ice sheet response in simulations using the Parallel Ice Sheet Model is comparable to the median response of models in three Antarctic Ice Sheet Intercomparison projects – initMIP, LARMIP-2 and ISMIP6 – conducted with a range of ice sheet models, the projected 21st century sea level contribution differs significantly depending on these two factors. For the highest emission scenario RCP8.5, this leads to projected ice loss ranging from 1.4 to 4.0 cm of sea level equivalent in simulations in which ISMIP6 ocean forcing drives the PICO ocean box model where parameter tuning leads to a comparably low sub-shelf melt sensitivity and in which no surface forcing is applied. This is opposed to a likely range of 9.1 to 35.8 cm using the exact same initial setup, but emulated from the LARMIP-2 experiments with a higher melt sensitivity, even though both projects use forcing from climate models and melt rates are calibrated with previous oceanographic studies. Furthermore, using two initial states, one with a previous historic simulation from 1850 to 2014 and one starting from a steady state, we show that while differences between the ice sheet configurations in 2015 seem marginal at first sight, the historic simulation increases the susceptibility of the ice sheet to ocean warming, thereby increasing mass loss from 2015 to 2100 by 5 % to 50 %. Hindcasting past ice sheet changes with numerical models would thus provide valuable tools to better constrain projections. Our results emphasize that the uncertainty that arises from the forcing is of the same order of magnitude as the ice dynamic response for future sea level projections.
Effects of extreme melt events on ice flow and sea level rise of the Greenland Ice Sheet
2023, Beckmann, Johanna, Winkelmann, Ricarda
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.
ISMIP6 Antarctica: A multi-model ensemble of the Antarctic ice sheet evolution over the 21st century
2020, Seroussi, Hélène, Nowicki, Sophie, Payne, Antony J., Goelzer, Heiko, Lipscomb, William H., Abe-Ouchi, Ayako, Agosta, Cécile, Albrecht, Torsten, Asay-Davis, Xylar, Barthel, Alice, Calov, Reinhard, Cullather, Richard, Dumas, Christophe, Galton-Fenzi, Benjamin K., Gladstone, Rupert, Golledge, Nicholas R., Gregory, Jonathan M., Greve, Ralf, Hattermann, Tore, Hoffman, Matthew J., Humbert, Angelika, Huybrechts, Philippe, Jourdain, Nicolas C., Kleiner, Thomas, Larour, Eric, Leguy, Gunter R., Lowry, Daniel P., Little, Chistopher M., Morlighem, Mathieu, Pattyn, Frank, Pelle, Tyler, Price, Stephen F., Quiquet, Aurélien, Reese, Ronja, Schlegel, Nicole-Jeanne, Shepherd, Andrew, Simon, Erika, Smith, Robin S., Straneo, Fiammetta, Sun, Sainan, Trusel, Luke D., Van Breedam, Jonas, van de Wal, Roderik S. W., Winkelmann, Ricarda, Zhao, Chen, Zhang, Tong, Zwinger, Thomas
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 presentday 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. © Author(s) 2020.
Global warming due to loss of large ice masses and Arctic summer sea ice
2020, Wunderling, Nico, Willeit, Matteo, Donges, Jonathan F., Winkelmann, Ricarda
Several large-scale cryosphere elements such as the Arctic summer sea ice, the mountain glaciers, the Greenland and West Antarctic Ice Sheet have changed substantially during the last century due to anthropogenic global warming. However, the impacts of their possible future disintegration on global mean temperature (GMT) and climate feedbacks have not yet been comprehensively evaluated. Here, we quantify this response using an Earth system model of intermediate complexity. Overall, we find a median additional global warming of 0.43 °C (interquartile range: 0.39−0.46 °C) at a CO2 concentration of 400 ppm. Most of this response (55%) is caused by albedo changes, but lapse rate together with water vapour (30%) and cloud feedbacks (15%) also contribute significantly. While a decay of the ice sheets would occur on centennial to millennial time scales, the Arctic might become ice-free during summer within the 21st century. Our findings imply an additional increase of the GMT on intermediate to long time scales.
Stabilizing the West Antarctic Ice Sheet by surface mass deposition
2019, Feldmann, Johannes, Levermann, Anders, Mengel, Matthias
There is evidence that a self-sustaining ice discharge from the West Antarctic Ice Sheet (WAIS) has started, potentially leading to its disintegration. The associated sea level rise of more than 3m would pose a serious challenge to highly populated areas including metropolises such as Calcutta, Shanghai, New York City, and Tokyo. Here, we show that the WAIS may be stabilized through mass deposition in coastal regions around Pine Island and Thwaites glaciers. In our numerical simulations, a minimum of 7400 Gt of additional snowfall stabilizes the flow if applied over a short period of 10 years onto the region (−2 mm year−1 sea level equivalent). Mass deposition at a lower rate increases the intervention time and the required total amount of snow. We find that the precise conditions of such an operation are crucial, and potential benefits need to be weighed against environmental hazards, future risks, and enormous technical challenges. Copyright © 2019 The Authors,
The tipping points and early warning indicators for Pine Island Glacier, West Antarctica
2021-3-25, Rosier, Sebastian H. R., Reese, Ronja, Donges, Jonathan F., De Rydt, Jan, Gudmundsson, G. Hilmar, Winkelmann, Ricarda
Mass loss from the Antarctic Ice Sheet is the main source of uncertainty in projections of future sea-level rise, with important implications for coastal regions worldwide. Central to ongoing and future changes is the marine ice sheet instability: once a critical threshold, or tipping point, is crossed, ice internal dynamics can drive a self-sustaining retreat committing a glacier to irreversible, rapid and substantial ice loss. This process might have already been triggered in the Amundsen Sea region, where Pine Island and Thwaites glaciers dominate the current mass loss from Antarctica, but modelling and observational techniques have not been able to establish this rigorously, leading to divergent views on the future mass loss of the West Antarctic Ice Sheet. Here, we aim at closing this knowledge gap by conducting a systematic investigation of the stability regime of Pine Island Glacier. To this end we show that early warning indicators in model simulations robustly detect the onset of the marine ice sheet instability. We are thereby able to identify three distinct tipping points in response to increases in ocean-induced melt. The third and final event, triggered by an ocean warming of approximately 1.2 ∘C from the steady-state model configuration, leads to a retreat of the entire glacier that could initiate a collapse of the West Antarctic Ice Sheet.
Future Sea Level Change Under Coupled Model Intercomparison Project Phase 5 and Phase 6 Scenarios From the Greenland and Antarctic Ice Sheets
2021, Payne, Antony J., Nowicki, Sophie, Abe‐Ouchi, Ayako, Agosta, Cécile, Alexander, Patrick, Albrecht, Torsten, Asay‐Davis, Xylar, Aschwanden, Andy, Barthel, Alice, Bracegirdle, Thomas J., Calov, Reinhard, Chambers, Christopher, Choi, Youngmin, Cullather, Richard, Cuzzone, Joshua, Dumas, Christophe, Edwards, Tamsin L., Felikson, Denis, Fettweis, Xavier, Galton‐Fenzi, Benjamin K., Goelzer, Heiko, Gladstone, Rupert, Golledge, Nicholas R., Gregory, Jonathan M., Greve, Ralf, Hattermann, Tore, Hoffman, Matthew J., Humbert, Angelika, Huybrechts, Philippe, Jourdain, Nicolas C., Kleiner, Thomas, Munneke, Peter Kuipers, Larour, Eric, Le clec'h, Sebastien, Lee, Victoria, Leguy, Gunter, Lipscomb, William H., Little, Christopher M., Lowry, Daniel P., Morlighem, Mathieu, Nias, Isabel, Pattyn, Frank, Pelle, Tyler, Price, Stephen F., Quiquet, Aurélien, Reese, Ronja, Rückamp, Martin, Schlegel, Nicole‐Jeanne, Seroussi, Hélène, Shepherd, Andrew, Simon, Erika, Slater, Donald, Smith, Robin S., Straneo, Fiammetta, Sun, Sainan, Tarasov, Lev, Trusel, Luke D., Van Breedam, Jonas, Wal, Roderik, Broeke, Michiel, Winkelmann, Ricarda, Zhao, Chen, Zhang, Tong, Zwinger, Thomas
Projections of the sea level contribution from the Greenland and Antarctic ice sheets (GrIS and AIS) rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous Coupled Model Intercomparison Project phase 5 (CMIP5) effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) – Part 1: Boundary conditions and climatic forcing
2020, Albrecht, Torsten, Winkelmann, Ricarda, Levermann, Anders
Simulations of the glacial–interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a “cookbook” for the growing community of PISM users and paleo-ice sheet modelers in general. For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
A simple equation for the melt elevation feedback of ice sheets
2016, Levermann, Anders, Winkelmann, Ricarda
In recent decades, the Greenland Ice Sheet has been losing mass and has thereby contributed to global sea-level rise. The rate of ice loss is highly relevant for coastal protection worldwide. The ice loss is likely to increase under future warming. Beyond a critical temperature threshold, a meltdown of the Greenland Ice Sheet is induced by the self-enforcing feedback between its lowering surface elevation and its increasing surface mass loss: the more ice that is lost, the lower the ice surface and the warmer the surface air temperature, which fosters further melting and ice loss. The computation of this rate so far relies on complex numerical models which are the appropriate tools for capturing the complexity of the problem. By contrast we aim here at gaining a conceptual understanding by deriving a purposefully simple equation for the self-enforcing feedback which is then used to estimate the melt time for different levels of warming using three observable characteristics of the ice sheet itself and its surroundings. The analysis is purely conceptual in nature. It is missing important processes like ice dynamics for it to be useful for applications to sea-level rise on centennial timescales, but if the volume loss is dominated by the feedback, the resulting logarithmic equation unifies existing numerical simulations and shows that the melt time depends strongly on the level of warming with a critical slowdown near the threshold: the median time to lose 10 % of the present-day ice volume varies between about 3500 years for a temperature level of 0.5 °C above the threshold and 500 years for 5 °C. Unless future observations show a significantly higher melting sensitivity than currently observed, a complete meltdown is unlikely within the next 2000 years without significant ice-dynamical contributions.