1990-2019年贡嘎山地区典型冰川表碛覆盖变化及其空间差异

表碛覆盖型冰川是中国西部较为常见的冰川类型。表碛层存在于大气—冰川冰界面,强烈影响大气圈与冰冻圈之间的热交换。表碛厚度的空间异质性可极大地改变冰川的消融率和物质平衡过程,进而影响冰川径流过程和下游水资源。基于Landsat TM/TIRS数据,运用能量平衡方程反演了贡嘎山地区冰川表碛厚度,研究了贡嘎山地区冰川在1990—2019年间表碛覆盖范围及厚度变化情况,同时对比了东西坡差异。结果表明：① 贡嘎山地区冰川表碛扩张总面积达43.824 km²。其中,海螺沟冰川扩张2.606 km²、磨子沟冰川1.959 km²、燕子沟冰川1.243 km²、大贡巴冰川0.896 km²、小贡巴冰川0.509 km²、南门关沟冰川2.264 km²,年均扩张率分别为3.2%、11.1%、1.5%、0.9%、1.0%和6.5%;② 海螺沟冰川、磨子沟冰川、燕子沟冰川、大贡巴冰川、小贡巴冰川、南门关沟冰川表碛平均增厚分别为5.2 cm、3.1 cm、3.7 cm、6.8 cm、7.3 cm和13.1 cm;③ 西坡冰川表碛覆盖度高,表碛覆盖年均扩张率低,冰川末端退缩量小;东坡冰川表碛覆盖年均扩张率高,但表碛覆盖度总体低于西坡,冰川末端退缩量大。     

表碛覆盖冰川信息提取方法——以波密县为例

针对表碛覆盖冰川提取困难的问题,本文结合TM影像、DEM和野外调查资料,基于地物光谱、地形和纹理特征,将研究区分为表碛覆盖冰川区、无表碛覆盖冰川区和阴影下冰川区,采用模糊逻辑方法建立定量化提取规则,综合利用了面向对象和波段比值两个方法的各自优点,建立面向对象-波段比值信息提取方法,进行了有无表碛物覆盖的冰川信息的分区提取,并以波密县为典型研究区,通过混淆矩阵进行精度分析。研究结果表明:1)提出的面向对象—波段比值法,能够较完整地提取研究区表碛覆盖冰川,精度可达93%,冰川信息提取总体精度为98%,kappa系数为0.97;2)研究区冰川总面积为7 655.27 km~2,占研究区总面积的46.18%,无表碛覆盖冰川、表碛覆盖冰川、阴影下冰川面积分别为5615.20 km~2、65.13 km~2、1 974.94 km~2,占总冰川73.35%、0.85%、25.80%。初步建立的表碛覆盖冰川自动信息提取方法可以为类似地区冰川遥感监测提供依据和参考。     

2000—2014年喀喇昆仑山音苏盖提冰川表面高程变化

喀喇昆仑山区冰川由于存在正物质平衡或跃动、前进现象,被称之为“喀喇昆仑异常”,不过该地区冰川变化差异显著,尤其是大型表碛覆盖冰川,呈现与其他类型冰川明显的差异性响应,为理解喀喇昆仑冰川异常的机理,冰川尺度的详细变化研究十分必要。音苏盖提冰川位于喀喇昆仑山乔戈里峰北坡,是中国面积最大的冰川,是典型的大型表碛覆盖冰川。通过应用TanDEM-X/TerraSAR-X（2014年2月）与SRTM-X DEM（2000年2月）的差分干涉测量方法计算音苏盖提冰川表面高程变化,并结合冰川表面流速对冰川表面高程变化和跃动进行分析和讨论。结果表明：2000—2014年音苏盖提冰川表面高程平均下降了1.68±0.94 m,即冰川整体厚度在减薄,年变化率为-0.12±0.07 m·a^-1。冰川表面高程变化分布不均,其中南分支（S）冰流冰川整体减薄较为显著,冰川南分支冰流运动速度较快,前进/跃动的末端占据了冰川的主干,阻滞原主干冰川物质的向下运移（跃动）,导致原主干冰舌表面高程上升;冰川厚度减薄随着海拔升高先下降后保持稳定,同时呈现一定的波动性;低海拔表碛区域消融大于裸冰区,可能存在较薄表碛,因热传导高、覆盖大量冰面湖塘和冰崖存在,加速了冰川消融;在坡度小于30 °的区域,冰川厚度减薄随着坡度的减小而加剧;坡向朝南冰川厚度略微增加（0.01 m）,西南坡向冰川厚度略微减薄（-0.03 m）,其他坡向冰川厚度减薄明显。近14 a来,表碛覆盖的音苏盖提冰川表面高程整体下降表明物质处于亏损状态,冰川跃动导致局部冰川表面高程的增加。     

1990—2015年喜马拉雅山冰川变化的遥感监测

基于Landsat系列遥感数据,运用比值阈值法（B3/B5）和目视解译,研究1990—2015年喜马拉雅山冰川面积的分布与变化特征。结果表明：25年间研究区冰川面积共减少2553.10 km²,年均退缩率为0.44%/a,研究时段冰川加速退缩。研究区冰川主要分布在西段地区,中段次之,东段最少,近25年来西段、东段和中段地区冰川均表现为退缩趋势,其中东段地区退缩最快,中段最慢。从地形分布和变化特征看,5°~25°范围内冰川的分布面积较多,近25年来各坡度等级冰川均在退缩,其中25°~30°之间冰川面积退缩较快,在极平缓/极陡峭地区退缩较慢。尽管8个坡向上冰川均表现为退缩趋势,但退缩幅度有所差异,北坡与西北坡冰川退缩较慢,其他坡向退缩较快。研究时段表碛物覆盖型与非表碛物覆盖型冰川均在退缩,但后者的退缩幅度较大,表明研究区表碛物在一定程度上抑制了冰川消融。     

基于无人机摄影测量的梅里雪山明永冰川末端表面高程动态监测

位于横断山区的梅里雪山地区,受地形与气候的影响,该地区光学遥感及微波遥感应用较为受限。论文基于大疆精灵4RTK无人机搭配中海达UBase,于2018年11月9日和2019年11月12日对梅里雪山地区明永冰川末端开展了2期摄影测量,估算冰川消融区表面高程变化,监测冰川动态变化。研究结果表明:冰川区表面形态在不同位置有显著差别,航测冰川区中下部大量表碛覆盖,有少量沿冰川主流线方向发育的裂隙;航测冰川区中上部表碛覆盖较少,大量横向裂隙发育。2018年11月至2019年11月,明永冰川末端表面高程平均变化-1.67 m。冰川表面高程变化空间分布差异显著,在航测冰川区上部高程差存在正负交替现象,航测冰川区中部以减薄为主,而航测冰川区下部以增厚为主。与高亚洲地区其他山系冰川变化对比,梅里雪山地区是冰川表面高程变化最显著的地区。     

喜马拉雅地区叶如藏布流域冰川和冰湖变化遥感监测研究

叶如藏布流域冰川和冰湖众多,冰川融水是当地重要的淡水资源,是冰湖扩张的重要补给,冰湖溃决是当地潜在的自然灾害,因此分析该区域冰川和冰湖的现状与变化特征具有重要的现实意义。基于Landsat系列遥感影像,分析1990—2020年叶如藏布流域冰川和冰湖的分布与变化特征。结果表明：（1） 近30 a来叶如藏布流域冰川面积整体呈退缩趋势,由1990年167.80 km²退缩到2020年128.92 km²,共退缩38.88 km²,年均退缩率为0.77%·a^-1,且研究区冰川主要分布在海拔5800~6400 m之间,集中分布在5°~20°的坡度上。（2） 与冰川变化趋势相反,研究时段冰湖整体表现为扩张趋势,由1990年5.72 km²增加到2020年8.81 km²,30 a共增加3.09 km²,年均增长率为1.80%·a^-1。（3） 冰湖主要分布在海拔5000~5600 m范围内,坡度在0~10°分布面积较多,表碛覆盖型冰川与非表碛覆盖型冰川对冰湖有着不同程度的影响。（4） 1990—2017年叶如藏布流域温度与降水波动较大,温度整体呈上升趋势,降水量则波动下降,导致叶如藏布流域的冰川消融,冰湖扩张。通过上述研究以期为叶如藏布流域地区提供详细的冰川和冰湖面积分布与变化特征基础数据,为防灾减灾提供一定的支撑。     

基于能量守恒和3S技术的分布式冰川融水径流过程模拟

面向对象的分布式冰川融水径流模型的建立,不仅可以提高模拟计算的精度,而且可清楚地反映冰川融水径流各影响因子的特征,对研究中国西部地区流域水文循环具有重要意义.本研究借助3S技术,根据能量守恒和水量平衡原理,建立了包括栅格冰川融水、融水产流以及汇流子系统且具有严格物理意义的分布式冰川融水径流过程模型.通过DEM数字高程技...     

基于多源遥感数据的玛纳斯河流域冰川物质平衡变化

冰川物质平衡变化是连接气候和水资源的重要纽带,对河川径流有重要的调节功能。本文采用MOD11C3和TRMM 3B43等多源遥感数据驱动度日模型,模拟了2000—2016年玛纳斯河（简称玛河）流域冰川物质平衡过程,并分析了冰川融水对径流的补给规律。结果表明： ① 通过构建气温及降水反演模型能有效校正气象遥感原数据的精度,且经降尺度后能较精细刻画冰川区气候变化特征。冰川区年均气温和降水量分别为-7.57 ℃和410.71 mm,海拔4200 m处为气候变化剧烈地带,气温直减率以其为界上下分别为-0.03 ℃/100 m和-0.57 ℃/100 m,降水梯度分别为-2.66 mm/100 m和4.8 mm/100 m,海拔大于4700 m后降水又以5.17 mm/100 m递增。② 研究期内流域冰川持续呈负平衡状态,累积物质平衡达-9811.19 mm w.e.,年均物质平衡介于-464.85~-632.19 mm w.e.之间。垂向物质平衡在消融区和积累区分别以244.83 mm w.e./100 m、18.77 mm w.e./100 m递增。2000—2002年、2008—2010年冰川消融减缓,2002—2008年、2010—2016年消融加剧,其中2005—2009年期间冰川亏损最为强烈。③ 年内河川径流对冰川物质平衡变化响应强烈,尤以7月、8月物质平衡亏损最为严重占全年总量的75.4%,使得同期河川径流量占全年径流总量的55.1%。年际冰川融水补给率波动于19%~31%之间,可能是不同年份降水和积雪融水补给率差异较大所致。玛河与天山北坡其他河流冰川融水贡献率非常接近,也进一步证实了本研究物质平衡估算结果的可靠性。本研究可为其他流域冰川物质平衡研究提供借鉴和参考。     

东帕米尔高原冰川运动特征分析

冰川运动控制着冰量输送变化,为冰川变化和冰川灾害研究提供重要信息.为了探讨东帕米尔高原冰川运动特征及其影响因素,基于ITS_LIVE和GoLIVE分析了不同规模、不同地形条件、表碛/非表碛区域的冰川运动速度状况.研究结果表明:(1)东帕米尔高原冰川平均运动速度为5.31 m·a-1,冰川运动速度与冰川规模相关,表现为大...     

基于无人机影像的冰面流速与高程变化提取方法

无人机遥感能够获取高精度、高分辨率的冰川表面三维形态,有助于揭示冰川的运动和消融变化,现已成为冰川变化监测的重要技术支撑。针对无人机高分辨率影像能够有效识别表碛覆盖型冰川表面细节特征点的优势,本文引入金字塔影像集和最小二乘匹配方法联合开展特征点的追踪,用以实现冰川表面特征点三维位移场提取和动态变化分析。以贡嘎山贡巴冰川冰舌区作为典型研究区域,利用高分辨率无人机影像生成消融期前后（2018年6月9日和10月17日）的正射影像数据和数字表面模型,基于2期正射影像提取冰川表面特征点并实施追踪以获得三维位移场,在此基础上初步探明贡巴冰川冰舌段表面位移速率最大为7.51 cm/d,垂向消融速率超过11 cm/d。本文相关研究数据和结果可为冰川及复杂地形山区地质灾害的监测提供参考。     

Spatial distribution of supraglacial debris thickness on glaciers of the China-Pakistan Economic Corridor and surroundings

Debris-covered glaciers, characterized by the presence of supraglacial debris mantles in their ablation zones, are widespread in the China-Pakistan Economic Corridor (CPEC) and surroundings. For these glaciers, thin debris layers accelerate the melting of underlying ice compared to that of bare ice, while thick debris layers retard ice melting, called debris-cover effect. Knowledge about the thickness and thermal properties of debris cover on CPEC glaciers is still unclear, making it difficult to assess the regional debris-cover effect. In this study, thermal resistance of the debris layer estimated from remotely sensed data reveals that about 54.0% of CPEC glaciers are debris-covered glaciers, on which the total debris-covered area is about 5,072 km², accounting for 14.0% of the total glacier area of the study region. We find that marked difference in the extent and thickness of debris cover is apparent from region to region, as well as the debris-cover effect. 53.3% of the total debris-covered area of the study region is concentrated in Karakoram, followed by Pamir with 30.2% of the total debris-covered area. As revealed by the thermal resistance, the debris thickness is thick in Hindu Kush on average, with the mean thermal resistance of 7.0×10^-2 ((m²?K)/W), followed by Karakoram, while the thickness in western Himalaya is thin with the mean value of 2.0×10^-2 ((m²?K)/W). Our findings provide a basis for better assessments of changes in debris-covered glaciers and their associated hydrological impacts in the CPEC and surroundings.     

Glacier mapping based on Chinese high-resolution remote sensing GF-1 satellite and topographic data

The precise glacier boundary is a fundamental requirement for glacier inventory,the assessment of climate change and water management in remote mountain areas.However,some glaciers in mountain areas are covered by debris.The high spatial resolution images bring opportunities in mapping debris-covered glaciers.To discuss the capability of Chinese GaoFen-1 satellite lacking the short wave infrared band and thermal infrared band in mapping glaciers,this study distinguished supraglacial terrain from surrounding debris by combining GaoFen-1(GF-1)wide-field-view(WFV)images,the ratio of the thermal infrared imagery and morphometric parameters(DEM and slope)with 30 m resolution.The overall accuracy of 90.94%indicated that this method was effective for mapping supraglacial terrain in mountain areas.Comparing this result with the combination of GF-1 WFV and low-resolution morphometric parameters shows that a high-quality DEM and the thermal infrared band enhanced the accuracy of glacier mapping especially debris-covered ice in steep terrain.The user's and producer's accuracies of glacier area were also improved from 89.67%and 85.95%to 92.83%and 90.34%,respectively.GF data is recommended for mapping heavily debris-covered glaciers and will be combined with SAR data for future studies.     

Study of thermal properties of supraglacial debris and degree-day factors on Lirung Glacier,Nepal

The extensive debris that covers glaciers in the ablation zone of the Himalayan region plays an important part in regulating ablation rates and water availability for the downstream region. The melt rate of ice is determined by the amount of heat conducted through debris material lying over the ice. This study presents the vertical temperature gradients, thermal properties in terms of thermal diffusivity and thermal conductivity, and positive degree-day factors for the debris-covered portion of Lirung Glacier in Langtang Valley, Nepal Himalaya using field-based measurements from three different seasons.Field measurements include debris temperatures at different debris thicknesses, air temperature, and ice melt during the monsoon(2013), winter(2013), and pre-monsoon(2014) seasons. We used a thermal equation to estimate thermal diffusivity and thermal conductivity, and degree-day factors(DDF) were calculated from cumulative positive temperature and ice melt of the measurement period. Our analysis of debris temperature profiles at different depths of debris show the daily linear gradients of-20.81 °C/m, 4.05 °C/m, and-7.79 °C/m in the monsoon, winter, and pre-monsoon seasons, respectively. The values of thermal diffusivity and thermal conductivity in the monsoon season were 10 times greater than in the winter season. The large difference in these values is attributed to surface temperature and moisture content within the debris. Similarly, we found higher values of DDFs at thinner debris for the pre-monsoon season than in the monsoon season although we observed less melting during the pre-monsoon season. This is attributed to higher cumulative temperature during the monsoon season than in the pre-monsoon season. Our study advances our understanding of heat conductivity through debris material in different seasons, which supports estimating ice melt and discharge from glacierized river basins with debris-covered glaciers in the Himalayan region.     

INFLUENCE OF SUPRAGLACIAL DEBRIS ON SUMMER ABLATION AND MASS BALANCE IN THE 24K GLACIER,SOUTHEAST TIBETAN PLATEAU

Many temperate glaciers in the southeast Tibetan Plateau are covered by supraglacial debris in the ablation area. To evaluate the effect of such debris on summer ablation and mass balance, the surface ablation on the 24K Glacier was measured in the summer of 2008. Mean ablation rates varied from 10 to 52 mm/day, strongly correlated to debris thickness. Synchronous observations of air temperature allowed application of a simple degree‐day model to calculate ablation rates. Maximum values of both ablation rate and degree‐day factors appeared in the middle area of the glacier where the debris layer thickness was about 1 cm. The simulated daily ablation obtained from the degree‐day approach showed that the debris layers significantly affected the total summer ablation. The calculated ablation would be increased by 36% of the total ablation with the actual surface debris cover if glacier surface had been assumed to be entirely debris free. If completely covered by 65 cm thick debris in the ablation area, the glacier would experience a 59% decrease in summer ablation. The presence of a debris cover also leads to a change in the ablation gradient in the ablation zone.     

中国西部冰川度日因子的空间变化特征

在冰川与积雪消融研究中,度日模型应用较为广泛,该模型是基于冰雪消融与气温,尤其是正积温之间的线性关系建立的。度日因子是该模型的重要参数,反映了单位正积温产生的冰雪消融量,其空间变化特征对于不同模型模拟冰雪消融过程的精度有较大影响。 本文根据中国西部不同地区数十条冰川的短期考察和观测资料,分析了西部冰川度日因子的空间变化特征,结果表明：由于青藏高原及其周围地区独特的气候和热量条件,西部冰川度日因子具有明显的区域特征。在同一冰川上,度日因子的空间变化较为明显。从冰川类型来看,与极大陆型及亚大陆型冰川相比,海洋型冰川的度日因子较大。总体看来,西部冰川的度日因子由西北向东南逐渐增大,这与中国西部冰川的气候环境变化趋势是一致的,即在干冷的气候条件下,度日因子较小;而在暖湿的气候条件下度日因子较大。     

Estimate of Ablation Rate of Glacier Ice under a Supraglacial Debris Layer

The ablation rate under a debris layer is very difficult to measure directly at debris-covered glaciers, because the surface is highly heterogeneous, and the ablation rate varies tremendously from place to place. Heat budget considerations with a debris layer on top of glacier ice suggested that 'thermal resistance' of the debris layer could be estimated from surface temperature and the heat fluxes at the debris surface, and the ablation rate of the underlying glacier ice from the thermal resistance and meteorological data. The method was tested at the Lirung Glacier in Langtang Valley, Nepal Himalayas, using the thermal band of LANDSAT satellite for estimating surface temperature distribution of the debris top surface. The amount of melt water thus estimated was compatible with the observed discharge data from the glacier basin for periods of the monsoon season in 1985 and the pre-monsoon to the monsoon season in 1996. The investigation also revealed that the amount of discharge was much larger than the amount of precipitation over the basin, and it was suggested that the melt water from the debris-covered glacier contributes significantly to the river flow as a result of the shrinkage of the glacier.     

Thermal Properties of a Supraglacial Debris Layer with Respect to Lithology and Grain Size

This paper focuses on the impacts of debris cover on ice melt with regards to lithology and grain size. Ten test plots were established with different debris grain sizes and debris thicknesses consisting of different natural material. For each plot, values of thermal conductivity were determined. The observations revealed a clear dependence of the sub‐debris ice melt on the layer thickness, grain size, porosity and moisture content. For the sand fraction the moisture content played a dominant role. These test fields were water saturated most of the time, resulting in an increased thermal conductivity. Highly porous volcanic material protected the ice much more effectively from melting than similar layer thicknesses of the local mica schist. However, the analysis of thermal diffusivities demonstrated that the vertical moisture distribution of the debris cover must be taken into consideration, with the diffusivity values being significantly lower in deeper layers.     

Glaciers in Svalbard: mass balance, runoff and freshwater flux

Gain or loss of the freshwater stored in Svalbard glaciers has both global implications for sea level and, on a more local scale, impacts upon the hydrology of rivers and the freshwater flux to fjords. This paper gives an overview of the potential runoff from the Svalbard glaciers. The freshwater flux from basins of different scales is quantified. In small basins (A < 10 km²), the extra runoff due to the negative mass balance of the glaciers is related to the proportion of glacier cover and can at present yield more than 20% higher runoff than if the glaciers were in equilibrium with the present climate. This does not apply generally to the ice masses of Svalbard, which are mostly much closer to being in balance. The total surface runoff from Svalbard glaciers due to melting of snow and ice is roughly 25 ± 5 km³ a⁻¹, which corresponds to a specific runoff of 680 ± 140 mm a⁻¹, only slightly more than the annual snow accumulation. Calving of icebergs from Svalbard glaciers currently contributes significantly to the freshwater flux and is estimated to be 4 ± 1 km³ a⁻¹ or about 110 mm a⁻¹.     

Assessment of the Evolution in Velocity of Two Debris‐Covered Valley Glaciers in Nepal and New Zealand

Feature tracking of orthorectified pairs of Advanced Spaceborne Thermal Emission and Reflection Radiometer satellite images is used to calculate velocities for the Tasman Glacier, New Zealand (2002–2014) and the Khumbu Glacier, Nepal (2001–2008). Velocities in the middle and upper ablation zones of both glaciers show a long‐term decrease of ～10–20%, while the terminus of Khumbu Glacier has remained near stagnation throughout the study period. In contrast, there has been a recent acceleration of the lower terminus of Tasman Glacier, from ～5 m a^–1 in 2002 to 40 m a^–1 in 2014. Both of these glaciers have an extensive supraglacial debris cover across their lower ablation regions, with the Khumbu Glacier terminating on land and the Tasman Glacier terminating in a proglacial lake. The rapid recent increase in velocity of the terminus of Tasman Glacier is closely correlated with the increase in size of its proglacial lake. These results indicate the complex dynamic changes that mountain valley glaciers may undergo in response to long‐term negative mass balance.