HUSCAP logo Hokkaido Univ. logo

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

Creative Commons License

Files in This Item:
Blatter_Greve_2015_PolarSci -.pdf1.17 MBPDFView/Open
Please use this identifier to cite or link to this item:http://hdl.handle.net/2115/65187

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
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/
Type: article (author version)
URI: http://hdl.handle.net/2115/65187
Appears in Collections:雑誌発表論文等 (Peer-reviewed Journal Articles, etc)

Submitter: Greve Ralf

 

Feedback - Hokkaido University