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

方振祥，王宁练，3，李想，张玉杰，邰雪楠

（1.陕西省地表系统与环境承载力重点实验室，陕西 西安 710127；2.西北大学城市与环境学院/地表系统与灾害研究院，陕西 西安 710127；3.中国科学院青藏高原研究所，北京 100101）

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（中分辨率成像光谱仪）温度产品已经广泛应用于极地地区的地表温度时空变化研究［12，26-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 格陵兰岛位置示意图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.shtml）［39］，GBI每月数据来自于美国国家海洋和大气管理局地球系统研究实验室（US NOAA’s ESRL）（https：//www.esrl.noaa.gov/psd/gcos_wgsp/Timeseries/GBI_UL/）［40］以及Hanna［41］。

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地表温度产品可能会产生一些冷偏差［27，35］，Hall等［12］将冰盖表面温度大于等于-1℃定义为融化发生，同样，也将-1℃作为评估格陵兰冰盖表面是否发生融化的阈值，并取每年夏季中格陵兰表面温度高于-1℃的最大范围作为当年冰盖表面的最大融化范围。值得注意的是，基于8天晴空温度数据所得到的结果是保守的，因为在8天期间内格陵兰表面可能会存在温度高于-1℃的更大范围的情况。

2.2.2 趋势分析

采用一元线性回归法拟合格陵兰冰盖表面温度时间空间年际变化。以线性函数T(t)=at+b拟合，按最小二乘法求得a和b。a值的正负表示该变量随时间变化的方向，正值表示随时间呈增长的趋势，负值表示有减小的趋势，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年格陵兰温暖的夏季以高压阻塞为特征，冰盖南部反气旋的大气条件导致一股相对温暖的西南气流平流到格陵兰冰盖的西侧，从而导致了冰盖更温暖更干燥，同时，较低的反照率也加剧了冰盖的高温［40，43］。

图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（unit：℃）

为进一步了解格陵兰冰盖夏季表面温度变化，基于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.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）-1（P＜0.05）的增长趋势，随后波动下降。且2000—2012年夏季格陵兰冰盖表面温度在六个区域都表现出较为显著的趋势，增长率介于1.16～1.93℃·（10a）-1之间，增长率最大的是北部，最小的是东南部。相反，2012—2020年格陵兰冰盖夏季表面温度除北部区域外都表现出不同程度的减少，减少率介于0.54～2.04℃·（10a）-1之间，减少率最大的为冰盖中西部，最小的为西北部。

表2 格陵兰冰盖夏季表面温度变化率［单位：℃·（10a）−1］Table 2 The summer surface temperature change rate of the Greenland Ice Sheet［unit：℃·（10a）−1］

图3 2000—2020年格陵兰冰盖夏季各月以及平均表面温度距平变化Fig.3 Monthly and mean surface temperature anomalies of the Greenland Ice Sheet during the summer of 2000—2020

图4 2000—2020年格陵兰冰盖夏季表面温度变化率（单位：℃·a-1）Fig.4 Summer surface temperature variability of the Greenland Ice Sheet from 2000 to 2020（unit：℃·a-1）

整个研究期间北部是增温幅度最大的区域，比其他任何区域的两倍都多。同样，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年，格陵兰冰盖夏季表面温度整体和各个部分都表现出巨大的差异性。值得注意的是，由于研究时段较短，每一年的温度都会较容易影响温度变化的整体趋势，特别是温度异常的年份，并非所有趋势在统计上都是显著的，但是基于气象站点气温数据以及气候模型都能得到相似的结果［12，40］。

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

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

图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 根据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 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 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 2000—2020年格陵兰冰盖夏季表面温度、NAO以及GBI距平趋势变化Fig.9 Trends of Greenland Ice Sheet summer surface temperature，NAO and GBI Anomaly from 2000 to 2020

5 结论

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

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

图10 格陵兰冰盖夏季表面温度与NAO以及GBI距平相关分析Fig.10 Correlation analysis of summer Greenland Ice Sheet surface temperature with NAO and GBI Anomaly

（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）。