冰川冻土, 2021, 43(1): 24-35 doi: 10.7522/j.issn.1000-0240.2019.1002

冰冻圈与全球变化

天山乌鲁木齐河源1号冰川东支能量-物质平衡模拟研究

李宏亮,1,2, 王璞玉,1,2, 李忠勤1,3, 金爽1, 徐春海1, 梁鹏斌3, 岳晓英1, 羊旻1

1.中国科学院 西北生态环境资源研究院 冰冻圈科学国家重点实验室/天山冰川观测试验站, 甘肃 兰州 730000

2.中国科学院大学, 北京 100049

3.西北师范大学 地理与环境科学学院, 甘肃 兰州 730070

Study on the energy-mass balance simulation of the east branch of the Urumqi Glacier No. 1, Tianshan Mountains

LI Hongliang,1,2, WANG Puyu,1,2, LI Zhongqin1,3, JIN Shuang1, XU Chunhai1, LIANG Pengbin3, YUE Xiaoying1, YANG Min1

1.State Key Laboratory of Cryospheric Science / Tianshan Glaciological Station, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China

2.University of Chinese Academy of Sciences, Beijing 100049, China

3.College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, China

通讯作者: 王璞玉,研究员,主要从事冰川变化及其机理研究. E-mail:wangpuyu@lzb.ac.cn

编委: 武俊杰

收稿日期: 2019-01-08   修回日期: 2019-05-30  

基金资助: 国家自然科学基金项目.  41771077
国家自然科学基金国际(地区)合作研究项目.  41761134093
中国科学院战略性先导科技专项(A类).  XDA20020102

Received: 2019-01-08   Revised: 2019-05-30  

作者简介 About authors

李宏亮,硕士研究生,主要从事冰川能量-物质平衡模拟研究.E-mail:lhl_0922@yeah.net , E-mail:lhl_0922@yeah.net

摘要

基于天山乌鲁木齐河源1号冰川东支海拔4 025 m处自动气象站的观测数据和同期物质平衡花杆观测数据,采用COSIMA模型,对该冰川东支2018年消融期单点能量-物质平衡进行了模拟。结果显示:物质平衡模拟值为(-0.67±0.03)m w.e.,与实测值有非常好的一致性,相关系数达0.96。造成冰川消融的能量来源于净短波辐射(84%)、感热通量(16%);冰川能量支出为净长波辐射(55%)、冰川消融耗热(32%)、潜热交换(7%)及地热通量(6%)。受能量收支影响,模拟物质平衡主要取决于表面消融和固态降水。与我国其他区域大陆型冰川研究结果比较发现,乌鲁木齐河源1号冰川物质损失较为显著,能量通量主要取决于海拔以及气候条件,再冻结和固态降水显著小于羌塘1号冰川和扎当冰川,推测与单条冰川所处的大气环流有关。

关键词: 乌鲁木齐河源1号冰川 ; 消融期 ; 能量平衡 ; 物质平衡 ; 模拟

Abstract

Glacier ablation can be described using energy-mass balance models. This paper conducts an energy budget experiment to quantify energy components in the Urumqi Glacier No.1, Tianshan Mountains. Based on the automatic weather station (4 025 m a.s.l.) observation data and mass balance field survey data by five individual ablation stakes on the east branch of the glacier, the single-point energy-mass balance was simulated by “COupled Snowpack and Ice surface energy and MAss balance model” (COSIMA) during the ablation season in 2018. Results showed that the modelled accumulative mass balance was (-0.67±0.03) m w.e. during the study period, which agreed well with the in-situ observation. The correlation coefficient between both reached to 0.96. The main energy component causing the glacier ablation were net shortwave radiation (84%) and sensible heat flux (16%). The main energy expenditure were the net longwave radiation (55%), the heat flux for snow/ice ablation (32%), latent heat flux (7%) and ground heat flux (6%). Affected by energy budget, modelled accumulative mass balance mainly depended on surface melt and snowfall. Comparing with the other continental glaciers in China, it is found that mass loss of the Urumqi Glacier No.1 was more significant and that energy fluxes mainly depends on altitude and climatic conditions. In addition, the refreezing and snowfall on the Urumqi Glacier No.1 were significantly less than that on the Qiangtang No.1 Glacier and Zhadang Glacier, probably due to the atmospheric circulation over the single glacier. These findings will enhance our understanding of the mechanisms of glacier changes.

Keywords: Urumqi Glacier No.1 ; ablation season ; energy balance ; mass balance ; simulation

PDF (2921KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李宏亮, 王璞玉, 李忠勤, 金爽, 徐春海, 梁鹏斌, 岳晓英, 羊旻. 天山乌鲁木齐河源1号冰川东支能量-物质平衡模拟研究[J]. 冰川冻土, 2021, 43(1): 24-35 doi:10.7522/j.issn.1000-0240.2019.1002

LI Hongliang, WANG Puyu, LI Zhongqin, JIN Shuang, XU Chunhai, LIANG Pengbin, YUE Xiaoying, YANG Min. Study on the energy-mass balance simulation of the east branch of the Urumqi Glacier No. 1, Tianshan Mountains[J]. Journal of Glaciology and Geocryology, 2021, 43(1): 24-35 doi:10.7522/j.issn.1000-0240.2019.1002

0 引言

冰川(包括冰盖)作为冰冻圈要素之一,因其对气候的高度敏感性和重要的反馈作用而倍受关注。随着气候变暖,冰川变化对海平面、淡水资源、生态、环境、气候等产生了深刻的影响。根据IPCC第五次评估报告,末次盛冰期以来,全球海平面平均上升了120 m,目前全球山地冰川对海平面上升的贡献为(0.76±0.37) mm⋅a-1,是仅次于海洋热膨胀的海平面上升第二大贡献者1。近50年来,我国天山山脉冰川面积缩小了18%2,冰川的退缩对区域生态环境、绿洲经济发展、工农业生产以及居民用水等造成至关重要的影响3-4

冰川物质平衡是表征冰川积累和消融量值的重要冰川学参数之一,主要受控于能量收支状况,对气候变化有敏感的响应。物质平衡及其动态变化是引起冰川规模和径流变化的物质基础,是连接冰川与气候、冰川与水资源的重要纽带5。传统的物质平衡监测主要利用花杆/雪坑法,但受限于高海拔、陡峭地势,耗时费力。近年来利用不同时期数字高程模型(DEM)差值来计算物质平衡,使得评估大尺度、长时间序列的冰川物质平衡成为可能6,但在认识冰川变化的物理机制及对气候变化响应过程和机理方面存在一定的限制。

目前,国内外冰川学聚焦于冰川物质平衡模拟。模型模拟主要有两种方法。一种是采用半经验的方法,建立冰川物质平衡和气象指标的关系,如线性回归模型7、度日因子模型8等。Rasmussen等7利用线性回归模型在斯瓦尔巴群岛的模拟结果指出,该模型对冬季平衡的模拟效果良好,同时该模型操作简单能够重建物质平衡,已经在挪威、冰岛和瑞典等地得到了应用。Huintjes等9利用修正的度日因子模型在乌鲁木齐河源1号冰川东支计算了物质平衡,发现在引入潜在太阳短波辐射修正模型后,模拟的物质平衡和消融速率有了显著的改善,但该方法只涉及表面物质平衡的计算,在气候条件复杂的大区域中是不适用的10。另一种方法是基于能量平衡描述冰川物理过程的物质平衡模型,如分布式能量-物质平衡模型11和表层能量-物质平衡模型(SEMB)12。Hock等11通过改进反照率参数化方案以及辐射分量计算方法,考虑地形因素如坡度、坡向等,建立的分布式能量-物质平衡模型能够对小时尺度的冰雪消融进行准确模拟。与度日因子模型相比,该方法则更为细致地体现了物质平衡的计算。此外,Yang等13、Li等14、Zhu等15利用表层能量-物质平衡模型在藏东南地区进行了系统的模拟研究,关注云量和反照率对物质平衡的影响,同时结合区域大气环流如南亚季风,探讨其与物质平衡的关系,进一步在区域尺度上分析了能量-物质平衡的差异性,从而揭示冰川对区域气候的响应过程。这些模型仅考虑了冰川表面的能量平衡分布,而由德国柏林洪堡大学研制开发的COSIMA模型,则充分考虑了大气、冰川表面以及表面以下10 m内雪层的物质和能量交换过程,通过建立物质平衡与气象要素(如气温、降水、辐射等)变化之间的关系,实现对整条冰川规模、小时尺度的冰川物质平衡模拟计算。目前该模型已经在青藏高原扎当冰川16和普若岗日冰帽17等地进行了能量-物质平衡模拟应用,取得了很好的模拟效果。

鉴于此,基于乌鲁木齐河源1号冰川东支2018年消融期冰面气象站和物质平衡花杆观测数据,利用COSIMA模型,开展了冰川能量-物质平衡模拟研究,以期为模拟预估由冰川变化引发的水资源与水循环、生态环境等变化,同时揭示冰川消融机理,为研发延缓冰川消融、保护冰川的措施提供科学依据。

1 研究区概况

乌鲁木齐河源1号冰川(43°06′ N、86°49′ E)位于我国天山中部喀拉乌成山脉主脉北坡乌鲁木齐河源上游(图1),属冰斗山谷冰川,由于冰川消融,于1993年分离为东、西两支。2015年1号冰川总面积为1.558 km2,其中东支面积为1.058 km2,海拔介于3 743~4 267 m18。受夏季气温升高、反照率下降、辐射增强及冰川温度升高影响,近50年来1号冰川在面积、末端位置及冰川温度等方面发生了显著变化,其中末端位置和面积变化率分别为4.2 m⋅a-1、0.006 km2⋅a-1[19-20。该区域主要受对流层上部西风急流、西伯利亚反气旋环流和西风环流气旋扰动三大动力要素控制2。大西沟气象站(海拔3 593 m)观测数据和物质平衡观测数据显示,1959—2017年乌鲁木齐河源1号冰川年均气温为-4.6 ℃,年均降水量为460 mm。5—8月是主要的消融期,降水主要来源于西风携带的水汽,其中5—8月降水占年降水量的78%,降水类型主要为固态降水,如雪、雹及雨夹雪21

图1

图1   乌鲁木齐河源1号冰川观测示意图

Fig.1   The map showing the location of the Urumqi Glacier No.1 (a) and the observation system on the east branch (b), together with the photos of automatic weather station (AWS) (c), weighing precipitation/snow gauge (d) and ablation stake (e)


2 数据与方法

2.1 气象数据

利用冰面自动气象站AWS1获取消融期(2018年4月29日14:00至9月1日10:00,北京时间)小时尺度气象数据(图1)。该气象站架设于海拔4 025 m处,接近东支主流线,地势相对平坦(坡度<2°)。观测的气象要素包括气温、相对湿度、气压、风速、降水量、入射和出射的长短波辐射等,数据集存储在Campbell CR1000中,每十分钟记录一个气象要素。具体气象站传感器类型和指标见表1

表1   自动气象站传感器信息

Table 1  Sensor information of the AWSs used in this study

气象站点传感器型号生产国/生产商测量项目及单位精度架设高度
AWS1,AWS2,AWS3HC2-S3瑞士/Rotronic气温/℃±0.1 ℃(23 ℃)2 m
相对湿度/%±0.8%(23 ℃)2 m
AWS1,AWS3Young 05103美国/RM Young风速/(m⋅s-1±0.3 m⋅s-12 m
AWS1,AWS3Geonor T-200B挪威/Geonor降水/mm±0.1 mm(-40~60 ℃)1.7 m
AWS1,AWS3PTB110美国/Campbell气压/hPa±0.3 hPa(20 ℃)2 m
AWS1,AWS3CNR4荷兰/Kipp&Zonen短波和长波辐射/(W⋅m-2±10%(-40~80 ℃,0%~100%RH1.5 m

新窗口打开| 下载CSV


然而,在复杂冰面气象环境的影响下,AWS1在观测期间出现了48 h的数据缺失,因此,为了填补这一数据空白来获得连续的气象数据集,本文利用AWS2和AWS3相应时间段内的数据来进行补充,具体方法见2.3节。AWS2是安装在海拔3 882 m花杆D2处距离冰面高度2 m的便携式气象仪。AWS3位于海拔3 835 m冰川监测塔处,距离1号冰川末端约100 m,该气象站具备同时段30 min尺度的入射及出射长短波辐射、气温、相对湿度、风速及降水量的观测记录。考虑到辐射传感器顶部会因频繁的降雪而形成积雪和霜,从而影响短波辐射的观测,因此,本文利用van den Broeke等22的方法对短波辐射进行修正,该方法则可以使净短波辐射的计算精度提高至小时尺度。基于临界气温的参数化方法来判断降水的状态,COSIMA中采用双临界气温指数法进行固液态降水分离:当气温高于5 ℃时为降雨,低于1 ℃时为降雪,气温介于两临界气温之间则按线性插值计算23。考虑到风速对固态降水的影响,采用T-200B原始降水数据的修正方法24,即

CE=exp(- 0.2Ua)

式中:CE为捕捉率(%);Ua为风速(m⋅s-1)。一般认为液态降水的捕捉率为90%,进行液态降水的修正。

2.2 冰川物质平衡数据

为了进行冰川物质平衡的验证,在AWS1同海拔附近布设了5个花杆点(图1)。2018年消融期,每隔半个月进行了详细的观测,选取消融期的7次,分别是6月20日、7月20日、8月1日、8月10日、8月14日、8月18日以及9月1日。具体每根花杆物质平衡观测内容包括:花杆至冰川表面的垂直高度、附加冰的厚度、粒雪层的厚度和密度及雪坑剖面结构。此外,密度测量必须要求雪坑垂直高度大于5 cm,否则冰雪密度采用野外观测经验值,冰和雪的密度分别取值为900 kg⋅m-3和300 kg⋅m-3。单点物质平衡计算公式为

bn=bs+bice+bsi

式中:bs、bice、bsi分别为雪、冰川冰以及附加冰的物质平衡,具体参见文献[25]。由此获得5根花杆的单点物质平衡值,取其平均值作为冰川物质平衡的验证数据。

然而,在野外观测过程中,由于气候和人为等因素会造成物质平衡的观测存在误差,如观测人员花杆读数、雪/粒雪密度的测量及假设、花杆的倾倒及下沉、雪/粒雪层的误判、附加冰难以辨别以及重复测量等26-27。综上所述,据Thibert等26提出的单点物质平衡误差计算方法,该方法综合考虑了单点物质平衡测量时因雪/粒雪层误判以及重复测量造成的误差,最终确定冰川物质平衡计算结果误差为±0.16 m w.e.。

2.3 COSIMA模型

基于COSIMA(COupled Snowpack and Ice surface energy and MAss balance model)开展1号冰川的单点物质平衡模拟,该模型是由表面能量平衡模型和表层下雪冰模型耦合形成的。其中,表面能量平衡模型主要基于如下能量平衡等式。

F=SWin(1-α)+LWin+LWout+Qsens+Qlat+QG

式中:SWinLWinLWoutQsensQlatQGα分别为入射短波辐射、入射与出射长波辐射、感热与潜热通量、地热通量以及反照率。需要说明的是,只有当冰川表面温度达到消融点时消融才会发生,此时F的数值就是各能量项的代数和,规定方向指向冰川表面的能量通量为正值。

COSIMA中利用块体空气动力学方法计算湍流热通量28,则冰川表面感热通量和潜热通量可表示为

Qsens=ρairK2lnhz/zo2uTair-Ts
Qlat=ρairLE/SK2lnhz/zo2uqair-qs
ρair=p100/287.058Tair1+qair0.608
qair/s=RHair/s0.622E/sp-E/s/100

式中:ρair为空气密度(kg⋅m-3);K为卡尔曼常数(0.4);hz为仪器架设高度,本文为2 m;zo为表面粗糙长度,新雪、粒雪和冰的粗糙长度分别为0.24 mm、4 mm及1.7 mm29-31u为风速(m⋅s-1);TairTs分别为2 m处和冰川表面气温(℃);LE为蒸发潜热(2.514×106 J⋅kg-1);LS为升华潜热(2.849×106 J⋅kg-1);qairqs分别为2 m处和冰川表面比湿(kg⋅kg-1);RHair/s为2 m处和冰川表面相对湿度(%),其中冰川表面的相对湿度假设为100%;p为气压(hPa);EEs分别为2 m处和冰川表面饱和水汽压(hPa)。

针对湍流通量计算的不确定性,COSIMA利用Braithwaite32描述方法对湍流通量进行稳定性修正,如下所示。

Qsens/lat=Q'sens/latRi≤0.01)
Qsens/lat=Q'sens/lat1-5Ri2 (0.01<Ri≤0.2)
Qsens/lat=0Ri>0.2)
Ri=gTair(c)hz/Tairu2

式中:Qsens/latQ'sens/lat分别为稳定性修正后和修正前的感热和潜热通量;Ri为总体理查逊数;g为重力加速度(9.81 m⋅s-2);Tair(c)Tair分别为摄氏度气温(℃)和开尔文气温(K)。

COSIMA模型表层下雪冰模型部分是将冰面以下10 m范围内的距离以0.2 m等分,由于每一层的水热组成存在差异,从而使得每一层的物理性质(冰温、密度以及液体水含量)不尽相同16。基于野外实测数据,初始雪深设置为0.32 m。初始温度剖面利用1号冰川实测的表面温度和冰川底部温度(-5.2 ℃)通过线性插值获得。初始的雪坑密度剖面也是基于1号冰川密度经验值以插值的方式得到,其中表层雪和粒雪密度分别采用经验值280 kg⋅m-3和265 kg⋅m-3;雪坑以下冰川冰的密度采用经验值900 kg⋅m-3

采用COSIMA模型进行1号冰川物质平衡模拟研究具有如下两方面优势:第一,在各能量通量作用下,物质平衡以物质流的形式发生变化,COSIMA模型可以将各物质流作用下导致的逐小时物质平衡变化进行区分,实现物质平衡精细化研究;第二,与以往只考虑表层物质平衡相比,COSIMA模型可以较为系统的呈现物质平衡在表层下的变化过程,比如融水渗浸、再冻结以及密实化过程16,使得物质平衡的计算更为全面、合理。

下面详细介绍模型中需要的数据,输入数据包括入射短波辐射、云因子、气温(2 m)、降水量、风速(2 m)、相对湿度(2 m)及气压。并且在复杂冰面气象环境的影响下,AWS1在8月20日12:00至8月22日12:00内的数据缺失,为了获得连续的气象数据集,本文利用AWS2及AWS3相应时间段内的各气象要素来进行补充。

(1) 入射短波辐射

COSIMA模型中入射短波辐射由CNR4传感器直接测得。虽然受大气透明度及复杂地形的影响,入射短波辐射在空间上呈现差异性1133,但是,入射短波辐射在较小海拔间仍具有强相关性,因此,尝试建立了AWS1与AWS3入射短波辐射的拟合关系,两者拟合相关系数达0.95,使用AWS3的入射短波值来代替缺失时段。

(2) 云因子

云对辐射的吸收和反射程度不同,其中10%~20%辐射量被云吸收,同时36%~80%的辐射量被云反射34,所以考虑云对短波辐射的影响是十分必要的。由于缺乏云量数据的直接观测,选择利用式(12)计算得到云因子,其中SWTOA为大气顶层太阳辐照度,取其常数为1 368 W⋅m-2[35

N=1.3-1.4SWin/SWTOA

(3) 气温和降水量

气温和降水量作为气候变化最为直接的参数,其变化程度深刻影响着大气与冰川区之间的能量交换,进而影响冰川区物质的积累和消融。由于存在缺失数据,因此,利用AWS1、AWS2和AWS3的气温和降水数据,考虑海拔梯度-1.1 ℃⋅(100m)-1、1.6 mm⋅(100m)-1,填补AWS1气温和降水量空白值。

(4) 风速、相对湿度及气压

近地面湍流交换受大气和冰川表面湿度差及风速等因素的影响36,利用Hock等11、Braithwaite等37的方法插补AWS1风速,其理论依据是假设整个冰川面上风速相同的情况下建立气象站点之间的关系,其中,Hock等11通过比较Storglacia¨ren冰川上分布在不同海拔的风速,发现各气象站间散点分布明显,风速的高低与各站点所处位置有关,因此,无法得出一个一般性的方案来定量分析风速在冰川面上的空间变化。所以,在考虑地形对风速的影响下,建立布设在冰舌和平衡线处风速的正相关关系,获得了良好的分布式模拟结果,因此,利用6月24日至8月20日AWS1和AWS3风速,其相关系数为0.41,虽然两站点风速相关性较低,但也可用于插补AWS1站点缺失值。相对湿度的插补同样依据Hock等11描述的方法,研究指出相对湿度在气象站点间的变化相对较小,所以,可以假设两站点间相对湿度是相同的,基于这一假设可利用气象站点相对湿度拓展为整个冰川面上的相对湿度。并且实际情况中AWS1与AWS2站点之间距离较近,因此,可利用AWS2同时段相对湿度获得AWS1连续数据序列。为了精确地获得气压与高度的对应关系,利用等温大气压高方程插补缺失的气压值,式中是把大气当成干空气处理的,并且实际大气并非等温大气,重力加速度也随纬度和高度而有变化,但在大气底层(<100 km)重力加速度随高度的变化不大,通过计算标准大气中气压与高度的对应值,发现气压仍是满足随高度按指数递减的规律,且等温大气与实际大气非常接近38

PAWS1=PAWS3/exp(HAWS1-HAWS3)/8000(1+αt)

式中:PAWS1PAWS3分别为AWS1和AWS3处气压值;t为气层间平均气温;HAWS1HAWS3为AWS1和AWS3处海拔,分别为4 025 m和3 835 m;α为常数,约为0.036。

2.4 误差评估

根据Huintjes等16研究,COSIMA模拟物质平衡的日误差(Uday)评估采用如下方法。

Uday=16Xj,max-Xj,min/nj

式中:Xj,maxXj,min分别为利用6次测量间隔获得的降水梯度运行COSIMA得到的物质平衡最大值和最小值;nj 为测量间隔j的天数。利用降水梯度[(1.6±0.3)mm⋅(100m)-1]分别运行COSIMA,其中,1.6 mm⋅(100m)-1作为参照降水梯度。综合评估表明,消融期冰川物质平衡的日误差为±0.03 m w.e.。

3 结果与分析

3.1 消融期冰面气象状况

自动气象站观测结果显示,2018年消融期,1号冰川冰面气温波动存在显著的差异性。消融期初(4月29日至6月8日),气温波动幅度剧烈,日均值为-3.2 ℃;6月9日至8月28日,气温大于0 ℃;消融期末(8月19日至9月1日)气温持续递减[图2(a)]。整个消融期相对湿度均值约为69%,其变化态势与气温较为类似[图2(b)],其中,5月相对湿度低,且低于均值的比重为91%。整个消融期,风速在0.6~5.3 m⋅s-1内波动变化,均值为2.3 m⋅s-1图2(c)],以6月27日为分界,4月29日至6月27日风速相对较高,之后变化小。气压在消融期内的变化相对稳定,均值为625 hPa[图2(d)]。消融期降水量为692.2 mm,主要降水集中在6—8月[图2(e)]。

图2

图2   2018年消融期气象状况日均值

Fig.2   Daily mean values of meteorological condition during the ablation season in 2018: air temperature (a), relative humidity (b), wind speed (c), air pressure (d), precipitation (e)


3.2 COSIMA模型验证

利用消融期物质平衡模拟结果对比实测的物质平衡,从散点图结果可以看出[图3(a)],模拟的物质平衡值与实测的物质平衡值有非常好的一致性,相关系数达0.96,离散程度(RMSE)为±0.22 m w.e.。4月29日至8月17日物质平衡先缓慢减小,然后受强烈消融影响快速减小,8月17日后物质平衡变化逐渐减小。与实测值相比,7月20日与8月1日两次模拟物质平衡被低估,小于实测物质平衡,可能源于模型计算的湍流通量被低估,导致冰川消融减少[图3(b)]。尽管如此,整体来看模拟结果能很好的呈现消融期1号冰川的物质平衡变化,为进一步模拟其能量通量及物质平衡的变化特征提供了有效支撑。

图3

图3   模拟与实测累积物质平衡

Fig.3   Modelled and observed accumulative mass balances: scatter plot (a), change trend (b)


3.3 能量通量特征

消融期各能量项月尺度能量组成如表2所示。消融期冰川表面能量收入项中净短波辐射、感热通量所占比重的均值为84%、16%,所以,净短波辐射是冰川表面主要的能量来源,其次是感热通量。净短波辐射(SWnet)在月尺度上表现出先减小后增加的趋势,6月和8月分别出现最小值(+30.48 W⋅m-2)和最大值(+108.43 W⋅m-2),影响因素除太阳倾斜角外,实测入射短波辐射在8月显著增加是其主要原因,而6月出现最小值主要是源于频繁的降雪事件使得冰川表面的反照率增加,冰川表面接收的太阳辐射减小。在扎当冰川上的研究也得到了类似的结果39,说明消融期冰川表面接收的太阳辐射受表面反照率反馈机制的调控。感热通量在消融期内始终是正值(+10.89 W⋅m-2),表明冰川上空2 m处气温高于冰川表面温度,大气向冰川表面传输能量39。而净长波辐射、潜热通量、用于冰雪消融的能量及地热通量作为能量支出项所占比重依次为55%、7%、32%及6%,可以看出,能量的支出项主要是净长波辐射和用于冰雪消融的能量。消融期,净长波辐射月均值为-39.3 W⋅m-2,6—7月相对湿度较高(RH=69%),并且气温和模拟冰川表面温度等于0 ℃,导致净长波辐射(LWnet)在6—7月相对较高且存在少数正值(图4)。观测期间,气温较高,空气水汽压小于冰川表面水汽压,使得多数潜热通量在消融期内是负值,表明冰川表面以升华和蒸发为主,同时向大气输出热量23,潜热通量从负值转变为正值的天数中,地气温差和风速较小(图2),不利于湍流的发展,此时大气向冰川表面传输能量,水汽便会发生凝结或凝华40。地热通量(QG)是透射短波辐射(Qps)和热传导通量(QC)的和,均值为 -4.85 W⋅m-2,此外,受模型中冰川最底部冰温为常数的影响,极有可能高估了地热通量。当冰温达到0 ℃时,透射短波辐射将会用于冰下消融。此外,消融期由升华或蒸发消耗的能量比例Qlat/Qsens+Qlat在46%(5月)至5%(8月)之间,均值为30%。

表2   消融期乌鲁木齐河源1号冰川东支AWS1处月尺度能量组成变化及占比

Table 2  Monthly variation and contributions of energy components at the AWS1 site on the east branch of the Urumqi Glacier No.1 during the ablation season

月份收入项/(W⋅m-2支出项/(W⋅m-2
净短波辐射感热通量入射长波辐射净长波辐射潜热通量消融能量地热通量出射长波辐射
542.5013.72219.15-50.14-11.75-1.44-11.39-269.30
630.489.17269.75-18.98-2.73-16.02-1.15-288.80
747.808.73264.49-25.33-2.77-26.18-2.53-289.80
8108.4311.93222.49-62.76-0.67-49.56-4.35-285.20
平均值57.3010.89242.83-39.30-4.48-23.30-4.85-283.30
占比84%16%55%7%32%6%

新窗口打开| 下载CSV


图4

图4   消融期乌鲁木齐河源1号冰川东支日均能量组成状况

Fig.4   Daily energy components on the east branch of the Urumqi Glacier No.1 during the ablation season


1号冰川是典型的大陆型冰川,为了进一步研究该冰川能量通量的变化特征,将其与我国其他大陆型冰川已有研究开展对比分析(表3)。虽然不同研究所涉及的年份的气象状况各异、研究时段和模拟方法存在差异,但仍然能够从整体上揭示大陆型冰川能量通量的变化特征。与1号冰川早期能量平衡研究相似41-42,造成冰川消融的能量主要来源于净短波辐射,消融冰面由于其逆温层结,感热通量占比次之;冰川能量支出主要为净长波辐射、冰川消融耗热。由于消融冰面的蒸发受制于表面水汽压,使得冰面蒸发受到一定的抑制作用,潜热通量负值较大。与1号冰川早期能量平衡参数化方案相比,COSIMA更具系统性,且考虑冰面以下10 m范围内的渗浸短波辐射和热传导通量,更为细致地反映能量物质平衡过程机理。同一冰川不同海拔位置,海拔越高,净辐射越小41-46,净辐射值也存在明显的年际变化,同时,不同时期各能量通量的差异涉及冰面气象状况、参数化方案及研究时段等因素。此外,李晶等43指出表碛表面净辐射小于雪面净辐射,说明冰川下垫面差异对净辐射的影响也是不容忽视的。利用扎当冰川与其他区域冰川做比较,发现净辐射是扎当冰川最大的能量收入项(+39 W⋅m-2),占总的能量收入项的82%,同时考虑到扎当冰川的研究点海拔显著高于其他冰川,所以,不能将扎当冰川与其余冰川进行对比,但是净辐射在不同类型冰川上的比重差异是明显的47。相似海拔下冰川所处气象状况可直接影响湍流通量的大小,其中七一冰川(+14.2 W⋅m-2)及科契卡尔巴西冰川(+14.4 W⋅m-2)的感热通量大于乌鲁木齐河源1号冰川(+10.89 W⋅m-2)及老虎沟12号冰川(+5.7 W⋅m-2),表明气温越高,风速越大,感热通量越大434648,而潜热通量的低值说明潜热交换剧烈消耗的能量多,与其余冰川相比,消融期科契卡尔巴西冰川潜热交换剧烈(-23 W⋅m-2),说明科契卡尔巴西冰川表面饱和水汽压和空气水汽压差异明显,水汽压梯度最大,蒸发最为强烈43。消融能量也具有随海拔增加而减小的特征,但只涉及乌鲁木齐河源1号冰川以及老虎沟12号冰川。此外,基于扎当冰川和帕隆4号冰川及其他区域典型监测冰川能量通量的对比结果表明单条冰川所处气候环境对消融能量的影响的确存在,但难以定量化,而要准确衡量对消融能量的影响程度还有赖于从控制研究时段和采用相同模型等两方面进行着手研究47

表3   乌鲁木齐河源1号冰川与其他大陆型冰川能量组成对比

Table 3  Comparison of energy components of the Urumqi Glacier No.1 with other continental glaciers

冰川名称纬度/N经度/E海拔/m月份

净辐射

/(W⋅m-2

感热通量

/(W⋅m-2

潜热通量

/(W⋅m-2

地热通量

/(W⋅m-2

消融能量

/(W⋅m-2

文献来源
乌鲁木齐河源1号冰川43.1°86.8°4 0255—817.210.9-4.5-4.9-23.3本文
乌鲁木齐河源1号冰川43.1°86.8°3 910

6—8

(1986—1987年)

73.013.0-5.0-81.041
乌鲁木齐河源1号冰川43.1°86.8°3 910

7—8

(1986年)

54.511.6-4.0-62.142
科契卡尔巴西冰川42.0°80.0°4 2006—963.314.4-23.0-54.043
科契卡尔巴西冰川42.0°80.0°3 300697.78.6-53.444
老虎沟12号冰川39.4°96.4°5 0406—927.310.3-11.9-7.6-18.245
老虎沟12号冰川39.4°96.4°4 5506—980.85.7-5.8-80.846
扎当冰川30.4°90.6°5 6555—939.08.0-11.0-11.0-25.047
七一冰川39.5°97.7°4 4737—1063.314.2-6.1-15.5-55.848

新窗口打开| 下载CSV


3.4 物质平衡特征

冰川物质平衡是冰川积累和消融的代数和,取决于降水、消融、升华、蒸发以及再冻结过程。1号冰川累积物质平衡及其组成显示(表4),模拟的消融期累积物质平衡为-0.671 m w.e.,其中,表面消融是最主要的物质损失项(-0.742 m w.e.),其次是冰下消融(-0.114 m w.e.),升华/蒸发也会引起物质损失,但值相对较小(-0.022 m w.e.)。固态降水是主要的物质收入项(0.196 m w.e.),再冻结带来的物质收入仅为0.007 m w.e.。总的来说,几乎全部表层和表层下的消融量都随着冰川径流而损失掉。由于地表水汽压小于空气水汽压,加之风速较小,产生了0.005 m w.e.水汽凝结量。Weidemann等49指出积雪场的能量收支会受到风吹雪的影响,COSIMA模型中缺少风吹雪的参数化方案,因此会导致上述升华作用被高估,但消融期1号冰川消融强烈,冰川表面积雪性状会因气温和融水等因素发现改变,因此,只有当风速较大且频繁有降雪发生时升华作用则被认为高估。此外,Oerlemans等50在西格林兰冰盖及瑞士的冰川上的研究指出冰川重力流在消融冰川中非常常见,能够增强湍流交换,而在应用COSIMA模型时忽略了冰川重力流,因此可能会造成表面消融量的低估。虽然会存在上述问题,但总体来看,消融期1号冰川物质平衡的变化主要取决于表面消融和固态降水。

表4   乌鲁木齐河源1号冰川累积物质平衡组成及与其他大陆型冰川对比

Table 4  Accumulative mass balance components of the Urumqi Glacier No.1 and comparison with other continental glaciers

冰川名称纬度/N经度/E月份

累积物质平衡

/(m w.e.)

固态降水

/(m w.e.)

再冻结

/(m w.e.)

凝结

/(m w.e.)

表面消融

/(m w.e.)

冰下消融

/(m w.e.)

升华及蒸发

/(m w.e.)

文献来源
乌鲁木齐河源1号冰川43.1°86.8°5—8-0.6710.1960.0070.005-0.742-0.114-0.022本文
羌塘1号冰川33.0°88.7°5—9-0.1760.3800.135-0.532-0.07114
扎当冰川30.4°91.0°5—9-0.4310.4240.104-0.867-0.09239

新窗口打开| 下载CSV


为进一步研究乌鲁木齐河源1号冰川物质平衡组成,将其与其他区域监测的大陆型冰川进行比较。结果表明,与羌塘1号冰川、扎当冰川相比,乌鲁木齐河源1号冰川物质损失明显(表4)。夏季积累型冰川深受降水的季节变化及固态降水比重的影响51,而乌鲁木齐河源1号冰川固态降水显著小于其余两者,主要原因是羌塘1号冰川和扎当冰川分别在消融期深受西风和印度季风的影响,使得固态降水量增加。另一方面,尽管表面消融和升华/蒸发作用差异很小,但乌鲁木齐河源1号冰川再冻结作用显著小于其余两者。已知再冻结作用主要发生在消融期40,取决于冰/积雪温度及其物理性质52,并且冰/积雪温度越高,再冻结过程越弱16,所以,从定性的角度上可以得出乌鲁木齐河源1号冰川的冰/积雪温度高于其余两者,而再冻结过程的物理机制还有待进一步实测验证。

4 结论与展望

针对乌鲁木齐河源1号冰川东支,以气象数据驱动COSIMA模型,结合冰川物质平衡观测数据,实现了该冰川2018年消融期单点能量-物质平衡模拟。得到如下结论:

(1)物质平衡模拟值为(-0.67±0.03) m w.e.,模拟与实测物质平衡较为一致,相关系数达0.96,因此,COSIMA应用在乌鲁木齐河源1号冰川上能够较为准确地计算单点能量-物质变化。

(2)净短波辐射是主要的能量来源(84%),其次为感热通量(16%)。净长波辐射是主要的能量支出项(55%),其次分别是冰川消融耗热(32%)、潜热消耗(7%)及地热通量(6%)。受能量收支的影响,模拟的冰川物质平衡主要取决于表面消融和固态降水量,同时在模型自身局限性的影响下,模拟升华和表面消融量存在不可避免的误差。

(3)与其他大陆型冰川已有研究结果对比发现,不同海拔以及气候条件下能量通量的差异明显,同时,物质平衡对比结果表明乌鲁木齐河源1号冰川固态降水量和再冻结显著小于羌塘1号冰川和扎当冰川,推测与单条冰川所处大气环流密切有关。

鉴于此,为了更为深入地研究1号冰川能量-物质平衡特征及其对气候变化的响应,深刻认识气候机制作用下的冰川变化过程及机理,下一步将开展单点到分布式的能量-物质平衡模拟研究。第一,完善观测数据。需要在1号冰川消融区补充气象和四分量辐射数据的观测仪器,以便获得更高精度的分布式输入数据。第二,完善观测技术。除常规观测数据外,国际上对于冰川能量-物质平衡的模拟验证通常将遥感影像和延迟摄影技术结合,提取不同时间同一时刻雪线的瞬时变化,来获得研究期观测及模拟雪线的时空变化,借以验证模型模拟效果,从而有助于提高模型模拟精度。第三,修正输入数据以驱动模型,获得模拟结果。基于1号冰川DEM数据提取掩膜,利用WRF模式修正气象要素以获得不同海拔上气象要素梯度变化,通过双曲正切函数计算固态降水量比重,最后结合入射短波辐射数据驱动二维模型,获得1号冰川分布式能量-物质平衡模拟结果。

参考文献

IPCC. Climate change 2013: the physical science basis: contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change[M]. New YorkCambridge University Press2013.

[本文引用: 1]

Farinotti DLonguevergne LMoholdt Get al.

Substantial glacier mass loss in the Tien Shan over the past 50 years

[J]. Nature Geoscience, 201589): 716-722.

[本文引用: 2]

Wang PuyuLi ZhongqinLi Huilinet al.

Analysis of the relation between glacier volume change and area change in the Tianshan Mountains

[J]. Journal of Glaciology and Geocryology, 2017391): 9-15.

[本文引用: 1]

王璞玉李忠勤李慧林.

天山冰储量变化和面积变化关系分析研究

[J]. 冰川冻土, 2017391): 9-15.

[本文引用: 1]

Ye WanhuaWang FeitengLi Zhongqinet al.

Temporal and spatial distributions of the equilibrium line altitudes of the monitoring glaciers in High Asia

[J]. Journal of Glaciology and Geocryology, 2016386): 1459-1469.

[本文引用: 1]

叶万花王飞腾李忠勤.

高亚洲定位监测冰川平衡线高度时空分布特征

[J]. 冰川冻土, 2016386): 1459-1469.

[本文引用: 1]

Liu ChaohaiXie ZichuWang Chunzu.

A research on the mass balance processes of Glacier No

.1 at the headwaters of the Urumqi River, Tianshan Mountains[J]. Journal of Glaciology and Geocryology, 1997191): 17-24.

[本文引用: 1]

刘潮海谢自楚王纯足.

天山乌鲁木齐河源1号冰川物质平衡过程研究

[J]. 冰川冻土, 1997191): 17-24.

[本文引用: 1]

Neckel NKropáček JBolch Tet al.

Glacier mass changes on the Tibetan Plateau

2003-2009 derived from ICESat laser altimetry measurements[J/OL]. Environmental Research Letters, 201491) [2019-04-30]. .

[本文引用: 1]

Rasmussen L AKohler J.

Mass balance of three Svalbard glaciers reconstructed back to 1948

[J]. Polar Research, 2010262): 168-174.

[本文引用: 2]

Hock R.

Temperature index melt modeling in mountain areas

[J]. Journal of Hydrology, 20032821): 104-115.

[本文引用: 1]

Huintjes ELi HuilinSauter Tet al.

Degree-day modelling of the surface mass balance of Urumqi Glacier No.1, Tian Shan, China

[J]. The Cryosphere, 20104207-232.

[本文引用: 1]

Gabbi JCarenzo MPellicciotti Fet al.

A comparison of empirical and physically based glacier surface melt models for long-term simulations of glacier response

[J]. Journal of Glaciology, 2014601140-1154.

[本文引用: 1]

Hock RHolmgren B.

A distributed surface energy balance model for complex topography and its application to Storglaciären, Sweden

[J]. Journal of Glaciology, 200551172): 25-36.

[本文引用: 6]

Plach ANisancioglu K HLe Clec’h Set al.

Eemian Greenland surface mass balance strongly sensitive to SMB model choice

[J]. Climate of the Past, 20181410): 1463-1485.

[本文引用: 1]

Yang WeiYao TandongGuo Xiaofenget al.

Mass balance of a maritime glacier on the southeast Tibetan Plateau and its climatic sensitivity

[J]. Journal of Geophysical Research: Atmospheres, 201311817): 9579-9594.

[本文引用: 1]

Li ShenghaiYao TandongYang Weiet al.

Glacier energy and mass balance in the inland Tibetan Plateau: seasonal and interannual variability in relation to atmospheric changes

[J]. Journal of Geophysical Research: Atmospheres, 20181236390-6409.

[本文引用: 2]

Zhu MeilinYao TandongYang Weiet al.

Differences in mass balance behavior for three glaciers from different climatic regions on the Tibetan Plateau

[J]. Climate Dynamics, 2017509/10): 3457-3484.

[本文引用: 1]

Huintjes ESauter TSchröter Bet al.

Evaluation of a coupled snow and energy balance model for Zhadang Glacier, Tibetan Plateau, using glaciological measurements and time-lapse photography

[J]. Arctic, Antarctic, and Alpine Research, 2015473): 573-590.

[本文引用: 5]

Huintjes ENeckel NHochschild Vet al.

Surface energy and mass balance at Purogangri ice cap, central Tibetan Plateau, 2001-2011

[J]. Journal of Glaciology, 201561230): 1048-1061.

[本文引用: 1]

Xu ChunhaiLi ZhongqinWang Feitenget al.

Using an ultra-long-range terrestrial laser scanner to monitor the net mass balance of Urumqi Glacier No.1, eastern Tien Shan, China, at the monthly scale

[J]. Journal of Glaciology, 201763241): 792-802.

[本文引用: 1]

Li ZhongqinShen YongpingWang Feitenget al.

Response of glacier melting to climate change: take Urumqi Glacier No.1 as an example

[J]. Journal of Glaciology and Geocryology, 2007293): 333-342.

[本文引用: 1]

李忠勤沈永平王飞腾.

冰川消融对气候变化的响应: 以乌鲁木齐河源1号冰川为例

[J]. 冰川冻土, 2007293): 333-342.

[本文引用: 1]

Wang PuyuLi ZhongqinLi Huilinet al.

Analyses of recent observations of Urumqi Glacier No.1, Chinese Tianshan Mountains

[J]. Environmental Earth Sciences, 2016758): 1-11.

[本文引用: 1]

Yue XiaoyingZhao JunLi Zhongqinet al.

Spatial and temporal variations of the surface albedo and other factors influencing Urumqi Glacier No.1 in Tien Shan, China

[J]. Journal of Glaciology, 201763241): 899-911.

[本文引用: 1]

van den Broeke Mvan As DReijmer Cet al.

Assessing and improving the quality of unattended radiation observations in Antarctica

[J]. Journal of Atmospheric and Oceanic Technology, 2004219): 1417-1431.

[本文引用: 1]

Mölg TMaussion FYang Weiet al.

The footprint of Asian monsoon dynamics in the mass and energy balance of a Tibetan glacier

[J]. The Cryosphere, 201261445-1461.

[本文引用: 2]

He XiaoboYe BaishengDing Yongjian.

Bias correction for precipitation measurement in Tanggula Mountain, Tibetan Plateau

[J]. Advances in Water Science, 2009203): 403-408.

[本文引用: 1]

何晓波叶柏生丁永建.

青藏高原唐古拉山区降水观测误差修正分析

[J]. 水科学进展, 2009203): 403-408.

[本文引用: 1]

Wang PuyuLi ZhongqinLi Huilinet al.

Characteristics of a partially debris-covered glacier and its response to atmospheric warming in Mt. Tomor, Tien Shan, China

[J]. Global and Planetary Change, 201715911-24.

[本文引用: 1]

Thibert EVincent CBlanc Ret al.

Instruments and methods glaciological and volumetric mass-balance measurements: error analysis over 51 years for Glacier de Sarennes, French Alps

[J]. Journal of Glaciology, 200854186): 522-532.

[本文引用: 2]

Zemp MJansson PHolmlund Pet al.

Reanalysis of multi-temporal aerial images of Storglaciären, Sweden (1959-99): Part 2: comparison of glaciological and volumetric mass balances

[J]. The Cryosphere, 2010204): 345-357.

[本文引用: 1]

Oerlemans J. Glaciers and climate change[M]. Lisse, the NetherlandsSwets and Zeitlinger2001148.

[本文引用: 1]

Brock B WWillis I CMartin M J.

Measurement and parameterization of aerodynamic roughness length variations at Haut Glacier d’Arolla, Switzerland

[J]. Journal of Glaciology, 200652281-297.

[本文引用: 1]

Cullen N JMölg TKaser Get al.

Energy-balance model validation on the top of Kilimanjaro, Tanzania, using eddy covariance data

[J]. Annals of Glaciology, 200746227-233.

Gromke CManes CWalter Bet al.

Aerodynamic roughness length of fresh snow

[J]. Boundary Layer Meteorology, 201114121-34.

[本文引用: 1]

Braithwaite R J.

Aerodynamic stability and turbulent sensible heat flux over a melting ice surface, the Greenland Ice Sheet

[J]. Journal of Glaciology, 199541562-571.

[本文引用: 1]

Xin HuijuanHe YuanqingNiu Hewenet al.

Near-surface meteorological characteristics on the Baishui Glacier No.1, Mt. Yulong

[J]. Journal of Glaciology and Geocryology, 2018404): 676-684.

[本文引用: 1]

辛惠娟何元庆牛贺文.

玉龙雪山白水1号冰川近地层气象要素变化特征

[J]. 冰川冻土, 2018404): 676-684.

[本文引用: 1]

Cess R DZhang M HMinnis Pet al.

Absorption of solar radiation by clouds: observations versus models

[J]. Science, 19952675197): 496-499.

[本文引用: 1]

Favier VWagnon PChazarin J Pet al.

One year measurements of surface heat budget on the ablation zone of Antizana Glacier

15, Ecuadorian Andes[J/OL]. Journal of Geophysical Research, 2004, 109(D18) [2019-04-30]. .

[本文引用: 1]

Sun WeijunQin XiangRen Jiawenet al.

Surface energy balance in the accumulation zone of the Laohugou Glacier No.12 in the Qilian Mountains during ablation period

[J]. Journal of Glaciology and Geocryology, 2011331): 38-46.

[本文引用: 1]

孙维君秦翔任贾文.

祁连山老虎沟12号冰川积累区消融期能量平衡特征

[J]. 冰川冻土, 2011331): 38-46.

[本文引用: 1]

Braithwaite R JZhang Yu.

Sensitivity of mass balance of five Swiss glaciers to temperature changes assessed by tuning a degree-day model

[J]. Journal of Glaciology, 200046152): 7-14.

[本文引用: 1]

Zhou ShuzhenZhang RuyiZhang Chao. Meteorology and climatology[M]. BeijingHigher Education Press199784-86.

[本文引用: 1]

周淑贞张如一张超. 气象学与气候学[M]. 北京高等教育出版社199784-86.

[本文引用: 1]

Zhang GuoshuaiKang ShichangFujita K.

Energy and mass balance of Zhadang Glacier surface, central Tibetan Plateau

[J]. Journal of Glaciology, 201359213): 137-148.

[本文引用: 3]

Fang XiaoyuLi ZhongqinWuennemann Bet al.

Physical energy-balance and statistical glacier melting models comparison and testing for Shiyi Glacier, Heihe River basin, Qilian Mountains, China

[J]. Journal of Glaciology and Geocryology, 2015372): 336-350.

[本文引用: 2]

方潇雨李忠勤Wuennemann B.

冰川物质平衡模式及其对比研究: 以祁连山黑河流域十一冰川研究为例

[J]. 冰川冻土, 2015372): 336-350.

[本文引用: 2]

Kang Ersi.

A parameterized energy balance model of glacier melting on the Tianshan Mountains

[J]. Acta Geographica Sinica, 1994495): 467-476.

[本文引用: 3]

康尔泗.

天山冰川消融参数化能量平衡模型

[J]. 地理学报, 1994495): 467-476.

[本文引用: 3]

Kang ErsiYang DaqingZhang Yinsheng.

An experiment study of the water and heat balance in the source area of the Urumqi River in the Tian Shan Mountains

[J]. Annals of Glaciology, 19921655-66.

[本文引用: 2]

Li JingLiu ShiyinZhang Yong.

Snow surface energy balance over the ablation period on the Keqicar Baxi Glacier in the Tianshan Mountains

[J]. Journal of Glaciology and Geocryology, 2007293): 366-374.

[本文引用: 4]

李晶刘时银张勇.

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

[J]. 冰川冻土, 2007293): 366-374.

[本文引用: 4]

Han HaidongDing YongjianLiu Shiyin.

Estimation and analysis of heat balance parameters in the ablation season of debris-covered Kerqikaer Glacier, Tianshan Mountains

[J]. Journal of Glaciology and Geocryology, 2005271): 88-94.

[本文引用: 1]

韩海东丁永建刘时银.

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

[J]. 冰川冻土, 2005271): 88-94.

[本文引用: 1]

Sun WeijunQin XiangRen Jiawenet al.

The surface energy budget in the accumulation zone of the Laohugou Glacier No.12 in the western Qilian Mountains, China, in summer 2009

[J]. Arctic, Antarctic, and Alpine Research, 2012443): 296-305.

[本文引用: 1]

Sun WeijunQin XiangWang Yetanget al.

The response of surface mass and energy balance of a continental glacier to climate variability, western Qilian Mountains, China

[J]. Climate Dynamics, 2018509/10): 3557-3570.

[本文引用: 3]

Zhu MeilinYao TandongYang Weiet al.

Energy- and mass-balance comparison between Zhadang and Parlung No.4 Glaciers on the Tibetan Plateau

[J]. Journal of Glaciology, 201561227): 595-607.

[本文引用: 3]

Jiang XiWang NinglianHe Jianqiaoet al.

A distributed surface energy and mass balance model and its application to a mountain glacier

[J]. Chinese Science Bulletin, 20105520): 2079-2087.

[本文引用: 2]

蒋熹王宁练贺建桥.

山地冰川表面分布式能量-物质平衡模型及其应用

[J]. 科学通报, 20105518): 1757-1765.

[本文引用: 2]

Weidemann S SSauter TMalz Pet al.

Glacier mass changes of lake-terminating Grey and Tyndall Glaciers at the Southern Patagonia Icefield derived from geodetic observations and energy and mass balance modeling

[J]. Frontiers of Earth Science, 2018681): 1-16.

[本文引用: 1]

Oerlemans JGrisogono B.

Glacier winds and parametrizations of the related surface heat fluxes

[J]. Tellus, 2002545): 440-452.

[本文引用: 1]

Kang ShichangChen FengGao Tanguanget al.

Early onset of rainy season suppresses glacier melt: a case study on Zhadang Glacier, Tibetan Plateau

[J]. Journal of Glaciology, 200955192): 755-758.

[本文引用: 1]

Illangasekare T HRodney J W JMeier M Fet al.

Modeling of meltwater infiltration in subfreezing snow

[J]. Water Resources Research, 1990265): 1001-1012.

[本文引用: 1]

/