一维海冰热力模式的守恒型差分格式和数值模拟

THE CONSERVATIVE DIFFERENCE SCHEME AND NUMERICAL SIMULATION OF ONE DIMENSIONAL THERMODYNAMIC SEA ICE MODEL

根据前人研究的—维海冰热力模式,考虑到方程的特点,采用积分插值法构造了守恒型差分格式,并对海冰热力模式所描述的冰内热扩散传导和冰底厚度的变化进行了数值模拟.当冰面覆盖积雪时,用该解法同样计算了雪层内的热扩散,考虑冰雪之间的热力相互作用后,得到冰厚在有雪盖情况下的变化.数值计算表明由积分插值法构造的数值格式计算稳定,对时间和空间步长的选择自由,可以根据需要进行不同时间段长的数值模拟,该解法得到的冰内温度分布与实际观测到的近似.当有大气变化状况的预报时,该模式结合大气和冰表面热交换过程可进行冰厚的预报.

This paper is based on the former research of the one dimensional thermodynamic sea ice model. And the conservative difference scheme of 1-D model is created here by using the integral interpolation method. Some numerical calculations including temperature profile within ice layer, ablation and accretion of ice bottom and the ice thickness change with snow layer covered are presented. The simulation results shows: 1)Numerical scheme has good stability. 2) Arbitrary time and vertical space step are allowed to do the simulation. 3) The distribution of temperature profile within ice are similar to the real observation data. If it is further considering the couple of atmosphere and ice surface, the 1-D model could be used to do the ice thickness forecast.