Hokkaido University Collection of Scholarly and Academic Papers >
Institute of Low Temperature Science >
Peer-reviewed Journal Articles, etc >
Comparison and verification of enthalpy schemes for polythermal glaciers and ice sheets with a one-dimensional model
This item is licensed under:Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
Title: | Comparison and verification of enthalpy schemes for polythermal glaciers and ice sheets with a one-dimensional model |
Authors: | Blatter, Heinz Browse this author | Greve, Ralf Browse this author →KAKEN DB |
Keywords: | Glacier | Ice sheet | Polythermal ice | Modeling | Enthalpy method |
Issue Date: | 30-Apr-2015 |
Publisher: | National Institute of Polar Research |
Journal Title: | Polar Science |
Volume: | 9 |
Issue: | 2 |
Start Page: | 196 |
End Page: | 207 |
Publisher DOI: | 10.1016/j.polar.2015.04.001 |
Abstract: | The enthalpy method for the thermodynamics of polythermal glaciers and ice sheets is tested and verified by a one-dimensional problem (parallel-sided slab). The enthalpy method alone does not include explicitly the transition conditions at the cold-temperate transition surface (CTS) that separates the upper cold from the lower temperate layer. However, these conditions are important for correctly determining the position of the CTS. For the numerical solution of the polythermal slab problem, we consider a two-layer front-tracking scheme as well as three di erent one-layer schemes (conventional one-layer scheme, one-layer melting CTS scheme, one-layer freezing CTS scheme). Computed steady-state temperature and water-content profiles are verified with exact solutions, and transient solutions computed by the one-layer schemes are compared with those of the two-layer scheme, considered to be a reliable reference. While the conventional one-layer scheme (that does not include the transition conditions at the CTS) can produce correct solutions for melting conditions at the CTS, it is more reliable to enforce the transition conditions explicitly. For freezing conditions, it is imperative to enforce them because the conventional one-layer scheme cannot handle the associated discontinuities. The suggested numerical schemes are suitable for implementation in three dimensionalglacier and ice-sheet models. |
Rights: | ©2015, 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/65187 |
Appears in Collections: | 低温科学研究所 (Institute of Low Temperature Science) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)
|
Submitter: Greve Ralf
|