Journal Article Comparison of thermodynamics solvers in the polythermal ice sheet model SICOPOLIS

10 ( 1 )  , pp.11 - 23 , 2016-03
In order to model the thermal structure of polythermal ice sheets accurately, energy-conserving schemes and correct tracking of the cold–temperate transition surface (CTS) are necessary. We compare four different thermodynamics solvers in the ice sheet model SICOPOLIS. Two exist already, namely a two-layer 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 melting-CTS 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.
ポリサーマル氷河の氷温度分布を正確に計算するためには、エネルギー保存の法則に従うとともに、融点にある氷と融点以下の氷の境界(CTS)の位置を正確に知ることが必要である。我々は4つの異なる熱力学方程式の解法を氷床モデル(SICOPOLIS)で比較した。2つの手法は既存の二層ポリサーマル法(POLY)と単一相寒冷氷法(COLD)で、他の2つは本研究で初めて導入した単一相エンタルピー法、すなわち従来法(ENTC)と融解CTS法(ENTM)である。これらの手法の比較には、第2期EISMINTの単純地形シナリオを用いた(Payne et al., 2000, J. Glaciol. 46, 227–238)。基準となるPOLYの計算結果と比較することで、他の手法について検討を行った。COLDおよびENTCはCTSにおいて温度勾配を連続に保つことができなかったが、ENTMでは強制によってその連続が常に成り立っている。ENTCとENTMのどちらの手法からも氷河内の氷温と含水率の分布を推定することができるが,ENTMはCTSの位置に関してはENTMよりも正確であった。以上の結果から、融点以下と融点以上の氷について領域が設定されるPOLYと比較して、単一相エンタルピー法であるENTCとENTMはより実用的で簡便な手法といえる。

Number of accesses :  

Other information