几十年来，格陵兰冰盖经历着日益严重的物质损失，格陵兰冰盖已成为全球海平面变化的主要驱动因素。而冰盖表面高程变化是反应气候变化最直接的指标之一，结合适当的冰雪密度模型，可由冰面高程变化估算冰盖的质量变化。

我中心王泽民教授团队针对冰盖高程变化研究中存在的卫星间系统偏差改正不完善，不同卫星任务衔接处存在阶跃误差，进而影响冰盖高程变化估计精度的问题，以及在构建长时序冰面高程产品时，因极地冰盖的地面轨迹分布不均匀，在地形复杂、陡峭或者变化剧烈的冰盖边缘区域观测精度差而需要进行插值所产生的过度插值误差问题，提出了引入任务间升降轨道偏差概念的新平差策略，结合基于经验正交函数重构插值的格陵兰长时序冰面高程产品构建方法，联合多源卫星测高数据建立了1991–2020年间近30年的月均5km格网的格陵兰冰面高程时间序列产品。

与欧洲空间局格陵兰冰盖气候变化行动计划的冰面高程变化产品进行对比评估发现，该方法能够更加有效地消除任务间的系统偏差，得到更加自洽和可靠的冰面高程时间序列。相比于欧洲空间局的产品，该产品的冰面高程变化率的平均精度（RMS）和误差的离散度（P₉₀−P₁₀）分别提高了19.3%和8.9%，在2006–2010到2010–2014五个时期，平均精度和误差的离散度更是分别提高了30.9%和19.0%。 

图2 格陵兰冰盖（a）Jakobshavn Isbræ，（b）Upernavik Isstrøm，（c）Zachariæ Isstrøm和Nioghalvfjerdsfjorden与（d）King Christian X Land等冰面高程变化特征区域的部分格网点的时间序列

通过冰面高程不同时间时期的变化趋势分析，发现格陵兰冰面高程的长期变化同样具有显著的空间变化特征（如图1和图2所示）。格陵兰冰盖的沿岸区域多展现出长期连续地变薄趋势，在多个溢出冰川流域（如Jakobshavn Isbræ、Upernavik Isstrøm、Zachariæ Isstrøm等）均出现了明显的加速变薄过程。格陵兰冰盖在1991-2020 年期间的体积损失速率为53.8±4.5km³yr^-1，加速度为 2.2±0.3km³yr^-2，八个冰流域中有六个显示出体积损失，最剧烈的体积发生在格陵兰冰盖的西海岸（如图3所示）。

图3 格陵兰冰盖8个冰流域的冰盖体积变化时间序列

无论是变薄还是变厚，速率都不会保持不变，总是有加速或减速的时期，这是气候变化对冰盖变化的驱动力的明显痕迹。例如，由北大西洋涛动（NAO）驱动的2012年异常融化和2016-2017年以来的积累在上述区域的时间序列中清晰可见。

该研究成果近日在地学顶刊《Earth System Science Data》（IF=11.333）上发表，武汉大学中国南极测绘研究中心副研究员张保军为第一作者，王泽民教授和刘婷婷副教授为通讯作者。