干旱区研究 ›› 2022, Vol. 39 ›› Issue (3): 676-683.doi: 10.13866/j.azr.2022.03.02

• 泛第三极环境与绿色丝路 • 上一篇    下一篇

中亚高山冰川表面高程变化时序重建

都伟冰1,2,3(),张世琼1(),李均力2,3,包安明2,3,王双亭1,史宁可1,许琳娟4,高鑫1,马丹丹1,郑岩超1   

  1. 1.河南理工大学测绘与国土信息工程学院,河南 焦作 454003
    2.中国科学院新疆生态与地理研究所荒漠与绿洲生态国家重点实验室,新疆 乌鲁木齐 830011
    3.新疆遥感与地理信息系统应用重点实验室,新疆 乌鲁木齐 830011
    4.黄河水利委员会黄河水利科学研究院,河南 郑州 450003
  • 收稿日期:2021-08-01 修回日期:2021-10-20 出版日期:2022-05-15 发布日期:2022-05-30
  • 通讯作者: 张世琼
  • 作者简介:都伟冰(1985-),男,博士,主要从事多源遥感冰川监测研究. E-mail: dwb@hpu.edu.cn
  • 基金资助:
    中国科学院创新交叉团队资助项目(JCTD-2019-20);新疆天山创新团队(2020D14016);国家自然基金项目(U2003201);国家自然基金项目(41671 034);国家自然基金项目(41975036);国家自然基金项目(42075132);国家自然基金项目(41601364);国家自然基金项目(42041004);智慧中原地理信息技术协同创新中心项目(2020C002);中央级公益性科研院所基本科研业务费专项(HKY-JBYW-2018-03);新疆维吾尔自治区优秀博士后资助;河南省科技攻关项目(222102320306);教育部产学合作协同育人项目:基于PIE的遥感信息智能提取课程群建设(202102245019);河南理工大学博士基金(B2022-8)

Temporal reconstruction of alpine glacier surface elevation variation in Central Asia

DU Weibing1,2,3(),ZHANG Shiqiong1(),LI Junli2,3,BAO Anming2,3,WANG Shuangting1,SHI Ningke1,XU Linjuan4,GAO Xin1,MA Dandan1,ZHENG Yanchao1   

  1. 1. College of Survey and Territory Information Engineering, Henan Polytechnic University, Jiaozuo 454003, Henan, China
    2. State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, Xinjiang, China
    3. Key Laboratory of GIS & RS Application Xinjiang Uygur Autonomous Region, Urumqi 830011, Xinjiang, China
    4. Yellow River Institute of Hydraulic Research, Zhengzhou 450003, Henan, China
  • Received:2021-08-01 Revised:2021-10-20 Online:2022-05-15 Published:2022-05-30
  • Contact: Shiqiong ZHANG

摘要:

中亚高山冰川区地形复杂,站点观测和传统实地测量范围十分有限。卫星激光测高技术可实现大范围冰川表面高程变化监测。以2003—2009年ICESat激光测高数据为数据源,参考2000年的SRTM高程数据,建立冰川区点云去噪及其精度优化算法和多尺度冰川区表面高程时空变化分析模型,并分析了2003—2009年间中亚山区整体与各分区冰川表面高程时序变化。结果表明:中亚高山冰川区的冰川表面平均高程整体呈下降趋势,表现出明显的区域差异。其中,2003—2009年中亚冰川表面高程总体下降了9.59±1.89 m;I区(即西藏和青海南部)的冰川表面高程下降了6.51±2.9 m;II区(即新疆、青海北部和甘肃部分地区)下降了7.87±5.03 m;III区(即中国境外,中亚地区的部分国家)下降了9.81±5.1 m,且监测到2004—2005年冰川表面高程上升。本研究方法对冰川区点云类高程脚点监测具有一定的通用性,但对基准DEM的依赖度较高。

关键词: 冰川表面高程, ICESat, 卫星激光测高, 时序重建, 高山冰川, 中亚

Abstract:

The topography of the alpine glacier area in Central Asia is complex. Satellite laser altimetry technology can be used to monitor large-scale glacier surface elevation change. This study used ICESat laser altimetry data from 2003 to 2009, and shuttle radar topography mission (SRTM) elevation data and Randolph glacier inventory (RGI) 6.0 cataloging data for 2000 to monitor the surface elevation change of alpine glaciers in Central Asia. The glacial region of Central Asia is divided into three regions according to geographical form. First, the data were tested for normality to eliminate error. When a normal distribution was not observed, outlier values from the standard deviation of the linear regression were selected as abnormal. Point cloud denoising and its accuracy optimization algorithm and a multi-scale analysis model of temporal and spatial variation in surface elevation in the glacier area were established. A cubic polynomial model was used to fit the glacier surface elevation data for different scales across the whole area and sub regions of the high Asian mountains. A time series reconstruction of changes in alpine glacier surface elevation in Central Asia was carried out based on the data for each region from 2003. Temporal changes in glacier surface elevation from 2003 to 2009 were analyzed. The results showed that the average elevation for the whole region decreased, with clear regional differences. From 2003 to 2009, the surface elevation of high Asian glaciers decreased by 9.59 ± 1.89 m. The glacier surface elevation in Areas I (i.e., Tibet and southern Qinghai), II (i.e., Xinjiang, Northern Qinghai, and some parts of Gansu), and III (i.e., outside China and some countries in Central Asia) decreased by 6.51 ± 2.9 m, 7.87 ± 5.03 m, and 9.81 ± 5.1 m, respectively. The glacier surface elevation increased from 2004 to 2005. Area I showed the slowest decline in glacier elevation, followed by Area II; Area III showed the fastest decline. This research method has universal application potential for monitoring point cloud elevation in glacier areas. The model evaluation parameter (R2) was >0.98, indicating that the cubic polynomial relationship between ICESat data and SRTM elevation data showed strong universality in this area. However, this will make the point cloud data more sparse, and the model is highly dependent on the benchmark DEM.

Key words: glacier surface elevation, ICESat, satellite laser altimetry, temporal reconstruction, alpine glacier, Central Asia