多热冰盖模型SICOPOLIS中热力学求解器的比较

为了准确模拟多热冰盖的热结构，需要能量守恒格式以及冷-温过渡面（CTS）的正确追踪。我们在冰盖模型SICOPOLIS中比较了四种不同的热力学求解器。其中两种已存在，即两层多热格式（POLY）和单相冷冰格式（COLD），另两种为新实现的单层焓格式，即常规格式（ENTC）和融化CTS格式（ENTM）。比较使用了EISMINT第二阶段简化几何实验（Payne et al., 2000, J. Glaciol. 46, 227-238）的场景。以POLY格式为参考检验其他格式的性能。COLD格式和ENTC格式均未能产生跨CTS的连续温度梯度，而该梯度由ENTM格式显式强制实现。在确定CTS位置方面，ENTM比ENTC更精确，同时两种格式在整个冰柱的温度/含水量剖面表现均良好。因此，单层焓格式ENTC和ENTM是POLY格式（需处理冷冰和温冰两个不同数值域）的可行且更易实现的替代方案。