HUSCAP logo Hokkaido Univ. logo

Hokkaido University Collection of Scholarly and Academic Papers >
Institute of Low Temperature Science >
Peer-reviewed Journal Articles, etc >

Comparison of thermodynamics solvers in the polythermal ice sheet model SICOPOLIS

Creative Commons License

Files in This Item:
Greve_Blatter_2016_PolarSci -.pdf2.28 MBPDFView/Open
Please use this identifier to cite or link to this item:http://hdl.handle.net/2115/68069

Title: Comparison of thermodynamics solvers in the polythermal ice sheet model SICOPOLIS
Authors: Greve, Ralf Browse this author
Blatter, Heinz Browse this author
Keywords: Ice sheet
Thermodynamics
Polythermal ice
Enthalpy method
Modeling
Issue Date: 2-Jan-2016
Publisher: Elsevier
Journal Title: Polar Science
Volume: 10
Issue: 1
Start Page: 11
End Page: 23
Publisher DOI: 10.1016/j.polar.2015.12.004
Abstract: In order to model the thermal structure of polythermal ice sheets accurately, energy-conserving schemes and correct tracking of the coldetemperate transition surface (CTS) are necessary. We compare four different thermodynamics solvers in the ice sheet model SICOPOLIS. Two exist already, namely a twolayer polythermal scheme (POLY) and a single-phase cold-ice scheme (COLD), while the other two are newly-implemented, one-layer enthalpy schemes, namely a conventional scheme (ENTC) and a meltingCTS scheme (ENTM). The comparison uses scenarios of the EISMINT Phase 2 Simplified Geometry Experiments (Payne et al., 2000, J. Glaciol. 46, 227-238). The POLY scheme is used as a reference against which the performance of the other schemes is tested. Both the COLD scheme and the ENTC scheme fail to produce a continuous temperature gradient across the CTS, which is explicitly enforced by the ENTM scheme. ENTM is more precise than ENTC for determining the position of the CTS, while the performance of both schemes is good for the temperature/water-content profiles in the entire ice column. Therefore, the one-layer enthalpy schemes ENTC and ENTM are viable, easier implementable alternatives to the POLY scheme with its need to handle two different numerical domains for cold and temperate ice. (C) 2016 Elsevier B.V. and NIPR. All rights reserved.
Rights: ©2016, Elsevier. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/
http://creativecommons.org/licenses/by-nc-nd/4.0/
Type: article (author version)
URI: http://hdl.handle.net/2115/68069
Appears in Collections:低温科学研究所 (Institute of Low Temperature Science) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)

Submitter: Greve Ralf

Export metadata:

OAI-PMH ( junii2 , jpcoar )


 

Feedback - Hokkaido University