1990—2020年朋曲流域冰川变化及其对气候变化的响应

汤远航， 李梦琦， 邓 铃， 王小丽

（GIS应用研究重庆市高校重点实验室，重庆师范大学地理与旅游学院，重庆 401331）

冰川是冰冻圈的重要组成部分，也是自然界最宝贵的淡水资源，具有重要的地理意义［1］。在气候变暖的趋势下，全球范围内的冰川大都表现出不同程度的退缩趋势［2］。一方面，由于冰雪对太阳辐射具有强烈的反射效应，冰川大规模消融会引起地表辐射以及热量的失衡，从而导致大气环流改变。此外其消融过程中产生的融水注入海洋将会对原有洋流格局产生影响、导致海平面上升，进而改变海洋与大气的相互作用状况，最终影响全球气候格局［3］。另一方面，高山冰川剧烈消融还会导致一系列自然灾害的发生，如冰川泥石流、洪水等，这类现象主要发生于青藏高原东南部和喜马拉雅山等地区，其中我国西藏地区的雅鲁藏布江、波曲及朋曲流域尤为突出［4］。

近期研究结果显示，青藏高原东南部冰川亏损及面积萎缩幅度为最大；其次为喜马拉雅山南缘；青藏高原腹地和末端冰川面积退缩幅度相对较小，冰川物质亏损相对较弱；而帕米尔-喀喇昆仑-西昆仑地区的冰川退缩程度最小，部分冰川甚至呈现前进状态，表现出微弱的物质盈余［5］。朋曲流域位于青藏高原南缘，目前关于朋曲流域的研究，大多集中在区域冰湖面积变化及潜在溃决冰湖分析，尽管有部分关于朋曲流域冰川变化的研究，但研究时段序列较短，且关于朋曲流域冰川与气候变化研究较少［6］，研究时段较早，冰川近况未能同步更新。基于此，本研究采用Landsat 系列遥感影像、数字高程模型（DEM）等数据，使用遥感图像处理及目视解译方法，提取了朋曲流域1990—2020 年冰川边界，研究近30 a来朋曲流域冰川面积的分布、变化及其与气候变化的响应关系，以期为朋曲流域冰川变化的预测及灾害防治提供理论基础，也可为朋曲流域合理利用冰川融水作水资源提供理论依据［7］。

1 研究区概况

喜马拉雅山脉发育着大量冰川，朋曲流域便是其中之一。朋曲流域（27°49′～29°05′ N、85°38′～88°57′ E，平均海拔4500 m 以上）位于青藏高原南部，地处西藏自治区西南面，与尼泊尔接壤，流域南起喜马拉雅山脉，北起冈底斯山脉，西临佩沽措［8-10］（图1）。流域大部分位于喜马拉雅山脉北侧，有明显的焚风效应，气温较低、降雨较少、蒸发量大；少部分位于喜马拉雅山脉南侧，受印度洋暖湿气流的影响，降雨充沛［11］。因区域地势较高、地形复杂，又位于喜马拉雅山脉的高寒地区，有利的地势条件和气候条件下朋曲流域发育了大量冰川，冰川总面积可达1400 km2左右，冰川较为密集且集中分布于沿国境线的喜马拉雅山区，成为了对于周围地区经济和社会发展而言非常重要的固体水库［12］。近年来在全球变暖、强太阳辐射等因素的影响下，朋曲流域冰川处于退缩并加速融化的状态，导致周围地区的冰川湖数量及冰川湖面积增大，在流域水系本就发达的前提下，可能会带来许多严重的自然灾害，如冰湖决堤引发的洪水［9］。

2 数据与方法

2.1 数据来源

2.1.1遥感影像数据和DEM 资料 本文数据源于美国地质勘探局（United States Geological Survey，USGS）、Landsat 系列卫星数据源于地理空间数据云（http://www.gscloud.cn/）［13］，该数据不仅经过了GLS2005 全球陆地控制点系统的正射校正，还通过DEM数据做出了高精度的地形校正，且本文数据的大地测量校正建立在精确的DEM 数据和地面控制点的基础上，满足研究冰川的遥感监测要求［14］。DEM 数据采用的是地理空间数据云提供的30 m 的ASTER GDEM 数据。由于积雪、云层等因素会遮挡遥感视线或造成冰川边界模糊不清，所以本研究选择的影像数据尽量都是夏季无云、少云的清晰图像。但在实际操作中，我们保留了10月左右的部分含云量较高但对冰川区影响较小、总体效果较好的影像数据［13］。本文中最终选用的遥感影像如表1所示。

表1 遥感影像列表Tab.1 Information of the remote sensing images

2.1.2气象格点数据及其他数据 由于当前朋曲流域地区附近的气象站点数量较少，本研究在数据分析并综合考虑到目前已搜集到的数据适用性后，最终选用了源于美国特拉华大学地理系气候研究中心的1970—2017 年气温和降水格点数据。由于该数据已通过地面站点观测的方法，对朋曲流域气象的描述更为精确，容易得出贴合实际的结论［13］。此外，为了验证格点数据的精准性，本研究还收集了朋曲流域附近3 个气象站点的数据与前者进行比较。除上述提到的数据外，本研究还根据第二次冰川编目数据验证了2010 年左右朋曲流域地区的冰川边界提取结果。

2.2 研究方法

本研究采用了精度较高的比值阈值法来提取冰川边界［15-18］，即计算波段比值后生成比值图像，再设定阈值进行决策树分类，在二值化处理比值图像后，依照阈值设定得到冰川边界初步提取结果［19］。经过实验数据的反复尝试后，最终Landsat TM/ETM+影像的阈值设为1.8，而Landsat OLI 数据的阈值设为1.0。为了确保数据的精确性，本研究对比值阈值法计算的结果逐条进行了目视解译［20］。

2.3 精度评价

为了验证冰川边界的精确度，本研究以15 m为距离对1990—2020 年的研究区边界进行了缓冲区处理（表2），由缓冲区面积除以总面积得到误差率。误差率较小，符合本研究要求的精度［21］。冰川与积雪的区分主要依据冰川形态，冰舌部分较为光滑，形态规整，而积雪较呈散状分布，且对于每一条冰川，均采用Google Earth历史数据进行目视解译。

表2 朋曲流域边界缓冲区面积及误差率Tab.2 Area and error rate of boundary buffer zone in Pumqu Basin

3 结果与分析

3.1 冰川分布面积及变化特征

1990、2000、2010 年和2020 年朋曲流域冰川面积分别为1248.70 km2、1209.19 km2、1152.25 km2和1068.10 km2（表3）。对比研究时段内不同时期的冰川分布面积，近30 a朋曲流域冰川整体呈现出退缩趋势，但不同时段冰川退缩率有所不同。1990—2000 年，冰川面积退缩量为39.52 km2，年均退缩率为0.32 %·a-1；2000—2010 年，冰川面积退缩量为56.94 km2，与前一时段相比，该时段冰川加速退缩（0.47 %·a-1）；2010—2020 年，冰川面积退缩量为84.15 km2，年均退缩率为0.73%·a-1，3个时段中，该时段冰川年均退缩率最大。综上可知，1990—2020年朋曲流域冰川面积加速退缩，尤其是2010—2020年来，朋曲流域加速退缩的趋势更为显著（图2）［22-23］。

图2 1990—2020年朋曲流域部分冰川分布特征Fig.2 Distribution characteristics of some glaciers in Pumqu Basin during 1990—2020

表3 1990—2020年朋曲流域冰川分布面积及变化特征Tab.3 Distribution area and variation characteristics of glaciers in Pumqu Basin during 1990—2020

3.2 不同海拔冰川分布特征

1990、2000、2010 年和2020 年朋曲流域冰川面积在不同海拔高度具有明显的差异性，与海拔的变化情况基本一致，具有明显的规律性（图3）。其变化规律表现为：冰川分布面积随海拔的升高不断增大，到达一定高度后开始逐渐缩小。不同时段的冰川分布面积最大及最小区域均出现在同一海拔范围内。海拔4700～4900 m处为冰川分布面积最小的区域，5900～6100 m 为冰川分布面积最大的区域。2020 年冰川分布面积最小值为8.44 km2，最大值为211.79 km2。

本流域冰川呈现出如上的分布特征主要是受到地形地势和气候条件的影响［24］。在对流层内，气温随海拔的升高不断降低，而水分变化情况却不同，根据前人研究，空气中的水汽含量随海拔的升高先增加后减少［25］。由此使得冰川分布面积先随海拔的升高而增加，到达一定海拔高度后开始减少，最终形成图3所示的分布特征。

图3 1990—2020年朋曲流域不同海拔高度冰川分布面积Fig.3 Distribution areas of glaciers at different altitudes in Pumqu Basin during 1990—2020

3.3 不同坡向冰川分布情况及变化特征

1990—2020 年朋曲流域北坡冰川分布面积最大，西南坡冰川分布面积最小（图4）。1990、2000、2010 年和2020 年朋曲流域北坡冰川面积分别为246.91 km2、241.66 km2、233.02 km2和222.87 km2，分别占研究区对应时段冰川总面积的19.77%、19.99%、20.22%和17.27%，冰川面积减小的趋势愈发明显，但冰川面积在总冰川面积中的比例呈现出先增加后减小的趋势；西南坡的冰川分布面积分别为91.56 km2、87.95 km2、82.35 km2和73.58 km2，分别占研究区对应时段冰川总面积的7.33%、7.27%、7.15%和5.70%，冰川面积和占比均呈现较为明显的减小趋势。总体而言，北坡向冰川分布面积较多，南坡向分布面积较少。

图4 1990—2020年朋曲流域不同坡向冰川分布面积Fig.4 Distribution areas of glaciers in different slope directions in Pumqu Basin during 1990—2020

1990—2020 年朋曲流域各坡向冰川均表现为退缩趋势，且退缩趋势逐渐加大（图5）。其中南坡冰川的退缩率为所有坡向中最大，为22.62%，西南坡、东南坡和东坡次之，分别为19.65%、19.33%和16.84%；北坡的冰川退缩率为所有坡向中最小，为9.74%，西北坡、东北坡和西坡次之，分别为10.41%、12.69%和14.34%。

图5 1990—2020年朋曲流域不同坡向冰川退缩率Fig.5 Glacier retreat rates in different slope directions in Pumqu Basin during 1990—2020

3.4 不同坡度冰川分布情况及变化特征

本研究将研究区冰面坡度以5°为间隔分为12个等级［13］（图6），1990—2020 年研究区冰川在0～5°坡度范围分布较少，5°～10°坡度范围为最大，此后依次递减，而坡度＞55°的冰川分布面积稍高于50°～55°坡度范围，但仍低于45°～50°坡度范围。综上可知，朋曲流域冰川大多分布在平缓地带（5°～25°），而分布在坡度＜5°的极平缓地带和较为陡峭地带的冰川较少。

图6 1990—2020年朋曲流域不同坡度冰川分布面积Fig.6 Distribution areas of glaciers in different slopes in Pumqu Basin during 1990—2020

近30 a 来，朋曲流域不同坡度等级的冰川均呈退缩趋势（图7）。整体而言，冰川退缩率随坡度的增加呈现出先增加后减少的趋势，冰川退缩主要发生在15°～40°坡度范围内，其中30°～35°坡度范围退缩率最大，为18.41%，＞55°坡度退缩率最小，仅为8.95%。1990—2000 年不同坡度的退缩率大致相同，中、高坡度（35°～55°）退缩率略高；2000—2010年冰川退缩率在0～35°坡度范围内逐渐增加，而后随坡度的增加而减少，在坡度＞55°后几乎为0；2010—2020 年冰川退缩率大体呈现2 个波峰，一个是在30°～35°坡度范围，退缩率为9.99%，一个是在50°～55°坡度范围，退缩率为9.10%，但在坡度＞55°后依旧减少。综上可知，1990—2020年朋曲流域冰川在15°～45°坡度范围内退缩较大，而在极平缓或极陡峭区域退缩较小。

图7 1990—2020年朋曲流域不同坡度冰川退缩率Fig.7 Glacial retreat rates in different slopes in Pumqu Basin during 1990—2020

3.5 不同规模等级冰川分布情况及变化特征

为了分析1990—2020 年朋曲流域不同规模等级冰川分布情况及变化特征，本研究以喜马拉雅山地区冰川面积大小为标准将研究区朋曲流域冰川分为9 个等级（≤0.2 km2、0.2～0.5 km2、0.5～1 km2、1～2 km2、2～5 km2、5～10 km2、10～20 km2、20～50 km2和50～100 km2），分别计算了1990 年和2020 年研究区不同规模等级冰川条数和冰川面积的分布情况（表4）。

表4 1990年和2020年朋曲流域不同规模冰川条数与面积统计Tab.4 Statistics on the number and area of glaciers of different sizes in Pumqu Basin in 1990 and 2020

3.5.1冰川条数的分布情况和变化特征 朋曲流域冰川以小规模冰川为主，1990 年和2020 年面积≤0.5 km2的冰川分别有489 条和596 条，占各时段冰川总数的59.13%和67.88%。其中，面积稍大的冰川（≥20 km2）仅有12 条和9 条，仅占各时段冰川总数的1.45%和1.03%。研究区内的冰川面积均未超过100 km2。研究发现，总体上冰川条数与冰川规模呈相反趋势，即冰川的规模越大，对应条数越少。

1990—2020年冰川条数整体上呈现增加趋势，从1990 年的827 条增加到2020 年的878 条，增加了51 条。但仅有面积≤0.2 km2和面积处于10～20 km2这2 个规模的冰川条数呈现出增加的趋势，其余规模等级的冰川条数均呈现减少的趋势，其中冰川规模为1～2 km2的冰川条数减少的最多，达24条，其次是规模为0.2～0.5 km2和0.5～1 km2的冰川，均为16条。冰川条数减少的最少的是冰川规模为5～10 km2和20～50 km2的冰川，仅为2 条和3 条。从冰川条数的波动中可以发现，冰川的规模越大，冰川条数的变化越小。

3.5.2冰川面积的分布情况和变化特征 朋曲流域冰川分布面积最大是规模等级为20～50 km2的冰川，1990 年 和2020 年 面 积 分 别 为326.99 km2和244.88 km2，占各时段冰川面积总数的26.19%和22.93%。其次是规模等级为2～5 km2的冰川，1990年和2020 年面积分别为251.95 km2和192.87 km2，占各时段冰川面积总数的20.18%和18.06%。冰川分布面积最小的是规模等级为≤0.2 km2的冰川，1990 年和2020 年冰川分布面积分别为29.86 km2和30.94 km2，占相应时段冰川面积总数的2.39%和2.90%。

分析近30 a 研究区不同规模冰川面积变化可知，研究区内9 个规模等级的冰川，规模等级为≤0.2 km2和10～20 km2的冰川分布面积呈现出增长趋势。其余规模等级的冰川面积均呈减少趋势。在呈现增长趋势的不同规模等级的冰川中，处于10～20 km2的冰川的面积增长最多，由1990 年的87.51 km2增加到2020 年的141.04 km2，冰川分布面积增加了53.53 km2，约占1990年该等级规模冰川总面积的61.16%。而规模≤0.2 km2的冰川面积增长最少，由1990 年的29.86 km2增加到2020 年的30.94 km2，冰川分布面积增加了1.08 km2，约占1990 年该等级规模冰川总面积的3.62%。在呈现减少趋势的不同规模等级的冰川中，处于20～50 km2的冰川的面积减少最多，由1990 年的326.99 km2减少到2020年的244.88 km2，减少了82.11 km2，约占1990年该等级规模冰川总面积的25.11%。处于0.2～0.5 km2的冰川的面积减少最少，由1990年的66.83 km2减少到2020 年的58.86 km2，减少了7.97 km2，约占1990 年该等级规模冰川总面积的11.94%。可以发现，与等级规模较大的冰川相比，规模等级较小的冰川，其冰川面积波动较小。

从冰川条数和冰川面积2个方面分析研究区内1990—2020年的冰川可以发现，冰川的规模等级与冰川的波动呈现相反趋势。即大冰川的波动小，小冰川的波动大。

4 讨论

冰川坡度及面积是影响冰川对气候变化响应的重要地形因子［2］。通过对不同坡向和坡度冰川面积的分析，可以发现冰川对气候变化响应具有一定滞后期，冰川平均规模越小，滞后期越短［26］，前人选择典型区域，利用中值面积指标研究了整个喜马拉雅山冰川与气候变化的滞后期为12 a［27］，喜马拉雅山冰川平均规模为7.78 km2，而朋曲流域冰川平均规模为1.51 km2。因此，以10 a为滞后期，采用的朋曲流域气候数据提前至1970 年。为了更充分地反映研究区气候变化特征，采用年均气温和年降水数据来代表研究年份内不同时段气候的变化情况（图8）。朋曲流域附近3 个站点气候的变化情况（图9）如下：

图8 1970—2017年朋曲流域气温和降水量Fig.8 Temperature and precipitation in Pumqu Basin during 1970—2017

图9 1990—2018年朋曲流域附近站点气温和降水量Fig.9 Temperature and precipitation of stations near Pumqu Basin during 1990—2018

（1）1980—1989 年，朋曲流域年均气温整体出现上下波动趋势，在1984、1985、1988 年时达到最高，为1.5 ℃，最低年均气温为0.7 ℃。而年降水量整体变化波动较大，1981—1984 年变化幅度较小，1985 年达到峰值1773.4 mm，而1986 年急降至谷值1330.8 mm，后大幅上升。在对应时段（1990—2000年），朋曲流域冰川呈现退缩趋势年均退缩率为0.32%·a-1。

（2）1990—1999 年，朋曲流域的年均气温高于上一时段，总体趋势是在波动下降后上升，最高值为2.1 ℃，最低值为0.6 ℃。该阶段的年降水量小于前一阶段，趋势为首先下降到1479.58 mm 的谷值，然后上升到1706.91 mm 的峰值。在对应时段（2000—2010年），朋曲流域冰川主要呈现出退缩趋势，朋曲流域冰川年退缩率为0.47%·a-1，高于第一时段。根据朋曲流域冰川2000—2010 年的气温变化以及降水变化特征，可以看出朋曲流域气温升高且降水减少，这可能是导致冰川加速退缩的主要气候因素［13］。

（3）2000—2009 年，朋曲流域的年均气温逐年上升至2.5 ℃的峰值，该时段的年均气温显著高于前2 个时段。年降水量的总体变化幅度相对较小，呈小幅下降趋势，最低值为1469.51 mm，但年降水量均高于前一时段。对应年份（2010—2020年），朋曲流域的冰川呈退缩趋势，年退缩率为0.73%·a-1，高于前2 个时段。原因分析表明，尽管年降水量高于前2个阶段，但气温却明显升高，这可能是该时段冰川加速退缩的主要气候原因。据此可以推测：与降水减少相比，气温上升可能对冰川退缩有更大的影响。

（4）2010—2017年，朋曲流域的年均气温、年降水量都高于前3 个时段。由此可见未来十几年内，朋曲流域冰川可能将持续处于退缩状态且年退缩率可能逐年增大。

总体而言，1990—2020年朋曲流域冰川加速退缩，尤其是2010—2020 年来，朋曲流域加速退缩趋势更为显著。1990—2020 年朋曲流域各坡向冰川均表现为退缩趋势，且退缩趋势逐渐加大。其中，南坡冰川的退缩率在所有坡向中最大，西南坡、东南坡和东坡次之，北坡最小。并且小冰川较大冰川受到气候等各要素的影响更大，退缩现象更为显著。

对于朋曲流域冰川变化的状况，已有学者进行相关研究，如魏红等［28］采用1970年左右地形图数据与2001 年的ASTER 遥感影像数据发现，研究时段内朋曲流域冰川呈退缩趋势，退缩率为0.29%·a-1；Che等［29］分析了1970—2013年朋曲流域冰川分布和变化状况，结果表明该区域冰川呈加速退缩状态，后一时段（2001—2013 年）冰川退缩速率远大于前一时段（1970—2001年），本文的研究结果与前人一致。与本区域同类型的研究而言，本研究的时间序列更长。并且对于小范围的冰川研究来说，人工目视解译、遥感和地理信息系统技术相结合的研究方法是比较科学且全面的。因此，本研究对区域冰川变化和冰湖预测具有一定的补充作用，对冰川变化预测和灾害防治有重要意义。

5 结论

本文基于Landsat系列遥感影像数据，通过比值阈值法与目视解译，提取了1990—2020年朋曲流域冰川分布及变化特征，同时结合气象格点数据，分析了研究区冰川变化及其对气候变化的响应，得出以下结论：

（1）1990—2020 年朋曲流域冰川加速退缩，且加速退缩趋势在2010—2020年中尤为明显。

（2）朋曲流域冰川分布面积随着海拔的升高而不断扩大，达到一定高度后逐渐缩小，原因是受地形、地势和气候条件的影响。

（3）1990—2020年朋曲流域冰川所有坡向冰川均呈退缩趋势。南坡冰川退缩率最大，北坡冰川退缩率最小。研究区北坡的冰川分布面积最大，冰川面积呈下降趋势，但在冰川总面积中的比例呈上升趋势；西南坡冰川分布面积最小，冰川面积和比例均呈下降趋势。总体上看，北坡的冰川分布面积较大，南坡的冰川分布面积较小。

（4）朋曲流域的冰川大多分布在平缓带（5°～25°），而极平缓带（＜5°）或较为陡峭的坡带冰川面积相对较小。近30 a中，朋曲流域不同坡度的冰川均呈现出退缩的趋势。退缩率呈现出随坡度的增加先增大后减小的趋势。在坡度15°～45°范围内退缩较大，而在极平缓或极陡峭区域退缩较小。

（5）以10 a为滞后期研究气候对冰川退缩的影响结果表明：气温升高和降水减少可能是冰川加速退缩的主要原因，而前者对冰川退缩的影响大于后者。可以预见，未来十几年内，朋曲流域冰川可能将继续处于退缩状态，年退缩率可能会继续增加。与大冰川相比，小冰川更易受到气候的影响。

致谢：十分感谢冀琴老师对本文提出大量宝贵的修改意见，同时对USGS、地理空间数据云、美国特拉华大学地理系气候研究中心提供的数据以及寒区旱区科学数据中心共享的中国第二次冰川编目数据一并表示衷心的感谢。