Towards a fully unstructured ocean model for ice shelf cavity environments: Model development and verification using the Firedrake finite element framework

Scott, William I.; Kramer, Stephan C.; Holland, Paul R. ORCID: https://orcid.org/0000-0001-8370-289X; Nicholls, Keith W. ORCID: https://orcid.org/0000-0002-2188-4509; Siegert, Martin J.; Piggott, Matthew D.. 2023 Towards a fully unstructured ocean model for ice shelf cavity environments: Model development and verification using the Firedrake finite element framework. Ocean Modelling, 182, 102178. 20, pp. https://doi.org/10.1016/j.ocemod.2023.102178

Before downloading, please read NORA policies.

Preview

Text (Open Access)
© 2023 The Author(s). Published by Elsevier Ltd.
1-s2.0-S1463500323000197-main.pdf - Published Version
Available under License Creative Commons Attribution 4.0.
Download (3MB) | Preview

Official URL: https://www.sciencedirect.com/science/article/pii/...

Abstract/Summary

Numerical studies of ice flow have consistently identified the grounding zone of outlet glaciers and ice streams (the region where ice starts to float) as crucial for predicting the rate of grounded ice loss to the ocean. Owing to the extreme environments and difficulty of access to ocean cavities beneath ice shelves, field observations are rare. Estimates of melt rates derived from satellites are also difficult to make near grounding zones with confidence. Therefore, numerical ocean models are important tools to investigate these critical and remote regions. The relative inflexibility of structured grid models means, however, that they can struggle to resolve these processes in irregular cavity geometries near grounding zones. To help solve this issue, we present a new nonhydrostatic unstructured mesh model for flow under ice shelves built using the Firedrake finite element framework. We demonstrate our ability to simulate full ice shelf cavity domains using the community standard ISOMIP+ Ocean0 test case and compare our results against those obtained with the popular MITgcm model. Good agreement is found between the two models, despite their use of different discretisation schemes and the sensitivity of the melt rate parameterisation to grid resolution. Verification tests based on the Method of Manufactured Solutions (MMS) show that the new model discretisation is sound and second-order accurate. A main driver behind using Firedrake is the availability of an automatically generated adjoint model. Our first adjoint calculations, of sensitivities of melt rate with respect to different inputs in an idealised grounding zone domain, are promising and point to the ability to address a number of important questions on ocean influence on ice shelf vulnerability in the future.

Item Type:	Publication - Article
Digital Object Identifier (DOI):	https://doi.org/10.1016/j.ocemod.2023.102178
ISSN:	14635003
Additional Keywords:	FEM, Ice shelf cavities, Firedrake, Adjoint, Verification, MMS
Date made live:	27 Feb 2023 10:56 +0 (UTC)
URI:	https://nora.nerc.ac.uk/id/eprint/534086

Item Type:

Publication - Article

Digital Object Identifier (DOI):

https://doi.org/10.1016/j.ocemod.2023.102178

ISSN:

14635003

Additional Keywords:

FEM, Ice shelf cavities, Firedrake, Adjoint, Verification, MMS

Date made live:

27 Feb 2023 10:56 +0 (UTC)

URI:

https://nora.nerc.ac.uk/id/eprint/534086