EGU21-6854, updated on 25 Jan 2023
https://doi.org/10.5194/egusphere-egu21-6854
EGU General Assembly 2021
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.

Implementation of higher-order advection schemes in a numerical ice-sheet model for ice-core studies

Fuyuki Saito1, Ayako Abe-Ouchi2, and Takashi Obase2
Fuyuki Saito et al.
  • 1JAMSTEC, Yokohama, Japan (saitofuyuki@jamstec.go.jp)
  • 2AORI, Univ. Tokyo, Kashiwa, Japan

Computation of temperature and age fields by numerical ice-sheet models is an important issue for ice-core related studies.  Generally the evolution of temperature and/or age in an ice-sheet model is formulated using an advection equation.  There are many variation of the formulation, which differ in numerical aspects such as stability, accuracy, numerical diffusivity, conservation and/or computational costs.  Saito et al (2020, GMD) implement Rational Constrained Interpolation Profile (RCIP) scheme on vertical 1-d age computation of ice sheet, and demonstrate its efficiency, in particular, to preserve surface mass balance properties recorded at the deposit in terms of annual layer thickness.  Successively, we have been extending the development using RCIP or similar higher-order advection schemes on 3-d age or temperature computation.  In this study, we demonstrate 1-d temperature computation by various numerical schemes including classical upwind schemes and compare the accuracy of those schemes.

How to cite: Saito, F., Abe-Ouchi, A., and Obase, T.: Implementation of higher-order advection schemes in a numerical ice-sheet model for ice-core studies, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6854, https://doi.org/10.5194/egusphere-egu21-6854, 2021.

Displays

Display file