Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18
Time-dependent simulations of ice sheets require two equations to be solved: the mass transport equation, derived from the conservation of mass, and the stress balance equation, derived from the conservation of momentum. The mass transport equation controls the advection of ice from the interior of...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , |
Format: | Article in Journal/Newspaper |
Language: | English |
Published: |
Copernicus Publications
2021
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-14-2545-2021 https://doaj.org/article/b06ba2cc09c24d29805bdd2b41c2cd08 |
id |
ftdoajarticles:oai:doaj.org/article:b06ba2cc09c24d29805bdd2b41c2cd08 |
---|---|
record_format |
openpolar |
spelling |
ftdoajarticles:oai:doaj.org/article:b06ba2cc09c24d29805bdd2b41c2cd08 2023-05-15T13:24:17+02:00 Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 T. D. dos Santos M. Morlighem H. Seroussi 2021-05-01T00:00:00Z https://doi.org/10.5194/gmd-14-2545-2021 https://doaj.org/article/b06ba2cc09c24d29805bdd2b41c2cd08 EN eng Copernicus Publications https://gmd.copernicus.org/articles/14/2545/2021/gmd-14-2545-2021.pdf https://doaj.org/toc/1991-959X https://doaj.org/toc/1991-9603 doi:10.5194/gmd-14-2545-2021 1991-959X 1991-9603 https://doaj.org/article/b06ba2cc09c24d29805bdd2b41c2cd08 Geoscientific Model Development, Vol 14, Pp 2545-2573 (2021) Geology QE1-996.5 article 2021 ftdoajarticles https://doi.org/10.5194/gmd-14-2545-2021 2022-12-31T06:40:45Z Time-dependent simulations of ice sheets require two equations to be solved: the mass transport equation, derived from the conservation of mass, and the stress balance equation, derived from the conservation of momentum. The mass transport equation controls the advection of ice from the interior of the ice sheet towards its periphery, thereby changing its geometry. Because it is based on an advection equation, a stabilization scheme needs to be employed when solved using the finite-element method. Several stabilization schemes exist in the finite-element method framework, but their respective accuracy and robustness have not yet been systematically assessed for glaciological applications. Here, we compare classical schemes used in the context of the finite-element method: (i) artificial diffusion, (ii) streamline upwinding, (iii) streamline upwind Petrov–Galerkin, (iv) discontinuous Galerkin, and (v) flux-corrected transport. We also look at the stress balance equation, which is responsible for computing the ice velocity that “advects” the ice downstream. To improve the velocity computation accuracy, the ice-sheet modeling community employs several sub-element parameterizations of physical processes at the grounding line, the point where the grounded ice starts to float onto the ocean. Here, we introduce a new sub-element parameterization for the driving stress, the force that drives the ice-sheet flow. We analyze the response of each stabilization scheme by running transient simulations forced by ice-shelf basal melt. The simulations are based on an idealized ice-sheet geometry for which there is no influence of bedrock topography. We also perform transient simulations of the Amundsen Sea Embayment, West Antarctica, where real bedrock and surface elevations are employed. In both idealized and real ice-sheet experiments, stabilization schemes based on artificial diffusion lead systematically to a bias towards more mass loss in comparison to the other schemes and therefore should be avoided or employed with a ... Article in Journal/Newspaper Amundsen Sea Antarc* Antarctica Ice Sheet Ice Shelf West Antarctica Directory of Open Access Journals: DOAJ Articles Amundsen Sea West Antarctica Geoscientific Model Development 14 5 2545 2573 |
institution |
Open Polar |
collection |
Directory of Open Access Journals: DOAJ Articles |
op_collection_id |
ftdoajarticles |
language |
English |
topic |
Geology QE1-996.5 |
spellingShingle |
Geology QE1-996.5 T. D. dos Santos M. Morlighem H. Seroussi Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
topic_facet |
Geology QE1-996.5 |
description |
Time-dependent simulations of ice sheets require two equations to be solved: the mass transport equation, derived from the conservation of mass, and the stress balance equation, derived from the conservation of momentum. The mass transport equation controls the advection of ice from the interior of the ice sheet towards its periphery, thereby changing its geometry. Because it is based on an advection equation, a stabilization scheme needs to be employed when solved using the finite-element method. Several stabilization schemes exist in the finite-element method framework, but their respective accuracy and robustness have not yet been systematically assessed for glaciological applications. Here, we compare classical schemes used in the context of the finite-element method: (i) artificial diffusion, (ii) streamline upwinding, (iii) streamline upwind Petrov–Galerkin, (iv) discontinuous Galerkin, and (v) flux-corrected transport. We also look at the stress balance equation, which is responsible for computing the ice velocity that “advects” the ice downstream. To improve the velocity computation accuracy, the ice-sheet modeling community employs several sub-element parameterizations of physical processes at the grounding line, the point where the grounded ice starts to float onto the ocean. Here, we introduce a new sub-element parameterization for the driving stress, the force that drives the ice-sheet flow. We analyze the response of each stabilization scheme by running transient simulations forced by ice-shelf basal melt. The simulations are based on an idealized ice-sheet geometry for which there is no influence of bedrock topography. We also perform transient simulations of the Amundsen Sea Embayment, West Antarctica, where real bedrock and surface elevations are employed. In both idealized and real ice-sheet experiments, stabilization schemes based on artificial diffusion lead systematically to a bias towards more mass loss in comparison to the other schemes and therefore should be avoided or employed with a ... |
format |
Article in Journal/Newspaper |
author |
T. D. dos Santos M. Morlighem H. Seroussi |
author_facet |
T. D. dos Santos M. Morlighem H. Seroussi |
author_sort |
T. D. dos Santos |
title |
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
title_short |
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
title_full |
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
title_fullStr |
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
title_full_unstemmed |
Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18 |
title_sort |
assessment of numerical schemes for transient, finite-element ice flow models using issm v4.18 |
publisher |
Copernicus Publications |
publishDate |
2021 |
url |
https://doi.org/10.5194/gmd-14-2545-2021 https://doaj.org/article/b06ba2cc09c24d29805bdd2b41c2cd08 |
geographic |
Amundsen Sea West Antarctica |
geographic_facet |
Amundsen Sea West Antarctica |
genre |
Amundsen Sea Antarc* Antarctica Ice Sheet Ice Shelf West Antarctica |
genre_facet |
Amundsen Sea Antarc* Antarctica Ice Sheet Ice Shelf West Antarctica |
op_source |
Geoscientific Model Development, Vol 14, Pp 2545-2573 (2021) |
op_relation |
https://gmd.copernicus.org/articles/14/2545/2021/gmd-14-2545-2021.pdf https://doaj.org/toc/1991-959X https://doaj.org/toc/1991-9603 doi:10.5194/gmd-14-2545-2021 1991-959X 1991-9603 https://doaj.org/article/b06ba2cc09c24d29805bdd2b41c2cd08 |
op_doi |
https://doi.org/10.5194/gmd-14-2545-2021 |
container_title |
Geoscientific Model Development |
container_volume |
14 |
container_issue |
5 |
container_start_page |
2545 |
op_container_end_page |
2573 |
_version_ |
1766378495837995008 |