人类活动与山地冰冻圈的交互日益增加加剧了冰冻圈灾害风险。陡峭的山坡和高地形起伏及地震活动使山区容易发生极具破坏性的物质运动。冰川和冻土对气候变化的敏感性正在加剧这些物质运动和灾害过程^［1］。青藏高原及周边地区是除南、北两极地区之外全球最重要的冰川资源富集地^［2］。作为“亚洲水塔”的重要组成部分，青藏高原及周边地区冰川对亚洲地区水资源压力的缓解具有重要意义，特别为我国西部地区的水资源安全、生态安全和经济社会发展提供重要保障^［2-3］。随着气候变暖的持续和人类活动的增加，青藏高原及周边地区冰川整体处于快速消融状态，不仅影响冰川水资源储备，增大了冰川自身的不稳定性，加剧了灾害发生风险^［4-5］。

在全球气候变暖大背景下，青藏高原过去相对稳定的冰川已经变得越来越不稳定，有的还出现了剧烈运动和强烈消融，进而导致各种冰川灾害的发生^［4］。冰崩等以往罕见的冰冻圈灾害时有发生。2016年和2018年，在西藏西北部阿里和东南部林芝先后发生了大规模冰崩事件。这两种类型完全不同的冰川连续发生冰崩灾害事件，很可能说明目前青藏高原的冰川整体上处于不稳定状态，灾害发生概率与空间范围正逐渐增大^［4-5］。冰川灾害有时还引发次生灾害，形成一个从冰冻圈开始，继而影响到岩石圈、水圈、生物圈、人类圈的灾害链，延长和放大灾害后果^［6］。

青藏高原冰雪快速消融和冰冻圈灾害风险不断加大的情况下，及时有效地监测和了解冰冻圈灾害发生地点、规模和动态变化成为制定防灾减灾决策和采取应对措施的首要任务和前提。然而，冰冻圈灾害区域一般地处高寒、高海拔偏远地区，经济欠发达，基础设施薄弱，抵御和应对自然灾害的能力弱，难以用常规的地面观测来完成灾害监测和调查，如青藏高原冰冻圈灾害多发生在海拔4 000 m以上的高山区，使得遥感技术成为冰冻圈灾害监测最为有效手段。尤其是，以卫星遥感为主的监测方法已在冰崩、冰川跃动和雪崩等冰冻圈灾害监测、应急响应和风险评价中得到了广泛应用，已成为全球山地冰冻圈灾害监测的主要手段。

李壮等^［7］通过实地调查，结合Landsat等多期遥感影像以及InSAR雷达数据，对色东普沟高位地质灾害发育特征进行分析，结果表明2018年的两次灾害事件物源来自于1处高位岩崩区和3处高位冰崩区。柴波等^［8］利用高分卫星资料和谷歌地球影像结合模型对西藏聂拉木县嘉龙错冰湖溃决危险性进行了评价，结果表明冰崩和雪崩是高山区冰湖溃决的重要致灾因子，嘉龙错上方冰川存在发生滑坡的可能性。汤明高等^［9］利用Google Earth、GF-2和ZY-3高分辨率遥感影像系统研究了青藏高原冰崩隐患数量、类型、发育规律及危险性，以期为应对青藏高原冰崩灾害及防灾减灾提供科学支撑。2002年9月发生在高加索山脉的Kolka冰崩是一起罕见的巨型冰崩灾害，冰崩体积达130×10⁶ m³，造成了130多人遇难^［10］。Kääb等^［11］利用ASTER卫星遥感数据通过提取数字高程等数据对Kolka冰崩进行了分析，重建了此次冰崩的动力学过程。Pandey^［12］利用谷歌地球高分辨率卫星影像对印度喜马拉雅喜马偕尔的岩石冰川进行第一次综合编目，结果表明研究区共有516个岩石冰川，所占面积大约为353 km²，其中59%的以冰川为源头，而41%为岩屑为源头。相关更多的研究工作关注基于Landsat等卫星遥感影像分析亚洲高山区中部冰川的跃动特征，从而为冰川跃动灾害监测和防治提供重要支撑^［13-15］。Kääb等^［16］建议综合多种遥感技术对冰崩等冰冻圈灾害过程的潜在威胁进行评估，能够从一定程度上减少冰川灾害造成的损失。

在全球变暖趋势下，冰崩、冰川跃动、雪崩、冰湖溃决等冰冻圈灾害在青藏高原及其周边地区呈逐渐增多趋势，受灾范围也在不断扩大^［4-7］。研究表明，在青藏高原，西藏的冰崩隐患点数量最多，占整个青藏高原的近60%，且大部分位于喜马拉雅山脉和念青唐古拉山及其东延高山^［9］。因此，对未来青藏高原特别是西藏冰冻圈灾害监测和应对措施提出了更高的要求。利用高分辨率卫星遥感图像对冰崩高危险区的冰川进行有针对性的动态监测，进而开展冰崩灾害风险评估和预警，是防灾减灾和减轻灾害损失的重要途径。Sentinel-2相比Landsat 等卫星遥感数据无论从空间分辨率还是从重访周期、波谱分辨率和图像宽幅都实现了革命性跨越。卫星传感器有专门用于冰雪监测和云雪识别的短波红外波段，使得Sentinel-2在冰雪和冰冻圈灾害监测中具有其他高分辨率卫星遥感数据无可比拟的优势。

为此，研究利用Sentinel-2卫星遥感数据监测了2016年西藏阿里阿汝错和2018年西藏林芝色东普的冰崩灾害发生、发展全过程，展示了Sentinel-2在冰冻圈灾害监测和应急响应中的巨大应用潜力。研究工作对于青藏高原冰冻圈灾害监测和应急响应以及全球其他类似的灾害监测和应急服务具有重要的借鉴意义。

欧盟委员会和欧洲航天局（ESA）于2003年共同推出地球环境与安全监测计划（Global Monitoring for Environment and Security，GMES）。该计划的主要目标是通过对欧洲及部分非欧洲国家现有和未来发射的卫星以及现场观测数据进行协调管理和集成，实现环境与安全的实时动态监测，为决策者提供数据，以帮助制定环境法案，同时对诸如自然灾害和人道主义危机等紧急状况作出响应，保证欧洲的可持续发展和提升国际竞争力。2012年12月，欧盟委员会将该计划更名为“哥白尼计划”。哥白尼计划是迄今最为宏大的地球观测计划，它提供精确、及时和方便获取的信息以改进环境的管理，了解和减轻气候变化的影响，确保民众安全。

Sentinel系列是哥白尼计划的核心组成部分，主要对陆地和海洋环境变化进行连续监测，为自然灾害预警及灾后应急提供帮助^［17-21］。Sentinel-2是一个极轨多光谱高分辨率光学成像卫星组，是到目前为止全球唯一能够免费获取的空间和光谱分辨率最高、重访周期最短、幅宽最大的光学遥感卫星。目前，Sentinel-2 由A、B两颗极地轨道卫星组成，分别于北京时间2015年6月23日和2017年3月7日发射升空，位于离地面786 km高度的太阳同步轨道上，彼此位置相距180°。Sentinel-2 通过A、B两颗卫星组成星座协同观测，大大提高了卫星对地观测周期，每5 d就能提供一次从地球56º S到83º N范围陆地表面影像数据。今后C星和D星陆续投入运行后，卫星重访周期提高一倍，其在环境和灾害监测、应急管理响应时间和时效性提升一倍。Sentinel-2已经成为哥白尼计划卫星系列中应用最广泛的卫星，60%的数据下载量来自Sentinel-2卫星。Sentinel-2卫星遥感数据已在全球环境变化、农作物生长、土地覆盖变化制图、自然灾害（如洪水、林火、滑坡、地震和火山喷发等）监测、应急管理及救援等方面得到了广泛应用。

Sentinel-2卫星的波段设计源自其设计阶段与广大用户广泛讨论的结果，同时继承了Landsat和SPOT卫星的成功经验^［22-24］。卫星上搭载了一个采用推扫式成像模式的多光谱成像仪（MSI），是新一代多光谱光学成像仪，共有13个波段，其中波段4（红）、3（绿）、2（蓝）和近红外波段8的地面分辨率为10 m，3个植被红边波段5、6、7和两个短波红外波段11、12 及窄近红外波段8 a的地面分辨率为20 m，其余波段1、9和10的地面分辨率为60 m，详见表1。

Sentinel-2卫星影像采用了与Landsat 8影像一致的UTM投影，提供L1C和L2A两种产品，其中L1C是经过传感器辐射定标和几何配准的大气上层表观反射率产品，L2A则是经过云检测和大气校正后的标准地表反射率产品。用户可以从欧空局（ESA）官网下载Sentinel-2大气校正插件Sen2Cor对L1C产品进行大气校正，生产L2A级别产品。Sentinle-2 L1C产品在UTM投影基础上按美国军用格网系统（MGRS）进行分幅，将每条UTM投影带按南北向细分为20个区域，再进一步分割为109.8 km×109.8 km的格网。用户直接可以从哥白尼计划官方网站（https：∥scihub.copernicus.eu/dhus/#/home）免费下载L1C和L2A产品。ESA还提供卫星遥感数据处理和分析工具SNAP，方便用户显示和处理Sentinel和Landsat系列等卫星遥感数据。

归一化积雪指数NDSI（Normalized Difference Snow Index）是卫星遥感定量识别冰雪的主要方法。计算NDSI的Sentinel-2 MSI传感器对应波段是绿波段3和短波红外波段11，计算公式如下：

Hall等^［25-27］发现，当NDSI大于等于0.40时能够较好地识别地表积雪，并且能够区分大部分云。在此基础上，利用水体在近红外波段反射率低而雪在该波段具有较高反射率的光谱特性，近红外波段的反射率大于0.11消除水体对积雪判识的影响^［28］。ESA提供的Sentinel-2大气校正插件Sen2Cor中有积雪检测算法。该算法中NDSI的积雪识别阈值下限为0.2，同时近红外波段8a的反射率大于0.15时，该像元识别为积雪像元。

在阿汝错冰崩体识别方法上，在ENVI5.6遥感图像处理软件中除了采用上述的NDSI计算方法之外，还用了可见光单波段阈值设定和目视解译相结合的方法对冰崩体面积进行定量计算；对色东普沟高位地质灾害影响范围识别上则以目视解译方法为主。

主要技术路线是，首先从哥白尼计划官方网站下载研究区L1C或L2A图像。若用于目视解译则下载10 m分辨率的L1C数据即可，如果用于计算NDSI等定量应用则需要下载L2A产品，但图像的分辨率也会降至20 m。下载好的L1C图像直接可以用ESA SNAP软件中打开，利用该工具可以进行Band4-3-2真彩色等波段合成和图像增强，以期达到最佳显示效果，最后导出后直接可以在ArcGIS 10.0中进行目视解译和制图等。

2016年7月17日，西藏阿里地区日土县阿汝错湖区的53号冰川发生了崩塌，冰崩体携带大量碎屑物质向偏东方向滑动约7 km后快速冲入阿汝错，产生20 m高的湖啸，导致湖泊水位上涨9 m，阿汝错湖水向对岸前推了约240 m。此次冰崩造成当地9名牧民遇难，数百头牲畜被埋，大量草场被毁^［29］。9月21日，同一冰川群的50号冰川再次发生冰崩^{［4-5，29-30］}。如此大规模的冰川崩塌事件连续发生在一直以来冰川活动较为稳定的青藏高原内陆地区，极为罕见^［29-31］。冰崩发生后，冰崩体的后期融化使其下游湖美玛错水位上涨迅速，2016—2019年湖泊水位升高3 m，冰体融化对湖泊扩张的贡献率达到23%^［32］。

Sentinel-2获取的卫星遥感图像很好地记录了此次冰崩事件发生的过程。图1（a）是53号冰川发生崩塌4 d后于2016年7月21日Sentinel-2卫星首次观测到的冰崩灾害实况卫星遥感影像图。经测量，冰崩堆积体长6.80 km、宽2.45km、面积为9.67 km²，平均厚度约7.5 m^［5］，冰崩体积约为68×10⁶ m^3［30］，冰崩扇淹没湖泊最长距离675 m，淹没湖泊面积0.75 km²。

2016年7月21日至2017年的春季由于研究区为冷季气温较低的时段，冰崩体消融速度较为缓慢（图2），2017年5月14日的冰崩体面积仍有7.01 km²，冰崩体面积减少了27.5%。这一期间冰崩体面积波动是由于春季降雪导致冰崩体被积雪覆盖引起的。研究发现，冰崩体的消融过程并非是像冰川那样从外围到里的逐渐萎缩消融过程，而是冰崩体不同区域破碎化的消融过程（图1（c））。2017年5月14日至9月26日是冰崩体快速消融的阶段，与研究区年内气温最高时段一致。2017年5月14日至7月8日消融达到2.98 km²（图1（d）），消融面积达到每10 d 0.54 km²，7月8日至8月15日共计减少面积2.10 km²，每10 d的消融面积达到0.55 km²。8月至9月的消融速度出现了明显减少，为0.21 km²/10 d，不及7月至8月消融速度的一半。2017年9月26日冰崩体面积出现冰崩发生以来的最小，仅1.06 km²。Sentinel-2卫星图像显示，之后残留的冰崩体一直存在，而且被积雪覆盖之后面积略有波动（图2）。2018年6月8日冰崩体面积降至1 km²以下，为0.95 km²，7月8日仅剩下0.14 km²（图1（e））。可以看出，53号冰崩残留体最终于2018年7月正式全部消融，从冰崩发生至完全消融先后持续了2 a整的时间。

同样，阿汝错50号冰崩发生于2016年9月21日，Sentinel-2监测到的最近一次晴空图像是2016年10月9日的图像（图1（b））。该图监测显示，50号冰崩体长5.80 km、宽1.90 km、面积为6.75 km²，平均厚度超过30 m^［5］，估算的冰崩体约为83×10⁶ m^3［30］。大部分冰崩体滑向偏东的阿汝错方向，但未到达阿汝错湖岸，距湖岸最近距离约1 km。少部分冰崩体滑向东北方向，面积0.84 km²，占50号冰崩体面积的12%。

2016年10月9日之后，冰崩体进入了逐渐消融和面积不断减少的过程（图2）。至2017年4月24日面积减少至5.78 km²，这一阶段由于湖区气温较低，不利于冰体的消融，面积减少速度缓慢（图1（c））。2017年5月两景图像上，由于积雪覆盖出现了冰崩体面积略有增加的情况。从图2中可以看出，5月中旬至8月中旬冰崩体消融速度明显加快，但是相比53号冰崩体其消融速度非常缓慢，主要原因是50号冰崩体很厚，达到30 m，需要更长的时间融化。冰崩残留体较薄的条件下，气温上升引起冰崩体的破碎化，加速其消融。11月之后消融速度出现了明显减少，直至下一年的春季，冰崩体面积减少幅度不大。2017年11月8日冰崩体面积为4.13 km²，过了5个月之后，2018年4月2日冰崩体面积仍保留有3.89 km²，仅减少了0.24 km²。2018年7月31日的面积为1.74 km²，之后冰崩体消融和面积减少也非常缓慢，2019年7月3日和2020年7月25日分别是1.18 km²和0.76 km²，至2021年6月22日，50号冰崩体的残留堆积体面积仍有0.58 km²（图1（f））。

2022年2月16日，对阿汝错53号和50号冰崩发生区域行了实地调查。地面调查同样表明，53号冰崩体中的冰体部分已全部消融，现场基本是大小岩石和碎石等冰崩残留体，个别大的石头直径在1 m以上（图3（a））。从现场可以看出，53号冰川是在山体斜坡上的冰川在重力作用下冰川前部断裂解体后发生了冰崩。由于山体坡度较大，冰川发生断裂后快速向东冲向阿汝错，约7 km之后最后冲入湖中。现场还发现，冰川发生断裂后53号冰川的后半部分保留完好，而且随着冰川物质的累积，冰川在不断发育，今后仍存在再次发生冰崩的可能性（图3（b））。50号冰崩体现场整体上与此类似，冰崩体中的冰体部分已基本消融，现场覆盖物基本上是由岩石、碎石和泥土等组成，个别大石头直径也在1 m以上。

从现场调查可以看出，阿汝错50号冰川的冰崩更有可能是一次冰川垮塌，是在山体斜坡上的冰川在重力作用下发生了垮塌，进而发生了冰崩（图3（c））。由于50号冰川的山体坡度小于53号冰川，冰川发生崩塌后向东冲向阿汝错方向，但冰崩体运动速度要小于53号冰崩体，所以冰崩体最终未到达湖泊。其冰崩体面积也小于53号冰崩体，两者相差2.92 km²，但是其估算的冰崩体积大于53号冰崩体15×10⁶m³，表明该冰崩体厚度远大于53号冰崩体，使得崩塌体后部至今仍有未消融的冰体。

冰崩是指冰川突然崩解而垮塌的现象^［2］。西藏冰川作为青藏高原和亚洲高山区冰川的主要组成部分，历史上冰崩灾害事件时有发生。例如，1940年亚东上康布地区的冰川突然滑动，1954年年楚河上游桑旺错后方冰川的滑动，造成了400多人死亡，卡若拉川曾先后于1963年7月、1975年9月和1978年8月都发生过冰崩，巨大冰块下滑到公路路面上，1950年8月藏东南南迦巴瓦峰则隆弄冰川因地震引发巨大冰崩，一度造成雅鲁藏布江断流^［33］。这些冰崩都曾引起冰川洪水、泥石流、公路和河流堵塞等，造成很大的生命和财产损失，但是这些冰川的冰舌所在地形部位很陡，几乎是在悬冰舌的情况下发生了灾难性的冰川崩塌和滑动。相反，阿汝错冰崩是一起历史罕见的低角度山地冰川的大规模冰崩灾害事件，冰崩体积合计达151×10⁶ m³。之前，2002年报道的体积为130×10⁶ m³的高加索山脉Kolka冰崩被认为是全球唯一的大量冰崩流动和跃动相伴的灾难性冰川崩塌事件^［34］。无论是低角度山地冰川还是悬冰川，冰川崩塌一般都发生在盛夏高温多雨时期。

2016年在阿汝错流域发生的两次冰崩灾害性事件迅速成为国内外焦点，国际上许多知名科学家纷纷研究此次冰崩发生的原因^［5，30］。冰崩是冰川不稳定的表现，主要发生因素包括气候、地形、冰体热力状况、基岩不稳定性及地震活动等^［4］。近年来，全球气候变暖逐渐被发现可能是引发冰崩灾害的主要深层次原因^［4-5］。Kääb等^［30］通过卫星遥感、数值模拟和野外调查对阿汝错冰崩事件的起因进行了系统研究，结果表明阿汝错两次冰崩事件不存在单一的触发因子，是多种因素综合作用的结果，是气候、天气驱动的外力强迫作用在特定的多温和软性基岩的冰川特性上引起的，是一起历史罕见的低角度冰川的巨大灾难性不稳定事件^［30］。

色东普沟位于西藏米林县雅鲁藏布江下游左岸加拉白垒峰下方，属于念青唐古拉山脉东延高山峡谷区。流域面积67 km²，近南北向发育，上游地形宽阔陡峭，冰川发育，中下游地形狭窄，呈阶梯状“Ｖ”字型。流域最高点为加拉白垒峰，海拔高度7 294 m，最低点为沟口与雅鲁藏布江交汇处，海拔2 746 m，高差4 548 m^［7，35］。色东普沟主沟长约7.6 km，入江处沟口宽度230 m。

色东普沟地处藏东南高山峡谷地带，由于受气温、强降水、地震以及特殊地形地貌条件的综合影响，该流域地质灾害频发，历史上发生过多起高位地质灾害，多次造成雅鲁藏布江堵江事件。Sentinel-2自发射以来记录了色东普沟冰崩—岩崩—碎屑流—堰塞湖—溃决洪水灾害链的发生发展过程。

色东普沟及其周边第一张晴空Sentinel-2卫星图像接收日期是2015年12月6日（图4（a））。该卫星图像上显示，色东普沟内除了主沟、冰川下部冰碛物及冰雪覆盖区之外的堆积体上都长有植被，主沟中部阶地的高地上也长有植被，水流从高地两旁通过。这些地表特征足以说明，色东普沟在过去较长一段时间内相对稳定，未发生过大的碎屑流事件。由于历史上曾发生过碎屑流堵塞河道事件，江内一直有堵江残留体造成的江心滩，其长约1.7 km，宽439 m，面积39.3×10⁴ m²。此处河道宽646 m，68.0%的河道处于堵塞状态，江水两侧过流，以左侧河道为主。已有记录和卫星遥感图像解译结果表明，1984年之前大规模堵江事件后水流从河道左侧自然漫顶过流，1984—2013年相对稳定，未发生大规模堵江事件，2014年色东普冰崩造成大规模堵江事件，最后水流仍从左侧自然漫顶过流^［35］。

2016年3月5日的Sentinel-2图像显示，河道堆积区由于水位下降而面积稍有变大之外，河道总体形态和水流状况与2015年12月6日基本一致（图5（a））。2016年5月24日，河道总体形态、江心滩面积和水流状况等与2015年12月6日基本一致，见图4（b）。可见，2015—2016 年色东普沟较稳定，未出现大规模碎屑流和堵江情况。

2017年11月5日的Sentinel-2图像显示，原有的江心滩基本被淹没，露出两处小块江心滩，合计面积仅3.9×10⁴ m²，同时在沟口出现了5.2×10⁴ m²的冲积扇，表明之前沟内出现了大规模碎屑流事件（图5（b））。据文献记录，2017 年10 月22 日色东普沟口出现堵江，崩滑冲击过程曾引起强烈地面震动，水位上涨约30 m后坝体自然溃决。此次堵江很可能是由加拉白垒峰西侧冰川整体断裂后形成的，是一起冰崩—岩崩—碎屑流造成的冰崩—堵溃链式灾害^［35-36］。

2017年12月30日的Sentinel-2卫星图像显示，沟口出现了56.4×10⁴ m²的大型冲积扇，水流从原来的河道左侧改为右侧通过（图4（c））。说明在此之前色东普沟发生过大规模碎屑流事件。根据文献，2017年12月21日色东普沟发生冰崩事件，导致冰崩—岩崩—碎屑物再次堵塞雅鲁藏布江形成堰塞湖，堵江72 h后堰塞体自然溃决，河流从河道右侧通过。对比事件前后Sentinel-2遥感影像显示，冰崩是触发此次滑坡事件的主要因素，也是一起典型的冰崩—堵溃链式灾害。此次冰崩与2017 年11月18日的林芝6.5级地震密切相关。该地震对流域造成严重扰动，破坏了冰川、冰碛物和沟谷堆积物的完整性和稳定性，变得更加松散，为之后的冰崩—碎屑流事件提供了环境条件和丰富的物质基础^［36］。

2018年夏季6月8日的Sentinel-2图像更为直观地显示了2017年底多次冰崩碎屑流和两次堵江事件之后色东普流域全景及沟口堵塞处的清晰情景（图4（d））。可以看出，从冰雪覆盖下端到沟口段整个沟道被冰崩—岩崩—碎屑流刮铲、深切和水流侵蚀痕迹清晰，原有沟内堆积体上的植被全部消失，成为裸露的地表，沟道两侧被碎屑流刮铲后裸露地表明显增大。相比2017年底，沟口冲积扇面积减少至43.6×10⁴ m²，右侧河道明显加宽，水流更为畅通。

2018年10月17日色东普沟发生高位冰崩—岩崩—碎屑流，崩滑体冲出沟口并堵塞雅鲁藏布江，形成堰塞湖，历时56 h后堰塞坝体溃决，水流通过堰塞湖右岸自然漫顶过流。针对此次冰崩—碎屑流发生过程的研究表明，色东普沟后部加拉白垒峰西北偏西侧海拔6 000 m左右近山顶部位发生冰崩岩崩，在重力作用下沿陡峻斜坡向南南西崩滑，断裂解体后迅速滑塌铲刮沟道并裹挟碎屑物转向南南东冲击，沿途再次解体，快速奔流，铲刮沟底和两侧，最终冲出沟口堵塞雅鲁藏布江，形成堰塞湖^［35］。冰崩岩崩最高点至沟口总长10.2 km。10月29日色东普沟再次发生冰崩—碎屑流—泥石流，冲出沟口，堆积在已有堰塞体之上，导致再次堵塞雅鲁藏布江，24 h后江水自然过顶漫流^［7，36］。针对灾害发生区域山高谷深，人烟稀少，交通困难的实际，提出了适应自然、全面避让和适当疏导的防灾减灾对策^［35］。两次冰崩—碎屑流—堰塞湖—溃决洪水灾害链造成20多个村庄被淹，近6 000 人受灾。2018年10月31日获取的Sentinel-2卫星图像很好地记录了此次大规模冰崩碎屑流事件。原有的沟口冲积扇被新的碎屑物覆盖后雅鲁藏布江被完全堵塞，沟口冲积扇总面积达72.4×10⁴ m²，堰塞湖水位上涨后沟口水域面积明显加大，流域后部中间高位山沟右侧南南西向的冰川岩石崩滑区极为明显（图4（e））。

2018年11月30日的Sentinel-2卫星图像显示了雅鲁藏布江堰塞湖自然溃决1个月之后的色东普沟沟口情景（图4（f））。河流水位显著下降后水流两侧漏出大面积冲积扇，面积达95.2×10⁴ m²，河流从沟口大型冲积扇右侧通过。河道处于严重堵塞状态，最窄处河道宽不到50 m。

2015年12月6日和2018年6月8日两张遥感影像（图4（a）和图4（d））对比可以发现，图4（d）上色东普沟中后段堆积体上的植被几乎全部消失，沟道下段两侧山体被刮铲的面积明显增大，沟口出现了大面积的冲积扇，2015年沟口对岸残留的洪积扇全部淹没消失，水流从原来的左侧改道为右侧。2018年10月17日和29日的两次大规模冰崩—岩崩—碎屑流使得沟内刮铲侵蚀下切程度和沟口堆积体面积进一步增大，沟口的堵塞程度进一步加剧，水流宽度进一步收窄。

2017年和2018年4次冰崩—碎屑流事件不仅对流域造成了重大灾害，对色东普沟沟口水流状况和地貌形态产生了重大影响。特别是2017年11月18日林芝6.5级地震之后发生的2017年12月21日冰崩—碎屑流对流域影响最大，沟口大量堆积体导致雅鲁藏布江水流长期从河道左侧改道从右侧流过，河道处于严重堵塞状态，最窄处不到50 m。

色东普沟冰川属于海洋型冰川，冰川物质累积和消融比大陆型冰川更为剧烈，气候变暖引起冰雪消融，在局地强降水和地震等因素作用在色东普沟有利地形地貌条件下，经常引发冰崩及相关灾害链^［37］。同时，色东普沟属于藏东南高山高寒地区，日夜温差大，沟谷后部冰雪所在高山冻融作用强烈，碎裂松动岩体分布较广，沟内时常发生岩崩^［38］。特别是，色东普沟特殊地形地貌下的“漏斗效应” 放大了冰川融化和集中降雨对松散堆积体的影响， 加剧了沟内滑坡的孕育、形成及后期演化^［39］。因此，在气温上升、强降水、地形地貌和地震等因素的综合作用下，色东普沟冰崩灾害呈多发频发态势，且具有一定的周期性，今后很长的时间内仍会多次发生^［40］。

除了此次发生在喜马拉雅东段的冰崩—岩崩—碎屑流灾害事件之外，2021年2月7日，喜马拉雅西段查莫里地区发生了特大规模岩崩—冰崩—泥石流事件，导致超过200人死亡或失踪和两座水电站被毁的重大灾害链事件^［1］。在喜马拉雅东西两侧的灾害事件中，高山冰雪和山体岩石的高位崩塌是灾害触发的直接原因。喜马拉雅西段南坡下游人口密集、水电开发集中，灾害链末端造成的损失更为严重^［41-46］。色东普沟相较于查莫里流域不存在水电站等大型基础设施，加之人类居住和活动有限，造成的灾害损失有限。然而，喜马拉雅流域高位崩塌造成的重大灾害链事件凸显了高寒山区冰冻圈灾害监测预警和可持续发展理念的极端重要性。

本研究利用Sentinel-2卫星遥感数据重现了西藏阿里和林芝的冰崩灾害发生全过程，展现了其在青藏高原冰冻圈灾害监测领域的应用情况和前景。得出的主要结论和建议如下：

（1） 阿汝错53号冰崩是在山体斜坡上的冰川在重力作用下冰川前部断裂解体后发生了冰崩，而且随着冰川物质的累积和发育，今后仍有可能再次发生；阿汝错50号冰崩则是一起典型的冰川前部垮塌产生的冰崩。

（2） Sentinel-2监测显示，阿汝错53号冰崩体最终于2018年7月全部消融，从崩塌至最后消融先后持续了2 a的时间，而50号冰崩体厚度远大于53号冰崩体，消融速度缓慢，2021年6月22日仍有0.58 km²的冰崩残留体。2017年和2018年色东普沟大规模冰崩—岩崩—碎屑流事件不仅对色东普沟地貌地表景观和河口水流状况产生了重要影响，对上下游流域造成了重大灾害。

（3） 在全球变暖持续和青藏高原环境总体变暖变湿的大背景下，高原冰冻圈灾害趋于活跃，及时监测和掌握这些灾害发生和变化特点是防灾减灾的第一步。喜马拉雅流域东西两次高位崩塌造成的重大灾害链事件凸显了山区灾害监测和可持续发展理念的重要性。因此，建议系统开展青藏高原和喜马拉雅高山区灾害风险调查，加强遥感和地面监测，逐步建立监测预警体系。未来，藏东南念青唐古拉山及其东延高山和喜马拉雅山区仍是青藏高原冰崩、雪崩和冰湖溃决等冰冻圈灾害发生的重点区域，应给予高度关注，进一步加强监测调查研究和预警预测，降低灾害风险，提高山区可持续发展能力。

（4） 目前国产GF-1等高分系列高分辨率卫星遥感传感器波段数较少，基本上是红绿蓝和近红外4波段加上全色波段，使得在地表植被参数和积雪等冰冻圈要素的定量遥感监测中的应用受到了限制。所以，建议未来国产高分辨率卫星遥感传感器中除了增加短波红外波段，收窄近红外波段宽度外，还需要增加气溶胶、水汽吸收和卷云检测等波段，从而用传感器自身波段观测值直接进行大气校正，进而更好地应用于地表植被和积雪等陆地生态系统参数的定量化遥感监测和反演。

（5） 在全球环境和灾害监测以及应急响应中除了目前可免费获取的Sentinel-2高分辨率卫星遥感数据之外，利用商业卫星提供的遥感数据也是重要的选项。其中，美国Planet公司拥有全世界最大的遥感卫星群，为用户提供快速更新的、定制化的高分辨率商用卫星影像。公司提出“利用空间帮助地球生命”的口号，其使命是“通过每日对地球所有陆地进行成像，使全球变化可见、可达和可行动”。Planet公司在世界上首次研发成功微卫星群技术，已拥有200多颗卫星，是世界上唯一具有全球高分辨率、高频次、全覆盖能力的遥感卫星系统。它们组成卫星群，保证在世界任何陆地区域每天获取3—5 m分辨率的影像，实现每天监测全球一次，对灾害监测和评估、应急救灾管理等具有不可比拟的优势。