【文献阅读】锡金地区的冰湖溃决风险预测（Remya S N等人，2019.01，Current science）

一、文献情况介绍

文章题目：《Volume estimation of existing and potential glacier lakes: A case study in Sikkim Himalaya》

这是2019年1月发表在《Current science》上的一篇文章，文章的质量一般，但还是有一些值得学习的地方。查了一下这个期刊，是印度办的，影响因子为0.88，中科院综合类期刊分区4区，网友评论说大概3-6周就能接收，审稿人有时候只有1位。

文献的引用格式：

Remya & Kulkarni, A & Srinivasalu, Pradeep & G Shrestha, D. (2019). Volume estimation of existing and potential glacier lakes: A case study in Sikkim Himalaya. Current science. 

文献的链接：

文章没有给出DOI

文献的下载地址：

https://www.researchgate.net/publication/330081428_Volume_estimation_of_existing_and_potential_glacier_lakes_A_case_study_in_Sikkim_Himalaya

二、文章导读

这里先放上作者的摘要部分吧，毕竟是整个文章的核心。

Glacial lake outburst floods (GLOFs) are a hazard commonly reported in the glaciated terrain of the Himalaya. Glacier lakes can form if the glaciers retreat and the bottom topography overdeepens. We have adopted a technique to estimate the depth and volume of lakes using parameters such as glacier surface velocity, slope and laminar flow of ice. The technique is now automated using python programming language. The method was applied in Sikkim Himalaya to map potential lake sites and also to predict further expansion of existing glacier lakes. The studies were carried out for ten glaciers in the Tista river basin. The analysis suggests nine potential lake sites including further expansion of four existing glacier lakes. To validate results, the model lake boundary in 2001 was compared with the satellite observed value of 2015 and field measurements were made at the South Lhonak lake. The volume of the South Lhonak lake (in 2015) was estimated as 60 10.8 million m^3 ; with prolonged retreat of the glacier, the lake is likely to expand to a maximum volume of 90 16.2 million m^3. The above technique can provide new information to planners, leading to a more realistic approach in understanding the disaster potential of glacier lakes.

作者在摘要里面说到了，这篇文章提出用于估算湖泊的体积和深度的方法，用到的参数包括冰川表面流速，坡度等；并将该方法用于锡金地区，进行潜在湖泊制图和冰湖扩张预测。

三、文章详细介绍

1. 背景

随着全球气候变暖，自然灾害频繁发生。比如，冰湖溃决现象，1981、1987、2013年均有发生。而基于遥感技术的宏观调查发现，以下这些地区的冰湖存在溃决风险：

Satluj, Chenab, Ganga, Tista, Chandra, Zanskar and Central Himalayan river basins where a large numbers of glacie.

现在的关键问题是如何评估冰湖的溃决危害，传统的面积-体积法估算湖泊体积的误差比较大，作者采用了一种计算湖泊体积的方法，并应用在了锡金地区的Tista流域（Tista river basin）。

2. 实验数据

1）研究区数据

首先是研究区，是锡金地区的Tista流域，覆盖面积为136平方千米，该流域内包含多条冰川：

所使用影像为Landsat ETM+ 2001年和2002年数据，Landsat 8 OLI 2015年数据，数据来源为USGS（链接：https://earthexplorer.usgs.gov/）；冰川边界数据为RGI（Randolph Glacier Inventory (RGI) ）和ASTER DEM数据，DEM数据的水平方向精度为30米，竖直方向为17米，这两个数据的来源为：http://reverb.echo.nasa.gov/reverb

2）相关数据（in-situ data）

深度数据（bathymetric data）是从一篇文献《South Lhonak glacial lake outburst flood, North Sikkim》中获得；2014年8月28-30号之间还进行了实地取样，实地调查数据包括湖泊深度，电阻测量（用于推测终碛的内部岩性和组成），电阻率用以表示沉淀结构和冰碛下的分布，200个GPS点的随机深度测量（使用仪器为Hondex PS-7 Depth Sounder，该仪器的最大测量值为80米），GPS测量点位的误差为1-3m，仪器测深误差为2-3米，2016年8月28日还进行了一次补充测量，使用仪器为（Hondex HE-770，测量范围为3-500m）。

3. 方法

先来看以下作者给出的流程图：

大致可以分为以下几步：1）初始化冰川；2）用COSI-Corr方法估算冰川表面流速，用冰层流动（laminar flow of ice
）来建模深度；3）用DEM和深度来估算底部地形和岩床的过侵蚀；4）通过比较2001和2015年数据来验证结论。

对于冰川表面流速，用15m的landsat全色波段相对进行配准，再计算亚像元的相关性（在COSI-Corr软件中进行）。

对于深度，可以使用流动方程来模拟：

                                                                                                                      ①

其中，和表示分别表示表面流速和基础流速，设定为的25%；是冰川厚度；是爬行参数；n设为3。

我们为了推导深度H，需要用到基础应力模型：

                                                                                                                               ②

其中，表示冰密度；表示重力加速度；是形状因子；是坡度，从DEM数据获得。

联立上述的两个方程，我们可以得到冰层的厚度H为：

                                                                                                                    ③

该方程的不确定性可以推导为：

对于RGI数据，只用到了大约5%的范围，它的不确定性大约为12.5%。

冰层密度的不确定性大约为10%，方差为830~923kgm^(-3)。

sina的不确定性为8.7%

冰体积的总体不确定性为18%。

最后利用python里面的tkinter模块开发GUI，实现该功能。

4. 实验结果

具体结果就不多介绍了，作者制作了一些图，关于冰的体积估算图：

湖泊的过侵蚀图：

还有湖泊体积估算图：

作者最后量化了各个参数及其误差。

四、小结

这篇文章挺普通的，方法上作者并没有自己的任何创新，连核心算法都是用的别人的，关键工作在于自己的分析，不过做的图还挺好看。