Articles | Volume 17, issue 16
https://doi.org/10.5194/gmd-17-6227-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-17-6227-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
22 Aug 2024
Development and technical paper |
| 22 Aug 2024
| 22 Aug 2024
Numerical stabilization methods for level-set-based ice front migration
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
Mathieu Morlighem
Department of Earth Sciences, Dartmouth College, Hanover, NH 03755, USA
G. Hilmar Gudmundsson
Department of Geography and Environmental Sciences, Northumbria University, Newcastle upon Tyne, NE1 8ST, UK
Related authors
Shfaqat A. Khan, Helene Seroussi, Mathieu Morlighem, William Colgan, Veit Helm, Gong Cheng, Danjal Berg, Valentina R. Barletta, Nicolaj K. Larsen, William Kochtitzky, Michiel van den Broeke, Kurt H. Kjær, Andy Aschwanden, Brice Noël, Jason E. Box, Joseph A. MacGregor, Robert S. Fausto, Kenneth D. Mankoff, Ian M. Howat, Kuba Oniszk, Dominik Fahrner, Anja Løkkegaard, Eigil Y. H. Lippert, and Javed Hassan
Earth Syst. Sci. Data Discuss., https://doi.org/10.5194/essd-2024-348, https://doi.org/10.5194/essd-2024-348, 2024
Revised manuscript under review for ESSD
Short summary
Short summary
The surface elevation of the Greenland Ice Sheet is changing due to surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from 2003 to 2023. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and useful for GIS ice sheet modelling.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
EGUsphere, https://doi.org/10.5194/egusphere-2024-1620, https://doi.org/10.5194/egusphere-2024-1620, 2024
Short summary
Short summary
Calving, the breaking of ice bodies from the terminus of a glacier, plays an important role in the mass losses of Greenland ice sheets. However, calving parameters have been poorly understood because of the intensive computational demands of traditional numerical models. To address this issue and find the optimal calving parameter that best represents real observations, we develop deep-learning emulators based on graph neural network architectures.
Joel A. Wilner, Mathieu Morlighem, and Gong Cheng
The Cryosphere, 17, 4889–4901, https://doi.org/10.5194/tc-17-4889-2023, https://doi.org/10.5194/tc-17-4889-2023, 2023
Short summary
Short summary
We use numerical modeling to study iceberg calving off of ice shelves in Antarctica. We examine four widely used mathematical descriptions of calving (
calving laws), under the assumption that Antarctic ice shelf front positions should be in steady state under the current climate forcing. We quantify how well each of these calving laws replicates the observed front positions. Our results suggest that the eigencalving and von Mises laws are most suitable for Antarctic ice shelves.
Gong Cheng, Nina Kirchner, and Per Lötstedt
The Cryosphere, 15, 715–742, https://doi.org/10.5194/tc-15-715-2021, https://doi.org/10.5194/tc-15-715-2021, 2021
Short summary
Short summary
We present an inverse modeling approach to improve the understanding of spatiotemporally variable processes at the inaccessible base of an ice sheet by determining the sensitivity of direct surface observations to perturbations of basal conditions. Time dependency is proved to be important in these types of problems. The effect of perturbations is analyzed based on analytical and numerical solutions.
Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 13, 2245–2258, https://doi.org/10.5194/gmd-13-2245-2020, https://doi.org/10.5194/gmd-13-2245-2020, 2020
Short summary
Short summary
A full Stokes subgrid scheme in two dimensions for the grounding line migration problem is presented in the open-source finite-element framework Elmer/ICE. This method can achieve comparable results to previous research using a more than 20 times larger mesh size, which can be used to improve the efficiency in marine ice sheet simulations.
Gong Cheng and Per Lötstedt
The Cryosphere, 14, 673–691, https://doi.org/10.5194/tc-14-673-2020, https://doi.org/10.5194/tc-14-673-2020, 2020
Short summary
Short summary
We present a time-dependent inverse method for ice sheet modeling. By investigating the sensitivity of the observations of the velocity and the height at the surface to the basal conditions of the ice, we show that if the basal parameters are time dependent, then time cannot be ignored in the inversion. By looking at the numerical features, we conclude that adding the height information of an ice sheet in the velocity inversion procedure could improve the robustness of the inference.
Sergey Marchenko, Gong Cheng, Per Lötstedt, Veijo Pohjola, Rickard Pettersson, Ward van Pelt, and Carleen Reijmer
The Cryosphere, 13, 1843–1859, https://doi.org/10.5194/tc-13-1843-2019, https://doi.org/10.5194/tc-13-1843-2019, 2019
Short summary
Short summary
Thermal conductivity (k) of firn at Lomonosovfonna, Svalbard, is estimated using measured temperature evolution and density. The optimized k values (0.2–1.6 W (m K)−1) increase downwards and over time and are most sensitive to systematic errors in measured temperature values and their depths, particularly in the lower part of the profile. Compared to the density-based parameterizations, derived k values are consistently larger, suggesting a faster conductive heat exchange in firn.
Eef C. H. van Dongen, Nina Kirchner, Martin B. van Gijzen, Roderik S. W. van de Wal, Thomas Zwinger, Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 11, 4563–4576, https://doi.org/10.5194/gmd-11-4563-2018, https://doi.org/10.5194/gmd-11-4563-2018, 2018
Short summary
Short summary
Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve. Therefore, approximations to the FS equations are used, especially when modeling an ice sheet on long time spans. Here, we report a combination of an approximation with the FS equations that allows simulating the dynamics of ice sheets over long time spans without introducing artifacts caused by application of approximations in parts of the domain where they are not valid.
Jowan M. Barnes, G. Hilmar Gudmundsson, Daniel N. Goldberg, and Sainan Sun
EGUsphere, https://doi.org/10.5194/egusphere-2025-328, https://doi.org/10.5194/egusphere-2025-328, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
Calving is where ice breaks off the front of glaciers. It has not been included widely in modelling as it is difficult to represent. We use our ice flow model to investigate the effects of calving floating ice shelves in West Antarctica. More calving leads to more ice loss and greater sea level rise, with local differences due to the shape of the bedrock. We find that ocean forcing and calving should be considered equally when trying to improve how models represent the real world.
Jamie Barnett, Felicity Alice Holmes, Joshua Cuzzone, Henning Åkesson, Mathieu Morlighem, Matt O'Regan, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2025-653, https://doi.org/10.5194/egusphere-2025-653, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
Understanding how ice sheets have changed in the past can allow us to make better predictions for the future. By running a state-of-the-art model of Ryder Glacier, North Greenland, over the past 12,000 years we find that both a warming atmosphere and ocean play a key role in the evolution of the Glacier. Our conclusions stress that accurately quantifying the ice sheet’s interactions with the ocean are required to predict future changes and reliable sea level rise estimates.
Francesca Baldacchino, Nicholas R. Golledge, Mathieu Morlighem, Huw Horgan, Alanna V. Alevropoulos-Borrill, Alena Malyarenko, Alexandra Gossart, Daniel P. Lowry, and Laurine van Haastrecht
The Cryosphere, 19, 107–127, https://doi.org/10.5194/tc-19-107-2025, https://doi.org/10.5194/tc-19-107-2025, 2025
Short summary
Short summary
Understanding how the Ross Ice Shelf flow is changing in a warming world is important for predicting ice sheet change. Field measurements show clear intra-annual variations in ice flow; however, it is unclear what mechanisms drive this variability. We show that local perturbations in basal melt can have a significant impact on ice flow speed, but a combination of forcings is likely driving the observed variability in ice flow.
Felicity Alice Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2024-3839, https://doi.org/10.5194/egusphere-2024-3839, 2024
Short summary
Short summary
Northern Greenland contains some of the ice sheet's last remaining glaciers with floating ice tongues. One of these is Ryder Glacier, which has been relatively stable in recent decades in contrast to nearby glaciers. Here, we use a computer model to simulate Ryder Glacier until 2300 under both a low and high emissions scenario. Very high levels of surface melt under a high emissions future leads to a sea-level rise contribution an order of magnitude higher than under a low emissions future.
Richard Parsons, Sainan Sun, G. Hilmar Gudmundsson, Jan Wuite, and Thomas Nagler
The Cryosphere, 18, 5789–5801, https://doi.org/10.5194/tc-18-5789-2024, https://doi.org/10.5194/tc-18-5789-2024, 2024
Short summary
Short summary
In 2022, multi-year landfast sea ice in Antarctica's Larsen B embayment disintegrated, after which time an increase in the rate at which Crane Glacier discharged ice into the ocean was observed. As the fast ice was joined to the glacier terminus, it could provide resistance against the glacier's flow, slowing down the rate of ice discharge. We used numerical modelling to quantify this resistive stress and found that the fast ice provided significant support to Crane prior to its disintegration.
Shfaqat A. Khan, Helene Seroussi, Mathieu Morlighem, William Colgan, Veit Helm, Gong Cheng, Danjal Berg, Valentina R. Barletta, Nicolaj K. Larsen, William Kochtitzky, Michiel van den Broeke, Kurt H. Kjær, Andy Aschwanden, Brice Noël, Jason E. Box, Joseph A. MacGregor, Robert S. Fausto, Kenneth D. Mankoff, Ian M. Howat, Kuba Oniszk, Dominik Fahrner, Anja Løkkegaard, Eigil Y. H. Lippert, and Javed Hassan
Earth Syst. Sci. Data Discuss., https://doi.org/10.5194/essd-2024-348, https://doi.org/10.5194/essd-2024-348, 2024
Revised manuscript under review for ESSD
Short summary
Short summary
The surface elevation of the Greenland Ice Sheet is changing due to surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from 2003 to 2023. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and useful for GIS ice sheet modelling.
Xianwei Wang, Hilmar Gudmundsson, and David Holland
EGUsphere, https://doi.org/10.5194/egusphere-2024-2790, https://doi.org/10.5194/egusphere-2024-2790, 2024
Short summary
Short summary
Understanding why iceberg calved during drifting in the ocean is important to understand the life cycle and the influence on the surrounding ocean of an iceberg. This study explains why iceberg A68a calved when approaching the South Georgia Island in late 2020 during its drifting in the Southern Ocean using satellite observation and modeling, which was caused by collision with seamount.
Brad Reed, J. A. Mattias Green, Adrian Jenkins, and G. Hilmar Gudmundsson
The Cryosphere, 18, 4567–4587, https://doi.org/10.5194/tc-18-4567-2024, https://doi.org/10.5194/tc-18-4567-2024, 2024
Short summary
Short summary
We use a numerical ice-flow model to simulate the response of a 1940s Pine Island Glacier to changes in melting beneath its ice shelf. A decadal period of warm forcing is sufficient to push the glacier into an unstable, irreversible retreat from its long-term position on a subglacial ridge to an upstream ice plain. This retreat can only be stopped when unrealistic cold forcing is applied. These results show that short warm anomalies can lead to quick and substantial increases in ice flux.
Joshua Cuzzone, Aaron Barth, Kelsey Barker, and Mathieu Morlighem
EGUsphere, https://doi.org/10.5194/egusphere-2024-2091, https://doi.org/10.5194/egusphere-2024-2091, 2024
Short summary
Short summary
We use an ice sheet model to simulate the Last Glacial Maximum conditions of the Laurentide Ice Sheet (LIS) across the Northeast United States. A complex thermal history existed for the (LIS), that allowed for high erosion across most of the NE USA, but prevented erosion across high elevation mountain peaks and areas where ice flow was slow. This has implications for geologic studies which rely on the erosional nature of the LIS to reconstruct its glacial history and landscape evolution.
Younghyun Koo, Gong Cheng, Mathieu Morlighem, and Maryam Rahnemoonfar
EGUsphere, https://doi.org/10.5194/egusphere-2024-1620, https://doi.org/10.5194/egusphere-2024-1620, 2024
Short summary
Short summary
Calving, the breaking of ice bodies from the terminus of a glacier, plays an important role in the mass losses of Greenland ice sheets. However, calving parameters have been poorly understood because of the intensive computational demands of traditional numerical models. To address this issue and find the optimal calving parameter that best represents real observations, we develop deep-learning emulators based on graph neural network architectures.
Sebastian H. R. Rosier, G. Hilmar Gudmundsson, Adrian Jenkins, and Kaitlin A. Naughten
EGUsphere, https://doi.org/10.5194/egusphere-2024-1838, https://doi.org/10.5194/egusphere-2024-1838, 2024
Short summary
Short summary
Glaciers in the Amundsen Sea region of Antarctica have been retreating and losing mass, but their future contribution to global sea level rise remains highly uncertain. We use an ice sheet model and uncertainty quantification methods to evaluate the probable range of mass loss from this region for two future climate scenarios and find that the rate of ice loss until 2100 will likely remain similar to present-day observations, with little sensitivity to climate scenario over this short timeframe.
J. Rachel Carr, Emily A. Hill, and G. Hilmar Gudmundsson
The Cryosphere, 18, 2719–2737, https://doi.org/10.5194/tc-18-2719-2024, https://doi.org/10.5194/tc-18-2719-2024, 2024
Short summary
Short summary
The Greenland Ice Sheet is one of the world's largest glaciers and is melting quickly in response to climate change. It contains fast-flowing channels of ice that move ice from Greenland's centre to its coasts and allow Greenland to react quickly to climate warming. As a result, we want to predict how these glaciers will behave in the future, but there are lots of uncertainties. Here we assess the impacts of two main sources of uncertainties in glacier models.
Cristina Gerli, Sebastian Rosier, G. Hilmar Gudmundsson, and Sainan Sun
The Cryosphere, 18, 2677–2689, https://doi.org/10.5194/tc-18-2677-2024, https://doi.org/10.5194/tc-18-2677-2024, 2024
Short summary
Short summary
Recent efforts have focused on using AI and satellite imagery to track crevasses for assessing ice shelf damage and informing ice flow models. Our study reveals a weak connection between these observed products and damage maps inferred from ice flow models. While there is some improvement in crevasse-dense regions, this association remains limited. Directly mapping ice damage from satellite observations may not significantly improve the representation of these processes within ice flow models.
Ole Richter, Ralph Timmermann, G. Hilmar Gudmundsson, and Jan De Rydt
EGUsphere, https://doi.org/10.5194/egusphere-2024-648, https://doi.org/10.5194/egusphere-2024-648, 2024
Short summary
Short summary
The new coupled ice sheet-ocean model addresses challenges related to horizontal resolution through advanced mesh flexibility, enabled by the use of unstructured grids. We describe the new model, verify its functioning in an idealised setting and demonstrate its advantages in a global-ocean/Antarctic ice sheet domain. The results of this study comprise an important step towards improving predictions of the Antarctic contribution to sea level rise over centennial time scales.
In-Woo Park, Emilia Kyung Jin, Mathieu Morlighem, and Kang-Kun Lee
The Cryosphere, 18, 1139–1155, https://doi.org/10.5194/tc-18-1139-2024, https://doi.org/10.5194/tc-18-1139-2024, 2024
Short summary
Short summary
This study conducted 3D thermodynamic ice sheet model experiments, and modeled temperatures were compared with 15 observed borehole temperature profiles. We found that using incompressibility of ice without sliding agrees well with observed temperature profiles in slow-flow regions, while incorporating sliding in fast-flow regions captures observed temperature profiles. Also, the choice of vertical velocity scheme has a greater impact on the shape of the modeled temperature profile.
Anjali Sandip, Ludovic Räss, and Mathieu Morlighem
Geosci. Model Dev., 17, 899–909, https://doi.org/10.5194/gmd-17-899-2024, https://doi.org/10.5194/gmd-17-899-2024, 2024
Short summary
Short summary
We solve momentum balance for unstructured meshes to predict ice flow for real glaciers using a pseudo-transient method on graphics processing units (GPUs) and compare it to a standard central processing unit (CPU) implementation. We justify the GPU implementation by applying the price-to-performance metric for up to million-grid-point spatial resolutions. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions.
Youngmin Choi, Helene Seroussi, Mathieu Morlighem, Nicole-Jeanne Schlegel, and Alex Gardner
The Cryosphere, 17, 5499–5517, https://doi.org/10.5194/tc-17-5499-2023, https://doi.org/10.5194/tc-17-5499-2023, 2023
Short summary
Short summary
Ice sheet models are often initialized using snapshot observations of present-day conditions, but this approach has limitations in capturing the transient evolution of the system. To more accurately represent the accelerating changes in glaciers, we employed time-dependent data assimilation. We found that models calibrated with the transient data better capture past trends and more accurately reproduce changes after the calibration period, even with limited observations.
Hélène Seroussi, Vincent Verjans, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Peter Van Katwyk, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 17, 5197–5217, https://doi.org/10.5194/tc-17-5197-2023, https://doi.org/10.5194/tc-17-5197-2023, 2023
Short summary
Short summary
Mass loss from Antarctica is a key contributor to sea level rise over the 21st century, and the associated uncertainty dominates sea level projections. We highlight here the Antarctic glaciers showing the largest changes and quantify the main sources of uncertainty in their future evolution using an ensemble of ice flow models. We show that on top of Pine Island and Thwaites glaciers, Totten and Moscow University glaciers show rapid changes and a strong sensitivity to warmer ocean conditions.
Joel A. Wilner, Mathieu Morlighem, and Gong Cheng
The Cryosphere, 17, 4889–4901, https://doi.org/10.5194/tc-17-4889-2023, https://doi.org/10.5194/tc-17-4889-2023, 2023
Short summary
Short summary
We use numerical modeling to study iceberg calving off of ice shelves in Antarctica. We examine four widely used mathematical descriptions of calving (
calving laws), under the assumption that Antarctic ice shelf front positions should be in steady state under the current climate forcing. We quantify how well each of these calving laws replicates the observed front positions. Our results suggest that the eigencalving and von Mises laws are most suitable for Antarctic ice shelves.
Emily A. Hill, Benoît Urruty, Ronja Reese, Julius Garbe, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Ricarda Winkelmann, Mondher Chekki, David Chandler, and Petra M. Langebroek
The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, https://doi.org/10.5194/tc-17-3739-2023, 2023
Short summary
Short summary
The grounding lines of the Antarctic Ice Sheet could enter phases of irreversible retreat or advance. We use three ice sheet models to show that the present-day locations of Antarctic grounding lines are reversible with respect to a small perturbation away from their current position. This indicates that present-day retreat of the grounding lines is not yet irreversible or self-enhancing.
Ronja Reese, Julius Garbe, Emily A. Hill, Benoît Urruty, Kaitlin A. Naughten, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, David Chandler, Petra M. Langebroek, and Ricarda Winkelmann
The Cryosphere, 17, 3761–3783, https://doi.org/10.5194/tc-17-3761-2023, https://doi.org/10.5194/tc-17-3761-2023, 2023
Short summary
Short summary
We use an ice sheet model to test where current climate conditions in Antarctica might lead. We find that present-day ocean and atmosphere conditions might commit an irreversible collapse of parts of West Antarctica which evolves over centuries to millennia. Importantly, this collapse is not irreversible yet.
Alice C. Frémand, Peter Fretwell, Julien A. Bodart, Hamish D. Pritchard, Alan Aitken, Jonathan L. Bamber, Robin Bell, Cesidio Bianchi, Robert G. Bingham, Donald D. Blankenship, Gino Casassa, Ginny Catania, Knut Christianson, Howard Conway, Hugh F. J. Corr, Xiangbin Cui, Detlef Damaske, Volkmar Damm, Reinhard Drews, Graeme Eagles, Olaf Eisen, Hannes Eisermann, Fausto Ferraccioli, Elena Field, René Forsberg, Steven Franke, Shuji Fujita, Yonggyu Gim, Vikram Goel, Siva Prasad Gogineni, Jamin Greenbaum, Benjamin Hills, Richard C. A. Hindmarsh, Andrew O. Hoffman, Per Holmlund, Nicholas Holschuh, John W. Holt, Annika N. Horlings, Angelika Humbert, Robert W. Jacobel, Daniela Jansen, Adrian Jenkins, Wilfried Jokat, Tom Jordan, Edward King, Jack Kohler, William Krabill, Mette Kusk Gillespie, Kirsty Langley, Joohan Lee, German Leitchenkov, Carlton Leuschen, Bruce Luyendyk, Joseph MacGregor, Emma MacKie, Kenichi Matsuoka, Mathieu Morlighem, Jérémie Mouginot, Frank O. Nitsche, Yoshifumi Nogi, Ole A. Nost, John Paden, Frank Pattyn, Sergey V. Popov, Eric Rignot, David M. Rippin, Andrés Rivera, Jason Roberts, Neil Ross, Anotonia Ruppel, Dustin M. Schroeder, Martin J. Siegert, Andrew M. Smith, Daniel Steinhage, Michael Studinger, Bo Sun, Ignazio Tabacco, Kirsty Tinto, Stefano Urbini, David Vaughan, Brian C. Welch, Douglas S. Wilson, Duncan A. Young, and Achille Zirizzotti
Earth Syst. Sci. Data, 15, 2695–2710, https://doi.org/10.5194/essd-15-2695-2023, https://doi.org/10.5194/essd-15-2695-2023, 2023
Short summary
Short summary
This paper presents the release of over 60 years of ice thickness, bed elevation, and surface elevation data acquired over Antarctica by the international community. These data are a crucial component of the Antarctic Bedmap initiative which aims to produce a new map and datasets of Antarctic ice thickness and bed topography for the international glaciology and geophysical community.
Sebastian H. R. Rosier, Christopher Y. S. Bull, Wai L. Woo, and G. Hilmar Gudmundsson
The Cryosphere, 17, 499–518, https://doi.org/10.5194/tc-17-499-2023, https://doi.org/10.5194/tc-17-499-2023, 2023
Short summary
Short summary
Future ice loss from Antarctica could raise sea levels by several metres, and key to this is the rate at which the ocean melts the ice sheet from below. Existing methods for modelling this process are either computationally expensive or very simplified. We present a new approach using machine learning to mimic the melt rates calculated by an ocean model but in a fraction of the time. This approach may provide a powerful alternative to existing methods, without compromising on accuracy or speed.
Bertie W. J. Miles, Chris R. Stokes, Adrian Jenkins, Jim R. Jordan, Stewart S. R. Jamieson, and G. Hilmar Gudmundsson
The Cryosphere, 17, 445–456, https://doi.org/10.5194/tc-17-445-2023, https://doi.org/10.5194/tc-17-445-2023, 2023
Short summary
Short summary
Satellite observations have shown that the Shirase Glacier catchment in East Antarctica has been gaining mass over the past 2 decades, a trend largely attributed to increased snowfall. Our multi-decadal observations of Shirase Glacier show that ocean forcing has also contributed to some of this recent mass gain. This has been caused by strengthening easterly winds reducing the inflow of warm water underneath the Shirase ice tongue, causing the glacier to slow down and thicken.
Jowan M. Barnes and G. Hilmar Gudmundsson
The Cryosphere, 16, 4291–4304, https://doi.org/10.5194/tc-16-4291-2022, https://doi.org/10.5194/tc-16-4291-2022, 2022
Short summary
Short summary
Models must represent how glaciers slide along the bed, but there are many ways to do so. In this paper, several sliding laws are tested and found to affect different regions of the Antarctic Ice Sheet in different ways and at different speeds. However, the variability in ice volume loss due to sliding-law choices is low compared to other factors, so limited empirical knowledge of sliding does not prevent us from making predictions of how an ice sheet will evolve.
Francesca Baldacchino, Mathieu Morlighem, Nicholas R. Golledge, Huw Horgan, and Alena Malyarenko
The Cryosphere, 16, 3723–3738, https://doi.org/10.5194/tc-16-3723-2022, https://doi.org/10.5194/tc-16-3723-2022, 2022
Short summary
Short summary
Understanding how the Ross Ice Shelf will evolve in a warming world is important to the future stability of Antarctica. It remains unclear what changes could drive the largest mass loss in the future and where places are most likely to trigger larger mass losses. Sensitivity maps are modelled showing that the RIS is sensitive to changes in environmental and glaciological controls at regions which are currently experiencing changes. These regions need to be monitored in a warming world.
Joshua K. Cuzzone, Nicolás E. Young, Mathieu Morlighem, Jason P. Briner, and Nicole-Jeanne Schlegel
The Cryosphere, 16, 2355–2372, https://doi.org/10.5194/tc-16-2355-2022, https://doi.org/10.5194/tc-16-2355-2022, 2022
Short summary
Short summary
We use an ice sheet model to determine what influenced the Greenland Ice Sheet to retreat across a portion of southwestern Greenland during the Holocene (about the last 12 000 years). Our simulations, constrained by observations from geologic markers, show that atmospheric warming and ice melt primarily caused the ice sheet to retreat rapidly across this domain. We find, however, that iceberg calving at the interface where the ice meets the ocean significantly influenced ice mass change.
Yannic Fischler, Martin Rückamp, Christian Bischof, Vadym Aizinger, Mathieu Morlighem, and Angelika Humbert
Geosci. Model Dev., 15, 3753–3771, https://doi.org/10.5194/gmd-15-3753-2022, https://doi.org/10.5194/gmd-15-3753-2022, 2022
Short summary
Short summary
Ice sheet models are used to simulate the changes of ice sheets in future but are currently often run in coarse resolution and/or with neglecting important physics to make them affordable in terms of computational costs. We conducted a study simulating the Greenland Ice Sheet in high resolution and adequate physics to test where the ISSM ice sheet code is using most time and what could be done to improve its performance for future computer architectures that allow massive parallel computing.
Tom Mitcham, G. Hilmar Gudmundsson, and Jonathan L. Bamber
The Cryosphere, 16, 883–901, https://doi.org/10.5194/tc-16-883-2022, https://doi.org/10.5194/tc-16-883-2022, 2022
Short summary
Short summary
We modelled the response of the Larsen C Ice Shelf (LCIS) and its tributary glaciers to the calving of the A68 iceberg and validated our results with observations. We found that the impact was limited, confirming that mostly passive ice was calved. Through further calving experiments we quantified the total buttressing provided by the LCIS and found that over 80 % of the buttressing capacity is generated in the first 5 km of the ice shelf downstream of the grounding line.
Thomas Frank, Henning Åkesson, Basile de Fleurian, Mathieu Morlighem, and Kerim H. Nisancioglu
The Cryosphere, 16, 581–601, https://doi.org/10.5194/tc-16-581-2022, https://doi.org/10.5194/tc-16-581-2022, 2022
Short summary
Short summary
The shape of a fjord can promote or inhibit glacier retreat in response to climate change. We conduct experiments with a synthetic setup under idealized conditions in a numerical model to study and quantify the processes involved. We find that friction between ice and fjord is the most important factor and that it is possible to directly link ice discharge and grounding line retreat to fjord topography in a quantitative way.
Thiago Dias dos Santos, Mathieu Morlighem, and Douglas Brinkerhoff
The Cryosphere, 16, 179–195, https://doi.org/10.5194/tc-16-179-2022, https://doi.org/10.5194/tc-16-179-2022, 2022
Short summary
Short summary
Projecting the future evolution of Greenland and Antarctica and their potential contribution to sea level rise often relies on computer simulations carried out by numerical ice sheet models. Here we present a new vertically integrated ice sheet model and assess its performance using different benchmarks. The new model shows results comparable to a three-dimensional model at relatively lower computational cost, suggesting that it is an excellent alternative for long-term simulations.
Emily A. Hill, Sebastian H. R. Rosier, G. Hilmar Gudmundsson, and Matthew Collins
The Cryosphere, 15, 4675–4702, https://doi.org/10.5194/tc-15-4675-2021, https://doi.org/10.5194/tc-15-4675-2021, 2021
Short summary
Short summary
Using an ice flow model and uncertainty quantification methods, we provide probabilistic projections of future sea level rise from the Filchner–Ronne region of Antarctica. We find that it is most likely that this region will contribute negatively to sea level rise over the next 300 years, largely as a result of increased surface mass balance. We identify parameters controlling ice shelf melt and snowfall contribute most to uncertainties in projections.
Matt O'Regan, Thomas M. Cronin, Brendan Reilly, Aage Kristian Olsen Alstrup, Laura Gemery, Anna Golub, Larry A. Mayer, Mathieu Morlighem, Matthias Moros, Ole L. Munk, Johan Nilsson, Christof Pearce, Henrieka Detlef, Christian Stranne, Flor Vermassen, Gabriel West, and Martin Jakobsson
The Cryosphere, 15, 4073–4097, https://doi.org/10.5194/tc-15-4073-2021, https://doi.org/10.5194/tc-15-4073-2021, 2021
Short summary
Short summary
Ryder Glacier is a marine-terminating glacier in north Greenland discharging ice into the Lincoln Sea. Here we use marine sediment cores to reconstruct its retreat and advance behavior through the Holocene. We show that while Sherard Osborn Fjord has a physiography conducive to glacier and ice tongue stability, Ryder still retreated more than 40 km inland from its current position by the Middle Holocene. This highlights the sensitivity of north Greenland's marine glaciers to climate change.
Thiago Dias dos Santos, Mathieu Morlighem, and Hélène Seroussi
Geosci. Model Dev., 14, 2545–2573, https://doi.org/10.5194/gmd-14-2545-2021, https://doi.org/10.5194/gmd-14-2545-2021, 2021
Short summary
Short summary
Numerical models are routinely used to understand the past and future behavior of ice sheets in response to climate evolution. As is always the case with numerical modeling, one needs to minimize biases and numerical artifacts due to the choice of numerical scheme employed in such models. Here, we assess different numerical schemes in time-dependent simulations of ice sheets. We also introduce a new parameterization for the driving stress, the force that drives the ice sheet flow.
Jowan M. Barnes, Thiago Dias dos Santos, Daniel Goldberg, G. Hilmar Gudmundsson, Mathieu Morlighem, and Jan De Rydt
The Cryosphere, 15, 1975–2000, https://doi.org/10.5194/tc-15-1975-2021, https://doi.org/10.5194/tc-15-1975-2021, 2021
Short summary
Short summary
Some properties of ice flow models must be initialised using observed data before they can be used to produce reliable predictions of the future. Different models have different ways of doing this, and the process is generally seen as being specific to an individual model. We compare the methods used by three different models and show that they produce similar outputs. We also demonstrate that the outputs from one model can be used in other models without introducing large uncertainties.
Sebastian H. R. Rosier, Ronja Reese, Jonathan F. Donges, Jan De Rydt, G. Hilmar Gudmundsson, and Ricarda Winkelmann
The Cryosphere, 15, 1501–1516, https://doi.org/10.5194/tc-15-1501-2021, https://doi.org/10.5194/tc-15-1501-2021, 2021
Short summary
Short summary
Pine Island Glacier has contributed more to sea-level rise over the past decades than any other glacier in Antarctica. Ice-flow modelling studies have shown that it can undergo periods of rapid mass loss, but no study has shown that these future changes could cross a tipping point and therefore be effectively irreversible. Here, we assess the stability of Pine Island Glacier, quantifying the changes in ocean temperatures required to cross future tipping points using statistical methods.
Gong Cheng, Nina Kirchner, and Per Lötstedt
The Cryosphere, 15, 715–742, https://doi.org/10.5194/tc-15-715-2021, https://doi.org/10.5194/tc-15-715-2021, 2021
Short summary
Short summary
We present an inverse modeling approach to improve the understanding of spatiotemporally variable processes at the inaccessible base of an ice sheet by determining the sensitivity of direct surface observations to perturbations of basal conditions. Time dependency is proved to be important in these types of problems. The effect of perturbations is analyzed based on analytical and numerical solutions.
Bertie W. J. Miles, Jim R. Jordan, Chris R. Stokes, Stewart S. R. Jamieson, G. Hilmar Gudmundsson, and Adrian Jenkins
The Cryosphere, 15, 663–676, https://doi.org/10.5194/tc-15-663-2021, https://doi.org/10.5194/tc-15-663-2021, 2021
Short summary
Short summary
We provide a historical overview of changes in Denman Glacier's flow speed, structure and calving events since the 1960s. Based on these observations, we perform a series of numerical modelling experiments to determine the likely cause of Denman's acceleration since the 1970s. We show that grounding line retreat, ice shelf thinning and the detachment of Denman's ice tongue from a pinning point are the most likely causes of the observed acceleration.
Jan De Rydt, Ronja Reese, Fernando S. Paolo, and G. Hilmar Gudmundsson
The Cryosphere, 15, 113–132, https://doi.org/10.5194/tc-15-113-2021, https://doi.org/10.5194/tc-15-113-2021, 2021
Short summary
Short summary
We used satellite observations and numerical simulations of Pine Island Glacier, West Antarctica, between 1996 and 2016 to show that the recent increase in its flow speed can only be reproduced by computer models if stringent assumptions are made about the material properties of the ice and its underlying bed. These assumptions are not commonly adopted in ice flow modelling, and our results therefore have implications for future simulations of Antarctic ice flow and sea level projections.
Kate Winter, Emily A. Hill, G. Hilmar Gudmundsson, and John Woodward
Earth Syst. Sci. Data, 12, 3453–3467, https://doi.org/10.5194/essd-12-3453-2020, https://doi.org/10.5194/essd-12-3453-2020, 2020
Short summary
Short summary
Satellite measurements of the English Coast in the Antarctic Peninsula reveal that glaciers are thinning and losing mass, but ice thickness data are required to assess these changes, in terms of ice flux and sea level contribution. Our ice-penetrating radar measurements reveal that low-elevation subglacial channels control fast-flowing ice streams, which release over 39 Gt of ice per year to floating ice shelves. This topography could make ice flows susceptible to future instability.
Xiangbin Cui, Hafeez Jeofry, Jamin S. Greenbaum, Jingxue Guo, Lin Li, Laura E. Lindzey, Feras A. Habbal, Wei Wei, Duncan A. Young, Neil Ross, Mathieu Morlighem, Lenneke M. Jong, Jason L. Roberts, Donald D. Blankenship, Sun Bo, and Martin J. Siegert
Earth Syst. Sci. Data, 12, 2765–2774, https://doi.org/10.5194/essd-12-2765-2020, https://doi.org/10.5194/essd-12-2765-2020, 2020
Short summary
Short summary
We present a topographic digital elevation model (DEM) for Princess Elizabeth Land (PEL), East Antarctica. The DEM covers an area of approximately 900 000 km2 and was built from radio-echo sounding data collected in four campaigns since 2015. Previously, to generate the Bedmap2 topographic product, PEL’s bed was characterised from low-resolution satellite gravity data across an otherwise large (>200 km wide) data-free zone.
Eric Larour, Lambert Caron, Mathieu Morlighem, Surendra Adhikari, Thomas Frederikse, Nicole-Jeanne Schlegel, Erik Ivins, Benjamin Hamlington, Robert Kopp, and Sophie Nowicki
Geosci. Model Dev., 13, 4925–4941, https://doi.org/10.5194/gmd-13-4925-2020, https://doi.org/10.5194/gmd-13-4925-2020, 2020
Short summary
Short summary
ISSM-SLPS is a new projection system for future sea level that increases the resolution and accuracy of current projection systems and improves the way uncertainty is treated in such projections. This will pave the way for better inclusion of state-of-the-art results from existing intercomparison efforts carried out by the scientific community, such as GlacierMIP2 or ISMIP6, into sea-level projections.
Martin Rückamp, Angelika Humbert, Thomas Kleiner, Mathieu Morlighem, and Helene Seroussi
Geosci. Model Dev., 13, 4491–4501, https://doi.org/10.5194/gmd-13-4491-2020, https://doi.org/10.5194/gmd-13-4491-2020, 2020
Short summary
Short summary
We present enthalpy formulations within the Ice-Sheet and Sea-Level System model that show better performance than earlier implementations. A first experiment indicates that the treatment of discontinuous conductivities of the solid–fluid system with a geometric mean produce accurate results when applied to coarse vertical resolutions. In a second experiment, we propose a novel stabilization formulation that avoids the problem of thin elements. This method provides accurate and stable results.
Heiko Goelzer, Sophie Nowicki, Anthony Payne, Eric Larour, Helene Seroussi, William H. Lipscomb, Jonathan Gregory, Ayako Abe-Ouchi, Andrew Shepherd, Erika Simon, Cécile Agosta, Patrick Alexander, Andy Aschwanden, Alice Barthel, Reinhard Calov, Christopher Chambers, Youngmin Choi, Joshua Cuzzone, Christophe Dumas, Tamsin Edwards, Denis Felikson, Xavier Fettweis, Nicholas R. Golledge, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Sebastien Le clec'h, Victoria Lee, Gunter Leguy, Chris Little, Daniel P. Lowry, Mathieu Morlighem, Isabel Nias, Aurelien Quiquet, Martin Rückamp, Nicole-Jeanne Schlegel, Donald A. Slater, Robin S. Smith, Fiamma Straneo, Lev Tarasov, Roderik van de Wal, and Michiel van den Broeke
The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, https://doi.org/10.5194/tc-14-3071-2020, 2020
Short summary
Short summary
In this paper we use a large ensemble of Greenland ice sheet models forced by six different global climate models to project ice sheet changes and sea-level rise contributions over the 21st century.
The results for two different greenhouse gas concentration scenarios indicate that the Greenland ice sheet will continue to lose mass until 2100, with contributions to sea-level rise of 90 ± 50 mm and 32 ± 17 mm for the high (RCP8.5) and low (RCP2.6) scenario, respectively.
Hélène Seroussi, Sophie Nowicki, Antony J. Payne, Heiko Goelzer, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Torsten Albrecht, Xylar Asay-Davis, Alice Barthel, Reinhard Calov, Richard Cullather, Christophe Dumas, Benjamin K. Galton-Fenzi, Rupert Gladstone, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Tore Hattermann, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Nicolas C. Jourdain, Thomas Kleiner, Eric Larour, Gunter R. Leguy, Daniel P. Lowry, Chistopher M. Little, Mathieu Morlighem, Frank Pattyn, Tyler Pelle, Stephen F. Price, Aurélien Quiquet, Ronja Reese, Nicole-Jeanne Schlegel, Andrew Shepherd, Erika Simon, Robin S. Smith, Fiammetta Straneo, Sainan Sun, Luke D. Trusel, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, Chen Zhao, Tong Zhang, and Thomas Zwinger
The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, https://doi.org/10.5194/tc-14-3033-2020, 2020
Short summary
Short summary
The Antarctic ice sheet has been losing mass over at least the past 3 decades in response to changes in atmospheric and oceanic conditions. This study presents an ensemble of model simulations of the Antarctic evolution over the 2015–2100 period based on various ice sheet models, climate forcings and emission scenarios. Results suggest that the West Antarctic ice sheet will continue losing a large amount of ice, while the East Antarctic ice sheet could experience increased snow accumulation.
Sophie Nowicki, Heiko Goelzer, Hélène Seroussi, Anthony J. Payne, William H. Lipscomb, Ayako Abe-Ouchi, Cécile Agosta, Patrick Alexander, Xylar S. Asay-Davis, Alice Barthel, Thomas J. Bracegirdle, Richard Cullather, Denis Felikson, Xavier Fettweis, Jonathan M. Gregory, Tore Hattermann, Nicolas C. Jourdain, Peter Kuipers Munneke, Eric Larour, Christopher M. Little, Mathieu Morlighem, Isabel Nias, Andrew Shepherd, Erika Simon, Donald Slater, Robin S. Smith, Fiammetta Straneo, Luke D. Trusel, Michiel R. van den Broeke, and Roderik van de Wal
The Cryosphere, 14, 2331–2368, https://doi.org/10.5194/tc-14-2331-2020, https://doi.org/10.5194/tc-14-2331-2020, 2020
Short summary
Short summary
This paper describes the experimental protocol for ice sheet models taking part in the Ice Sheet Model Intercomparion Project for CMIP6 (ISMIP6) and presents an overview of the atmospheric and oceanic datasets to be used for the simulations. The ISMIP6 framework allows for exploring the uncertainty in 21st century sea level change from the Greenland and Antarctic ice sheets.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 13, 2245–2258, https://doi.org/10.5194/gmd-13-2245-2020, https://doi.org/10.5194/gmd-13-2245-2020, 2020
Short summary
Short summary
A full Stokes subgrid scheme in two dimensions for the grounding line migration problem is presented in the open-source finite-element framework Elmer/ICE. This method can achieve comparable results to previous research using a more than 20 times larger mesh size, which can be used to improve the efficiency in marine ice sheet simulations.
Donald A. Slater, Denis Felikson, Fiamma Straneo, Heiko Goelzer, Christopher M. Little, Mathieu Morlighem, Xavier Fettweis, and Sophie Nowicki
The Cryosphere, 14, 985–1008, https://doi.org/10.5194/tc-14-985-2020, https://doi.org/10.5194/tc-14-985-2020, 2020
Short summary
Short summary
Changes in the ocean around Greenland play an important role in determining how much the ice sheet will contribute to global sea level over the coming century. However, capturing these links in models is very challenging. This paper presents a strategy enabling an ensemble of ice sheet models to feel the effect of the ocean for the first time and should therefore result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.
Gong Cheng and Per Lötstedt
The Cryosphere, 14, 673–691, https://doi.org/10.5194/tc-14-673-2020, https://doi.org/10.5194/tc-14-673-2020, 2020
Short summary
Short summary
We present a time-dependent inverse method for ice sheet modeling. By investigating the sensitivity of the observations of the velocity and the height at the surface to the basal conditions of the ice, we show that if the basal parameters are time dependent, then time cannot be ignored in the inversion. By looking at the numerical features, we conclude that adding the height information of an ice sheet in the velocity inversion procedure could improve the robustness of the inference.
Anders Levermann, Ricarda Winkelmann, Torsten Albrecht, Heiko Goelzer, Nicholas R. Golledge, Ralf Greve, Philippe Huybrechts, Jim Jordan, Gunter Leguy, Daniel Martin, Mathieu Morlighem, Frank Pattyn, David Pollard, Aurelien Quiquet, Christian Rodehacke, Helene Seroussi, Johannes Sutter, Tong Zhang, Jonas Van Breedam, Reinhard Calov, Robert DeConto, Christophe Dumas, Julius Garbe, G. Hilmar Gudmundsson, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, William H. Lipscomb, Malte Meinshausen, Esmond Ng, Sophie M. J. Nowicki, Mauro Perego, Stephen F. Price, Fuyuki Saito, Nicole-Jeanne Schlegel, Sainan Sun, and Roderik S. W. van de Wal
Earth Syst. Dynam., 11, 35–76, https://doi.org/10.5194/esd-11-35-2020, https://doi.org/10.5194/esd-11-35-2020, 2020
Short summary
Short summary
We provide an estimate of the future sea level contribution of Antarctica from basal ice shelf melting up to the year 2100. The full uncertainty range in the warming-related forcing of basal melt is estimated and applied to 16 state-of-the-art ice sheet models using a linear response theory approach. The sea level contribution we obtain is very likely below 61 cm under unmitigated climate change until 2100 (RCP8.5) and very likely below 40 cm if the Paris Climate Agreement is kept.
Sebastian H. R. Rosier and G. Hilmar Gudmundsson
The Cryosphere, 14, 17–37, https://doi.org/10.5194/tc-14-17-2020, https://doi.org/10.5194/tc-14-17-2020, 2020
Short summary
Short summary
The flow of ice shelves is now known to be strongly affected by ocean tides, but the mechanism by which this happens is unclear. We use a viscoelastic model to try to reproduce observations of this behaviour on the Filchner–Ronne Ice Shelf in Antarctica. We find that tilting of the ice shelf explains the short-period behaviour, while tidally induced movement of the grounding line (the boundary between grounded and floating ice) explains the more complex long-period response.
Jan De Rydt, Gudmundur Hilmar Gudmundsson, Thomas Nagler, and Jan Wuite
The Cryosphere, 13, 2771–2787, https://doi.org/10.5194/tc-13-2771-2019, https://doi.org/10.5194/tc-13-2771-2019, 2019
Short summary
Short summary
Two large icebergs are about to break off from the Brunt Ice Shelf in Antarctica. Rifting started several years ago and is now approaching its final phase. Satellite data and computer simulations show that over the past 2 decades, growth of the ice shelf has caused a build-up of forces within the ice, which culminated in its fracture. These natural changes in geometry coincided with large variations in flow speed, a process that is thought to be relevant for all Antarctic ice shelf margins.
Sergey Marchenko, Gong Cheng, Per Lötstedt, Veijo Pohjola, Rickard Pettersson, Ward van Pelt, and Carleen Reijmer
The Cryosphere, 13, 1843–1859, https://doi.org/10.5194/tc-13-1843-2019, https://doi.org/10.5194/tc-13-1843-2019, 2019
Short summary
Short summary
Thermal conductivity (k) of firn at Lomonosovfonna, Svalbard, is estimated using measured temperature evolution and density. The optimized k values (0.2–1.6 W (m K)−1) increase downwards and over time and are most sensitive to systematic errors in measured temperature values and their depths, particularly in the lower part of the profile. Compared to the density-based parameterizations, derived k values are consistently larger, suggesting a faster conductive heat exchange in firn.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Tyler Pelle, Mathieu Morlighem, and Johannes H. Bondzio
The Cryosphere, 13, 1043–1049, https://doi.org/10.5194/tc-13-1043-2019, https://doi.org/10.5194/tc-13-1043-2019, 2019
Short summary
Short summary
How ocean-induced melt under floating ice shelves will change as ocean currents evolve remains a big uncertainty in projections of sea level rise. In this study, we combine two of the most recently developed melt models to form PICOP, which overcomes the limitations of past models and produces accurate ice shelf melt rates. We find that our model is easy to set up and computationally efficient, providing researchers an important tool to improve the accuracy of their future glacial projections.
Joshua K. Cuzzone, Nicole-Jeanne Schlegel, Mathieu Morlighem, Eric Larour, Jason P. Briner, Helene Seroussi, and Lambert Caron
The Cryosphere, 13, 879–893, https://doi.org/10.5194/tc-13-879-2019, https://doi.org/10.5194/tc-13-879-2019, 2019
Short summary
Short summary
We present ice sheet modeling results of ice retreat over southwestern Greenland during the last 12 000 years, and we also test the impact that model horizontal resolution has on differences in the simulated spatial retreat and its associated rate. Results indicate that model resolution plays a minor role in simulated retreat in areas where bed topography is not complex but plays an important role in areas where bed topography is complex (such as fjords).
Mathieu Morlighem, Michael Wood, Hélène Seroussi, Youngmin Choi, and Eric Rignot
The Cryosphere, 13, 723–734, https://doi.org/10.5194/tc-13-723-2019, https://doi.org/10.5194/tc-13-723-2019, 2019
Short summary
Short summary
Many glaciers along the coast of Greenland have been retreating. It has been suggested that this retreat is triggered by the presence of warm water in the fjords, and surface melt at the top of the ice sheet is exacerbating this problem. Here, we quantify the vulnerability of northwestern Greenland to further warming using a numerical model. We find that in current conditions, this sector alone will contribute more than 1 cm to sea rise level by 2100, and up to 3 cm in the most extreme scenario.
Thiago Dias dos Santos, Mathieu Morlighem, Hélène Seroussi, Philippe Remy Bernard Devloo, and Jefferson Cardia Simões
Geosci. Model Dev., 12, 215–232, https://doi.org/10.5194/gmd-12-215-2019, https://doi.org/10.5194/gmd-12-215-2019, 2019
Short summary
Short summary
The reduction of numerical errors in ice sheet modeling increases the results' accuracy reliability. We improve numerical accuracy by better capturing grounding line dynamics, while maintaining a low computational cost. We implement an adaptive mesh refinement (AMR) technique in the Ice Sheet System Model and compare AMR simulations with uniformly refined meshes. Our results show that the computational time with AMR is significantly shorter than for uniformly refined meshes for a given accuracy.
Emily A. Hill, G. Hilmar Gudmundsson, J. Rachel Carr, and Chris R. Stokes
The Cryosphere, 12, 3907–3921, https://doi.org/10.5194/tc-12-3907-2018, https://doi.org/10.5194/tc-12-3907-2018, 2018
Short summary
Short summary
Floating ice tongues in Greenland buttress inland ice, and their removal could accelerate ice flow. Petermann Glacier recently lost large sections of its ice tongue, but there was little glacier acceleration. Here, we assess the impact of future calving events on ice speeds. We find that removing the lower portions of the ice tongue does not accelerate flow. However, future iceberg calving closer to the grounding line could accelerate ice flow and increase ice discharge and sea level rise.
Hongju Yu, Eric Rignot, Helene Seroussi, and Mathieu Morlighem
The Cryosphere, 12, 3861–3876, https://doi.org/10.5194/tc-12-3861-2018, https://doi.org/10.5194/tc-12-3861-2018, 2018
Short summary
Short summary
Thwaites Glacier, West Antarctica, has experienced rapid grounding line retreat and mass loss in the past decades. In this study, we simulate the evolution of Thwaites Glacier over the next century using different model configurations. Overall, we estimate a 5 mm contribution to global sea level rise from Thwaites Glacier in the next 30 years. However, a 300 % uncertainty is found over the next 100 years, ranging from 14 to 42 mm, depending on the model setup.
Youngmin Choi, Mathieu Morlighem, Michael Wood, and Johannes H. Bondzio
The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, https://doi.org/10.5194/tc-12-3735-2018, 2018
Short summary
Short summary
Calving is an important mechanism that controls the dynamics of Greenland outlet glaciers. We test and compare four calving laws and assess which calving law has better predictive abilities. Overall, the calving law based on von Mises stress is more satisfactory than other laws, but new parameterizations should be derived to better capture the detailed processes involved in calving.
Eef C. H. van Dongen, Nina Kirchner, Martin B. van Gijzen, Roderik S. W. van de Wal, Thomas Zwinger, Gong Cheng, Per Lötstedt, and Lina von Sydow
Geosci. Model Dev., 11, 4563–4576, https://doi.org/10.5194/gmd-11-4563-2018, https://doi.org/10.5194/gmd-11-4563-2018, 2018
Short summary
Short summary
Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve. Therefore, approximations to the FS equations are used, especially when modeling an ice sheet on long time spans. Here, we report a combination of an approximation with the FS equations that allows simulating the dynamics of ice sheets over long time spans without introducing artifacts caused by application of approximations in parts of the domain where they are not valid.
Nicole-Jeanne Schlegel, Helene Seroussi, Michael P. Schodlok, Eric Y. Larour, Carmen Boening, Daniel Limonadi, Michael M. Watkins, Mathieu Morlighem, and Michiel R. van den Broeke
The Cryosphere, 12, 3511–3534, https://doi.org/10.5194/tc-12-3511-2018, https://doi.org/10.5194/tc-12-3511-2018, 2018
Short summary
Short summary
Using NASA supercomputers and a novel framework, in which Sandia National Laboratories' statistical software is embedded in the Jet Propulsion Laboratory's ice sheet model, we run a range of 100-year warming scenarios for Antarctica. We find that 1.2 m of sea level contribution is achievable, but not likely. Also, we find that bedrock topography beneath the ice drives potential for regional sea level contribution, highlighting the need for accurate bedrock mapping of the ice sheet interior.
Edward C. King, Jan De Rydt, and G. Hilmar Gudmundsson
The Cryosphere, 12, 3361–3372, https://doi.org/10.5194/tc-12-3361-2018, https://doi.org/10.5194/tc-12-3361-2018, 2018
Short summary
Short summary
Ice shelves are thick sheets of ice floating on the ocean off the coasts of Antarctica and Greenland. They help regulate the flow of ice off the continent. Ice shelves undergo a natural cycle of seaward flow, fracture, iceberg production and regrowth. The Brunt Ice Shelf recently developed two large cracks. We used ground-penetrating radar to find out how the internal structure of the ice might influence the present crack development and the future stability of the ice shelf.
Ronja Reese, Ricarda Winkelmann, and G. Hilmar Gudmundsson
The Cryosphere, 12, 3229–3242, https://doi.org/10.5194/tc-12-3229-2018, https://doi.org/10.5194/tc-12-3229-2018, 2018
Short summary
Short summary
Accurately representing grounding-line flux is essential for modelling the evolution of the Antarctic Ice Sheet. Currently, in some large-scale ice-flow modelling studies a condition on ice flux across grounding lines is imposed using an analytically motivated parameterisation. Here we test this expression for Antarctic grounding lines and find that it provides inaccurate and partly unphysical estimates of ice flux for the highly buttressed ice streams.
Emily A. Hill, J. Rachel Carr, Chris R. Stokes, and G. Hilmar Gudmundsson
The Cryosphere, 12, 3243–3263, https://doi.org/10.5194/tc-12-3243-2018, https://doi.org/10.5194/tc-12-3243-2018, 2018
Short summary
Short summary
The dynamic behaviour (i.e. acceleration and retreat) of outlet glaciers in northern Greenland remains understudied. Here, we provide a new long-term (68-year) record of terminus change. Overall, recent retreat rates (1995–2015) are higher than the last 47 years. Despite region-wide retreat, we found disparities in dynamic behaviour depending on terminus type; grounded glaciers accelerated and thinned following retreat, while glaciers with floating ice tongues were insensitive to recent retreat.
Veronika Emetc, Paul Tregoning, Mathieu Morlighem, Chris Borstad, and Malcolm Sambridge
The Cryosphere, 12, 3187–3213, https://doi.org/10.5194/tc-12-3187-2018, https://doi.org/10.5194/tc-12-3187-2018, 2018
Short summary
Short summary
The paper includes a model that can be used to predict zones of fracture formation in both floating and grounded ice in Antarctica. We used observations and a statistics-based model to predict fractures in most ice shelves in Antarctica as an alternative to the damage-based approach. We can predict the location of observed fractures with an average success rate of 84% for grounded ice and 61% for floating ice and mean overestimation error of 26% and 20%, respectively.
Hélène Seroussi and Mathieu Morlighem
The Cryosphere, 12, 3085–3096, https://doi.org/10.5194/tc-12-3085-2018, https://doi.org/10.5194/tc-12-3085-2018, 2018
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2637–2652, https://doi.org/10.5194/tc-12-2637-2018, https://doi.org/10.5194/tc-12-2637-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula. The model was also used to simulate the distribution of temperature within the ice, which governs the rate at which the ice can deform. This is especially important for glaciers like the Fleming Glacier, which has both regions of rapid deformation and regions of rapid sliding at the bed.
Chen Zhao, Rupert M. Gladstone, Roland C. Warner, Matt A. King, Thomas Zwinger, and Mathieu Morlighem
The Cryosphere, 12, 2653–2666, https://doi.org/10.5194/tc-12-2653-2018, https://doi.org/10.5194/tc-12-2653-2018, 2018
Short summary
Short summary
A combination of computer modelling and observational data were used to infer the resistance to ice flow at the bed of the Fleming Glacier on the Antarctic Peninsula in both 2008 and 2015. The comparison suggests the grounding line retreated by ~ 9 km from 2008 to 2015. The retreat may be enhanced by a positive feedback between friction, melting and sliding at the glacier bed.
Aleah Sommers, Harihar Rajaram, and Mathieu Morlighem
Geosci. Model Dev., 11, 2955–2974, https://doi.org/10.5194/gmd-11-2955-2018, https://doi.org/10.5194/gmd-11-2955-2018, 2018
Short summary
Short summary
Meltwater drainage beneath glaciers and ice sheets influences how fast they move and is complicated and constantly changing. Most models distinguish between
fastand
slowdrainage with different equations for each system. The SHAKTI model allows for the ice–water drainage arrangement to transition naturally between different types of flow. This model can be used to understand how drainage affects glacier speeds and the associated ice loss to further inform predictions of sea level rise.
Sebastian H. R. Rosier and G. Hilmar Gudmundsson
The Cryosphere, 12, 1699–1713, https://doi.org/10.5194/tc-12-1699-2018, https://doi.org/10.5194/tc-12-1699-2018, 2018
Short summary
Short summary
Ocean tides cause strong modulation of horizontal ice shelf flow, most notably at a fortnightly frequency that is absent in the vertical tidal forcing. We propose that tidal bending in the margins of the ice shelf produces sufficiently large stresses that the effective viscosity of ice in these regions is reduced during high and low tide. This effect can explain many features of the observed behaviour and implies that ice shelves in areas with strong tides move faster than they otherwise would.
Joshua K. Cuzzone, Mathieu Morlighem, Eric Larour, Nicole Schlegel, and Helene Seroussi
Geosci. Model Dev., 11, 1683–1694, https://doi.org/10.5194/gmd-11-1683-2018, https://doi.org/10.5194/gmd-11-1683-2018, 2018
Short summary
Short summary
This paper details the implementation of higher-order vertical finite elements in the Ice Sheet System Model (ISSM). When using higher-order vertical finite elements, fewer vertical layers are needed to accurately capture the thermal structure in an ice sheet versus a conventional linear vertical interpolation, therefore greatly improving model runtime speeds, particularly in higher-order stress balance ice sheet models. The implications for paleoclimate ice sheet simulations are discussed.
Konstanze Haubner, Jason E. Box, Nicole J. Schlegel, Eric Y. Larour, Mathieu Morlighem, Anne M. Solgaard, Kristian K. Kjeldsen, Signe H. Larsen, Eric Rignot, Todd K. Dupont, and Kurt H. Kjær
The Cryosphere, 12, 1511–1522, https://doi.org/10.5194/tc-12-1511-2018, https://doi.org/10.5194/tc-12-1511-2018, 2018
Short summary
Short summary
We investigate the effect of neglecting calving on Upernavik Isstrøm, West Greenland, between 1849 and 2012.
Our simulation is forced with observed terminus positions in discrete time steps and is responsive to the prescribed ice front changes.
Simulated frontal retreat is needed to obtain a realistic ice surface elevation and velocity evolution of Upernavik.
Using the prescribed terminus position change we gain insight to mass loss partitioning during different time periods.
Heiko Goelzer, Sophie Nowicki, Tamsin Edwards, Matthew Beckley, Ayako Abe-Ouchi, Andy Aschwanden, Reinhard Calov, Olivier Gagliardini, Fabien Gillet-Chaulet, Nicholas R. Golledge, Jonathan Gregory, Ralf Greve, Angelika Humbert, Philippe Huybrechts, Joseph H. Kennedy, Eric Larour, William H. Lipscomb, Sébastien Le clec'h, Victoria Lee, Mathieu Morlighem, Frank Pattyn, Antony J. Payne, Christian Rodehacke, Martin Rückamp, Fuyuki Saito, Nicole Schlegel, Helene Seroussi, Andrew Shepherd, Sainan Sun, Roderik van de Wal, and Florian A. Ziemen
The Cryosphere, 12, 1433–1460, https://doi.org/10.5194/tc-12-1433-2018, https://doi.org/10.5194/tc-12-1433-2018, 2018
Short summary
Short summary
We have compared a wide spectrum of different initialisation techniques used in the ice sheet modelling community to define the modelled present-day Greenland ice sheet state as a starting point for physically based future-sea-level-change projections. Compared to earlier community-wide comparisons, we find better agreement across different models, which implies overall improvement of our understanding of what is needed to produce such initial states.
Hafeez Jeofry, Neil Ross, Hugh F. J. Corr, Jilu Li, Mathieu Morlighem, Prasad Gogineni, and Martin J. Siegert
Earth Syst. Sci. Data, 10, 711–725, https://doi.org/10.5194/essd-10-711-2018, https://doi.org/10.5194/essd-10-711-2018, 2018
Short summary
Short summary
Accurately characterizing the complexities of the ice-sheet dynamic specifically close to the grounding line across the Weddell Sea (WS) sector in the ice-sheet models provides challenges to the scientific community. Our main objective is to comprehend these complexities, adding accuracy to the projection of future ice-sheet dynamics. Therefore, we have developed a new bed elevation digital elevation model across the WS sector, which will be of value to ice-sheet modelling experiments.
Felicity S. Graham, Mathieu Morlighem, Roland C. Warner, and Adam Treverrow
The Cryosphere, 12, 1047–1067, https://doi.org/10.5194/tc-12-1047-2018, https://doi.org/10.5194/tc-12-1047-2018, 2018
Short summary
Short summary
Ice sheet flow is anisotropic, depending on the nature of the stress applied. However, most large-scale ice sheet models rely on the Glen flow relation, which ignores anisotropic effects. We implement a flow relation (ESTAR) for anisotropic ice in a large-scale ice sheet model. In ice shelf simulations, the Glen flow relation overestimates velocities by up to 17 % compared with ESTAR. Our results have implications for ice sheet model simulations of paleo-ice extent and sea level rise prediction.
Jan De Rydt, G. Hilmar Gudmundsson, Thomas Nagler, Jan Wuite, and Edward C. King
The Cryosphere, 12, 505–520, https://doi.org/10.5194/tc-12-505-2018, https://doi.org/10.5194/tc-12-505-2018, 2018
Short summary
Short summary
We provide an unprecedented view into the dynamics of two active rifts in the Brunt Ice Shelf through a unique set of field observations, novel satellite data products, and a state-of-the-art ice flow model. We describe the evolution of fracture width and length in great detail, pushing the boundaries of both spatial and temporal coverage, and provide a deeper insight into the process of iceberg formation, which exerts an important control over the mass balance of the Antarctic Ice Sheet.
Werner M. J. Lazeroms, Adrian Jenkins, G. Hilmar Gudmundsson, and Roderik S. W. van de Wal
The Cryosphere, 12, 49–70, https://doi.org/10.5194/tc-12-49-2018, https://doi.org/10.5194/tc-12-49-2018, 2018
Short summary
Short summary
Basal melting of ice shelves is a major factor in the decline of the Antarctic Ice Sheet, which can contribute significantly to sea-level rise. Here, we investigate a new basal melt model based on the dynamics of meltwater plumes. For the first time, this model is applied to all Antarctic ice shelves. The model results in a realistic melt-rate pattern given suitable data for the topography and ocean temperature, making it a promising tool for future simulations of the Antarctic Ice Sheet.
Eric Larour, Daniel Cheng, Gilberto Perez, Justin Quinn, Mathieu Morlighem, Bao Duong, Lan Nguyen, Kit Petrie, Silva Harounian, Daria Halkides, and Wayne Hayes
Geosci. Model Dev., 10, 4393–4403, https://doi.org/10.5194/gmd-10-4393-2017, https://doi.org/10.5194/gmd-10-4393-2017, 2017
Short summary
Short summary
This work presents a new way of carrying out simulations using the C++ based Ice Sheet System Model (ISSM) within a web page. This allows for a new generation of websites that can rely on the entire code of a climate model, without compromising or simplifying the physics implemented in such a model. We believe this approach will enable better education/outreach websites as well as improve access to complex climate models without compromising their integrity.
Sebastian H. R. Rosier, G. Hilmar Gudmundsson, Matt A. King, Keith W. Nicholls, Keith Makinson, and Hugh F. J. Corr
Earth Syst. Sci. Data, 9, 849–860, https://doi.org/10.5194/essd-9-849-2017, https://doi.org/10.5194/essd-9-849-2017, 2017
Short summary
Short summary
Tides can affect the flow of ice at hourly to yearly timescales. In some cases the ice responds at a different frequency than is found in the tidal forcing; for example, on Rutford Ice Stream the strongest response is at a fortnightly period. A new compilation of GPS data across the Ronne Ice Shelf and adjoining ice streams shows that this fortnightly modulation in ice flow is found across the entire region. Measurements of this kind can provide insights into ice rheology and basal processes.
Hongju Yu, Eric Rignot, Mathieu Morlighem, and Helene Seroussi
The Cryosphere, 11, 1283–1296, https://doi.org/10.5194/tc-11-1283-2017, https://doi.org/10.5194/tc-11-1283-2017, 2017
Short summary
Short summary
We combine 2-D ice flow model with linear elastic fracture mechanics (LEFM) to model the calving behavior of Thwaites Glacier, West Antarctica. We find the combination of full-Stokes (FS) model and LEFM produces crevasses that are consistent with observations. We also find that calving is enhanced with pre-existing surface crevasses, shorter ice shelves or undercut at the ice shelf front. We conclude that the FS/LEFM combination is capable of constraining crevasse formation and iceberg calving.
Daniel Farinotti, Douglas J. Brinkerhoff, Garry K. C. Clarke, Johannes J. Fürst, Holger Frey, Prateek Gantayat, Fabien Gillet-Chaulet, Claire Girard, Matthias Huss, Paul W. Leclercq, Andreas Linsbauer, Horst Machguth, Carlos Martin, Fabien Maussion, Mathieu Morlighem, Cyrille Mosbeux, Ankur Pandit, Andrea Portmann, Antoine Rabatel, RAAJ Ramsankaran, Thomas J. Reerink, Olivier Sanchez, Peter A. Stentoft, Sangita Singh Kumari, Ward J. J. van Pelt, Brian Anderson, Toby Benham, Daniel Binder, Julian A. Dowdeswell, Andrea Fischer, Kay Helfricht, Stanislav Kutuzov, Ivan Lavrentiev, Robert McNabb, G. Hilmar Gudmundsson, Huilin Li, and Liss M. Andreassen
The Cryosphere, 11, 949–970, https://doi.org/10.5194/tc-11-949-2017, https://doi.org/10.5194/tc-11-949-2017, 2017
Short summary
Short summary
ITMIX – the Ice Thickness Models Intercomparison eXperiment – was the first coordinated performance assessment for models inferring glacier ice thickness from surface characteristics. Considering 17 different models and 21 different test cases, we show that although solutions of individual models can differ considerably, an ensemble average can yield uncertainties in the order of 10 ± 24 % the mean ice thickness. Ways forward for improving such estimates are sketched.
Henning Åkesson, Kerim H. Nisancioglu, Rianne H. Giesen, and Mathieu Morlighem
The Cryosphere, 11, 281–302, https://doi.org/10.5194/tc-11-281-2017, https://doi.org/10.5194/tc-11-281-2017, 2017
Short summary
Short summary
We present simulations of the history of Hardangerjøkulen ice cap in southern Norway using a dynamical ice sheet model. From mid-Holocene ice-free conditions 4000 years ago, Hardangerjøkulen grows nonlinearly in response to a linear climate forcing, reaching maximum extent during the Little Ice Age (~ 1750 AD). The ice cap exhibits spatially asymmetric growth and retreat and is highly sensitive to climate change. Our results call for reassessment of glacier reconstructions from proxy records.
Feras Habbal, Eric Larour, Mathieu Morlighem, Helene Seroussi, Christopher P. Borstad, and Eric Rignot
Geosci. Model Dev., 10, 155–168, https://doi.org/10.5194/gmd-10-155-2017, https://doi.org/10.5194/gmd-10-155-2017, 2017
Short summary
Short summary
This work presents the results from testing a suite of numerical solvers on a standard ice sheet benchmark test. We note the relevance of this test to practical simulations and identify the fastest solvers for the transient simulation. The highlighted solvers show significant speed-ups in relation to the default solver (~1.5–100 times faster) and enable a new capability for solving massive, high-resolution models that are critical for improving projections of ice sheets and sea-level change.
Eric Larour, Jean Utke, Anton Bovin, Mathieu Morlighem, and Gilberto Perez
Geosci. Model Dev., 9, 3907–3918, https://doi.org/10.5194/gmd-9-3907-2016, https://doi.org/10.5194/gmd-9-3907-2016, 2016
Short summary
Short summary
We present an approach to derive the adjoint state of the C++ coded Ice Sheet System Model. The approach enables data assimilation of observations to improve projections of polar ice sheet mass balance and contribution to sea-level rise. It is applicable to other Earth science frameworks relying on C++ and parallel computing, is non-intrusive, and enables computation of transient adjoints for any type of physics, hence providing insights into the sensitivities of any model to its inputs.
Janin Schaffer, Ralph Timmermann, Jan Erik Arndt, Steen Savstrup Kristensen, Christoph Mayer, Mathieu Morlighem, and Daniel Steinhage
Earth Syst. Sci. Data, 8, 543–557, https://doi.org/10.5194/essd-8-543-2016, https://doi.org/10.5194/essd-8-543-2016, 2016
Short summary
Short summary
The RTopo-2 data set provides consistent maps of global ocean bathymetry and ice surface topographies for Greenland and Antarctica at 30 arcsec grid spacing. We corrected data from earlier products in the areas of Petermann, Hagen Bræ, and Helheim glaciers, incorporated original data for the floating ice tongue of Nioghalvfjerdsfjorden Glacier, and applied corrections for the geometry of Getz, Abbot, and Fimbul ice shelf cavities. The data set is available from the PANGAEA database.
Xylar S. Asay-Davis, Stephen L. Cornford, Gaël Durand, Benjamin K. Galton-Fenzi, Rupert M. Gladstone, G. Hilmar Gudmundsson, Tore Hattermann, David M. Holland, Denise Holland, Paul R. Holland, Daniel F. Martin, Pierre Mathiot, Frank Pattyn, and Hélène Seroussi
Geosci. Model Dev., 9, 2471–2497, https://doi.org/10.5194/gmd-9-2471-2016, https://doi.org/10.5194/gmd-9-2471-2016, 2016
Short summary
Short summary
Coupled ice sheet–ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of ice sheets and glaciers, including assessing their contributions to sea level change. Here we describe the idealized experiments that make up three interrelated Model Intercomparison Projects (MIPs) for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities.
Hongju Yu, Eric Rignot, Mathieu Morlighem, and Helene Seroussi
The Cryosphere Discuss., https://doi.org/10.5194/tc-2016-101, https://doi.org/10.5194/tc-2016-101, 2016
Revised manuscript not accepted
Short summary
Short summary
We performed a 2D Full-Stokes (FS) modeling study of grounding line dynamics and calving of Thwaites Glacier, West Antarctica. We compare FS with simplified models on grounding line migration and we combine FS with Linear Elastic Fracture Mechanics to simulate crevasse propagation. We find that only FS is able to provide reliable grounding line migration and to explain observed crevasse. We conclude that it may be essential to employ FS in the grounding line region for 2D simulations.
Pierre Rampal, Sylvain Bouillon, Einar Ólason, and Mathieu Morlighem
The Cryosphere, 10, 1055–1073, https://doi.org/10.5194/tc-10-1055-2016, https://doi.org/10.5194/tc-10-1055-2016, 2016
Short summary
Short summary
The Arctic sea ice cover has changed drastically over the last decades and undergone a shift in its dynamical regime, as seen by the increase of extreme fracturing events and the acceleration of sea ice drift. In this paper we present a new sea ice model, neXtSIM, that is capable of simulating both sea ice drift and deformation as observed from satellites, with similar spatial and temporal scaling properties. At the same time, the model reproduces sea ice area, extent, and volume correctly.
David H. Jones, Carl Robinson, and G. Hilmar Gudmundsson
Geosci. Instrum. Method. Data Syst., 5, 65–73, https://doi.org/10.5194/gi-5-65-2016, https://doi.org/10.5194/gi-5-65-2016, 2016
Short summary
Short summary
Long-term records from high-precision GPS receivers are essential for studying geophysical movement. Existing, commercially available, precision GPS receivers are not intended for long-term, autonomous deployment. We have designed a GPS receiver that is better suited for this application. In this paper, we discuss the receiver design and compare its performance with that of some of the commercially available receivers.
Johannes H. Bondzio, Hélène Seroussi, Mathieu Morlighem, Thomas Kleiner, Martin Rückamp, Angelika Humbert, and Eric Y. Larour
The Cryosphere, 10, 497–510, https://doi.org/10.5194/tc-10-497-2016, https://doi.org/10.5194/tc-10-497-2016, 2016
Short summary
Short summary
We implemented a level-set method in the ice sheet system model. This method allows us to dynamically evolve a calving front subject to user-defined calving rates. We apply the method to Jakobshavn Isbræ, West Greenland, and study its response to calving rate perturbations. We find its behaviour strongly dependent on the calving rate, which was to be expected. Both reduced basal drag and rheological shear margin weakening sustain the acceleration of this dynamic outlet glacier.
Laura A. Stevens, Fiamma Straneo, Sarah B. Das, Albert J. Plueddemann, Amy L. Kukulya, and Mathieu Morlighem
The Cryosphere, 10, 417–432, https://doi.org/10.5194/tc-10-417-2016, https://doi.org/10.5194/tc-10-417-2016, 2016
Short summary
Short summary
Here we pair detailed hydrographic measurements collected with an autonomous underwater vehicle as close as 150 m from the ice–ocean interface of the Saqqarliup sermia–Sarqardleq Fjord system, West Greenland, with modeled and observed subglacial discharge locations and magnitudes. We find evidence of two main types of subsurface glacially modified water localized in space that are consistent with runoff discharged at two locations along the grounding line.
K. Le Morzadec, L. Tarasov, M. Morlighem, and H. Seroussi
Geosci. Model Dev., 8, 3199–3213, https://doi.org/10.5194/gmd-8-3199-2015, https://doi.org/10.5194/gmd-8-3199-2015, 2015
Short summary
Short summary
A long-term challenge for any model of complex large-scale processes
is accounting for the impact of unresolved sub-grid (SG) processes.
We quantify the impact of SG mass-balance and ice fluxes on glacial
cycle ensemble results for North America. We find no easy solutions to
accurately capture these impacts. We show that SG process
representation and associated parametric uncertainties can have
significant impact on coarse resolution model results for glacial
cycle ice sheet evolution.
S. H. R. Rosier, G. H. Gudmundsson, and J. A. M. Green
The Cryosphere, 9, 1649–1661, https://doi.org/10.5194/tc-9-1649-2015, https://doi.org/10.5194/tc-9-1649-2015, 2015
Short summary
Short summary
We use a full-Stokes model to investigate the long period modulation of Rutford Ice Stream flow by the ocean tide. We find that using a nonlinear sliding law cannot fully explain the measurements and an additional mechanism, whereby tidally induced subglacial pressure variations are transmitted upstream from the grounding line, is also required to match the large amplitude and decay length scale of the observations.
D. H. Jones and G. H. Gudmundsson
Nat. Hazards Earth Syst. Sci., 15, 1243–1250, https://doi.org/10.5194/nhess-15-1243-2015, https://doi.org/10.5194/nhess-15-1243-2015, 2015
Short summary
Short summary
Icebergs are a natural hazard to maritime operations in polar regions. Iceberg populations are increasing, as is the demand for access to both Arctic and Antarctic seas. Soon the ability to reliably track icebergs may become a necessity for continued operational safety. In this paper we describe the design of a tracking sensor that can be deployed from an aircraft during surveys of Antarctic icebergs, and detail the results of its first deployment operation on iceberg B-31.
J. Wuite, H. Rott, M. Hetzenecker, D. Floricioiu, J. De Rydt, G. H. Gudmundsson, T. Nagler, and M. Kern
The Cryosphere, 9, 957–969, https://doi.org/10.5194/tc-9-957-2015, https://doi.org/10.5194/tc-9-957-2015, 2015
Short summary
Short summary
We present new analysis of satellite data showing the variability of glacier velocities in the Larsen B area, Antarctic Peninsula, back to 1995. Velocity data and estimates of ice thickness are used to derive ice discharge at different epochs. Velocities of the glaciers remain to date well above the velocities of the pre-collapse period. The response of individual glaciers differs, and velocities show significant temporal fluctuations, implying major variations in ice discharge and mass balance.
C. Martín, R. Mulvaney, G. H. Gudmundsson, and H. Corr
Clim. Past, 11, 547–557, https://doi.org/10.5194/cp-11-547-2015, https://doi.org/10.5194/cp-11-547-2015, 2015
E. Larour, J. Utke, B. Csatho, A. Schenk, H. Seroussi, M. Morlighem, E. Rignot, N. Schlegel, and A. Khazendar
The Cryosphere, 8, 2335–2351, https://doi.org/10.5194/tc-8-2335-2014, https://doi.org/10.5194/tc-8-2335-2014, 2014
Short summary
Short summary
We present a temporal inversion of surface mass balance and basal friction for the Northeast Greenland Ice Sheet between 2003 and 2009, using the altimetry record from ICESat. The inversion relies on automatic differentiation of ISSM and demonstrates the feasibility of assimilating altimetry records into reconstructions of the Greenland Ice Sheet. The boundary conditions provide a snapshot of the state of the ice for this period and can be used for further process studies.
H. Seroussi, M. Morlighem, E. Larour, E. Rignot, and A. Khazendar
The Cryosphere, 8, 2075–2087, https://doi.org/10.5194/tc-8-2075-2014, https://doi.org/10.5194/tc-8-2075-2014, 2014
J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick
The Cryosphere, 8, 2031–2045, https://doi.org/10.5194/tc-8-2031-2014, https://doi.org/10.5194/tc-8-2031-2014, 2014
S. H. R. Rosier, G. H. Gudmundsson, and J. A. M. Green
The Cryosphere, 8, 1763–1775, https://doi.org/10.5194/tc-8-1763-2014, https://doi.org/10.5194/tc-8-1763-2014, 2014
H. Seroussi, M. Morlighem, E. Rignot, J. Mouginot, E. Larour, M. Schodlok, and A. Khazendar
The Cryosphere, 8, 1699–1710, https://doi.org/10.5194/tc-8-1699-2014, https://doi.org/10.5194/tc-8-1699-2014, 2014
S. Adhikari, E. R. Ivins, E. Larour, H. Seroussi, M. Morlighem, and S. Nowicki
Solid Earth, 5, 569–584, https://doi.org/10.5194/se-5-569-2014, https://doi.org/10.5194/se-5-569-2014, 2014
R. Anderson, D. H. Jones, and G. H. Gudmundsson
Nat. Hazards Earth Syst. Sci., 14, 917–927, https://doi.org/10.5194/nhess-14-917-2014, https://doi.org/10.5194/nhess-14-917-2014, 2014
G. H. Gudmundsson
The Cryosphere, 7, 647–655, https://doi.org/10.5194/tc-7-647-2013, https://doi.org/10.5194/tc-7-647-2013, 2013
J. De Rydt, G. H. Gudmundsson, H. F. J. Corr, and P. Christoffersen
The Cryosphere, 7, 407–417, https://doi.org/10.5194/tc-7-407-2013, https://doi.org/10.5194/tc-7-407-2013, 2013
G. H. Gudmundsson, J. Krug, G. Durand, L. Favier, and O. Gagliardini
The Cryosphere, 6, 1497–1505, https://doi.org/10.5194/tc-6-1497-2012, https://doi.org/10.5194/tc-6-1497-2012, 2012
Related subject area
Numerical methods
Hydro-geomorphological modelling of leaky wooden dam efficacy from reach to catchment scale with CAESAR-Lisflood 1.9j
Enhancing single precision with quasi-double precision: achieving double-precision accuracy in the Model for Prediction Across Scales – Atmosphere (MPAS-A) version 8.2.1
Advances in land surface forecasting: a comparison of LSTM, gradient boosting, and feed-forward neural networks as prognostic state emulators in a case study with ecLand
Subgrid corrections for the linear inertial equations of a compound flood model – a case study using SFINCS 2.1.1 Dollerup release
Introducing Iterative Model Calibration (IMC) v1.0: a generalizable framework for numerical model calibration with a CAESAR-Lisflood case study
Development of a high-order global dynamical core using the discontinuous Galerkin method for an atmospheric large-eddy simulation (LES) and proposal of test cases: SCALE-DG v0.8.0
A joint reconstruction and model selection approach for large-scale linear inverse modeling (msHyBR v2)
Assimilation of snow water equivalent from AMSR2 and IMS satellite data utilizing the local ensemble transform Kalman filter
The Paleochrono-1.1 probabilistic model to derive a common age model for several paleoclimatic sites using absolute and relative dating constraints
asQ: parallel-in-time finite element simulations using ParaDiag for geoscientific models and beyond
Explicit stochastic advection algorithms for the regional-scale particle-resolved atmospheric aerosol model WRF-PartMC (v1.0)
Potential based Thermodynamics with Consistent Conservative Cascade Transport for Implicit Large Eddy Simulation: PTerodaC3TILES version 1.0
The Measurement Error Proxy System Model: MEPSM v0.2
Modelling chemical advection during magma ascent
Consistent point data assimilation in Firedrake and Icepack
A computationally efficient parameterization of aerosol, cloud and precipitation pH for application at global and regional scale (EQSAM4Clim-v12)
Assessing the benefits of approximately exact step sizes for Picard and Newton solver in simulating ice flow (FEniCS-full-Stokes v.1.3.2)
Assessing effects of climate and technology uncertainties in large natural resource allocation problems
VISIR-2: ship weather routing in Python
Incremental analysis update (IAU) in the Model for Prediction Across Scales coupled with the Joint Effort for Data assimilation Integration (MPAS–JEDI 2.0.0)
Decision-making strategies implemented in SolFinder 1.0 to identify eco-efficient aircraft trajectories: application study in AirTraf 3.0
CLAQC v1.0 – Country Level Air Quality Calculator. An Empirical Modeling Approach
Developing meshing workflows in Gmsh v4.11 for the geologic uncertainty assessment of high-temperature aquifer thermal energy storage
Development and preliminary validation of a land surface image assimilation system based on the Common Land Model
NorSand4AI: a comprehensive triaxial test simulation database for NorSand constitutive model materials
ParticleDA.jl v.1.0: a distributed particle-filtering data assimilation package
HETerogeneous vectorized or Parallel (HETPv1.0): an updated inorganic heterogeneous chemistry solver for the metastable-state NH4+–Na+–Ca2+–K+–Mg2+–SO42−–NO3−–Cl−–H2O system based on ISORROPIA II
Three-dimensional geological modelling of igneous intrusions in LoopStructural v1.5.10
Estimating volcanic ash emissions using retrieved satellite ash columns and inverse ash transport modeling using VolcanicAshInversion v1.2.1, within the operational eEMEP (emergency European Monitoring and Evaluation Programme) volcanic plume forecasting system (version rv4_17)
Accounting for uncertainties in forecasting tropical-cyclone-induced compound flooding
An automatic mesh generator for coupled 1D–2D hydrodynamic models
Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme
Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 2: A semi-discrete error analysis framework for assessing coupling schemes
jsmetrics v0.2.0: a Python package for metrics and algorithms used to identify or characterise atmospheric jet streams
P3D-BRNS v1.0.0: a three-dimensional, multiphase, multicomponent, pore-scale reactive transport modelling package for simulating biogeochemical processes in subsurface environments
MinVoellmy v1: a lightweight model for simulating rapid mass movements based on a modified Voellmy rheology
Scalable Feature Extraction and Tracking (SCAFET): a general framework for feature extraction from large climate data sets
Sweep interpolation: a cost-effective semi-Lagrangian scheme in the Global Environmental Multiscale model
CHONK 1.0: landscape evolution framework: cellular automata meets graph theory
Perspectives of physics-based machine learning strategies for geoscientific applications governed by partial differential equations
Calibration of absorbing boundary layers for geoacoustic wave modeling in pseudo-spectral time-domain methods
GeoINR 1.0: an implicit neural network approach to three-dimensional geological modelling
Parameterization and tuning of the Bay of Biscay Atlantis model v1
A comparison of Eulerian and Lagrangian methods for vertical particle transport in the water column
AutoQS v1: automatic parametrization of QuickSampling based on training images analysis
Implementation and application of ensemble optimal interpolation on an operational chemistry weather model for improving PM2.5 and visibility predictions
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling
GStatSim V1.0: a Python package for geostatistical interpolation and conditional simulation
Leveraging Google's Tensor Processing Units for tsunami-risk mitigation planning in the Pacific Northwest and beyond
An improved subgrid channel model with upwind-form artificial diffusion for river hydrodynamics and floodplain inundation simulation
Joshua M. Wolstenholme, Christopher J. Skinner, David Milan, Robert E. Thomas, and Daniel R. Parsons
Geosci. Model Dev., 18, 1395–1411, https://doi.org/10.5194/gmd-18-1395-2025, https://doi.org/10.5194/gmd-18-1395-2025, 2025
Short summary
Short summary
Leaky wooden dams are a type of natural flood management intervention that aims to reduce flood risk downstream by temporarily holding back water during a storm event and releasing it afterwards. These structures alter the river hydrology, and therefore the geomorphology, yet often this is excluded from numerical models. Here we show that by not simulating geomorphology, we are currently underestimating the efficacy of these structures to reduce the flood peak and store water.
Jiayi Lai, Lanning Wang, Qizhong Wu, Yizhou Yang, and Fang Wang
Geosci. Model Dev., 18, 1089–1102, https://doi.org/10.5194/gmd-18-1089-2025, https://doi.org/10.5194/gmd-18-1089-2025, 2025
Short summary
Short summary
High-performance computing limitations often hinder numerical model development. Traditional models use double precision for accuracy, which is computationally expensive. Lower precision reduces costs but can introduce errors. The quasi-double-precision (QDP) algorithm helps mitigate these errors. This study applies the QDP algorithm to the Model for Prediction Across Scales – Atmosphere, showing reduced errors and computational time, making it an efficient solution for large-scale simulations.
Marieke Wesselkamp, Matthew Chantry, Ewan Pinnington, Margarita Choulga, Souhail Boussetta, Maria Kalweit, Joschka Bödecker, Carsten F. Dormann, Florian Pappenberger, and Gianpaolo Balsamo
Geosci. Model Dev., 18, 921–937, https://doi.org/10.5194/gmd-18-921-2025, https://doi.org/10.5194/gmd-18-921-2025, 2025
Short summary
Short summary
We compared spatiotemporal forecasts of three machine learning models that learned water and energy
states on the land surface from a physical model scheme. The forecasting models were developed with reanalysis data and simulations on a European scale and transferred to the globe. We found that all approaches deliver highly accurate approximations of the physical dynamic at long time horizons, implying their usefulness to advance land surface forecasting with synthetic data.
states on the land surface from a physical model scheme. The forecasting models were developed with reanalysis data and simulations on a European scale and transferred to the globe. We found that all approaches deliver highly accurate approximations of the physical dynamic at long time horizons, implying their usefulness to advance land surface forecasting with synthetic data.
Maarten van Ormondt, Tim Leijnse, Roel de Goede, Kees Nederhoff, and Ap van Dongeren
Geosci. Model Dev., 18, 843–861, https://doi.org/10.5194/gmd-18-843-2025, https://doi.org/10.5194/gmd-18-843-2025, 2025
Short summary
Short summary
Accurate flood risk assessments are crucial for storm protection. To achieve efficiency, computational costs must be minimized. This paper introduces a novel subgrid approach for linear inertial equations (LIEs) with bed level and friction variations, implemented in the Super-Fast INundation of CoastS (SFINCS) model. Pre-processed lookup tables enhance simulation precision with lower costs. Validations show significant accuracy improvement even at coarser resolutions.
Chayan Banerjee, Kien Nguyen, Clinton Fookes, Gregory Hancock, and Thomas Coulthard
Geosci. Model Dev., 18, 803–818, https://doi.org/10.5194/gmd-18-803-2025, https://doi.org/10.5194/gmd-18-803-2025, 2025
Short summary
Short summary
In geosciences, the reliance on numerical models necessitates the precise calibration of their parameters to effectively translate information from observed to unobserved settings. We introduce a generalizable framework for calibrating numerical models, with a case study of the geomorphological model CAESAR-Lisflood. This approach efficiently identifies the optimal set of parameters for a given numerical model, enabling retrospective and prospective analyses at various temporal resolutions.
Yuta Kawai and Hirofumi Tomita
Geosci. Model Dev., 18, 725–762, https://doi.org/10.5194/gmd-18-725-2025, https://doi.org/10.5194/gmd-18-725-2025, 2025
Short summary
Short summary
When considering future high-resolution atmospheric simulations with O(10–100 m) grid spacing, such as global large-eddy simulations, numerical errors with the conventional low-order dynamical cores can contaminate the effects of turbulent parameterization. To achieve high-order discretization, we developed a global atmospheric dynamical core using the discontinuous Galerkin method (DGM). By conducting several validation tests, we discussed the impact of a high-order DGM on atmospheric flows.
Malena Sabaté Landman, Julianne Chung, Jiahua Jiang, Scot M. Miller, and Arvind K. Saibaba
Geosci. Model Dev., 17, 8853–8872, https://doi.org/10.5194/gmd-17-8853-2024, https://doi.org/10.5194/gmd-17-8853-2024, 2024
Short summary
Short summary
Making an informed decision about what prior information to incorporate or discard in an inverse model is important yet very challenging, as it is often not straightforward to distinguish between informative and non-informative variables. In this study, we develop a new approach for incorporating prior information in an inverse model using predictor variables, while simultaneously selecting the relevant predictor variables for the estimation of the unknown quantity of interest.
Joonlee Lee, Myong-In Lee, Sunlae Tak, Eunkyo Seo, and Yong-Keun Lee
Geosci. Model Dev., 17, 8799–8816, https://doi.org/10.5194/gmd-17-8799-2024, https://doi.org/10.5194/gmd-17-8799-2024, 2024
Short summary
Short summary
We developed an advanced snow water equivalent (SWE) data assimilation framework using satellite data based on a land surface model. The results of this study highlight the beneficial impact of data assimilation by effectively combining land surface model and satellite-derived data according to their relative uncertainty, thereby controlling not only transitional regions but also the regions with heavy snow accumulation that are difficult to detect by satellite.
Frédéric Parrenin, Marie Bouchet, Christo Buizert, Emilie Capron, Ellen Corrick, Russell Drysdale, Kenji Kawamura, Amaëlle Landais, Robert Mulvaney, Ikumi Oyabu, and Sune Olander Rasmussen
Geosci. Model Dev., 17, 8735–8750, https://doi.org/10.5194/gmd-17-8735-2024, https://doi.org/10.5194/gmd-17-8735-2024, 2024
Short summary
Short summary
The Paleochrono-1.1 probabilistic dating model allows users to derive a common and optimized chronology for several paleoclimatic sites from various archives (ice cores, speleothems, marine cores, lake cores, etc.). It combines prior sedimentation scenarios with chronological information such as dated horizons, dated intervals, stratigraphic links and (for ice cores) Δdepth observations. Paleochrono-1.1 is available under an open-source license.
Joshua Hope-Collins, Abdalaziz Hamdan, Werner Bauer, Lawrence Mitchell, and Colin Cotter
EGUsphere, https://doi.org/https://doi.org/10.48550/arXiv.2409.18792, https://doi.org/https://doi.org/10.48550/arXiv.2409.18792, 2024
Short summary
Short summary
Effectively using modern supercomputers requires massively parallel algorithms. Time-parallel algorithms calculate the system state (e.g. the atmosphere) at multiple times simultaneously and have exciting potential, but are tricky to implement and still require development. We have developed software to simplify implementing and testing the ParaDiag algorithm on supercomputers. We show that for some atmospheric problems it can enable faster or more accurate solutions than traditional techniques.
Jeffrey H. Curtis, Nicole Riemer, and Matthew West
Geosci. Model Dev., 17, 8399–8420, https://doi.org/10.5194/gmd-17-8399-2024, https://doi.org/10.5194/gmd-17-8399-2024, 2024
Short summary
Short summary
This paper introduces a numerical method for simulating particle-based aerosol transport in atmospheric models. We detail the various numerical properties of the advection order method and demonstrate its implementation in a 3D weather prediction model (WRF) for the first time. Particle-based techniques improve the accuracy of aerosol size and composition predictions, which are key for aerosol–cloud and aerosol–radiation interactions.
John Thuburn
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-153, https://doi.org/10.5194/gmd-2024-153, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
A new computational fluid dynamics code for simulating the atmospheric boundary layer and convection is presented. Moist thermodynamics is formulated via thermodynamic potentials, avoiding inconsistencies that can be introduced with conventional approaches. Numerical methods typical of weather and climate models are used, with no explicit subgrid scheme. Results highlight some advantages (e.g., large time steps) and disadvantages (e.g., weak vertical fluxes near the surface) of this approach.
Matt J. Fischer
Geosci. Model Dev., 17, 6745–6760, https://doi.org/10.5194/gmd-17-6745-2024, https://doi.org/10.5194/gmd-17-6745-2024, 2024
Short summary
Short summary
In paleoclimate research, proxy system models (PSMs) model chemical or biological systems which receive environmental input and output (e.g., geochemical signals). The environmental inputs are rarely noiseless, which causes problems when calibrating multi-input PSMs. Here a PSM is developed which includes generalized noise in both model inputs and outputs and prior information. A quasi-Bayesian method enhances the stability of the solution of the Measurement Error Proxy System Model.
Hugo Dominguez, Nicolas Riel, and Pierre Lanari
Geosci. Model Dev., 17, 6105–6122, https://doi.org/10.5194/gmd-17-6105-2024, https://doi.org/10.5194/gmd-17-6105-2024, 2024
Short summary
Short summary
Predicting the behaviour of magmatic systems is important for understanding Earth's matter and heat transport. Numerical modelling is a technique that can predict complex systems at different scales of space and time by solving equations using various techniques. This study tests four algorithms to find the best way to transport the melt composition. The "weighted essentially non-oscillatory" algorithm emerges as the best choice, minimising errors and preserving system mass well.
Reuben W. Nixon-Hill, Daniel Shapero, Colin J. Cotter, and David A. Ham
Geosci. Model Dev., 17, 5369–5386, https://doi.org/10.5194/gmd-17-5369-2024, https://doi.org/10.5194/gmd-17-5369-2024, 2024
Short summary
Short summary
Scientists often use models to study complex processes, like the movement of ice sheets, and compare them to measurements for estimating quantities that are hard to measure. We highlight an approach that ensures accurate results from point data sources (e.g. height measurements) by evaluating the numerical solution at true point locations. This method improves accuracy, aids communication between scientists, and is well-suited for integration with specialised software that automates processes.
Swen Metzger, Samuel Rémy, Jason E. Williams, Vincent Huijnen, and Johannes Flemming
Geosci. Model Dev., 17, 5009–5021, https://doi.org/10.5194/gmd-17-5009-2024, https://doi.org/10.5194/gmd-17-5009-2024, 2024
Short summary
Short summary
EQSAM4Clim has recently been revised to provide an accurate and efficient method for calculating the acidity of atmospheric particles. It is based on an analytical concept that is sufficiently fast and free of numerical noise, which makes it attractive for air quality forecasting. Version 12 allows the calculation of aerosol composition based on the gas–liquid–solid and the reduced gas–liquid partitioning with the associated water uptake for both cases, including the acidity of the aerosols.
Niko Schmidt, Angelika Humbert, and Thomas Slawig
Geosci. Model Dev., 17, 4943–4959, https://doi.org/10.5194/gmd-17-4943-2024, https://doi.org/10.5194/gmd-17-4943-2024, 2024
Short summary
Short summary
Future sea-level rise is of big significance for coastal regions. The melting and acceleration of glaciers plays a major role in sea-level change. Computer simulation of glaciers costs a lot of computational resources. In this publication, we test a new way of simulating glaciers. This approach produces the same results but has the advantage that it needs much less computation time. As simulations can be obtained with fewer computation resources, higher resolution and physics become affordable.
Jevgenijs Steinbuks, Yongyang Cai, Jonas Jaegermeyr, and Thomas W. Hertel
Geosci. Model Dev., 17, 4791–4819, https://doi.org/10.5194/gmd-17-4791-2024, https://doi.org/10.5194/gmd-17-4791-2024, 2024
Short summary
Short summary
This paper applies a cutting-edge numerical method, SCEQ, to show how uncertain climate change and technological progress affect the future utilization of the world's scarce land resources. The paper's key insight is to illustrate how much global cropland will expand when future crop yields are unknown. The study finds the range of outcomes for land use change to be smaller when using this novel method compared to existing deterministic models.
Gianandrea Mannarini, Mario Leonardo Salinas, Lorenzo Carelli, Nicola Petacco, and Josip Orović
Geosci. Model Dev., 17, 4355–4382, https://doi.org/10.5194/gmd-17-4355-2024, https://doi.org/10.5194/gmd-17-4355-2024, 2024
Short summary
Short summary
Ship weather routing has the potential to reduce CO2 emissions, but it currently lacks open and verifiable research. The Python-refactored VISIR-2 model considers currents, waves, and wind to optimise routes. The model was validated, and its computational performance is quasi-linear. For a ferry sailing in the Mediterranean Sea, VISIR-2 yields the largest percentage emission savings for upwind navigation. Given the vessel performance curve, the model is generalisable across various vessel types.
Soyoung Ha, Jonathan J. Guerrette, Ivette Hernández Baños, William C. Skamarock, and Michael G. Duda
Geosci. Model Dev., 17, 4199–4211, https://doi.org/10.5194/gmd-17-4199-2024, https://doi.org/10.5194/gmd-17-4199-2024, 2024
Short summary
Short summary
To mitigate the imbalances in the initial conditions, this study introduces our recent implementation of the incremental analysis update (IAU) in the Model for Prediction Across Scales – Atmospheric (MPAS-A) component coupled with the Joint Effort for Data assimilation Integration (JEDI) through the cycling system. A month-long cycling run demonstrates the successful implementation of the IAU capability in the MPAS–JEDI cycling system.
Federica Castino, Feijia Yin, Volker Grewe, Hiroshi Yamashita, Sigrun Matthes, Simone Dietmüller, Sabine Baumann, Manuel Soler, Abolfazl Simorgh, Maximilian Mendiguchia Meuser, Florian Linke, and Benjamin Lührs
Geosci. Model Dev., 17, 4031–4052, https://doi.org/10.5194/gmd-17-4031-2024, https://doi.org/10.5194/gmd-17-4031-2024, 2024
Short summary
Short summary
We introduce SolFinder 1.0, a decision-making tool to select trade-offs between different objective functions for optimal aircraft trajectories, including fuel use, flight time, NOx emissions, contrail distance, and climate impact. The module is included in the AirTraf 3.0 submodel and uses weather conditions simulated by the EMAC atmospheric model. This paper focuses on the ability of SolFinder to identify eco-efficient trajectories, reducing a flight's climate impact at limited cost penalties.
Stefania Renna, Francesco Granella, Lara Aleluia Reis, and Paulina Estela Schulz-Antipa
EGUsphere, https://doi.org/10.2139/ssrn.4591252, https://doi.org/10.2139/ssrn.4591252, 2024
Short summary
Short summary
The Country Level Air Quality Calculator (CLAQC) is a new fast modeling tool that predicts globally country-level monthly and annual concentrations of two major air pollutants, fine particulate matter (PM2.5) and tropospheric ozone (O3). It was designed to inform national and regional climate and pollution mitigation policies. It is easy to use and computationally efficient, allowing the simulation of a large number of policy scenarios for policy assessments and optimization frameworks.
Ali Dashti, Jens C. Grimmer, Christophe Geuzaine, Florian Bauer, and Thomas Kohl
Geosci. Model Dev., 17, 3467–3485, https://doi.org/10.5194/gmd-17-3467-2024, https://doi.org/10.5194/gmd-17-3467-2024, 2024
Short summary
Short summary
This study developed new meshing workflows to enable the automatic generation of meshes that follow geological models. The workflow allows for importing several geological models as input for Gmsh and later exporting the same number of high-quality meshes. This way, geological uncertainty is directly included in the numerical simulations. This study evaluates the impact of the geological uncertainty on thermohydraulic performance of two reservoirs for high-temperature heat storage applications.
Wangbin Shen, Zhaohui Lin, Zhengkun Qin, and Juan Li
Geosci. Model Dev., 17, 3447–3465, https://doi.org/10.5194/gmd-17-3447-2024, https://doi.org/10.5194/gmd-17-3447-2024, 2024
Short summary
Short summary
In this study, a land surface image assimilation system capable of optimizing the spatial structure of the background field is constructed by introducing the curvelet analysis method and taking the similarity of image structure as a weak constraint. The findings demonstrate that the assimilation of surface soil moisture observation images effectively and reasonably enhances the spatial structure of soil moisture analysis field.
Luan Carlos de Sena Monteiro Ozelim, Michéle Dal Toé Casagrande, and André Luís Brasil Cavalcante
Geosci. Model Dev., 17, 3175–3197, https://doi.org/10.5194/gmd-17-3175-2024, https://doi.org/10.5194/gmd-17-3175-2024, 2024
Short summary
Short summary
The paper addresses synthetic dataset challenges in nonlinear constitutive modeling of soils, providing two datasets: one with 2000 soil types, 40 test conditions each (160 000 triaxial tests), and a second with 2048 soil types, 42 test conditions each (172 032 triaxial tests). Each dataset is a 4000×10 matrix applicable for multivariate forecasting and geotechnical simulations. In addition, a new Python code is introduced, empowering researchers to leverage Python packages for NorSand analyses.
Daniel Giles, Matthew M. Graham, Mosè Giordano, Tuomas Koskela, Alexandros Beskos, and Serge Guillas
Geosci. Model Dev., 17, 2427–2445, https://doi.org/10.5194/gmd-17-2427-2024, https://doi.org/10.5194/gmd-17-2427-2024, 2024
Short summary
Short summary
Digital twins of physical and human systems informed by real-time data are becoming ubiquitous across a wide range of settings. Progress for researchers is currently limited by a lack of tools to run these models effectively and efficiently. A key challenge is the optimal use of high-performance computing environments. The work presented here focuses on a developed open-source software platform which aims to improve this usage, with an emphasis placed on flexibility, efficiency, and scalability.
Stefan J. Miller, Paul A. Makar, and Colin J. Lee
Geosci. Model Dev., 17, 2197–2219, https://doi.org/10.5194/gmd-17-2197-2024, https://doi.org/10.5194/gmd-17-2197-2024, 2024
Short summary
Short summary
This work outlines a new solver written in Fortran to calculate the partitioning of metastable aerosols at thermodynamic equilibrium based on the forward algorithms of ISORROPIA II. The new code includes numerical improvements that decrease the computational speed (compared to ISORROPIA II) while improving the accuracy of the partitioning solution.
Fernanda Alvarado-Neves, Laurent Ailleres, Lachlan Grose, Alexander R. Cruden, and Robin Armit
Geosci. Model Dev., 17, 1975–1993, https://doi.org/10.5194/gmd-17-1975-2024, https://doi.org/10.5194/gmd-17-1975-2024, 2024
Short summary
Short summary
Previous work has demonstrated that adding geological knowledge to modelling methods creates more accurate and reliable models. Following this reasoning, we added constraints from magma emplacement mechanisms into existing modelling frameworks to improve the 3D characterisation of igneous intrusions. We tested the method on synthetic and real-world case studies, and the results show that our method can reproduce intrusion morphologies with no manual processing and using realistic datasets.
André R. Brodtkorb, Anna Benedictow, Heiko Klein, Arve Kylling, Agnes Nyiri, Alvaro Valdebenito, Espen Sollum, and Nina Kristiansen
Geosci. Model Dev., 17, 1957–1974, https://doi.org/10.5194/gmd-17-1957-2024, https://doi.org/10.5194/gmd-17-1957-2024, 2024
Short summary
Short summary
It is vital to know the extent and concentration of volcanic ash in the atmosphere during a volcanic eruption. Whilst satellite imagery may give an estimate of the ash right now (assuming no cloud coverage), we also need to know where it will be in the coming hours. This paper presents a method for estimating parameters for a volcanic eruption based on satellite observations of ash in the atmosphere. The software package is open source and applicable to similar inversion scenarios.
Kees Nederhoff, Maarten van Ormondt, Jay Veeramony, Ap van Dongeren, José Antonio Álvarez Antolínez, Tim Leijnse, and Dano Roelvink
Geosci. Model Dev., 17, 1789–1811, https://doi.org/10.5194/gmd-17-1789-2024, https://doi.org/10.5194/gmd-17-1789-2024, 2024
Short summary
Short summary
Forecasting tropical cyclones and their flooding impact is challenging. Our research introduces the Tropical Cyclone Forecasting Framework (TC-FF), enhancing cyclone predictions despite uncertainties. TC-FF generates global wind and flood scenarios, valuable even in data-limited regions. Applied to cases like Cyclone Idai, it showcases potential in bettering disaster preparation, marking progress in handling cyclone threats.
Younghun Kang and Ethan J. Kubatko
Geosci. Model Dev., 17, 1603–1625, https://doi.org/10.5194/gmd-17-1603-2024, https://doi.org/10.5194/gmd-17-1603-2024, 2024
Short summary
Short summary
Models used to simulate the flow of coastal and riverine waters, including flooding, require a geometric representation (or mesh) of geographic features that exhibit a range of disparate spatial scales from large, open waters to small, narrow channels. Representing these features in an accurate way without excessive computational overhead presents a challenge. Here, we develop an automatic mesh generation tool to help address this challenge. Our results demonstrate the efficacy of our approach.
Hui Wan, Kai Zhang, Christopher J. Vogl, Carol S. Woodward, Richard C. Easter, Philip J. Rasch, Yan Feng, and Hailong Wang
Geosci. Model Dev., 17, 1387–1407, https://doi.org/10.5194/gmd-17-1387-2024, https://doi.org/10.5194/gmd-17-1387-2024, 2024
Short summary
Short summary
Sophisticated numerical models of the Earth's atmosphere include representations of many physical and chemical processes. In numerical simulations, these processes need to be calculated in a certain sequence. This study reveals the weaknesses of the sequence of calculations used for aerosol processes in a global atmosphere model. A revision of the sequence is proposed and its impacts on the simulated global aerosol climatology are evaluated.
Christopher J. Vogl, Hui Wan, Carol S. Woodward, and Quan M. Bui
Geosci. Model Dev., 17, 1409–1428, https://doi.org/10.5194/gmd-17-1409-2024, https://doi.org/10.5194/gmd-17-1409-2024, 2024
Short summary
Short summary
Generally speaking, accurate climate simulation requires an accurate evolution of the underlying mathematical equations on large computers. The equations are typically formulated and evolved in process groups. Process coupling refers to how the evolution of each group is combined with that of other groups to evolve the full set of equations for the whole atmosphere. This work presents a mathematical framework to evaluate methods without the need to first implement the methods.
Tom Keel, Chris Brierley, and Tamsin Edwards
Geosci. Model Dev., 17, 1229–1247, https://doi.org/10.5194/gmd-17-1229-2024, https://doi.org/10.5194/gmd-17-1229-2024, 2024
Short summary
Short summary
Jet streams are an important control on surface weather as their speed and shape can modify the properties of weather systems. Establishing trends in the operation of jet streams may provide some indication of the future of weather in a warming world. Despite this, it has not been easy to establish trends, as many methods have been used to characterise them in data. We introduce a tool containing various implementations of jet stream statistics and algorithms that works in a standardised manner.
Amir Golparvar, Matthias Kästner, and Martin Thullner
Geosci. Model Dev., 17, 881–898, https://doi.org/10.5194/gmd-17-881-2024, https://doi.org/10.5194/gmd-17-881-2024, 2024
Short summary
Short summary
Coupled reaction transport modelling is an established and beneficial method for studying natural and synthetic porous material, with applications ranging from industrial processes to natural decompositions in terrestrial environments. Up to now, a framework that explicitly considers the porous structure (e.g. from µ-CT images) for modelling the transport of reactive species is missing. We presented a model that overcomes this limitation and represents a novel numerical simulation toolbox.
Stefan Hergarten
Geosci. Model Dev., 17, 781–794, https://doi.org/10.5194/gmd-17-781-2024, https://doi.org/10.5194/gmd-17-781-2024, 2024
Short summary
Short summary
The Voellmy rheology has been widely used for simulating snow and rock avalanches. Recently, a modified version of this rheology was proposed, which turned out to be able to predict the observed long runout of large rock avalanches theoretically. The software MinVoellmy presented here is the first numerical implementation of the modified rheology. It consists of MATLAB and Python classes, where simplicity and parsimony were the design goals.
Arjun Babu Nellikkattil, Danielle Lemmon, Travis Allen O'Brien, June-Yi Lee, and Jung-Eun Chu
Geosci. Model Dev., 17, 301–320, https://doi.org/10.5194/gmd-17-301-2024, https://doi.org/10.5194/gmd-17-301-2024, 2024
Short summary
Short summary
This study introduces a new computational framework called Scalable Feature Extraction and Tracking (SCAFET), designed to extract and track features in climate data. SCAFET stands out by using innovative shape-based metrics to identify features without relying on preconceived assumptions about the climate model or mean state. This approach allows more accurate comparisons between different models and scenarios.
Mohammad Mortezazadeh, Jean-François Cossette, Ashu Dastoor, Jean de Grandpré, Irena Ivanova, and Abdessamad Qaddouri
Geosci. Model Dev., 17, 335–346, https://doi.org/10.5194/gmd-17-335-2024, https://doi.org/10.5194/gmd-17-335-2024, 2024
Short summary
Short summary
The interpolation process is the most computationally expensive step of the semi-Lagrangian (SL) approach. In this paper we implement a new interpolation scheme into the semi-Lagrangian approach which has the same computational cost as a third-order polynomial scheme but with the accuracy of a fourth-order interpolation scheme. This improvement is achieved by using two third-order backward and forward polynomial interpolation schemes in two consecutive time steps.
Boris Gailleton, Luca C. Malatesta, Guillaume Cordonnier, and Jean Braun
Geosci. Model Dev., 17, 71–90, https://doi.org/10.5194/gmd-17-71-2024, https://doi.org/10.5194/gmd-17-71-2024, 2024
Short summary
Short summary
This contribution presents a new method to numerically explore the evolution of mountain ranges and surrounding areas. The method helps in monitoring with details on the timing and travel path of material eroded from the mountain ranges. It is particularly well suited to studies juxtaposing different domains – lakes or multiple rock types, for example – and enables the combination of different processes.
Denise Degen, Daniel Caviedes Voullième, Susanne Buiter, Harrie-Jan Hendricks Franssen, Harry Vereecken, Ana González-Nicolás, and Florian Wellmann
Geosci. Model Dev., 16, 7375–7409, https://doi.org/10.5194/gmd-16-7375-2023, https://doi.org/10.5194/gmd-16-7375-2023, 2023
Short summary
Short summary
In geosciences, we often use simulations based on physical laws. These simulations can be computationally expensive, which is a problem if simulations must be performed many times (e.g., to add error bounds). We show how a novel machine learning method helps to reduce simulation time. In comparison to other approaches, which typically only look at the output of a simulation, the method considers physical laws in the simulation itself. The method provides reliable results faster than standard.
Carlos Spa, Otilio Rojas, and Josep de la Puente
Geosci. Model Dev., 16, 7237–7252, https://doi.org/10.5194/gmd-16-7237-2023, https://doi.org/10.5194/gmd-16-7237-2023, 2023
Short summary
Short summary
This paper develops a calibration methodology of all absorbing techniques typically used by Fourier pseudo-spectral time-domain (PSTD) methods for geoacoustic wave simulations. The main contributions of the paper are (a) an implementation and quantitative comparison of all absorbing techniques available for PSTD methods through a simple and robust numerical experiment, and (b) a validation of these absorbing techniques in several 3-D seismic scenarios with gradual heterogeneity complexities.
Michael Hillier, Florian Wellmann, Eric A. de Kemp, Boyan Brodaric, Ernst Schetselaar, and Karine Bédard
Geosci. Model Dev., 16, 6987–7012, https://doi.org/10.5194/gmd-16-6987-2023, https://doi.org/10.5194/gmd-16-6987-2023, 2023
Short summary
Short summary
Neural networks can be used effectively to model three-dimensional geological structures from point data, sampling geological interfaces, units, and structural orientations. Existing neural network approaches for this type of modelling are advanced by the efficient incorporation of unconformities, new knowledge inputs, and improved data fitting techniques. These advances permit the modelling of more complex geology in diverse geological settings, different-sized areas, and various data regimes.
Ane Lopez de Gamiz-Zearra, Cecilie Hansen, Xavier Corrales, Iñaki Quincoces, Izaskun Preciado, and Eider Andonegi
EGUsphere, https://doi.org/10.5194/egusphere-2023-1368, https://doi.org/10.5194/egusphere-2023-1368, 2023
Short summary
Short summary
This paper describes the development of the first calibrated end-to-end Atlantis model for the Bay of Biscay that will help us improve the comprehension of the spatial functioning of the Bay of Biscay ecosystem and help establishing management measures of human activities. Our results highlighted the importance of lower trophic levels to the pelagic system and demonstrate the importance of having accurate and precise data for biological processes.
Tor Nordam, Ruben Kristiansen, Raymond Nepstad, Erik van Sebille, and Andy M. Booth
Geosci. Model Dev., 16, 5339–5363, https://doi.org/10.5194/gmd-16-5339-2023, https://doi.org/10.5194/gmd-16-5339-2023, 2023
Short summary
Short summary
We describe and compare two common methods, Eulerian and Lagrangian models, used to simulate the vertical transport of material in the ocean. They both solve the same transport problems but use different approaches for representing the underlying equations on the computer. The main focus of our study is on the numerical accuracy of the two approaches. Our results should be useful for other researchers creating or using these types of transport models.
Mathieu Gravey and Grégoire Mariethoz
Geosci. Model Dev., 16, 5265–5279, https://doi.org/10.5194/gmd-16-5265-2023, https://doi.org/10.5194/gmd-16-5265-2023, 2023
Short summary
Short summary
Multiple‐point geostatistics are widely used to simulate complex spatial structures based on a training image. The use of these methods relies on the possibility of finding optimal training images and parametrization of the simulation algorithms. Here, we propose finding an optimal set of parameters using only the training image as input. The main advantage of our approach is to remove the risk of overfitting an objective function.
Siting Li, Ping Wang, Hong Wang, Yue Peng, Zhaodong Liu, Wenjie Zhang, Hongli Liu, Yaqiang Wang, Huizheng Che, and Xiaoye Zhang
Geosci. Model Dev., 16, 4171–4191, https://doi.org/10.5194/gmd-16-4171-2023, https://doi.org/10.5194/gmd-16-4171-2023, 2023
Short summary
Short summary
Optimizing the initial state of atmospheric chemistry model input is one of the most essential methods to improve forecast accuracy. Considering the large computational load of the model, we introduce an ensemble optimal interpolation scheme (EnOI) for operational use and efficient updating of the initial fields of chemical components. The results suggest that EnOI provides a practical and cost-effective technique for improving the accuracy of chemical weather numerical forecasts.
Thomas Richter, Véronique Dansereau, Christian Lessig, and Piotr Minakowski
Geosci. Model Dev., 16, 3907–3926, https://doi.org/10.5194/gmd-16-3907-2023, https://doi.org/10.5194/gmd-16-3907-2023, 2023
Short summary
Short summary
Sea ice covers not only the pole regions but affects the weather and climate globally. For example, its white surface reflects more sunlight than land. The oceans around the poles are therefore kept cool, which affects the circulation in the oceans worldwide. Simulating the behavior and changes in sea ice on a computer is, however, very difficult. We propose a new computer simulation that better models how cracks in the ice change over time and show this by comparing to other simulations.
Emma J. MacKie, Michael Field, Lijing Wang, Zhen Yin, Nathan Schoedl, Matthew Hibbs, and Allan Zhang
Geosci. Model Dev., 16, 3765–3783, https://doi.org/10.5194/gmd-16-3765-2023, https://doi.org/10.5194/gmd-16-3765-2023, 2023
Short summary
Short summary
Earth scientists often have to fill in spatial gaps in measurements. This gap-filling or interpolation can be accomplished with geostatistical methods, where the statistical relationships between measurements are used to inform how these gaps should be filled. Despite the broad utility of these methods, there are few freely available geostatistical software applications. We present GStatSim, a Python package for performing different geostatistical interpolation methods.
Ian Madden, Simone Marras, and Jenny Suckale
Geosci. Model Dev., 16, 3479–3500, https://doi.org/10.5194/gmd-16-3479-2023, https://doi.org/10.5194/gmd-16-3479-2023, 2023
Short summary
Short summary
To aid risk managers who may wish to rapidly assess tsunami risk but may lack high-performance computing infrastructure, we provide an accessible software package able to rapidly model tsunami inundation over real topography by leveraging Google's Tensor Processing Unit, a high-performance hardware. Minimally trained users can take advantage of the rapid modeling abilities provided by this package via a web browser thanks to the ease of use of Google Cloud Platform.
Youtong Rong, Paul Bates, and Jeffrey Neal
Geosci. Model Dev., 16, 3291–3311, https://doi.org/10.5194/gmd-16-3291-2023, https://doi.org/10.5194/gmd-16-3291-2023, 2023
Short summary
Short summary
A novel subgrid channel (SGC) model is developed for river–floodplain modelling, allowing utilization of subgrid-scale bathymetric information while performing computations on relatively coarse grids. By including adaptive artificial diffusion, potential numerical instability, which the original SGC solver had, in low-friction regions such as urban areas is addressed. Evaluation of the new SGC model through structured tests confirmed that the accuracy and stability have improved.
Cited articles
Benn, D. I. and Astrom, J. A.: Calving glaciers and ice shelves, Adv. Phys.-X, 3, 1513819, https://doi.org/10.1080/23746149.2018.1513819, 2018. a
Biswas, R., Devine, K. D., and Flahert, J. E.: Parallel, adaptive finite element methods for conservation laws, Appl. Numer. Math., 14, 255–283, https://doi.org/10.1016/0168-9274(94)90029-9, 1994. a
Black, T. E. and Joughin, I.: Weekly to monthly terminus variability of Greenland's marine-terminating outlet glaciers, The Cryosphere, 17, 1–13, https://doi.org/10.5194/tc-17-1-2023, 2023. a, b
Bondzio, J., Morlighem, M., Seroussi, H., Kleiner, T., Ruckamp, M., Mouginot, J., Moon, T., Larour, E., and Humbert, A.: The mechanisms behind Jakobshavn Isbræ's acceleration and mass loss: A 3-D thermomechanical model study, Geophys. Res. Lett., 44, 6252–6260, https://doi.org/10.1002/2017GL073309, 2017. a
Bondzio, J. H., Seroussi, H., Morlighem, M., Kleiner, T., Rückamp, M., Humbert, A., and Larour, E. Y.: Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, The Cryosphere, 10, 497–510, https://doi.org/10.5194/tc-10-497-2016, 2016. a, b, c
Cheng, D., Hayes, W., Larour, E., Mohajerani, Y., Wood, M., Velicogna, I., and Rignot, E.: Calving Front Machine (CALFIN): glacial termini dataset and automated deep learning extraction method for Greenland, 1972–2019, The Cryosphere, 15, 1663–1675, https://doi.org/10.5194/tc-15-1663-2021, 2021. a
Cheng, G.: Animations of the evolution of misfits, Zenodo [video], https://doi.org/10.5281/zenodo.10454554, 2024. a
Cheng, G., Morlighem, M., Mouginot, J., and Cheng, D.: Helheim Glacier's Terminus Position Controls Its Seasonal and Inter-Annual Ice Flow Variability, Geophys. Res. Lett., 49, e2021GL097085, https://doi.org/10.1029/2021GL097085, 2022. a
Choi, Y., Morlighem, M., Wood, M., and Bondzio, J. H.: Comparison of four calving laws to model Greenland outlet glaciers, The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, 2018. a
Choi, Y., Morlighem, M., Rignot, E., and Wood, M.: Ice dynamics will remain a primary driver of Greenland ice sheet mass loss over the next century, Nature Commun. Earth Environ., 2, 26, https://doi.org/10.1038/s43247-021-00092-z, 2021. a, b
Crawford, A. J., Benn, D. I., Todd, J., Astrom, J. A., Bassis, J. N., and Zwinger, T.: Marine ice-cliff instability modeling shows mixed- mode ice-cliff failure and yields calving rate parameterization, Nat. Commun., 12, 2701, https://doi.org/10.1038/s41467-021-23070-7, 2021. a
dos Santos, T. D., Morlighem, M., and Seroussi, H.: Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18, Geosci. Model Dev., 14, 2545–2573, https://doi.org/10.5194/gmd-14-2545-2021, 2021. a, b, c
Eriksson, K., Estep, D., Hansbo, P., and Johnson, C.: Computational differential equations, Cambridge University Press, ISBN 9780521567381, 1996. a
Gibou, F., Fedkiw, R., and Osher, S.: A review of level-set methods and some recent applications, J. Comput. Phys., 353, 82–109, https://doi.org/10.1016/j.jcp.2017.10.006, 2018. a
Goelzer, H., Nowicki, S., Payne, A., Larour, E., Seroussi, H., Lipscomb, W. H., Gregory, J., Abe-Ouchi, A., Shepherd, A., Simon, E., Agosta, C., Alexander, P., Aschwanden, A., Barthel, A., Calov, R., Chambers, C., Choi, Y., Cuzzone, J., Dumas, C., Edwards, T., Felikson, D., Fettweis, X., Golledge, N. R., Greve, R., Humbert, A., Huybrechts, P., Le clec'h, S., Lee, V., Leguy, G., Little, C., Lowry, D. P., Morlighem, M., Nias, I., Quiquet, A., Rückamp, M., Schlegel, N.-J., Slater, D. A., Smith, R. S., Straneo, F., Tarasov, L., van de Wal, R., and van den Broeke, M.: The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6, The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, 2020. a
Gong, C.: enigne/Levelset: Level-set Stabilization (v1.1), Zenodo [code and data set], https://doi.org/10.5281/zenodo.10454657, 2024. a
Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Advances in Geophysical and Environmental Mechanics and Mathematics, Springer Science & Business Media, https://doi.org/10.1007/978-3-642-03415-2, 2009. a
Gross, S. and Reusken, A.: Numerical Methods for Two-phase Incompressible Flows, vol. 40 of Springer Series in Computational Mathematics, 1–472, https://doi.org/10.1007/978-3-642-19686-7, 2011. a
Gudmundsson, G. H., Paolo, F. S., Adusumilli, S., and Fricker, H. A.: Instantaneous Antarctic ice sheet mass loss driven by thinning ice shelves, Geophys. Res. Lett., 46, 13903–13909, https://doi.org/10.1029/2019GL085027, 2019. a, b
Gudmundsson, H.: GHilmarG/UaSource: Ua2019b, Zenodo [code], https://doi.org/10.5281/zenodo.3706624, 2020. a, b, c
Hartmann, D., Meinke, M., and Schroeder, W.: The constrained reinitialization equation for level set methods, J. Comput. Phys., 229, 1514–1535, https://doi.org/10.1016/j.jcp.2009.10.042, 2010. a, b
Henri, F., Coquerelle, M., and Lubin, P.: Geometrical level set reinitialization using closest point method and kink detection for thin filaments, topology changes and two-phase flows, J. Comput. Phys., 448, 110704, https://doi.org/10.1016/j.jcp.2021.110704, 2022. a, b
IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, vol. In Press, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, https://doi.org/10.1017/9781009157896, 2021. a
ISSM Team: Ice-sheet and Sea-level System Model source code, v4.23 r27696, Zenodo [code], https://doi.org/10.5281/zenodo.7850841, 2023. a, b
King, M. D., Howat, I. M., Candela, S. G., Noh, M. J., Jeong, S., Noel, B. P. Y., van den Broeke, M. R., Wouters, B., and Negrete, A.: Dynamic ice loss from the Greenland Ice Sheet driven by sustained glacier retreat, Commun. Earth Environ., 1, 1, https://doi.org/10.1038/s43247-020-0001-2, 2020. a
Larour, E., Seroussi, H., Morlighem, M., and Rignot, E.: Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM), J. Geophys. Res., 117, 1–20, https://doi.org/10.1029/2011JF002140, 2012. a
Larson, M. G. and Bengzon, F.: The Finite Element Method: Theory, Implementation, and Applications, Springer Publishing Company, Incorporated, https://doi.org/10.1007/978-3-642-33287-6, 2013. a, b, c, d
Li, C., Xu, C., Gui, C., and Fox, M. D.: Level set evolution without re-initialization: A new variational formulation, in: 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Vol 1, Proceedings, edited by: Schmid, C., Soatto, S., and Tomasi, C., PROCEEDINGS – IEEE Computer Society Conference On Computer Vision And Pattern Recognition, 430–436, IEEE Comp Soc, conference on Computer Vision and Pattern Recognition, 20–25 June 2005, San Diego, CA, 2005. a
MacAyeal, D. R.: Large-scale ice flow over a viscous basal sediment: Theory and application to Ice Stream B, Antarctica, J. Geophys. Res., 94, 4071–4087, 1989. a
Mouginot, J., Rignot, E., Bjørk, A. A., van den Broeke, M., Millan, R., Morlighem, M., Noël, B., Scheuchl, B., and Wood, M.: Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018, P. Natl. Acad. Sci. USA., 116, 9239–9244, https://doi.org/10.1073/pnas.1904242116, 2019. a, b
Pattyn, F. and Morlighem, M.: The uncertain future of the Antarctic Ice Sheet, Science, 367, 1331–1335, https://doi.org/10.1126/science.aaz5487, 2020. a
Pollard, D., DeConto, R. M., and Alley, R. B.: Potential Antarctic Ice Sheet retreat driven by hydrofracturing and ice cliff failure, Earth Planet Sc. Lett., 412, 112–121, https://doi.org/10.1016/j.epsl.2014.12.035, 2015. a
Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070, https://doi.org/10.5194/tc-14-3033-2020, 2020. a
Sussman, M., Smereka, P., and Osher, S.: A Level set Aspproach for Computing Solutions to Incompressible Two-Phase Flow, J. Comput. Phys., 114, 146–159, https://doi.org/10.1006/jcph.1994.1155, 1994. a
Toure, M. K. and Soulaimani, A.: Stabilized finite element methods for solving the level set equation without reinitialization, Comput. Math. Appl., 71, 1602–1623, https://doi.org/10.1016/j.camwa.2016.02.028, 2016. a, b
Wood, M., Rignot, E., Fenty, I., An, L., Bjørk, A., van den Broeke, M., Cai, C., Kane, E., Menemenlis, D., Millan, R., Morlighem, M., Mouginot, J., Noël, B., Scheuchl, B., Velicogna, I., Willis, J. K., and Zhang, H.: Ocean forcing drives glacier retreat in Greenland, Sci. Adv., 7, 1315–1332, https://doi.org/10.1126/sciadv.aba7282, 2021. a
Short summary
We conducted a comprehensive analysis of the stabilization and reinitialization techniques currently employed in ISSM and Úa for solving level-set equations, specifically those related to the dynamic representation of moving ice fronts within numerical ice sheet models. Our results demonstrate that the streamline upwind Petrov–Galerkin (SUPG) method outperforms the other approaches. We found that excessively frequent reinitialization can lead to exceptionally high errors in simulations.
We conducted a comprehensive analysis of the stabilization and reinitialization techniques...
