English
 
Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT
 
 
DownloadE-Mail
  A comparison of the stability and performance of depth-integrated ice-dynamics solvers

Robinson, A., Goldberg, D., Lipscomb, W. H. (2022): A comparison of the stability and performance of depth-integrated ice-dynamics solvers. - The Cryosphere, 16, 2, 689-709.
https://doi.org/10.5194/tc-16-689-2022

Item is

Files

show Files
hide Files
:
28102oa.pdf (Publisher version), 3MB
Name:
28102oa.pdf
Description:
-
Visibility:
Public
MIME-Type / Checksum:
application/pdf / [MD5]
Technical Metadata:
Copyright Date:
-
Copyright Info:
-

Locators

show

Creators

show
hide
 Creators:
Robinson, Alexander1, Author              
Goldberg, Daniel2, Author
Lipscomb, William H.2, Author
Affiliations:
1Potsdam Institute for Climate Impact Research, Potsdam, ou_persistent13              
2External Organizations, ou_persistent22              

Content

show
hide
Free keywords: -
 Abstract: In the last decade, the number of ice-sheet models has increased substantially, in line with the growth of the glaciological community. These models use solvers based on different approximations of ice dynamics. In particular, several depth-integrated dynamics solvers have emerged as fast solvers capable of resolving the relevant physics of ice sheets at the continental scale. However, the numerical stability of these schemes has not been studied systematically to evaluate their effectiveness in practice. Here we focus on three such solvers, the so-called Hybrid, L1L2-SIA and DIVA solvers, as well as the well-known SIA and SSA solvers as boundary cases. We investigate the numerical stability of these solvers as a function of grid resolution and the state of the ice sheet for an explicit time discretization scheme of the mass conservation step. Under simplified conditions with constant viscosity, the maximum stable time step of the Hybrid solver, like the SIA solver, has a quadratic dependence on grid resolution. In contrast, the DIVA solver has a maximum time step that is independent of resolution as the grid becomes increasingly refined, like the SSA solver. A simple 1D implementation of the L1L2-SIA solver indicates that it should behave similarly, but in practice, the complexity of its implementation appears to restrict its stability. In realistic simulations of the Greenland Ice Sheet with a nonlinear rheology, the DIVA and SSA solvers maintain superior numerical stability, while the SIA, Hybrid and L1L2-SIA solvers show markedly poorer performance. At a grid resolution of Δx=4 km, the DIVA solver runs approximately 20 times faster than the Hybrid and L1L2-SIA solvers as a result of a larger stable time step. Our analysis shows that as resolution increases, the ice-dynamics solver can act as a bottleneck to model performance. The DIVA solver emerges as a clear outlier in terms of both model performance and its representation of the ice-flow physics itself.

Details

show
hide
Language(s): eng - English
 Dates: 2022-01-122022-02-252022-02-25
 Publication Status: Finally published
 Pages: 21
 Publishing info: -
 Table of Contents: -
 Rev. Type: Peer
 Identifiers: DOI: 10.5194/tc-16-689-2022
PIKDOMAIN: RD1 - Earth System Analysis
Organisational keyword: RD1 - Earth System Analysis
MDB-ID: Entry suspended
OATYPE: Gold Open Access
 Degree: -

Event

show

Legal Case

show

Project information

show

Source 1

show
hide
Title: The Cryosphere
Source Genre: Journal, SCI, Scopus, p3, oa
 Creator(s):
Affiliations:
Publ. Info: -
Pages: - Volume / Issue: 16 (2) Sequence Number: - Start / End Page: 689 - 709 Identifier: CoNE: https://publications.pik-potsdam.de/cone/journals/resource/140507
Publisher: Copernicus