巴托拉等冰川的新近变化及对中巴公路的影响

李志杰， 王宁练,3， 常佳雯

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

喀喇昆仑公路（Karakoram Highway，KKH，以下简称公路）北起新疆喀什，南至巴基斯坦塔科特，是中巴两国实现互联互通，建设中巴经济走廊的战略通道。公路穿越的东帕米尔高原、喀喇昆仑山是高亚洲冰川分布最集中的地区之一，冰川灾害曾多次威胁和中断公路通行，造成重大损失［1-3］。位于巴基斯坦北部罕萨河（洪扎河）中游的巴托拉段是整个公路沿线距离冰川最近、冰川灾害最集中的路段，其中以巴托拉（Batura）、帕苏（Pasu）、固尔金（Ghulkin）等3 条冰川对公路的威胁最大也最为直接［4-6］。

1973年，巴托拉冰川主排水道（巴托拉河）改道引发的冰川洪水冲毁了建设中的公路和公路桥。1974—1975 年，为提供合理的修复设计，中国考察组对巴托拉冰川开展了包括末端进退、流速、厚度、融水径流等全方面、详尽的科学考察［4］。最终考察组提出按照原线修复的科学建议，确保了公路的长期稳定通行［7］。此后的1978 年、1980 年、1994 年等，考察组再次对巴托拉地区开展科学考察，积累了丰富的地质地貌、冰川物理、水文、气象等科学资料［7-8］。在此基础上，学者们揭示了巴托拉、帕苏、固尔金等冰川在20世纪近百年的进退变化，并预测巴托拉冰川将于2010s前后再次发生前进［4-5,7-8］。

进入21 世纪以来，随着气候波动的加剧，巴托拉、帕苏、固尔金等冰川发生了许多新的变化，这些变化对公路的正常通行产生了新的威胁［9］。例如巴托拉冰川的主排水洞位置不断北移，最终导致融水通道被冰舌挤压切断，时隔近50 a 后再次引发巴托拉河的大改道；帕苏冰川末端的冰碛阻塞湖频繁溃决，2008 年1 月和4 月2 次引发冰湖溃决洪水，中断公路通行［10］；固尔金冰川末端的频繁波动导致冰崖下的排水通道不断迁移，仅在2008年就引发4次泥石流，致使公路桥底严重淤积［6,10-11］。因此，鉴于灾害风险的日渐凸显和公路无可取代的地位，有必要对巴托拉、帕苏、固尔金冰川的历史与新近变化进行揭示和分析，为公路的安全运行和灾害防治提供参考。

1 研究区概况

巴托拉、帕苏、固尔金冰川发源于喀喇昆仑山西段主脊——洪扎喀喇昆仑山北坡（图1）。洪扎喀喇昆仑山平均海拔超过6000 m，最高峰巴托拉峰海拔高达7795 m，为冰川发育提供了极为有利的地形条件，使得巴托拉冰川群成为了公路沿线规模最大、末端海拔最低、进退最为频繁的冰川群之一。巴托拉、帕苏、固尔金冰川上游粒雪盆面积开阔，补给物质充足，中下游则深入了干热河谷，消融强烈，因此属于典型的多温型冰川，并表现出较高的物质平衡水平［4］。本研究的解译结果表明，2021 年巴托拉冰川长约56.5 km，面积307.12±4.32 km2，末端海拔约2620 m；帕苏冰川长约26.4 km，面积64.81±0.72 km2，末端海拔约2700 m；固尔金冰川长约17.8 km，面积28.26±0.47 km2，末端海拔约2540 m。

图1 巴托拉、帕苏、固尔金冰川的地理位置Fig.1 Geographic location of the Batura,Pasu,and Ghulkin Glacier

2 数据与方法

2.1 数据源

2.1.1 Landsat 影像 本研究使用了获取于1972—2022 年的25 景Landsat 系列影像，以揭示研究冰川的相关变化（表1）。所用影像均为L1T 级产品，由USGS（United States Geological Survey，https://earthex⁃plorer.usgs.gov/）进行了系统辐射校正和几何校正，并结合DEM进行了地形校正。由于Landsat系列影像有相当高的正射校正精度，达到1/2 像元左右，因此本研究直接利用Landsat 影像开展冰川变化研究［12］。为提高冰川边界、出水洞、融水径流等地物的解译或识别精度，本研究选取了消融季末期，冰川区无云、山体阴影覆盖较小的影像［13］。

表1 遥感影像数据Tab.1 List of the Landsat images used in this study

2.1.2 冰面高程变化数据集 Brun等［14］利用ASTER立体像对提取DEM，基于大地测量学方法，通过DEM的空间匹配与误差校正，重建了2000—2016年高亚洲冰川的表面高程变化与物质平衡，并生成了高亚洲冰面高程变化数据集。数据集中的误差主要源于DEM的生成精度及其空间匹配偏差，因此以非冰川区等地形稳定区的高程残差来校正并评价冰面高程变化的误差［14］。该数据集的空间分辨率为30 m，可获取于PANGAEA(https://doi.pangaea.de/10.1594/PANGAEA.876545)。

2.1.3 冰面流速数据集 NASA 的ITS_LIVE（Intermission time series of land ice velocity and elevation）产品包含了1985—2020 年高亚洲主要冰川的流速数据。该产品基于Landsat 影像，采用Gardner 等提出的自动裂缝特征跟踪处理链方法，通过进行局部归一化、过采样、特征跟踪等提取冰川流速［15-16］。影像的空间匹配偏差是该产品的主要误差来源，尽管在生成过程中已基于基岩等稳定地形对其进行校正，但无法彻底消除。因此，ITS_LIVE产品同样基于非冰川区的流速残差评价冰川流速的不确定性［15］。1999年之前（TM影像）冰川流速估算的不确定性超过10 m·a-1，1999年之后（ETM+影像）降低为约2 m·a-1，2013 年之后（OLI 影像）降低为约0.5 m·a-1［16-17］。本研究选用了其中1985—2018 年的冰川年平均流速合成数据，空间分辨率为120 m（https://its-live.jpl.nasa.gov/）。

2.2 方法

2.2.1 冰川边界解译 由于研究冰川表面存在大范围表碛覆盖，且进退变化频繁，本研究采用目视解译提取冰川边界［13］。在中上游裸冰区，直接根据Landsat 假彩色影像（热红外、红、绿波段组合）中显著的色彩差异进行冰川边界的数字化（图1）［18］。在表碛覆盖的冰川下游，则根据末端冰崖与出水洞的组合、高程变化等地形和水文特征识别冰川边界［19］。假彩色影像中，裸冰区与表碛区、表碛区与基岩区、冰崖与水流等的色彩差异均得到了显著增强，由此可准确识别冰舌、出水洞等的位置。野外考察资料表明，研究冰川冰舌消融强烈，因此末端冰崖与出水口的组合即代表了冰川运动的最前缘［4,6,20］。在实际解译中，本研究首先完成2021年的冰川边界提取。在此基础上结合不同时期的影像，修改冰舌等发生变化的部位，以获取其他年份的冰川边界。2.2.2 解译精度评价 基于光学影像解译冰川边界的精度主要受影像空间分辨率的影响，因此本研究采用统计冰川轮廓线经过的像元数量的方法，评价目视解译提取冰川边界的不确定性［21］。公式如下：

E=N·λ/2

式中：E为冰川面积的不确定性；N为冰川边界经过的像元数量；λ为像元面积。

2.2.3 出水洞与融水通道识别 出水洞指位于冰川末端，冰内和冰下水系的排水出口，本研究结合野外考察记录和遥感影像识别冰川出水洞与融水通道的迁移变化。如图2所示，在消融季末期的Land⁃sat 和Google Earth 影像中，得益于巴托拉冰川较大的排水量，冰川出水洞与融水通道的准确位置可获得有效识别。此外，在1973—2021 年2 次大改道期间，巴托拉河下游保持稳定，上游随冰舌波动而延伸，因此结合不同时期的遥感影像即可准确识别出水洞与融水通道的变迁过程。

图2 巴托拉冰川出水洞与融水径流的识别Fig.2 Identification of drainage system of the Batura Glacier

3 结果与分析

3.1 冰川进退变化

近150 a 来，巴托拉冰川发生过多次进退波动。1873年巴托拉冰川发生前进，但此后近50 a间冰舌位置基本保持稳定。1930s 冰川开始衰退，1944—1966年退缩827 m。1960s末转入前进，1966—1980年前进约136 m。1980s初再次衰退，至1994年累计后退了150—200 m［4-5,7,22-23］。影像资料表明（图3a），1994—2008 年巴托拉冰舌继续后退约410 m，2008年冰舌再次前进，至2021 年7 月累计前进约230 m。综上所述，近150 a 来巴托拉冰舌每隔50～60 a就会发生一轮前进，但每轮前进一般只可维持10～20 a，总体而言，近150 a 来巴托拉冰川以退缩为主导［23］。

图3 1972—2021年巴托拉、帕苏、固尔金冰川进退变化Fig.3 Terminus positions at different periods for the Batura,Pasu,and Ghulkin Glacier from 1972 to 2021

帕苏冰川虽曾于1907—1913 年发生过大规模跃动，导致冰舌前进约1200 m［24］，但此后却长期处于退缩状态，仅在2000 年前后发生过一次前进。1972—1998年，帕苏冰舌累计后退约180 m，1998—2002 年前进220 m，2002—2021 年再次后退约570 m。因此近百年来，帕苏冰川总体处于加速退缩状态（图3b）。

近百年来，固尔金冰川进退频繁，堪称公路沿线最活跃的冰川之一。1885—1980年固尔金冰川发生了3 轮进退，期间还伴随着多次小幅度波动［6］。1966—1978 年，固尔金南冰舌后退260 m，1978—1980年却又前进260 m，此后前进速率减缓，1980—1994 年累计前进100 m［6］。之后再度衰退，至2021年累计退缩约390 m。北冰舌更加活跃（图3c），1972—1990 年后退200 m，1990—1998 年又前进170 m。之后，冰舌开始在小范围内频繁波动，最终在1998—2018 年累计前进120 m。2018 年起冰舌再次衰退，至2021年已后退约36 m。

3.2 融水通道变化

1973年6—7月，巴托拉河改道后冰川主出水洞由中流线附近转移至冰舌南侧，但此后主出水洞位置开始缓慢迁移。2000年之前迁移速率较慢，基本稳定在冰舌南侧。2000 年之后，迁移速率明显加快，逐步越过中流线到达冰舌北侧。至2021 年6月，主出水洞已处于冰舌西北侧，排水道随之延伸完全横亘于冰舌前（图4a）。2021年6—7月，冰舌的突然前进挤压切断了主排水道，引发了近50 a 来巴托拉河的又一次大改道（图4b）。由于冰坝阻挡，融水径流被迫在冰碛丘陵中侵蚀切割出一条新的通道，直接导致巴托拉河由2.9 km缩短为2.5 km，纵比降由65‰增大至76‰，因此流速显著增快。

图4 巴托拉、固尔金冰川融水径流排水通道变化Fig.4 Variations of the meltwater runoff drainage channel of the Batura and Ghulkin Glacier

受地形影响，帕苏冰川的排水通道十分稳定。但在1981年之后，冰舌附近在融水补给下形成了一个小型终碛阻塞湖，2011年实测面积约0.12 km2［10］。在水量增加、冰舌崩塌、谷坡落石等的共同作用下，阻塞湖曾于2007年7月、2008年1月、4月等数次发生溃决，引发冰湖溃决洪水［3,10］。但经过数次溃决以及融水径流的长期侵蚀，近年来帕苏冰湖终碛堤的完整性已遭到破坏，冰湖发育规模受到很大限制。

固尔金段是公路沿线受冰川融水影响最大的路段之一。由于冰舌分为南北两支且融水径流流程极短，不足1 km，因此冰川融水呈辐散状，极不稳定且频繁发生改道。固尔金冰川主要排水通道约有6 条，小型通道则难以计数［6,10］。其中，北冰舌融水集中在M1通道，南冰舌的M2通道也已基本废弃（图4c）。长期以来，M5为南冰舌的主排水道，但近年来M4 的水量也在增大。M3 曾在2014 年前后短暂成为主排水道，但之后水量又迅速减少，M6 水量也较为有限。因此，近年来南冰舌融水通道以M5与M4为主，M3与M6为辅。

3.3 冰面高程与流速

近年来“喀喇昆仑异常”现象备受关注，但就处于喀喇昆仑山西北边缘的巴托拉地区而言，2000年以来的冰川物质亏损却有所加速［20,25］。基于KH-9 DEM和SRTM DEM的研究发现，1973—1999年巴托拉和帕苏冰川物质平衡分别为0.00±0.10 m w.e.a-1和0.05±0.11 m w.e.a-1，表明2000年以前冰川物质收支基本保持稳定［20］。但此后的2000—2016年，巴托拉、帕苏、固尔金冰川表面高程均发生了下降（图5a），下降速率分别为0.09±0.29 m·a-1、0.15±0.17 m·a-1和0.08±0.26 m·a-1，尤以冰川末端表现最为明显［14］。尽管下降速率较为微弱，但这依旧表明2000年以来巴托拉地区冰川发生了物质亏损。

图5 巴托拉、帕苏、固尔金冰川表面高程变化（a）与年平均流速（b）Fig.5 Surface elevation change(a)and annual average flow velocity(b)of the Batura,Pasu,and Ghulkin Glacier

长期以来，3 条研究冰川均维持着较高的表面流速。1985—2018年，巴托拉、帕苏、固尔金冰川的平均流速分别为59.60±0.28 m·a-1、106.44±0.44 m·a-1和19.50±0.35 m·a-1。就流速的空间分布而言（图5b），高值区域均出现在冰川上游的粒雪线附近。巴托拉冰川的最大流速出现在粒雪盆出口，第一、第二冰流构成的冰瀑布处，接近1200 m·a-1。帕苏冰川的最大流速出现在相同的位置，约940 m·a-1。固尔金冰川流速相对平缓，最大流速也出现在了冰瀑布以下的附近位置。相较于巴托拉和固尔金冰川，帕苏冰川流速较快，意味着其物质平衡水平更高，这是帕苏冰舌消融更强烈的主要原因。

4 讨论

4.1 冰川变化的影响因素

山地冰川是气候与地形共同作用的产物，但在十年或百年尺度上，冰川变化主要受气候条件的影响［26］。冰川物质平衡直接由当年气候状态决定，末端进退却更多地受气候长周期波动的影响，因为从气候变化到冰舌波动需要一个能量积聚和动力反映的过程，即存在滞后性。滞后期的长短主要与冰川规模有关，小规模冰川为3～30 a，大规模冰川则会更长［4,27］。

就巴托拉、帕苏、固尔金这种大型复式山地冰川而言，气温的长周期波动是影响冰舌进退的关键因素［8］。如图5所示，2000年以来，冰川积累区基本保持稳定，消融区高程则有所降低，表明近几十年来罕萨河谷的升温直接导致了冰川物质亏损的加剧［28-29］。而巴托拉、帕苏、固尔金冰舌在减薄的同时发生进退波动，则充分体现了冰舌进退对气候响应的滞后性以及冰川规模对于滞后性的影响。有研究表明，近150 a来巴托拉冰舌进退相较于气温波动存在60 a 左右的滞后期，冰川的数次前进均与气温的阶段性变化存在很好的对应关系［4,7-8］。基于这一判断，不仅成功预测了巴托拉冰舌在1970s—2000年的进退，还成功预测了冰舌在2010s 的最新一轮的前进。

跃动冰川指周期性发生快速运动的冰川，当跃动发生时其速度往往是平时的数倍至上百倍［9］。有研究将巴托拉第一、第二冰流等部位的高流速、大量裂隙与弧拱发育视为跃动的表现［24］，但实际上这对于大型山地冰川而言并不罕见。对巴托拉这种巨型冰川而言，支冰流跃动对末端进退的影响也十分有限，如克拉牙依拉克、费德琴科等冰川支冰流的跃动并未改变整体的退缩趋势［30-31］。帕苏冰川曾于1907—1913 年，固尔金冰川曾于1913—1925 年、1966—1978 年发生跃动［24］，但近50 a 来2 条冰川的积累区并未发生显著增厚，因此在近期发生大规模跃动的可能性较低。

4.2 融水径流的影响因素

较高的物质平衡水平导致3条冰川的融水径流量明显高于其他地区的同等规模冰川，巴托拉冰川甚至高于近710 km2的费德琴科冰川［4］。巴托拉河多年平均流量为40.1 m3·s-1，瞬时最大流量可超过600 m3·s-1，其中近80%的径流量集中在5—9 月，近60%集中在7—8 月［4,32］。巴托拉河90%的水量来源于冰川融水，因此径流量主要受气温控制，随气温上升而增加，在夏季降水条件下有时甚至会因气温下降而减少［4］。此外，即便在7—8月巴托拉河的径流量波动也十分剧烈，常随极端高温天气的出现而陡涨［32］。帕苏、固尔金冰川融水径流的年内变化与巴托拉河类似，但变化过程远不及前者剧烈。

巴托拉河径流年际变化较为平缓，主要由于：（1）干旱年份降水偏少而融水偏多，湿润年份降水偏多而融水偏少，起到了天然的调节作用［31］。（2）冰舌区年平均气温接近10 ℃，导致消融十分剧烈，存在热融喀斯特现象，因此冰内裂隙、孔穴，冰碛物的滞水能力也能起到一定的径流调节作用［4,32］。（3）1970s 以来巴托拉流域冰川基本保持稳定，即便2000年以来物质亏损有所增加，但对融水径流的影响依旧有限。实际上不仅是巴托拉河，1980s 以来整个罕萨河流域的年径流量也无明显变化［33］。

巴托拉河在1973 年和2021 年的2 次改道有一定的相似性，改道均发生在6—7月融水径流快速上升的时期、改道前主排水洞均已发生长期缓慢的迁移、改道的发生均与冰舌突然前进有关，因此2次改道的原因也是基本相同的。如图1 所示，巴托拉冰川的物质来源极不对称，冰舌北侧主要由尤克苏戈兹、北、西冰流等构成、南侧主要由第一、第二冰流构成。由于南侧冰流拥有更大的高差与更开阔的粒雪盆，因此物质供给量更大，直接导致冰舌南侧厚度与流速均高于北侧［4,34］。冰川中下游，南侧弧拱与裸冰的延伸范围明显大于北侧就是二者存在流速差的直观表现。在这种流速差的长期作用下，主排水洞受到推挤而缓慢北移，排水道随之延长并横亘于冰舌前，最终被突然前进的冰舌切断而被迫改道。因此，冰舌南北两侧物质供给量的差异是巴托拉河改道的根本原因，冰舌的突然前进则是直接原因［4］。

4.3 对公路的影响

本研究在公路沿线选取了4 个稳定的参考点，分别为巴托拉桥、帕苏桥、M1桥、M4桥，以参考点与冰舌的最近距离来衡量公路与冰川的距离。如图6所示，近50 a来巴托拉冰舌与公路的距离总体上在增大，当前尚有1 km 以上，因此在短期内不会对公路产生直接威胁。帕苏冰舌自2000 年以来不断加速后退，与公路的距离也在增大，因此短期内也不会直接威胁公路。尽管固尔金南北冰舌仅在小范围内波动，但由于距离过近，且在将来不能排除发生跃动的可能，因此对公路的直接威胁将长期存在。

图6 1972—2021年巴托拉、帕苏、固尔金冰舌与公路的距离变化Fig.6 Variations of the distances from glacier tongues of the Batura,Pasu,and Ghulkin Glacier to the Karakoram highway from 1972 to 2021

巴托拉冰川对公路更现实的威胁在于巴托拉河新的大改道。2021年6—7月改道后，巴托拉河流程缩短，比降增大，下切侵蚀能力增强，近1 a年来，新河道明显拓宽加深（图7）。更关键的是新河道下游折向南流，与公路基本平行，因此水流将在转弯处直接冲击路基。甚至完全可以推测，正是由于路基坚固才迫使水流南流返回故道。巴托拉河具有流量大、径流高度集中、陡涨陡落等特点［32］，因此在未来的强烈消融季中，随着水量陡涨，超过600 m3·s-1的水流裹挟冰块砂石奔腾下泄，将对路基产生强烈侵蚀，对公路构成重大威胁。

图7 2022年5月27日巴托拉冰舌区概况Fig.7 Overview of Batura Glacier tongue on May 27,2022

帕苏冰川近年来消融强烈导致融水径流增大，但由于冰湖发育规模受限，大规模冰湖溃决洪水对于下游公路的威胁趋于减弱。固尔金冰川融水轻则在路面积水结冰，影响通行，重则侵蚀路基、淤积桥底，甚至引发泥石流洪水。但在2008年采取工程措施后，固尔金段的融水灾害已得到了有效治理［3］。然而，固尔金冰川融水径流的迁移摆动问题仍不容忽视，一方面它可能会导致过水工程废弃，但更重要的是径流量可能会超过单一通道的排水能力，进而引发新的灾害。

5 结论

中巴喀喇昆仑公路巴托拉段受冰川灾害威胁最为严重，是影响公路畅通运行的关键路段。基于野外考察与遥感影像等资料，分析了巴托拉、帕苏、固尔金冰川近百年来的变化、变化原因及对中巴公路的影响，主要结论如下：

（1）近百年来，特别是2000 年以来，巴托拉、帕苏、固尔金冰川总体处于微弱退缩状态，冰舌发生后退减薄或基本保持稳定，对公路的直接威胁趋于减弱。

（2）在冰川消融加剧与冰舌波动的共同影响下，巴托拉段冰川融水径流频繁迁移改道，冰川融水洪水、冰湖溃决洪水、冰川泥石流等次生灾害对公路的威胁正日渐凸显。

（3）巴托拉段冰川变化对公路最直接的威胁在于2021 年6—7 月巴托拉河新的改道，改道后水流将直接冲击侵蚀路基，影像公路安全通行。