Articles | Volume 13, issue 3
https://doi.org/10.5194/gmd-13-1155-2020
© Author(s) 2020. 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-13-1155-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper
| 
12 Mar 2020
Development and technical paper |
| 12 Mar 2020
| 12 Mar 2020
A new open-source viscoelastic solid earth deformation module implemented in Elmer (v8.4)
CSC–IT Center for Science Ltd., Espoo, Finland
Grace A. Nield
Surveying and Spatial Sciences, School of Technology, Environments and Design, University of Tasmania, Hobart, Australia
Department of Geography, Durham University, Durham, UK
Juha Ruokolainen
CSC–IT Center for Science Ltd., Espoo, Finland
Matt A. King
Surveying and Spatial Sciences, School of Technology, Environments and Design, University of Tasmania, Hobart, Australia
Related authors
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Ben Galton-Fenzi, and Poul Christoffersen
EGUsphere, https://doi.org/10.5194/egusphere-2024-1005, https://doi.org/10.5194/egusphere-2024-1005, 2024
Short summary
Short summary
Our research delves into the future ice loss in Antarctica’s Wilkes Subglacial Basin (WSB) and its impact on sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. According to our best model results, under high-emission scenarios, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300, potentially raising global sea levels by up to 0.34 m by 2500.
Iain Wheel, Douglas I. Benn, Anna J. Crawford, Joe A. Todd, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2023-2778, https://doi.org/10.5194/egusphere-2023-2778, 2024
Short summary
Short summary
Calving, the detachment of large icebergs from glaciers, is one of the largest uncertainties in future sea level rise projections. This process is poorly understood and there is an absence of detailed models capable of simulating calving. A new 3D calving model has been developed to better understand calving at glaciers where detailed modelling was previously limited. Importantly, the new model is very flexible. By allowing for unrestricted calving geometries, it can be applied at any location.
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.
Denis Cohen, Guillaume Jouvet, Thomas Zwinger, Angela Landgraf, and Urs H. Fischer
E&G Quaternary Sci. J., 72, 189–201, https://doi.org/10.5194/egqsj-72-189-2023, https://doi.org/10.5194/egqsj-72-189-2023, 2023
Short summary
Short summary
During glacial times in Switzerland, glaciers of the Alps excavated valleys in low-lying regions that were later filled with sediment or water. How glaciers eroded these valleys is not well understood because erosion occurred near ice margins where ice moved slowly and was present for short times. Erosion is linked to the speed of ice and to water flowing under it. Here we present a model that estimates the location of water channels beneath the ice and links these locations to zones of erosion.
André Löfgren, Thomas Zwinger, Peter Råback, Christian Helanow, and Josefin Ahlkrona
EGUsphere, https://doi.org/10.5194/egusphere-2023-1507, https://doi.org/10.5194/egusphere-2023-1507, 2023
Short summary
Short summary
This paper investigates a stabilization method for free-surface flows in the context of glacier simulations. Previous applications of the stabilization on ice flows have only considered simple ice-sheet benchmark problems; in particular the method has not been tested on real-world glacier domains. This work addresses this shortcoming by demonstrating that the stabilization works well also in this case, and increases stability and robustness without negatively impacting computation times.
Maryam Zarrinderakht, Christian Schoof, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2023-807, https://doi.org/10.5194/egusphere-2023-807, 2023
Short summary
Short summary
We used a model to study how crevasses propagate in ice shelves. Our model combines a viscous model and a fracture mechanics model. We studied periodic crevasses on an ice shelf being stretched. We show that existing models based only on stress cannot fully explain how crevasses grow and lead to iceberg calving. This model can be a useful tool to train a low-dimensional representation calving law for an ice sheet model.
Douglas I. Benn, Adrian Luckman, Jan A. Åström, Anna J. Crawford, Stephen L. Cornford, Suzanne L. Bevan, Thomas Zwinger, Rupert Gladstone, Karen Alley, Erin Pettit, and Jeremy Bassis
The Cryosphere, 16, 2545–2564, https://doi.org/10.5194/tc-16-2545-2022, https://doi.org/10.5194/tc-16-2545-2022, 2022
Short summary
Short summary
Thwaites Glacier (TG), in West Antarctica, is potentially unstable and may contribute significantly to sea-level rise as global warming continues. Using satellite data, we show that Thwaites Eastern Ice Shelf, the largest remaining floating extension of TG, has started to accelerate as it fragments along a shear zone. Computer modelling does not indicate that fragmentation will lead to imminent glacier collapse, but it is clear that major, rapid, and unpredictable changes are underway.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev., 14, 889–905, https://doi.org/10.5194/gmd-14-889-2021, https://doi.org/10.5194/gmd-14-889-2021, 2021
Short summary
Short summary
Retreat of the Antarctic ice sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
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.
Shahbaz Memon, Dorothée Vallot, Thomas Zwinger, Jan Åström, Helmut Neukirchen, Morris Riedel, and Matthias Book
Geosci. Model Dev., 12, 3001–3015, https://doi.org/10.5194/gmd-12-3001-2019, https://doi.org/10.5194/gmd-12-3001-2019, 2019
Short summary
Short summary
Scientific workflows enable complex scientific computational scenarios, which include data intensive scenarios, parametric executions, and interactive simulations. In this article, we applied the UNICORE workflow management system to automate a formerly hard-coded coupling of a glacier flow model and a calving model, which contain many tasks and dependencies, ranging from pre-processing and data management to repetitive executions on heterogeneous high-performance computing (HPC) resources.
Joe Todd, Poul Christoffersen, Thomas Zwinger, Peter Råback, and Douglas I. Benn
The Cryosphere, 13, 1681–1694, https://doi.org/10.5194/tc-13-1681-2019, https://doi.org/10.5194/tc-13-1681-2019, 2019
Short summary
Short summary
The Greenland Ice Sheet loses 30 %–60 % of its ice due to iceberg calving. Calving processes and their links to climate are not well understood or incorporated into numerical models of glaciers. Here we use a new 3-D calving model to investigate calving at Store Glacier, West Greenland, and test its sensitivity to increased submarine melting and reduced support from ice mélange (sea ice and icebergs). We find Store remains fairly stable despite these changes, but less so in the southern side.
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.
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.
Sue Cook, Jan Åström, Thomas Zwinger, Benjamin Keith Galton-Fenzi, Jamin Stevens Greenbaum, and Richard Coleman
The Cryosphere, 12, 2401–2411, https://doi.org/10.5194/tc-12-2401-2018, https://doi.org/10.5194/tc-12-2401-2018, 2018
Short summary
Short summary
The growth of fractures on Antarctic ice shelves is important because it controls the amount of ice lost as icebergs. We use a model constructed of multiple interconnected blocks to predict the locations where fractures will form on the Totten Ice Shelf in East Antarctica. The results show that iceberg calving is controlled not only by fractures forming near the front of the ice shelf but also by fractures which formed many kilometres upstream.
Yongmei Gong, Thomas Zwinger, Jan Åström, Bas Altena, Thomas Schellenberger, Rupert Gladstone, and John C. Moore
The Cryosphere, 12, 1563–1577, https://doi.org/10.5194/tc-12-1563-2018, https://doi.org/10.5194/tc-12-1563-2018, 2018
Short summary
Short summary
In this study we apply a discrete element model capable of simulating ice fracturing. A microscopic-scale discrete process is applied in addition to a continuum ice dynamics model to investigate the mechanisms facilitated by basal meltwater production, surface meltwater and ice crack opening, for the surge in Basin 3, Austfonna ice cap. The discrete element model is used to locate the ice cracks that can penetrate though the full thickness of the glacier and deliver surface water to the bed.
Dorothée Vallot, Jan Åström, Thomas Zwinger, Rickard Pettersson, Alistair Everett, Douglas I. Benn, Adrian Luckman, Ward J. J. van Pelt, Faezeh Nick, and Jack Kohler
The Cryosphere, 12, 609–625, https://doi.org/10.5194/tc-12-609-2018, https://doi.org/10.5194/tc-12-609-2018, 2018
Short summary
Short summary
This paper presents a new perspective on the role of ice dynamics and ocean interaction in glacier calving processes applied to Kronebreen, a tidewater glacier in Svalbard. A global modelling approach includes ice flow modelling, undercutting estimation by a combination of glacier energy balance and plume modelling as well as calving by a discrete particle model. We show that modelling undercutting is necessary and calving is influenced by basal friction velocity and geometry.
Hakime Seddik, Ralf Greve, Thomas Zwinger, and Shin Sugiyama
The Cryosphere, 11, 2213–2229, https://doi.org/10.5194/tc-11-2213-2017, https://doi.org/10.5194/tc-11-2213-2017, 2017
Short summary
Short summary
The Shirase Glacier in Antarctica is studied by means of a computer model. This model implements two physical approaches to represent the glacier flow dynamics. This study finds that it is important to use the more precise and sophisticated method in order to better understand and predict the evolution of fast flowing glaciers. This may be important to more accurately predict the sea level change due to global warming.
Rupert Michael Gladstone, Roland Charles Warner, Benjamin Keith Galton-Fenzi, Olivier Gagliardini, Thomas Zwinger, and Ralf Greve
The Cryosphere, 11, 319–329, https://doi.org/10.5194/tc-11-319-2017, https://doi.org/10.5194/tc-11-319-2017, 2017
Short summary
Short summary
Computer models are used to simulate the behaviour of glaciers and ice sheets. It has been found that such models are required to be run at very high resolution (which means high computational expense) in order to accurately represent the evolution of marine ice sheets (ice sheets resting on bedrock below sea level), in certain situations which depend on sub-glacial physical processes.
T. Zwinger, T. Malm, M. Schäfer, R. Stenberg, and J. C. Moore
The Cryosphere, 9, 1415–1426, https://doi.org/10.5194/tc-9-1415-2015, https://doi.org/10.5194/tc-9-1415-2015, 2015
Short summary
Short summary
By deploying a large-scale high-resolution turbulent CFD simulation using the present-day topography of the Scharffenbergbotnen (SBB) valley, we show how the surrounding topography redirects incoming easterly katabatic storm fronts to impact the blue ice areas (BIA) inside the valley, where the snow cover frequently is removed. A further simulation of a reconstructed topography at the Late Glacial Maximum further reveals that the BIA at SBB must have formed after this period.
M. Schäfer, F. Gillet-Chaulet, R. Gladstone, R. Pettersson, V. A. Pohjola, T. Strozzi, and T. Zwinger
The Cryosphere, 8, 1951–1973, https://doi.org/10.5194/tc-8-1951-2014, https://doi.org/10.5194/tc-8-1951-2014, 2014
R. Gladstone, M. Schäfer, T. Zwinger, Y. Gong, T. Strozzi, R. Mottram, F. Boberg, and J. C. Moore
The Cryosphere, 8, 1393–1405, https://doi.org/10.5194/tc-8-1393-2014, https://doi.org/10.5194/tc-8-1393-2014, 2014
B. Sun, J. C. Moore, T. Zwinger, L. Zhao, D. Steinhage, X. Tang, D. Zhang, X. Cui, and C. Martín
The Cryosphere, 8, 1121–1128, https://doi.org/10.5194/tc-8-1121-2014, https://doi.org/10.5194/tc-8-1121-2014, 2014
S. Cook, I. C. Rutt, T. Murray, A. Luckman, T. Zwinger, N. Selmes, A. Goldsack, and T. D. James
The Cryosphere, 8, 827–841, https://doi.org/10.5194/tc-8-827-2014, https://doi.org/10.5194/tc-8-827-2014, 2014
T. Zwinger, M. Schäfer, C. Martín, and J. C. Moore
The Cryosphere, 8, 607–621, https://doi.org/10.5194/tc-8-607-2014, https://doi.org/10.5194/tc-8-607-2014, 2014
T. Sato, T. Shiraiwa, R. Greve, H. Seddik, E. Edelmann, and T. Zwinger
Clim. Past, 10, 393–404, https://doi.org/10.5194/cp-10-393-2014, https://doi.org/10.5194/cp-10-393-2014, 2014
B. de Fleurian, O. Gagliardini, T. Zwinger, G. Durand, E. Le Meur, D. Mair, and P. Råback
The Cryosphere, 8, 137–153, https://doi.org/10.5194/tc-8-137-2014, https://doi.org/10.5194/tc-8-137-2014, 2014
J. A. Åström, T. I. Riikilä, T. Tallinen, T. Zwinger, D. Benn, J. C. Moore, and J. Timonen
The Cryosphere, 7, 1591–1602, https://doi.org/10.5194/tc-7-1591-2013, https://doi.org/10.5194/tc-7-1591-2013, 2013
O. Gagliardini, T. Zwinger, F. Gillet-Chaulet, G. Durand, L. Favier, B. de Fleurian, R. Greve, M. Malinen, C. Martín, P. Råback, J. Ruokolainen, M. Sacchettini, M. Schäfer, H. Seddik, and J. Thies
Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, https://doi.org/10.5194/gmd-6-1299-2013, 2013
A. S. Drouet, D. Docquier, G. Durand, R. Hindmarsh, F. Pattyn, O. Gagliardini, and T. Zwinger
The Cryosphere, 7, 395–406, https://doi.org/10.5194/tc-7-395-2013, https://doi.org/10.5194/tc-7-395-2013, 2013
L. Zhao, L. Tian, T. Zwinger, R. Ding, J. Zong, Q. Ye, and J. C. Moore
The Cryosphere Discuss., https://doi.org/10.5194/tcd-7-145-2013, https://doi.org/10.5194/tcd-7-145-2013, 2013
Revised manuscript not accepted
F. Gillet-Chaulet, O. Gagliardini, H. Seddik, M. Nodet, G. Durand, C. Ritz, T. Zwinger, R. Greve, and D. G. Vaughan
The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, https://doi.org/10.5194/tc-6-1561-2012, 2012
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Ben Galton-Fenzi, and Poul Christoffersen
EGUsphere, https://doi.org/10.5194/egusphere-2024-1005, https://doi.org/10.5194/egusphere-2024-1005, 2024
Short summary
Short summary
Our research delves into the future ice loss in Antarctica’s Wilkes Subglacial Basin (WSB) and its impact on sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. According to our best model results, under high-emission scenarios, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300, potentially raising global sea levels by up to 0.34 m by 2500.
Iain Wheel, Douglas I. Benn, Anna J. Crawford, Joe A. Todd, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2023-2778, https://doi.org/10.5194/egusphere-2023-2778, 2024
Short summary
Short summary
Calving, the detachment of large icebergs from glaciers, is one of the largest uncertainties in future sea level rise projections. This process is poorly understood and there is an absence of detailed models capable of simulating calving. A new 3D calving model has been developed to better understand calving at glaciers where detailed modelling was previously limited. Importantly, the new model is very flexible. By allowing for unrestricted calving geometries, it can be applied at any location.
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.
Denis Cohen, Guillaume Jouvet, Thomas Zwinger, Angela Landgraf, and Urs H. Fischer
E&G Quaternary Sci. J., 72, 189–201, https://doi.org/10.5194/egqsj-72-189-2023, https://doi.org/10.5194/egqsj-72-189-2023, 2023
Short summary
Short summary
During glacial times in Switzerland, glaciers of the Alps excavated valleys in low-lying regions that were later filled with sediment or water. How glaciers eroded these valleys is not well understood because erosion occurred near ice margins where ice moved slowly and was present for short times. Erosion is linked to the speed of ice and to water flowing under it. Here we present a model that estimates the location of water channels beneath the ice and links these locations to zones of erosion.
André Löfgren, Thomas Zwinger, Peter Råback, Christian Helanow, and Josefin Ahlkrona
EGUsphere, https://doi.org/10.5194/egusphere-2023-1507, https://doi.org/10.5194/egusphere-2023-1507, 2023
Short summary
Short summary
This paper investigates a stabilization method for free-surface flows in the context of glacier simulations. Previous applications of the stabilization on ice flows have only considered simple ice-sheet benchmark problems; in particular the method has not been tested on real-world glacier domains. This work addresses this shortcoming by demonstrating that the stabilization works well also in this case, and increases stability and robustness without negatively impacting computation times.
Maryam Zarrinderakht, Christian Schoof, and Thomas Zwinger
EGUsphere, https://doi.org/10.5194/egusphere-2023-807, https://doi.org/10.5194/egusphere-2023-807, 2023
Short summary
Short summary
We used a model to study how crevasses propagate in ice shelves. Our model combines a viscous model and a fracture mechanics model. We studied periodic crevasses on an ice shelf being stretched. We show that existing models based only on stress cannot fully explain how crevasses grow and lead to iceberg calving. This model can be a useful tool to train a low-dimensional representation calving law for an ice sheet model.
Douglas I. Benn, Adrian Luckman, Jan A. Åström, Anna J. Crawford, Stephen L. Cornford, Suzanne L. Bevan, Thomas Zwinger, Rupert Gladstone, Karen Alley, Erin Pettit, and Jeremy Bassis
The Cryosphere, 16, 2545–2564, https://doi.org/10.5194/tc-16-2545-2022, https://doi.org/10.5194/tc-16-2545-2022, 2022
Short summary
Short summary
Thwaites Glacier (TG), in West Antarctica, is potentially unstable and may contribute significantly to sea-level rise as global warming continues. Using satellite data, we show that Thwaites Eastern Ice Shelf, the largest remaining floating extension of TG, has started to accelerate as it fragments along a shear zone. Computer modelling does not indicate that fragmentation will lead to imminent glacier collapse, but it is clear that major, rapid, and unpredictable changes are underway.
Ole Richter, David E. Gwyther, Matt A. King, and Benjamin K. Galton-Fenzi
The Cryosphere, 16, 1409–1429, https://doi.org/10.5194/tc-16-1409-2022, https://doi.org/10.5194/tc-16-1409-2022, 2022
Short summary
Short summary
Tidal currents may play an important role in Antarctic ice sheet retreat by changing the rate at which the ocean melts glaciers. Here, using a computational ocean model, we derive the first estimate of present-day tidal melting that covers all of Antarctica. Our results suggest that large-scale ocean models aiming to accurately predict ice melt rates will need to account for the effects of tides. The inclusion of tide-induced friction at the ice–ocean interface should be prioritized.
Grace A. Nield, Matt A. King, Rebekka Steffen, and Bas Blank
Geosci. Model Dev., 15, 2489–2503, https://doi.org/10.5194/gmd-15-2489-2022, https://doi.org/10.5194/gmd-15-2489-2022, 2022
Short summary
Short summary
We present a finite-element model of post-seismic solid Earth deformation built in the software package Abaqus for the purpose of calculating post-seismic deformation in the far field of major earthquakes. The model is benchmarked against an existing open-source post-seismic model demonstrating good agreement. The advantage over existing models is the potential for simple modification to include 3-D Earth structure, non-linear rheologies and alternative or multiple sources of stress change.
Steven J. Phipps, Jason L. Roberts, and Matt A. King
Geosci. Model Dev., 14, 5107–5124, https://doi.org/10.5194/gmd-14-5107-2021, https://doi.org/10.5194/gmd-14-5107-2021, 2021
Short summary
Short summary
Simplified schemes, known as parameterisations, are sometimes used to describe physical processes within numerical models. However, the values of the parameters are uncertain. This introduces uncertainty into the model outputs. We develop a simple approach to identify plausible ranges for model parameters. Using a model of the Antarctic Ice Sheet, we find that the value of one parameter can depend on the values of others. We conclude that a single optimal set of parameter values does not exist.
Rupert Gladstone, Benjamin Galton-Fenzi, David Gwyther, Qin Zhou, Tore Hattermann, Chen Zhao, Lenneke Jong, Yuwei Xia, Xiaoran Guo, Konstantinos Petrakopoulos, Thomas Zwinger, Daniel Shapero, and John Moore
Geosci. Model Dev., 14, 889–905, https://doi.org/10.5194/gmd-14-889-2021, https://doi.org/10.5194/gmd-14-889-2021, 2021
Short summary
Short summary
Retreat of the Antarctic ice sheet, and hence its contribution to sea level rise, is highly sensitive to melting of its floating ice shelves. This melt is caused by warm ocean currents coming into contact with the ice. Computer models used for future ice sheet projections are not able to realistically evolve these melt rates. We describe a new coupling framework to enable ice sheet and ocean computer models to interact, allowing projection of the evolution of melt and its impact on sea level.
Bogdan Matviichuk, Matt King, and Christopher Watson
Solid Earth, 11, 1849–1863, https://doi.org/10.5194/se-11-1849-2020, https://doi.org/10.5194/se-11-1849-2020, 2020
Short summary
Short summary
The Earth deforms as the weight of ocean mass changes with the tides. GPS has been used to estimate displacements of the Earth at tidal periods and then used to understand the properties of the Earth or to test models of ocean tides. However, there are important inaccuracies in these GPS measurements at major tidal periods. We find that combining GPS and GLONASS gives more accurate results for constituents other than K2 and K1; for these, GLONASS or ambiguity resolved GPS are preferred.
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.
Shahbaz Memon, Dorothée Vallot, Thomas Zwinger, Jan Åström, Helmut Neukirchen, Morris Riedel, and Matthias Book
Geosci. Model Dev., 12, 3001–3015, https://doi.org/10.5194/gmd-12-3001-2019, https://doi.org/10.5194/gmd-12-3001-2019, 2019
Short summary
Short summary
Scientific workflows enable complex scientific computational scenarios, which include data intensive scenarios, parametric executions, and interactive simulations. In this article, we applied the UNICORE workflow management system to automate a formerly hard-coded coupling of a glacier flow model and a calving model, which contain many tasks and dependencies, ranging from pre-processing and data management to repetitive executions on heterogeneous high-performance computing (HPC) resources.
Joe Todd, Poul Christoffersen, Thomas Zwinger, Peter Råback, and Douglas I. Benn
The Cryosphere, 13, 1681–1694, https://doi.org/10.5194/tc-13-1681-2019, https://doi.org/10.5194/tc-13-1681-2019, 2019
Short summary
Short summary
The Greenland Ice Sheet loses 30 %–60 % of its ice due to iceberg calving. Calving processes and their links to climate are not well understood or incorporated into numerical models of glaciers. Here we use a new 3-D calving model to investigate calving at Store Glacier, West Greenland, and test its sensitivity to increased submarine melting and reduced support from ice mélange (sea ice and icebergs). We find Store remains fairly stable despite these changes, but less so in the southern side.
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.
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.
Lenneke M. Jong, Rupert M. Gladstone, Benjamin K. Galton-Fenzi, and Matt A. King
The Cryosphere, 12, 2425–2436, https://doi.org/10.5194/tc-12-2425-2018, https://doi.org/10.5194/tc-12-2425-2018, 2018
Short summary
Short summary
We used an ice sheet model to simulate temporary regrounding of a marine ice sheet retreating across a retrograde bedrock slope. We show that a sliding relation incorporating water-filled cavities and the ice overburden pressure at the base allows the temporary regrounding to occur. This suggests that choice of basal sliding relation can be important when modelling grounding line behaviour of regions where potential ice rises and pinning points are present and regrounding could occur.
Sue Cook, Jan Åström, Thomas Zwinger, Benjamin Keith Galton-Fenzi, Jamin Stevens Greenbaum, and Richard Coleman
The Cryosphere, 12, 2401–2411, https://doi.org/10.5194/tc-12-2401-2018, https://doi.org/10.5194/tc-12-2401-2018, 2018
Short summary
Short summary
The growth of fractures on Antarctic ice shelves is important because it controls the amount of ice lost as icebergs. We use a model constructed of multiple interconnected blocks to predict the locations where fractures will form on the Totten Ice Shelf in East Antarctica. The results show that iceberg calving is controlled not only by fractures forming near the front of the ice shelf but also by fractures which formed many kilometres upstream.
Yongmei Gong, Thomas Zwinger, Jan Åström, Bas Altena, Thomas Schellenberger, Rupert Gladstone, and John C. Moore
The Cryosphere, 12, 1563–1577, https://doi.org/10.5194/tc-12-1563-2018, https://doi.org/10.5194/tc-12-1563-2018, 2018
Short summary
Short summary
In this study we apply a discrete element model capable of simulating ice fracturing. A microscopic-scale discrete process is applied in addition to a continuum ice dynamics model to investigate the mechanisms facilitated by basal meltwater production, surface meltwater and ice crack opening, for the surge in Basin 3, Austfonna ice cap. The discrete element model is used to locate the ice cracks that can penetrate though the full thickness of the glacier and deliver surface water to the bed.
Dorothée Vallot, Jan Åström, Thomas Zwinger, Rickard Pettersson, Alistair Everett, Douglas I. Benn, Adrian Luckman, Ward J. J. van Pelt, Faezeh Nick, and Jack Kohler
The Cryosphere, 12, 609–625, https://doi.org/10.5194/tc-12-609-2018, https://doi.org/10.5194/tc-12-609-2018, 2018
Short summary
Short summary
This paper presents a new perspective on the role of ice dynamics and ocean interaction in glacier calving processes applied to Kronebreen, a tidewater glacier in Svalbard. A global modelling approach includes ice flow modelling, undercutting estimation by a combination of glacier energy balance and plume modelling as well as calving by a discrete particle model. We show that modelling undercutting is necessary and calving is influenced by basal friction velocity and geometry.
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.
Wolfgang Rack, Matt A. King, Oliver J. Marsh, Christian T. Wild, and Dana Floricioiu
The Cryosphere, 11, 2481–2490, https://doi.org/10.5194/tc-11-2481-2017, https://doi.org/10.5194/tc-11-2481-2017, 2017
Short summary
Short summary
Predicting changes of the Antarctic Ice Sheet involves fully understanding ice dynamics at the transition between grounded and floating ice. We map tidal bending of ice by satellite using InSAR, and we use precise GPS measurements with assumptions of tidal elastic bending to better interpret the satellite signal. It allows us to better define the grounding-line position and to refine the shape of tidal flexure profiles.
Hakime Seddik, Ralf Greve, Thomas Zwinger, and Shin Sugiyama
The Cryosphere, 11, 2213–2229, https://doi.org/10.5194/tc-11-2213-2017, https://doi.org/10.5194/tc-11-2213-2017, 2017
Short summary
Short summary
The Shirase Glacier in Antarctica is studied by means of a computer model. This model implements two physical approaches to represent the glacier flow dynamics. This study finds that it is important to use the more precise and sophisticated method in order to better understand and predict the evolution of fast flowing glaciers. This may be important to more accurately predict the sea level change due to global warming.
Riccardo E. M. Riva, Thomas Frederikse, Matt A. King, Ben Marzeion, and Michiel R. van den Broeke
The Cryosphere, 11, 1327–1332, https://doi.org/10.5194/tc-11-1327-2017, https://doi.org/10.5194/tc-11-1327-2017, 2017
Short summary
Short summary
The reduction of ice masses stored on land has made an important contribution to sea-level rise over the last century, as well as changed the Earth's shape. We model the solid-earth response to ice mass changes and find significant vertical deformation signals over large continental areas. We show how deformation rates have varied strongly throughout the last century, which affects the interpretation and extrapolation of recent observations of vertical land motion and sea-level change.
Rupert Michael Gladstone, Roland Charles Warner, Benjamin Keith Galton-Fenzi, Olivier Gagliardini, Thomas Zwinger, and Ralf Greve
The Cryosphere, 11, 319–329, https://doi.org/10.5194/tc-11-319-2017, https://doi.org/10.5194/tc-11-319-2017, 2017
Short summary
Short summary
Computer models are used to simulate the behaviour of glaciers and ice sheets. It has been found that such models are required to be run at very high resolution (which means high computational expense) in order to accurately represent the evolution of marine ice sheets (ice sheets resting on bedrock below sea level), in certain situations which depend on sub-glacial physical processes.
T. Zwinger, T. Malm, M. Schäfer, R. Stenberg, and J. C. Moore
The Cryosphere, 9, 1415–1426, https://doi.org/10.5194/tc-9-1415-2015, https://doi.org/10.5194/tc-9-1415-2015, 2015
Short summary
Short summary
By deploying a large-scale high-resolution turbulent CFD simulation using the present-day topography of the Scharffenbergbotnen (SBB) valley, we show how the surrounding topography redirects incoming easterly katabatic storm fronts to impact the blue ice areas (BIA) inside the valley, where the snow cover frequently is removed. A further simulation of a reconstructed topography at the Late Glacial Maximum further reveals that the BIA at SBB must have formed after this period.
M. Schäfer, F. Gillet-Chaulet, R. Gladstone, R. Pettersson, V. A. Pohjola, T. Strozzi, and T. Zwinger
The Cryosphere, 8, 1951–1973, https://doi.org/10.5194/tc-8-1951-2014, https://doi.org/10.5194/tc-8-1951-2014, 2014
R. Gladstone, M. Schäfer, T. Zwinger, Y. Gong, T. Strozzi, R. Mottram, F. Boberg, and J. C. Moore
The Cryosphere, 8, 1393–1405, https://doi.org/10.5194/tc-8-1393-2014, https://doi.org/10.5194/tc-8-1393-2014, 2014
B. Sun, J. C. Moore, T. Zwinger, L. Zhao, D. Steinhage, X. Tang, D. Zhang, X. Cui, and C. Martín
The Cryosphere, 8, 1121–1128, https://doi.org/10.5194/tc-8-1121-2014, https://doi.org/10.5194/tc-8-1121-2014, 2014
S. Cook, I. C. Rutt, T. Murray, A. Luckman, T. Zwinger, N. Selmes, A. Goldsack, and T. D. James
The Cryosphere, 8, 827–841, https://doi.org/10.5194/tc-8-827-2014, https://doi.org/10.5194/tc-8-827-2014, 2014
B. C. Gunter, O. Didova, R. E. M. Riva, S. R. M. Ligtenberg, J. T. M. Lenaerts, M. A. King, M. R. van den Broeke, and T. Urban
The Cryosphere, 8, 743–760, https://doi.org/10.5194/tc-8-743-2014, https://doi.org/10.5194/tc-8-743-2014, 2014
T. Zwinger, M. Schäfer, C. Martín, and J. C. Moore
The Cryosphere, 8, 607–621, https://doi.org/10.5194/tc-8-607-2014, https://doi.org/10.5194/tc-8-607-2014, 2014
T. Sato, T. Shiraiwa, R. Greve, H. Seddik, E. Edelmann, and T. Zwinger
Clim. Past, 10, 393–404, https://doi.org/10.5194/cp-10-393-2014, https://doi.org/10.5194/cp-10-393-2014, 2014
B. de Fleurian, O. Gagliardini, T. Zwinger, G. Durand, E. Le Meur, D. Mair, and P. Råback
The Cryosphere, 8, 137–153, https://doi.org/10.5194/tc-8-137-2014, https://doi.org/10.5194/tc-8-137-2014, 2014
J. A. Åström, T. I. Riikilä, T. Tallinen, T. Zwinger, D. Benn, J. C. Moore, and J. Timonen
The Cryosphere, 7, 1591–1602, https://doi.org/10.5194/tc-7-1591-2013, https://doi.org/10.5194/tc-7-1591-2013, 2013
O. Gagliardini, T. Zwinger, F. Gillet-Chaulet, G. Durand, L. Favier, B. de Fleurian, R. Greve, M. Malinen, C. Martín, P. Råback, J. Ruokolainen, M. Sacchettini, M. Schäfer, H. Seddik, and J. Thies
Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, https://doi.org/10.5194/gmd-6-1299-2013, 2013
A. S. Drouet, D. Docquier, G. Durand, R. Hindmarsh, F. Pattyn, O. Gagliardini, and T. Zwinger
The Cryosphere, 7, 395–406, https://doi.org/10.5194/tc-7-395-2013, https://doi.org/10.5194/tc-7-395-2013, 2013
L. Zhao, L. Tian, T. Zwinger, R. Ding, J. Zong, Q. Ye, and J. C. Moore
The Cryosphere Discuss., https://doi.org/10.5194/tcd-7-145-2013, https://doi.org/10.5194/tcd-7-145-2013, 2013
Revised manuscript not accepted
F. Gillet-Chaulet, O. Gagliardini, H. Seddik, M. Nodet, G. Durand, C. Ritz, T. Zwinger, R. Greve, and D. G. Vaughan
The Cryosphere, 6, 1561–1576, https://doi.org/10.5194/tc-6-1561-2012, https://doi.org/10.5194/tc-6-1561-2012, 2012
Related subject area
Solid Earth
REHEATFUNQ (REgional HEAT-Flow Uncertainty and aNomaly Quantification) 2.0.1: a model for regional aggregate heat flow distributions and anomaly quantification
A new temperature–photoperiod coupled phenology module in LPJ-GUESS model v4.1: optimizing estimation of terrestrial carbon and water processes
High-precision 1′ × 1′ bathymetric model of Philippine Sea inversed from marine gravity anomalies
Deciphering past earthquakes from the probabilistic modeling of paleoseismic records – the Paleoseismic EArthquake CHronologies code (PEACH, version 1)
Benchmarking the accuracy of higher order particle methods in geodynamic models of transient flow
Modelling detrital cosmogenic nuclide concentrations during landscape evolution in Cidre v2.0
IMEX_SfloW2D v2: a depth-averaged numerical flow model for volcanic gas–particle flows over complex topographies and water
Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle
AdaHRBF v1.0: gradient-adaptive Hermite–Birkhoff radial basis function interpolants for three-dimensional stratigraphic implicit modeling
PySubdiv 1.0: open-source geological modeling and reconstruction by non-manifold subdivision surfaces
Reconstructing tephra fall deposits via ensemble-based data assimilation techniques
ClinoformNet-1.0: stratigraphic forward modeling and deep learning for seismic clinoform delineation
Addressing challenges in uncertainty quantification: the case of geohazard assessments
DeepISMNet: three-dimensional implicit structural modeling with convolutional neural network
Towards automatic finite-element methods for geodynamics via Firedrake
MagmaFOAM-1.0: a modular framework for the simulation of magmatic systems
A global, spherical finite-element model for post-seismic deformation using Abaqus
SMAUG v1.0 – a user-friendly muon simulator for the imaging of geological objects in 3-D
CliffDelineaTool v1.2.0: an algorithm for identifying coastal cliff base and top positions
Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm
Structural, petrophysical, and geological constraints in potential field inversion using the Tomofast-x v1.0 open-source code
Spatial agents for geological surface modelling
RHEA v1.0: Enabling fully coupled simulations with hydro-geomechanical heterogeneity
Modelling of faults in LoopStructural 1.0
PALEOSTRIPv1.0 – a user-friendly 3D backtracking software to reconstruct paleo-bathymetries
LoopStructural 1.0: time-aware geological modelling
Sub3DNet1.0: a deep-learning model for regional-scale 3D subsurface structure mapping
Analytical solutions for mantle flow in cylindrical and spherical shells
Towards a model for structured mass movements: the OpenLISEM hazard model 2.0a
GO_3D_OBS: the multi-parameter benchmark geomodel for seismic imaging method assessment and next-generation 3D survey design (version 1.0)
PLUME-MoM-TSM 1.0.0: a volcanic column and umbrella cloud spreading model
HydrothermalFoam v1.0: a 3-D hydrothermal transport model for natural submarine hydrothermal systems
Synthetic seismicity distribution in Guerrero–Oaxaca subduction zone, Mexico, and its implications on the role of asperities in Gutenberg–Richter law
CobWeb 1.0: machine learning toolbox for tomographic imaging
pygeodyn 1.1.0: a Python package for geomagnetic data assimilation
IMEX_SfloW2D 1.0: a depth-averaged numerical flow model for pyroclastic avalanches
A multilayer approach and its application to model a local gravimetric quasi-geoid model over the North Sea: QGNSea V1.0
Development of an automatic delineation of cliff top and toe on very irregular planform coastlines (CliffMetrics v1.0)
Bayesian inference of earthquake rupture models using polynomial chaos expansion
Geodynamic diagnostics, scientific visualisation and StagLab 3.0
SaLEM (v1.0) – the Soil and Landscape Evolution Model (SaLEM) for simulation of regolith depth in periglacial environments
SILLi 1.0: a 1-D numerical tool quantifying the thermal effects of sill intrusions
The SPACE 1.0 model: a Landlab component for 2-D calculation of sediment transport, bedrock erosion, and landscape evolution
Ellipsoids (v1.0): 3-D magnetic modelling of ellipsoidal bodies
Global-scale modelling of melting and isotopic evolution of Earth's mantle: melting modules for TERRA
pynoddy 1.0: an experimental platform for automated 3-D kinematic and potential field modelling
Open-source modular solutions for flexural isostasy: gFlex v1.0
FPLUME-1.0: An integral volcanic plume model accounting for ash aggregation
PyXRD v0.6.7: a free and open-source program to quantify disordered phyllosilicates using multi-specimen X-ray diffraction profile fitting
r.randomwalk v1, a multi-functional conceptual tool for mass movement routing
Malte Jörn Ziebarth and Sebastian von Specht
Geosci. Model Dev., 17, 2783–2828, https://doi.org/10.5194/gmd-17-2783-2024, https://doi.org/10.5194/gmd-17-2783-2024, 2024
Short summary
Short summary
Thermal energy from Earth’s active interior constantly dissipates through Earth’s surface. This heat flow is not spatially uniform, and its exact pattern is hard to predict since it depends on crustal and mantle properties, both varying across scales. Our new model REHEATFUNQ addresses this difficulty by treating the fluctuations of heat flow within a region statistically. REHEATFUNQ estimates the regional distribution of heat flow and quantifies known structural signals therein.
Shouzhi Chen, Yongshuo H. Fu, Mingwei Li, Zitong Jia, Yishuo Cui, and Jing Tang
Geosci. Model Dev., 17, 2509–2523, https://doi.org/10.5194/gmd-17-2509-2024, https://doi.org/10.5194/gmd-17-2509-2024, 2024
Short summary
Short summary
It is still a challenge to achieve an accurate simulation of vegetation phenology in the dynamic global vegetation models (DGVMs). We implemented and coupled the spring and autumn phenology models into one of the DGVMs, LPJ-GUESS, and substantially improved the accuracy in capturing the start and end dates of growing seasons. Our study highlights the importance of getting accurate phenology estimations to reduce the uncertainties in plant distribution and terrestrial carbon and water cycling.
Dechao An, Jinyun Guo, Xiaotao Chang, Zhenming Wang, Yongjun Jia, Xin Liu, Valery Bondur, and Heping Sun
Geosci. Model Dev., 17, 2039–2052, https://doi.org/10.5194/gmd-17-2039-2024, https://doi.org/10.5194/gmd-17-2039-2024, 2024
Short summary
Short summary
Seafloor topography, as fundamental geoinformation in marine surveying and mapping, plays a crucial role in numerous scientific studies. In this paper, we focus on constructing a high-precision seafloor topography and bathymetry model for the Philippine Sea (5° N–35° N, 120° E–150° E), based on shipborne bathymetric data and marine gravity anomalies, and evaluate the reliability of the model's accuracy.
Octavi Gómez-Novell, Bruno Pace, Francesco Visini, Joanna Faure Walker, and Oona Scotti
Geosci. Model Dev., 16, 7339–7355, https://doi.org/10.5194/gmd-16-7339-2023, https://doi.org/10.5194/gmd-16-7339-2023, 2023
Short summary
Short summary
Knowing the rate at which earthquakes happen along active faults is crucial to characterize the hazard that they pose. We present an approach (Paleoseismic EArthquake CHronologies, PEACH) to correlate and compute seismic histories using paleoseismic data, a type of data that characterizes past seismic activity from the geological record. Our approach reduces the uncertainties of the seismic histories and overall can improve the knowledge on fault rupture behavior for the seismic hazard.
Rene Gassmöller, Juliane Dannberg, Wolfgang Bangerth, Elbridge Gerry Puckett, and Cedric Thieulot
EGUsphere, https://doi.org/10.5194/egusphere-2023-2765, https://doi.org/10.5194/egusphere-2023-2765, 2023
Short summary
Short summary
Numerical models that use simulated particles are a powerful tool for investigating flow in the interior of the Earth, but the accuracy of these models is not fully understood. We here present two new benchmarks, that allow to measure the model accuracy. We then document that better accuracy matters for applications like convection beneath an oceanic plate. Our benchmarks and methods are freely available to help the community develop better models.
Sébastien Carretier, Vincent Regard, Youssouf Abdelhafiz, and Bastien Plazolles
Geosci. Model Dev., 16, 6741–6755, https://doi.org/10.5194/gmd-16-6741-2023, https://doi.org/10.5194/gmd-16-6741-2023, 2023
Short summary
Short summary
We present the development of a code to simulate simultaneously the dynamics of landscapes over geological time and the evolution of the concentration of cosmogenic isotopes in grains throughout their transport from the slopes to the river outlets. This new model makes it possible to study the relationship between the detrital signal of cosmogenic isotope concentration measured in sediment and the erosion--deposition processes in watersheds.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, and Samantha Engwell
Geosci. Model Dev., 16, 6309–6336, https://doi.org/10.5194/gmd-16-6309-2023, https://doi.org/10.5194/gmd-16-6309-2023, 2023
Short summary
Short summary
We present version 2 of the numerical code IMEX-SfloW2D. With this version it is possible to simulate a wide range of volcanic mass flows (pyroclastic avalanches, lahars, pyroclastic surges), and here we present its application to transient dilute pyroclastic density currents (PDCs). A simulation of the 1883 Krakatau eruption demonstrates the capability of the numerical model to face a complex natural case involving the propagation of PDCs over the sea surface and across topographic obstacles.
Caroline J. van Calcar, Roderik S. W. van de Wal, Bas Blank, Bas de Boer, and Wouter van der Wal
Geosci. Model Dev., 16, 5473–5492, https://doi.org/10.5194/gmd-16-5473-2023, https://doi.org/10.5194/gmd-16-5473-2023, 2023
Short summary
Short summary
The waxing and waning of the Antarctic ice sheet caused the Earth’s surface to deform, which is stabilizing the ice sheet and mainly determined by the spatially variable viscosity of the mantle. Including this feedback in model simulations led to significant differences in ice sheet extent and ice thickness over the last glacial cycle. The results underline and quantify the importance of including this local feedback effect in ice sheet models when simulating the Antarctic ice sheet evolution.
Baoyi Zhang, Linze Du, Umair Khan, Yongqiang Tong, Lifang Wang, and Hao Deng
Geosci. Model Dev., 16, 3651–3674, https://doi.org/10.5194/gmd-16-3651-2023, https://doi.org/10.5194/gmd-16-3651-2023, 2023
Short summary
Short summary
We propose a Hermite–Birkhoff radial basis function (HRBF) formulation, AdaHRBF, with an adaptive gradient magnitude for continuous 3D stratigraphic potential field (SPF) modeling of multiple stratigraphic interfaces. In the linear system of HRBF interpolants constrained by the scattered on-contact attribute points and off-contact attitude points of a set of strata in 3D space, we add a novel optimization term to iteratively obtain the true gradient magnitude.
Mohammad Moulaeifard, Simon Bernard, and Florian Wellmann
Geosci. Model Dev., 16, 3565–3579, https://doi.org/10.5194/gmd-16-3565-2023, https://doi.org/10.5194/gmd-16-3565-2023, 2023
Short summary
Short summary
In this work, we propose a flexible framework to generate and interact with geological models using explicit surface representations. The essence of the work lies in the determination of the flexible control mesh, topologically similar to the main geological structure, watertight and controllable with few control points, to manage the geological structures. We exploited the subdivision surface method in our work, which is commonly used in the animation and gaming industry.
Leonardo Mingari, Antonio Costa, Giovanni Macedonio, and Arnau Folch
Geosci. Model Dev., 16, 3459–3478, https://doi.org/10.5194/gmd-16-3459-2023, https://doi.org/10.5194/gmd-16-3459-2023, 2023
Short summary
Short summary
Two novel techniques for ensemble-based data assimilation, suitable for semi-positive-definite variables with highly skewed uncertainty distributions such as tephra deposit mass loading, are applied to reconstruct the tephra fallout deposit resulting from the 2015 Calbuco eruption in Chile. The deposit spatial distribution and the ashfall volume according to the analyses are in good agreement with estimations based on field measurements and isopach maps reported in previous studies.
Hui Gao, Xinming Wu, Jinyu Zhang, Xiaoming Sun, and Zhengfa Bi
Geosci. Model Dev., 16, 2495–2513, https://doi.org/10.5194/gmd-16-2495-2023, https://doi.org/10.5194/gmd-16-2495-2023, 2023
Short summary
Short summary
We propose a workflow to automatically generate synthetic seismic data and corresponding stratigraphic labels (e.g., clinoform facies, relative geologic time, and synchronous horizons) by geological and geophysical forward modeling. Trained with only synthetic datasets, our network works well to accurately and efficiently predict clinoform facies in 2D and 3D field seismic data. Such a workflow can be easily extended for other geological and geophysical scenarios in the future.
Ibsen Chivata Cardenas, Terje Aven, and Roger Flage
Geosci. Model Dev., 16, 1601–1615, https://doi.org/10.5194/gmd-16-1601-2023, https://doi.org/10.5194/gmd-16-1601-2023, 2023
Short summary
Short summary
We discuss challenges in uncertainty quantification for geohazard assessments. The challenges arise from limited data and the one-off nature of geohazard features. The challenges include the credibility of predictions, input uncertainty, and assumptions’ impact. Considerations to increase credibility of the quantification are provided. Crucial tasks in the quantification are the exhaustive scrutiny of the background knowledge coupled with the assessment of deviations of assumptions made.
Zhengfa Bi, Xinming Wu, Zhaoliang Li, Dekuan Chang, and Xueshan Yong
Geosci. Model Dev., 15, 6841–6861, https://doi.org/10.5194/gmd-15-6841-2022, https://doi.org/10.5194/gmd-15-6841-2022, 2022
Short summary
Short summary
We present an implicit modeling method based on deep learning to produce a geologically valid and structurally compatible model from unevenly sampled structural data. Trained with automatically generated synthetic data with realistic features, our network can efficiently model geological structures without the need to solve large systems of mathematical equations, opening new opportunities for further leveraging deep learning to improve modeling capacity in many Earth science applications.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
Federico Brogi, Simone Colucci, Jacopo Matrone, Chiara Paola Montagna, Mattia De' Michieli Vitturi, and Paolo Papale
Geosci. Model Dev., 15, 3773–3796, https://doi.org/10.5194/gmd-15-3773-2022, https://doi.org/10.5194/gmd-15-3773-2022, 2022
Short summary
Short summary
Computer simulations play a fundamental role in understanding volcanic phenomena. The growing complexity of these simulations requires the development of flexible computational tools that can easily switch between sub-models and solution techniques as well as optimizations. MagmaFOAM is a newly developed library that allows for maximum flexibility for solving multiphase volcanic flows and promotes collaborative work for in-house and community model development, testing, and comparison.
Grace A. Nield, Matt A. King, Rebekka Steffen, and Bas Blank
Geosci. Model Dev., 15, 2489–2503, https://doi.org/10.5194/gmd-15-2489-2022, https://doi.org/10.5194/gmd-15-2489-2022, 2022
Short summary
Short summary
We present a finite-element model of post-seismic solid Earth deformation built in the software package Abaqus for the purpose of calculating post-seismic deformation in the far field of major earthquakes. The model is benchmarked against an existing open-source post-seismic model demonstrating good agreement. The advantage over existing models is the potential for simple modification to include 3-D Earth structure, non-linear rheologies and alternative or multiple sources of stress change.
Alessandro Lechmann, David Mair, Akitaka Ariga, Tomoko Ariga, Antonio Ereditato, Ryuichi Nishiyama, Ciro Pistillo, Paola Scampoli, Mykhailo Vladymyrov, and Fritz Schlunegger
Geosci. Model Dev., 15, 2441–2473, https://doi.org/10.5194/gmd-15-2441-2022, https://doi.org/10.5194/gmd-15-2441-2022, 2022
Short summary
Short summary
Muon tomography is a technology that is used often in geoscientific research. The know-how of data analysis is, however, still possessed by physicists who developed this technology. This article aims at providing geoscientists with the necessary tools to perform their own analyses. We hope that a lower threshold to enter the field of muon tomography will allow more geoscientists to engage with muon tomography. SMAUG is set up in a modular way to allow for its own modules to work in between.
Zuzanna M. Swirad and Adam P. Young
Geosci. Model Dev., 15, 1499–1512, https://doi.org/10.5194/gmd-15-1499-2022, https://doi.org/10.5194/gmd-15-1499-2022, 2022
Short summary
Short summary
Cliff base and top lines that delimit coastal cliff faces are usually manually digitized based on maps, aerial photographs, terrain models, etc. However, manual mapping is time consuming and depends on the mapper's decisions and skills. To increase the objectivity and efficiency of cliff mapping, we developed CliffDelineaTool, an algorithm that identifies cliff base and top positions along cross-shore transects using elevation and slope characteristics.
Holly Kyeore Han, Natalya Gomez, and Jeannette Xiu Wen Wan
Geosci. Model Dev., 15, 1355–1373, https://doi.org/10.5194/gmd-15-1355-2022, https://doi.org/10.5194/gmd-15-1355-2022, 2022
Short summary
Short summary
Interactions between ice sheets, sea level and the solid Earth occur over a range of timescales from years to tens of thousands of years. This requires coupled ice-sheet–sea-level models to exchange information frequently, leading to a quadratic increase in computation time with the number of model timesteps. We present a new sea-level model algorithm that allows coupled models to improve the computational feasibility and precisely capture short-term interactions within longer simulations.
Jérémie Giraud, Vitaliy Ogarko, Roland Martin, Mark Jessell, and Mark Lindsay
Geosci. Model Dev., 14, 6681–6709, https://doi.org/10.5194/gmd-14-6681-2021, https://doi.org/10.5194/gmd-14-6681-2021, 2021
Short summary
Short summary
We review different techniques to model the Earth's subsurface from geophysical data (gravity field anomaly, magnetic field anomaly) using geological models and measurements of the rocks' properties. We show examples of application using idealised examples reproducing realistic features and provide theoretical details of the open-source algorithm we use.
Eric A. de Kemp
Geosci. Model Dev., 14, 6661–6680, https://doi.org/10.5194/gmd-14-6661-2021, https://doi.org/10.5194/gmd-14-6661-2021, 2021
Short summary
Short summary
This is a proof of concept and review paper of spatial agents, with initial research focusing on geomodelling. The results may be of interest to others working on complex regional geological modelling with sparse data. Structural agent-based swarming behaviour is key to advancing this field. The study provides groundwork for research in structural geology 3D modelling with spatial agents. This work was done with NetLogo, a free agent modelling platform used mostly for teaching complex systems.
José M. Bastías Espejo, Andy Wilkins, Gabriel C. Rau, and Philipp Blum
Geosci. Model Dev., 14, 6257–6272, https://doi.org/10.5194/gmd-14-6257-2021, https://doi.org/10.5194/gmd-14-6257-2021, 2021
Short summary
Short summary
The hydraulic and mechanical properties of the subsurface are inherently heterogeneous. RHEA is a simulator that can perform couple hydro-geomechanical processes in heterogeneous porous media with steep gradients. RHEA is able to fully integrate spatial heterogeneity, allowing allocation of distributed hydraulic and geomechanical properties at mesh element level. RHEA is a valuable tool that can simulate problems considering realistic heterogeneity inherent to geologic formations.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, Guillaume Caumon, Mark Jessell, and Robin Armit
Geosci. Model Dev., 14, 6197–6213, https://doi.org/10.5194/gmd-14-6197-2021, https://doi.org/10.5194/gmd-14-6197-2021, 2021
Short summary
Short summary
Fault discontinuities in rock packages represent the plane where two blocks of rock have moved. They are challenging to incorporate into geological models because the geometry of the faulted rock units are defined by not only the location of the discontinuity but also the kinematics of the fault. In this paper, we outline a structural geology framework for incorporating faults into geological models by directly incorporating kinematics into the mathematical framework of the model.
Florence Colleoni, Laura De Santis, Enrico Pochini, Edy Forlin, Riccardo Geletti, Giuseppe Brancatelli, Magdala Tesauro, Martina Busetti, and Carla Braitenberg
Geosci. Model Dev., 14, 5285–5305, https://doi.org/10.5194/gmd-14-5285-2021, https://doi.org/10.5194/gmd-14-5285-2021, 2021
Short summary
Short summary
PALEOSTRIP has been developed in the framework of past Antarctic ice sheet reconstructions for periods when bathymetry around Antarctica differed substantially from today. It has been designed for users with no knowledge of numerical modelling and allows users to switch on and off the processes involved in backtracking and backstripping. Applications are broad, and it can be used to restore any continental margin bathymetry or sediment thickness and to perform basin analysis.
Lachlan Grose, Laurent Ailleres, Gautier Laurent, and Mark Jessell
Geosci. Model Dev., 14, 3915–3937, https://doi.org/10.5194/gmd-14-3915-2021, https://doi.org/10.5194/gmd-14-3915-2021, 2021
Short summary
Short summary
LoopStructural is an open-source 3D geological modelling library with a model design allowing for multiple different algorithms to be used for comparison for the same geology. Geological structures are modelled using structural geology concepts and techniques, allowing for complex structures such as overprinted folds and faults to be modelled. In the paper, we demonstrate automatically generating a 3-D model from map2loop-processed geological survey data of the Flinders Ranges, South Australia.
Zhenjiao Jiang, Dirk Mallants, Lei Gao, Tim Munday, Gregoire Mariethoz, and Luk Peeters
Geosci. Model Dev., 14, 3421–3435, https://doi.org/10.5194/gmd-14-3421-2021, https://doi.org/10.5194/gmd-14-3421-2021, 2021
Short summary
Short summary
Fast and reliable tools are required to extract hidden information from big geophysical and remote sensing data. A deep-learning model in 3D image construction from 2D image(s) is here developed for paleovalley mapping from globally available digital elevation data. The outstanding performance for 3D subsurface imaging gives confidence that this generic novel tool will make better use of existing geophysical and remote sensing data for improved management of limited earth resources.
Stephan C. Kramer, D. Rhodri Davies, and Cian R. Wilson
Geosci. Model Dev., 14, 1899–1919, https://doi.org/10.5194/gmd-14-1899-2021, https://doi.org/10.5194/gmd-14-1899-2021, 2021
Short summary
Short summary
Computational models of Earth's mantle require rigorous verification and validation. Analytical solutions of the underlying Stokes equations provide a method to verify that these equations are accurately solved for. However, their derivation in spherical and cylindrical shell domains with physically relevant boundary conditions is involved. This paper provides a number of solutions. They are provided in a Python package (Assess) and their use is demonstrated in a convergence study with Fluidity.
Bastian van den Bout, Theo van Asch, Wei Hu, Chenxiao X. Tang, Olga Mavrouli, Victor G. Jetten, and Cees J. van Westen
Geosci. Model Dev., 14, 1841–1864, https://doi.org/10.5194/gmd-14-1841-2021, https://doi.org/10.5194/gmd-14-1841-2021, 2021
Short summary
Short summary
Landslides, debris flows and other types of dense gravity-driven flows threaten livelihoods around the globe. Understanding the mechanics of these flows can be crucial for predicting their behaviour and reducing disaster risk. Numerical models assume that the solids and fluids of the flow are unstructured. The newly presented model captures the internal structure during movement. This important step can lead to more accurate predictions of landslide movement.
Andrzej Górszczyk and Stéphane Operto
Geosci. Model Dev., 14, 1773–1799, https://doi.org/10.5194/gmd-14-1773-2021, https://doi.org/10.5194/gmd-14-1773-2021, 2021
Short summary
Short summary
We present the 3D multi-parameter synthetic geomodel of the subduction zone, as well as the workflow designed to implement all of its components. The model contains different geological structures of various scales and complexities. It is intended to serve as a tool for the geophysical community to validate imaging approaches, design acquisition techniques, estimate uncertainties, benchmark computing approaches, etc.
Mattia de' Michieli Vitturi and Federica Pardini
Geosci. Model Dev., 14, 1345–1377, https://doi.org/10.5194/gmd-14-1345-2021, https://doi.org/10.5194/gmd-14-1345-2021, 2021
Short summary
Short summary
Here, we present PLUME-MoM-TSM, a volcanic plume model that allows us to quantify the formation of aggregates during the rise of the plume, model the phase change of water, and include the possibility to simulate the initial spreading of the tephra umbrella cloud intruding from the volcanic column into the atmosphere. The model is first applied to the 2015 Calbuco eruption (Chile) and provides an analytical relationship between the upwind spreading and some characteristic of the volcanic column.
Zhikui Guo, Lars Rüpke, and Chunhui Tao
Geosci. Model Dev., 13, 6547–6565, https://doi.org/10.5194/gmd-13-6547-2020, https://doi.org/10.5194/gmd-13-6547-2020, 2020
Short summary
Short summary
We present the 3-D hydro-thermo-transport model HydrothermalFoam v1.0, which we designed to provide the marine geosciences community with an easy-to-use and state-of-the-art tool for simulating mass and energy transport in submarine hydrothermal systems. HydrothermalFoam is based on the popular open-source platform OpenFOAM, comes with a number of tutorials, and is published under the GNU General Public License v3.0.
Marisol Monterrubio-Velasco, F. Ramón Zúñiga, Quetzalcoatl Rodríguez-Pérez, Otilio Rojas, Armando Aguilar-Meléndez, and Josep de la Puente
Geosci. Model Dev., 13, 6361–6381, https://doi.org/10.5194/gmd-13-6361-2020, https://doi.org/10.5194/gmd-13-6361-2020, 2020
Short summary
Short summary
The Mexican subduction zone along the Pacific coast is one of the most active seismic zones in the world, where every year larger-magnitude earthquakes shake huge inland cities such as Mexico City. In this work, we use TREMOL (sThochastic Rupture Earthquake ModeL) to simulate the seismicity observed in this zone. Our numerical results reinforce the hypothesis that in some subduction regions single asperities are responsible for producing the observed seismicity.
Swarup Chauhan, Kathleen Sell, Wolfram Rühaak, Thorsten Wille, and Ingo Sass
Geosci. Model Dev., 13, 315–334, https://doi.org/10.5194/gmd-13-315-2020, https://doi.org/10.5194/gmd-13-315-2020, 2020
Short summary
Short summary
We present CobWeb 1.0, a graphical user interface for analysing tomographic images of geomaterials. CobWeb offers different machine learning techniques for accurate multiphase image segmentation and visualizing material specific parameters such as pore size distribution, relative porosity and volume fraction. We demonstrate a novel approach of dual filtration and dual segmentation to eliminate edge enhancement artefact in synchrotron-tomographic datasets and provide the computational code.
Loïc Huder, Nicolas Gillet, and Franck Thollard
Geosci. Model Dev., 12, 3795–3803, https://doi.org/10.5194/gmd-12-3795-2019, https://doi.org/10.5194/gmd-12-3795-2019, 2019
Short summary
Short summary
The pygeodyn package is a geomagnetic data assimilation tool written in Python. It gives access to the Earth's core flow dynamics, controlled by geomagnetic observations, by means of a reduced numerical model anchored to geodynamo simulation statistics. It aims to provide the community with a user-friendly and tunable data assimilation algorithm. It can be used for education, geomagnetic model production or tests in conjunction with webgeodyn, a set of visualization tools for geomagnetic models.
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, Giacomo Lari, and Alvaro Aravena
Geosci. Model Dev., 12, 581–595, https://doi.org/10.5194/gmd-12-581-2019, https://doi.org/10.5194/gmd-12-581-2019, 2019
Short summary
Short summary
Pyroclastic avalanches are a type of granular flow generated at active volcanoes by different mechanisms, including the collapse of steep pyroclastic deposits (e.g., scoria and ash cones) and fountaining during moderately explosive eruptions. We present IMEX_SfloW2D, a depth-averaged flow model describing the granular mixture as a single-phase granular fluid. Benchmark cases and preliminary application to the simulation of the 11 February pyroclastic avalanche at Mt. Etna (Italy) are shown.
Yihao Wu, Zhicai Luo, Bo Zhong, and Chuang Xu
Geosci. Model Dev., 11, 4797–4815, https://doi.org/10.5194/gmd-11-4797-2018, https://doi.org/10.5194/gmd-11-4797-2018, 2018
Short summary
Short summary
A multilayer approach is parameterized for model development, and the multiple layers are located at different depths beneath the Earth’s surface. This method may be beneficial for gravity/manget field modeling, which may outperform the traditional single-layer approach.
Andres Payo, Bismarck Jigena Antelo, Martin Hurst, Monica Palaseanu-Lovejoy, Chris Williams, Gareth Jenkins, Kathryn Lee, David Favis-Mortlock, Andrew Barkwith, and Michael A. Ellis
Geosci. Model Dev., 11, 4317–4337, https://doi.org/10.5194/gmd-11-4317-2018, https://doi.org/10.5194/gmd-11-4317-2018, 2018
Short summary
Short summary
We describe a new algorithm that automatically delineates the cliff top and toe of a cliffed coastline from a digital elevation model (DEM). The algorithm builds upon existing methods but is specifically designed to resolve very irregular planform coastlines with many bays and capes, such as parts of the coastline of Great Britain.
Hugo Cruz-Jiménez, Guotu Li, Paul Martin Mai, Ibrahim Hoteit, and Omar M. Knio
Geosci. Model Dev., 11, 3071–3088, https://doi.org/10.5194/gmd-11-3071-2018, https://doi.org/10.5194/gmd-11-3071-2018, 2018
Short summary
Short summary
One of the most important challenges seismologists and earthquake engineers face is reliably estimating ground motion in an area prone to large damaging earthquakes. This study aimed at better understanding the relationship between characteristics of geological faults (e.g., hypocenter location, rupture size/location, etc.) and resulting ground motion, via statistical analysis of a rupture simulation model. This study provides important insight on ground-motion responses to geological faults.
Fabio Crameri
Geosci. Model Dev., 11, 2541–2562, https://doi.org/10.5194/gmd-11-2541-2018, https://doi.org/10.5194/gmd-11-2541-2018, 2018
Short summary
Short summary
Firstly, this study acts as a compilation of key geodynamic diagnostics and describes how to automatise them for a more efficient scientific procedure. Secondly, it outlines today's key pitfalls of scientific visualisation and provides means to circumvent them with, for example, a novel set of fully scientific colour maps. Thirdly, it introduces StagLab 3.0, a software that applies such fully automated diagnostics and state-of-the-art visualisation in the blink of an eye.
Michael Bock, Olaf Conrad, Andreas Günther, Ernst Gehrt, Rainer Baritz, and Jürgen Böhner
Geosci. Model Dev., 11, 1641–1652, https://doi.org/10.5194/gmd-11-1641-2018, https://doi.org/10.5194/gmd-11-1641-2018, 2018
Short summary
Short summary
We introduce the Soil and
Landscape Evolution Model (SaLEM) for the prediction of soil parent material evolution following a lithologically differentiated approach. The GIS tool is working within the software framework SAGA GIS. Weathering, erosion and transport functions are calibrated using extrinsic and intrinsic parameter data. First results indicate that our approach shows evidence for the spatiotemporal prediction of soil parental material properties.
Karthik Iyer, Henrik Svensen, and Daniel W. Schmid
Geosci. Model Dev., 11, 43–60, https://doi.org/10.5194/gmd-11-43-2018, https://doi.org/10.5194/gmd-11-43-2018, 2018
Short summary
Short summary
Igneous intrusions in sedimentary basins have a profound effect on the thermal structure of the hosting sedimentary rocks. In this paper, we present a user-friendly 1-D FEM-based tool, SILLi, that calculates the thermal effects of sill intrusions on the enclosing sedimentary stratigraphy. The motivation is to make a standardized numerical toolkit openly available that can be widely used by scientists with different backgrounds to test the effects of magmatic bodies in a wide variety of settings.
Charles M. Shobe, Gregory E. Tucker, and Katherine R. Barnhart
Geosci. Model Dev., 10, 4577–4604, https://doi.org/10.5194/gmd-10-4577-2017, https://doi.org/10.5194/gmd-10-4577-2017, 2017
Short summary
Short summary
Rivers control the movement of sediment and nutrients across Earth's surface. Understanding how rivers change through time is important for mitigating natural hazards and predicting Earth's response to climate change. We develop a new computer model for predicting how rivers cut through sediment and rock. Our model is designed to be joined with models of flooding, landslides, vegetation change, and other factors to provide a comprehensive toolbox for predicting changes to the landscape.
Diego Takahashi and Vanderlei C. Oliveira Jr.
Geosci. Model Dev., 10, 3591–3608, https://doi.org/10.5194/gmd-10-3591-2017, https://doi.org/10.5194/gmd-10-3591-2017, 2017
Short summary
Short summary
Ellipsoids are the only bodies for which the self-demagnetization can be treated analytically. This property is useful for modelling compact orebodies having high susceptibility. We present a review of the magnetic modelling of ellipsoids, propose a way of determining the isotropic susceptibility above which the self-demagnetization must be considered, and discuss the ambiguity between confocal ellipsoids, as well as provide a set of routines to model the magnetic field produced by ellipsoids.
Hein J. van Heck, J. Huw Davies, Tim Elliott, and Don Porcelli
Geosci. Model Dev., 9, 1399–1411, https://doi.org/10.5194/gmd-9-1399-2016, https://doi.org/10.5194/gmd-9-1399-2016, 2016
Short summary
Short summary
Currently, extensive geochemical databases of surface observations exist, but satisfying explanations of underlying mantle processes are lacking. We have implemented a new way to track both bulk compositions and concentrations of trace elements in a mantle convection code. In our model, chemical fractionation happens at evolving melting zones. We compare our results to a semi-analytical theory relating observed arrays of correlated Pb isotope compositions to melting age distributions.
J. Florian Wellmann, Sam T. Thiele, Mark D. Lindsay, and Mark W. Jessell
Geosci. Model Dev., 9, 1019–1035, https://doi.org/10.5194/gmd-9-1019-2016, https://doi.org/10.5194/gmd-9-1019-2016, 2016
Short summary
Short summary
We often obtain knowledge about the subsurface in the form of structural geological models, as a basis for subsurface usage or resource extraction. Here, we provide a modelling code to construct such models on the basis of significant deformational events in geological history, encapsulated in kinematic equations. Our methods simplify complex dynamic processes, but enable us to evaluate how events interact, and finally how certain we are about predictions of structures in the subsurface.
A. D. Wickert
Geosci. Model Dev., 9, 997–1017, https://doi.org/10.5194/gmd-9-997-2016, https://doi.org/10.5194/gmd-9-997-2016, 2016
Short summary
Short summary
Earth's lithosphere bends beneath surface loads, such as ice, sediments, and mountain belts. The pattern of this bending, or flexural isostatic response, is a function of both the loads and the spatially variable strength of the lithosphere. gFlex is an easy-to-use program to calculate flexural isostastic response, and may be used to better understand how ice sheets, glaciers, large lakes, sedimentary basins, volcanoes, and other surface loads interact with the solid Earth.
A. Folch, A. Costa, and G. Macedonio
Geosci. Model Dev., 9, 431–450, https://doi.org/10.5194/gmd-9-431-2016, https://doi.org/10.5194/gmd-9-431-2016, 2016
Short summary
Short summary
We present FPLUME-1.0, a steady-state 1-D cross-section-averaged eruption column model based on the buoyant plume theory (BPT). The model accounts for plume bending by wind, entrainment of ambient moisture, effects of water phase changes, particle fallout and re-entrainment, a new parameterization for the air entrainment coefficients and a model for wet aggregation of ash particles in presence of liquid water or ice.
M. Dumon and E. Van Ranst
Geosci. Model Dev., 9, 41–57, https://doi.org/10.5194/gmd-9-41-2016, https://doi.org/10.5194/gmd-9-41-2016, 2016
Short summary
Short summary
This paper presents a FOSS model called PyXRD used to improve the quantification of complex mixed-layer phyllosilicate assemblages using X-ray diffraction. The novelty of this model is the ab initio incorporation of the multi-specimen method, making it possible to share phases and their parameters across multiple specimens. We present results from a comparison of PyXRD with Sybilla v2.2.2 and a number of theoretical experiments illustrating the use of the multi-specimen set-up.
M. Mergili, J. Krenn, and H.-J. Chu
Geosci. Model Dev., 8, 4027–4043, https://doi.org/10.5194/gmd-8-4027-2015, https://doi.org/10.5194/gmd-8-4027-2015, 2015
Short summary
Short summary
r.randomwalk is a flexible and multi-functional open-source GIS tool for simulating the propagation of mass movements. Mass points are routed from given release pixels through a digital elevation model until a defined break criterion is reached. In contrast to existing tools, r.randomwalk includes functionalities to account for parameter uncertainties, and it offers built-in functions for validation and visualization. We show the key functionalities of r.randomwalk for three test areas.
Cited articles
Adhikari, S., Ivins, E. R., Larour, E., Seroussi, H., Morlighem, M., and Nowicki, S.: Future Antarctic bed topography and its implications for ice sheet dynamics, Solid Earth, 5, 569–584, https://doi.org/10.5194/se-5-569-2014, 2014. a
Australian National Computational Infrastructure: Online description of
raijin system, available at:
https://nci.org.au/systems-services/peak-system/raijin/ (last access:
19 July 2019), 2017. a
Barletta, V., Bevis, M., Smith, B., Wilson, T., Brown, A., Bordoni, A., Willis,
M., Khan, S., Rovira-Navarro, M., Dalziel, I., Smalley, R., Kendrick, E.,
Konfal, S., Caccamise, D., Aster, R., Nyblade, A., and Wiens, D.: Observed
rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability,
Science, 360, 1335–1339, https://doi.org/10.1126/science.aao1447, 2018. a
Bueler, E., Lingle, C., and Brown, J.: Fast computation of a viscoelastic
deformable Earth model for ice-sheet simulations, Ann. Glaciol., 46,
97–105, https://doi.org/10.3189/172756407782871567, 2007. a
Byckling, M., Kataja, J., Klemm, M., and Zwinger, T.: OpenMP SIMD
Vectorization and Threading of the Elmer Finite Element Software, in:
Proceedings 13th International Workshop on OpenMP, Springer
Lecture Notes, 123–137, https://doi.org/10.1007/978-3-319-65578-9_9, 2017. a, b
de Boer, B., Stocchi, P., and van de Wal, R. S. W.: A fully coupled 3-D ice-sheet–sea-level model: algorithm and applications, Geosci. Model Dev., 7, 2141–2156, https://doi.org/10.5194/gmd-7-2141-2014, 2014. a
Dziewonski, A. and Anderson, D.: Preliminary reference Earth model, Phys. Earth
Planet. In., 25, 297–356, https://doi.org/10.1016/0031-9201(81)90046-7, 1981. a
Eisenstat, S., Elman, H., and Schultz, M.: Variational iterative methods for
nonsymmetric systems of linear equations, SIAM J. Numer. Anal., 20, 345–357,
1983. a
Farrell, W. and Clark, J.: On Postglacial Sea Level, Geophys. J. Roy. Astr.
S., 46, 647–667, https://doi.org/10.1111/j.1365-246X.1976.tb01252.x, 1976. a
Gagliardini, O., Zwinger, T., Gillet-Chaulet, F., Durand, G., Favier, L., de Fleurian, B., Greve, R., Malinen, M., Martín, C., Råback, P., Ruokolainen, J., Sacchettini, M., Schäfer, M., Seddik, H., and Thies, J.: Capabilities and performance of Elmer/Ice, a new-generation ice sheet model, Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, 2013. a
Geuzaine, C. and Remacle, J.-F.: Gmsh: A 3-D finite element mesh generator
with built-in pre- and post-processing facilities, Int. J.
Numer. Meth. Eng., 11, 1309–1331, https://doi.org/10.1002/nme.2579,
2009. a
Gomez, N., Mitrovica, J., Huybers, P., and Clark, P.: Sea level as a
stabilizing factor for marine-ice-sheet grounding lines, Nat. Geosci., 3,
850–853, https://doi.org/10.1038/ngeo1012, 2010. a
Gomez, N., Pollard, D., and Mitrovica, J.: A 3-D coupled ice sheet – sea
level model applied to Antarctica through the last 40 ky, Earth
Planet. Sc. Lett., 384, 88–99, https://doi.org/10.1016/j.epsl.2013.09.042,
2013. a, b
Gomez, N., Pollard, D., and Holland, D.: Sea-level feedback lowers projections
of future Antarctic Ice-Sheet mass loss, Nat. Communi., 6, 8798,
https://doi.org/10.1038/ncomms9798, 2015. a
Greve, R.: Glacial Isostasy: Models for the Response of the Earth to Varying
Ice Loads, in: Continuum Mechanics and Applications in Geophysics and the
Environment Springer, Berlin, Germany, 307–325,
https://doi.org/10.1007/978-3-662-04439-1_16, 2001. a
Hibbitt, D., Karlsson, B., and Sorensen, P.: Getting Started with ABAQUS,
Version (6.14), hibbitt, Karlsson & Sorensen, Inc. Edn., 2016. a
Karypis, G. and Kumar, V.: Multilevel k-way Partitioning Scheme for Irregular
Graphs, J. Parallel Distr. Com., 48, 96–129,
https://doi.org/10.1006/jpdc.1997.1404, 1998. a
Larour, E., Seroussi, H., Adhikari, S., Ivins, E., Caron, L., Morlighem, M.,
and Schlegel, N.: Slowdown in Antarctic mass loss from solid Earth and
sea-level feedbacks, Science, 364, 7908, https://doi.org/10.1126/science.aav7908,
2019. a
Le Meur, E. and Huybrechts, P.: A comparison of different ways of dealing with
isostasy: examples from modeling the Antarctic ice sheet during the last
glacial cycle, Ann. Glaciol., 23, 309–317,
https://doi.org/10.3189/S0260305500013586, 1996. a, b, c
Nield, G., Barletta, V., Bordoni, A., King, M., Whitehouse, P., Clarke, P.,
Domack, E., Scambos, T., and Berthier, E.: Rapid bedrock uplift in the
Antarctic Peninsula explained by viscoelastic response to recent ice
unloading, Earth Planet. Sc. Lett., 397, 32–41,
https://doi.org/10.1016/j.epsl.2014.04.019, 2014. a
Peltier, W.: The impulse response of a Maxwell earth, Rev. Geophys. Space
Ge., 12, 649–669, https://doi.org/10.1029/RG012i004p00649, 1974. a, b
Råback, P., Malinen, M., Ruokolainen, J., Pursula, A., and Zwinger, T.: Elmer
Models Manual,
available at: http://www.nic.funet.fi/pub/sci/physics/elmer/doc/ElmerModelsManual.pdf (last access: 5 March 2020),
2019. a
Rutt, I., Hagdorn, M., Hulton, N., and Payne, A.: The Glimmer community ice
sheet model, J. Geophys. Res., 114, F02004, https://doi.org/10.1029/2008JF001015,
2009. a
Schoof, C.: Ice Sheet Grounding Line Dynamics: Steady States, Stability, and
Hysteresis, J. Geophys. Res., 112, F03S28,
https://doi.org/10.1029/2006JF000664, 2007. a
Spada, G.: The theory behind TABOO, available at: https://github.com/danielemelini/TABOO/blob/master/DOC/TABOO-theory.pdf (last access: 5 March 2020), 2003. a
Spada, G., Antonioli, A., Boschi, L., Brandi, V., Cianetti, S., Galvani, G.,
Giunchi, C., Perniola, B., Piana Agostinetti, N., Piersanti, A., and Stocchi,
P.: TABOO, User Guide, available at: https://github.com/danielemelini/TABOO/blob/master/DOC/TABOO_User_Guide.pdf (last access: 5 March 2020), 2003. a
Steffen, H., Kaufmann, G., and Wu, P.: Three-dimensional finite-element
modeling of the glacial isostatic adjustment in Fennoscandia, Earth Planet.
Sc. Lett., 250, 358–375, https://doi.org/10.1016/j.epsl.2006.08.003, 2006. a, b
Whitehouse, P. L.: Glacial isostatic adjustment modelling: historical perspectives, recent advances, and future directions, Earth Surf. Dynam., 6, 401–429, https://doi.org/10.5194/esurf-6-401-2018, 2018. a
Wu, P. and Johnston, P.: Validity of Using Flat-Earth Finite Element Models in
the Study of Postglacial Rebound, in: Dynamics of the Ice Age Earth,
Trans Tech Publications Ltd, Switzerland, 191–202, 1998. a
Zwinger, T.: Deformation in flat-earth model Elmer/Earth, https://doi.org/10.5446/44086, 2019. a
Short summary
We present a newly developed flat-earth model, Elmer/Earth, for viscoelastic treatment of solid earth deformation under ice loads. Unlike many previous approaches with proprietary software, this model is based on the open-source FEM code Elmer, with the advantage for scientists to apply and alter the model without license constraints. The new-generation full-stress ice-sheet model Elmer/Ice shares the same code base, enabling future coupled ice-sheet–glacial-isostatic-adjustment simulations.
We present a newly developed flat-earth model, Elmer/Earth, for viscoelastic treatment of solid...
