论文部分内容阅读
受全球气候变暖影响,近年来多数冰川物质平衡处于持续亏损状态。而我国的山地冰川储备丰富,且大多集中在藏东南地区。这些冰川不仅是我国淡水资源的主要源头之一,对于藏东南地区甚至整个西南地区的气候有着重要影响,而且冰川运动可能引发一系列的自然灾害。因此,研究山地冰川的运动状态及影响其动态变化的相关因素,可为冰川退化、气候分析以及冰川灾害预警等提供理论和数据支撑。为此,本文主要做了以下工作:首先,本文系统的分析了目前冰川运动监测常用的方法,根据各方法的优缺点及山地冰川的特点,最终本文选择像素偏移跟踪技术(Pixel Offset Tracking,POT)提取冰川流速场。然而,冰川表面极易发生匹配失相干,影像对之间的相干性过低会导致失匹配像元数量增加,从而引起监测空间覆盖度的下降,同时误匹配像元的数量也会增加,进而影响POT解算结果的精度及可靠性。为了克服该问题,本文提出了PO-SBAS模型,通过限制影像对的时间基线长度来保持影像对之间的相关性,降低匹配失相干的影响。另外,通过影像对组合获取了较多的冗余观测量,解算时可降低误匹配噪声干扰,从而提高解算精度。其次,考虑到冰面地物散射特性与季节变化的关系,而地面散射特性与匹配失相干息息相关。因此,本文在PO-SBAS模型的基础上,提出一种分季节的时间基线优化组合方法,使PO-SBAS模型更适用于冰川运动监测。基于所提出的方法,对藏东地区布加岗日冰川进行分析。获取2015~2019年间覆盖研究区域的76景C波段Sentinel-1卫星的降轨SAR影像,首先,利用PO-SBAS模型提取研究区域的冰川流速场,并以冰川区域内失匹配像元的数量为评价指标,分析不同季节以及相同季节不同时间基线长度对匹配失相干影响。结果表明,消融前期冰川表面地物散射特性变化剧烈,随时间基线长度增加匹配质量下降,消融后期48天及以上的基线组合匹配质量迅速下降,冻结期内对96天及以内的基线长度变化不敏感,跨越冻结期和消融期的影像组合随时间基线长度增加匹配质量下降。然后,基于上述分析,提出分季节的时间基线优化组合方法,即将每年冻结期、消融前期和消融后期的SAR影像分别构成子基线集。其中,冻结期内,基线长度控制在96天及以内;消融前期以及跨越冻结期和消融期内,则选取最短时间基线;而消融后期内,基线长度控制在48天以内。以这种基线组合方式优化PO-SBAS模型,解算时序冰川流速场。第三,基于基线组合优化后的PO-SBAS时序解算模型,本文首次提取了布加岗日冰川运动场,并分析了运动场的时空演变过程与冰川规模、地形、海拔、气候等因素的联系。研究表明,(1)布加岗日冰川平均流速局部地区可达到0.94 m/d,从冰川运动场与冰川空间分布的角度看,冰川运动速度与冰川面积呈正相关关系,冰川面积越大则运动速度越快,积累区的运动速度高于消融区,越接近冰舌末端运动速度越慢,经过陡坡区域运动速度加快。冰舌末端是冰川对于气候变化反应最敏感的区域,冰川流速的时序波动最剧烈,随着高程上升冰川流速的时序差异逐渐降低。受气候影响,该地区迎风西南坡冰川的运动速度快于背风东北坡,速度波动也更剧烈,对于气候变化反应更加敏感。(2)从布加岗日冰川运动场的时序变化角度看,冰川流速呈季节性波动,春夏季节的运动速度大于秋冬季节。冰舌末端流速随着温度升高而升高,对气候变化的响应具有同步性,随着高程上升冰川流速对气候变化的响应时间逐渐滞后,季节性变化幅度也逐渐减小。(3)从布加岗日冰川表面高程变化情况来看,冰川厚度在2000~2014年间整体下降了1.07 m,冰川物质处于亏损状态,冰舌末端存在不同程度的厚度减薄现象,随着海拔上升冰川表面高程逐渐稳定。其中,南坡足学会(S1)冰川的厚度减薄程度最大,可达到-76.7 m,北坡的东碑加(N4)冰川呈厚度减薄的区域海拔位置更高,中轴线上的平均厚度变化为-15.70 m。结合布加岗日冰川运动场和高程变化情况分析,当冰川物质长期处于亏损状态,冰舌厚度逐渐减薄,冰川流速减缓。根据局部区域水热条件分布情况的不同,水汽来源更加丰富的南坡冰川相比北坡冰川的退缩量相对较小,冰川流速也更快。