Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14)
<p>Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution in ice sheet models to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of...
Main Authors: | , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Copernicus Publications
2019-01-01
|
Series: | Geoscientific Model Development |
Online Access: | https://www.geosci-model-dev.net/12/215/2019/gmd-12-215-2019.pdf |
id |
doaj-f5751877efa4409792959422d4637ec2 |
---|---|
record_format |
Article |
spelling |
doaj-f5751877efa4409792959422d4637ec22020-11-24T21:50:09ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032019-01-011221523210.5194/gmd-12-215-2019Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14)T. D. dos Santos0M. Morlighem1H. Seroussi2P. R. B. Devloo3J. C. Simões4Department of Structures, School of Civil Engineering, Architecture and Urban Design, University of Campinas (UNICAMP), Campinas, SP, BrazilDepartment of Earth System Science, University of California, Irvine, CA, USAJet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USADepartment of Structures, School of Civil Engineering, Architecture and Urban Design, University of Campinas (UNICAMP), Campinas, SP, BrazilPolar and Climate Center, Geosciences Institute, Federal University of Rio Grande do Sul (UFRGS), Porto Alegre, RS, Brazil<p>Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution in ice sheet models to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of the grounding line indeed plays a major role in ice sheet dynamics, as it is a fundamental control on marine ice sheet stability. Numerical modeling of a grounding line requires significant computational resources since the accuracy of its position depends on grid or mesh resolution. A technique that improves accuracy with reduced computational cost is the adaptive mesh refinement (AMR) approach. We present here the implementation of the AMR technique in the finite element Ice Sheet System Model (ISSM) to simulate grounding line dynamics under two different benchmarks: MISMIP3d and MISMIP<span class="inline-formula">+</span>. We test different refinement criteria: (a) distance around the grounding line, (b) a posteriori error estimator, the Zienkiewicz–Zhu (ZZ) error estimator, and (c) different combinations of (a) and (b). In both benchmarks, the ZZ error estimator presents high values around the grounding line. In the MISMIP<span class="inline-formula">+</span> setup, this estimator also presents high values in the grounded part of the ice sheet, following the complex shape of the bedrock geometry. The ZZ estimator helps guide the refinement procedure such that AMR performance is improved. Our results show that computational time with AMR depends on the required accuracy, but in all cases, it is significantly shorter than for uniformly refined meshes. We conclude that AMR without an associated error estimator should be avoided, especially for real glaciers that have a complex bed geometry.</p>https://www.geosci-model-dev.net/12/215/2019/gmd-12-215-2019.pdf |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
T. D. dos Santos M. Morlighem H. Seroussi P. R. B. Devloo J. C. Simões |
spellingShingle |
T. D. dos Santos M. Morlighem H. Seroussi P. R. B. Devloo J. C. Simões Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) Geoscientific Model Development |
author_facet |
T. D. dos Santos M. Morlighem H. Seroussi P. R. B. Devloo J. C. Simões |
author_sort |
T. D. dos Santos |
title |
Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) |
title_short |
Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) |
title_full |
Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) |
title_fullStr |
Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) |
title_full_unstemmed |
Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14) |
title_sort |
implementation and performance of adaptive mesh refinement in the ice sheet system model (issm v4.14) |
publisher |
Copernicus Publications |
series |
Geoscientific Model Development |
issn |
1991-959X 1991-9603 |
publishDate |
2019-01-01 |
description |
<p>Accurate projections of the evolution of ice sheets in a changing climate
require a fine mesh/grid resolution in ice sheet models to correctly capture
fundamental physical processes, such as the evolution of the grounding line,
the region where grounded ice starts to float. The evolution of the grounding
line indeed plays a major role in ice sheet dynamics, as it is a fundamental
control on marine ice sheet stability. Numerical modeling of a grounding line
requires significant computational resources since the accuracy of its
position depends on grid or mesh resolution. A technique that improves
accuracy with reduced computational cost is the adaptive mesh refinement
(AMR) approach. We present here the implementation of the AMR technique in
the finite element Ice Sheet System Model (ISSM) to simulate grounding line
dynamics under two different benchmarks: MISMIP3d and MISMIP<span class="inline-formula">+</span>. We test
different refinement criteria: (a) distance around the grounding line, (b) a
posteriori error estimator, the Zienkiewicz–Zhu (ZZ) error estimator, and
(c) different combinations of (a) and (b). In both benchmarks, the ZZ error
estimator presents high values around the grounding line. In the MISMIP<span class="inline-formula">+</span> setup,
this estimator also presents high values in the grounded
part of the ice sheet, following the complex shape of the bedrock geometry.
The ZZ estimator helps guide the refinement procedure such that AMR
performance is improved. Our results show that computational time with AMR
depends on the required accuracy, but in all cases, it is significantly
shorter than for uniformly refined meshes. We conclude that AMR without an
associated error estimator should be avoided, especially for real glaciers
that have a complex bed geometry.</p> |
url |
https://www.geosci-model-dev.net/12/215/2019/gmd-12-215-2019.pdf |
work_keys_str_mv |
AT tddossantos implementationandperformanceofadaptivemeshrefinementintheicesheetsystemmodelissmv414 AT mmorlighem implementationandperformanceofadaptivemeshrefinementintheicesheetsystemmodelissmv414 AT hseroussi implementationandperformanceofadaptivemeshrefinementintheicesheetsystemmodelissmv414 AT prbdevloo implementationandperformanceofadaptivemeshrefinementintheicesheetsystemmodelissmv414 AT jcsimoes implementationandperformanceofadaptivemeshrefinementintheicesheetsystemmodelissmv414 |
_version_ |
1725885023786106880 |