冰川冻土, 2022, 44(3): 872-884 doi: 10.7522/j.issn.1000-0240.2022.0084

冰冻圈与全球变化

2000—2020年格陵兰冰盖夏季表面温度变化及其对物质平衡的影响

方振祥,1,2, 王宁练,1,2,3, 李想1,2, 张玉杰1,2, 邰雪楠1,2

1.陕西省地表系统与环境承载力重点实验室,陕西 西安 710127

2.西北大学 城市与环境学院/地表系统与灾害研究院,陕西 西安 710127

3.中国科学院 青藏高原研究所,北京 100101

Summer surface temperature change of Greenland Ice Sheet from 2000 to 2020 and its effect on mass balance

FANG Zhenxiang,1,2, WANG Ninglian,1,2,3, LI Xiang1,2, ZHANG Yujie1,2, TAI Xuenan1,2

1.Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity,Xi’an 710127,China

2.Institute of Earth Surface System and Hazards,College of Urban and Environmental Sciences,Northwest University,Xi’an 710127,China

3.Institute of Tibetan Plateau Research,Chinese Academy of Sciences,Beijing 100101,China

通讯作者: 王宁练,教授,主要从事冰冻圈与全球变化研究. E-mail: nlwang@nwu.edu.cn

收稿日期: 2022-04-13   修回日期: 2022-06-14  

基金资助: 中国科学院战略先导科技专项(A类).  XDA19070302
国家自然科学基金项目.  41971088

Received: 2022-04-13   Revised: 2022-06-14  

作者简介 About authors

方振祥,硕士研究生,主要从事冰川温度研究.E-mail:fzxdhh@163.com , E-mail:fzxdhh@163.com

摘要

基于MODIS温度产品,着重分析了2000—2020年格陵兰冰盖夏季表面温度和表面融化范围的年际变化趋势;联合IMBIE(冰盖物质平衡对比实验)数据分析表面温度对于冰盖物质平衡的影响;进一步讨论了大气环流对于格陵兰冰盖表面温度变化的影响。结果表明:格陵兰冰盖夏季表面温度和融化范围趋势较为一致,2000年初期呈现出显著的上升趋势,2012年达到峰值,随后波动下降;整个研究阶段北部区域是增温速率最大的区域,高于其他任何区域两倍,东南部和西南部是温度最高的区域却具有最小的增长率;格陵兰冰盖夏季表面温度、融化范围以及物质平衡之间都具有显著的相关性,同时格陵兰冰盖夏季表面温度每上升1 ℃,会导致其物质损失增加74.29 Gt·a-1;最后,经过对北大西洋涛动(NAO)和格陵兰阻塞指数(GBI)指数的分析得到,格陵兰冰盖夏季表面温度受到GBI的影响要强于NAO的影响,冰盖夏季表面温度和NAO呈现出负相关(r=-0.64,P<0.05),和GBI呈现出正相关(r=0.77,P<0.05)。

关键词: 表面温度 ; 融化 ; 物质平衡 ; 北大西洋涛动 ; 格陵兰阻塞指数 ; 格陵兰冰盖

Abstract

The monitoring of the Greenland Ice Sheet is of great significance to global and Arctic climate change. As the largest Ice Sheet in the Northern Hemisphere, the Greenland Ice Sheet is critical to the climate response and has always been the focus of global attention. Based on the MODIS temperature product data, this paper analyzes the interannual trend of the Greenland Ice Sheet summer surface temperature and surface melting range from 2000 to 2020. The effect of surface temperature on the mass balance of the Greenland Ice Sheet is analyzed with IMBIE data, and the effect of atmospheric circulation on the surface temperature of the Greenland Ice Sheet is discussed. The results show that the trend of summer surface temperature and melting range of the Greenland Ice Sheet is consistent, showing a positive upward trend from 2000 to 2012, reaching a peak in 2012, and then showing a downward trend of fluctuation. The northern part of the Greenland Ice Sheet had the highest rate of warming during the study period, more than twice as fast as any other region, the south-east and south-west had the highest temperatures but the lowest growth rates; And there was a strong correlation between the Greenland Ice Sheet’s summer surface temperature, the extent of melting, and the mass balance, at the same time, an increase of 1 °C in the summer surface temperature of the Greenland Ice Sheet will lead to an increase of 74.29 Gt·a-1 in mass loss. Finally, the NAO (North Atlantic Oscillation) and GBI (Greenland blocking index) indices are analyzed, the influence of GBI on the surface temperature of Greenland Ice Sheet in summer is stronger than that of NAO. Greenland Ice Sheet surface temperature is negatively correlated with NAO (r=-0.64, P<0.05) and positively correlated with GBI (r=0.77, P<0.05).

Keywords: Surface Temperature ; Melt ; Mass balance ; NAO ; GBI ; Greenland Ice Sheet

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

本文引用格式

方振祥, 王宁练, 李想, 张玉杰, 邰雪楠. 2000—2020年格陵兰冰盖夏季表面温度变化及其对物质平衡的影响[J]. 冰川冻土, 2022, 44(3): 872-884 doi:10.7522/j.issn.1000-0240.2022.0084

FANG Zhenxiang, WANG Ninglian, LI Xiang, ZHANG Yujie, TAI Xuenan. Summer surface temperature change of Greenland Ice Sheet from 2000 to 2020 and its effect on mass balance[J]. Journal of Glaciology and Geocryology, 2022, 44(3): 872-884 doi:10.7522/j.issn.1000-0240.2022.0084

0 引言

IPCC(政府间气候变化专门委员会)第六次评估报告提出,全球变暖背景下,全球范围内几乎所有的冰川都处于退缩状态1。极地冰盖对气候变化的响应及其变化对海平面与大洋洋流的影响,长期以来都是全球变化研究关注的焦点问题2-4。自20世纪90年代以来,诸多研究已经证实格陵兰冰盖发生着显著的融化以及物质损失5-8,且融化表现出加速的趋势9-12。Zwally等13对比格陵兰冰盖1992—2002年和2003—2007年两个时期的物质平衡,结果表明冰盖物质损失显著增加;陈国栋等14利用ICESat数据,揭示出2000—2009年期间格陵兰冰盖高程和体积的减少趋势;卢飞等15基于GRACE数据反演了近年来格陵兰冰盖物质变化,证实格陵兰冰盖加速融化且消融区域主要集中在南部和西北区域;进一步研究表明16,假如格陵兰冰盖全部融化,将会导致全球海平面上升7.4 m,会对海岸经济社会造成严重的影响17。因此,对于格陵兰冰盖表面温度、融化范围以及物质平衡的研究对于认识格陵兰冰盖变化及其机理具有重要意义。

格陵兰冰盖表面温度(Land Surface Temperature, LST)对其辐射收支与物质平衡变化具有重要的表征意义12。相关模型预测18-20,格陵兰冰盖表面温度每上升1 ℃会产生20%~50%融冰,且夏季温度的升高会使融化的速度加倍,从而导致更多的物质损失。目前对于格陵兰冰盖温度研究的数据主要来自于自动气象站、区域气候模型以及遥感21。然而,格陵兰冰盖上气象观测点很少,大多不均匀地分布在海岸地区,对于整个格陵兰冰盖的表面温度的测量很难获取,且很难进行长时间序列的观测,同时单点测量的代表性较差,不能完整反应格陵兰冰盖表面温度的变化情况。区域气候模型往往过于粗糙不能捕捉局地尺度气象过程的空间变异性22。因此,利用遥感数据成为近年来获取格陵兰冰盖表面温度的研究趋势。自20世纪70年代起,热红外卫星遥感因具有长时间序列、多波段、大尺度等优点,已成为区域乃至全球范围内解决实地观测资料不足的重要的手段23-25。热红外卫星遥感观测不仅可以重复观测,而且可以大范围的获取地表温度信息,具有较高数据质量及时空分辨率的MODIS(中分辨率成像光谱仪)温度产品已经广泛应用于极地地区的地表温度时空变化研究1226-28。因此,利用高时空分辨率卫星数据可对格陵兰冰盖表面温度的时空变化特征进行系统研究。

准确监测格陵兰冰盖表面温度的时空变化以及物质平衡之间的响应关系,不仅可以更好地认识到格陵兰冰盖的变化及其机理,对于评估区域气候乃至全球尺度的气候变化也具有重要的科学意义和现实意义。最近研究表明29-31,由于增强的表面融化以及异常的气候条件,21世纪格陵兰冰盖物质相比20世纪损失更多。为了更好地理解本世纪格陵兰冰盖表面融化、物质平衡以及表面过程,基于MODIS温度产品数据,首先分析了2000—2020年夏季格陵兰冰盖整体及六个主要部分的表面温度以及表面融化范围的年际变化趋势;结合IMBIE(Ice Sheet Mass Balance Inter-comparison Exercise)物质平衡数据分析了冰盖表面温度对于物质平衡的影响;最后结合大气环流指数,分析了2000—2020年夏季格陵兰冰盖表面温度变化的原因。

1 研究区概况

格陵兰岛(59.46~83.39° N,73.08°~11.39° W)是世界上面积最大的岛屿,面积约为216万平方千米,约83.7%被冰雪覆盖。与南极冰盖相比,格陵兰冰盖对于气候的响应更为敏感。图1是格陵兰冰盖的分布范围,Rignot等32-33在开展冰盖物质平衡对比实验(IMBIE 2016)研究时曾将格陵兰冰盖划分为北部(NO)、东北部(NE)、东南部(SE)、西南部(SW)、中西部(CW)以及西北部(NW)六个区域,如图1所示。关于格陵兰冰盖夏季表面温度变化的相关分析也按照这六个区域进行分析。

图1

图1   格陵兰岛位置示意图

Fig. 1   Map of the location of Greenland. contour lines are also represented


2 数据和方法

2.1 数据获取

2.1.1 MODIS温度产品

中分辨率成像光谱仪MODIS搭载于Terra和Aqua卫星上,Terra与Aqua卫星数据最早分别始于2000年和2002年。MODIS数据产品极为丰富,本文使用的MODIS温度产品来源于美国国家航空航天局NASA的数据发布服务器(https://ladsweb.modaps.eosdis.nasa.gov/),已经广泛应用于格陵兰冰盖表面温度的研究,且具有较好的适用性和准确性:Hall等34比较了MODIS LST与格陵兰冰盖现场数据,表明冰雪地区产品误差约为2 ℃,且与ASTER等温度产品具有较好的一致性;Koenig等35的分析表明,在格陵兰冰盖MOD11资料与实测数据的RMSE为3.1 ℃;Stroeve等36利用MOD11A1数据与格陵兰自动气象站自动收集的原始数据进行比较,均方根误差为0.07,相关系数高达0.79。MOD11A2是等级为3的MODIS地表温度和发射率产品,时间分辨率为8天,由存储在1公里正弦网格上的每日1公里LST产品(MOD11A1)取8天期间晴空LST的平均值合成(有较好的代表性且可以很好地减轻云量造成的数据减少)。本文使用其温度数据集(LST_DAY)和质量控制集(QC_DAY),选用白天数据,来分析2000—2020年期间格陵兰冰盖的夏季(6—8月)表面温度的变化。

2.1.2 ERA5气温数据

ERA5是第五代ECMWF(欧洲中期天气预报中心)大气在分析全球气候数据,提供了1950年至今的大气、陆地和海洋气候变量的估计值。本文使用的数据为ERA5每日2 m平均气温数据(2000—2019年)。格陵兰冰盖表面温度和近地表气温具有较高的相关性,因此对比其距平变化,可以有效验证MODIS温度数据的可用性,进一步支撑我们的研究结论。

2.1.3 物质平衡数据

IMBIE成立于2011年,旨在协调冰盖物质平衡的卫星测量结果。IMBIE主要是由欧洲空间局(ESA)以及美国航空航天局(NASA)支持的国际科学家之间的合作,主要是作为对IPCC的贡献,其提高了冰盖物质平衡测量和相关全球海平面贡献的可信度。本文使用的数据为IMBIE 2019年发布的关于格陵兰冰盖物质平衡的最新报告,使用的数据为2000—2018年格陵兰冰盖的物质平衡数据和相应误差数据 (IMBIE 2019 Greenland Dataset;http://imbie.org/data-downloads/)。

2.1.4 大气环流指数

长期以来,格陵兰气候受到NAO和GBI的高度影响,用北大西洋涛动(NAO)和格陵兰阻塞指数(GBI)来分析大气环流对格陵兰冰盖表面温度的影响37-38。NAO每月数据来自于(https://www.cpc.ncep.noaa.gov/products/precip/CWlink/pna/nao.shtml39,GBI每月数据来自于美国国家海洋和大气管理局地球系统研究实验室(US NOAA’s ESRL)(https://www.esrl.noaa.gov/psd/gcos_wgsp/Timeseries/GBI_UL/40以及Hanna41

2.2 研究方法
2.2.1 数据预处理

使用了可完整覆盖格陵兰冰盖的2000—2020年6—8月MOD11A2数据,条幅号为H15V02、H16V00、H16V01、H16V02、H17V00、H17V01、H17V02,大约共有1 764幅。数据预处理主要包括影像拼接、投影转换以及裁剪。采用MODIS Reprojection Tools (MRT)工具将不同轨道号数据进行拼接,提取LST_DAY和QC_DAY数据,然后将正弦投影转换为横轴墨卡托投影,最后对拼接后的区域进行裁剪得到研究区域的完整图幅。MODIS温度产品的每一个网格都有一个对应的质量控制值(quality control),为0、1、2、3四个等级,表示误差值分为0~1 K、1~2 K、2~3 K以及>3 K(K为开尔文,热力学温度)。本研究使用等级为0和1的像元,剔除等级为2和3的误差较大的像元以减少数据本身带来的误差影响,即MODIS温度产品自身评估误差被认为>2 K的像元被剔除并用Null值替代,以确保更接近晴空条件下的冰盖实际表面温度。温度数据集有效数据(DN)范围为7 500~65 535,将经过质量筛选的数据转化为摄氏度(LST=0.02×DN-273.15),由所有的温度数据计算得到平均夏季表面温度,进行分析计算。夏季平均温度是由12个影像平均得到,月平均温度由相应的4个影像计算得到。

Välisuo等42总结了15个关于格陵兰冰盖表面融化研究的定义,主要使用融化程度、融化天数以及持续时间等来定义融化,且利用MODIS温度产品对格陵兰冰盖融化范围的评估已经广泛适用。相关研究表明,MODIS地表温度产品可能会产生一些冷偏差2735,Hall等12将冰盖表面温度大于等于 -1 ℃定义为融化发生,同样,也将-1 ℃作为评估格陵兰冰盖表面是否发生融化的阈值,并取每年夏季中格陵兰表面温度高于-1 ℃的最大范围作为当年冰盖表面的最大融化范围。值得注意的是,基于8天晴空温度数据所得到的结果是保守的,因为在8天期间内格陵兰表面可能会存在温度高于-1 ℃的更大范围的情况。

2.2.2 趋势分析

采用一元线性回归法拟合格陵兰冰盖表面温度时间空间年际变化。以线性函数T(t)=at+b拟合,按最小二乘法求得aba值的正负表示该变量随时间变化的方向,正值表示随时间呈增长的趋势,负值表示有减小的趋势,a绝对值的大小表示变化的快慢程度。

2.2.3 相关性分析

采用皮尔逊相关系数分析格陵兰冰盖表面温度变化与大气环流指数变化之间的相关性,相关系数越接近1表明两个变量之间正相关越强,越接近 -1表明两个变量之间负相关越强,0表示两个变量之间无关。

3 结果

3.1 格陵兰冰盖夏季表面温度总体变化特征

2000—2020年格陵兰冰盖夏季(6—8月)表面平均温度以及六个主要区域的温度如图2表1所示。结果表明,整个冰盖,近21年夏季平均表面温度为(-7.72±4.34) ℃。其中,2012年夏季表面温度最高,达到(-5.92±4.01) ℃;2000年夏季冰盖表面温度最低,为(-9.24±4.91) ℃。研究期间,2012年和2019年六个区域都表现出较高的温度。2012年,除冰盖北部区域,其他区域表面温度均上升至最高值,特别是东南部和西南部,温度分别达到(-2.74±3.39) ℃和(-1.72±1.87) ℃;2019年冰盖北部区域表面温度达到最高,同时,冰盖表面温度整体上升至较高值(-6.51±3.93) ℃,然而2013年和2020年表面温度都有所下降。2000年,冰盖夏季平均表面温度最低,北部、东北部以及西北部都呈现出近21年的最低温度,而其余流域最低温度分别出现在2018年和2020年。其中,2012年格陵兰温暖的夏季以高压阻塞为特征,冰盖南部反气旋的大气条件导致一股相对温暖的西南气流平流到格陵兰冰盖的西侧,从而导致了冰盖更温暖更干燥,同时,较低的反照率也加剧了冰盖的高温4043

图2

图2   格陵兰冰盖夏季表面平均温度,绿色的为陆地及非连续的冰川

Fig. 2   Mean summer surface temperature of the Greenland Ice Sheet, green are terrestrial and discontinuous glaciers


表1   冰盖整体以及各流域夏季平均表面温度和标准差 (℃)

Table 1  Mean summer surface temperature and standard deviation of the Greenland Ice Sheet as a whole and basins

NONESESWCWNW平均
2000年-10.22±3.65-10.80±4.14-4.51±4.39-3.60±2.94-10.07±4.22-10.57±5.06-9.24±4.91
2001年-8.98±3.72-10.39±4.27-5.35±4.94-4.24±3.24-8.91±4.16-8.24±4.85-8.46±4.65
2002年-7.72±3.28-8.86±3.92-5.97±5.02-4.57±3.51-9.79±4.22-8.49±4.39-7.96±4.30
2003年-7.04±3.37-8.18±4.07-4.04±4.09-2.82±2.68-8.04±4.13-8.63±4.51-7.08±4.36
2004年-8.20±3.23-8.51±3.66-4.38±4.36-3.21±2.70-8.19±4.04-8.79±4.20-7.55±4.18
2005年-7.29±3.32-9.10±3.94-5.08±4.49-4.01±2.92-8.63±4.18-7.72±4.50-7.54±4.28
2006年-9.31±3.14-9.63±3.99-4.91±4.62-3.83±3.01-9.56±4.30-9.06±4.65-8.43±4.48
2007年-7.63±3.60-8.93±3.99-4.05±4.35-2.66±2.33-8.04±4.34-7.41±4.65-7.21±4.49
2008年-6.45±3.24-8.82±4.58-4.89±4.65-3.48±3.01-8.27±4.13-7.42±4.38-7.12±4.48
2009年-7.85±3.18-10.18±3.93-5.41±4.53-4.07±3.22-9.24±4.56-8.27±4.43-8.20±4.43
2010年-7.42±3.57-7.96±3.84-3.33±3.96-2.27±2.52-7.34±4.10-7.44±4.71-6.68±4.36
2011年-6.82±3.23-8.96±3.84-4.56±4.18-3.26±2.71-8.08±3.81-7.55±4.14-7.18±4.15
2012年-6.48±3.42-7.24±3.53-2.74±3.39-1.72±1.87-6.34±3.70-6.74±4.26-5.92±4.01
2013年-9.83±3.44-10.04±4.03-5.52±4.36-4.09±3.03-9.32±4.37-9.49±4.68-8.81±4.51
2014年-8.16±3.55-8.82±4.22-3.89±4.23-2.53±2.69-8.78±4.32-9.48±4.75-7.72±4.68
2015年-7.51±3.16-9.22±3.79-4.96±4.31-4.44±3.17-8.67±3.71-7.76±4.18-7.65±4.09
2016年-7.56±3.38-8.73±4.11-3.56±4.05-2.36±2.57-7.60±4.00-8.07±4.76-7.12±4.50
2017年-9.01±3.60-9.29±4.09-4.87±4.65-3.85±3.31-9.49±4.25-9.31±4.58-8.30±4.56
2018年-9.25±3.61-10.21±3.81-5.42±4.38-4.58±3.06-10.09±4.47-10.00±4.95-8.94±4.54
2019年-6.35±3.02-7.81±3.66-3.84±3.79-3.33±2.66-7.48±3.89-7.12±4.15-6.51±3.93
2020年-7.58±3.62-10.43±4.46-6.08±5.43-4.42±3.44-9.29±4.44-8.69±4.69-8.47±4.80
平均-7.94±3.34-9.13±3.98-4.64±4.31-3.49±2.80-8.66±4.11-8.93±4.49-7.72±4.34

注:保留小数点后两位,下划波浪线的为各个区域温度最低对应的年份,双下划线为最高温度对应的年份。

新窗口打开| 下载CSV


为进一步了解格陵兰冰盖夏季表面温度变化,基于2000—2020年平均表面温度,计算了格陵兰冰盖月平均温度距平以及夏季平均温度距平如图3所示。同时,将由ERA5再分析数据计算得到的格陵兰冰盖夏季平均近地表气温年际距平也显示在图中。通过对温度距平进行线性拟合分析,结果得到,夏季平均、6月和7月表现出不同的增长趋势,并分别为0.16 ℃·(10a)-1、0.03 ℃·(10a)-1和0.79 ℃·(10a)-1,其中仅有7月表面温度的上升趋势率[0.79 ℃·(10a)-1]通过90%置信度检验。

图3

图3   2000—2020年格陵兰冰盖夏季各月以及平均表面温度距平变化

Fig. 3   Monthly and mean surface temperature anomalies of the Greenland Ice Sheet during the summer of 2000—2020


3.2 格陵兰冰盖夏季表面温度时空变化特征

采用一元线性回归方法对格陵兰冰盖夏季表面温度进行了计算(图4表2)。由图3可以看出,格陵兰冰盖表面温度在2012年前后表现出大致相反的变化趋势,于是本文将分别对2000—2012年和2012—2020年格陵兰冰盖表面温度变化趋势进行分析。整个研究期间而言,格陵兰冰盖夏季表面温度呈现出微弱的增长趋势,增长率为0.16 ℃·(10a)-1,且在六个区域都表现出不同的增长趋势,分别为0.34 ℃·(10a)-1、0.14 ℃·(10a)-1、0.09 ℃·(10a)-1、0.01 ℃·(10a)-1、0.12 ℃·(10a)-1、0.13 ℃·(10a)-1。冰盖东南部和西南部是表面温度最高的区域,平均表面温度分别达到(-4.64±4.31) ℃和 (-3.49±2.80) ℃,但两个区域却有最小的变化率,变化率均小于0.1 ℃·(10a)-1。Hall等34对于格陵兰冰盖的研究也表明,东南部和西南部区域较于其他流域具有较高的温度但是却具有较小的变化率,这可能是由于东南部和西南部区域面积较小,但与较大比例的海岸地区相邻所导致的25,海岸地区温度往往高于冰盖内陆温度,同时,东南部的位势高度也较低。2012年前后冰盖表面温度的趋势表现出巨大的差异性,2012年之前冰盖整体表现出1.61 ℃·(10a)-1P<0.05)的增长趋势,随后波动下降。且2000—2012年夏季格陵兰冰盖表面温度在六个区域都表现出较为显著的趋势,增长率介于1.16~1.93 ℃·(10a)-1之间,增长率最大的是北部,最小的是东南部。相反,2012—2020年格陵兰冰盖夏季表面温度除北部区域外都表现出不同程度的减少,减少率介于0.54~2.04 ℃·(10a)-1之间,减少率最大的为冰盖中西部,最小的为西北部。

图4

图4   2000—2020年格陵兰冰盖夏季表面温度变化率(单位:℃·a-1

Fig. 4   Summer surface temperature variability of the Greenland Ice Sheet from 2000 to 2020 (unit: ℃·a-1


表2   格陵兰冰盖夏季表面温度变化率[单位:℃·(10a)-1

Table 2  The summer surface temperature change rate of the Greenland Ice Sheet [unit: ℃·(10a)-1

区域2000—2012年2012—2020年2000—2020年
NO1.930.400.34
NE1.41-1.490.14
SE1.16-1.880.09
SW1.23-2.000.01
CW1.82-2.040.12
NW1.91-0.540.13
总计1.61-1.070.16

注:2000—2012年趋势(东北和东南部除外)通过95%置信度检验,其他趋势未通过。

新窗口打开| 下载CSV


整个研究期间北部是增温幅度最大的区域,比其他任何区域的两倍都多。同样,2012年之后仅有北部区域呈现微弱的升温[0.40 ℃·(10a)-1]。大尺度的环流变化可以有效解释这一现象。格陵兰冰盖地表径流的增加与格陵兰阻塞频率的增加有关。格陵兰岛上空更频繁和更强的高压阻塞意味着格陵兰岛上空反气旋条件的增加。20世纪90年代以来,反气旋阻塞事件的数量增加了一倍44,导致来自南方的暖气团频繁沿格陵兰冰盖西侧向北平流45,于是夏季格陵兰冰盖西北部区域比南部区域变暖要快。同时,Preece等46确定了三种主要的阻塞类型:(1)一个高振幅的欧米茄块(a high-amplitude Omega block);(2)一个低振幅的、静止的夏季山脊(a lower-amplitude, stationary summer ridge);(3)一个气旋破波模式(a cyclonic wave breaking pattern),并记录了每种阻塞类型的表面能量和物质平衡的时空进程。相对于所有的阻塞期,欧米茄块(Omega block)在北部冰盖产生更多的融化,低振幅的、静止的夏季山脊(cyclonic wave breaking pattern)在格陵兰岛东北部产生更多的融化。且结果表明夏季格陵兰岛最近的阻塞趋势主要是由欧米茄块(Omega block)的增加驱动的。相比格陵兰冰盖南部区域,两种阻塞类型在北部产生更多的融化,降低了表面反照率,同样可能会导致冰盖北部区域升温高于南部。以上原因可能导致了北部区域较大的增温幅度。总体来说,2000—2012年和2012—2020年,格陵兰冰盖夏季表面温度整体和各个部分都表现出巨大的差异性。值得注意的是,由于研究时段较短,每一年的温度都会较容易影响温度变化的整体趋势,特别是温度异常的年份,并非所有趋势在统计上都是显著的,但是基于气象站点气温数据以及气候模型都能得到相似的结果1240

3.3 格陵兰冰盖夏季表面温度、融化范围与物质平衡

基于MODIS温度数据构建了格陵兰冰盖表面每年的融化范围情况。2000—2020年格陵兰冰盖夏季最大融化范围以及各个日期内的融化范围如图5~6所示。其中,最大融化范围由当年夏季8天期间内的最大融化面积所占百分比所得。对于整个研究期间,2000年初期,格陵兰冰盖融化范围呈现出显著上升的趋势(1.20%·a-1P<0.05),在2012年达到峰值,随后波动下降(-2.02%·a-1)。夏季最大融化范围较大的年份为2012年、2019年、2020年,较小的年份为2001年、2018年。经研究发现,2012年夏季格陵兰冰盖经历了21世纪以来最大程度的融化。2012年夏季年每一个8天期间都达到较大融化,其中最大融化记录发生在193天,有52.22%的区域发生融化现象。2018年融化范围最小,在各个流域也都相比其余年份较小,2018年193天得到了格陵兰冰盖最大融化范围唯一低于20%的融化,远低于其他年份。

图5

图5   2000—2020年夏季(153~241天)格陵兰冰盖每8天内发生融化的范围,百分比是相对于整个冰盖融化所占面积,虚线显示的是研究期间7月份

Fig. 5   The extent to which the Greenland Ice Sheet melted every eight days in the summer of 2000—2020 (153~241 days), the percentage is relative to the total melting area of the Ice Sheet, the dotted line shows the month of July during the study period


图6

图6   根据MODIS温度数据绘制得到:2000—2020年夏季最大融化年值趋势线(a);2012年夏季最大融化范围(b);2018年夏季最大融化范围(c)

Fig. 6   trend line of summer maximum annual melting from 2000 to 2020(a), maximum melting range in 2012 summer (b), and maximum melting range in 2018 summer (c), are plotted based on MODIS temperature data


图中虚线内显示的是7月格陵兰冰盖融化情况,7月是格陵兰冰盖融化最大的月份,研究期间所有年份所记录的最大融化均发生在7月。同时,2000年以来,7月格陵兰冰盖平均融化范围最大,6月次之,8月最低,融化范围分别为20.63%、12.32%、8.32%。值得注意的是,由于使用的温度数据为8天晴空表面温度平均值,所以得到的融化范围是较为保守的,但在8天晴空表面温度更高的区域,融化可能更久更剧烈。

IMBIE近年来广泛用于评估格陵兰冰盖的物质平衡47-48图7显示了格陵兰冰盖夏季表面温度、平均融化范围以及物质平衡变化的时间序列。格陵兰冰盖表面温度与融化范围具有较好的一致性,而二者都与物质平衡表现出大致相反的趋势,即表面温度越高融化范围越大物质损失越大。研究期间冰盖物质损失速率总体呈现倒“V”形(即先升后降),增长率和减少率分别为24.63 Gt·a-1和28.77 Gt·a-1,且均通过95%置信度检验。其中,2010—2012年是冰盖物质损失速率最大的阶段,平均损失率达到了(335±65) Gt·a-1,其表面温度也较高,高于任何其他连续三年的平均温度。

图7

图7   2000—2018年格陵兰冰盖物质平衡(灰色表示误差值)和温度变化时间序列

Fig. 7   Time series of Greenland Ice Sheet mass balance (grey error) and temperature change from 2000 to 2018


图8为格陵兰冰盖夏季表面平均温度、平均融化范围以及物质平衡相关性分析,结果显示,冰盖表面平均温度和平均融化范围都与物质平衡具有良好的负相关性,相关系数分别为-0.71和-0.75。夏季表面温度与融化范围具有显著正相关,相关系数高达0.91,且均通过置信度检验。同时,由线性拟合关系可知,格陵兰冰盖夏季表面温度每升高1 ℃,会导致其物质损失增加74.29 Gt·a-1

图8

图8   2000—2018年格陵兰冰盖夏季表面温度、融化范围、物质平衡相关性分析

Fig.8   Correlation analysis of surface temperature, melting extent and mass balance of Greenland Ice Sheet from 2000 to 2018


4 讨论

4.1 大气环流对格陵兰冰盖夏季表面温度的影响

多年来,格陵兰冰盖气候受到各种大气环流的影响,其中,格陵兰冰盖表面温度受到北大西洋涛动(NAO)以及格陵兰阻塞指数(GBI)的影响最为明显37-38。NAO相当于是一个气压场,当NAO指数越大,亚热带的高压比通常更高,冰岛的低压相对更低,气压差的增加导致更多更强的风暴沿更偏北的方向穿越大西洋,这会导致格陵兰岛和加拿大北部更加寒冷干燥,反之,NAO指数越小,亚热带高压相对低,冰岛低压相对高,气压梯度的减小导致穿越由西到东的风暴越来越少,越来越弱,格陵兰岛会更加温暖25。GBI是格陵兰地区(60°~80° N,20°~80° W)的平均500 hPa位势高度,它衡量的是格陵兰上空的阻塞情况,对北半球的气候和天气都有影响。于是,本文利用这两个指数来代表大气环流对于格陵兰冰盖表面温度的影响,分析格陵兰冰盖表面温度与格陵兰地区大气环流变化之间的联系。基于之前的研究,本文计算了近年来夏季格陵兰冰盖NAO和GBI指数,然后将其与格陵兰冰盖表面温度距平进行相关性分析,并对其进行显著性分析(图9~10)。结果表明,格陵兰冰盖表面温度和NAO以及GBI均有显著的相关性,且冰盖表面温度受GBI的影响程度要强于NAO的影响,冰盖表面温度和NAO呈现负相关(r=-0.64),和GBI表现出正相关(r=0.77),且两者都达到95%的置信度。图8显示了格陵兰冰盖表面温度以及两种指数的季节距平时间序列变化,NAO和GBI之间也具有很好的负相关性,NAO在2013年和2018年达到较高值,而GBI这两年也达到最低值。对于温度较高的2012年和2019年,GBI达到最高,NAO也更小。格陵兰冰盖表面温度变化较大的年份,NAO和GBI也均表现出较好的一致性。异常的NAO以及GBI指数导致了近年来出现的两次极端的高温情况。较低的NAO指数使得气压梯度的降低导致更少更弱的风暴由西到东,促使格陵兰冰盖表面的反照率降低,反照率的变化会引起冰盖吸收的太阳辐射发生较大的变化,进一步促进了格陵兰冰盖表面温度的升高。

图9

图9   2000—2020年格陵兰冰盖夏季表面温度、NAO以及GBI距平趋势变化

Fig. 9   Trends of Greenland Ice Sheet summer surface temperature, NAO and GBI Anomaly from 2000 to 2020


图10

图10   格陵兰冰盖夏季表面温度与NAO以及GBI距平相关分析

Fig. 10   Correlation analysis of summer Greenland Ice Sheet surface temperature with NAO and GBI Anomaly


5 结论

基于MODIS温度产品数据,分析了2000—2020年格陵兰冰盖夏季整体及六个主要部分的表面温度以及冰盖表面融化范围的年际变化趋势,并结合IMBIE物质平衡数据分析了表面温度与物质平衡之间的关系,最后探讨了大气环流对于冰盖表面温度的影响,结果表明:

(1)2000—2020年,格陵兰冰盖夏季表面温度表现出0.16 ℃·(10a)-1的微弱的增长趋势。其中,2012年是温度最高的年份,且2012年前格陵兰冰盖呈现出显著的升温趋势[1.61 ℃·(10a)-1P<0.05],随后波动下降[-1.07 ℃·(10a)-1],2012年前后格陵兰冰盖整体和六个部分都表现出较大差异性。

(2)7月是格陵兰冰盖表面温度波动最大的月份,也是温度最高的月份。6月和8月温度不明显的变化趋势导致格陵兰冰盖夏季平均表面温度变化相对平缓。研究期间所有六个区域都有不同程度的升温,其中东南部[(-4.64±4.31) ℃]和西南部[(-3.49±2.80) ℃]是温度最高的部分,却具有最小的增长率,正北部分冰盖表面温度增长率[0.34℃·(10a)-1]最大,高于其他任何区域的两倍。

(3)格陵兰冰盖夏季表面温度、平均融化范围以及物质平衡之间都呈现出显著的相关性,同时线性拟合结果表明格陵兰冰盖夏季表面温度每上升1 ℃,会导致其物质损失增加74.29 Gt·a-1

(4)近20年格陵兰冰盖夏季表面温度与NAO和GBI表现出显著的相关性。且冰盖夏季表面温度受GBI的影响要强于NAO的影响,冰盖夏季表面温度和NAO呈现出负相关(r=-0.64,P<0.05),和GBI呈现出正相关(r=0.77,P<0.05)。

参考文献

Allan R PHawkins EBellouin Net al. IPCC, 2021: summary for Policymakers[M]. CambridgeCambridge University Press2021.

[本文引用: 1]

Greve R.

On the response of the Greenland ice sheet to greenhouse climate change

[J]. Climatic Change, 2000463): 289-303.

[本文引用: 1]

Velicogna I.

Increasing rates of ice mass loss from the Greenland and Antarctic Ice Sheets revealed by GRACE

[J]. Geophysical Research Letters, 20093619): L19503.

Zwally H JGiovinetto M BLi Jet al.

Mass changes of the Greenland and Antarctic Ice Sheets and shelves and contributions to sea-level rise: 1992—2002

[J]. Journal of Glaciology, 200551175): 509-527.

[本文引用: 1]

Abdalati WSteffen K.

Update on the Greenland Ice Sheet melt extent: 1979—1999

[J]. Journal of Geopgysical Research Atmospheres, 20001-16.

[本文引用: 1]

Steffen KNghiem S VHuff Ret al.

The melt anomaly of 2002 on the Greenland Ice Sheet from active and passive microwave satellite observations

[J]. Geophysical Research Letters, 20043120): L20402.

Box J EFettweis XStroeve J Cet al.

Greenland Ice Sheet albedo feedback: thermodynamics and atmospheric drivers

[J]. The Cryosphere, 201264): 821-839.

Neumann G.

The extreme melt across the Greenland Ice Sheet in 2012

[J]. Geophysical Research Letters, 20123920): L20502.

[本文引用: 1]

Mote T L.

Greenland surface melt trends 1973-2007: Evidence of a large increase in 2007

[J]. Geophysical Research Letters, 20073422).

[本文引用: 1]

Fettweis XMabille GErpicum Met al.

The 1958-2009 Greenland Ice Sheet surface melt and the mid-tropospheric atmospheric circulation

[J]. Climate Dynamics, 2011361): 139-159.

Hanna EHuybrechts PSteffen Ket al.

Increased runoff from melt from the Greenland Ice Sheet: a response to global warming

[J]. Journal of Climate, 2008212): 331-341.

Hall D KComiso J CDiGirolamo N Eet al.

Variability in the surface temperature and melt extent of the Greenland Ice Sheet from MODIS

[J]. Geophysical Research Letters, 20134010): 2114-2120.

[本文引用: 5]

Zwally H JLi JBrenner A Cet al.

Greenland Ice Sheet mass balance: distribution of increased mass loss with climate warming; 2003–07 versus 1992—2002

[J]. Journal of Glaciology, 201157201): 88-102.

[本文引用: 1]

Chen GuodongZhang Shengjun.

Elevation and Volume Change Determination of Greenland Ice Sheet Based on ICESat Observations

[J]. Chinese Journal of Geophysics, 2019627): 2417-2428.

[本文引用: 1]

陈国栋张胜军.

利用ICESat数据确定格陵兰冰盖高程和体积变化

[J]. 地球物理学报, 2019627): 2417-2428.

[本文引用: 1]

Lu FeiYou WeiFan Dongming.

Analysis of Greenland Ice mass change based on GRACE

[J]. Journal of Geodesy and Geodynamics. 2015354): 640-644.

[本文引用: 1]

卢飞游为范东明.

基于GRACE的格陵兰冰盖质量变化分析

[J]. 大地测量与地球动力学, 2015354): 640-644.

[本文引用: 1]

Bamber J LAspinall W P.

An expert judgement assessment of future sea level rise from the Ice Sheets

[J]. Nature Climate Change, 201334): 424-427.

[本文引用: 1]

Rowley R JKostelnick J CBraaten Det al.

Risk of rising sea level to population and land area

[J]. Eos, Transactions American Geophysical Union, 2007889): 105-107.

[本文引用: 1]

Oerlemans J.

The mass balance of the Greenland Ice Sheet: sensitivity to climate change as revealed by energy-balance modelling

[J]. The Holocene, 199111): 40-48.

[本文引用: 1]

Hanna EHuybrechts PJanssens Iet al.

Runoff and mass balance of the Greenland Ice Sheet: 1958—2003

[J]. Journal of Geophysical Research: Atmospheres, 2005, 110: (D13108).

Janssens IHuybrechts P.

The treatment of meltwater retardation in mass balance parameterisations of the Greenland Ice Sheet

[C] //. International Symposium on the Verification of Cryospheric Models. Zürich (CH)August 1999.

[本文引用: 1]

Yu XWang TDing Met al.

Assessment of MODIS surface temperature products of greenland Ice Sheet using in-situ measurements

[J]. Land, 2022115): 593.

[本文引用: 1]

Leeson A AVan Wessem J MLigtenberg S R Met al.

Regional climate of the Larsen B embayment 1980—2014

[J]. Journal of Glaciology, 201763240): 683-690.

[本文引用: 1]

Zou DefuZhao LinWu Tonghuaet al.

Assessing the applicability of MODIS land surface temperature products in continuous permafrost regions in the central Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2015372): 308-317.

[本文引用: 1]

邹德富赵林吴通华.

MODIS 地表温度产品在青藏高原连续多年冻土区的适用性分析

[J]. 冰川冻土, 2015372): 308-317.

[本文引用: 1]

Wang ZhixiaZhuotong NanZhao Lin.

The applicability of MODIS land surface temperature products to simulating the permafrost distribution over the Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2011331): 132-143.

王之夏南卓铜赵林.

MODIS地表温度产品在青藏高原冻土模拟中的适用性评价

[J]. 冰川冻土, 2011331): 132-143.

Qin DaheDing Yongjin.

Crysopheric changes and their impacts, present, trends and key issues

[J]. Climate Change Research, 200954):187-195.

[本文引用: 3]

秦大河丁永建.

冰冻圈变化及其影响研究——现状、趋势及关键问题

[J]. 气候变化研究进展, 200954):187-195.

[本文引用: 3]

Mortimer C ASharp MWouters B.

Glacier surface temperatures in the Canadian High Arctic, 2000—15

[J]. Journal of Glaciology, 201662235): 963-975.

[本文引用: 1]

Hall D KWilliams Jr R SCasey K Aet al.

Satellite‐derived, melt‐season surface temperature of the Greenland Ice Sheet 2000—2005) and its relationship to mass balance

[J]. Geophysical Research Letters, 20063311L11501.

[本文引用: 1]

Hall D KWilliams R SLuthcke S Bet al.

Greenland Ice Sheet surface temperature, melt and mass loss: 2000–06

[J]. Journal of Glaciology, 200854184): 81-93.

[本文引用: 1]

Tedesco MSerreze MFettweis X.

Diagnosing the extreme surface melt event over southwestern Greenland in 2007

[J]. The Cryosphere, 200822): 159-166.

[本文引用: 1]

Pörtner H ORoberts D CMasson-Delmotte Vet al.

The ocean and cryosphere in a changing climate

[J]. IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, 2019.

Wouters BBamber J LVan den Broeke M Ret al.

Limits in detecting acceleration of Ice Sheet mass loss due to climate variability

[J]. Nature Geoscience, 201368): 613-616.

[本文引用: 1]

Rignot EVelicogna Ivan den Broeke M Ret al.

Acceleration of the contribution of the Greenland and Antarctic Ice Sheets to sea level rise

[J]. Geophysical research letters, 2011385): L05503.

[本文引用: 1]

Peng ZhenranHu ZhengwangWang Linsonget al.

Spatio-temporal characteristics of Ice Sheet melting in Greenland and contributions to sea level rise from 2003 to 2015

[J]. Earth Science, 2021462): 743-758.

[本文引用: 1]

彭桢燃胡正旺王林松.

格陵兰岛冰盖消融时空特征(2003~2015年)及其对海平面上升的贡献

[J]. 地球科学, 2021462): 743-758.

[本文引用: 1]

Hall D KBox J ECasey K Aet al.

Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland

[J]. Remote Sensing of Environment, 200811210): 3739-3749.

[本文引用: 2]

Koenig L SHall D K.

Comparison of satellite, thermochron and air temperatures at Summit, Greenland, during the winter of 2008/09

[J]. Journal of Glaciology, 201056198): 735-741.

[本文引用: 2]

Stroeve J CBox J EHaran T.

Evaluation of the MODIS (MOD10A1) daily snow albedo product over the Greenland Ice Sheet

[J]. Remote Sensing of Environment, 20061052): 155-171.

[本文引用: 1]

Previdi MVeron D E.

North Atlantic Oscillation–related climate variability in a regional atmospheric model

[J]. Journal of Geophysical Research: Atmospheres, 2005110D16106.

[本文引用: 2]

Hanna EJones J MCappelen Jet al.

The influence of North Atlantic atmospheric and oceanic forcing effects on 1900—2010 Greenland summer climate and ice melt/runoff

[J]. International Journal of Climatology, 2013334): 862-880.

[本文引用: 2]

Peng ZWang LChen C.

The short-term mass change of Greenland Ice Sheet and the atmospheric forcing

[C]//IOP Conference Series: Earth and Environmental Science. IOP Publishing20216601): 012091.

[本文引用: 1]

Hanna ECappelen JFettweis Xet al.

Greenland surface air temperature changes from 1981 to 2019 and implications for ice‐sheet melt and mass‐balance change

[J]. International Journal of Climatology, 202141E1336-E1352.

[本文引用: 3]

Hanna ECropper T EHall R Jet al.

Greenland blocking index 1851—2015: a regional climate change signal

[J]. International Journal of Climatology, 20163615): 4847-4861.

[本文引用: 1]

Välisuo IVihma TPirazzini Ret al.

Interannual variability of atmospheric conditions and surface melt in Greenland in 2000—2014

[J]. Journal of Geophysical Research: Atmospheres, 201812318): 10443-10463.

[本文引用: 1]

Hanna EFettweis XMernild S Het al.

Atmospheric and oceanic climate forcing of the exceptional Greenland Ice Sheet surface melt in summer 2012

[J]. International Journal of Climatology, 2014344): 1022-1037.

[本文引用: 1]

Fettweis XHanna ELang Cet al.

Brief communication “Important role of the mid-tropospheric atmospheric circulation in the recent surface melt increase over the Greenland Ice Sheet

[J]. The Cryosphere, 201371): 241-248.

[本文引用: 1]

Tedesco MFettweis X.

Unprecedented atmospheric conditions (1948—2019) drive the 2019 exceptional melting season over the Greenland Ice Sheet

[J]. The Cryosphere, 2020144): 1209-1223.

[本文引用: 1]

Preece J RWachowicz L JMote T Let al.

Summer Greenland Blocking Diversity and Its Impact on the Surface Mass Balance of the Greenland Ice Sheet

[J]. Journal of Geophysical Research: Atmospheres, 20221274): e2021JD035489.

[本文引用: 1]

Mankoff K DFettweis XLangen P Let al.

Greenland Ice Sheet mass balance from 1840 through next week

[J]. Earth System Science Data, 20211310): 5001-5025.

[本文引用: 1]

Simonsen S BBarletta V RColgan W Tet al.

Greenland Ice Sheet mass balance (1992—2020) from calibrated radar altimetry

[J]. Geophysical Research Letters, 2021483): e2020GL091216.

[本文引用: 1]

/