Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG
The thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheets may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs...
Main Authors: | , , , , |
---|---|
Format: | Text |
Language: | English |
Published: |
2020
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-2020-78 https://gmd.copernicus.org/preprints/gmd-2020-78/ |
id |
ftcopernicus:oai:publications.copernicus.org:gmdd84602 |
---|---|
record_format |
openpolar |
spelling |
ftcopernicus:oai:publications.copernicus.org:gmdd84602 2023-05-15T16:40:00+02:00 Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG Rückamp, Martin Humbert, Angelika Kleiner, Thomas Morlighem, Mathieu Seroussi, Helene 2020-04-24 application/pdf https://doi.org/10.5194/gmd-2020-78 https://gmd.copernicus.org/preprints/gmd-2020-78/ eng eng doi:10.5194/gmd-2020-78 https://gmd.copernicus.org/preprints/gmd-2020-78/ eISSN: 1991-9603 Text 2020 ftcopernicus https://doi.org/10.5194/gmd-2020-78 2020-07-20T16:22:13Z The thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheets may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs). Additionally, the highly anisotropic geometry of the 3D elements in ice sheet modelling poses a problem for stabilization approaches in advection dominated problems. Here, we present extended enthalpy formulations within the finite-element Ice Sheet System Model (ISSM) that show a better performance to earlier implementations. In a first polythermal-slab experiment, we found that the treatment of the discontinuous conductivities at the CTS with a geometric mean produce more accurate results compared to the arithmetic or harmonic mean. This improvement is particularly efficient when applied to coarse vertical resolutions. In a second ice dome experiment, we find that the numerical solution is sensitive to the choice of stabilization parameters in the well-established Streamline Upwind Petrov–Galerkin (SUPG) method. As standard literature values for the SUPG stabilization parameter are not accounting for the highly anisotropic geometry of the 3D elements in ice sheet modelling, we propose a novel Anisotropic SUPG (ASUPG) formulation. This formulation circumvents the problem of high aspect-ratio by treating the horizontal and vertical directions separately in the stabilization coefficients. The ASUPG method provides accurate results for the thermodynamic equation on geometries with very small aspect ratios like ice sheets. Text Ice Sheet Copernicus Publications: E-Journals |
institution |
Open Polar |
collection |
Copernicus Publications: E-Journals |
op_collection_id |
ftcopernicus |
language |
English |
description |
The thermal state of an ice sheet is an important control on its past and future evolution. Some parts of the ice sheets may be polythermal, leading to discontinuous properties at the cold–temperate transition surface (CTS). These discontinuities require a careful treatment in ice sheet models (ISMs). Additionally, the highly anisotropic geometry of the 3D elements in ice sheet modelling poses a problem for stabilization approaches in advection dominated problems. Here, we present extended enthalpy formulations within the finite-element Ice Sheet System Model (ISSM) that show a better performance to earlier implementations. In a first polythermal-slab experiment, we found that the treatment of the discontinuous conductivities at the CTS with a geometric mean produce more accurate results compared to the arithmetic or harmonic mean. This improvement is particularly efficient when applied to coarse vertical resolutions. In a second ice dome experiment, we find that the numerical solution is sensitive to the choice of stabilization parameters in the well-established Streamline Upwind Petrov–Galerkin (SUPG) method. As standard literature values for the SUPG stabilization parameter are not accounting for the highly anisotropic geometry of the 3D elements in ice sheet modelling, we propose a novel Anisotropic SUPG (ASUPG) formulation. This formulation circumvents the problem of high aspect-ratio by treating the horizontal and vertical directions separately in the stabilization coefficients. The ASUPG method provides accurate results for the thermodynamic equation on geometries with very small aspect ratios like ice sheets. |
format |
Text |
author |
Rückamp, Martin Humbert, Angelika Kleiner, Thomas Morlighem, Mathieu Seroussi, Helene |
spellingShingle |
Rückamp, Martin Humbert, Angelika Kleiner, Thomas Morlighem, Mathieu Seroussi, Helene Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
author_facet |
Rückamp, Martin Humbert, Angelika Kleiner, Thomas Morlighem, Mathieu Seroussi, Helene |
author_sort |
Rückamp, Martin |
title |
Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
title_short |
Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
title_full |
Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
title_fullStr |
Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
title_full_unstemmed |
Extended enthalpy formulations in the ice flow model ISSM version 4.17: discontinuous conductivity and anisotropic SUPG |
title_sort |
extended enthalpy formulations in the ice flow model issm version 4.17: discontinuous conductivity and anisotropic supg |
publishDate |
2020 |
url |
https://doi.org/10.5194/gmd-2020-78 https://gmd.copernicus.org/preprints/gmd-2020-78/ |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_source |
eISSN: 1991-9603 |
op_relation |
doi:10.5194/gmd-2020-78 https://gmd.copernicus.org/preprints/gmd-2020-78/ |
op_doi |
https://doi.org/10.5194/gmd-2020-78 |
_version_ |
1766030358392864768 |