冰川冻土, 2022, 44(3): 863-871 doi: 10.7522/j.issn.1000-0240.2022.0083

冰冻圈与全球变化

中国极地冰盖数值模拟研究进展

赵励耘,1, 唐学远2, 张通3, 冷伟4, 艾松涛5, 刘成彦6, 王玉哲7, 刘岩1, 岳超1

1.北京师范大学 全球变化与地球系统科学研究院，北京 100875

2.自然资源部 极地科学重点实验室 中国极地研究中心，上海 200136

3.北京师范大学 地表过程与资源生态国家重点实验室，北京 100875

4.中国科学院 数学与系统科学研究院，北京 100190

5.武汉大学 中国南极测绘研究中心，湖北 武汉 430079

6.南方海洋科学与工程广东省实验室（珠海），广东 珠海 519080

7.山东师范大学 地理与环境学院，山东 济南 250014

Progress in numerical simulation of polar ice sheets in China

ZHAO Liyun,1, TANG Xueyuan2, ZHANG Tong3, LENG Wei4, AI Songtao5, LIU Chengyan6, WANG Yuzhe7, LIU Yan1, YUE Chao1

1.College of Global Change and Earth System Science，Beijing Normal University，Beijing 100875，China

2.Key Laboratory of Polar Science of Ministry of Natural Resources （MNR），Polar Research Institute of China，Shanghai 200136，China

3.State Key Laboratory of Earth Surface Processes and Resource Ecology，Beijing Normal University，Beijing 100875，China

4.Academy of Mathematics and Systems Science，Chinese Academy of Sciences，Beijing 100190，China

5.Chinese Antarctic Center of Surveying and Mapping，Wuhan University，Wuhan 430079，China

6.Southern Marine Science and Engineering Guangdong Laboratory （Zhuhai），Zhuhai 519080，Guangdong，China

7.College of Geography and Environment，Shandong Normal University，Jinan 250014，China

收稿日期: 2021-12-27   修回日期: 2022-04-02  

基金资助: 国家自然科学基金极地基础科学前沿专项.  41941006
国家重点研发计划项目.  2018YFC1406104

Received: 2021-12-27   Revised: 2022-04-02  

作者简介 About authors

赵励耘，副教授，主要从事冰川变化和冰川数值模拟研究.E-mail：zhaoliyun@bnu.edu.cn , E-mail：zhaoliyun@bnu.edu.cn

摘要

冰盖数值模拟是一种基于多源观测数据，通过构建并求解冰流动力学方程组，理解冰流运动物理机制以及诊断和预估其演化过程的方法，目前已被广泛应用于冰盖变化研究。本文简要介绍了极地冰盖数值模拟方法，归纳综述了近十余年我国学者在极地冰盖数值模拟方面的研究进展，厘清我国在冰盖数值模拟领域遇到的瓶颈和关键问题。阐述了如何与我国的极地冰盖科考优势区域深度结合，协同多源强化观测和数值模拟，研发和改进冰盖模式，提高冰盖模拟能力，对定量估算极地冰盖的物质平衡及其对未来海平面上升的影响做出实质贡献。通过逐步发展冰盖模式的研究能力，有望将来在冰盖关键动力过程和机制的科学认识上有所突破。

关键词： 极地冰盖 ; 数值模拟 ; 模式发展 ; 圈层耦合

Abstract

Numerical simulation of ice sheets is a method to understand the physical mechanism of ice flow movement， diagnose and predict its evolution by constructing and solving the dynamic equations of ice flow based on multi-source observation data. Numerical simulation has been widely used in the study of ice sheet change. In this review， we briefly introduce the numerical simulation methods of polar ice sheets， summarize the research progress of Chinese scholars in polar ice sheets numerical simulation in the last decade， and identify the bottlenecks and key problems encountered in the field of ice sheet numerical simulation in China. We expound how to combine the advantageous areas of polar ice sheet scientific research in China， cooperate with multi-source enhanced observation and numerical simulation， develop and improve the ice sheet model， enhance the ice sheet simulation ability， and contribute to the quantitative estimation of polar ice sheets mass balance and its impact on future sea level rise. By gradually developing the research ability of ice sheet model， it is expected to make a breakthrough in the scientific understanding of key dynamic processes and mechanisms of ice sheets in the future.

Keywords： polar ice sheets ; numerical simulation ; model development ; sphere coupling

PDF (2078KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

赵励耘, 唐学远, 张通, 冷伟, 艾松涛, 刘成彦, 王玉哲, 刘岩, 岳超. 中国极地冰盖数值模拟研究进展[J]. 冰川冻土, 2022, 44(3): 863-871 doi:10.7522/j.issn.1000-0240.2022.0083

ZHAO Liyun, TANG Xueyuan, ZHANG Tong, LENG Wei, AI Songtao, LIU Chengyan, WANG Yuzhe, LIU Yan, YUE Chao. Progress in numerical simulation of polar ice sheets in China[J]. Journal of Glaciology and Geocryology, 2022, 44(3): 863-871 doi:10.7522/j.issn.1000-0240.2022.0083

0 引言

极地冰盖（图1）对全球热量平衡和海平面上升具有重要的影响。IPCC第六次评估报告指出，气候变化评估中最大的不确定性之一来自对极地冰盖对未来海平面上升贡献量的估计，格陵兰冰盖物质损失对20世纪全球海平面变化贡献了29%，南极冰盖物质损失区域差异性大，其中西南极冰盖底部由于受到海洋环流影响融化速率相对较快［1-2］。此外，东南极尤其是威尔克斯盆地（Wilkes Basin）也发现潜在的不稳定性［3］。

由于极地冰盖气候恶劣、难以接近，评估当前极地冰盖变化的主要手段是基于遥感监测，在部分区域结合野外实地测量。遥感技术可以对冰盖实现全覆盖、连续性监测。地面实地测量难以大范围实施，但是可以更精确地探测冰盖内部构造和底部地形等关键物理信息。遥感技术和地面实地测量等多源观测为研究极地冰盖变化提供了必备的数据资料（表1）。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   南极冰盖（a）；北极的格陵兰冰盖、斯瓦尔巴群岛和冰岛（b）［在图中标注了我国建立的台站（长城站、中山站、昆仑站、泰山站、罗斯海新站、黄河站）以及本文提到的一些研究区域或冰川位置（威尔克斯盆地、伊丽莎白公主地、埃默里冰架、普利兹湾、Totten冰川、雅各布港冰川、瓦特纳冰帽）。Dome A与昆仑站非常近，Austre Lovénbreen冰川和Pedersenbreen冰川与黄河站距离非常近，不再标注］

Fig. 1   Antarctic Ice Sheet （a）； Greenland Ice Sheet， Svalbard Islands and Iceland in the Arctic （b） ［The stations established by China （Great Wall Station， Zhongshan Station， Kunlun Station， Taishan Station， Ross Sea New Station and Yellow River Station） and some study regions or glacier locations mentioned in this paper （Wilkes Basin， Princess Elizabeth Land， Amery ice shelf， Prydz Bay， Totten Glacier， Jakobshavn Glacier， Vatnajökull ice cap） are marked in the figure. Dome A， very close to Kunlun Station， Austre Lovénbreen Glacier and Pedersenbreen Glacier almost next to Yellow River Station， are not marked］


表1   冰盖关键动力过程、模型中的相关变量和主要观测手段

Table 1  Key dynamic processes of ice sheet， relevant variables in the model and main observational approaches

关键动力过程模型中的相关变量主要观测手段
冰面水文过程冰面径流卫星遥感、无人机
冰面高程变化冰面物质平衡、冰流速卫星遥感、GPS、花杆、雪坑
冰下环境和冰下水文过程冰温、冰流速、冰底滑动系数、冰盖底部融化速率、有效压力、水流通量冰雷达、回声测深仪、钻孔
冰架崩解冰架末端退缩速率卫星遥感
冰架底部物质平衡冰架底部消融/冻结速率回声测深仪、钻孔
内部构造冰温、冰流速、冰晶组构、冰龄、杂质冰雷达、钻孔
触地线迁移触地线位置卫星遥感

新窗口打开| 下载CSV


基于多源观测数据作为模式边界条件和外部驱动，可以通过数值求解冰流动力学模型，详细描述冰流的运动或预测冰盖的未来演化。冰盖数值模拟通常包括两种：一种是诊断模拟（diagnostic modeling），用以诊断冰盖在稳态假设和特定边界条件及参数下的动态平衡状态；也可以用来测试冰盖对某些参数或边界条件的敏感性以及微小变化对模型行为的影响，从而确定可能反馈机制的性质。诊断模型可用于提高对控制特定冰流行为的参数的理解，或研究一般冰盖中一个或多个物理参数与边界条件的重要性。另一种是预估模拟（prognostic modeling），用以描述冰盖运动的过程，通过调整模型并使其包括不同的物理过程，以便更好地将模拟预估与实际情况相匹配，从而更好地预测冰盖的未来演变。预估模型既可直接用于预估（或预测）冰盖的未来变化；也可用于重建冰盖的历史演化过程，从而有助于准确理解冰盖未来可能发生的变化。

如果要考虑到作用在冰上的所有应力分量，full-Stokes模型可能是最合适的工具。在某些情形下忽略部分应力分量，可以得到简化近似方程。基于full-Stokes模型或者简化近似的冰流模型，各国研究团队陆续研发了一些用于模拟冰盖状态或者预估其未来变化的数值求解方案，即冰流模式［4］。虽然我国在冰盖数值模拟方向的研究起步较晚，但是目前也形成了三个自主研发的冰流模式——多温型陆地冰模型PoLIM［5］、三维并行冰盖演化模型［6］和高精度并行有限元冰盖模型［7-9］。PoLIM为二维Blatter-Pattyn高阶近似热-力耦合冰流模型，以焓方法计算冰川热力场，包含冰下水文模型，采用有限差分方法求解计算。张怀等开发了三维并行冰盖演化模型，进行了格陵兰冰盖的模拟演化［6］。冷伟等研发的三维full-Stokes和高阶近似冰盖模式求解器，使用了适用于大规模数值模拟的快速求解方法［7-8］，并采用了保持局部质量守恒的数值离散方法以提高计算精度［9］，具有热-力耦合和处理触地线进退过程的功能［10-11］。

冰盖与大气、基岩、海洋接触，在边界处存在复杂的过程，包括冰盖表面物质平衡、冰盖水文过程、冰下水热环境、冰架崩解、冰架底部消融和再冻结等［12-13］（表1、图2）。这些过程都会对冰流系统的物质平衡和稳定性产生重要的影响，是最近冰盖动力学研究领域的热点和难点。目前在理论和观测上对上述关键动力过程和机制的认识都存在很大不足，严重制约了模式的发展和模拟结果的置信度。尽管国际冰盖模式发展时间较长，但是其对上述过程多采取简单的参数化处理或者尚未考虑，其模拟结果仍存在很大的不确定性［4］。特别的，近年来开展的冰盖-海洋耦合模式比较计划，不同耦合模式对海洋强迫的模拟结果差异很大。因此，模式的改进需要以冰盖典型区域长期持续的高精度精细测量作为基础支撑以提升对冰盖关键过程和机制的理解。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   冰盖关键动力过程示意图（本示意图的地形更适合西南极冰盖，而不适用于东南极和格陵兰冰盖大部分地区）

Fig. 2   Schematic diagram of key dynamic process of ice sheet （The topography in this diagram is more suitable for West Antarctic Ice Sheet than most region in East Antarctic and Greenland Ice Sheet）


前些年已有数篇针对极地冰盖数值模拟的评述文章［14-15］。近十余年来，我国冰川学家在极地冰盖数值模拟方面产出了不少新的研究成果。本文简要介绍冰盖数值模拟方法，进而梳理了我国在极地冰盖数值模拟方面的新进展，并对未来的研究前景进行了展望。

1 冰盖数值模拟方法简述

开发冰盖动力学数值或解析模型的主要用途是为了更好地了解冰盖的流动行为以及它们如何对外部环境的强迫做出响应。这些模型通常基于一些描述冰流的基本定律或假设，经过简化后可以得到解析或数值解。冰盖动力学模型包括动力学方程和初始条件以及边界条件。动力学方程主要包括冰盖的质量、动量、能量守恒定律和本构方程Glen流动定律。目前最为完整的动力学方程是热-力耦合的三维full-Stokes方程，它将描述质量、动量守恒的Stokes方程和描述能量守恒的热传导方程进行耦合，适用于各种情形的冰盖模拟。在Stokes方程的基础上进行不同程度的简化近似，可以得到静水压力近似（hydrostatic approximation）、一阶近似（first order approximation）、浅冰盖近似（shallow ice approximation）、浅冰架近似（shallow shelf approximation）等［16-17］。在full-Stokes方程中忽略垂直切应力τzx 和τzy 沿着运动方向的导数，可以得到静水压力近似。在静水压力近似基础上，忽略速度竖直分量的水平导数，即假设速度竖直分量vz 仅依赖于竖直坐标z和时间t，vz =vz （z，t），得到一阶近似。在一阶近似基础上进一步简化，把冰盖看作是平行底面的剪切流，水平切应力τxz 和τyz 对水平方向的动量平衡起到重要作用，忽略法向应力偏量和竖直切应力τxy，可以得到浅冰盖近似模型；而如果应用在冰架上，假设冰架底面切应力为零，流速的水平分量不依赖于深度，即vx =vx （x，y，t），vy =vy （x，y，t），则得到浅冰架近似模型。浅冰盖近似和浅冰架近似可以结合成混合模型（hybrid model），应用到整个冰盖的模拟［18-20］。

黏性系数是动力学方程的关键参数，直接影响着冰流的速度大小，并且敏感地依赖于冰温和冰晶组构等因素。因此，热-力耦合的动力学方程的模拟结果相对更接近实际情形。full-Stokes方程及其近似方程是建立在冰各向同性的假设之上，由于实际冰盖冰晶组构各向异性等原因，冰流动力学方程中的重要参数（例如冰的黏性增强因子、冰底摩擦系数）可以通过求解反问题来反演，使得表面冰流速的模拟值与观测值误差最小［21-22］。原则上，对于full-Stokes冰盖模型而言，随着计算机技术的发展使得完整应力分布的充分计算已成为可能并得到逐步实现。但对于百年以上的积分，某些大陆尺度的冰盖演化模拟在数值计算上仍不可行。

技术上，冰盖动力学模型求解的空间域是基于地形数据建立的研究对象的数字高程模型。基于冰厚、海平面高程和冰底高程数据，可以判断出冰盖触地部分和冰架部分。在研究区域的边界（冰面、冰底、冰崖、侧面等）要给出合理的边界条件，主要包括：（1）冰面温度、表面物质平衡；（2）冰底触地部分的热力平衡方程（包括地热通量等）、运动状态等；（3）冰盖底部的压融点温度和物质平衡等；（4）冰崖处的海水压强和崩解过程等。其中，表面物质平衡和冰架底面物质平衡体现的是大气和海洋对冰盖作用的强迫场。冰盖表面物质平衡的获取途径较多，包括实地测量、再分析资料和区域气候模式的模拟结果等。冰架底面物质平衡的直接监测条件受限（主要是回声测深仪［23］和钻孔），缺乏基础数据。以往和目前的研究多采用参数化方案或者海洋模式的模拟结果进行估算。如何给出较为精确的冰架底面物质平衡是目前研究的主要难点之一。另外，冰盖底部的运动状态很难监测，而且对冰下水热环境变化非常敏感等，具有较大的不确定性。以往研究中的常用方法是在假定的滑动定律形式下通过数据同化反演底部滑动系数［24-25］。冰崖处的崩解过程对南极冰盖的稳定性非常重要，也是南极冰盖物质损失很大的不确定性来源。但是崩解过程如何在冰盖模型中精确地刻画是迄今尚未较好解决的难题。在以往的冰盖模拟研究中常见的处理方式有：冰崖的位置保持固定，或者利用冰裂隙深度、冰面径流、应变率、应力、冰厚等对崩解速率进行参数化［4，26-28］。

冰盖模型的初始条件是确定起始时间冰盖的状态。如何得到符合冰盖实际情况的初始状态，即模式的初始化，是一个研究难题。前两年国际专门针对模式初始化开展了模式比较计划initMIP［29-30］。常见的方法有长期的预热模拟（spin-up）和基于数据同化方法的稳态模拟，即融合观测数据和数值模拟得到同化的模式参数和较为接近实际情形的初始状态。同化的参数在预估冰盖未来演化状态中是否仍然适用，是有待商榷的问题。

根据基本物理定律而构造的用来模拟冰盖状态或者预估其未来变化的一套流体力学和热力学偏微分方程组及其求解方案，称之为冰流模式。针对冰盖研发的冰流模式也称为冰盖模式。冰盖模式通过数值求解带有初始条件和边界条件的物理模型，模拟冰盖的速度场、温度场、应力场等主要物理量的时空分布以及冰盖形态和触地线位置的变化，从而有效地理解冰盖动力演化过程。冰盖模式计算的大致流程见图3。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   冰盖动力学数值模式计算流程图

Fig. 3   Flow chart of numerical simulation of ice sheet model


2 中国在南极冰盖数值模拟领域的主要进展

我国早期开展的南极冰盖数值模拟研究主要是一些理想化的冰盖诊断模拟试验，探索冰盖运动变化特征。唐学远等［31］利用基于浅冰近似模型的三维冰盖模式GLIMMER，计算了方形冰盖积分区域的流场特征量，研究了冰盖在大尺度长时间序列条件下对气候变化的反馈，考察了稳定态下冰盖演化的周期性行为。唐学远等［32］采用三维冰盖模式GLIMMER，对南极冰盖开展了理想化冰盖数值试验，给出了冰盖处于稳定状态下的各种演化特征曲线，并将该模式从三维情形简化为二维情形，在耦合温度场条件下模拟了二维冰流的演化，探讨了简化后的二维模型在模拟Dome A地区年代-深度-古积累率关系的应用前景。Zhang等［11］比较了三维full-Stokes冰盖模式求解器FELIX-S［7-9］和冰流模式Elmer/Ice两个不同Stokes冰盖模式对简化的海洋性冰盖的模拟结果，研究表明：对于诊断模拟，使用相同的几何形状和网格时，两个模型给出了相似的结果（差异<2%），原因是两个模型的结构是相似的；对于预测模拟，发现FELIX-S（Elmer/Ice）接地线相对更后退（前进），结果与在诊断试验结果中观察到的微小差异一致，这是由于两个模型中底部边界条件的不同选择造成的［11］。

在过去四十年，我国南极科学考察以中山站、昆仑站、泰山站等为支撑平台，在东南极伊丽莎白公主地和Dome A地区逐步实施了冰芯钻探和强化观测［33-34］。沿中山站—Dome A断面，已逐步建设了星-空-地一体化的业务监测空间网络。基于这些极地考察区域的观测优势，我国在南极内陆冰盖典型区域开展数值模拟研究，取得了一些国际水平的研究成果。

Dome A在昆仑站附近，位于东南极冰盖分冰岭的中心，被认为是获取古老冰芯记录的理想地点［35］。Sun等［36］利用冰流模式Elmer/Ice，基于各向同性和简单的各向异性的冰晶组构的假设条件，对Dome A深冰芯所在区域冰温和冰龄的空间分布进行诊断模拟，模拟结果表明昆仑站可能提供过去60~70万年垂直钻探的高分辨率记录。Zhao等［37］采用冰流模式Elmer/Ice，将Wang等［38］发现的各向异性冰晶组构随深度交替现象进行参数化，并利用雷达等时层的年代信息作为约束条件，重新模拟了Dome A地区的冰温和冰龄空间分布，估算了昆仑站冰芯位置的底部融化速率，指出超过一百万年的古老冰有望在距离昆仑站5~6 km以内、冰底以上200 m深度处找到。Wang等［39］利用内部等时层和Vostok冰芯年代资料，得到Dome A六条等时层位上的深度-年代关系，并将其输入到Dansgaard-Johnsen冰流模式，揭示了Dome A地区和昆仑站深冰钻探点在16万年以来6个不同历史时期的古积累率的时空变化。张良甫等［40］采用基于full-Stokes方程的三维热力-动力耦合冰流模式Elmer/Ice，模拟了甘布尔采夫山脉的速度场和温度场。泰山站是中国于2014年新建立的科考站。Tang等［41］对泰山站的冰川与气象条件开展了综合的分析，并采用一维垂向模型模拟了冰龄和冰底温度。

在过去十几年，我国对东南极的埃默里冰架（Amery Ice Shelf）和普里兹湾（Prydz Bay）开展了长期的观测，对以这两个区域为代表的南极冰盖边缘开展了数值模拟研究。基于高分辨率的海洋-海冰-冰架耦合数值模式，Liu等［42］首次命名了普里兹湾东岸流（PBECC），继而诊断了PBECC的海洋热输送能力及其对Amery冰架底部质量平衡的影响［43］。Cheng等［44］利用自主研发的冰架海洋边界层模式对Amery冰架底部的边界层冰晶羽流过程开展了深入的研究和数值模拟，提高了我国自主研发冰架海洋边界层模式的技术水平。季青原等［45］利用PISM冰盖模式模拟了Amery冰架的速度场，模拟结果总体上与观测结果较为接近。Li等［46］利用Úa冰流模式模拟和分析了Amery冰架在三次崩解事件前后的变化，发现崩解事件对Amery冰架的稳定性影响很小；也对Amery冰架开展了反演和预测试验，模拟了冰架解体极端情景下的冰流加速状况［47］。Sun等［48］基于冰盖模式BISICLES，采取多种冰架底面融化率参数化方案，对南极包含Totten冰川的奥罗拉（Aurora）流域开展动力学模拟和预估研究。

3 中国在北极的冰盖数值模拟领域主要进展

我国在北极地区的冰盖模拟研究主要集中在格陵兰冰盖和斯瓦尔巴群岛、冰岛等地区的冰川。在针对整个格陵兰冰盖的数值模拟方面，Yan等［49］利用冰盖模式SICOPOLIS，探讨了不同现代气候强迫和初始化方案对于现代格陵兰冰盖模拟结果的影响，预估了21世纪格陵兰冰盖物质平衡的变化以及冰盖融化对全球海平面上升的贡献。

表面融水和溢出冰川动力过程的物质输送是造成格陵兰冰盖物质损失的两大主要因素。格陵兰冰盖表面消融产生大量融水，形成冰面水系并输送到冰裂隙、竖井、冰面湖，下渗到冰盖内部和底部，通过促进底部滑动和提升冰温，对冰盖运动速率与稳定性产生影响。格陵兰冰盖溢出冰川末端大多延伸到狭长的峡湾，通过峡湾连通到开阔的海域。海洋环流的变化带动峡湾水体的变化，对冰川的物质平衡和动力过程起到重要的影响。运用数值模拟预估格陵兰冰盖表面径流和溢出冰川的演化是理解和评估格陵兰冰盖未来物质平衡变化的重要方面。

在针对格陵兰冰盖表面径流的模拟研究中，我国在不同的时空尺度上开展了研究。Smith和Yang等［50］利用卫星遥感与无人机影像，提取了一个典型集水区的冰面水系，结合实地测量的冰面河出口位置径流量，提出了一种测算冰面径流量的新方法，研究结果表明，目前常用的区域气候模式过高地估算了冰面实际径流量，冰面融水输送的方式和效率对于冰盖物质平衡以及动力过程具有重要影响。Moore等［51］利用冰盖表面能量-物质平衡模型，研究平流层注入硫酸盐气溶胶的地球工程试验G4对格陵兰冰盖表面径流的影响，模拟结果显示在G4情景下冰盖表面径流比RCP4.5情景下减少了20%。Yue等［52］利用度日模型，开展了21世纪和22世纪不同气候情景下格陵兰冰盖表面径流的模拟和预估研究。

在溢出冰川动力过程的数值模拟方面，我国主要对格陵兰冰盖最快的溢出冰川——雅各布港冰川进行了数值模拟。Guo和Zhao等［53］采用冰盖模式BISICLES和同化的参数方案，重现了该冰川在2004—2013年间由于海水增暖造成的冰流速度加速等行为过程，在季节尺度上模拟的冰流速度、冰架末端位置变化、触地线位置变化与观测值具有很好的一致性，并且在RCP4.5情景下预估了冰川的动力学演化。

除了格陵兰冰盖，我国科研工作者依托黄河站为支撑平台，对北极斯瓦尔巴群岛的冰川等开展了实地观测和数值模拟有效结合的研究。Ai等［54］利用Elmer/Ice冰流模式对北极黄河站附近的Austre Lovénbreen冰川开展了数值模拟，发现该冰川除了以前实测发现的快速冰流区，还有另一个速度更快的区域。基于该模拟发现，2016年起开展了高精度GPS加密测量，证实了该发现，新发现的快速冰流区比以前测得的快速流速区流速要快8%左右。此外，以历史观测记录和Elmer/Ice冰流模式为参照，艾松涛等［55］在忽略动力学因素的条件下，利用ArcGIS软件在冰川轮廓内的每个像素处计算冰面高程变化，进而模拟山地冰川面积、冰厚和体积的长期变化，并探讨了该方法的可行性和适用条件。在2015/2016年的北极异常升温背景下，黄河站Austre Lovénbreen冰川和Pedersenbreen冰川的流速都有较大增加。Wang等［56］通过模拟发现，底部滑动是冰流增速的主要影响因素，直接原因是冰下水量（subglacial water）的增加。随着冰川融水渗透到冰川内部，还可能出现冬季缓慢释放至冰川末端形成冰湖的现象［57］。通过冰川生命周期的模拟，斯瓦尔巴群岛上类似Austre Lovénbreen冰川的小冰川很可能在一百余年后彻底消失［58］。另外，Yue等［59］应用冰盖模式PISM，模拟和预估了冰岛冰帽在21世纪物质平衡和动力过程的演化。

4 总结与展望

回顾过去十几年，我国在极地冰盖数值模拟领域的研究取得了明显的进步。我国已经自主研发了三个冰盖模式，并在持续改进。同时，在极地冰盖及其典型区域上应用冰盖模式模拟，得到了一些重要的科学认识，取得了具备国际水平的研究成果。经过数十年的极地科考活动，我国在极地冰盖的监测方面形成了强化监测优势区域，其中包括南极的Dome A、Amery冰架、伊丽莎白公主地和北极的黄河站等，积累了大量宝贵的数据资源。近年来在格陵兰冰盖典型区域也开展了空-天-地一体化监测［60］。这些都为研究冰盖变化关键过程提供了良好的数据平台。

与此同时，我国在极地冰盖数值模拟领域也面临着一些问题，制约着当前和未来的工作，其中包括：研究人员和人才很少，没有形成实力较强的团队；拥有自主知识产权的冰盖模式在研发上处于起步阶段，需要持续投入；我国在国际冰盖模式比较计划中的参与度不足；研究成果的数量和质量尽管有很大进步，但是涉及入海冰川或冰架崩解和消融的模拟研究较少；国际上已有研究团队探索了大气环流模式与冰盖模式的耦合［61-64］以及海洋模式与冰盖模式的耦合［65-67］，提高模式模拟的精度，而我国目前鲜有海-冰-气耦合模式方面的研究。

针对这些问题，需要提高对冰盖模式模拟研究领域的重视程度，持续加大投入力度，加强国内研究院所的交流合作，有效结合实地观测、遥感监测和数值模拟，形成研究合力；加强与大气、海洋等其他圈层领域互动，推进海-冰-气多圈层耦合和相互作用研究，提升极地冰盖模拟在气候变化科学领域的重要性；在我国传统观测优势区域，通过观测、模拟等综合研究手段开展联合攻关，有望在冰盖关键动力过程和机制的科学认识上有所突破，提升冰盖模式对物理过程的刻画能力，产出具有中国特色的世界先进水平的研究成果。展望未来，可能需要在以下几个方面加强研究力量：

（1）充分利用多源数据和区域气候模式，构建长时间高分辨率时空连续的冰盖表面物质平衡，准确刻画区域尺度空间变化特征；开展格陵兰冰盖表面融水和冰面水系的观测和模拟，关联冰面-冰下融水，探索冰面融水对冰盖运动的影响。

（2）开展冰盖典型流域的动力学模拟，研究海洋强迫对入海冰川末端或者冰架的作用机理，改进冰架消融和崩解在冰盖模式中的实现，探索入海冰川或者冰架与海洋之间的相互作用和耦合模拟，加深对关键动力学机制的理解。

（3）研究冰盖冰下水系分布、冰下水文过程及其对冰盖不稳定性的影响，开展冰下水文过程模拟，在冰盖模式中考虑冰下水文过程的影响。

参考文献

Fox-Kemper B， Hewitt H T， Xiao C， et al.

Ocean， cryosphere and sea level change

［M］// Climate change 2021： the physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press， 2021.

[本文引用: 1]

Rignot E， Mouginot J， Scheuchl B， et al.

Inaugural article： four decades of Antarctic Ice Sheet mass balance from 1979—2017

［J］. Proceedings of the National Academy of Sciences of the United States of America， 2019， 116（4）： 1095-1103.

[本文引用: 1]

Shen Q， Wang H， Shum C， et al.

Recent high-resolution Antarctic ice velocity maps reveal increased mass loss in Wilkes Land， East Antarctica

［J］. Scientific reports， 2018， 8（1）： 1-8.

[本文引用: 1]

Seroussi H， Nowicki S， Payne A， et al.

ISMIP6 Antarctica： a multi-model ensemble of the Antarctic Ice Sheet evolution over the 21st century

［J］. The Cryosphere， 2020， 14（9）： 3033-3070.

[本文引用: 3]

Wang Y， Zhang T， Xiao C， et al.

A two-dimensional， higher-order， enthalpy-based thermomechanical ice flow model for mountain glaciers and its benchmark experiments

［J］. Computers & Geosciences， 2020， 141： 104526.

[本文引用: 1]

Zhang H， Ju L， Gunzburger M， et al.

Coupled models and parallel simulations for three-dimensional Full-Stokes ice sheet modeling

［J］. Numerical Mathematics： Theory， Methods and Application， 2011， 3： 396-418.

[本文引用: 2]

Leng W， Ju L， Gunzburger M， et al.

A parallel high-order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments

［J］. Journal of Geophysical Research： Earth Surface， 2012， 117： F01001.

[本文引用: 3]

Leng W， Ju L， Gunzburger M， et al.

Manufactured solutions and the verification of three-dimensional Stokes ice-sheet models

［J］. The Cryosphere， 2013（7）： 19-29.

[本文引用: 1]

Leng W， Ju L， Gunzburger M， et al.

A parallel computational model for three-dimensional， thermo-mechanical Stokes flow simulations of glaciers and ice sheets

［J］. Communications in Computational Physics， 2014， 16（4）： 1056-1080.

[本文引用: 3]

Leng W， Ju L， Xie Y， et al.

Finite element three dimensional Stokes ice sheet dynamics model with enhanced local mass conservation

［J］. Journal of Computational Physics， 2014， 274： 299-311.

[本文引用: 1]

Zhang T， Price S， Ju L， et al.

A comparison of two Stokes ice sheet models applied to the Marine Ice Sheet Model Intercomparison Project for plan view models （MISMIP3d）

［J］. The Cryosphere， 2017， 11， 179-190.

[本文引用: 3]

Bell R， Ferraccioli F， Creyts T， et al.

Widespread persistent thickening of the east Antarctic Ice Sheet by freezing from the base

［J］. Science， 2011， 331（6024）： 1592-1595.

[本文引用: 1]

Siegert M， Ross N， Le Brocq A.

Recent advances in understanding Antarctic subglacial lakes and hydrology

［J］. Philosophical Transactions of the Royal Society A： Mathematical， Physical and Engineering Sciences， 2016， 374（2059）： 20140306.

[本文引用: 1]

Tang Xueyuan， Zhang Zhanhai， Sun Bo.

Progress and prospect in numerical modelling of the Antarctic Ice Sheet

［J］. Chinese journal of polar research， 2006， 18（4）： 290-300.

[本文引用: 1]

唐学远， 张占海， 孙波.

南极冰盖数值模拟研究进展与展望

［J］. 极地研究， 2006， 18（4）： 290-300.

[本文引用: 1]

Tang Xueyuan， Sun Bo， Li Yuansheng， et al.

Some recent progress of Antarctic Ice Sheet research

［J］. Advances in Earth Science， 2009， 24（11）： 1210-1218.

[本文引用: 1]

唐学远， 孙波， 李院生， 等.

南极冰盖研究最新进展

［J］. 地球科学进展， 2009， 24（11）： 1210-1218.

[本文引用: 1]

Greve R， Blatter H. Dynamics of Ice Sheets and Glaciers［M］. Springer Science & Business Media， Berlin， Germany， 2009.

[本文引用: 1]

Schoof C， Hindmarsh R.

Thin-film flows with wall slip： an asymptotic analysis of higher order glacier flow models

［J］. Quarterly Journal of Mechanics and Applied Mathematics， 2010， 63（1）： 73-114.

[本文引用: 1]

Hulbe C and MacAyeal D R.

A new numerical model of coupled inland ice sheet， ice stream， and ice shelf flow and its application to the West Antarctic Ice Sheet

［J］. Journal of Geophysical Research， 1999， 104（B11）： 25349-25366.

[本文引用: 1]

Pollard D and DeConto R M.

Description of a hybrid ice sheet shelf model， and application to Antarctica

［J］. Geoscientific Model Development， 2012， 5： 1273-1295.

Pattyn F.

Sea-level response to melting of Antarctic ice shelves on multi-centennial timescales with the fast Elementary Thermomechanical Ice Sheet model （f.ETISh v1.0）

［J］. The Cryosphere， 2017， 11： 1851-1878.

[本文引用: 1]

Gillet-Chaulet F， Gagliardini O， Seddik H， et al.

Greenland Ice Sheet contribution to sea-level rise from a new-generation ice-sheet model

［J］. The Cryosphere， 2012， 6： 1561-1576.

[本文引用: 1]

Cornford S L， Martin D F， Payne A J， et al.

Century scale simulations of the response of the West Antarctic Ice Sheet to a warming climate

［J］. The Cryosphere， 2015， 9： 1579-1600.

[本文引用: 1]

Nicholls K， Corr H， Stewart C.

A ground-based radar for measuring vertical strain rates and time-varying basal melt rates in ice sheets and shelves

［J］. Journal of Glaciology， 2015， 61（230）： 1079-1087.

[本文引用: 1]

Zhao C， Gladstone R M， Warner R C， et al.

Basal friction of Fleming Glacier， Antarctica-Part 1： sensitivity of inversion to temperature and bedrock uncertainty

［J］. The Cryosphere， 2018， 12（8）： 2637-2652.

[本文引用: 1]

Gong Y， Zwinger T， Cornford S， et al.

Importance of basal boundary conditions in transient simulations： case study of a surging marine-terminating glacier on Austfonna， Svalbard

［J］. Journal of Glaciology， 2016， 63（237）： 106-117.

[本文引用: 1]

Martin M A， Winkelmann R， Haseloff M， et al.

The Potsdam Parallel Ice Sheet Model （PISM-PIK）–Part 2： dynamic equilibrium simulation of the Antarctic Ice Sheet

［J］. The Cryosphere， 2011， 5： 727-740.

[本文引用: 1]

DeConto R M， Pollard D.

Contribution of Antarctica to past and future sea-level rise

［J］. Nature， 2016， 531： 591-597.

Goelzer H， Nowicki S， Payne A.

The future sea-level contribution of the Greenland Ice Sheet： a multi-model ensemble study of ISMIP6

［J］. The Cryosphere， 2020， 14（9）： 3071-3096.

[本文引用: 1]

Goelzer H， Nowicki S， Edwards T， et al.

Design and results of the ice sheet model initialisation experiments initMIP-Greenland： an ISMIP6 intercomparison

［J］. The Cryosphere， 2018， 12： 1433-1460.

[本文引用: 1]

Seroussi H， Nowicki S， Simon E， et al.

initMIP-Antarctica： an ice sheet model initialization experiment of ISMIP6

［J］. The Cryosphere， 2019， 13（5）： 1441-1471.

[本文引用: 1]

Tang Xueyuan， Sun Bo， Li Yuansheng， et al.

GLIMMER Antarctic Ice Sheet model， an experimental research of moving boundary condition

［J］. Journal of Glaciology and Geocryology， 2007， 29（6）： 905-913.

[本文引用: 1]

唐学远， 孙波， 张占海， 等.

南极冰盖GLIMMER模式移动边界试验研究

［J］. 冰川冻土， 2007， 29（6）： 905-913.

[本文引用: 1]

Tang Xueyuan， Sun Bo， Li Yuansheng， et al.

GLIMMER Antarctic Ice Sheet model， model simplification of the two-dimensional ice stream

［J］. Progress in Natural Science， 2007， 17（4）： 480-487.

[本文引用: 1]

唐学远， 张占海， 孙波， 等.

南极冰盖GLIMMER模式试验及其对二维冰流情形的模型简化

［J］. 自然科学进展， 2007， 17（4）： 480-487.

[本文引用: 1]

Sun Bo， Siegert M， Mudd S， et al.

The Gamburtsev mountains and the origin and early evolution of the Antarctic Ice Sheet

［J］. Nature， 2009， 459（7247）： 690-693.

[本文引用: 1]

Cui X， Jeofry H， Greenbaum J， et al.

Bed topography of Princess Elizabeth Land in East Antarctica

［J］. Earth System Science Data， 2020， 12（4）： 2765-2774.

[本文引用: 1]

Tang Xueyuan， Sun Bo， Li Yuansheng， et al.

Review of the glaciological research progress and future development of deep ice core plan at Dome A， East Antarctica

［J］. Chinses Journal of Polar Research， 2012， 24（1）： 77-86.

[本文引用: 1]

唐学远 孙波， 李院生， 等.

冰穹A冰川学研究进展及深冰芯计划展望

［J］. 极地研究， 2012， 24（1）： 77-86.

[本文引用: 1]

Sun B， Moore J， Zwinger T， et al.

How old is the ice beneath Dome A， Antarctica？

［J］. The Cryosphere， 2014， 8（3）： 1121-1128.

[本文引用: 1]

Zhao L， Moore J， Sun Bo， et al.

Where is the 1-million-year-old ice at Dome A？

［J］. The Cryosphere， 2018， 12（5）： 1651-1663.

[本文引用: 1]

Wang B， Sun B， Martin C， et al.

Summit of the East Antarctic Ice Sheet underlain by thick ice-crystal fabric layers linked to glacial-interglacial environmental change

［J］. Geological Society， London， Special Publications， 2018， 461（1）： 131-143.

[本文引用: 1]

Wang T， Sun B， Tang X， et al.

Spatio-temporal variability of past accumulation rates inferred from isochronous layers at Dome A， East Antarctica

［J］. Annals of Glaciology， 2016， 57（73）： 87-93.

[本文引用: 1]

Zhang Liangpu， Tang Xueyuan， Yang Shuhu， et al.

Numerical simulations of east Antarctic Ice Sheet based on the Elmer/Ice model

［J］. Chinese Journal of Polar Research， 2017， 29（3）： 390-398.

[本文引用: 1]

张良甫， 唐学远， 杨树瑚， 等.

基于Elmer/Ice冰盖模式的南极Gamburtsev山脉Lambert冰流区域的数值模拟研究

［J］. 极地研究， 2017， 29（3）： 390-398.

[本文引用: 1]

Tang X， Guo J， Dou Y， et al.

Glaciological and meteorological conditions at the Chinese Taishan Station， East Antarctica

［J］. Frontiers in Earth Science， 2020， 8： 250.

[本文引用: 1]

Liu C， Wang Z， Cheng C， et al.

Modeling modified Circumpolar Deep Water intrusions onto the Prydz Bay continental shelf， East Antarctica

［J］. Journal of Geophysical Research： Oceans， 2017， 122（7）： 5198-5217.

[本文引用: 1]

Liu C， Wang Z， Cheng C， et al.

On the modified Circumpolar Deep Water upwelling over the Four Ladies Bank in Prydz Bay， East Antarctica

［J］. Journal of Geophysical Research： Oceans， 2018， 123（11）： 7819-7838.

[本文引用: 1]

Cheng C， Jenkins A， Wang Z， et al.

Modeling the vertical structure of the ice shelf-ocean boundary current under supercooled condition with suspended frazil ice processes： a case study underneath the Amery Ice Shelf， East Antarctica

［J］. Ocean Modelling， 2020， 156： 101712.

[本文引用: 1]

Ji Qingyuan， Wang Bangbing， Sun Bo.

Applicability PISM for velocity analysis of the Amery ice shelf， East Antarctica

［J］. Chinese journal of Polar Research， 2015， 27（3）： 229-236.

[本文引用: 1]

季青原， 王帮兵， 孙波.

PISM冰盖模式对Amery冰架流速场模拟的适用性

［J］. 极地研究， 2015， 27（3）： 229-236.

[本文引用: 1]

Li T， Liu Y， Cheng X.

Recent and imminent calving events do little to impair Amery ice shelf's stability

［J］. Acta Oceanologica Sinica， 2020， 39（5）： 168-170.

[本文引用: 1]

Li Teng， Chen Zhuoqi， Li Huilin， et al.

A new-generation ice flow model Úa and its application in Amery Ice Shelf， Antarctica

［J］. Journal of Glaciology and Geocryology， 2020， 42（1）： 254-264.

[本文引用: 1]

李腾， 陈卓奇， 李慧林， 等.

新一代冰流模式乌阿及其在南极埃默里冰架的应用

［J］. 冰川冻土， 2020， 42（1）： 254-264.

[本文引用: 1]

Sun S， Cornford S， Gwyther D， et al.

Impact of ocean forcing on the Aurora Basin in the 21st and 22nd centuries

［J］. Annals of Glaciology， 2016， 57（73）： 79-86.

[本文引用: 1]

Yan Q， Zhang Z， Gao Y， et al.

Sensitivity of the modeled present-day Greenland Ice Sheet to climatic forcing and spin-up methods and its influence on future sea level projections

［J］. Journal of Geophysical Research Earth Surface， 2013， 118（4）： 2174-2189.

[本文引用: 1]

Smith L， Yang Kang， Pitcher L， et al.

Direct measurements of meltwater runoff on the Greenland Ice Sheet surface

［J］. Proceedings of the National Academy of Sciences of the United States of America， 2017， 114（50）： E10622-E10631.

[本文引用: 1]

Moore J C， Yue C， Zhao L， et al.

Greenland Ice Sheet response to stratospheric aerosol injection geoengineering

［J］. Earth’s Future， 2019， 7.

[本文引用: 1]

Yue C， Zhao L， Wolovick M， et al.

Greenland Ice Sheet surface runoff projections to 2200 using degree-day methods

［J］. Atmosphere， 2021， 12（12）： 1569.

[本文引用: 1]

Guo X， Zhao L， Gladstone R， et al.

Simulated retreat of Jakobshavn Isbræ during the 21st century

［J］. The Cryosphere， 2019， 13（11）： 3139-3153.

[本文引用: 1]

Ai S， Ding X， An J， et al.

Discovery of the fastest ice flow along the central flow line of Austre Lovénbreen， a poly-thermal valley glacier in Svalbard

［J］. Remote Sensing， 2019， 11（12）： 1488.

[本文引用: 1]

Ai Songtao， Yan Boya， Wang Zemin， et al.

Projection of long-term changes of mountain glacier based on GIS grid operation method

［J］. Polar Research， 2019， 31（3）： 267-275.

[本文引用: 1]

艾松涛， 严博雅， 王泽民， 等.

基于GIS栅格操作方法模拟山地冰川长期变化

［J］. 极地研究， 2019， 31（3）： 267-275.

[本文引用: 1]

Wang Z， Yan B， Ai S， et al.

Quantitative analysis of Arctic ice flow acceleration with increasing temperature

［J］. Acta Oceanologica Sinica， 2021， 40（1）： 22-32.

[本文引用: 1]

Ai S， Yan B， Wang Z， Yan M.

A decadal record of inter-annual surface ice flow from Pedersenbreen， Svalbard （2005-15）

［J］. Polar Science， 2019， 22： 100485.

[本文引用: 1]

Wang Z， Lin G， Ai S.

How long will an Arctic mountain glacier survive？ A case study of Austre Lovénbreen， Svalbard

［J］. Polar Research， 2019， 38： 3519.

[本文引用: 1]

Yue C， Schmidt L， Zhao L， et al.

Insensitivity of mass loss of Icelandic Vatnajökull ice cap to solar geoengineering

［J］. The Cryosphere Discussions， 2021： 1-20.

[本文引用: 1]

Xiao Cunde， Chen Zhuoqi， Jiang Liming， et al.

A study of monitoring， simulation and climate impact of Greenland Ice Sheet

［J］. Advances in Earth Science， 2019， 34（8）： 781-786.

[本文引用: 1]

效存德， 陈卓奇， 江利明， 等.

格陵兰冰盖监测、模拟及气候影响研究

［J］. 地球科学进展， 2019， 34（8）： 781-786.

[本文引用: 1]

Huybrechts P， Janssens I， Poncin C， et al.

The response of the Greenland Ice Sheet to climate changes in the 21st century by interactive coupling of an AOGCM with a thermomechanical ice-sheet model

［J］. Annual Glaciology， 2002， 35： 409-415.

[本文引用: 1]

Helsen M M， van de Wal R S W， van den Broeke M R， et al.

Coupling of climate models and ice sheet models by surface mass balance gradients： application to the Greenland Ice Sheet

［J］. The Cryosphere， 2012， 6： 255-272.

Vizcaino M， Mikolajewicz U， Ziemen F， et al.

Coupled simulations of Greenland Ice Sheet and climate change up to A.D. 2300

［J］. Geophysical Research Letters， 2015， 42， 3927-3935.

Le clec’h S， Charbit S， Quiquet A， et al.

Assessment of the Greenland Ice Sheet-atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model

［J］. The Cryosphere， 2019， 13： 373-395.

[本文引用: 1]

Seroussi H， Nakayama Y， Larour E， et al.

Continued retreat of Thwaites Glacier， West Antarctica， controlled by bed topography and ocean circulation

［J］. Geophysical Research Letters， 2017， 44： 6191-6199.

[本文引用: 1]

Beckmann J， Perrette M， Beyer S， et al.

Modeling the response of Greenland outlet glaciers to global warming using a coupled flow line-plume model

［J］. The Cryosphere， 2019， 13， 2281-2301.

Favier L， Jourdain N C， Jenkins A， et al.

Assessment of sub-shelf melting parameterisations using the ocean-ice-sheet coupled model NEMO （v3.6）–Elmer/Ice （v8.3）

［J］. Geoscientific Model Development， 2019， 12： 2255-2283.

[本文引用: 1]

/

〈 〉

陇ICP备2021001824号-10
　 甘公网安备 62010202000676号
版权所有 © 《冰川冻土》编辑部
地址：兰州市天水中路8号　
邮编：730000　电话：0931-8260767　
E-mail：edjgg@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司