根据我国第22次南极科学考察现场观测数据以及现场采集的冰芯分析数据,考虑太阳辐射、云量等因素对南极海冰厚度变化的影响,利用最小二乘法辨识气温与冰表面温度之间的关系,并以辨识结果所得到的冰表面温度作为上边界条件,依据实测情况取结冰点-1.81℃作为下边界条件,对一维热传导方程采用Crank-Nicolson格式进行离散,对南极中山站内拉湾附近海域海冰厚度变化进行了数值模拟,并与现场观测的海冰厚度数据进行了对比。
This paper studies the numerical simulation of Antarctic sea ice thickness, adopting implicit Crank-Nicolson format to discrete the one-dimensional heat transfer equation. We have based on the observation data of the 22nd CHINARE(2006) and ice core analysis data collected at the scene and considered Solar radiation, cloud cover and other factors on the Antarctic sea ice thickness changes. We take the recognition results by the method of least squares which is the relationship between air temperature and ice surface temperature as the upper boundary conditions, and according to the actual situation to take freezing point(--1.81 ~C ) as the lower boundary condition. This paper studies the numerical simulation of sea ice thickness around Nella Bay of Zhongshan Station in Antarctica, adopting implicit Crank-Nicolson format to discrete the one-dimensional heat transfer equation, and compares simulative ice thickness and measured ice thickness.