基于GF-1光谱数据的青藏地区冰川资源现状遥感调查
【类型】期刊
【作者】安国英,韩磊,黄树春,谷延群,支瑞荣,郭兆成,童立强(中国国土资源航空物探遥感中心;四川省地质调查院;中国地质大学(武汉);湖南省遥感中心;河北省地矿局水文工程地质勘查院)
【作者单位】中国国土资源航空物探遥感中心;四川省地质调查院;中国地质大学(武汉);湖南省遥感中心;河北省地矿局水文工程地质勘查院
【刊名】现代地质
【关键词】 冰川资源;冰川分布;遥感;冰川编目;青藏地区
【资助项】中国地质调查局“青藏冰川变化与冰湖溃决灾害遥感综合调查”项目(121201203000160012)
【ISSN号】1000-8527
【页码】P584-594
【年份】2019
【期号】第3期
【期刊卷】1;|7;|8;|4
【摘要】以2014—2015年的GF-1为主、少量OLI影像为基础,参考第二次中国冰川目录等文献资料,修编完成青海省和西藏自治区两省区的现代冰川编目,查明青藏两省区目前共有冰川24 796条,总面积约2.624×10~4km~2,约占青藏两省区区域面积的1.37%,冰川储量为2.027×1~03~2.121×10~3km~3。调查区冰川数量以面积<1.0 km~2、冰川面积介于1.0~10.0 km~2之间的冰川为主,其中面积<1.0 km~2的冰川有19 983条,占总数量的80.59%,面积介于1.0~10.0km~2之间的冰川面积为11 962.40 km~2,占总面积的45.59%;面积最大的中锋冰川的面积达237.37 km~2。调查区内的山系(高原)均有冰川分布,念青唐古拉山冰川数量最多,其次是喜马拉雅山和冈底斯山,这3座山系冰川数量占调查区内冰川总数量的63.33%;念青唐古拉山、喜马拉雅山和昆仑山的冰川面积和冰储量位列前3位,其冰川面积和冰储量分别占总数的68.09%和73.44%;然而昆仑山和羌塘高原的单条冰川的平均面积大于念青唐古拉山和喜马拉雅山的平均面积。从冰川海拔分布来看,海拔5 000~6 500 m之间是冰川集中发育区域,约占调查区冰川数量和冰川总面积的85%以上。调查区的冰川在各流域的分布差异显著,恒河流域是冰川分布数量最多、面积最大的一级外流区,其数量占冰川总量的47%以上,面积占总面积的52%以上;青藏高原内陆流域的冰川数量、面积次之,其冰川数量占总数量的21%,面积占总面积的24%以上,并且内流区单条冰川的平均面积略大于外流区的平均面积。总体上,西藏的冰川数量、面积和冰储量分别占西藏和青海两省区的84.92%、84.92%、86.68%,单条冰川的平均面积两省区相近。
【全文】 文献传递
基于GF-1光谱数据的青藏地区冰川资源现状遥感调查
摘要:以2014—2015年的GF-1为主、少量OLI影像为基础,参考第二次中国冰川目录等文献资料,修编完成青海省和西藏自治区两省区的现代冰川编目,查明青藏两省区目前共有冰川24 796条,总面积约2.624×104 km2,约占青藏两省区区域面积的1.37%,冰川储量为2.027×103~2.121×103 km3。调查区冰川数量以面积<1.0 km2、冰川面积介于1.0~10.0 km2之间的冰川为主,其中面积<1.0 km2的冰川有19 983条,占总数量的80.59%,面积介于1.0~10.0 km2之间的冰川面积为11 962.40 km2,占总面积的45.59%;面积最大的中锋冰川的面积达237.37 km2。调查区内的山系(高原)均有冰川分布,念青唐古拉山冰川数量最多,其次是喜马拉雅山和冈底斯山,这3座山系冰川数量占调查区内冰川总数量的63.33%;念青唐古拉山、喜马拉雅山和昆仑山的冰川面积和冰储量位列前3位,其冰川面积和冰储量分别占总数的68.09%和73.44%;然而昆仑山和羌塘高原的单条冰川的平均面积大于念青唐古拉山和喜马拉雅山的平均面积。从冰川海拔分布来看,海拔5 000~6 500 m之间是冰川集中发育区域,约占调查区冰川数量和冰川总面积的85%以上。调查区的冰川在各流域的分布差异显著,恒河流域是冰川分布数量最多、面积最大的一级外流区,其数量占冰川总量的47%以上,面积占总面积的52%以上;青藏高原内陆流域的冰川数量、面积次之,其冰川数量占总数量的21%,面积占总面积的24%以上,并且内流区单条冰川的平均面积略大于外流区的平均面积。总体上,西藏的冰川数量、面积和冰储量分别占西藏和青海两省区的84.92%、84.92%、86.68%,单条冰川的平均面积两省区相近。
关键词:冰川资源;冰川分布;遥感;冰川编目;青藏地区
0 引 言
中国是世界上中低纬度山岳冰川最发育的国家。2012年完成的第二次冰川编目中代表21世纪初期中国冰川的分布状况,显示中国共有冰川48 571条,总面积约5.18×104 km2,冰储量4.3×104~4.7×104 km3[1],与代表20世纪60—80年代第一次冰川编目结果(中国有冰川46 377条,总面积约5.94×104 km2,冰储量5 600 km3[2])相比,已经发生较大变化。第二次冰川编目统计显示,青海和西藏两省区的冰川数量、冰川面积和冰储量约占全国冰川总数量、面积及冰储量的1/2以上,其中西藏的冰川数量和面积均列第一,而冰储量略低于新疆;青海的冰川数量和规模仅次于西藏和新疆而位居第3。冰川作为重要的固体水资源,主要的大江、大河都有冰川融水补给,尤其是西北干旱区的水资源很大程度上依赖于冰川融水,因此冰川与冰冻圈的研究一直以来备受中国政府和学术界关注[2-5]。
从20世纪70年代以来,随着航空摄影、遥感光学影像、雷达、激光测高、数字高程模型等技术的进步和发展,采用遥感和地理信息相结合的技术已成为冰川研究的重要方法。例如,对于现代冰川的考察研究,由过去的单条冰川深入到区域性全面的编目与冰川的变化监测,从冰川雪线变化、冰心、冰川能量平衡等到冰川变化与水文、气象关系等多个方面观测[6-11];对冰川变化的研究,从单体冰川定位监测、区域考察到不同时相的遥感影像或航片资料对比;对冰川资源的评价由定性、半定量到定量,从面积变化到高程变化进而确定厚度变化及冰储量等[12-13]。冰川研究所使用的遥感信息,由航摄的黑白航片到彩色航片,及Landsat MSS、TM/ETM+、ASTER、SPOT、IKONOS、ALOS PALSAR、SRTM、GDEM、InSAR等多种数据类型[6,12-18],目前国内外的多家机构免费共享的光学和微波遥感数据,如TM/ETM+,ASTER、HJ-1A/B、ERS-1/2、ENVISAT、PALSAR等,为冰川变化监测提供了重要的数据基础[16-17]。而美国2013年2月11日发射的Land-Sad 8 OLI卫星数据,弥补自2003年5月31日LandSat 7机载扫描校正仪发生故障而中断的大面积冰川卫星监测共享数据。中国同年4 月 26 日发射的GF-1卫星,除用于精细探测局部的环境污染状况,还可用于水、大气环境和生态环境质量等大范围的宏观监测与评价[19],截至目前还没见使用该数据对冰冻圈进行调查分析的成果发表。
在全球气候变暖,高海拔地区特别是热带高海拔地区变暖趋势更为明显的背景下,冰川普遍退缩,由冰川变化引发的灾害事件增多,冰湖溃决造成突发性洪水和泥石流等灾害也成为高海拔地区危害严重的地质灾害类型之一[4,20-21]。为落实防灾、减灾和生态环境文明建设,2016年,中国地质调查局实施“青藏冰川变化与冰湖溃决灾害遥感综合调查”项目,旨在充分认识青藏地区冰川与冰湖分布及变化状况,为国土空间管控提供基础数据支持,项目任务之一就是完成三期青藏地区冰川与冰湖编目,以认清青藏地区冰川、冰湖动态变化。本文基于上述项目工作成果,选用国产卫星GF-1光谱数据开展遥感调查,以青海和西藏的冰川分布为例,对该区最新冰川分布情况进行分析和总结,提供国产卫星在冰川学领域的应用实例,可以为同行和其他研究者提供参考。
1 调查区概况
调查区位于西藏和青海2省区,地理坐标为78°24′—103°55′E,26°48′—39°12′N,面积约192万km2。其地势总体是西高东低,由北西向南东倾斜,西藏自治区平均海拔在4 000 m以上[22],青海省平均海拔在3 000 m以上[23]。区内主体为高原山地地形,在高原山地中夹有宽谷盆地和山间谷地,山地包括昆仑山山系、祁连山山系、唐古拉山山系、念青唐古拉山山系、冈底斯山山系、喜马拉雅山山系等多条长度在1 000 km以上的山脉,其海拔在5 000 m以上大都终年积雪,冰川广布,是中国主要的冰川分布区。中国三类冰川,海洋型、亚大陆型和极大陆型冰川[2],从调查区南部、东南部向北部、西北部均有发育。
2 使用数据与方法
2.1 数据源简介
GF-1卫星搭载了具有60 km幅宽的两台2 m分辨率全色/8 m分辨率多光谱相机以及800 km幅宽的四台16 m分辨率多光谱相机,是典型的高时间分辨率(宽覆盖)遥感卫星[19],其中16 m分辨率的数据有4个波段,覆盖周期4天,可用于大范围内的冰川监测(图1(a))。
Landsat-8卫星,是美国陆地探测卫星系列的后续卫星。装备有陆地成像仪(Operational Land Imager,简称“OLI”)和热红外传感器(Thermal Infrared Sensor,简称“TIRS”)。OLI有9个波段的感应器,覆盖了从红外到可见光的波长范围。与Landsat-7卫星的ETM+传感器相比,OLI增加的蓝色波段(0.433~0.453 μm)和短波红外波段(band 9,1.360~1.390 μm),还用于海岸带观测和云检测[24]。Landsat 系列卫星数据,曾是全球冰冻圈冰川监测的主要数据[14-20](图1(b))。

图1 不同遥感影像对比(局部)图
Fig.1 Comparison of GF-1 image and OLI image of glaciers (in part)
(a)GF-1影像图(由波段4、3、2合成);(b)OLI影像图 (由波段7、4、1合成)
表1 青藏地区调查工作使用的影像时间
Table 1 Image schedules used in the survey of the Qinghai-Tibet region

获取时间数量/景获取时间数量/景获取时间数量/景获取时间数量/景获取时间数量/景2014-01142014-0442014-0752014-11122015-0112014-0222014-05162014-0832014-12232015-0312014-0322014-06222014-09102015-111
2.2 使用数据
本次调查工作使用的数据有卫星影像、数字高程模型数据及“中国第二次冰川编目”矢量数据等。其中GF-1和OLI卫星遥感影像,用于提取冰川边界;SRTM V4数据作为DEM模型数据,用于复合冰川面积的分割与冰川几何参数提取;冰川编码依据“中国第二次冰川编目”。而本次调查使用GF-1影像114景、OLI影像3景,包括2014年114景,2015年3景。数据选取原则:(1)以空间分辨率为16 m的GF-1影像为主,以OLI影像为补充,并以在线的Google Earth影像、OLI影像等作为参考;(2)数据时相,首选夏秋时节,由于工作区范围大,实际使用影像各季节均有(表1);(3)总体数据云、雾覆盖率<10%,阴影较少、色调层次分明,经处理后满足工作需要。本次数据,总体上可解译程度较高, GF-1影像质量尚好,可清晰识别冰川,色彩层次丰富,但细节较OLI模糊,藏东南个别地区积雪覆盖率<15%;OLI数据质量较好,纹理清晰,云、雪覆盖少,积雪覆盖普遍<10%(图1)。
2.3 研究方法
首先,使用软件系统对卫星影像数据进行自动化处理,输出用于解译的数字高程模型、正射影像等产品,并将投影设定为Albers正轴等积双标准纬线割圆锥投影,WGS84 坐标系,以便准确计算各条冰川面积参数。其次,提取冰川边界,采用人工目视解译,勾绘冰川分布区边界(图2);对于含有表碛覆盖的冰川,在有经验的冰川专家指导下,消除内部岛状冰碛物,辨别冰川末端和侧碛垄位置。再其次,借助GIS技术利用流域边界水文地质特征和坡向差自动提取山脊线,并经过人工修正获取冰川区山脊线矢量作为分冰岭(图3),划分出单条冰川(图4),完成冰川矢量边界提取。

图2 影像图上勾绘冰川边界
Fig.2 Glacial boundaries delineated on remote sensing image
在提取单条冰川后,参考“中国第二次冰川编目”中的冰川编码,逐条进行属性信息的编录,包括冰川编码、山系、流域编码、行政区域信息及由GIS软件提取的几何参数和计算结果,如面积、冰储量、周长、朝向、坡度、平均高程、最大高程、最小高程、中值面积高程等数据,还填写了数据源、影像采集时间和提取人、审核人等字段信息。其中,冰储量计算方法采用“中国第二次冰川编目”使用的由Radi和Hock[25]、Grinsted[26]提出的计算方法,本文中冰储量为上述两种算法的统计结果。
V=0.0365A1.375
(1)
V=0.0433A1.29
(2)
式中:V为冰川冰储量,km3;A为冰川面积,km2。

图3 分冰岭自动提取与人工修正结果比对
Fig.3 Comparison between the ice ridge via automatic extraction and the artificial revision
(a)自动提取分冰岭;(b)经人工修正分冰岭

图4 分冰岭分割单条冰川过程
Fig.4 Images of a single glacier divided by ice ridges
(a)复合冰川边界;(b)以山脊线作为分冰岭;(c)划分成单条冰川
最后,完成冰川边界提取与分割、属性信息录入后,经过外业查证后,修订疑似的侧碛垄冰川边界,完善冰川遥感调查工作。本次调查工作采用的最小面积基于影像数据分辨率16 m/像元计算,提取最小冰川面积为0.000 2 km2。
3 结果与讨论
3.1 青藏地区冰川分布与规模
据本次调查统计,青藏地区冰川条数共24 796条,面积26 238.25 km2,冰储量(2 074.13±47.07) km3。其中,青海有冰川3 738条,冰川总面积为3 957.92 km2,冰储量为(276.19±0.31) km3;西藏有冰川21 058条,冰川总面积为22 280.33 km2,冰储量为(1 797.94±46.76) km3;青海冰川数量和面积占总数量的15.08%,而冰储量占总数的13.32%,略低于冰川数量、面积占比。就单体冰川而言,这两个地区的平均面积接近,最小面积也接近,最大面积的冰川出现在西藏的昆仑山。与第二次冰川编目相比,两地区冰川总量呈减少趋势;所占比例上,西藏的冰川总体上有所减少,青海略有增加(表2)。
从调查区冰川面积规模等级组成(表3)可以看出,随着冰川面积等级的增大,冰川数量、面积和冰储量呈先增加后减少趋势。青藏地区冰川数量以面积<1.0 km2的冰川为主,共19 983条,约占冰川总数的80.59%,其面积占总面积的19.88%,冰储量不足总储量的8%,其中面积<0.1 km2的冰川数量占总数的26.19%,冰川面积占总面积的1.27%,冰储量仅占总储量的0.26%;面积介于1.0~10.0 km2的冰川面积最大,占总面积的45.59%,冰储量占总储量的34.95%,位列第二;面积介于10~100 km2冰川数量有366条,仅占总数的1.48%,面积和冰储量分别占总面积的30.43%和总冰储量的46.22%;冰川面积≥100 km2的冰川数量只有7条,占总面积的4.10%,其冰储量占到总冰储量的11.28%。与“中国第二次冰川编目”结果相比,各等级面积的冰川数量所占比例近似;从各等级的冰川面积所占比例得知,面积≥100 km2 的冰川显著减少,面积介于1.0~10.0 km2之间的冰川比例有所增加,其他等级冰川占比近似;而冰储量方面,面积≥100 km2冰川的冰储量比例在减少,面积<100 km2冰川的冰储量比例却在增加。
表2 青藏地区现代冰川信息统计
Table 2 Statistics of modern glaciers in the Qinghai-Tibet region

冰川数量数量/条占比/%冰川面积面积/km2占比/%冰储量冰储量/km3占比/%单体冰川面积/km2平均值最大值最小值青海3 73815.083 957.9215.08276.19±0.3113.321.058 887.070.004 8西藏21 05884.9222 280.3384.921 797.94±46.7686.681.058 0237.370.005 2总计24 796100.0026 238.25100.002 074.13±47.07 100.001.058 2237.370.004 8青海*3 80314.813 936.7514.21274.77±0.32 12.161.035 283.940.010 3西藏*21 87285.1923 795.4185.791 984.70±61.2187.841.087 9237.46 0.010 0总计*25 675100.0027 732.16100.002 259.47±61.48100.001.080 1237.460.010 0
注:“*”数据根据“中国第二次冰川编目”数据库信息统计。
表3 青藏地区不同规模冰川信息统计
Table 3 Statistics of glacier areas in the Qinghai-Tibet region

规模/km2冰川数量数量/条占比/%占比*/%冰川面积面积/km2占比/%占比*/%冰储量冰储量/km3占比/%占比*/%<0.16 49526.1926.06333.821.271.245.33±1.07 0.260.220.1~1.013 48854.4055.374 882.7118.6118.27151.27±17.467.296.481.0~104 44017.9117.0811 962.4045.5942.36724.99±20.8334.9529.3710~1003661.481.447 983.9330.4330.44958.59±55.9046.2243.24>10070.030.041 075.384.107.69233.94±30.5211.2820.69总计24 796100.00100.0026 238.25100.00100.002 074.13±47.07 100.00100.00
注:“占比*”的数据引自刘时银等[1],不同冰川面积分级的数据所占百分比合并计算。
表4 青藏地区面积大于86 km2的现代冰川信息对比
Table 4 Statistics of modern glaciers with area over 86 km2 in the Qinghai-Tibet region

名称所处山系长度/km面积/km2冰储量/km3物理类型所属流域末端高程/m冰川形态所在省区中峰冰川昆仑山26.675237.3758.79±8.58极大陆型阿克赛钦湖5 350.7冰帽 西藏雅弄冰川念青唐古拉山31.116170.9137.88±5.01海洋型 雅鲁藏布江3 987.4山谷冰川西藏恰青冰川念青唐古拉山37.454167.8436.97±4.86海洋型 雅鲁藏布江2 682.5山谷冰川西藏崇测冰川昆仑山25.642166.2636.50±4.79极大陆型阿克赛钦湖5 305.3冰帽 西藏夏曲冰川念青唐古拉山22.874113.3421.86±2.52海洋型 雅鲁藏布江3 150.9山谷冰川西藏古里雅冰川昆仑山16.712111.8821.49±2.46极大陆型阿克赛钦湖5 487.3冰帽 西藏弓形冰川昆仑山18.291107.7820.44±2.31极大陆型阿克赛钦湖5 360.8冰帽 西藏喜日弄浦冰川念青唐古拉山19.56996.1017.54±1.90海洋型 雅鲁藏布江4 277.5山谷冰川西藏莫诺马哈冰川昆仑山17.67087.0715.37±1.60极大陆型柴达木内陆4 898.6山谷冰川青海增冰川喜马拉雅山25.32386.4915.23±1.58亚大陆型雅鲁藏布江5 312.5山谷冰川西藏
如表4所示,调查区内面积最大的10条冰川,单条冰川面积≥86.0 km2,总面积为1 345.04 km2,占区内冰川总面积的5.13%;冰储量为(282.08±35.61) km3,占总冰储量的13.60%。这几条巨大的冰川集中分布在昆仑山、念青唐古拉山,其中分布在昆仑山南坡的中峰冰川是调查区中面积最大的冰川,面积为237.37 km2。这10条冰川中有7条冰川面积超过100 km2,与“中国第二次冰川编目”成果中调查区内面积≥100 km2的8条冰川相比[1],已经减少了1条,即那龙冰川面积发生巨大的变化,已经小于100 km2。
3.2 各山系(高原)冰川分布
调查区自北向南依次发育有祁连山、喀喇昆仑山、昆仑山、羌塘高原、唐古拉山、冈底斯山、念青唐古拉山、喜马拉雅山和横断山等9座山系(高原),这些山脉或山峰的绝对海拔高度及冰川平衡线以上的相对高差,为冰川形成提供了广阔的积累空间和水热条件,从而成为冰川集中分布的区域。以500 m为间距,将冰川平均高程划分为不同的高程区间,分别对冰川数量、面积和冰储量等进行统计(表5和图5(a)),查明5 000~5 500 m海拔区间的冰川数量、面积和冰储量分别占总数的32.48%、33.29%和31.67%;而5 500~6 000 m海拔区间的冰川数量、面积和冰储量分别占总数的45.09%、39.06%和36.72%,冰川规模最大;而6 000~6 500 m海拔区间冰川数量、面积和冰储量分别占12.06%、14.64%和19.21%。这三个海拔区间,冰川数量、面积和冰储量分别占总量的89.63%、86.99%和87.60%。可见,青藏地区5 000~6 500 m之间是冰川集中发育区域,占冰川总量的85%以上。海拔4 000 m以下的地区受气温影响不利于冰川发育,而6 500 m以上由于山地面积较小,冰川面积相应也较少,这两个区域的冰川数量之和占总数量的0.25%,冰川面积总和也仅占冰川总面积的0.10%。从单条冰川看,在4 500~5 000 m的海拔区间,平均面积最大,为1.38 km2;>6 500 m海拔区间,平均面积最小,仅为0.27 km2。以中值面积高程为例,统计各山系(高原)中冰川的分布范围,如图5(b)所示。喜马拉雅山和念青唐古拉山海拔分布范围最大,冰川数量众多;其次是冈底斯山、横断山和昆仑山,羌塘高原和横断山近似;唐古拉山和昆仑山分布范围最小。冰川发育和规模受多重因素制约,山地海拔和高海拔地区的宽展程度是冰川存在的必要条件,热量和水分条件组合是冰川发育的物质基础。低纬度地区山脉的冰川,只能在较高的海拔范围存在,有丰富的降水能维持冰川向较低海拔延伸;纬度高、温度低和较丰富的降水极有利于冰川发育。
表5 青藏地区现代冰川的平均高程分布信息统计
Table 5 Statistics of modern glaciers at different altitudes in the Qinghai-Tibet region

平均高程/m数量数量/条占比/%面积面积/km2占比/%冰储量冰储量/km3占比/%单条冰川平均值/km2最大值/km2最小值/km23 500~4 000200.08 16.880.061.12±0.010.050.844 210.420.086 94 000~4 5001990.80184.750.7011.02±0.230.530.928 418.190.030 24 500~5 0002 3129.323 201.3412.20244.59±3.6011.791.384 7113.340.005 45 000~5 5008 05332.488 734.8833.29656.84±11.2531.671.084 7170.910.004 85 500~6 00011 18045.0910 247.7539.06761.61±8.5336.720.916 687.070.005 26 000~6 5002 99012.063 841.3414.64398.53±23.9619.211.284 7237.370.007 4>6 500 420.17 11.300.04 0.43±0.030.020.269 0 3.010.016 1总计24 796100.0026 238.25100.002 074.13±47.07100.001.058 2237.370.004 8

图5 青藏地区不同高程冰川面积分布(a)及山系冰川高程分布范围(b)
Fig.5 Distribution of glacial sizes at different altitudes (a) and elevations of modern glaciers in mountain ranges (b) in the Qinghai-Tibet region
表6 青藏地区现代冰川在各山系的分布信息统计
Table 6 Statistics of modern glaciers on the Qinghai-Tibetan mountains

行标签数量数量/条占比/%面积面积/km2占比/%冰储量储量/km3 占比/%单条冰川平均值/km2最大值/km2最小值/km2冈底斯山3 86515.591 296.254.9455.85±3.48 2.690.335 421.540.007 4横断山1 1984.83829.663.1641.35±1.941.990.692 515.940.011 2喀喇昆仑山1 7467.041 529.005.83105.25±0.30 5.070.875 747.890.005 2昆仑山2 0598.303 214.3512.25353.67±24.08 17.051.561 1237.370.010 2念青唐古拉山6 36425.678 394.5431.99691.37±19.68 33.331.319 1170.910.010 3祁连山1 1364.58831.293.1746.44±1.45 2.240.731 820.020.004 8羌塘高原1 4685.922 056.137.84162.79±2.63 7.851.400 668.750.010 3唐古拉山1 4886.001 829.416.97139.27±1.636.711.229 454.150.012 4喜马拉雅山5 47222.076 257.6323.85478.15±6.2423.051.143 686.490.009 5总计24 796100.0026 238.25100.00 2 074.13±47.07100.001.058 2237.370.004 8
调查区各山系的冰川统计(表6)显示,喜马拉雅山和念青唐古拉山冰川数量均在5 000条以上,这2座山系共分布冰川11 836条,面积总计14 652.17 km2,分别占调查区冰川数量、面积相应总量的47.74%和55.84%,冰储量占总量的56.39%。冈底斯山中的冰川数量为3 865条,占总数量的15.59%,居第三位,但是冰川面积不足总面积的5%,平均面积最小,仅有0.34 km2,是冰川平均面积最小的山系;昆仑山冰川面积居第三位,占总面积的18.92%,且平均面积最大可达1.56 km2。羌塘高原冰川数量和面积不如前几者,但深居青藏高原腹地,该区域分布若干海拔6 000 m以上的较为平坦的山峰,致使冰川平均面积达1.40 km2,位居第二。喜马拉雅山虽然非常高峻,但山脊较狭窄而限制了冰川扩展,冰川平均面积只有1.14 km2,比总体平均规模略高。冰川数量或面积最少的3座山系分别为冈底斯山、横断山和祁连山,没有冰川面积大于25 km2的冰川,且冰川平均规模均在0.75 km2以下。
3.3 各流域冰川数量与分布
调查区覆盖了8个一级流域(表7),分属内流区和外流区。根据统计,外流区和内流区的冰川数量分别为17 225条和7 571条,相应的面积为17 876.92 km2和8 361.33 km2,分别占冰川总面积的68.13%和31.87%,相应的冰储量为(1 358.10±22.28) km3和(716.03±24.79) km3,分别占总储量的65.48%和34.52%。中国境内的恒河流域,冰川数量最多,面积、冰储量亦最大,分别占冰川总量的47.87%、52.75%和52.29%。其次是内流区的青藏高原内陆冰川数量、面积和冰储量,分别占冰川总量的21.89%、24.75%和28.81%;冰川分布数量最少和冰川规模最小的黄河流域,仅有冰川 161条,面积128.24 km2,冰储量(8.67±0.03) km3。从冰川平均面积来看,长江流域平均面积1.27 km2,为最大;其次是青藏高原内流区,平均面积约1.20 km2,与恒河流域平均面积1.17 km2接近;中国境内印度河上游冰川平均面积最小,分别为0.44 km2,其他流域介于1.17~0.44 km2之间。
从一级流域冰川规模等级组成(表8)得知,除印度河流域冰川数量以面积<0.1 km2居多之外,其他流域均以面积为0.1~1.0 km2的冰川数量最多;除黄河流域冰川面积以10.0~100.0 km2最多外,其他流域冰川面积以1.0~10.0 km2为最多;冰川面积≥100 km2的7条冰川,其中青藏高原内流区分布有4条,均分布在阿克赛钦湖二级流域,恒河流域分布有3条,也是集中分布在雅鲁藏布江二级流域内(表4)。黄河、湄公河和印度河3个流域,没有面积≥50.0 km2的冰川分布,其中湄公河流域仅有3条冰川的面积>10 km2,单条冰川面积均小于16.0 km2,共森龙巴冰川最大,面积仅为15.94 km2;黄河流域有3条面积≥10.0 km2的冰川,哈龙冰川最大,面积为20.74 km2。萨尔温江流域有2条面积>50 km2的冰川,即若果冰川,其面积为57.38 km2,以及麻果龙冰川,其面积为51.76 km2;恒河流域的冰川,有5条面积>86 km2,15条面积>50 km2;印度河流域有3条面积>50 km2的冰川。
表7 青藏地区现代冰川在不同流域的分布信息统计
Table 7 Statistics of modern glaciers in different drainage basins of the Qinghai-Tibet region

流域 数量数量/条占比/%面积面积/km2占比/%冰储量冰储量/km3占比/%单条冰川平均值/km2最大值/km2最小值/km2外流区黄河1610.65128.240.49 8.67±0.030.420.796 520.740.011 1长江8913.591 132.814.3282.39±0.463.971.271 454.150.013 0湄公河4681.89269.181.0313.96±0.540.670.575 215.940.010 2萨尔温江2 1338.601 762.456.72122.10±0.635.890.826 357.380.011 2恒河1 187047.8713 839.6552.751 084.54±22.1452.291.165 9170.910.009 5印度河1 7026.86744.582.8446.44±0.382.240.437 548.670.008 1小计17 22569.4717 876.9268.131 358.10±22.2865.481.037 8170.910.008 1内流区东亚内流域 2 1438.641 866.307.11118.40±1.16 5.710.870 987.070.004 8青藏高原内陆5 42821.896 495.0324.75597.64±25.9528.811.196 6237.370.005 2小计7 57130.538 361.3331.87716.03±24.7934.521.104 4237.370.004 8总计24 796100.00 26 238.25100.00 2 074.13±47.07 100.00 1.058 2237.370.004 8
表8 青藏地区一级流域冰川规模分布信息统计
Table 8 Statistics of glaciers in the first-level drainage basins in the Qinghai-Tibet region

规模/km2黄河数量/条面积/km2长江数量/条面积/km2湄公河数量/条面积/km2萨尔温江数量/条面积/km2<0.1482.671438.721227.0344727.760.1~1.08727.23528204.0427884.461 345458.301.0~102346.88201540.5265137.41321773.6610~100351.4719379.54340.2820502.74>10000000000总计161128.248911 132.81468269.182 1331 762.46规模/km2恒河数量/条面积/km2印度河数量/条面积/km2东亚内流域数量/条面积/km2青藏高原内陆数量/条面积/km2<0.12 544137.4088038.7153828.661 77382.870.1~1.06 7252 509.24700216.021 187431.122 638952.301.0~102 4066 495.74111293.794011 078.809122 595.6210~1001924 245.1911196.0617327.711012 240.95>1003452.0900004623.29总计11 87013 839.661 702744.582 1431 866.305 4286 495.03
3.4 各行政区中冰川数量与分布
本次调查显示(表9),西藏7个市(地区)均有冰川分布,阿里、林芝、日喀则地区是冰川主要分布区域,占青藏地区冰川规模的50%以上,其中阿里地区冰川数量最多,为5 964条,冰川面积为5 274.38 km2,冰储量为(479.30±21.80) km3,但冰川面积略小于林芝地区的5 640.47 km2;其次为日喀则,冰川数量为4 135条,冰川面积为3 855.64 km2,冰储量(287.54±3.24) km3,略少于林芝地区的面积和冰储量;拉萨有754条冰川,冰川面积为466.53 km2,冰储量为(24.90±0.87) km3,是西藏境内最小的冰川发育区。就单体冰川而言,林芝地区的平均面积最大,约为1.42 km2;其次为昌都地区,为1.28 km2;拉萨的平均面积最小,仅为0.62 km2。
青海境内的冰川相对较少,最大的冰川是发育在布喀达坂峰(又称新青峰)南坡的莫诺马哈冰川,面积为87.07 km2。青海省境内冰川主要分布在海西、玉树、海北、果洛和海南5个自治州,其中海西州冰川数量2 239条,冰川面积2 628.83 km2,冰储量(183.34±0.13) km3,分别占青海冰川总量的59.88%、66.42%和66.38%;其次是玉树州,冰川数量为1 084条,冰川面积1 129.47 km2,冰储量(81.76±0.71) km3,分别占青海冰川总量的29.0%、28.54%和29.60%;海南州冰川规模最小,仅分布有6条小冰川,总面积也仅为0.71 km2。就单条冰川平均面积而言,海西州冰川面积最大,为1.17km2;果洛州次之,为1.15 km2;海南州的平均面积最小,仅为0.12 km2。中国的长江和黄河均发源于青海省的冰川区,冰川融水对于补给江河上游径流具有重要意义,而环绕柴达木盆地山脉上的冰川对于该区农业和石油化工的发展具有重要的经济价值[1]。
表9 青藏地区不同地区现代冰川信息统计
Table 9 Statistics of modern glaciers in the districts of the Qinghai-Tibet region

行政区划冰川数量数量/条 占比/%冰川面积面积/km2占比/%冰储量 冰储量/km3 占比/%单条冰川平均值/km2最大值/km2最小值/km2青海西藏果洛州940.38108.030.417.96±0.04 0.381.149 220.7380.011 1海北州3151.2790.880.353.12±0.300.150.288 52.9780.009 2海南州60.020.7100.0200.118 50.2110.049 2海西州2 2399.032 628.8310.02183.34±0.138.841.174 154.1500.004 8玉树州1 0844.371 129.474.3081.76±0.713.941.041 987.0730.010 2小计 3 73815.073 957.9215.08276.19±0.3113.321.058 887.0730.004 8阿里地区5 96424.055 274.3820.10479.30±21.8023.110.884 4237.3740.005 2昌都1 8967.652 420.719.23240.06±13.0411.571.276 7170.9120.011 2拉萨7543.04466.531.7824.90±0.871.200.618 727.4130.010 2林芝3 97916.055 640.4721.50423.26±5.5620.411.417 6167.8410.014 5那曲2 58710.432 698.9510.29203.63±2.939.821.043 365.7420.010 3日喀则4 13516.683 855.6414.69287.54±3.2413.860.932 473.1570.007 4山南地区1 7437.031 923.657.33139.26±1.076.711.103 686.4920.009 5小计 21 05884.9222 280.3384.921 797.94±46.7686.681.058 0237.3740.005 2总计24 796100.0026 238.25100.002 074.13±47.07100.001.058 2237.3740.004 8
4 结 论
(1)截至2014年底,青藏地区有面积≥0.000 2 km2的冰川24 796条,总面积约2.624×104 km2,约占青藏地区总面积的1.37%,冰川储量为2.027~2.121×103 km3,其中分布西藏的冰川数量和面积占84.92%、84.92%,冰储量约占86.68%,基本反映了青藏地区的现代冰川现状。
(2)青藏地区冰川数量以面积<1.0 km2的冰川为主,约占冰川总数量的80.59%;面积以介于1.0~10.0 km2的冰川为主,占冰川总面积的45.59%;其次面积为10~100 km2的冰川,所占比例为30.43%,位居第二;面积>100.0 km2的冰川共7条,其中面积最大的中峰冰川面积为237.37 km2,分布在昆仑山峰的南坡。
(3)青藏地区9座山系(高原)均有冰川分布,海拔5 000~6 500 m之间是冰川集中发育区域,占冰川总面积的85%以上。喜马拉雅山和念青唐古拉山2座山系发育了11 836条冰川,面积为14 652.17 km2,冰储量(1 169.52±25.92) km3,分别占冰川相应总量的47.74%、55.84%和56.39%。冰川数量、面积和冰储量最少的山系是横断山和祁连山,分别占调查区冰川数量、面积和冰储量相应总量的9.41%、6.33%和4.23%。昆仑山冰川平均面积最大,为1.56 km2;冈底斯山冰川平均面积最小,仅为0.34 km2。
(4)青藏地区外流区冰川数量和面积均多于内流区,分别约占冰川总面积的68.13%和31.87%。恒河外流区是冰川分布数量最多、面积和冰储量最大的一级流域;其次是青藏高原内流区。这两个区域的冰川数量、面积、冰储量分别占调查区冰川总量的69.76%、77.50%和81.10%。黄河流域是冰川数量最少、规模最小的一级流域,仅有冰川161条,面积128.24 km2。长江流域冰川平均面积最大,为1.27 km2。印度河流域冰川平均面积最小,为0.44 km2。
(5)冰川在西藏的7个地区(市)和青海的5个州中均有分布,西藏阿里、林芝和日喀则地区冰川数量和面积最多,占青藏地区的50%以上,单条冰川平均面积林芝地区最大,其次是昌都地区,青海省海南州的平均冰川面积最小,仅为0.12 km2;青海无论是冰川数量还是冰川面积都远少于西藏。
(6) 本次调查使用的GF-1影像可清晰识别冰川,色彩层次亦丰富,质量尚好但细节较OLI模糊,总体上能够满足冰川调查需要;与国际上冰川监测常用的ETM或者OLI相比,GF-1有4个波段相对较少,4天的覆盖周期时间分辨率较高,可以作为冰川调查的数据源。
致谢:本文是“青藏冰川变化与冰湖溃决灾害遥感综合调查”二级项目中“青藏地区冰川冰川湖动态变化中分辨率遥感调查”子项目的部分调查成果,谨对参与完成该调查工作的其他人员表示感谢。
参考文献:
[1] 刘时银,姚晓军,郭万钦,等.基于第二次冰川编目的中国冰川现状[J].地理学报,2015,70(1):3-16.
[2] 施雅风.简明中国冰川编目[M].上海:上海科学普及出版社,2005:1-194.
[3] 杨针娘.中国冰川水资源[M]. 兰州:甘肃科学技术出版社,1991:1-158.
[4] 施雅风.中国冰川与环境——现在、过去和未来[M].北京:科学出版社,2000:1-426.
[5] 王宗太, 苏宏超.世界和中国的冰川分布及其水资源意义[J].冰川冻土, 2003, 25(5): 198-202.
[6] 方洪宾,赵福岳,张振德,等.青藏高原现代生态地质环境遥感调查与演变研究[M]. 北京:地质出版社,2009:78-105.
[7] 姚檀栋, 杨志红, 刘景寿.冰芯记录所揭示的青藏高原升温[J].科学通报,1994,39(5):438-441.
[8] 张东启,效存德,秦大河.近几十年来喜马拉雅山冰川变化及其对水资源的影响[J].冰川冻土,2009,31(5):892-895.
[9] 姚檀栋,刘时银,蒲健辰,等.高亚洲冰川的近期退缩及其对西北水资源的影响[J].中国科学(D辑),2004,34(6) :535-543.
[10] 赵刚,杨太保,田洪阵.1990-2011年南天山地区冰川面积变化对气候的响应[J].水土保持研究,2014,21(2):257-268.
[11] 姜珊,杨太保,田洪阵.1973—2010年基于RS和GIS的马兰冰川退缩与气候变化关系研究[J].冰川冻土,2012,34(3):522-529.
[12] 谢自楚,刘潮海.冰川学导论[M].上海:上海科学普及出版社,2010:1-486.
[13] 刘时银,刘潮海,谢自楚.冰川观测与研究方法[M].北京:科学出版社,2012:1-211.
[14] RAUP B,KAAB A,KARGEL J,et al. Remote sensing and GIS technology in the Global Land Ice Measurement for Space (GLIMS) Project[J]. Computers & Geosciences,2007,33(1):104-125.
[15] 李成秀,杨太保,田洪阵.1990-2011 年西昆仑峰区冰川变化的遥感监测[J].地理科学进展,2013,32(4):548-558.
[16] 周文明.基于光学影像的冰川面积参数提取研究[D].长沙:中南大学,2013.
[17] 金珊珊.基于多时相遥感影像和多源DEM的冰川变化研究[D].青岛:山东科技大学,2013.
[18] XIAO C D, LIU S Y, ZHAO L, et al. Observed changes of cryosphere in China over the second half of the 20th century: an overview [J]. Annals of Glaciology, 2007, 46(1): 382-390.
[19] 赵少华,王桥,杨一鹏,等.高分一号卫星环境遥感监测应用示范研究[J].卫星应用,2015(3):38-40.
[20] 安国英,郭兆成,童立强,等.青藏地区冰川、冰川湖动态变化中分辨率遥感调查项目成果报告[R].北京:中国地质调查局国土资源航空物探遥感中心,2017.
[21] 童立强,祁生文,安国英,等.喜马拉雅山地区重大地质灾害遥感调查[M]. 北京:科技出版社,2010:148-180.
[22] 西藏自治区地质矿产局.西藏自治区区域地质志[M].北京:地质出版社,1994:1-20.
[23] 青海省地质矿产局.青海省区域地质志[M].北京:地质出版社,1991:1-5.
[24] 中国遥感数据.Landsat-8卫星数据已开通共享服务[EB/OL].[2013-10-17]. http://rs.ceode.ac.cn/landsat 8 -2013-10-17.jsp.
[25] RADI V, HOCK R. Regional and global volumes of glaciers derived from statistical upscaling of glacier inventory data[J].Journal of Geophysical Research, 2010, 115: F01010, doi: 10.1029/2009JF001373.
[26] GRINSTED A. An estimate of global glacier volume[J]. The Cryosphere, 2013, 7: 141-151.
Remote Sensing Survey of Glaciers Based on GF-1 Spectral Data in the Qinghai-Tibet Region
Abstract: Based on remote sensing images, including mainly GF-1 and minor OLI in 2014/2015, the glacier inventory of Qinghai Province and Tibet Autonomous Region is completed with reference to the Second Chinese Glacier Inventory and other literatures. The results show that there are currently 24,796 glaciers with a total area of about 2.624×104 km2, covering about 1.37% of the land area of the surveyed region.The ice volume is estimated to be about 2.027×103 to 2.121×103 km3. These 24,796 glaciers are dominated by those with sizes below 1.0 km2(19,983 glaciers, 80.59% of the total number of glaciers) or 1.0 to 10.0 km2(11,962.40 km2, 45.59% of the total glacial size). The Zhongfeng Glacier is the largest glacier (237.37 km2). The glaciers are spatially distributed in nine mountains or plateaus: The Nyainqentanglula Mountains contain the largest number of glaciers, followed by the Himalayas and the Gangdese Mountains, which altogether account for 63.33% of the total glacier number in the region.The sizes and ice volumes of the Nyainqentanglula, Himalayas and Kunlun Mountains rank among the top three, accounting for 68.09% of the total size and 73.44% of the total ice volume. The average single glacier size in the Kunlun Mountains and Qiangtang plateau is greater than that in the Nyainqentanglula and the Himalayas. Over 85% of the number or the area of the glaciers is concentrated in the 5,000 m to 6,500 m elevation. Analogously, there is an obvious difference of glacial resources in different drainage basins. On one hand, the Ganges basin is the largest first-level drainage in terms of the number (47% of total) or size (52% of total) of glaciers. On the other hand, the number (21% of total) and size (24% of total) of glaciers in the inner basin of the Qinghai-Tibet Plateau are smaller than those of the Ganges, and the average single glacier size in the inner stream is slightly larger than the average of the outflow basin. In general, the Tibetan glacial resources (in terms of number and size of glacier and ice volume) account for nearly 85% of the total glacial resources in Tibet and Qinghai, and the average single glacial size is similar between the glaciers in Tibet and Qinghai.
Key words:glacier resource;glacier distribution;remote sensing;glacier inventory;Qinghai-Tibet region
中图分类号:P627;P642.5
文献标志码:A
文章编号:1000-8527(2018)03-0584-11
DOI:10.19657/j.geoscience.1000-8527.2018.03.16
收稿日期:2017-10-10;改回日期:2017-12-20;责任编辑:潘令枝。
基金项目:中国地质调查局 “青藏冰川变化与冰湖溃决灾害遥感综合调查”项目(121201203000160012)。
引用格式:安国英,韩磊,黄树春,等.基于GF-1光谱数据的青藏地区冰川资源现状遥感调查[J]. 现代地质, 2018,32(3): 584-594.
AN Guoying, HAN Lei, HUANG Shuchun,et al.Remote Sensing Survey of Glaciers Based on GF-1 Spectral Data in the Qinghai-Tibet Region[J].Geoscience,2018,32(3): 584-594.