北京大学地球与空间科学学院遥感与地理信息系统研究所18级本科生徐玥，在本科期间参与北京大学全球变化水循环遥感大数据实验室科研项目，近期以第一作者身份在《Geophysical Research Letters》(GRL)发表关于利用激光测高卫星数据反演南极海冰厚度的研究成果。GRL是地学领域的顶级期刊之一(IF:4.72)，发表地球科学领域方面高质量的创新研究成果。

该工作针对国际上激光测高卫星南极海冰厚度反演精度有待提高、新发射的Ice, Cloud and Land Elevation Satellite-2 (ICESat-2，2018至今)与ICESat(2003-2009)数据之间的一致性亟待检验等问题，提出了改进的逐点动态冰雪比单层法One-Layer Method-improved(OLMi)，给出全面的误差分析，分析了2018-2020年逐月南极海冰厚度时空分布规律，同时探讨了ICESat和ICESat-2在空间覆盖度、几何采样和高度反演算法等方面的数据一致性。该工作大大推进了高精度、长时间序列环南极海冰关键参数数据集的生产，为全球气候变化研究提供了重要支撑。

背景

南极海冰是地球冰雪圈的重要组成部分之一，不仅影响着南半球的海气热量交换，同时也是气候变化的“温度计”。目前，使用被动微波卫星对南极海冰面积进行的监测已经持续了四十多年，但是长时间序列高精度环南极海冰厚度产品依然缺乏。

新一代激光雷达测高卫星ICESat-2在2018年9月成功发射和运行，发布了时空分辨率更高的激光测高数据，在理论上已经可以获取逐月尺度南极海冰的厚度分布情况。结合第一代ICESat卫星在2003-2008年正常运行期间的激光测高数据，可以联合分析、探究全球变化对南极海冰潜在的影响。

数据和方法

本研究使用了2003-2008年Lowest Level Elevation算法生产的ICESat海冰总干舷高度（Total Freeboard, TF）数据、2018-2020年ICESat-2的ATL10海冰干舷产品数据。

由于南极的积雪厚度较深，本研究改进了OLM单层法（OLMi），将雪层厚度从总中减去，即，更好地反映南极海冰的厚度特性。

反演公式的误差主要来源于海水、海冰、冰上积雪的密度、海冰总干舷高度和经验公式中的独立参数()，根据高斯误差传播公式，海冰厚度的反演误差可以被初步估计。

结果和分析

1、月尺度时空分布

第一代激光测高卫星ICESat的数据受激光点数量和工作时间所限，只能支撑全南极海域分季度的海冰时序分析，ICESat-2的发射频率、光束数量都有了显著的提升，能够支撑逐月尺度全南极海域海冰情况的展示和分析，由此对南极海冰厚度在一年内生长和融化的情况有了更加精细的刻画（图 2）：

2、一致性评估

ICESat和ICESat-2的数据差异主要体现在空间覆盖度、几何采样和高度反演算法三个方面。

ICESat的数据覆盖了14.3%的海冰浓度大于50%的格网，而ICESat-2能够覆盖39.6%，大幅提高了数据的空间覆盖度。ICESat的足迹点大小约70m，而ICESat-2只有17m。将ICESat-2的数据重采样到70m之后，计算的厚度结果比原结果低2.0%（0.005m）左右。

Freeboard反演算法是两代数据差异最大的地方，在将ICESat的Lowest Level Elevation（LLE）算法应用在ICESat-2的海冰高度数据后，发现LLE算法比ICESat-2的ATL10算法计算结果低估厚度14%（0.036m）左右。

作者信息：

Xu, Y., Li, H., Liu, B., Xie, H., Ozsoy‐Cicek, B., 2021.Deriving Antarctic Sea‐Ice Thickness From Satellite Altimetry and EstimatingConsistency for NASA's ICESat/ICESat‐2 Missions. Geophysical Research Letters48, e2021GL093425. doi: https://doi.org/10.1029/2021GL093425