Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11)
Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using advanced numerical solvers for transient ice-flow simulations conducted with the JP...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , , , |
Format: | Text |
Language: | English |
Published: |
2018
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-10-155-2017 https://gmd.copernicus.org/articles/10/155/2017/ |
id |
ftcopernicus:oai:publications.copernicus.org:gmd51207 |
---|---|
record_format |
openpolar |
spelling |
ftcopernicus:oai:publications.copernicus.org:gmd51207 2023-05-15T16:40:13+02:00 Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) Habbal, Feras Larour, Eric Morlighem, Mathieu Seroussi, Helene Borstad, Christopher P. Rignot, Eric 2018-09-27 application/pdf https://doi.org/10.5194/gmd-10-155-2017 https://gmd.copernicus.org/articles/10/155/2017/ eng eng doi:10.5194/gmd-10-155-2017 https://gmd.copernicus.org/articles/10/155/2017/ eISSN: 1991-9603 Text 2018 ftcopernicus https://doi.org/10.5194/gmd-10-155-2017 2020-07-20T16:23:52Z Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using advanced numerical solvers for transient ice-flow simulations conducted with the JPL–UCI Ice Sheet System Model (ISSM). We identify optimal numerical solvers by testing a broad suite of readily available solvers, ranging from direct sparse solvers to preconditioned iterative methods, on the commonly used Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models benchmark tests. Three types of analyses are considered: mass transport, horizontal stress balance, and incompressibility. The results of the fastest solvers for each analysis type are ranked based on their scalability across mesh size and basal boundary conditions. We find that the fastest iterative solvers are ∼ 1.5–100 times faster than the default direct solver used in ISSM, with speed-ups improving rapidly with increased mesh resolution. We provide a set of recommendations for users in search of efficient solvers to use for transient ice-flow simulations, enabling higher-resolution meshes and faster turnaround time. The end result will be improved transient simulations for short-term, highly resolved forward projections (10–100 year time scale) and also improved long-term paleo-reconstructions using higher-order representations of stresses in the ice. This analysis will also enable a new generation of comprehensive uncertainty quantification assessments of forward sea-level rise projections, which rely heavily on ensemble or sampling approaches that are inherently expensive. Text Ice Sheet Copernicus Publications: E-Journals Geoscientific Model Development 10 1 155 168 |
institution |
Open Polar |
collection |
Copernicus Publications: E-Journals |
op_collection_id |
ftcopernicus |
language |
English |
description |
Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using advanced numerical solvers for transient ice-flow simulations conducted with the JPL–UCI Ice Sheet System Model (ISSM). We identify optimal numerical solvers by testing a broad suite of readily available solvers, ranging from direct sparse solvers to preconditioned iterative methods, on the commonly used Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models benchmark tests. Three types of analyses are considered: mass transport, horizontal stress balance, and incompressibility. The results of the fastest solvers for each analysis type are ranked based on their scalability across mesh size and basal boundary conditions. We find that the fastest iterative solvers are ∼ 1.5–100 times faster than the default direct solver used in ISSM, with speed-ups improving rapidly with increased mesh resolution. We provide a set of recommendations for users in search of efficient solvers to use for transient ice-flow simulations, enabling higher-resolution meshes and faster turnaround time. The end result will be improved transient simulations for short-term, highly resolved forward projections (10–100 year time scale) and also improved long-term paleo-reconstructions using higher-order representations of stresses in the ice. This analysis will also enable a new generation of comprehensive uncertainty quantification assessments of forward sea-level rise projections, which rely heavily on ensemble or sampling approaches that are inherently expensive. |
format |
Text |
author |
Habbal, Feras Larour, Eric Morlighem, Mathieu Seroussi, Helene Borstad, Christopher P. Rignot, Eric |
spellingShingle |
Habbal, Feras Larour, Eric Morlighem, Mathieu Seroussi, Helene Borstad, Christopher P. Rignot, Eric Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
author_facet |
Habbal, Feras Larour, Eric Morlighem, Mathieu Seroussi, Helene Borstad, Christopher P. Rignot, Eric |
author_sort |
Habbal, Feras |
title |
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
title_short |
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
title_full |
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
title_fullStr |
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
title_full_unstemmed |
Optimal numerical solvers for transient simulations of ice flow using the Ice Sheet System Model (ISSM versions 4.2.5 and 4.11) |
title_sort |
optimal numerical solvers for transient simulations of ice flow using the ice sheet system model (issm versions 4.2.5 and 4.11) |
publishDate |
2018 |
url |
https://doi.org/10.5194/gmd-10-155-2017 https://gmd.copernicus.org/articles/10/155/2017/ |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_source |
eISSN: 1991-9603 |
op_relation |
doi:10.5194/gmd-10-155-2017 https://gmd.copernicus.org/articles/10/155/2017/ |
op_doi |
https://doi.org/10.5194/gmd-10-155-2017 |
container_title |
Geoscientific Model Development |
container_volume |
10 |
container_issue |
1 |
container_start_page |
155 |
op_container_end_page |
168 |
_version_ |
1766030596891475968 |