1967–2014年科西河流域面积≥0.05 km2冰湖数据集
| 1967–2014年科西河流域面积≥0.05 km2冰湖数据集 作者：张大弘 宫鹏 姚晓军 刘娟 2019年11月28日 |
 张大弘, 宫鹏, 姚晓军, 刘娟. 1967–2014年科西河流域面积≥0.05 km2冰湖数据集[J/OL]. 中国科学数据, 2019. (2019-11-27). DOI: 10.11922/csdata.2019.0064.zh. |
摘要&关键词
[编辑]摘要：科西河流域发育有众多的冰川和冰湖，冰川覆盖率达9.78%，冰湖数量221个，冰湖分布及变化不仅是预测该区域冰湖溃决洪水灾害的基础数据，而且可以在一定程度上反映高山寒区气候及冰川水资源的变化。本数据集应用1960s–2010s多源遥感影像数据（Corona、Landsat MSS/TM/ETM+/OLI）及地形图进行目视解译获得了科西河流域1960s、1970s、1980s、1990s、2000s和2010s六期冰湖数据。研究区空间范围为26.86°–29.14°N，85.39°–88.95°E，覆盖印德拉瓦迪、孙科西、绒辖曲、利库科拉、牛奶河、阿润和塔木尔等子流域。本数据集包括流域内面积≥0.05 km2的冰川湖，数据集时序较长且包含冰湖类别属性信息，可作为科西河流域冰湖时空变化、冰湖水量变化评估的基础数据，也可作为喜马拉雅山地区冰湖溃决危险性评估、水资源变化、气候变化及冰湖与冰川间变化响应关系等研究的重要依据。
关键词：冰湖；科西河；目视解译；青藏高原
Abstract & Keywords
[编辑]Abstract: There are many glaciers and glacial lakes in the Koshi River basin. The coverage rate of glaciers reaches 9.78% and the number of glacial lakes is 221. The distribution and change of glacial lakes in this area are the main factors for predicting the disaster of glacial lake outburst floods (GLOFs), and can also reflect the changes of climate and glacial water resources in Alpine frigid areas to a certain extent. The dataset includes the glacial lake data of six periods in the Koshi River basin, i.e., 1960s, 1970s, 1980s, 1990s, 2000s and 2010s, which was obtained by visual interpretation using multi-source remote sensing image (Corona, Landsat MSS/TM/ETM+/OLI) and topographic maps during 1960s to 2010s. Lying between 26.86°–29.14°N and 85.39°–88.95°E , the study area is composed of several sub basins including Indrawati, Sun Koshi, Tama Koshi, Likhu, Dudh Koshi, Arun and Tamor. The dataset includes the vector of glacial lakes with an area above 0.05 km2 in the Koshi River basin. The dataset has a long time series, which can be used as the basic data for the assessment of the temporal and spatial changes of glacial lakes and the changes of water volume of glacial lakes in the Koshi River basin. In addition, it can also be used as an important basis for the risk assessment of glacial lake outburst, the changes of water resources, climate change and the response relationship between glacial lakes and glaciers in this region.
Keywords: glacial lake; Koshi River; visual interpretation; Tibetan Plateau
数据库（集）基本信息简介
[编辑]| 数据库（集）名称 | 1967–2014年科西河流域面积≥0.05 km2冰湖数据集 |
| 数据作者 | 张大弘，宫鹏，姚晓军，刘娟 |
| 数据通信作者 | 姚晓军（yaoxj_nwnu@163.com） |
| 数据时间范围 | 1967–2014年 |
| 地理区域 | 地理范围包括26.86°–29.14°N，85.39°–88.95°E，北至雅鲁藏布江南部，南到尼泊尔与印度交界处，西临加德满都，东达中、尼、印三国接壤地区。 |
| 空间分辨率 | 全色1.8 m / 15 m、多光谱78 m / 30 m、地形图1:50000/1:100000 |
| 数据量 | 0.85 MB |
| 数据格式 | ESRI Shapefile文件（压缩为.zip格式） |
| 数据服务系统网址 | http://www.sciencedb.cn/dataSet/handle/885 |
| 基金项目 | 国家自然科学基金（41561016、41861013、41801052）；西北师范大学青年教师科研能力提升计划项目（NWNU-LKQN-14-4）。 |
| 数据库（集）组成 | 数据集共包括2个数据文件：（1）Basin_Koshi.zip是科西河流域边界数据，数据量0.1 MB；（2）Glacial_lake_Koshi.zip是1967–2014年科西河流域面积≥0.05 km2的冰湖数据，数据量0.75 MB。 |
Dataset Profile
[编辑]| Title | A dataset of glacial lakes in the Koshi River basin with an area above 0.05 km2 (1967 – 2014) |
| Data corresponding author | Yao Xiaojun (yaoxj_nwnu@163.com) |
| Data authors | Zhang Dahong, Gong Peng, Yao Xiaojun, Liu Juan |
| Time range | 1967 – 2014 |
| Geographical scope | Latitude & longitude (26.86°–29.14°N, 85.39°–88.95°E);located in the Central Himalayas; The specific areas: north to the south of the Yarlung Zangbo River, south to the border of Nepal and India, West to Kathmandu in India, East to the contiguous areas of China, Nepal and India. |
| Spatial resolution | Pan 1.8 m / 15 m 、Multispectral 78 m / 30m、Topographic maps 1:50000 / 1:100000 |
| Data volume | 0.85 MB |
| Data format | ESRI Shapefile file (compressed in *.zip) |
| Data service system | <http://www.sciencedb.cn/dataSet/handle/885> |
| Source(s) of funding | National Natural Science Foundation of China (41561016, 41861013, 41801052); Earlier Career Research Promotion Program of Northwest Normal University (NWNU-LKQN-14-4). |
| Dataset composition | The dataset consists of 2 subsets in total: (1) Basin_Koshi.zip is made up of boundary data of regions in the Koshi River Basin, with a data volume of 0.1 MB; (2) Glacial_lake_Koshi.zip is made up of boundary data of glacial lakes in the Koshi River Basin with an area above 0.05 km2 from 1967 to 2014, with a data volume of 0.75 MB. |
引 言
[编辑]冰湖溃决洪水（Glacial Lake Outburst Flood，GLOF）是高山冰川作用区常见的自然灾害之一[1][2]，主要类型是冰川阻塞湖和冰碛阻塞湖。就中国西藏自治区来说，20世纪以来发生的冰湖溃决洪水（泥石流）事件近30起[3]；欧洲的阿尔卑斯山、南美洲的安第斯山、中亚的天山以及青藏高原周边高山地区也是冰川溃决洪水的主要分布区[4]。冰湖溃决有可能造成大量的人员伤亡、牲畜死亡、毁坏大面积的草场或作物以及桥梁建筑等，尽可能预测并减少冰湖溃决的损失是当前重要任务，而冰湖分布及变化正是预测冰湖溃决洪水灾害的基础数据。冰湖扩张或萎缩可以灵敏地反映气候变化[5]，其面积或水量的变化能从一定程度上反映区域水资源分布的变化。因此，冰湖数据集是研究冰湖分布与变化、预测冰湖溃决灾害、区域水资源分布及气候变化的基础数据。
关于冰湖的概念，学者们给出的定义不尽相同[6][7][8][9][10][11][12]，但都强调冰湖是冰川作用所形成的，区别主要是冰川作用时期与湖泊物质来源。例如，百度百科（https://baike.baidu.com）将冰湖定义为由冰川挖蚀形成的洼地（冰斗）和冰碛物堵塞冰川槽谷积水而成的一类湖泊。姚治君等[6]认为冰湖是末次冰期以来冰川运动或退缩产生的融水在冰川表面、前部或侧部汇集形成的湖泊，此定义强调冰湖与冰川间的地理位置关系。姚晓军等[9]从冰湖编目和冰湖溃决洪水角度提出冰湖是以现代冰川融水为主要补给源或在冰碛垄洼地内积水形成的天然水体，并给出了冰湖类型及遥感判识指标。考虑冰湖识别的易识别性、可操作性及一致性，本文采用姚晓军等[9]提出的定义及冰湖分类体系。喜马拉雅山中段是冰湖溃决最频繁的地区之一[13]，历史上该地区曾发生过的冰湖堰塞湖溃决事件最为集中的区域是科西河流域上游地区[14][15]。作为恒河的一个重要支流，科西河流经中国、尼泊尔和印度3个国家，流域内具有世界上最大的高差，最为完整的气候、土壤和植被带谱。流域内分布有多座8000 m以上山峰，现代冰川发育，冰川覆盖率达9.78%。据统计，科西河流域历史上发生过17次较大规模的冰湖溃决事件，其中中国境内高达12次[16]。2015年尼泊尔博克拉地区发生8.1级地震，Kargel等[17]对震后灾区分析时指出，未来当该区域发生更大震级的地震或震源靠近冰湖分布区时，发生严重冰湖溃决洪水灾害的可能性依然存在。除地震外，雪（冰）崩、强降水、冰川强烈消融、冰碛坝内死冰消融等亦是引发冰湖溃决的诱因[18][14]。因此，建立长时间序列的冰湖数据集可为分析科西河流域内冰湖时空变化特征、冰湖危险性评估及冰湖溃决洪水（泥石流）灾害等研究提供基础数据支撑和科学依据。
1 数据采集和处理方法
[编辑]1.1 数据来源与预处理
[编辑]为获取科西河流域不同年代的冰湖数据，本研究选用的遥感影像数据源主要包括 USGS/NASA（http://earthexplore.usgs.gov）提供的14景Corona卫星影像和56景Landsat MSS/TM/ETM+/OLI遥感影像（表1）。其中，Corona卫星影像成像时间集中在1960s，是该时段可获得的质量较好的遥感影像；Landsat MSS/TM/ETM+遥感影像时间为1970s、1980s、1990s、2000s，Landsat OLI遥感影像则反映了2013年以来该流域冰湖状况。以Landsat MSS/TM/OLI数据作为主要数据源，可以免费获取大量的数据，且回访周期短，利用成像时间接近的影像通过交叉验证解译的方法可以避免云和山体阴影等因素的干扰，提高湖泊解译精度。为保证数据的准确性，本研究同时选用了该流域1966–1980年的83幅地形图，其中包括4幅1∶50万、22幅1∶10万、48幅1∶5万和7幅1∶2.5万地形图。受遥感影像空间分辨率限制，按Landsat MSS影像空间分辨率（78 m）计算，0.006 km2是像元解译的理论阈值，因此本数据集仅考虑面积≥0.05 km2的冰湖。此外，研究区所用的冰川数据从科技部科技基础性工作专项“中国冰川资源及其变化调查”项目组和国际山地中心（ICIMOD）获取[19]。
表1 数据集所采用的遥感影像列表
| 影像编号/轨道号 | 卫星/传感器 | 分辨率/m | 影像获取时间 |
| DS1037-1039DA026 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1037-1039DA027 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1037-1039DA028 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1044-1070DA095 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1044-1070DA096 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1044-1070DA097 | Corona KH-4 | 全色1.8 m | 1967-11-11 |
| DS1108-2217DA063 | Corona KH-4 | 全色1.8 m | 1969-12-18 |
| DS1108-2217DA064 | Corona KH-4 | 全色1.8 m | 1969-12-18 |
| DS1108-2217DA065 | Corona KH-4 | 全色1.8 m | 1969-12-18 |
| DS1108-2217DA066 | Corona KH-4 | 全色1.8 m | 1969-12-18 |
| DS1049-2136DA148 | Corona KH-4 | 全色1.8 m | 1968-12-21 |
| DS1108-2233DA055 | Corona KH-4 | 全色1.8 m | 1968-12-21 |
| DS1014-2218DA192 | Corona KH-4 | 全色1.8 m | 1968-12-21 |
| DS1014-2218DA068 | Corona KH-4 | 全色1.8 m | 1972-04-21 |
| 139/040 | Landsat MSS | 78 | 1976-12-18 |
| Landsat TM | 30 | 1989-12-28; 1990-01-13; 2000-10-07; 2000-11-08; 2011-08-27 | |
| Landsat OLI | 30 | 2013-06-13; 2014-01-07 | |
| 139/041 | Landsat MSS | 78 | 1976-12-18 |
| Landsat TM | 30 | 1989-10-25; 1989-11-10; 1999-09-19; 2000-12-26; 2011-08-27 | |
| Landsat OLI | 30 | 2013-06-13; 2013-11-20; 2013-12-06 | |
| 140/040 | Landsat MSS | 78 | 1976-12-19 |
| Landsat TM | 30 | 1989-01-17; 1992-09-30; 1999-11-29; 2000-09-12;2000-09-28; 2010-10-26; 2010-11-11; 2010-04-09; 2010-04-25 | |
| Landsat OLI | 30 | 2013-11-11; 2014-05-22 | |
| 140/041 | Landsat MSS | 78 | 1976-12-19; 1977-03-19 |
| Landsat TM | 30 | 1989-01-17; 1992-09-22; 1999-10-28; 2000-10-30; 2009-10-15 | |
| Landsat OLI | 30 | 2014-01-05; 2014-05-22 | |
| 141/040 | Landsat MSS | 78 | 1976-12-02 |
| Landsat TM | 30 | 1991-11-30; 1993-03-08; 1999-12-06; 2000-10-05;2000-11-22; 2009-06-08; 2009-10-30; 2010-10-01 | |
| Landsat OLI | 30 | 2013-12-04 | |
| 141/041 | Landsat MSS | 78 | 1976-12-02 |
| Landsat TM | 30 | 1991-11-30; 2001-02-26; 2001-12-27; 2009-05-07; 2010-11-02 | |
| Landsat OLI | 30 | 2014-04-11; 2014-05-13 |
冰湖提取会受到云、雪及山体阴影的影响，在选取遥感影像时遵循少云多期原则，以达到多时段多角度的交叉解译目的，从而尽量避免漏画、错画问题。除Corona遥感影像需进行几何校正和图幅拼接外，获取的Landsat卫星系列遥感影像在数据生产时已经过初步预处理（系统辐射校正、地面控制点几何校正及基于DEM的地形校正），经检查影像几何位置匹配较好，目视解译又属于定性判别地类，故这些遥感影像未做其他预处理工作。
1.2 研究区概况
[编辑]科西河流域位于青藏高原南部、喜马拉雅山中段（26.86°–29.14°N，85.39°–88.95°E），北抵雅鲁藏布江，南至尼泊尔与印度交界处，西临加德满都附近，东达中、尼、印三国接壤地区。地势北高南低，海拔落差极大，主河全长255 km，上游分Arun（中国境内称为朋曲）、Sun Koshi（中国境内称为波曲）和Tamor三个支流（图1）。流域东西长340 km，南北宽230 km，总面积约5.5×104 km2，位于中国和尼泊尔境内的面积分别占流域总面积的51.96%和48.04%[20]。流域内分布有5座海拔8000 m 以上山峰，围绕这些山峰发育了众多的冰川及冰湖。受地形影响，研究区气候类型差异明显，自北向南依次分布着喜马拉雅山过渡带、高山带、亚高山带、温带、亚热带和热带等，后三种气候类型分布在河流中下游，比较适宜农业发展，分布着大量人口。据统计，尼泊尔境内总人口434.59万人，占流域总人口的96.57%[18]。
图1 研究区概况图
1.3 数据处理过程
[编辑]冰湖作为一类自然水体，目前已有一些比较成熟的基于遥感影像的自动提取算法，如水体指数法、波段比值法、监督分类法及面向对象分割法等。由于本数据集采用的本底数据存在全色影像、地形图等多源数据，加之自动提取算法对影像质量要求较高，为保证数据采集精度，本数据集统一采用目视解译方法，具体步骤如下：
（1）获取解译底图：下载、收集并整理科西河流域的遥感影像、地形图等可提供解译对象信息的资料，按时间、类型、位置等归类，确保每个时段每个位置至少有一期影像覆盖。
（2）确定解译方式：因冰湖具有季节性变化特征，一年之内最稳的时段为9–12月，且少量影像（如Landsat ETM+）局部存在坏带，故以其他相邻时相的影像作为补充，通过交叉检验解译方式获取精度较高的湖泊信息。
（3）湖泊边界目视解译：以交叉检验解译的方式逐时期进行，显示比例尺控制在1:2000左右，沿水陆、冰陆或者冰雪界线像元中心数字化冰湖边界，精度控制在一个像元以内。
（4）确定冰湖与非冰湖：采用姚晓军等[4]提出的冰湖定义、冰湖与非冰湖的界定方法，逐一对不同时期提取到的研究区湖泊数据进行判定，剔除由大气降水、地下水等形式补给的非冰湖，包括研究区内大量存在的古冰蚀湖，最终得到研究区不同时期的冰湖数据集。
（5）湖泊数据质量检查：逐一检查各湖泊边界图形数据、对应的影像数据及属性数据等，确保无漏画错画、各数据项无遗漏和录入正确等。
（6）冰湖类型识别：将最终获得的冰湖数据导入Google Earth软件，根据冰湖分类体系[4]，逐一识别出所有冰湖的类型并写入属性表。
2 数据样本描述
[编辑]图2展示了科西河流域冰湖整体分布情况。冰湖多边形矢量数据以gl_XXXX形式命名，其中XXXX表示年份；冰湖点矢量数据gl_LocationCode给出了流域内所有的冰湖位置、编码、类型等信息。至2014年，流域内共发育冰川湖221个，其中冰川侵蚀湖83个（冰斗湖10个，其他冰川侵蚀湖73个），冰碛阻塞湖128个（终碛阻塞湖115个，侧碛阻塞湖13个），冰面湖4个，其他冰川湖6个。以面积持续扩张的冰湖GL085838E28322N为样本，图3展现了该湖泊在研究时序内的边界、周长及面积等的变化情况。图4给出了数据集中面积先减小后持续扩张的龙巴萨巴湖各时序内的边界分布情况。
图2 冰湖分布数据集示意图
图3 冰湖GL085838E28322N
图4 龙巴萨巴湖
在本数据集中，每个冰湖对象由图形数据和属性数据两部分组成，其中属性数据包括FID、形状类型（Shape）、冰湖编码（Code）、冰湖名称（Name_gl）、冰湖周长（Perimeter_gl）、冰湖面积（Area_gl）以及冰湖类型（Type_gl）等。冰湖编码结构为：GLnnnnnnEmmmmmN，其中，GL表示冰湖，n表示冰湖质心经度，m表示冰湖质心纬度，E、N分别代表东经和北纬，n的前三位和m的前两位分别表示经度、纬度的整数部分，后三位表示经纬度小数前三位。为保证冰川编码的一致性，以2010s时期冰湖质心为准。如冰湖在2010s时期不存在，则以第一次出现年代对应的冰湖质心进行编码。
3 数据质量控制和评估
[编辑]解译时依据姚晓军[4]提出的冰湖定义、冰湖与非冰湖的界定方法，遵循马荣华[21]、宫鹏[20]等所使用的湖泊边界遥感判译原则，精度控制在一个像元以内。参与数据集解译的长时序多源遥感数据精度存在差异，就空间分辨率而言，Corona KH-4数据1.8 m分辨率（全色，1.8 m），Landsat MSS数据（多光谱，78 m）、Landsat TM数据（多光谱，30 m）、Landsat ETM+（多光谱30 m，全色15 m，部分数据存在坏带）、Landsat OLI（多光谱30 m，全色15 m）以及参考的地形图（1:50000、1:100000）。因此本数据集的精度虽不及其他高分辨率、单一数据源生产的数据，整体精度仅能保证在最低分辨率78 m左右，但从数据时序长度来看，本数据集可以满足大多数相关领域的研究。
4 数据价值
[编辑]研究区分布有大量湖泊，面积小而数量多，故通过区分冰湖与非冰湖以及识别各个冰湖的类型[4]所得到的冰湖数据及湖泊类别属性是本数据集的价值所在，可以为当地相关部门以及从事冰湖溃决危险性评价工作的学者们提供宝贵的参考。
本数据集是依据不同时期多源遥感影像数据（1967–2014）生产的6期冰湖边界矢量数据，特点在于较长时间序列，可作为科西河流域冰湖时空变化、冰湖水量变化评估的基础数据，也可作为该区域冰湖溃决危险性评估、水资源变化、气候变化及冰湖与冰川间变化响应关系等研究的重要依据。
5 数据使用方法和建议
[编辑]本数据集解压后为ESRI Shapefile文件格式，可使用支持SHP格式的GIS与遥感软件读取与操作。数据集直接反映了1967–2014年科西河流域的冰湖分布特征，进一步分析处理可获得流域内冰湖的时空变化特征等信息，结合相关资料，亦可用作冰湖危险性评估及冰湖溃决洪水（泥石流）灾害的基础数据。另外，如果还有其他来源的数据，可对本数据做更新修改，也可用最新的数据更新2010s数据，以延长数据时序。
参考文献
[编辑]- ↑ OSTI R, EGASHIRA S. Hydrodynamic characteristics of the Tam Pokhari glacial lake outburst flood in the Mt. Everest region, Nepal[J]. Hydrological Processes, 2009, 23(20): 2943-2955.
- ↑ RICHARDSON S D, REYNOLDS J M. An overview of glacial hazards in the Himalayas[J]. Quaternary International, 2000, 65/66: 31-47
- ↑ 姚晓军, 刘时银, 孙美平, 等. 20 世纪以来西藏冰湖溃决灾害事件梳理[J]. 自然资源学报, 2014, 29(8): 1377-1390.
- ^ 4.0 4.1 4.2 4.3 4.4 CUI P, DANG C, CHENG Z, et al. Debris flows resulting from glacial-lake outburst floods in Tibet, China[J]. Physical Geography, 2010, 31(6): 508-527.
- ↑ 王欣, 刘时银, 莫宏伟, 等. 我国喜马拉雅山区冰湖扩张特征及其气候意义[J]. 地理学报, 2011, 66(7):895-904.
- ^ 6.0 6.1 姚治君, 段瑞, 董晓辉, 等. 青藏高原冰湖研究进展及趋势[J]. 地理科学进展, 2010, 29(1): 10-14.
- ↑ 陈晨, 郑江华, 刘永强, 等. 近20年中国阿尔泰山区冰川湖泊对区域气候变化响应的时空特征[J]. 地理研究, 2015, 34(2): 270-284.
- ↑ 唐升贵, 刘发祥, 赵振远. 西藏念青唐古拉山东段冰湖的时空分布特征及其影响因素探讨[J]. 科技资讯, 2014, 12(16): 35-36.
- ^ 9.0 9.1 9.2 姚晓军, 刘时银, 韩磊, 等. 冰湖的界定与分类体系——面向冰湖编目和冰湖灾害研究[J]. 地理学报, 2017, 72(7): 1173-1183.
- ↑ 曹学诚, 刘周周, 李维胜. 尼泊尔4·25地震震前冰湖制图与潜在危险性分析[J]. 冰川冻土, 2016, 38(3): 573-583.
- ↑ 中国科学院-水利部成都山地灾害与环境研究所, 西藏自治区交通厅科学研究所. 西藏泥石流与环境[M]. 成都: 成都科技大学出版社, 1999: 79.
- ↑ 秦大河, 姚檀栋, 丁永建, 等. 冰冻圈科学辞典[M]. 北京: 气象出版社, 2014: 2-54.
- ↑ REYNOLDS J M. Glacial-lake outburst floods（GLOF）in the Himalayas: an example of hazard mitigation from Nepal[J]. Geoscience and Development, 1995, 2: 6-8.
- ^ 14.0 14.1 RICHARDSON S D, REYNOLDS J M. An overview of glacial hazards in the Himalayas[J]. Quaternary International, 2000, 65/66, 31-47.
- ↑ LIU C H, SHARMA C K. Report on First Expedition toGlaciers and Glacier Lakes in the Pumqu (Arun) and Poqu (Bhote-SunKosi)River basins, XiZang (Tibet), China[M]. Beijing: Science Press, 1988.
- ↑ 邓伟, 张镱锂. 气候变化下Koshi河流域资源、环境与发展[M]. 成都: 四川科学技术出版社, 2014:82.
- ↑ KARGEL J S, LEONARD G J, SHUGAR D H, et al. Geomorphic and geologic controls of geohazards inducedby Nepal's 2015 Gorkha earthquake[J]. Science, 2016, 351(6269): aac8353．
- ^ 18.0 18.1 YAMADA T, SHARMA C K. Glacier lakes and outburst floods in the Nepal Himalaya[J]. IAHS, 1993, 218: 319-330.
- ↑ 刘时银, 姚晓军, 郭万钦, 等. 基于第二次冰川编目的中国冰川现状[J]. 地理学报, 2015, 70(1): 3-16.
- ^ 20.0 20.1 宫鹏, 姚晓军, 孙美平. 1967–2014年科西河流域冰湖时空变化[J]. 生态学报, 2017, 37(24): 8422-8432．
- ↑ 马荣华, 杨桂山, 段洪涛, 等. 中国湖泊的数量、面积与空间分布[J]. 中国科学: 地球科学, 2011, 41(3): 394-401.
数据引用格式
[编辑]张大弘, 宫鹏, 姚晓军, 刘娟. 1967–2014年科西河流域面积≥0.05 km2冰湖数据集[DB/OL]. Science Data Bank, 2019. (2019-09-05). DOI: 10.11922/sciencedb.885.
