《Remote Sensing Applications: Society and Environment》：Glacier dynamics monitoring of the Southern Patagonia Icefield using optical satellite images
编辑推荐：
冰川是全球气候变化的关键指示器。南巴塔哥尼亚冰原（Southern Patagonian Icefield, SPI）是世界第三大陆地冰体，估计总面积约13,000 km2，是南美洲最大的淡水储备库。鉴于研究区域的广阔性和难以到达性，卫星遥感被确定为通过欧洲空
冰川是全球气候变化的关键指示器。南巴塔哥尼亚冰原（Southern Patagonian Icefield, SPI）是世界第三大陆地冰体，估计总面积约13,000 km2，是南美洲最大的淡水储备库。鉴于研究区域的广阔性和难以到达性，卫星遥感被确定为通过欧洲空间局（ESA）的Sentinel-2星座分析冰川表面位移和速度的最合适方法。研究结果显示，该区域存在显著的冰体移动，速度值最高可达2 m/天，尤其是在消融区，且在东-西和北-南两个方向均有分布。研究人员对别德马（Viedma）冰川和廷德尔（Tyndall）冰川在2016年至2024年期间进行了全面的季节性趋势分析。结果显示位移速度有所增加，最高达2.6 m/天，凸显了由于温度变化引起的速度季节性趋势，这一趋势也得到了ESA的ERA5再分析气候数据集的印证。监测冰川动力学对于理解气候变化对环境的影响至关重要，因为它提供了关于这些冰体加速和运动的见解。因此，该方法为冰川位移提供了宝贵的信息，使得在难以进入和偏远地区进行数据收集成为可能。
**一、研究背景与科学问题**
冰川是全球重要的自然资源，也是全球和区域尺度上气候变化的关键指示器。温度波动和降雪模式在不同时间尺度上的变化——包括季节性、年际和年代际——都会对这些冰体产生影响。冰川表面速度是研究和分析冰川演化与变化以及缓解相关风险最重要的参数之一。然而，由于缺乏时间连续的流速数据，预测潜在灾害以及确定冰川变薄、加速和退缩之间的关系面临巨大挑战。这些信息对于计算水文平衡、物质平衡、侵蚀和冰川稳定性等关键参数至关重要。在全球和详细尺度上评估冰川运动，由于后勤挑战、可达性问题以及往往是先于科学挑战的技术和实际困难，并不容易实现。包括温度和降水模式、年降雪量和冰川速度在内的原位数据收集，由于研究区域难以进入以及数据有限或根本不存在而面临重重困难。
南巴塔哥尼亚冰原（SPI）位于阿根廷和智利边境，是继南极和格陵兰之后的世界第三大冰体，覆盖超过13,000平方公里，由48条主要冰川组成。由于其规模、恶劣天气和偏远位置，遥感技术被广泛用于研究该区域，使研究人员能够追踪雪冰分布、冰厚度以及质量和海拔的长期变化。卫星遥感为非常大且相对难以到达的区域进行定期观测提供了有效工具。本研究利用欧空局Sentinel-2星座获取的2021年光学影像数据集，使用COSI-Corr工具分析冰川的位移和速度，并将这些数据与包括温度在内的气候记录（ESA的ERA5卫星数据）进行比较。该研究的目的是利用光学影像估算冰川运动的日平均速度，并将该结果与外部参数相关联，以全面表征SPI的动力学特征，特别是在当前气候变化情景下。
**二、研究内容、技术方法与主要结论**
研究人员开展了以下工作：首先，选取2021年2月至4月（南半球夏末至初秋）期间的8对Sentinel-2光学影像，对整个SPI进行了大范围的冰川位移速度制图分析；其次，选取别德马（Viedma）冰川和廷德尔（Tyndall）冰川作为局部尺度研究对象，分析了2016年至2024年间夏冬两季的季节性形变速度；最后，将COSI-Corr获取的速度结果与Glims和ITS等已发表的参考数据集进行对比验证，并结合ERA5温度数据分析速度与温度之间的关系。研究表明，SPI整体最大形变速度可达±5 m/天，主要消融区速度超过2-3 m/天；Viedma冰川速度最高达2.62 m/天且无明显季节性差异，Tyndall冰川则表现出中等季节性，夏季速度约1.56 m/天高于冬季。与参考数据集对比显示两者高度一致（R2介于0.58至0.91之间），验证了光学影像相关方法在SPI区域冰川表面速度估算中的可靠性。然而，统计分析显示冰川速度与温度之间的相关性较低且不显著（p=0.38），表明短期冰川运动受内部冰动力学、基底条件和冰下水文等复杂因素的综合控制，而非温度单独主导。
**三、主要技术方法**
本研究采用的主要技术方法可概括为以下三项：第一，COSI-Corr光学影像相关模块，该模块集成于ENVI软件中，通过正交校正、配准和相关运算对光学卫星影像进行精确几何处理，采用频率相关器提取冰川位移的东-西（E/W）和北-南（N/S）两个分量；第二，Sentinel-2光学影像数据集，来自哥白尼计划的双卫星星座（Sentinel-2A于2015年发射，Sentinel-2B于2017年发射），其空间分辨率为10 m，本研究选用云量低于10%的影像以确保相关性质量；第三，ERA5再分析气候数据集，来自欧空局哥白尼项目，用于分析1950年至2020年研究区域的温度变化趋势及季节性温度异常。研究中还使用了信噪比（Signal-to-Noise Ratio, SNR）进行质量控制，过滤低相干性像元，并通过与Glims和ITS（基于Sentinel-1合成孔径雷达数据）数据集的对比进行交叉验证。
**四、研究结果**
**（1）南巴塔哥尼亚冰原——大范围分析**
对8对Sentinel-2影像进行COSI-Corr分析，研究人员生成了整个SPI区域的日形变速度图，包含东-西（E/W）和北-南（N/S）两个位移分量。最大形变速度记录在±5 m/天范围内，速度值超过2-3 m/天的区域主要集中在流域中部的消融区和以冰崩解为特征的冰川前端（calving glaciers）。沿冰川流向的形变模量图综合显示了受形变影响最严重的区域。分析主要冰川速度发现，朝向西侧的冰川（如Europa 4.48 m/天、Penguin 5.21 m/天、HPS-13 5.34 m/天、Pio XI 2.95 m/天）普遍比朝向东侧的冰川（如Upsala 2.46 m/天、Perito Moreno 3.4 m/天、Mayo 2.6 m/天）表现出更高的速度，但Tempeno冰川（1.91 m/天）和O'Higgins冰川（6.1 m/天）构成例外。
**（2）别德马（Viedma）冰川和廷德尔（Tyndall）冰川——局部尺度分析**
研究人员选取位于SPI中东部和南部的Viedma冰川和Tyndall冰川进行了2016至2024年间夏冬两季的详细分析，各自生成了18张日形变速度模量图。Viedma冰川在2016至2024年间夏季与冬季速度保持相对恒定，无明显季节性波动，最高速度出现在冰川中段和靠近湖泊的末端区域，速度值最高达2.62 m/天。Tyndall冰川则在大多数分析年份中记录到夏季位移率明显高于冬季，表明存在强烈的季节性影响，这与暖月较高温度相关，速度约为1.56 m/天，但2019年和2022年出现了夏冬位移率相对接近的例外情况。
**五、讨论总结与研究结论**
在讨论部分，研究人员将COSI-Corr获取的冰川表面速度与Glims和ITS两个参考数据集进行了对比验证。对SPI主要出口冰川的平均表面速度比较显示出强烈的一致性，大多数冰川的差异仅为十分之几米/天，在光学特征跟踪方法对中等分辨率卫星影像的典型不确定度范围内。通过100×100 m规则网格进行的空间相关性分析显示，散点图呈明显正相关，决定系数（R2）范围约为0.58至0.91，快速流动冰川（如Pio XI，R2=0.91）呈现较高相关性，而动动力学更复杂的冰川（如Tyndall，R2=0.58-0.60）相关性略低。不确定性评估考虑Sentinel-2影像10 m空间分辨率，假设位移不确定性约为1.5 m，估算速度不确定性范围约为0.02至0.06 m/天，显著低于实测冰川速度（0.4至1.2 m/天），表明导出的位移场可靠地代表了冰川运动。温度趋势分析显示，尽管存在显著季节性温度差异（Tyndall冬季温度约-5至-7°C，Viedma约-7至-12°C，两者夏季均稳定在约0°C），但2016-2024年间速度测量与相应温度值之间的统计比较显示低且不显著的相关性（p=0.38），表明温度变化并未直接影响冰川速度的短期变化，内部冰动力学、基底条件或冰下水文等其他因素可能在控制冰川运动中起主导作用。
研究结论指出：冰川代表重要的全球自然资源，是全球和区域气候变化的重要指标。本研究通过Sentinel-2光学影像结合COSI-Corr模块对SPI进行了分析，最高形变速度在±5 m/天范围内，最大速度超过2-3 m/天，主要在消融区和崩解冰川前端观测到。Viedma冰川在2016-2024年间无显著季节性差异，最高模量速度达2.62 m/天；Tyndall冰川则表现出中等季节性，冬季速度低于夏季，模量速度约1.56 m/天。包含速度和形变值的数据集以及显示该区域过去70年温度上升的图表，为冰川融化和气候变化影响提供了有力证据。随着全球温度持续上升，可以合理假设冰川速度和融化速率将会增加，冰川因其对温度变化的快速响应而作为气候变化的敏感指标。因此，持续监测冰川动力学对于深入理解气候变化的更广泛环境影响十分必要。
