表碛下冰面消融模型的改进

针对已有的表碛下冰面消融模型的缺点,通过对模型算法、模型假设和模型结构等方面进行改进,在模型的易用性和模拟精度上有了较大的提高.科其喀尔冰川的实例研究表明,改进模型对于表碛区不同试验点上地温的模拟是较好的,变化趋势基本一致,变化幅度也相当.但也存在模拟值与实测值之间变化相位差异的突出问题,实际操作可以通过适当的校正程序得以较好的解决.模型中对表碛热参数的常数化处理不会产生较大的模拟误差,此外计算开始时均一表碛温度的假设不会对表碛温度的整体模拟效果产生较大的影响.通过地温模拟得到的冰面消融速率与利用实测表碛热通量计算的冰面消融速率的比较说明,模型思路合理且算法可行.     

西天山托木尔峰南麓大型山谷冰川冰舌区消融特征分析

基于对托木尔峰南麓托木尔型山谷冰川的野外考察和典型冰川的定位观测,对冰面被表碛广泛覆盖的所谓“托木尔型”冰川冰舌区表碛与冰面消融的关系进行了研究. 结果表明：表碛对冰面消融、冰川水文过程、冰川变化等均具有重要影响,当表碛厚度超过3 cm时,表碛对冰面消融就产生明显抑制作用,且随着厚度增加,冰面消融显明减弱. 科其喀尔冰川表面的观测表明,由末端向上,表碛厚度逐渐减薄. 受表碛影响,科其喀尔冰川区最大的消融量出现在海拔3 800~3 900 m之间、表碛物厚度小于10 cm的区域内;冰川消融强度由此向上随着海拔的升高而下降,向下随表碛厚度的增大而减弱. 冰面湖的发育是表碛覆盖冰川的又一主要特征,湖水对冰面的融蚀和快速排泄成为冰面产汇流的主要过程. 科其喀尔冰川研究表明,两三个冰面湖排泄形成的融蚀冰量就相当于冰川末端退缩造成的冰量损失. 因此,冰面湖等热喀斯特地形的形成、扩张融蚀、融穿排泄、形成湖区低地,这一周而复始的过程不仅是其主要消融方式之一,而且也强烈的影响着冰川水文及冰川变化. 托木尔峰南麓地区大型冰川变化主要以厚度减薄为主,而不是像大多数冰川显著的变化主要表现在末端和面积减少方面.     

天山南坡科契卡尔巴西冰川消融期雪面能量平衡研究

根据2005年6～9月的野外观测资料,计算了天山南坡科契卡尔巴西冰川积雪表面的能量平衡各分量,其中感热和潜热采用空气动力学-梯度方法计算得到,净辐射由观测获得.结果表明,消融期净辐射是雪面最主要的能量来源,占能量收入的81.4%,平均值为63.3 W·m^-2;其次为感热供热,占18.6%,为14.4 W·m^-2.吸收的热量主要通过融化和蒸发两种方式消耗,融化和蒸发耗热分别为54.0 W·m^-2和23.0 W·m^-2,占能量总支出的69.5%和29.7%,剩下的0.8%由感热消耗.在积雪表面的能量组成中,感热值在6月和9月较大而8月较小;对于潜热来说,6月潜热交换绝对值最大,蒸发最强烈,这与6月风速大且天气晴好有关.另一方面,净辐射值较大的6月和9月融化热较多,说明此时段冰川消融较为强烈.     

喜马拉雅山珠穆朗玛峰北坡绒布冰川消融速率特征分析

根据1959年和2009年在喜马拉雅山珠穆朗玛峰北坡绒布冰川获得的冰川消融数据, 分析了该冰川消融速率变化特征.结果表明: 1) 在珠峰绒布冰川表碛覆盖区, 表碛厚度随海拔升高而降低. 2) 不同厚度表碛下的冰川消融速率差别较大; 当表碛厚度>8.5 cm时, 消融速率随表碛厚度的增加而减小; 促进冰川消融的表碛厚度阈值大于5 cm. 3) 从冰川消融速率的空间分布看, 绒布冰川大部分区域的消融速率<20 mm·d^-1, 最大消融速率出现在海拔5 400~5 450 m处. 4) 绒布冰川消融速率受表碛厚度和气温综合影响, 低海拔处表碛太厚, 高海拔处气温较低, 冰川消融在上述两海拔处均受抑制, 冰川消融速率较小; 在中海拔处, 表碛相对较薄, 气温相对较高, 冰川消融速率最大; 冰川日均消融速率与日均正积温正相关. 5) 喜马拉雅山南坡冰川消融速率大于北坡冰川消融速率.     

青藏公路沿线多年冻土区土壤热通量参数化方案的优化和检验

土壤热通量是地表能量平衡的重要分量,其估算方案在研究地表能量平衡研究中必不可少。利用青藏公路沿线5个站点0~20 cm的实测土壤层温、湿度及5 cm土壤热通量资料,以翁笃鸣气候学计算方案为基础建立了优化的5 cm土壤热通量计算方案。通过唐古拉和西大滩两个独立站点的检验结果表明,优化方案的结果相对于原方案有较大的改善,唐古拉和西大滩5 cm土壤热通量均方根误差值分别减小了3.2 W·m^-2和4.8 W·m^-2,而相对误差分别减小了61.9%和36.1%,即新方案能够较好地估算出青藏公路沿线多年冻土区5 cm土壤热通量。使用优化方案模拟了青藏公路沿线11个站点5 cm土壤热通量变化,结果显示,近十年青藏公路沿线土壤热通量呈现出增大的趋势,其中,5 cm土壤热通量增大了近1.0 W·m^-2,而且各观测场的年平均土壤热通量值均大于0.0 W·m^-2,表明就年尺度而言,热量有盈余,盈余热量用于加热下层土壤,引起活动层厚度增加,平均状况下土壤热通量每增大1.0 W·m^-2,活动层厚度增大约21.0 cm。     

慕士塔格峰洋布拉克冰川消融的观测分析

2001年7月4日至8月8日,在慕士塔格峰西侧的洋布拉克冰川海拔4600～4460m区间的冰舌段,进行了短期的冰面消融观测.慕士塔格峰冰川区暖期短,冰面强消融时期比较集中.观测期间,冰面纯消融厚度为640～1260mm水层,日平均消融厚度达26～39.6mm,推算冰舌区年消融量不低于1700～2000mm,比青藏高原内部的冰川消融强烈的多.7月21-22日出现最大消融值,在海拔4460m和4600m,日消融量分别为144.5mm和59.5mm.冰面消融随海拔上升而减小,日平均消融梯度:在裸露冰区为0.40～0.55mm·10m-1;在表碛覆盖区为0.21～3.53mm·10m-1,变幅较裸露冰区大.按裸露冰区的消融梯度计算出海拔4800m处的日平均消融量,和过去的研究资料比较,2001年冰面日平均消融量较1987年和1960年的消融量大,反映出慕士塔格峰区影响冰川消融的气候与全球气候变暖的特点是一致的.     

天山南坡科其卡尔巴契冰川度日因子变化特征研究

度日模型是估算冰川消融的一种简单而有效的方法.根据科其卡尔巴契冰川2003年的观测资料,分析了该冰川度日因子的空间变化规律及其影响因素.研究表明:各高度上的度日因子,介于2.0～9.7mm·℃^-1·d^-1之间变化,平均值为5.7mm·℃^-1·d^-1,与青藏高原各冰川及其它地区冰川相比较小;随着海拔的增高,度日因子随之递增;随平均气温的升高而随之递减.由于冰面状况复杂,度日因子变化幅度较大,裸冰区的度日因子明显大于表碛覆盖区.人为测量误差、反照率、地形等对度日因子的影响也不容忽视.     

藏东南嘎隆拉冰川表碛冻融过程与零点幕效应

冰川表碛区冻融过程的观测分析有助于厘清各层之间能量和水分传输的关系,从而为构建相应的能量平衡模型以及冰川径流模型提供理论支持。基于2015年10月至2016年11月嘎隆拉冰川表碛区自动气象站气象资料,对表碛区冻融过程进行分析,结果表明：表碛中出现冻土中常见的秋季和春季零点幕效应,秋季零点幕持续时间32天,春季零点幕持续58天;春季零点幕期间,气温和地表温度均呈现清晨低,午后高的日变化特征,而表碛层内温度没有明显的日变化特征,其体积含水率呈“脉冲式”变化;春季零点幕效应结束,底层冻结层破坏后,冰川冰消融才开始;嘎隆拉冰川表碛对冰川冰消融的影响不仅是因为表碛较厚,在消融量上抑制了冰川冰的消融,而且更重要的是表碛产生零点幕效应,延迟了表碛下冰川冰开始消融的日期,在时间上抑制了冰川冰消融。     

基于KNN机器学习方法对青藏高原唐古拉地区表层土壤水热状况的模拟

利用唐古拉站2004—2012年气象观测资料,基于KNN算法,结合机器学习思想,建立了一个气象回归模型,模拟了2005年唐古拉地区表层土壤水热变化趋势,结合实测数据,将模拟值与观测值进行对比,并对模型模拟效果进行了评估。结果表明：KNN模型能够较好地模拟活动层表层土壤水热状况,各层土壤温度的模拟值与观测值的相关系数均在0.99以上,均方根误差在1.25 ℃以内;不同深度土壤水分的模拟值与观测值的相关系数均在0.95以上,均方根误差在0.02 m³?m^-3以内。总体上,KNN模型能够对青藏高原多年冻土区唐古拉地区表层土壤水热状况进行较为精确地模拟,该模型对于青藏高原其他地区的适用性有待进一步研究验证。     

黄河源鄂陵湖地区辐射收支和地表能量平衡特征研究

利用2010年6-7月鄂陵湖野外试验的近地层观测数据,分析了在不同天气条件下黄河源鄂陵湖地区辐射分量、地表能量分量、土壤温度和反照率的变化特征. 结果表明：不同天气条件下,辐射和地表能量各分量日变化差异较大,晴天、阴天和雨天的地表反照率依次递减,平均反照率约为0.21;观测期内,平均辐射贡献从大到小依次为向上长波、向下长波、向下短波、向上短波,日积分值分别为31.4 MJ·m^-2、25.6 MJ·m^-2、22.4 MJ·m^-2、4.2 MJ·m^-2,净辐射(12.5 MJ·m^-2)占向下短波辐射的55.7%;平均地表能量和土壤温度的变化幅度较晴天小,感热、潜热、0 cm土壤热通量的平均日积分值分别占净辐射的21.2%、43.1%、8.2%;平均土壤温度变化幅度随深度增加逐渐减小,浅层土壤温度峰值较晴天低2 ℃,深层土壤温度相差不大. 云和降水的扰动削弱了向下短波辐射,导致平均感热通量和0 cm土壤热通量的峰值比晴天小,而平均潜热通量的峰值大于晴天. 由于湖泊水体巨大的热容量和水分供应,鄂陵湖地区的气温日较差较小,地表温度变化幅度变小,附近地表温度升高缓慢. 鄂陵湖区的地表能量平衡中,潜热通量占主导,感热和地表土壤热通量次之. 研究结果有助于理解气候变化背景下黄河源区湖泊的能量水分循环过程,为促进该地区光热资源的合理利用和畜牧业的可持续发展提供数据支持.     

科奇喀尔冰川夏季表碛区热量平衡参数的估算分析

利用能量平衡原理、热传导理论和通量传输理论建立了一个热量平衡参数的估算模型,对西天山的科奇喀尔冰川夏季消融区中部表碛区的热量平衡参数进行估算与分析.结果表明:净辐射是表碛面热量收支的主要热源,吸收的热量主要以潜热和感热的形式向大气输送水汽和热量,剩余部分用于表碛增温耗热.与消融区上部的冰面和表碛面相比,在消融区中部表碛面热量收入中感热输送减小,同时向上的地热输送增加.热平衡支出项中,感热交换、蒸发耗热和地热通量的比例分别为39.1%、39.9%和21%,其中感热通量与蒸发耗热的比例比消融区上部有所提高,蒸发耗热的增加比较显著.在总的热量支出中,平均只有7.8%的热量可以用于表碛下部的增温和向深层传导.     

Thermal regime of a supraglacial lake on the debris-covered Koxkar Glacier, southwest Tianshan, China

A supraglacial lake was surveyed on the Koxkar Glacier in southwest Tianshan from July to September 2007 and July to September 2008, and the temperature variation characteristics of the lake, debris and debris-free ice were analyzed at different depths to determine the thermal regimes. In addition, the discrepancies of temperature variation characteristics were investigated for different geomorphic units of the ablation zone of the Koxkar Glacier. It was found that daily temperature variation curves for deep water are V-shaped because meltwater from the glacier surface at temperatures of around 0°C feeds the lake and mixes with the relatively high-temperature surface water during the day. As the water temperature rises to approximately 4°C, the mixed water sinks and forms a low-temperature trough in the deep water of the lake in the middle of the day. The vertical lapse rate of the lake water temperature against depth (?0.33°C/m) has a magnitude lower than that of the debris (?4.29°C/m) and that of the debris-free ice (?0.38°C/m) in the Koxkar Glacier??s ablation zone. The temperature curve for the surface water largely varied between the temperature curves for the debris at depths of 0.2 and 0.5?m. The surface thermal condition of the ablation zone is significantly affected by the daily weather, and there is a limited influence in debris at a depth of 1?m and in the lake at a depth of 5?m.     

Imaging the debris internal structure and estimating the effect of debris layer on ablation of Glacier ice

This study is a comparative investigation of the debris layer and underlying ice of the Koxkar Glacier using multi-frequency GPR with antennae having different frequencies. Together with analysis of the fluctuation of the radar signal amplitude and polarity, the debris layer and underlying ice were analyzed on the basis of high-resolution GPR images. It was found that the optimal average velocity in the shallow layer (0–4 m) is 0.06 m/ns. Images obtained with different frequency antennas have different characteristics; and the performance of the 200 MHz antenna for a debriscovered glacier is the best. The interpretation of typical GPR image is validated by using FDTD numerical model. Combining the debris layer thickness and the underlying ice structure, the effect of debris layer on ablation of glacier ice and forecast of the glacier change in the aspect of thickness-thinning and glacier retreat can be estimated. This study can provide as a reference to the formation mechanisms and estimation of the ice volume of glaciers covered by debris.     

Estimate of heat flux at the ice-water interface in Lake Baikal from experimental data

An original system for measuring temperature in the ice cover and subglacial water and an increase in the ice thickness provides data necessary for calculation of the heat flux at the ice-water interface. Successive freezing of 1-mm temperature sensors during the ice growth allows us to measure temperature gradients in the vicinities of the ice-water interface for the first time. An analytical equation derived from the Stefan condition allows calculations of the heat flux at the phase boundary on the basis of the experimental data, which agree with independent estimates that have been made on the basis of the subglacial temperature gradients and are within the 4–39 W/m² range. The flow at the ice-water interface is comparable with the heat flux inside the ice depth and significantly affects the dynamics of the ice cover thickness.     

Debris control on glacier thinning—a case study of the Batal glacier,Chandra basin,Western Himalaya

Complex factors such as climate, glacial geometry, topographical features and debris covers have significant influence on the dynamics of the Himalayan glaciers. Presence of debris covers on the surface of glaciers can significantly alter the surface energy balance and influence the climatic response of glaciers. In this study, the influence of debris covers and its impact on the ablation processes were analyzed from the in situ data collected over the surface of the Batal glacier in Chandra Basin, Western Himalaya. Almost 90 % of the ablation zone of the Batal glacier is covered by debris, 35 % of which is thick debris (>10 cm). Fourteen stakes (depth ～10 m) with increasing altitude and with varying debris thicknesses were installed to cover the whole ablation zone. Among them, four stakes represent thin debris (<2 cm), two stakes represent 2–5 cm debris thickness, two stakes represent 5–25 cm debris thickness, three stakes represent 25–50 cm debris thickness and three stakes represent >50 cm debris thickness. Our study has revealed high surface melting (?2.0 cm. w.e.d^?1) in the debris free glacier while low surface melting observed in thick debris covered ice (?0.6 cm. w.e.d^?1). Although limited to one season, this observation revealed a significant difference in the rate of surface melting as per the increasing debris thickness. Contrasting to normal ablation pattern over glaciers, Batal has experienced inverse retreat rate of ablation along with increasing altitude. A high degree of negative correlation (r = ?0.82, p < 0.05) between ablation rate and debris thickness in Batal suggest a significant control of debris thickness over ablation rate.     

Role of debris cover to control specific ablation of adjoining Batal and Sutri Dhaka glaciers in Chandra Basin (Himachal Pradesh) during peak ablation season

As part of the on-going annual mass balance measurements on Batal and Sutri Dhaka glaciers, observations were made during peak ablation (August–September) season in 2013 to understand the response of debris covered and clean-ice (debris free) glacier surface to melting processes. Though, both the Batal and Sutri Dhaka glaciers have almost similar geographical disposition, Batal shows extensive debris cover (90% of the ablation area), while the latter is free from debris (only 5% of the ablation area). The thickness of debris in Batal glacier is inversely proportional to altitude, whereas Sutri Dhaka mostly experienced debris-free zone except snout area. Observation revealed that the vertical gradient of ablation rate in ablation area is contrastingly opposite in these two glaciers, reflecting significant control of debris thickness and their distribution over glacier surface on the ablation rates. While different thickness (2–100 cm) of debris have attenuated melting rates up to 70% of total melting, debris cover of <2 cm thickness has accelerated melting up to 10% of the total melting. Estimated melt ratio reveals that about 90% of the ablation area has experienced inhibited melting in Batal glacier, whereas only less than 5% ablation area of Sutri Dhaka has undergone inhibited melting. Comparison of topographical maps of 1962 with successive satellite images of the area demonstrates a terminus retreat of 373 ± 33.5 m and 579 ± 33.5 m for Batal and Sutri Dhaka glaciers for the period 1962–2013, respectively.