泸定Ms 6.8地震对海螺沟冰川的影响应急分析 Emergency Analysis of the Impact of the Luding Ms 6.8 Earthquake on Hailuogou Glacier 1 first-author 李为乐，博士，教授，主要从事地质灾害隐患识别与监测研究。liweile08@mail.cdut.edu.cn 李为乐，博士，教授，主要从事地质灾害隐患识别与监测研究。liweile08@mail.cdut.edu.cn liweile08@mail.cdut.edu.cn 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 corresp 许强，博士，教授。xq@cdut.edu.cn 许强，博士，教授。xq@cdut.edu.cn xq@cdut.edu.cn 成都理工大学地质灾害防治与地质环境保护国家重点实验室，四川 成都，610059 成都理工大学地质灾害防治与地质环境保护国家重点实验室，四川 成都，610059 State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China 2022-09-05四川省甘孜州泸定县发生Ms 6.8地震，震中距离海螺沟冰川约10 km，海螺沟冰川受此次地震影响如何，受到社会广泛关切。利用多时相Landsat 8和Sentinel-2光学卫星影像，基于归一化雪覆盖指数对海螺沟冰川面积进行监测，并利用Sentinel-1雷达卫星影像，基于像素偏移追踪技术对海螺沟地震前和震中位移进行监测。结果表明，2016―2022年海螺沟冰川8月份面积呈波动变化趋势，与日平均最高气温具有负相关性，冰川运动速度与地形坡度和日平均最高气温具有正相关性。泸定地震未造成海螺沟冰川运动速度大范围显著增大，但对冰瀑布前缘区域扰动明显。震后冰崩直接致灾可能性较小，但增加了泥石流致灾风险。 Objectives 2 On 5th September 2022，an Ms 6.8 earthquake struck Luding County, Ganzi Prefecture, Sichuan Province, China, with the epicenter about 10 km away from Hailuogou Glacier. The affect of this earthquake on Hailuogou Glacier has been widely concerned by the society. Methods 2 Firstly, the glacier area is monitored based on normalized difference snow index using multi-temporal Landsat 8 and Sentinel-2 optical satellite images. Secondly, Sentinel-1 synthetic aperture radar satellite image is used to monitor the displacement before and during the Luding earthquake based on pixel offset tracking technology. Results 2 The results indicate that the area of Hailuogou Glacier shows a shaking trend from 2016 to 2022, which have a negative correlation with the daily average maximum temperature. While the velocity has a positively correlation with the slope gradient and the daily average maximum temperature. The Luding earthquake did not cause a significant increase in the velocity of Hailuogou Glaciers in a large range, but it significantly disturbed the front area of the ice waterfall. Conclusions 2 The possibility of direct disaster caused by ice avalanches after the earthquake was low, but which increased the risk of disaster caused by debris flow. 泸定地震 海螺沟冰川 遥感监测 像素偏移追踪 归一化雪覆盖指数 应急分析 Luding earthquake Hailuogou Glacier remote sensing monitoring pixel offset tracking normalized difference snow index emergency analysis 2022-09-05 12时52分，四川省甘孜藏族自治州泸定县境内发生了Ms 6.8地震，震中位于泸定县磨西镇海螺沟冰川森林公园内（29.59°N，102.08°E），震源深度约16 km，地震持续时间约20 s，最高地震烈度达IX度。此次地震为主震-余震型地震，截至9月8日15时共记录到Ms 2.8及以上余震20次，最大为9月7日发生在石棉县田湾乡的Ms 4.5余震 ［ 1 1 - 2 2 1-2 ］ 。此次地震造成了严重的人员伤亡和财产损失，国务院抗震救灾指挥部办公室、应急管理部紧急启动了国家地震应急二级响应。截至9月13日国家二级应急响应结束，地震及次生地质灾害共造成93人死亡，25人失联，270余人受伤 ［ 2 2 ］ 。同时，地震还触发了大量的崩塌、滑坡等地质灾害，造成大量居民房屋和道路被毁，王岗坪乡、草科乡等地一度成为孤岛，部分区域道路短期内未能抢修通行。 卫星和航空遥感技术在此次地震应急响应中发挥了重要作用。地震发生后，相关部门紧急协调了国产雷达卫星（高分三号、巢湖一号、海丝一号、齐鲁一号、泰景四号等）、光学卫星（高分一号、高分二号、高分六号、高分七号、资源一号、资源三号、北京二号、北京三号、四维02号等）、高光谱卫星（珠海一号）、夜光卫星（吉林一号）等多源卫星资源对地震灾区进行拍摄。由于震后灾区主要为云雾天气，应急响应前期获取到有效数据的主要是雷达卫星。9月5日晚便成功获取了震中区域3 m分辨率的高分三号雷达卫星数据，各单位主要利用其强度影像进行次生地质灾害、震毁房屋和道路等灾情的目视判译。笔者团队利用该影像在海螺沟景区内识别出9处大型滑坡隐患，第一时间提供给现场地质队伍进行排查，经与后期无人机航拍影像对比，识别出的滑坡隐患均有变形迹象或已在地震中失稳，证明了高分辨率雷达卫星强度影像在识别大型滑坡隐患的可行性。直到9月10日高分六号、高分二号、北京三号等光学卫星才获取到有效的影像数据，相关单位迅速利用这些影像完成了重灾区同震地质灾害、损毁房屋、道路等的解译。其中，笔者团队利用高分六号和高分二号解译同震地质灾害2 717处，地质灾害总面积约23 km 2 。泸定地震地貌和地质概况如 图1 图1 所示，其中断层数据来源于四川省地震局提供的1∶50万地质图，地震烈度来源于应急管理部，地震数据来源于中国地震台网中心。此外，自2022-09-06开始，相关单位陆续利用无人机航空平台先后获取了磨西镇、湾东河沟口、大渡河石棉县王岗坪乡-新民乡段、石棉县草科乡等重灾区高分辨率影像，笔者所在团队第一时间利用这些影像完成了同震地质灾害、损毁房屋、道路等的解译，先后向抗震救灾应急指挥部、应急管理部、自然资源部等部门提交解译简报6份，为抗震救灾工作提供了及时支撑。 此次地震震中位于贡嘎山海螺沟冰川森林公园内，距离海螺沟冰川直线距离约10 km，海螺沟冰川受此次地震影响如何，受到社会广泛关切。为此，笔者团队第一时间收集了海螺沟冰川区域地震前后多时相光学和雷达卫星影像，通过归一化雪覆盖指数（normalized difference snow index，NDSI）计算、像素偏移追踪（pixel offset tracking， POT）、目视判译等方法，对海螺沟冰川地震前后面积变化、强形变区分布、运动速度变化等进行分析，为震后潜在冰崩等灾害风险评价提供参考。 1　研究区概况 1 1.1　区域地质背景 2 研究区地处青藏高原东南缘的横断山脉，属典型高山峡谷地貌，海拔高程最低约980 m，最高约7 550 m，海拔高差达6 570 m。大渡河自北向南穿过研究区，河水落差大，水电资源丰富，目前已修建多座大型水利发电站。 本次泸定地震区域位于鲜水河断裂带、龙门山断裂带和安宁河断裂带构成的川西“Y字形”断裂带的交汇处附近（ 图1 图1 ），初步判断此次地震发震断裂为鲜水河断裂带东南段的磨西断裂。鲜水河断裂带为左旋走滑断裂带，全长约350 km，是南北活动构造带的重要组成部分，也是中内陆地震活动性最强的断裂之一，自1700年以来已发生7级以上地震9次 ［ 3 3 - 9 9 3-9 ］ 。在此次地震发生前，已有学者预测未来30 a鲜水河断裂磨西段有发生Ms 7.0强震的危险 ［ 3 3 - 4 4 3-4 ］ 。 研究区地层岩性主要为早震旦系花岗岩、第三纪二长花岗岩、二叠系雷口坡组白云岩等，受构造长期挤压和风化作用，研究区岩土体破碎，是中国典型的地质灾害高易发区 ［ 1 1 - 2 2 1-2 ］ 。 1.2　海螺沟概况 2 海螺沟流域位于四川省甘孜州泸定县磨西镇贡嘎山东坡，面积约195 km 2 ，主沟长约15 km，沟口高程约1 500 m，流域顶部海拔7 556 m，高差高达6 056 m。此次地震在流域内触发了大量崩塌、滑坡地质灾害，造成通往海螺沟冰川的唯一道路多处被掩埋，直至2022-09-15才被抢修通行。利用2022-09-14拍摄的从沟口至冰川舌末端低海拔区，面积约43 km 2 的无人机航空影像，共解译同震地质灾害656处，总面积3.2 km 2 ，占影像覆盖区面积的7.4%。海螺沟一号营地附近同震滑坡分布如 图2 图2 所示。从 图2 图2 可以看出，流域内同震地质灾害最为严重的区域为一号营地上下游1.5 km范围，上下游各有一处大型滑坡将道路完全阻断。 海螺沟冰川卫星影像和照片如 图3 图3 所示。海螺沟冰川是典型的季风海洋型温性现代冰斗-山谷冰川，长约13.1 km，面积约25.7 km 2 ［ 10 10 - 12 12 10-12 ］ 。在全球气候变暖影响下，冰川消融严重，近100多年来海螺沟冰川末端海拔上升了约300 m，年均升高1.64 m ［ 11 11 ］ 。海螺沟冰川自上而下分为粒雪盆（4 980~7 556 m）、冰瀑布（3 700~4 980 m）和冰舌（2 900~3 700 m）3个部分（ 图3 图3 （a）），冰瀑布历史上曾与冰舌相连，但近年来由于冰川持续变薄，已与冰舌完全断开，形成垂直高差高约1 080 m的陡崖（ 图3 图3 （b）），上部粒雪盆地形相对宽缓，堆积的大量冰川从不同方向不断向下运动，行至冰瀑布上游狭窄沟谷段汇集收敛，地形坡度变陡且冰川前缘临空，导致该段流速显著升高，冰川被拉裂，裂隙发育，冰崩事件频发 ［ 11 11 - 12 12 11-12 ］ 。通过网络资料检索和历史Planet卫星影像对比，冰瀑布位置至少发生了2018-04-03、2018-08-18、2022-04-30、2022-07-29等多次大规模冰崩事件（ 图4 图4 ）。 2　光学卫星影像冰川面积监测 1 2.1　遥感数据 2 卫星影像是否适合冰川测绘取决于云量、采集日期、季节性降雪的存在和采集时的照明几何形状（太阳高度和方位角） ［ 13 13 - 14 14 13-14 ］ 。本文从美国地质调查局（ http：//earthexplorer.usgs.gov http://earthexplorer.usgs.gov ）下载2016—2022年Landsat 8和Sentinel-2卫星影像数据。为了确保冰川面积提取的准确性，选取研究区内无云且采集时间在8月—9月的遥感影像，数据信息见 表1 表1 。 年份 卫星 分辨率/m 采集日期 2016 Landsat 8 30 2016-08-25 2017 Landsat 8 30 2017-09-13 2018 Sentinel-2 10 2018-08-28 2019 Landsat 8 30 2019-08-25 2020 Sentinel-2 10 2020-08-27 2021 Sentinel-2 10 2021-08-02 2022 Sentinel-2 10 2022-08-12 2.2　冰川面积提取方法 2 NDSI是常用的植被指数方法在冰川提取中的扩展和延伸，能够很好区分冰川与背景信息，其计算式为： <math display="block" id="M1"><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">I</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">G</mi></mrow></msub><mo>-</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">W</mi><mi mathvariant="normal">I</mi><mi mathvariant="normal">R</mi></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">G</mi></mrow></msub><mo>+</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">W</mi><mi mathvariant="normal">I</mi><mi mathvariant="normal">R</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117261&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117259&type= 54.27133560 4.82600021 式中， <math id="M2"><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">G</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117251&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117263&type= 3.72533321 4.57200003 、 <math id="M3"><msub><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">W</mi><mi mathvariant="normal">I</mi><mi mathvariant="normal">R</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117270&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117268&type= 7.36600018 4.57200003 分别为可见光绿波段和短波红外（short-wave length infrared，SWIR）波段的反射率。 阈值分割方法是遥感冰川信息提取的有效方法，根据文献［ 15 15 ］确定海螺沟冰川的阈值为0.2，即NDSI > 0.2的像素被划分为冰川，NDSI≤0.2被划分为非冰川。 2.3　2016—2022年冰川面积变化 2 利用上述遥感数据和方法，得到2016—2022年海螺沟冰川8月的面积变化曲线，如 图5 图5 所示（2017年为9月数据）。从 图5 图5 可以看出，海螺沟冰川面积在6年间呈波动变化趋势，2016年为19.3 km 2 ，2017年增加至19.4 km 2 ，2018年剧减至18.2 km 2 ，为近6年内最小，2019年与2018年相当，为18.3 km 2 ，2020—2021年冰川面积呈显著增高趋势，2021年达到20.7 km 2 ，为近6年内最大，2022年又显著降低至18.9 km 2 。 由于目前尚未收集到海螺沟冰川附近气象站气象数据，本文从咕呱天气网站（ http：//www.ggtqw.com/2016luding/8yue.html http://www.ggtqw.com/2016luding/8yue.html ）下载了泸定县2016—2022年8月逐日最高和最低气温数据，用于应急状态下气温与冰川面积变化的分析，计算得到的每年8月平均最高气温和最低气温如 图5 图5 所示。从 图5 图5 可以看出，2016—2022年海螺沟冰川8月面积变化趋势与泸定县8月日平均最高气温具有负相关性，即平均最 高气温增高则冰川面积减少。如2019年8月平均最高气温由28.55 ℃逐渐降低至2021年的26.84 ℃，冰川面积则由2019年的18.3 km 2 增加至2021年的20.7 km 2 。2022年8月四川省多次出现有记录以来最高温，泸定县8月平均最高气温也上升至30.55 ℃，冰川面积也减少至18.9 km 2 。 3　雷达卫星影像冰川形变监测 1 3.1　雷达卫星数据 2 为了探测海螺沟冰川受地震影响情况，本文收集了欧洲空间局提供的地震前后6景Sentinel-1雷达卫星影像数据。Sentinel-1雷达卫星搭载了C波段合成孔径雷达，不同于遥感影像受天气影响的特征，雷达卫星可以全天时全天候提供连续影像 ［ 16 16 ］ 。本文选取通过干涉宽幅工作模式（interferometric wide swath， IWS）获取的空间分辨率为5 m×20 m的合成孔径雷达（synthetic aperture radar，SAR）影像，对海螺沟冰川进行数据处理，影像获取时间分别为2021-08-31、2021-09-12、2022-08-02、2022-08-14、2022-08-26、2022-09-07。 3.2　形变监测方法 2 由于合成孔径雷达差分干涉测量（differential interferometric SAR，DInSAR）对监测目标相干性要求较高，在类似冰川运动等大梯度地表形变场的监测中容易造成失相干和相位解缠错误 ［ 17 17 ］ ，因此本文选用POT对海螺沟冰川运动情况进行监测。POT技术主要通过对前后两景SAR影像的强度信息或相位信息精确配准后，获取同名像素坐标偏移量 ［ 18 18 ］ 。所获得的偏移量主要包括地表形变引起的偏移量、地形起伏引起的偏移量以及飞行轨道和成像姿态引起的偏移量 ［ 10 10 ］ ，计算式为： <math display="block" id="M4"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">t</mi></mrow></msub><mo>=</mo><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">f</mi></mrow></msub><mo>+</mo><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">b</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi></mrow></msub><mo>+</mo><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">M</mi><mtext> </mtext></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117288&type= 42.16400146 3.80999994 式中， <math id="M5"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">f</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117277&type= 7.53533268 3.72533321 表示配准后获取的偏移量； <math id="M6"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">f</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117310&type= 5.58799982 3.72533321 为地表形变引起的偏移量； <math id="M7"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">b</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117327&type= 6.85799980 3.72533321 为飞行轨道和成像姿态引起的偏移量； <math id="M8"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">M</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117308&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117306&type= 7.70466709 3.72533321 表示地形起伏引起的偏移量。 联合无形变区的轨道偏移数据以及外部数字高程模型（digital elevation model，DEM），采用最小二乘准则构建出系统偏移模型，去除地形起伏度（ <math id="M9"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">D</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">M</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117323&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117317&type= 7.70466709 3.72533321 ）、飞行轨道及成像姿态（ <math id="M10"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">b</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117331&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117327&type= 6.85799980 3.72533321 ）引起的偏移量，从而获得海螺沟冰川地表形变引起的偏移量（ <math id="M11"><msub><mrow><mi>D</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">f</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=45117310&type= 5.58799982 3.72533321 ）。 3.3　监测结果分析 2 利用雷达卫星数据和POT方法，得到了4期贡嘎山主峰及海螺沟周边冰川的结果，分别是2021-08-31—2021-09-12、2022-08-02—2022-08-14、2022-08-14—2022-08-26、2022-08-26—2022-09-07（以下分别简称Ⅰ、Ⅱ、Ⅲ、Ⅳ期），其中泸定地震前后（Ⅳ期）的POT结果如 图6 图6 所示。由 图6 图6 可以看出，此次地震导致贡嘎山主峰东侧海螺沟、磨子沟、燕子沟、南门关沟等冰川及其后缘山体均探测到较显著形变信号，而主峰西侧冰川受地震影响相对较小。其中海螺沟冰川冰瀑布区形变最为显著，形变显著区面积约 1.5 km 2 ，位移量最大约30 m，最大平均运动速度约2.5 m/d。 4个期次的结果图上均显示该区域有较显著的形变信号，本文将重点对该区域POT结果进行对比分析。首先，以2022-07-29 Planet卫星影像为底图，根据各期次形变渲染图上不同形变量的色彩差异（位移量 > 4.0 m）圈定形变显著区范围，结果如 图7 图7 所示。分析发现，4个期次的海螺沟冰川强形变区面积呈现线性增大的趋势，如 图8 图8 所示，由Ⅰ期的0.49 km 2 逐渐增加至Ⅳ期的1.51 km 2 。在空间上，强形变区依次向沟谷左岸和下游扩展。推测原因一方面在于2021年8月研究区平均最高气温（26.84 ℃）比2022年同期（30.55 ℃）降低了3.71 ℃，因此冰川运动性2021年比2022年低；另一方面，2022年8月研究区气温持续升高导致冰川的运动性逐渐增强，强形变区的面积逐渐增加。值得注意的是，从 图7 图7 可以看出，此次地震并未导致强变形区的面积增加速率增大。 进一步在强形变区选取了1-1′和2-2′两条剖面线，对比分析剖面线上不同期次的形变特 征，其中1-1′剖面位于Ⅳ期强形变区中轴线位置，2-2′剖面位于Ⅰ期强形变区中轴线位置（ 图7 图7 ），分别绘制其位移量曲线和平均运动速度曲线（ 图9 图9 ）。具体分析如下。 1）1-1′剖面 在此剖面上，Ⅰ、Ⅱ期位移不显著，剖面线从其强变形区的边缘穿过。Ⅲ、Ⅳ期在剖面线0~1 300 m段位移基本一致，运动速度呈缓慢增加趋势，运动速度由0.5 m/d增至1.0 m/d。Ⅲ期运动速度从1 300 m后快速增加，到1 880 m处运动速度增加至最大2.7 m/d；之后运动速度陡降，到2 200 m处运动速度降至0.1 m/d。Ⅳ期运动速度从1 500 m处显著增加，到1 810 m处运动速度由1.2 m/d增加至最大2.3 m/d；之后运动速度呈现剧烈波动下降的趋势。 2）2-2′剖面 在此剖面上，Ⅱ、Ⅲ、Ⅳ期在0~1 300 m段运动速度基本一致，均呈缓慢增加趋势；Ⅰ期也呈现缓慢增加趋势，但速率小于Ⅱ、Ⅲ、Ⅳ期。Ⅰ期运动速度在剖面线0~1 680 m段呈缓慢增加的趋势，由0.25 m/d增加至1.0 m/d；之后运动速度陡降，到1 820 m处降至0.1 m/d。Ⅱ期运动速度在剖面线0~1 850 m段呈缓慢增加的趋势，由0.32 m/d增加至1.7 m/d；之后运动速度陡降，到2 000 m处降至0.1 m/d。Ⅲ期运动速度在剖面线0~1 420 m段呈缓慢增加的趋势，由0.36 m/d增加至1.2 m/d；在1 420 m后运动速度快速增加，至1 620 m处增加至最大2.4 m/d；之后运动速度陡降，到2 000 m处降至0.1 m/d。Ⅳ期运动速度在剖面线0~1 300 m段呈缓慢增加的趋势，由0.25 m/d增加至1.1 m/d；在1 300~2 000 m段运动速度呈缓慢降低趋势，降低至0.7 m/d，其中 在1 700m附近运动速度有所波动；2 000 m之后运动速度陡增至2.3 m/d，在2 150 m之后开始剧降，到2 280 m处降至0.05 m/d。 3.4　监测结果讨论 2 1）监测结果可靠性分析 文献［ 10 10 ］利用日本ALOS雷达卫星2007—2018年共38景PALSAR-1/2影像，基于POT技术提取了海螺沟冰川形变速率，结果显示冰瀑布及上游区域多年平均速度为1.2 m/d，与本文探测到的结果吻合（ 图9 图9 （b）、 图9 图9 （d））；同时探测到冰瀑布位置为运动速度最快的部位，运动速度最大为2.5 m/d，与本文泸定地震前一期次在1-1′剖面探测到的最大速度2.7 m/d接近（ 图9 图9 （b））。因此，本文在地震应急状态下虽只利用Sentinel-1雷达计算了4个期次冰川位移数据，但计算结果总体上应该是可靠的。后续笔者团队将收集更多期次雷达影像对该冰川运动特征进行深入分析。 2）地形坡度对运动速度的影响分析 从Ⅲ、Ⅳ期的1-1′剖面平均运动速度曲线可以看出，冰川运动速度从上至下可分为速度缓增段（0~1 300 m）、速度陡增段（1 300~1 880 m）和速度陡降段（1 880~2 600 m）3个区段（见 图10 图10 ）。通过地震前后Planet卫星影像（分辨率为3.0 m）对比发现，速度缓增段冰川表面裂缝细长、排列紧密有序；速度陡增段冰川表面裂缝宽大、排列较凌乱；速度陡降段则正好对应冰崩频发区。 图11 图11 为1-1′剖面坡度曲线和速度分段图。从 图11 图11 可以看出，1 300 m和1 880 m速度分界点分别对应地形上的裂点，其中速度缓增段平均坡度仅11.6°，速度陡增段平均坡度增加至20.5°，速度陡降段平均坡度进一步增加至34.6°。由此可见，地形坡度与冰川运动速度有强相关性，总体上冰川运动速度随地形坡度增高而增大；当平均坡度达到20.5°时，运动速度显著增大；当平均坡度达到34.6°时，运动速度急剧降低。对照光学影像（ 图10 图10 ）发现，运动速度剧降区在影像上呈高亮度显示，与周围冰川对比明显，正好与冰瀑布前端的冰崩高发区对应。 运动速度急剧降低可能有两方面的原因：一方面，由于该区域地形坡度陡，前缘临空条件好，不同规模的冰崩频繁发生。对于冰崩频繁发生的区域，地表破坏严重，前后两期雷达影像强度图上同名像素点丢失，导致POT技术偏移估计信噪比较低 ［ 10 10 ］ ，计算结果可能失真，因此该区域测得的地表速度并不一定能真实反映冰川主体的运动速度。另一方面，冰瀑布三维遥感影像如 图12 图12 所示，呈上下游窄、中间宽的地貌特征，上游宽约850 m，中部最宽约1 100 m，下游宽约580 m，只有中部宽度的一半，有可能对冰川运动起到阻挡作用，使冰川运动速度迅速降低，后缘冰川不断在该处堆积增高，达到一定程度，发生冰崩失稳。 3）气温对运动速度的影响分析 如上文所述，冰川面积与日平均最高气温具 有负相关性，那么气温与冰川运动速度是否相关呢？如 图9 图9 （b）所示，在1-1′剖面速度缓增段，Ⅲ、Ⅳ期的平均运动速度相当；在速度陡增段，Ⅲ期的平均运动速度要显著大于Ⅳ期。如 图9 图9 （d）所示，在2-2′剖面速度缓增段，Ⅱ、Ⅲ、Ⅳ期的平均运动速度相当，Ⅰ期的平均运动速度略小于Ⅱ、Ⅲ、Ⅳ期；在速度陡增段，运动速度从高到低依次是Ⅲ、Ⅱ、Ⅰ、Ⅳ期。利用从咕呱天气网站下载的气温数据计算4个期次对应的日平均最高气温，分别为25.46 ℃、31.38 ℃、32.08 ℃、24.5 ℃。可见日平均最高气温从高到低依次也是Ⅲ、Ⅱ、Ⅰ、Ⅳ期，与2-2′剖面速度陡增段的速度大小关系正好对应，满足日平均最高气温越大，平均运动速度也越大的特点；在1-1′剖面上速度陡增段，如果排除Ⅰ、Ⅱ期（剖面线从强形变区边缘通过），则Ⅲ期平均运动速度大于Ⅳ期，也满足日平均最高气温越大、平均运动速度也越大的特点。由此可见，冰川平均运动速度与日平均最高气温具有正相关性。 4）泸定地震对冰川影响分析 如上文1-1′和2-2′剖面平均运动速度分析，在1-1′剖面上，Ⅳ期在1 810 m之前的平均运动曲线与震前相比无明显异常，1 810 m之后呈现剧烈波动下降的趋势，与震前连续下降趋势有明显差异；在2-2′剖面上，Ⅳ期在1 700 m之前的平均运动速度与震前相比无明显异常，1 700 m之后出现波动，并在2 000 m之后陡增，2 150 m后开始剧降，与震前连续下降趋势有明显差异。推测上述运动速度异常主要受此次地震扰动影响。对比地震前后3.0 m分辨率的Planet卫星影像，发现地震疑似在对应运动速度异常区域触发了冰崩。但由于现有震后Planet卫星影像分辨率有限且有云雾影响，尚无法判断发生冰崩的具体规模和位置，需要后续进行现场调查或获取震后高分辨率遥感影像进行对比判断。 通过对地震前后冰瀑布位置运动速度的对比分析，可以初步判断此次地震对冰川的直接影响区域主要是冰瀑布前缘区域，对后缘冰川直接影响较小，震后发生大规模冰崩的可能性较小，而且冰瀑布下方较长范围的宽缓冰舌区对冰崩物质具有较好的缓冲作用，初步判断震后冰崩事件直接致灾的可能性较小。但此次地震在流域内触发了大量崩塌、滑坡灾害（ 图2 图2 ），并使沟内冰碛物和周边山体震裂松动，为泥石流灾害的发生提供了丰富的物源，未来发生暴雨泥石流的风险增加，后续应加强相关调查和研究。 4　结　语 1 本文综合利用多时相光学卫星影像NDSI计算和雷达卫星影像POT技术对海螺沟冰川面积和运动速度进行定量监测和详细分析，得到以下初步认识： 1）基于Landsat 8和Sentinel-2光学卫星影像NDSI 计算结果，发现2016—2022年海螺沟冰川面积呈波动变化趋势，2018年面积最小，2021年面积最大，与泸定县日平均最高气温具有较好的负相关性。 2）基于泸定地震前后Sentinel-1雷达卫星影像POT计算结果，发现海螺沟冰川强形变区主要位于冰瀑布范围，面积约1.5 km 2 ，位移量最大约30 m，最大平均运动速度约2.5 m/d。强形变区从上往下可分为速度缓增段、速度陡增段和速度陡降段。不同段冰川表面具有明显的形态特征差异，速度缓增段冰川表面裂缝细长、排列紧密有序，速度陡增段冰川表面裂缝宽大、排列较凌乱，速度陡降段则正好对应冰崩频发区。 3）地形坡度与冰川运动速度有强相关性，总体上冰川运动速度随地形坡度增高而增大，当平均坡度达到20.5°时，运动速度显著增大，当平均坡度达到34.6°时，运动速度急剧降低。推测运动速度急剧降低的原因一方面可能为冰崩作用使得偏移估计信噪比低，计算结果失真；另一方面可能为冰瀑布前缘沟道显著变窄，对冰川运动起到阻挡作用。 4）冰川平均运动速度与日平均最高气温具有正相关性，2022-08-14—2022-08-26期间冰川运动速度最大，正好与该段时间川渝地区出现大范围极高温异常天气对应。 5）泸定地震未造成海螺沟冰川运动速度大范围显著增大，但对冰瀑布前缘区域扰动明显。震后发生大规模冰崩的可能性较小，而且冰瀑布下方较长范围的宽缓冰舌区对冰崩物质具有较好的缓冲作用，初步判断震后冰崩事件直接致灾的可能性较小。但震后沟内泥石流物源显著增多，未来发生暴雨泥石流的风险增加，后续应加强泥石流物源调查和致灾风险评价。 本文是地震应急状态下开展的初步研究，利用的Sentinel-1雷达影像数据期次较少，后续将收集更多期次雷达影像对该冰川运动特征进行长时间序列分析，进一步揭示该冰川运动规律和地震影响。 参 考 文 献 Wang Xin ， Fang Chengyong ， Tang Xiaochuan ， et al . Research on Emergency Evaluation of Landslides Induced by Luding Ms 6.8 Earthquake ［J］. Geomatics and Information Science of Wuhan University ， 2022 ， DOI： 10.13203/J.whugis20220586 http://dx.doi.org/10.13203/J.whugis20220586 . Research on Emergency Evaluation of Landslides Induced by Luding Ms 6.8 Earthquake 王欣 ， 方成勇 ， 唐小川 ， 等 . 泸定Ms 6.8级地震诱发滑坡应急评价研究 ［J］. 武汉大学学报（信息科学版） ， 2022 ， DOI： 10.13203/J.whugis20220586 http://dx.doi.org/10.13203/J.whugis20220586 . 泸定Ms 6.8级地震诱发滑坡应急评价研究 Fan Xuanmei ， Wang Xin ， Dai Lanxin ， et al . Characteristics and Spatial Distribution Pattern of Ms 6.8 Luding Earthquake Occurred on September 5， 2022 ［J］. Journal of Engineering Geology ， 2022 . http：//www.gcdz.org/article/id/77d603ac-db72-4aac-822e-f37abe547b9c http://www.gcdz.org/article/id/77d603ac-db72-4aac-822e-f37abe547b9c . Characteristics and Spatial Distribution Pattern of Ms 6.8 Luding Earthquake Occurred on September 5， 2022 范宣梅 ， 王欣 ， 戴岚欣 ， 等 . 2022年Ms 6.8级泸定地震诱发地质灾害特征与空间分布规律研究 ［J］. 工程地质学报 ， 2022 . http：//www.gcdz.org/article/id/77d603ac-db72-4aac-822e-f37abe547b9c http://www.gcdz.org/article/id/77d603ac-db72-4aac-822e-f37abe547b9c . 2022年Ms 6.8级泸定地震诱发地质灾害特征与空间分布规律研究 Cheng Jia ， Xu Xiwei . Features of Earthquake Clustering from Calculation of Coulomb Stress Around the Bayan Har Block， Tibetan Plateau ［J］. Seismology and Geology ， 2018 ， 40 （ 1 ）： 133 - 154 . Features of Earthquake Clustering from Calculation of Coulomb Stress Around the Bayan Har Block， Tibetan Plateau 程佳 ， 徐锡伟 . 巴颜喀拉块体周缘强震间应力作用与丛集活动特征初步分析 ［J］. 地震地质 ， 2018 ， 40 （ 1 ）： 133 - 154 . 巴颜喀拉块体周缘强震间应力作用与丛集活动特征初步分析 Bai Mingkun ， MarieLuce C ， Li Haibing ， et al . Late Quaternary Slip Rate and Earthquake Hazard Along the Qianning Segment， Xianshuihe Fault ［J］. Acta Geologica Sinica ， 2022 ， 96 （ 7 ）： 2312 - 2332 . Late Quaternary Slip Rate and Earthquake Hazard Along the Qianning Segment， Xianshuihe Fault 白明坤 ， Chevalier Marie-Luce ， 李海兵 ， 等 . 鲜水河断裂带乾宁段晚第四纪走滑速率及区域强震危险性研究 ［J］. 地质学报 ， 2022 ， 96 （ 7 ）： 2312 - 2332 . 鲜水河断裂带乾宁段晚第四纪走滑速率及区域强震危险性研究 Xu Xiwei ， Wu Xiyan ， Yu Guihua ， et al . Seismo-Geological Signatures for Identifying M≥7.0 Earthquake Risk Areas and Their Premilimary Application in China's Mainland ［J］. Seismology and Geology ， 2017 ， 39 （ 2 ）： 219 - 275 . Seismo-Geological Signatures for Identifying M≥7.0 Earthquake Risk Areas and Their Premilimary Application in China's Mainland 徐锡伟 ， 吴熙彦 ， 于贵华 ， 等 . 中国大陆高震级地震危险区判定的地震地质学标志及其应用 ［J］. 地震地质 ， 2017 ， 39 （ 2 ）： 219 - 275 . 中国大陆高震级地震危险区判定的地震地质学标志及其应用 Deng Qidong ， Cheng Shaoping ， Ma Ji ， et al . Seismic Activities and Earthquake Potential in the Tibetan Plateau ［J］. Chinese Journal of Geophysics ， 2014 ， 57 （ 7 ）： 2025 - 2042 . Seismic Activities and Earthquake Potential in the Tibetan Plateau 邓起东 ， 程绍平 ， 马冀 ， 等 . 青藏高原地震活动特征及当前地震活动形势 ［J］. 地球物理学报 ， 2014 ， 57 （ 7 ）： 2025 - 2042 . 青藏高原地震活动特征及当前地震活动形势 Wang Minjie ， Li Tianbin ， Meng Lubo ， et al . Back Analysis of Stress Field in the Intersection Region of Y Shaped Fault， Sichuan ［J］. Journal of Railway Science and Engineering ， 2015 ， 12 （ 5 ）： 1088 - 1095 . Back Analysis of Stress Field in the Intersection Region of Y Shaped Fault， Sichuan 王敏杰 ， 李天斌 ， 孟陆波 ， 等 . 四川“Y字形”断裂交汇部应力场反演分析 ［J］. 铁道科学与工程学报 ， 2015 ， 12 （ 5 ）： 1088 - 1095 . 四川“Y字形”断裂交汇部应力场反演分析 Xu Jing ， Shao Zhigang ， Ma Hongsheng ， et al . Evolution of Coulomb Stress and Stress Interaction Among Strong Earthquakes Along the Xianshuihe Fault Zone ［J］. Chinese Journal of Geophysics ， 2013 ， 56 （ 4 ）： 1146 - 1158 . Evolution of Coulomb Stress and Stress Interaction Among Strong Earthquakes Along the Xianshuihe Fault Zone 徐晶 ， 邵志刚 ， 马宏生 ， 等 . 鲜水河断裂带库仑应力演化与强震间关系 ［J］. 地球物理学报 ， 2013 ， 56 （ 4 ）： 1146 - 1158 . 鲜水河断裂带库仑应力演化与强震间关系 Xiong Wei ， Tan Kai ， Yu Pengfei ， et al . Triggering of Mw 5.9 Kangding Earthquake by Coulomb Stress Evolution Along Xianshuihe Fault Zone Since 1955 ［J］. Journal of Geodesy and Geodynamics ， 2016 ， 36 （ 2 ）： 95 - 100 . Triggering of Mw 5.9 Kangding Earthquake by Coulomb Stress Evolution Along Xianshuihe Fault Zone Since 1955 熊维 ， 谭凯 ， 余鹏飞 ， 等 . 鲜水河断裂近期库仑应力演化及其与康定Mw 5.9地震的关系 ［J］. 大地测量与地球动力学 ， 2016 ， 36 （ 2 ）： 95 - 100 . 鲜水河断裂近期库仑应力演化及其与康定Mw 5.9地震的关系 Liu Guoxiang ， Zhang Bo ， Zhang Rui ， et al . Monitoring Dynamics of Hailuogou Glacier and the Secondary Landslide Disasters Based on Combination of Satellite SAR and Ground-Based SAR ［J］. Geomatics and Information Science of Wuhan University ， 2019 ， 44 （ 7 ）： 980 - 995 . Monitoring Dynamics of Hailuogou Glacier and the Secondary Landslide Disasters Based on Combination of Satellite SAR and Ground-Based SAR 刘国祥 ， 张波 ， 张瑞 ， 等 . 联合卫星SAR和地基SAR的海螺沟冰川动态变化及次生滑坡灾害监测 ［J］. 武汉大学学报（信息科学版） ， 2019 ， 44 （ 7 ）： 980 - 995 . 联合卫星SAR和地基SAR的海螺沟冰川动态变化及次生滑坡灾害监测 Liu Qiao ， Liu Shiyin ， Zhang Yong ， et al . Surface Ablation Features and Recent Variation of the Lower Ablation Area of the Hailuogou Glacier， Mt.Gongga ［J］. Journal of Glaciology and Geocryology ， 2011 ， 33 （ 2 ）： 227 - 236 . Surface Ablation Features and Recent Variation of the Lower Ablation Area of the Hailuogou Glacier， Mt.Gongga 刘巧 ， 刘时银 ， 张勇 ， 等 . 贡嘎山海螺沟冰川消融区表面消融特征及其近期变化 ［J］. 冰川冻土 ， 2011 ， 33 （ 2 ）： 227 - 236 . 贡嘎山海螺沟冰川消融区表面消融特征及其近期变化 Liu Qiao ， Zhang Yong . Studies on the Dynamics of Monsoonal Temperate Glaciers in Mt.Gongga： A Review ［J］. Mountain Research ， 2017 ， 35 （ 5 ）： 717 - 726 . Studies on the Dynamics of Monsoonal Temperate Glaciers in Mt.Gongga： A Review 刘巧 ， 张勇 . 贡嘎山海洋型冰川监测与研究： 历史、现状与展望 ［J］. 山地学报 ， 2017 ， 35 （ 5 ）： 717 - 726 . 贡嘎山海洋型冰川监测与研究： 历史、现状与展望 Wang L C ， Yu K ， Chang L ， et al . Response of Glacier Area Variation to Climate Change in the Kaidu-Kongque River Basin， Southern Tianshan Mountains During the last 20 Years ［J］. China Geology ， 2021 ， 4 （ 3 ）： 389 - 401 . Response of Glacier Area Variation to Climate Change in the Kaidu-Kongque River Basin， Southern Tianshan Mountains During the last 20 Years Wu Renzhe ， Huang Xinru ， Li Yong ， et al . Analysis of the Snow Line Change and its Influencing Factors in the Glaciers in Bugyai Kangri Region from 1994 to 2016 ［J］. Journal of Geomatics ， 2020 ， DOI： 10.14188/j.2095-6045.2020075 http://dx.doi.org/10.14188/j.2095-6045.2020075 . Analysis of the Snow Line Change and its Influencing Factors in the Glaciers in Bugyai Kangri Region from 1994 to 2016 吴仁哲 ， 黄心茹 ， 李勇 ， 等 . 布加岗日地区冰川1994至2016年雪线变化及其影响因素分析 ［J］. 测绘地理信息 ， 2020 ， DOI： 10.14188/j.2095-6045.2020075 http://dx.doi.org/10.14188/j.2095-6045.2020075 . 布加岗日地区冰川1994至2016年雪线变化及其影响因素分析 Pratibha S ， Kulkarni A V . Decadal Change in Supraglacial Debris Cover in Baspa Basin， Western Himalaya ［J］. Current Science ， 2018 ， 114 （ 4 ）： 792 . Decadal Change in Supraglacial Debris Cover in Baspa Basin， Western Himalaya Lu Huiyan ， Li Weile ， Xu Qiang ， et al . Early Detection of Landslides in the Upstream and Downstream Areas of the Baige Landslide， the Jinsha River Based on Optical Remote Sensing and InSAR Technologies ［J］. Geomatics and Information Science of Wuhan University ， 2019 ， 44 （ 9 ）： 1342 - 1354 . Early Detection of Landslides in the Upstream and Downstream Areas of the Baige Landslide， the Jinsha River Based on Optical Remote Sensing and InSAR Technologies 陆会燕 ， 李为乐 ， 许强 ， 等 . 光学遥感与InSAR结合的金沙江白格滑坡上下游滑坡隐患早期识别 ［J］. 武汉大学学报（信息科学版） ， 2019 ， 44 （ 9 ）： 1342 - 1354 . 光学遥感与InSAR结合的金沙江白格滑坡上下游滑坡隐患早期识别 Deng Fanghui ， Zhou Chunxia ， Wang Zemin ， et al . Ice-Flow Velocity Derivation of the Confluence Zone of the Amery Ice Shelf Using Offset-Tracking Method ［J］. Geomatics and Information Science of Wuhan University ， 2015 ， 40 （ 7 ）： 901 - 906 . Ice-Flow Velocity Derivation of the Confluence Zone of the Amery Ice Shelf Using Offset-Tracking Method 邓方慧 ， 周春霞 ， 王泽民 ， 等 . 利用偏移量跟踪测定Amery冰架冰流汇合区的冰流速 ［J］. 武汉大学学报（信息科学版） ， 2015 ， 40 （ 7 ）： 901 - 906 . 利用偏移量跟踪测定Amery冰架冰流汇合区的冰流速 Wasowski J ， Bovenga F . Investigating Landslides and Unstable Slopes with Satellite Multi Temporal Interferometry： Current Issues and Future Perspectives ［J］. Engineering Geology ， 2014 ， 174 ： 103 - 138 . Investigating Landslides and Unstable Slopes with Satellite Multi Temporal Interferometry： Current Issues and Future Perspectives 感谢四川测绘地理信息局、四川省国防科技情报研究所和深圳飞马机器人科技有限公司提供震后遥感数据，同时感谢欧洲空间局提供的Sentinel-1雷达卫星数据和美国Planet卫星公司提供的Planet光学卫星数据。本文撰写过程中参考了中国科学院、水利部成都山地灾害与环境研究所刘巧研究员震后发表的关于地震对海螺沟冰川影响的相关科普论文，在此表示感谢。 10.13203/j.whugis20220593 P237 A 1671-8860（2023）01-0047-11 李为乐,陈俊伊,陆会燕，等.泸定Ms 6.8地震对海螺沟冰川的影响应急分析[J].武汉大学学报（信息科学版）,2023,48(1):47-57.DOI: 10.13203/j.whugis20220593 http://dx.doi.org/10.13203/j.whugis20220593 李为乐,陈俊伊,陆会燕，等.泸定Ms 6.8地震对海螺沟冰川的影响应急分析[J].武汉大学学报（信息科学版）,2023,48(1):47-57.DOI:10.13203/j.whugis20220593 LI Weile,CHEN Junyi,LU Huiyan,et al.Emergency Analysis of the Impact of the Luding Ms 6.8 Earthquake on Hailuogou Glacier[J].Geomatics and Information Science of Wuhan University,2023,48(1):47-57.DOI: 10.13203/j.whugis20220593 http://dx.doi.org/10.13203/j.whugis20220593 LI Weile,CHEN Junyi,LU Huiyan,et al.Emergency Analysis of the Impact of the Luding Ms 6.8 Earthquake on Hailuogou Glacier[J].Geomatics and Information Science of Wuhan University,2023,48(1):47-57.DOI:10.13203/j.whugis20220593 国家重点研发计划（2021YFC3000401）；国家自然科学基金（41941019）；第二次青藏高原科学考察项目（2019QZKK0201）。 2022-09-18 2023-01-05 2023-06-29 33226352 武汉大学学报（信息科学版） 1 48