喜马拉雅山脉呈近EW向弧形展布于西藏高原南侧, 是地球上最年轻和最雄伟的褶皱山系, 其同位素年龄只有10~20 Ma。喜马拉雅山脉随青藏高原自新近纪开始大幅度抬升以来, 形成了世界上最年轻而又最高的高原^{[1, 2, 3, 4]}, 发育了许多规模巨大的现代冰川, 冰斗、角峰、刃脊、冰碛垅岗、冰碛丘陵和冰碛台地等冰川地貌广泛发育^{[5, 6, 7, 8, 9, 10, 11]}。

冰川湖(简称“ 冰湖” )多分布于现代冰川的前缘、侧缘以及古冰川谷、古冰斗内, 分布位置海拔高, 湖面面积较小。在地貌上, 冰湖需要有山间谷地、冰斗、古冰川谷等积水洼地才能形成。涉及现代冰川的有冰面湖、冰内湖、冰坝湖; 涉及冰川退缩的有冰川终碛湖、冰蚀槽谷-冰碛湖、冰斗湖、冰川侧碛阻塞支谷湖等^[12]。多数冰湖后缘与现代冰川相连或距离现代冰川的冰舌较近, 在现代冰川前进或跃动、冰舌断裂、冰湖岸坡出现崩塌或滑坡、温度骤然增加导致冰川融化加速、湖口溯源侵蚀加剧、坝体下部管涌引起塌陷等诸多可能因素的影响下, 易造成冰湖溃决灾害, 其中部分会演化发展, 引发泥石流而使灾害进一步扩大^{[12, 13]}。

根据笔者的遥感调查结果, 分布在喜马拉雅山脉、面积大于0.014 km²的冰湖大约有1 680个, 约占西藏全区冰湖总数的67.7%。冰川终碛湖占冰湖总数的40%~45%, 其中面积大于1 km²的冰川终碛湖只有12个, 最大的是康马县冲巴雍水库, 水面面积约11 km²(人工筑坝); 其次是面积为5.7 km²的康马县桑旺错。据迄今为止所了解到的情况, 喜马拉雅山的冰湖溃决都是冰川终碛湖溃决^{[14, 15]}。冰川终碛湖受形成条件和自然环境的影响, 往往会发生溃决, 造成洪水和泥石流灾害, 严重威胁人类生产和生活、生存与发展, 成为喜马拉雅山地区一重大地质灾害隐患^{[12, 13, 14, 15, 16, 17]}。

从20世纪30年代中期起, 很多学者就对西藏冰湖溃决进行过研究, 例如崔鹏等^[17]、李吉均^[18]、徐道明等^{[19, 20]}和李德基等^[21]。1999年吕儒仁等系统总结了1995年以前西藏发生的12次冰湖溃决事件, 并对每次事件进行了详细描述^[12]。但因研究区自然环境极度恶劣、高寒缺氧、地形切割极大, 大部分地区人迹罕至, 常规地面调查工作难以开展, 导致地质研究调查程度极低; 加之调查技术手段落后, 对冰湖溃决的研究仍限于溃决事件发生后才能开展, 不能及时指出冰湖溃决隐患, 难以满足防灾减灾要求。

本文发挥遥感技术优势, 采用QuickBird和ETM等卫星遥感数据, 结合野外考察, 研究了喜马拉雅山地区冰湖溃决灾害隐患, 并对其影响因素进行了分析。研究结果对西藏地区冰湖溃决灾害的预警与防治具有重要意义。

经研究发现, 喜马拉雅山地区的冰湖在区域分布上具有由东南向西北逐步递减的趋势^[14]。冰湖分布最多的地区是山南、日喀则、林芝和阿里地区, 分布密度最大的地区是山南地区的洛扎一带。从实地考察和文献资料来看, 高度在海拔5 000 m以上的冰湖约占冰湖总数的50%~60%, 海拔4 000~5 000 m的冰湖占40%左右, 海拔低于4 000 m的冰湖仅占0.5%~2%。海洋性冰川分布区的冰湖的海拔高度普遍低于大陆性冰川分布区的冰湖。同一地区的冰斗湖的海拔常常高于冰碛湖(冰斗湖在雪线附近, 冰碛湖则在雪线以下)。在地理位置上, 林芝地区冰湖的高度明显低于日喀则和阿里地区冰湖的高度(高差甚至能达到1 500 m)。由于林芝地区地形切割强烈, 多处于峡谷区, 冰湖的数量多, 但规模一般较小; 日喀则地区冰湖的数量不多, 但规模较大。

冰面湖、冰内湖及冰坝湖不仅体积都非常小, 而且分布在大型现代冰川体上, 湖深很浅, 形态很不稳定。分布数量较多、规模较大的是冰川终碛湖, 其次是冰斗湖和侧碛阻塞支谷湖。在冰川融化季节, 冰川终碛湖和侧碛阻塞支谷湖的湖水多数是浑浊的, 湖水中的泥沙来自于冰川运动沿程所携带的岩石风化碎屑物质, 湖面上漂浮着大小不同的冰块。发生溃决最多的冰湖是冰川终碛湖, 多沿冰川谷分布, 一般呈长条状; 其后缘距现代冰川冰舌距离为0~数百m, 近年来随着全球气候的变暖, 这种距离的大部分在进一步增加^[15]。

由于冰湖的补给来源主要是冰雪融水, 其分布往往与冰川分布区域紧密相连, 尤其是冰碛湖。冰川在运动过程中不仅强烈掏蚀下垫冰床和沟床, 而且作用于沿程的2侧谷坡, 不断改造槽谷形态, 并在不同位置形成不同类型的冰湖。冰川的大小与数量取决于积雪区面积, 积雪区的面积愈大, 冰川的规模就愈大, 汇集到冰舌地段和冰湖内的消融水量也就越多, 冰湖的规模也就越大。同样的消融水补给面积, 海洋性冰川分布区的冰川和冰湖规模要大于大陆性冰川分布区, 其主要原因是海洋性冰川具有补给快、消融快的特点^[14]。本次遥感调查发现, 一般补给区的面积只有大于2 km²才能形成具一定规模的冰川和冰湖。冰湖分布的数量与冰川分布的数量成正比, 冰湖的规模大小与冰川的大小成正相关, 目前发生溃决的冰湖补给区面积都大于3 km²(多数在2~40 km²之间)。

喜马拉雅山脉的措美以东为海洋性冰川分布区。该段山脉的南北2侧地形皆十分陡峭, 部分地段形成再生冰川; 冰湖分布的海拔有高有低, 海拔最低的冰湖高程在3 500 m, 最高可达5 500 m; 冰斗湖分布在山脉上部, 冰碛湖分布在山脉下部。喜马拉雅山脉的措美— 吉隆段为海洋性冰川与大陆性冰川的过渡区, 降水由东向西逐步递减, 植被也由森林覆盖区变成稀疏植被覆盖区; 分布位置最低的冰湖高程从3 500 m增加到4 200 m, 最高可达5 800 m; 该段也是冰湖分布最多的区段。喜马拉雅山脉的吉隆以西为大陆性冰川分布区。由于降水明显逐步递减, 该段山脉南北坡的植被发育差异很大, 北侧基本无植被, 雪线高度大于5 000 m, 冰湖数量也由东向西大幅度递减; 冰湖分布位置的最低高度从4 200 m增加到4 500 m, 最高可达5 500 m。

冰川终碛湖的补给来源有冰川融化补给、降雨直接补给、降雪直接补给、沟谷地表径流补给、地下水补给和其他冰川湖补给。冰川终碛湖主要是通过天然泄水口进行排泄, 也存在地下渗漏和蒸发排泄。天然径流排泄的泄水口往往分布在终碛垄上, 由于天然终碛垄堆积相对较紧密, 有一定的宽度, 不均一, 渗透性相对较差, 一般情况下的稳定性较好; 但终碛垄多数起伏较大, 最凹的地方就变成了天然泄水口, 当水位上升较快、最低的天然泄水口达不到泄流要求时, 就会在较低的凹口形成新的天然泄水口。因此, 1个冰碛湖通常有1个到多个天然泄水口, 用以满足正常的排泄。还有1种排泄方式是非正常突然超量排泄, 这种排泄即为冰湖溃决排泄。

中、高分辨率的ETM和QuickBird等卫星遥感图像忠实地纪录了冰湖溃决形成的印迹。冰湖溃决过程中形成的水体变化和地形地貌改变都可以通过遥感图像的纹理、色调、形态等影像特征表达出来。从遥感图像中的典型影像特征分析和实地踏勘验证来看, 已发生过溃决的冰湖一般都会在沟口部位有明显的扇状堆积体, 在遥感影像中呈现扇状暗色调区域, 堆积体有较强的浮雕般凸起感, 花纹呈斑点状, 花纹结构较为粗糙, 色调不均一; 冰湖的湖岸有阶地; 终碛上有明显的溃口和“ U” 型冲沟, 冰湖下游的老冰碛物上也有明显的“ U” 型冲沟; 冰湖下游河道冰碛物和坡积物发育, 物源丰富, 河道多呈明显的“ U” 型冲沟。图1示出冲巴吓错溃决冰湖的ETM和QuickBird卫星遥感影像图, 从图像中可以明显地看到溃决冰湖的影像特征。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 冲巴吓错冰湖ETM和QuickBird影像图Fig.1 ETM and QuickBird images covering Chongbaxia glacier lake

冰湖溃决灾害隐患是指冰湖存在溃决的危险性, 并有灾害风险。其判断依据是冰湖周围有高陡的边坡、破碎的岩体, 冰湖上游有高位冰川, 冰湖内有足量的湖水存储, 冰湖下游有明确的致灾体。据此, 本次遥感调查最终确定喜马拉雅山地区存在19个溃决灾害隐患冰湖, 其中13个存在重大溃决灾害隐患, 6个存在较大溃决灾害隐患(表1、图2)。

表1

Tab.1

表1(Tab.1)

表1 喜马拉雅山地区冰湖溃决灾害隐患 Tab.1 Glacier lakes burst potential in Himalayan area 序号 冰川湖名称 所在地区 湖面海 拔/m 冰川湖面 积/万m2 成灾 风险 序号 冰川湖名称 所在地区 湖面海 拔/m 冰川湖面 积/万m2 成灾 风险 1 扎日朗普错 洛扎县 5 421 39.0 重大 11 曲典错 岗巴县 5 270 51.6 较大 2 桑旺错 康马县 5 125 569.6 重大 12 直习错 定结县 5 350 113.7 重大 3 白湖 康马县 4 924 134.6 重大 13 金错 定结县 5 365 47.3 重大 4 冲巴吓错 康马县 4 912 82.4 重大 14 虾错 定结县 5 410 72.7 重大 5 钟勤错 康马县 4 898 22.5 较大 15 玛朗错 定日县 5 340 384.0 重大 6 错嘎布 康马县 4 898 20.6 重大 16 查玛曲旦登错 聂拉木县 5 420 51.2 重大 7 藏玛桑错 康马县 5 294 41.6 较大 17 帕曲错 聂拉木县 5 310 54.9 重大 8 养莫可朗错 亚东县 4 990 14.0 较大 18 郭如错 聂拉木县 5 277 436.0 重大 9 朗错 亚东县 5 000 58.3 较大 19 郭骆强错 聂拉木县 5 350 506.4 重大 10 倬莫拉日神湖 亚东县 5 210 55.3 较大

表1 喜马拉雅山地区冰湖溃决灾害隐患 Tab.1 Glacier lakes burst potential in Himalayan area

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 喜马拉雅山地区冰湖溃决灾害隐患分布图Fig.2 Distribution of glacier lakes burst potential in Himalayan area

以冲巴吓错溃决隐患为例, 对喜马拉雅山地区冰川终碛湖溃决灾害隐患进行分析。

冲巴吓错位于康马县萨马达乡雅拉山口北坡, 湖面面积82.4万m², 湖口高程4 912 m, 终碛高程4 938 m。冲巴吓错的遥感影像特征显示, 在终碛下有明显的“ 扇状” 堆积体, 堆积体的色调饱和度较高, 色调均一, 植被不发育。冰川湖周围植被很少。湖中见冰川碎块, 湖水浑浊, 湖后缘为陡崖(如图1所示)。终碛上有明显的“ U” 型冲沟, 下游河道有明显过水冲刷痕迹。湖口下游的沟谷为“ U” 型谷, 冰碛物发育, 物源丰富。湖口下游4.5 km处为冲巴芒错, 河床平均比降48‰ ; 下游14.4 km处为冲巴水库, 河床平均比降24‰ 。冰川前缘为一陡崖, 进入冰川湖后, 平均坡度13.6° , 山峰高程5 571 m。

本次遥感调查的结果(图3)显示, 1990年11月14日―2000年11月17日期间, 冲巴吓错曾发生过溃决事件, 前端湖水后退320 m, 最大湖宽从500 m减至365 m; 溃口深约30 m, 宽约35 m(图4), 总洪水量约375万m³。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 冲巴吓错冰湖遥感解译图Fig.3 Remote sensing interpretation for Chongbaxia glacier lake

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 冲巴吓错2000年溃决口实景照片Fig.4 Field photo for burst outlet of Chongbaxia glacier glake in 2000

2008年, 笔者对冲巴吓错进行了实地调查与走访, 得知该湖于2000年8月6日下午因后缘冰川滑坡引发涌浪溃坝洪水和泥石流, 使冰湖下游的萨马达、康马、少岗和南尼等4个乡受灾, 共冲毁农田222亩（1亩=666.7 m²）、草场8 160亩、林地80亩、水渠4 600 m、公路328 m、桥梁10座、护河堤140 m、羊圈7个。从遥感影像及实地调查中看到, 冲巴吓错坝体2侧高山耸立, 岩体破碎, 崩塌、落石很多。据此, 判断该湖有再次溃决的高危险; 但因现在的溢流口与终碛最低高程基本一样, 湖水较浅, 故预测其灾害风险会低于2000年溃决形成的灾害。

本文采用遥感方法调查了喜马拉雅山地区冰湖分布情况, 并以冲巴吓错为例, 对冰川终碛湖溃决灾害隐患进行了分析; 在此基础上, 对冰川终碛湖溃决隐患影响因素进行了研究。得出以下结论:

1)喜马拉雅山地区的冰湖在区域分布上具有由东南向西北逐步递减的趋势。冰湖分布最多的地区是山南、日喀则、林芝和阿里地区, 分布密度最大的地区是山南地区的洛扎一带。林芝地区地形切割强烈, 多处于峡谷区, 冰湖的数量多, 但规模一般较小。日喀则地区冰川湖的数量不多, 但规模较大。

2)喜马拉雅山地区存在19个溃决灾害隐患的冰湖均是冰川终碛湖, 其中13个存在重大溃决灾害隐患, 6个存在较大溃决灾害隐患。

3)引起冰川终碛湖发生溃决灾害的条件和因素较多。其中, 冰湖周围的冰崩、雪崩、基岩崩塌和滑坡入湖涌浪溃坝等是最重要的诱发因素。