高亚洲冰川区度日因子空间分布数据集

张勇，刘时银，王欣

1. 湖南科技大学资源环境与安全工程学院，湖南湘潭 411201

2. 云南大学国际河流与生态安全研究院，昆明 650500

数据库（集）基本信息简介

引 言

高亚洲泛指包括青藏高原、帕米尔高原、天山和阿尔泰山等构成的海拔超过2000 m的地区，是除极地之外最大的冰川发育区，冰川总面积约占全球冰川总面积的13.6%[1-2]。冰川是高亚洲地区极其重要的固体水资源，其融水是周边地区数条大江大河重要的补给来源[3-5]，尤其是我国西北干旱区的水资源在很大程度上依赖于冰川融水[6-8]。在全球变暖背景下，高亚洲地区冰川呈现不同程度的退缩与强烈的物质亏损[2,5,9-10]，其融水在不同时空尺度上的变化对区域水资源形成与变化、冰湖溃决洪水、冰川洪水及其诱发泥石流等灾害的影响十分显著[6-7,11-14]，尤其影响下游地区数亿人口的水供给和食品安全[3]。因此，冰川融水模拟是评估高亚洲地区水资源变化及影响的重要组成部分之一，其研究已成为社会各界关注的热点问题[7,15]。

在冰川融水计算方法中，度日模型是最为简单、应用最为广泛的方法。该方法是基于冰川消融与气温之间的线性关系建立的，能够给出类似能量平衡模型的输出结果[16-18]。在度日模型中，度日因子是其重要的参数，是冰雪表面能量传递与转化过程的简化描述[17-19]，其空间变化特征对于冰川消融模拟的精度有较大影响[16-18]。在冰川尺度上，度日因子可通过冰川区消融与气温观测计算获取；在区域/全球尺度上，度日因子具有较大的时空变异性[17-18]。在众多应用度日模型进行区域/全球冰川模拟的研究中，基于典型冰川物质平衡观测数据[20]对冰川冰和雪的度日因子进行率定，进而推广应用于每条冰川[21-26]。目前，全球开展长期物质平衡观测的冰川数量有限，而高亚洲地区具有长期物质平衡观测的冰川数量更少[20]。基于高亚洲现有冰川观测数据率定度日因子，进而开展整个区域冰川融水及其影响的评估，可能会引起一定的误差。为此，本文收集了过去几十年来高亚洲地区不同地区40余条冰川的短期考察和观测数据（图1），获取了每条观测冰川的度日因子值，包括冰川冰与雪的度日因子。以此数据为基础，构建了冰川冰与雪度日因子的转换公式，获取了整个高亚洲冰川区度日因子的空间分布格网数据。本数据精度已在典型区域的冰川物质平衡模拟应用中得到了验证[27-28]，获取数据可靠，可为进一步研究高亚洲区域冰川融水及其水资源效应和灾害效应提供必要的数据支撑。

图1 高亚洲地区观测冰川分布图

1 数据采集和处理方法

1.1 数据源

本文采用的数据包括冰川区消融和气温观测数据、0.5°格网气温和降水数据以及冰川编目数据。其中，冰川区消融和气温数据来源于过去几十年来高亚洲不同地区40余条冰川不同时期的考察和观测数据[17,20,29-31]（图1和表1）。这些冰川的消融数据主要是通过花杆实测获取的[17,20]，大部分冰川消融观测时长超过1个月。1951-2011年期间的0.5°格网日气温和降水数据[32]是选取了距离冰川最近的网格数据。高亚洲冰川信息来自于Randolph冰川编目[1]，主要包括冰川末端海拔高度和冰川经纬度。

表1 高亚洲地区观测冰川信息[17,20,29-31]

3 2 D u n a g i r i 3 0.5 5 7 9.9 0 1 9 8 6-1 9 9 0年3 3 S h a u n e G a r a n g 3 1.2 8 7 8.3 3 1 9 8 2-1 9 9 0年3 4 G o r G a r a n g 3 1.6 2 7 8.8 2 1 9 7 7-1 9 8 3年3 5 N e h N a r 3 4.2 7 7 5.8 7 1 9 8 0-1 9 8 4年3 6 L e v i y A k t r u 5 0.0 8 8 7.7 3 1 9 7 7-2 0 1 1年3 7 M a l i y A k t r u 5 0.0 8 8 7.7 5 1 9 6 2-2 0 1 1年3 8 P r a v i y A k t r u 5 0.0 8 8 7.7 3 1 9 8 0-1 9 9 0年3 9 V o d o p a d n i y 5 0.1 0 8 7.7 0 1 9 7 7-2 0 1 1年4 0 K o s m o d e m a y a 4 3.0 0 7 7.1 0 1 9 6 5-1 9 9 0年4 1 P a r t i z a n 4 3.0 0 7 7.1 0 1 9 6 5-1 9 9 0年4 2 V i s y a c h i e 1-2 4 3.0 0 7 7.1 0 1 9 6 5-1 9 9 0年4 3 O r d z h o n i k i d z e 4 3.0 0 7 7.1 0 1 9 6 5-1 9 9 0年

1.2 数据处理步骤

本文基于收集的冰川区消融和气温观测数据，采用张勇等[16]方法计算了冰川冰与雪的度日因子值，其计算公式如下：

式中，DDF为冰川冰或雪的度日因子（mm d-1℃-1），M是某时段内冰川冰或雪的消融水当量（mm w.e.），PDD是同一时段内的正积温（℃），一般由下式计算获取：

式中，Tt为某天（t）的日平均气温，Ht为逻辑变量。当Tt≥0℃时，Ht= 1.0，当Tt＜0℃时，Ht= 0.0。大部分冰川存在同期观测的气温数据，对于仅有消融资料而没有气温资料的观测冰川区，使用距离该冰川区最近的格网气温数据。

基于公式（1）和（2）获取了高亚洲不同地区40条冰川的冰川冰与雪度日因子值。同时基于格网日气温和降水数据分别计算了 40条冰川区的年平均气温和年平均降水量。基于计算的度日因子值、冰川海拔高度、冰川经纬度、冰川区年平均气温和降水数据，采用多元线性回归方法，建立了冰川冰与雪度日因子的转换公式（表2），其相关系数（r）分别为0.62和0.92（显著性水平p＜0.001）。进而基于冰川海拔高度、冰川经纬度、冰川区年平均气温和降水数据，应用表1中的转换公式获取了高亚洲地区每一条冰川的度日因子值。

表2 冰川冰与雪的度日因子转换公式

2 数据样本描述

本数据集的数据存储格式为 GeoTIFF格式，数据类型为 32位浮点型，命名为Spatial_distribution_of_DDFice.tif和 Spatial_distribution_of_DDFsnow.tif，单位是 mm d-1℃-1。样本展示如图2和图3。

图2 高亚洲地区冰川冰度日因子空间分布图

图3 高亚洲冰川区雪度日因子空间分布图

高亚洲冰川冰的度日因子值介于3.0-15.0 mm d-1℃-1之间，而雪的度日因子值与冰川冰相比较小，其介于1.5-10.0 mm d-1℃-1之间。青藏高原东南部和喜马拉雅山东段地区冰川冰的度日因子值相对其他地区较大，而喀喇昆仑山、昆仑山、帕米尔等地区的度日因子值相对较小（图2）。对于雪的度日因子来说，青藏高原中部、喜马拉雅山区相对较大，而天山东部和青藏高原东部地区相对较小（图3）。度日因子这一显著的区域特征，可能与由青藏高原及其周边地区独特的气候和热量条件密切相关。对于青藏高原来说，冰川的气候环境从西北向东南逐渐由干冷向温湿转变，同时，热量平衡组成类型则由蒸发型向凝结-蒸发型转变[33]。不同地区典型冰川的能量平衡研究亦证明了这一变化趋势[34]。在干冷的气候条件下，冰面蒸发（主要是升华作用）耗热较大，抑制了冰川消融，度日因子较小；而在暖湿的气候条件下，用于冰川消融的能量在热量平衡中所占比例较大，促进了冰川消融，度日因子较大[16]。

3 数据质量控制和评估

本文基于冰川编目数据和格网气象数据，应用度日因子转换公式（表 2）获取了高亚洲冰川区冰川冰和雪的度日因子值（图2和图3），其平均值分别为6.6±2.3和5.2±2.7 mm d-1℃-1（±表示置信水平95%的误差区间）。同时，通过高亚洲北部阿尔泰山地区（图1）不同典型冰川物质平衡观测数据对获取的度日因子的精度进行了综合评估[27-28]。基于度日因子转换公式，获取了阿尔泰山冰川区的度日因子值。该地区冰川冰的度日因子介于6.0-11.0 mm d-1℃-1之间，平均值为8.6 mm d-1℃-1；雪的度日因子介于1.5-6.1 mm d-1℃-1之间，平均值为4.1 mm d-1℃-1。

图4 高亚洲北部阿尔泰山不同冰川物质变化观测值与模拟值对比

从图4a可以看出，基于上述度日因子值模拟的不同典型冰川物质平衡值与观测值较为接近。总体上，不同冰川物质平衡观测值与模拟值的相关系数和平均绝对误差分别为 0.77（显著性水平 p＜0.001）和0.21 mw.e.。从该地区典型冰川流域1977-2005年间物质平衡变化序列来看（图4b），流域物质平衡模拟值变化趋势与观测值较为一致，其平均值分别为-0.16 mw.e.和-0.12 mw.e.。此外，基于遥感数据，应用大地测量法获取了阿尔泰地区14条典型冰川1999-2005年的物质变化速率[35]，同时应用上述方法获取的度日因子值模拟了这14条冰川同一时段内的物质变化速率。比较发现（图4c），这两种不同方法获取的结果较为一致。总体上，大地测量法获取的结果与模型模拟值的相关系数和平均绝对误差分别为0.83（显著性水平p＜0.001）和0.14 mw.e.。由此可见，应用转换公式获取的度日因子能够较理想地模拟出冰川消融和积累的变化过程。

4 数据价值

全球气候变暖所导致的高亚洲地区冰川加速消融退缩，不仅引起了区域水资源变化[3,6,13-14]，同时引起了冰湖溃决洪水和冰川泥石流等灾害发生频率的加剧和影响程度加大[7,11-13]。开展高亚洲地区冰川消融过程及其对水资源与灾害的影响系统评估，不仅对建立高亚洲区域气候变化与冰川变化、水资源变化之间的关系研究具有重要科学意义，而且对我国推动的“一带一路”计划陆上丝绸之路关键区的经济建设与可持续发展具有重大现实意义[15,36]。本数据集是度日模型的重要参数，该模型在冰川物质平衡研究中应用广泛，且模拟冰川消融过程能够给出类似能量平衡模型的输出结果[16-18]。然而，对于高亚洲地区广泛分布的冰川来说，有长期观测的冰川数量相对较少，无法为度日模型开展区域冰川消融模拟提供必要的参数。本数据集是在不同地区冰川区观测数据的基础上获取的，能够反映出高亚洲冰川区度日因子的空间分布特征，且参数精度在典型区域进行了综合评估。因此，本数据集可为高亚洲区域冰川消融模拟及其对区域水资源与相关冰川灾害的研究奠定参数基础。

5 数据使用方法和建议

高亚洲冰川度日因子空间分布数据集所有数据存储格式均为GeoTIFF格式。ArcGIS、SuperMap、ENVI、ERDAS等常用GIS与遥感软件均可支持该数据的读取和操作。度日因子值以像元值表示，单位是mm d-1℃-1，空间分辨率为0.5°。本数据集可靠性高，代表性强，可作为高亚洲区域冰川消融模拟的重要输入数据，为进一步研究区域冰川变化对气候的响应及其影响等研究奠定基础。