Implementation of higher-order vertical finite elements in ISSM v4.13 for improved ice sheet flow modeling over paleoclimate timescales

Cuzzone, Joshua K.; Morlighem, Mathieu; Larour, Eric; Schlegel, Nicole; Seroussi, Helene

doi:https://doi.org/10.5194/gmd-11-1683-2018

Articles | Volume 11, issue 5

https://doi.org/10.5194/gmd-11-1683-2018

© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

https://doi.org/10.5194/gmd-11-1683-2018

© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Articles | Volume 11, issue 5

Development and technical paper

|

03 May 2018

Development and technical paper |

| 03 May 2018

Implementation of higher-order vertical finite elements in ISSM v4.13 for improved ice sheet flow modeling over paleoclimate timescales

Joshua K. Cuzzone, Mathieu Morlighem, Eric Larour, Nicole Schlegel, and Helene Seroussi

Download

Final revised paper (published on 03 May 2018)
Supplement to the final revised paper
Preprint (discussion started on 11 Jan 2018)
Supplement to the preprint

Interactive discussion

Status: closed

AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment

- Printer-friendly version

- Supplement

RC1: 'Review', Anonymous Referee #1, 17 Feb 2018
RC2: 'Review of Cuzzone et al. for Geoscientific Model Development, 2018', Anonymous Referee #2, 22 Feb 2018
AC1: 'Response to reviewer #1', Joshua Cuzzone, 05 Apr 2018
AC2: 'Response to reviewer #2', Joshua Cuzzone, 05 Apr 2018

Peer-review completion

AR: Author's response | RR: Referee report | ED: Editor decision

AR by Joshua Cuzzone on behalf of the Authors (09 Apr 2018) Manuscript

ED: Publish as is (14 Apr 2018) by Jeremy Fyke

AR by Joshua Cuzzone on behalf of the Authors (19 Apr 2018) Manuscript

Short summary

This paper details the implementation of higher-order vertical finite elements in the Ice Sheet System Model (ISSM). When using higher-order vertical finite elements, fewer vertical layers are needed to accurately capture the thermal structure in an ice sheet versus a conventional linear vertical interpolation, therefore greatly improving model runtime speeds, particularly in higher-order stress balance ice sheet models. The implications for paleoclimate ice sheet simulations are discussed.