open access

https://hdl.handle.net/2115/65187

Comparison and verification of enthalpy schemes for polythermal glaciers and ice sheets with a one-dimensional model

英語

Enthalpy method

2015-04-30

National Institute of Polar Research

Polar Science

9(2)

196-207

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.

学術雑誌論文

Accepted Manuscript

Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International

©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/

1873-9652