中国极地研究中心（中国极地研究所）极地多圈层研究团队在深度学习自动提取冰盖内部层位结构方面取得新进展

       联合国政府间气候变化专门委员会第六次评估报告（IPCC，AR6，2021年）指出，在温室气体高排放情景下，极地冰盖的贡献将成为估计未来海平面上升的最大不确定性来源。然而，通过数值模式量化冰盖对海平面的贡献受到极地冰盖演变的不确定性的阻碍。当前的冰盖模拟仍然受到诸如基岩地形、地热通量、冰层结构等不完整边界（或初始）条件以及冰流过程过度参数化的限制。20世纪60年代以来，在南极和格陵兰冰盖收集了大量的冰雷达数据。研究表明，冰盖内部和底部的冰雷达成像是一种获取冰底冰岩界面（冰下地形）和冰内等时层（代表冰盖古表面）的有效方式。它们已被广泛用于获取冰层年龄-深度关系与寻找老冰、估算冰底粗糙度与冰下温度、反演古冰流、推断冰下水分布、反推古冰盖表面积累率和物质平衡时空变化等。近几十年来，虽然对冰盖内部结构的理解取得了巨大的进步，然而即便在一些最基础的问题上仍远未得到解决。主要表现在冰盖几何层位（特别是内部等时冰层）的形成机理及其快速自动提取等方面。目前，无论是在技术、理论分析还是数值模拟等层面都对此缺乏基本的理解。迄今为止，大多数提取冰盖冰岩界面和内部等时层仍采用人工半经验追踪方法。人工半经验追踪内部冰层极其费时费力。例如，沿着20000千米长的冰雷达断面追踪20层需要熟练的操作员花至少10年时间才能完成。

      鉴于人工半经验追踪方法的极低效率，迫切需要开发自动的方法来追踪冰下基岩与冰层。然而，自动追踪内部等时层是极具挑战的研究难题，目前仍不存在一个具备对不同冰盖区域都有效的内部等时层自动提取数学物理模型。为此国际南极科学研究委员会（SCAR）于2018年批准了一个南极冰盖等时冰层模型构建（AntArchitecture）行动小组。AntArchitecture的目标是利用雷达探测成像的内部等时层和表面，通过发展新的技术开发一个覆盖整个大陆、年龄深度的南极冰模型，以便用于确定南极冰盖在过去演化周期中的稳定性。 

      最近，中国极地研究中心（中国极地研究所）（简称“极地中心”）自然资源部极地科学重点实验室冰盖多圈层相互作用与气候变化研究团队联合南方科技大学地球与空间科学系陈晓非院士团队，在地球科学和遥感领域期刊IEEE Transaction on Geoscience and Remote Sensing（IF=8.125，以下简称IEEE TGRS）和Remote Sensing（IF=5.349）上发表了3篇相关研究论文，尝试解决自动快速提取冰盖几何层位和冰层连续性指数、冰底粗糙度定量计算等问题。该系列研究通过构建基于卷积神经网络的人工智能深度学习模型，结合不同的地球物理数据预处理方法，从雷达图像去噪提升几何层位成像质量入手，进行了基岩界面和等时层自动提取的相关研究。在此基础上，提出了一种基于地球物理雷达数据去噪、深度学习提取冰盖内/底部几何轮廓自动提取后快速计算冰层连续性指数、冰底粗糙度的新方法。

图1：自动提取冰盖基岩界面和冰层结构的深度学习模型EisNet网络结构

       在Dong 等人 （2021）中，我们基于卷积神经网络构建了一种新的可应用于雷达图像回波噪声剔除的EisNet冰层自动提取网络模型（图1）。该模型具有以下特征：（1）使用了一种新的雷达图像数据合成方法用以解决因训练数据与对应的层位标签较少带来的训练困难，适用于EisNet不同分支路径的两类数据集；（2）训练后的EisNet经过合成数据与实际观测数据的双重测试，对冰层与冰岩界面两种目标的提取效果较好，具有较强的泛化能力；（3）EisNet的出现极大地提升了冰岩界面提取的效率，并为提取冰岩界面提供了统一稳定的定量化标准，同时为内部等时层的提取提供了有效的实现方法。

      在此基础上，针对机载冰雷达探测数据的复杂去噪与冰盖几何层位的特征提取速率问题，Tang等人（2022a）建议将EisNet方法引入到数据处理过程中。在采用传统的F-K滤波和偏移相结合的方法去除雷达数据中的特征噪声和干扰后，聚焦雷达图像分层界面（图2）。通过比较滤波和深度学习两种方法的处理效果，提出了一种结合滤波和深度学习的融合方法来优化机载雷达数据的几何特征提取。现场数据检验表明，在不同的雷达数据噪声条件下，采用不同的方法和策略可以显著提高层位特征提取的效果。

图2. 雷达数据和处理结果实例.(a-c）分别是雷达数据按顺序进行的零时刻校正(二维插值以及去空气层信号）、KL滤波、F-K偏移处理后的结果

      内部层连续性指数和底部粗糙度是刻画冰盖内部冰层或底部地形与冰流关系的主要指标。在过去的二十年，已经发展出几种方法来定量计算冰层连续性指数和底部粗糙度，即利用回波信号强度的垂向梯度绝对值平均和基岩界面的水平起伏标准方差表示。然而直到最近，这些方法仍受限于未经处理的雷达数据集的规模和质量，而未得到进一步的应用。Tang等人（2022b）提出了一种基于地球物理雷达图像去噪、EisNet快速提取冰盖内/底部几何轮廓基础上的冰层连续性指数、冰底粗糙度计算新方法。使用地面和机载雷达数据进行了东南极中山站至Dome A车载雷达和机载雷达的数据成像、计算结果的对比（图3）。通过从去除噪声后的雷达数据中提取有关冰川分层连续性的信息，并利用以往人工半经验提取的层位信息进行了相关计算验证（图4）。研究认为该方法有潜力成为未来从冰盖内部结构中提取更多冰盖动力学历史信息的有效途径。

图3. (a)人工半经验和EisNet提取的中山站至Dome A机载雷达断面的表面高程和底部高程图；其中(b)和(c)分别是(a)中绿色矩形框对应的冰表面和冰岩界面附近局部雷达图像；(d)人工半经验与EisNet插值平滑处理后的基岩高程差直方图

图4. 中山站至Dome A机载雷达断面连续指数随深度变化结果