1980～2014年岗日嘎布地区冰川高程变化数据集

高亚洲冰、雪和环境专题 I 区论文（评审中） • 版本 ZH4

A dataset of glacier elevation changes in the Kangri Karpo Mountains during 1980 – 2014

吴坤鹏，刘时银，蒋宗立，魏俊锋

Wu Kunpeng, Liu Shiyin, Jiang Zongli, Junfeng Wei

DOI: 10.11922/csdata.2018.0011.zh
: 10.11922/sciencedb.574

： 2018 - 04 - 18

： 2018 - 04 - 28

1107 17 0

摘要&关键词

摘要：由于印度洋季风的影响，位于藏东南的岗日嘎布是青藏高原最湿润的地区，海洋性冰川在该地区集中发育。该地区的冰川物质损失对海平面上升、调节河川径流及冰湖溃决灾害等有重要影响。本文基于1980年航测地形图、2000年2月11–22日SRTM数字高程模型（DEM）和2014年的X波段TerraSAR/TanDEM雷达影像得到岗日嘎布地区冰川高程变化。大致空间范围在29°N～30°N、96°E～98°E内，包括帕隆藏布以南、贡日嘎布曲以北区域，覆盖面积达3600平方公里。数据集内包含两类文件：1）基于TerraSAR/TanDEM雷达影像获取的2014年数字高程模型数据；2）基于地形图、SRTM和TerraSAR/TanDEM雷达影像获取的1980～2014年冰川高程变化栅格数据。为克服岗日嘎布地区多云雨对光学遥感影像的影响，对大地测量法中数字高程模型及冰川高程变化的提取采用合成孔径雷达差分干涉测量（DInSAR）。与ICESat-GLAS数据相比较，2014年数字高程模型的垂直精度达到2.65米。利用DInSAR获取的冰川高程变化转化为物质平衡，与野外测量获取的冰川物质平衡有较好的一致性。本数据集为目前岗日嘎布地区最新的高程变化数据，有效减少了云、季节性冰雪的影响，可作为基础数据集，应用到各种冰川研究中。

关键词：数字高程模型；冰川高程变化；岗日嘎布；DInSAR；TerraSAR/TanDEM

Abstract & Keywords

Abstract: Due to the influence of the Indian monsoon, the Kangri Karpo Mountains, located in southeast Tibet, is the most humid and maritime glacier-concentrated region on the Tibetan Plateau. Glacier mass loss in Kangri Karpo importantly contributes to the rise of global mean sea level, the adjustment of run-off distribution, and the increase of glacial lake outburst floods (GLOFs). Based on topographic maps (1980), SRTM (2000), and TerraSAR-X/TanDEM-X (2014), the study obtained the latest glacier elevation changes on the Kangri Karpo Mountains. Within a spatial extent of 29°N ~ 30°N, 96°E ~ 98°E totaling an area of 3600 km², the Mountains cover the region south of the Purlung Tsangpo River and north of the Gongri–Gabo River. The dataset consists of two subsets: ① 2014 digital elevation model (DEM) data obtained from TerraSAR-X/TanDEM-X, and ② raster data of glacier elevation changes during 1980 – 2014 based on topographic maps, SRTM and TerraSAR-X/TanDEM-X acquisitions. To mitigate influences of cloud cover on optical images, Differential Interferometric Synthetic Aperture Radar (DInSAR) was used to extract DEM data and acquire glacier elevation changes. Compared with ICESat-GLAS, the vertical precision of TerraSAR/TanDEM DEM reached 2.65 m. Glacier mass balance resulting from glacier elevation changes showed a high consistency with field measurements. The dataset records the latest glacier elevation changes in the Kangri Karpo Mountains, with minimum influences of cloud and seasonal snow cover on the glacier maps. It can be used as a fundamental dataset in future glacier-related studies.

Keywords:

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

数据库（集）名称	1980～2014年岗日嘎布地区冰川高程变化数据集
数据作者	吴坤鹏，刘时银，蒋宗立，魏俊锋
数据通信作者	吴坤鹏（wukunpeng2008@lzb.ac.cn），刘时银（liusy@ lzb.ac.cn）
数据时间范围	1980～2014年
地理区域	地理范围为北纬28°54′16″~ 29°31′42″，东经96°22′59″~ 97°8′23″，包括帕隆藏布以南、贡日嘎布曲以北区域。
空间分辨率	10米
数据量	142 MB
数据格式	.tif (Geo TIFF, 32bit float)（压缩为.rar格式）
数据服务系统网址	http://www.sciencedb.cn/dataSet/handle/574
基金项目	科技部专项（2013FY111400），国家自然科学基金项目（41190084，41471067，41671057，41671075），云南大学人才引进项目，国际合作项目（131C11KYSB20160061-4）。
数据库（集）组成	数据集由2部分数据组成，其一为2014年数字高程模型数据，其二是1980～2014年冰川高程变化栅格数据集。2个数据文件分别为： 1. 数字高程模型文件夹中是2014年数字高程模型数据； 2. 冰川表面高程变化文件夹中是1980～2014年冰川高程变化栅格数据集。

Dataset Profile

Title	A dataset of glacier elevation changes in the Kangri Karpo Mountains during 1980 – 2014
Data authors	Wu Kunpeng, Liu Shiyin, Jiang Zongl, Wei Junfeng
Data corresponding authors	Wu Kunpeng (wukunpeng2008@lzb.ac.cn), Liu Shiyin (liusy@ lzb.ac.cn)
Time range	1980～2014
Geographical scope	A spatial extent of 28°54′16″N ~ 29°31′42″N, 96°22′59″E ~ 97°8′23″E, which covers the region south of the Purlung Tsangpo River and north of the Gongri–Gabo River.
Spatial resolution	10 m
Data volume	142 MB
Data format	.tif (Geo TIFF, 32 bit float)
Data service system	http://www.sciencedb.cn/dataSet/handle/574
Sources of funding	Fundamental Programme of the Ministry of Science and Technology of China (MOST) (Grant No. 2013FY111400), National Natural Science Foundation of China (Grant Nos. 41190084, 41471067, 41671057 and 41671075) and Talent Introduction Program of Yunnan University , International Collaboration Project (131C11KYSB20160061-4).
Dataset composition	The dataset consists of two subsets: 2014 digital elevation model (DEM) data, and raster data of glacier elevation changes during 1980 – 2014.

引言

冰川是冰冻圈系统的重要组成部分^[1]，对气候变化具有强烈的指示作用^[2]。在当今全球气候变暖的背景下，冰川普遍退缩^{[3, 4]}。这对区域水资源有重要影响^{[5, 6]}，并引起不同程度的冰川灾害^{[7, 8]}。随着中国第二次冰川编目的完成^{[9, 10]}，高亚洲地区的冰川面积、长度变化研究成果涌现^[11-15]，但仅从面积变化难以理解冰川变化差异。冰川表面高程变化及物质平衡^[16-24]，日益成为当前冰川变化研究领域关注的热点。遥感技术具有方便快捷、空间覆盖范围大等特点，已逐渐取代了GPS等传统观测手段对冰川等恶劣气候环境区域的地面监测任务，如基于光学遥感立体像对提取DEM、激光雷达测高等方法获取冰川表面高程变化^{[16-21, 23, 24]}。基于光学遥感的冰川物质平衡研究，关键问题是获取无云、无季节性雪盖的影像^[25]，但在大量表碛覆盖、云活动频繁的高山地区，获取冰川物质平衡仍是一个有挑战性的问题。特别是青藏高原东南部岗日嘎布地区，受印度洋季风影响的湿润气候，冰川消融期有频繁的云雨活动，加上季节性积雪和表碛覆盖，获得高质量的光学遥感影像较为困难。因此，对岗日嘎布地区的冰川物质平衡研究，不同数据及方法得到不同的物质平衡水平，甚至相同数据得到的结果相差也较大^{[17, 19, 26, 27]}。不同物质平衡水平，会给后续冰川变化研究、水文建模等带来很大的不确定性。

随着遥感技术的发展，合成孔径雷达（SAR）技术已经具备了多波段、多极化、多模式等特点，能够全天候、不受天气影响地获取较高时空分辨率的地面雷达反射数据^[28]，十分有利于藏东南地区开展大范围的冰川表面高程变化监测。雷达干涉技术包括重复观测交轨干涉和双基站Bistatic交轨干涉，其中双基站Bistatic交轨干涉图相位包含平地相位与地形相位，常用于冰川高程变化研究，如SRTM DEM和TanDEM影像。基于TerraSAR-X/TanDEM-X遥感影像，对青藏高原及其周边地区的冰川表面高程变化均取得了较好的研究成果^[29-32]，可见该方法十分可靠。本文获取了2014年2月18日和2014年3月13日覆盖岗日嘎布地区的TerraSAR-X/TanDEM-X数据，利用合成孔径雷达差分干涉测量（DInSAR），对研究区内冰川进行了表面高程变化分析，得到了岗日嘎布地区1980～2014年冰川表面高程变化数据。

1 数据采集和处理方法

1.1 数据源

岗日嘎布地区冰川表面高程变化数据集主要基于4种数据：基于1980年航测地形图得到的数字高程模型（TOPO DEM）、2000年2月的数字高程模型（SRTM DEM）、基于双基站模式的TerraSAR-X/TanDEM-X数据和中国冰川编目数据（表1，图1）。

其中，共选取了12幅1:50000的航测地形图，其航摄时间为1980年10月，采用北京1954坐标系和黄海1985高程系。首先将原始地形图按照300 dpi在大幅面扫描仪上扫描，将扫描的地形图进行几何纠正，保证纠正精度在半个像元以内。然后对地形图的等高线进行矢量化，通过生成不规则三角网（TIN）模型，转化为DEM栅格文件，通过七参数法将北京坐标系/黄海高程系转换到WGS1984/EGM96空间坐标系下，空间分辨率为30米。根据《国家基本比例尺地图编绘规范》的要求^[33]，所采用的地形图平坦区域（坡度小于2°）和丘陵区域（坡度2～6°）垂直精度优于3～5米，对于山区和高山区（坡度大于6°）垂直精度优于8～14米。

SRTM DEM为美国地质勘探局提供的、未经“空洞”填补处理的SRTM1数据产品，其获取时间为2000年2月，分辨率为30米，采用WGS84/EGM96空间坐标系。在90%的置信区间内，SRTM1垂直精度优于16米，其高程精度受地形的影响较大，在平坦区域的高程精度可达10米，且随地表坡度变化^[34]。

TerraSAR-X/TanDEM-X数据为德国空间局提供的双星分布式X波段SAR数据，获取时间2014年2月18日和2014年3月13日。该数据获取模式为收发分置式（bistatic mode），即一颗卫星发射脉冲，两颗卫星同时接收，通过对两个回波信号相位差进行解缠来获取地面高程信息，其垂直方向的绝对精度优于10米，相对精度在局部坡度角大于20%时达到4米，小于等于20%时达到2米^[28]。

中国冰川编目数据由寒区旱区科学数据中心提供（http://westdc.westgis.ac.cn/data），主要用于冰川区与非冰川区的划分。另外，为了对TerraSAR-X/TanDEM-X数据提取的数字高程模型进行评估，本文还获取了美国国家冰雪数据中心（National Snow and Ice Data Center，NSIDC，http://nsidc.org/data/icesat/）发布的ICESat-GLAS数据，使用的是GLA14数据产品（Release 633），为全球陆地地面高程数据（图1）。

表1 数据集用到的地形图和TerraSAR-X/TanDEM-X影像列表

时间	数据源	ID	分辨率（m）
1980/10	地形图/ TOPO DEM	H47e016002/H47e016003/ H47e016004 H47e017002/H47e017003/ H47e017004 H47e017005/H47e018003/ H47e018004 H47e018005	30
2000/02/11–2000/02/22	SRTM C-band	-	30
2014/02/18	TSX/TDX	20140218T232940_20140218T232948	12
2014/03/13	TSX/TDX	20140313T113609_20140313T113617	12

图1 研究区范围以及所用地形图和TerraSAR-X/TanDEM-X影像的覆盖范围

1.2 数据处理步骤

详尽的数据处理方法可见Wu等人发表的文献^[22]。冰川表面高程变化获取主要由3部分组成：获取2014年数字高程模型，获取1980～2014年和2000～2014年两个时间段的冰川表面高程变化，获取1980～2000年的冰川表面高程变化。本文采用瑞士GAMMA软件平台和美国ESRI公司的ArcGIS软件平台进行数据的处理和分析，图2显示了主要的处理步骤。

合成孔径雷达差分干涉测量，首先基于TerraSAR-X将TOPO DEM和SRTM模拟转换为SAR成像数据，TerraSAR-X与TanDEM-X为同一成像时间的SAR图像，TerraSAR-X与TOPO DEM/SRTM模拟的SAR数据为不同成像时间的SAR图像，TerraSAR-X分别与TanDEM-X和TOPO DEM/SRTM模拟的SAR数据进行配准、共轭相乘、生成两幅干涉图，其中同一成像时间的干涉图相位包含地形高程信息，而不同成像时间的干涉图相位包含地形高程信息和两次成像期间地表形变信息。对两幅干涉图去平地效应、并对低相干及相位异常区域进行掩膜，通过相位解缠，将相位转换为高程差。最后利用非冰川区域双线性拟合进行整体改正，得到优化的高程差，即得到研究区高程变化信息。利用TerraSAR-X与TanDEM-X生成的干涉图，进行相位解缠、基线估计，得到反映地面高程信息的真实相位，根据InSAR成像基本原理和坐标转换，得到最新的DEM数据。

TOPO DEM与SRTM之间的处理采用大地测量法，两幅DEM数据在没有配准的情况下直接高程相减得到高程差值图，高程差值图与坡度坡向之间具有以下余弦函数关系：

\[\frac{\mathrm{d}\mathrm{h}}{\mathrm{tan}\left(\mathrm{\alpha }\right)}=\mathrm{a}\mathrm{ }\mathrm{*}\mathrm{ }\mathrm{cos}\left(\mathrm{b}-\mathrm{\phi }\right) +\mathrm{ }\frac{\stackrel{-}{\mathrm{d}\mathrm{h}}}{\mathrm{tan}\left(\mathrm{\alpha }\right)}\]

式中：dh为不同DEM数据之间各个像元的高程差，a是平移矢量的大小，b是平移矢量的方向，α和φ分别为像元对应的坡度和坡向，\(\stackrel{-}{\mathrm{d}\mathrm{h}}\)为不同DEM数据之间的整体高程差异（垂直矢量的大小）。坡度的大小对DEM数据匹配的精度有很大影响，因此要剔除地形平坦的区域和坡度太大的区域，本文剔除了地面坡度小于5°和大于80°的区域。不同DEM数据之间的高程差会存在异常值，影响余弦曲线的拟合，本文对高程差样本进行统计，认为低于5%或高于95%分位数的高程差是异常值，进行剔除。然后通过最小二乘法拟合dh/tan(α)与坡向的余弦曲线，得到系数a、b、c，计算不同DEM数据在X、Y和Z方向的平移量。最后利用非冰川区高程差与最大曲率的关系，校正由不同分辨率造成的冰川区的高程差异残差。最终得到的高程差值即为研究区高程变化信息。

图2 岗日嘎布地区最新数字高程模型及冰川表面高程变化主要流程示意图

2 数据样本描述

岗日嘎布地区冰川高程变化数据集包含两部分数据：一是基于TerraSAR/TanDEM雷达影像获取的2014年数字高程模型数据，命名为TanDEM.tif，存放在“数字高程模型”文件夹中；二是基于地形图、SRTM和TerraSAR/TanDEM雷达影像获取的1980～2014年冰川高程变化栅格数据，命名为SRTM_TOPO_DIFF.tif、TDX_SRTM_DIFF.tif和TDX_TOPO_DIFF.tif，表示不同DEM之间获取的研究区高程变化，单位是米，存放在“冰川表面高程变化”文件夹中。样本展示如图3所示。

图3 岗日嘎布地区数字高程模型与冰川表面高程变化

a：2014年数字高程模型数据；b：SRTM与TOPO DEM之间获取的1980～2000年冰川表面高程变化；c：TanDEM与SRTM之间获取的2000～2014年冰川表面高程变化；d：TanDEM与TOPO DEM之间获取的1980～2014年冰川表面高程变化

3 数据质量控制和评估

TerraSAR-X/TanDEM-X遥感影像为X波段数据，在获取数字高程模型的过程中，X波段对冰雪穿透深度较小。但由于山区其特殊的地形，对微波影像有前向收缩、顶底叠置的影响，地形阴影对数据精度也有影响。在相位解缠的过程中，通过掩膜（mask）剔除影响数据精度的区域，使得获取的数字高程模型存在空洞。由于空洞所占比例较小，不影响数据质量，可以忽略不计。在数据处理的过程中，影像配准也会产生误差，本文使用的GAMMA雷达数据处理平台的影像配准算法误差控制在0.01个像素范围内。

从美国国家冰雪数据中心获取的ICESat-GLAS数据，其中GLA14数据产品（Release 633）为全球陆地地面高程数据。基于TerraSAR-X/TanDEM-X获取的数字高程模型与GLA14点相比较，剔除超出±100米的异常值，在非冰川区的平均值与标准差为2.65±1.48米。

本文采用的SRTM1为C波段数据，而C波段能够穿透冰雪，在测定积雪覆盖的冰川表面地形时得不到真实地面信息。相比C波段数据，X波段对冰雪的穿透效应要小，利用邻近的、且配准和校正之后的SRTM C波段和X波段相比较，得到SRTM1在岗日嘎布地区平均穿透1.24米。由于SRTM1为未经“空洞”填补处理的数据产品，与TOPO DEM、TerraSAR-X/TanDEM-X相比较得到的冰川表面高程变化数据存在空洞。

利用TOPO DEM、SRTM和TerraSAR-X/TanDEM-X数据获取地面高程变化信息后可利用非冰川区的高程差异残差对DEM数据间相对误差进行评估（图4）。无冰区高程差的均方根误差（RMSE）或标准差（STDV）可以作为高程差不确定性的初步估计，但是会高估高程变化的不确定性。根据空间邻近相似性原理，相邻的像元间存在自相关，因此需要除去自相关造成的影响。为此，引入标准平均误差（SE）对高程差不确定性进行评估。最终不同DEM数据之间高程差的精度（σ）可以用无冰区高程差的均值MED和标准平均误差SE表示：

\[\mathrm{\sigma }=\sqrt{{\mathrm{M}\mathrm{E}\mathrm{D}}^{2}+{\mathrm{S}\mathrm{E}}^{2}}\]

\[\mathrm{S}\mathrm{E}=\mathrm{S}\mathrm{T}\mathrm{D}\mathrm{V}/\sqrt{\mathrm{N}}\]

式中：N为去自相关后的像元个数。SRTM与TOPO DEM之间的去空间自相关距离选择600米；TerraSAR-X/TanDEM-X与SRTM/TOPO DEM之间的去空间自相关距离选择200米。结果显示（表2），误差校正后，残差平均值（MED）减小到了−0.53～0.67米之间，趋近于0；TOPO、SRTM和TSX/TDX数据间相对误差介于0.43～0.77米之间。

图4 岗日嘎布地区SRTM与TOPO DEM差值前（a）、差值后（b）在非冰川区的高程差

表2 岗日嘎布地区DEM差值后非冰川区误差特征分布

数据集	残差平均值/MED (m)	残差标准差/STDV (m)	去自相关的像元个数/N	标准误差/SE (m)	高程差的不确定性/\(\mathrm{\sigma }\) (m)
SRTM - TOPO	0.67	16.41	1829	0.38	0.77
TanDEM - SRTM	-0.42	9.93	16356	0.08	0.43
TanDEM - TOPO	-0.53	5.36	16356	0.04	0.53

4 数据价值

TerraSAR /TanDEM采用X波段的合成孔径雷达传感器，在测算冰川表面高程变化时对大气干扰及冰雪穿透有较好的屏蔽效果。利用光学遥感影像或激光测高数据得到的本研究区冰川表面高程变化，与本文的结果相差较大。通过对比不同研究成果及研究方法发现，本文结果与其他研究成果相差较大的主要原因是SRTM C波段对冰雪穿透的估算差异。Kääb等^[27]在测算帕米尔–喀喇昆仑山–喜马拉雅山冰川物质平衡时认为，SRTM C波段在念青唐古拉山平均穿透7～10米。本文利用利用邻近的、且配准和校正之后的SRTM C波段和X波段相比较，得到SRTM1在岗日嘎布地区平均穿透1.24米。综合前人研究发现，在喜马拉雅山及念青唐古拉山中西段地区^{[17, 18, 29, 30]}，SRTM对冰雪的平均穿透1.24～2.5米，因此在念青唐古拉山东段岗日嘎布地区1.24米的穿透深度是可信的。

野外实测冰川物质平衡是衡量冰川变化的最佳指标，本文搜集了岗日嘎布地区两条典型冰川（帕隆藏布4号冰川和10号冰川）的实测物质平衡数据^{[4, 35]}，通过与野外实测数据相比较可知，基于SRTM和TerraSAR-X/TanDEM-X数据获取的帕隆藏布4号和10号冰川表面高程变化转换的物质平衡，与野外实测具有较高的一致性，基于遥感数据获取的岗日嘎布地区冰川表面高程变化精度较高。

冰川表面高程变化数据是研究冰川物质平衡、冰川运动及区域水文变化特征等重要基础数据。对冰川表面高程变化数据集的充分利用，结合气候、水文数据的综合探讨，对区域乃至全球水资源分布和变化的研究、以及地理环境变化的特点和趋势分析有着重要意义。

5 数据使用方法和建议

岗日嘎布地区2014年数字高程模型及1980～2014年冰川表面高程变化数据集通过CERN综合中心数据资源服务网站提供全部数据下载服务。本数据集提供了ArcGIS软件平台下生产的32位浮点型GeoTIFF栅格数据，除了*.tif数据本身之外，还有包括了影像坐标文件*.tfw和投影、金字塔信息*.tfw.aux.xml。可以在ArcGIS、SuperMap等GIS平台软件直接调用，也可以使用Matlab等数据分析软件编码进行批量处理和分析。

致谢

感谢德国空间局（DLR）提供的TerraSAR-X/TanDEM-X数据，感谢美国地质调查局（USGS）提供的SRTM1数据，感谢中国科学院寒区旱区研究所科学数据中心提供的冰川编目数据集。

[1]

LI X, CHENG G, JIN H, et al. Cryospheric change in China[J]. Global and Planetary Change, 2008, 62(3): 210-218.