Nepal Earthquake Triggered Landslides Recognition and Deformation Analysis of Avalanches' Region
-
摘要: 2015-04-25尼泊尔Gorkha附近发生了Mw 7.8级强震，引发尼泊尔及中国西藏地区大量的山体滑坡和雪崩灾害，造成了严重的人员伤亡与经济损失。基于Landsat 8影像和亚像素相关性匹配技术识别形变失相关区域，并结合目视解译提取震害滑坡与雪崩沉积区域，统计分析了所提取的滑坡数目、面积等与地形因子的关系，同时对珠峰雪崩高发区的冰川形变进行了定量分析。实验结果表明，提取形变场中的形变失相关区域并结合目视解译验证是一种识别震害滑坡与雪崩沉积区域的有效方法；统计结果显示滑坡体随坡度分布呈现高斯特性，多发生在坡度大于30°的地形，易发于河谷两侧，且高坡度更容易触发单体面积大的滑坡体。该研究为利用Landsat 8影像数据进行震害滑坡与雪崩沉积区域提取、雪崩高发区的冰川流速监测预警工作提供了新视角。
-
关键词:
- 尼泊尔Gorkha地震 /
- 亚像素相关性匹配 /
- 滑坡 /
- 雪崩 /
- 冰川流速
Abstract: On April 25th 2015, an Mw 7.8 earthquake occurred near Gorkha, Nepal, which triggered substantial landslides and avalanches in Tibet and Nepal causing severe casualties and economic losses. In this study, based on Landsat 8 images and sub-pixel correlation technology, we acquire horizontal deformation fields of the earthquake region. Combined with the visual interpretation method, we detected the landslides and debris deposition areas of avalanches by extracting the corresponding decorrelation noises. We determined the correlation between the landslides and topography factors. Experimental results show that the extraction of decorrelation noises caused by geohazards from the deformation field, when combined with validation of visual interpretation method, is an effective method for locating the regions effected by landslides and avalanches. Furthermore, results show that landslides are more likely to occur at slope > 30° terrain, both sides of the valley, and high slope is more easily trigger larger single landslides. This study provides a new perspective on detection of earthquake-triggered landslides and debris deposition areas of avalanches, monitoring and early warning of glaciers' velocity based on the Landsat 8 satellite imagery.-
Keywords:
- Nepal Gorkha earthquake /
- sub-pixel correlation /
- landslides /
- avalanches /
- glacier velocity
-
作为最为广泛使用的电离层经验模型之一，国际参考电离层模型（international reference ionosphere, IRI）对电离层研究和应用具有十分重要的意义[1]。IRI模型利用了垂测仪、非相干散射雷达、卫星观测及探空火箭数据，可以输出任意指定时间和经纬度的60~2 000 km的电子密度和总电子含量（total electron content, TEC）等多种电离层参数信息[2]。自1972年以来，随着电离层观测数据的积累、观测技术的革新和建模方法的改进，IRI模型已发布了多个版本，目前最新版本为IRI-2016 [3-4]。
峰值密度NmF2和峰值高度hmF2是表征电子密度廓线形态的重要参量。对于NmF2的计算，IRI-2016模型中包含IRI_CCIR和IRI_URSI两种模型选项。对于hmF2的计算，在旧版本模型使用的IRI_CCIR hmF2模型的基础上，IRI-2016模型中另外新增了IRI_AMTB和IRI_Shubin两种模型选项。hmF2的IRI_CCIR模型是基于1954—1958年全球分布的约150个垂测仪观测资料，利用hmF2与参数M(3 000)F2的强相关关系由M(3 000)F2计算得到。其中，M(3 000)F2为3 000 km高度处最大可用频率和F2层临界频率（foF2）的比值。这种间接计算hmF2的方法有诸多缺陷。IRI-2016模型中新增的两种计算hmF2的模型选项直接对hmF2进行建模，不依赖M(3 000)F2。其中，IRI_AMTB模型是根据全球分布的26个垂测仪（其中包括中国区域的武汉站、安阳站和三亚站）在1998—2006年的hmF2观测值，纯粹利用地基垂测仪观测值进行建模[5]；IRI_Shubin模型的建立则使用了挑战性小卫星有效载荷（Challenging Minisatellite Payload，CHAMP）、重力恢复与气候实验（Gravity Recovery and Climate Experiment，GRACE）和气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology, Ionosphere and Climate，COSMIC)3个掩星任务的观测资料以及分布全球的62个垂测仪于1987—2012年的观测数据[6]。
COSMIC由美国国防部和中国台湾地区联合研制并于2006-04-15发射升空，在掩星观测任务中具有里程碑的意义[7]。COSMIC由6颗低轨卫星组成，轨道高度在800 km左右，轨道倾角为72°，并配备能够接收GPS信号的接收机。掩星事件发生时，GPS信号穿过电离层时受电子密度影响而发生折射，产生附加相位延迟。在信号直线传播和球对称假设等前提下，由传播路径上的总电子含量使用Abel积分变换方法可得到电子密度廓线。COSMIC掩星反演结果已经得到了多种独立观测手段的验证[8-10]。
中国区域（70°E~140°E，15°N~55°N）跨越中低纬地区，电离层变化幅度较大。文献[11]利用中国区域的4个垂测仪台站对IRI-2016模型在中国区域的性能进行了评估，但只能对垂测仪台站周围的电离层特征进行验证分析。COSMIC掩星观测具有全球覆盖、高精度、垂直分辨率高等特点，可用于对IRI-2016模型在整个中国区域的表现进行分析。本文的研究内容包括使用2008—2014年的COSMIC掩星数据对地磁平静条件下IRI-2016模型在中国区域的性能进行验证，并对IRI-2016模型在中国区域使用时不同NmF2和hmF2模型选项的输出结果进行统计分析，以确定IRI-2016模型在中国区域输出值的最优选项。
1 数据与模型
1.1 COSMIC掩星数据
COSMIC掩星二级产品“ionPrf”数据可从COSMIC数据存档与分析中心（COSMIC Data Analysis and Archive Center, CDAAC）的网站http://cdaac-www.cosmic.ucar.edu下载得到。该数据包含了掩星事件发生的时间、位置、NmF2、hmF2和电子密度廓线等信息。在分析地磁平静条件下IRI模型的表现时，需要先剔除3 h ap指数大于15的时间段的COSMIC掩星数据[8]。另一方面，尽管CDAAC在反演电子密度廓线时采取了一定的质量检核措施，但通过其检核的产品中仍然存在问题廓线[8]。文献[12]提出根据顶部电离层电子密度廓线的斜率和电子密度波动的平均偏差（mean deviation, MD）来剔除这些问题廓线。本文选择COSMIC电子密度廓线的标准为在420~490 km的斜率小于-0.02；平均偏差不超过1.5。此外，峰值高度hmF2范围应在200~500 km之间。
1.2 IRI-2016模型
IRI-2016模型输出值可从http://irimodel.org/网站获取。该模型中计算NmF2的选项包括IRI_CCIR模型和IRI_URSI模型。其中，IRI_URSI模型为推荐选项。计算峰值高度hmF2的选项包括原有的IRI_CCIR模型以及新增的IRI_AMTB模型和IRI_Shubin模型。其中，IRI_AMTB模型为推荐选项。模型的输入参数包括经纬度和时间。由于掩星观测具有水平漂移，本文选用COSMIC掩星观测的F2层峰值高度处的经纬度和发生时刻作为输入量，通过控制模型选项得到不同选项下输出的NmF2和hmF2。图 1(a)、1(b)分别为当地时间（local time, LT）2018-02-03T13:35和2008-02-08T23:46的COSMIC掩星电子密度廓线和IRI-2016模型不同选项下输出结果的比较图。可以看出，COSMIC电子密度廓线和IRI-2016模型不同选项以及选项之间的输出结果存在着明显差异。
![]()
图 1 COSMIC反演的电子密度廓线和IRI-2016模型不同选项输出值的比较Figure 1. Comparison Among Electron Density Profiles Derived from COSMIC RO and IRI-2016 Model with Different Model Options1.3 F10.7指数
F10.7指数是指10.7 cm（2 800 MHz）波长的太阳辐射通量，单位为sfu，与太阳黑子数之间有强相关性，能很好地表征太阳活动强度。图 2是F10.7指数的月均值在第24太阳活动周2008—2014年期间的变化。可以看出，2008—2010年为太阳活动低年，尤其是2008、2009年，太阳活动达到了极低的水平；之后太阳活动水平逐渐提高，在2011—2014年太阳活动维持在较高水平，为太阳活动高年。
2 结果分析
2.1 IRI-2016模型在中国区域的性能分析
为了了解IRI-2016模型在中国区域的表现，首先采用模型推荐的选项计算NmF2和hmF2，将结果与COSMIC掩星观测进行比较。图 3为2008年和2014年COSMIC掩星观测与IRI-2016模型输出的峰值密度NmF2与峰值高度hmF2的散点图。其中，红色实线代表COSMIC观测值等于IRI预测值，黑色虚线代表使用最小二乘方法得到的拟合线，R和S分别表示两种方法得到的相关系数和拟合线斜率。由图 3(a)和3(b)可以看出，尽管在太阳活动低年2008年和太阳活动高年2014年，COSMIC掩星观测和IRI模型输出的NmF2相关系数非常接近（分别为0.86和0.87），但在2014年数据对较为均匀地围绕红色实线分布，且拟合线和等值线非常接近，斜率高达0.90。而在2008年，数据对较为分散，且较多的数据对分布在等值线下侧，拟合直线与等值线之间存在较大偏差。拟合直线的斜率（0.70）表明在太阳活动低年2008年，相较于COSMIC掩星观测值，IRI-2016模型明显高估了NmF2。与峰值密度NmF2情形类似，图 3(c)和图 3(d)表明, 在2014年IRI-2016模型峰值高度hmF2输出值与COSMIC掩星观测值更为接近，拟合直线几乎与等值线重合。而在2008年，拟合直线的斜率为0.76，远小于2014年的斜率0.97，说明2008年IRI-2016模型的hmF2输出值相对于COSMIC掩星观测存在明显高估。同时，由相关系数可以看出，无论在2008年还是2014年，NmF2的相关系数均高于hmF2，表明IRI-2016模型和COSMIC掩星观测输出的NmF2一致性优于hmF2。
![]()
图 3 COSMIC掩星观测与IRI-2016模型采用推荐选项输出的NmF2与hmF2散点图Figure 3. Scatter Plots of NmF2 and hmF2 from COSMIC RO and IRI-2016 Model with Recommended Model Options为了进一步定量分析太阳活动低年和高年COSMIC掩星观测值和IRI-2016模型输出值的差异，定义了相对中误差绝对值的均值（mean absolute relative error，MARE）M，计算公式为：
$$ M=\frac{\mathop{\sum }^{}\left( \left| I-C \right|/C \right)}{N} $$ 式中，C表示由掩星数据反演的NmF2或者hmF2；I表示对应的由模型计算的NmF2或者hmF2；N为观测数据个数。
表 1给出了2008—2014年的NmF2和hmF2的MARE数值。可以看出，无论是NmF2还是hmF2，太阳活动低年（2008—2010年）的MARE均高于太阳活动高年（2011—2014年），在太阳活动极小年（2008—2009年），IRI-2016模型输出值与COSMIC掩星观测值的差异最大。
表 1 2008—2014年NmF2和hmF2的MARE/%Table 1. MARE of NmF2 and hmF2 During 2008―2014/%统计项 2008年 2009年 2010年 2011年 2012年 2013年 2014年 NmF2 31.112 2 31.481 8 27.836 7 26.7882 21.888 1 24.443 1 25.204 4 hmF2 11.689 0 11.944 2 9.747 9 8.460 8 7.071 9 8.021 6 7.344 6 因此，总体而言，对中国区域，无论是NmF2还是hmF2，以COSMIC掩星数据为参照，IRI-2016模型在太阳活动高年的表现都优于太阳活动低年。
2.2 太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布
为了进一步分析中国区域COSMIC掩星反演的峰值参数与IRI-2016模型输出值间的差异，对比了差异更为明显的太阳活动低年（2008—2010年）峰值参数的空间分布。考虑到电离层活跃水平，分析的时间段选为最具代表性的春秋季的中午时分（LT 12:00—14:00）。在本文中，季节的划分标准为：3月、4月、9月、10月为春秋季，5月—8月为夏季，11月、12月、1月和2月为冬季。图 4为2008—2010年春秋季中午（LT 12:00—14:00）由COSMIC掩星与IRI-2016模型得到的峰值参数NmF2与hmF2的均值在中国区域的空间分布。图 4(a)~4(d)显示，COSMIC掩星观测与IRI模型输出的峰值参数空间分布形态均反映出低纬地区电离层更为活跃，同时观测值与模型输出值之间存在明显差异。图 4(e)、4(g) 表明，对于NmF2，与COSMIC掩星反演结果对比，IRI模型输出值在15°N~30°N的低纬区域明显高估，相对偏差最高可达100%以上；而在30°N~55°N的中纬区域，IRI模型输出的NmF2普遍表现为低估。图 4(f)、4(h)表明，对于hmF2，情况恰恰相反：在30°N~55°N中纬区域，IRI模型的hmF2输出值普遍高于COSMIC观测值，而在15°N~30°N低纬区域，模型输出值普遍低于COSMIC观测值。由图 4（f）可见，在太阳活动低年，IRI-2016模型与COSMIC在中国区域得到的hmF2的绝对偏差一般不超过50 km。
![]()
图 4 2008—2010年春秋季LT 12：00—14：00 COSMIC掩星与IRI-2016模型的NmF2与hmF2在中国区域的空间分布Figure 4. Spatial Distributions of NmF2 and hmF2 over China Region at LT 12:00—14:00 During Equinox Months of 2008—2010 Derived from COSMIC RO and IRI-2016 Model2.3 IRI-2016模型峰值参数不同选项输出结果的比较分析
以上分析中所使用的IRI-2016模型输出值均基于模型的推荐选项得到。除了推荐选项，IRI-2016对NmF2和hmF2的计算均提供其他模型选项，而使用不同选项的模型输出结果存在差异。为了了解不同选项的模型在中国区域的表现，进一步以COSMIC掩星数据为参照，对IRI-2016模型中两种计算NmF2的选项（IRI_CCIR和IRI_URSI模型选项）和3种计算hmF2的选项（IRI_CCIR、IRI_AMTB和IRI_Shubin模型选项）的输出结果进行了比较。图 5给出了不同季节COSMIC掩星与IRI-2016模型两种不同选项得到的NmF2均值的日变化情况，其中图 5(a)~5(c)和图 5(d)~5(f)分别为太阳活动低年和高年的结果。统计中如果掩星事件发生时间与某一整点时刻相差不超过30 min，则该时刻的掩星数据和IRI模型输出结果被用于计算该整点时刻的平均值。可以看出，无论采用哪个选项，NmF2在约LT 04：00达到最小值，之后快速升高，在约LT 14：00达到全天最大值，之后逐渐下降。春秋季NmF2高于其他季节，电离层最为活跃。基于两种不同选项的IRI-2016模型输出值与COSMIC观测值均有相似的日变化趋势，但分别基于两种选项的模型输出结果之间存在一定差异。图 5(a)~5(c)表明采用两种选项的IRI-2016模型输出结果非常接近，但在中午时分，IRI_CCIR模型和IRI_URSI模型均高估了峰值密度，其中IRI_URSI的高估程度更严重。在春秋季和冬季，IRI_CCIR模型和IRI_URSI模型在LT 17：00—05：00则低估了峰值密度，其中IRI_URSI低估程度更严重。图 5(d)~5(f)与图 5(a)~5(c)的对比表明，随着太阳活动水平的增加，NmF2随之增大，但中午时分采用两种选项的IRI-2016模型输出值与COSMIC观测值之间的差异均有减小。特别是在太阳活动高年的夏季，IRI-2016模型和COSMIC观测值全天都符合得很好。而在春秋季LT 16：00—04：00，两种选项的IRI-2016模型输出值均小于COSMIC观测值，其中IRI_URSI的偏离度相对更大。
![]()
图 5 不同季节COSMIC掩星与IRI-2016模型两种不同选项下NmF2的日变化图Figure 5. Diurnal Variations of NmF2 from COSMIC RO and IRI-2016 Model Using Two Different Options for Different Seasons图 6给出了不同季节COSMIC掩星观测与IRI-2016模型3种不同选项输出的hmF2日变化情况。可以看出，3种选项对应的IRI-2016模型输出与COSMIC观测的hmF2呈现相似的日变化趋势，即在LT 06：00―07：00达到最小值，在LT 12：00达到小高峰，并在约LT 24：00LT达到全天最大值。在太阳活动低年，3种计算hmF2的模型选项中，IRI_Shubin模型输出值和COSMIC观测值最为接近，其他两种模型输出值相对于COSMIC观测值均存在高估。IRI_AMTB模型输出值相对于COSMIC观测值的高估现象更为突出，冬季最显著。与太阳活动低年相比较，在太阳活动高年，3种选项对应的hmF2模型输出值与COSMIC观测值的一致性均相对更好。但IRI_AMTB模型输出值在夜间的偏差仍较大。
![]()
图 6 不同季节COSMIC掩星与IRI-2016模型不同选项下hmF2的日变化图Figure 6. Diurnal Variations of hmF2 from COSMIC RO and IRI-2016 Model Using Different Options for Different Seasons3 讨论
本文以COSMIC掩星数据为参照，对IRI-2016模型在中国区域的电离层峰值参数NmF2和hmF2输出值随太阳活动强度、季节和昼夜的变化进行了对比分析。整体而言，模型输出值与COSMIC掩星观测值的相关性较高，但太阳活动低年的相关性低于太阳活动高年。两组数据拟合直线的斜率表明，IRI-2016模型输出值与COSMIC掩星观测值的一致性在太阳活动低年低于太阳活动高年。
文献[13]使用位于磁赤道附近的Jicamarca站垂测仪数据与IRI模型输出值进行验证分析，结果发现模型输出值与垂测仪观测值在太阳活动高年的相关性高于太阳活动低年，这与本文的研究结果一致。这主要是由于太阳活动低年，特别是2008年、2009年为太阳活动非寻常极小年，太阳活动水平极低且持续时间较长，IRI模型输出值会产生较大误差，甚至出现不合理的输出结果[3]。
太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布比较情况显示，虽然两者的空间分布形态基本一致，但仍存在明显差异。在15°N~30°N范围内，IRI-2016模型往往高估NmF2而低估hmF2。在30°N~55°N范围内，IRI-2016模型则存在低估NmF2而高估hmF2的现象。COSMIC掩星与IRI-2016模型之间的差异一方面是由于IRI-2016模型本身的误差，另一方面则是由于COSMIC掩星反演结果存在的误差。在球对称假设的前提下，COSMIC掩星反演的电子密度实际上是掩星切点附近整个掩星区域的平均值。因此，在电子密度梯度较大的区域，COSMIC掩星反演的电子密度与真实大气状态相比可能存在较大差异[14]。中国区域低纬度地区位于赤道电离层异常的驼峰区域附近，电子密度梯度很大，导致COSMIC反演的电子密度廓线相较于实际值可能存在较大偏差。
文献[15]使用中国区域位于同一经度带上的4个垂测仪台站（漠河站、北京站、武汉站和三亚站）数据对COSMIC掩星峰值参数进行了比较验证，发现在漠河站两者符合程度较高，但在武汉站和北京站中午时分COSMIC NmF2观测值高于垂测仪观测值，在三亚站中午时分COSMIC NmF2观测值低于垂测仪观测值。反演存在的误差可能使得在中纬度的武汉站（31.0°N，114.6°E）COSMIC反演的NmF2高于IRI-2016模型值，在低纬度的三亚站（18.3°N，109.6°E）COSMIC反演的NmF2低于IRI-2016模型值。文献[11]使用这4个垂测仪台站数据对IRI-2016模型计算hmF2的不同选项输出值结果进行比较分析，发现IRI-2016模型推荐的IRI_AMTB模型选项在4个台站均高估hmF2。本文结果也表明IRI-2016模型的hmF2输出值在30°N~55°N高于COSMIC掩星观测值。对于15°N~30°N范围，IRI-2016模型的hmF2输出值相对于掩星资料存在低估，原因可能是由于该范围位于赤道电离层异常驼峰附近，掩星反演的hmF2高于实际值[11]。
4 结语
本文对2008—2014年IRI-2016模型在中国区域的电离层峰值密度NmF2和峰值高度hmF2输出值与COSMIC掩星数据进行了对比，并对IRI-2016模型中计算这两个峰值参数的不同选项进行了比较分析。结论如下：
1）IRI-2016模型输出值与COSMIC掩星观测数据的一致性在太阳活动高年比太阳活动低年更好。模型NmF2输出值与COSMIC掩星观测值的相关性高于hmF2输出值与COSMIC掩星观测值的相关性。
2）在太阳活动低年春秋季的中午时段，与COSMIC掩星数据相比，IRI-2016模型在30°N~55°N对NmF2和hmF2分别存在低估和高估现象；在15°N~30°N区域内则恰恰相反。
3）以COSMIC掩星观测值为参照，IRI_CCIR模型和IRI_URSI模型均对NmF2在中午时分存在高估现象，该现象在太阳活动低年更为显著。而在春秋季的LT 17：00—04：00时间段，采用这两个模型选项得到的IRI-2016模型输出值均低于COSMIC观测值。
4）在太阳活动低年，无论任何季节，IRI_CCIR模型和IRI_AMTB模型的hmF2输出值均高于COSMIC掩星观测值，其中IRI_AMTB模型输出值相对于掩星观测的偏差更大，尤其是在冬季。IRI_Shubin模型的hmF2输出值与COSMIC观测值最为接近，即使在太阳活动高年，两者的偏差也小于其他两种模型。这表明COSMIC掩星资料在IRI_Shubin模型改进hmF2输出值中发挥了重要作用。
5）根据本文的研究结果，在中国区域采用IRI-2016计算NmF2和hmF2时，分别推荐使用IRI_CCIR和IRI_Shubin模型，而非官方推荐的IRI_URSI和IRI_AMTB模型。
-
![]()
图 1 震区构造背景与Landsat 8影像覆盖图
Figure 1. Background Setting and Coverage of Landsat 8 for the 2015 Nepal Earthquake
![]()
图 5 吉隆、聂拉木地区不同坡度水平滑坡体数目、总面积与平均面积分布直方图
Figure 5. Landslides' Quantity, Gross and Average Area Distribution Histogram of Jilong and Nyalam in Different Slope Levels
![]()
图 6 郎塘形变场和雪崩沉积区提取图
Figure 6. Extraction of Deformation Fields and Avalanches' Debris Areas in Langtang
表 1 本文实验所采用的Landsat 8影像信息
Table 1 Information of Landsat 8 Images Used in This Study
影像对 时段 轨道号 分幅号 获取时间 云量/% 影像对1 震前 141 41 2014-03-26 2.4 震后 141 41 2015-04-30 30.2 影像对2 震前 141 40 2014-01-05 4.4 震后 141 40 2015-04-30 30.2 影像对3 震前 140 41 2015-03-22 11.8 震后 140 41 2015-05-25 25.6 影像对4 震前 140 41 2014-03-19 17.4 震后 140 41 2015-05-22 13.1 
下载: 导出CSV
表 2 吉隆山体滑坡编目表
Table 2 Landslides Catalog of Jilong
序号 滑坡位置 滑坡体面积/m2 技术监测手段 高程/m 坡向/(°) 坡度/(°) 东经/(°) 北纬/(°) SPC 目视判读 J01 85.363 28.360 163 646 √ √ 2 615 289.6 52.8 J02 85.359 28.354 29 052 √ √ 2 370 268.7 43.5 J03 85.337 28.376 37 580 √ √ 2 490 80.4 30.2 J04 85.324 28.389 49 045 √ √ 2 737 312.1 45.9 J05 85.305 28.401 5 599 √ √ 2 837 359.1 37.3 J06 85.304 28.406 33 891 √ √ 2 809 28.9 40.3 J07 85.350 28.364 2 603 √ √ 2 543 10.6 47.5 J08 85.345 28.369 11 552 √ √ 2 415 165.9 38.9 J09 85.337 28.374 5 967 √ √ 2 450 102.1 33.3 J10 85.340 28.381 7 431 线状滑坡× √ 2 839 208.5 53.9 J11 85.337 28.374 1 295 小滑坡× √ 2 603 24.7 27.1 J12 85.355 28.364 680 小滑坡× √ 2 472 120.0 45.7 J13 85.370 28.376 6 789 √ √ 2 925 212.4 44.8 J14 85.358 28.370 42 480 √ √ 2 439 196.8 52.6 J15 85.355 28.372 1 057 小滑坡× √ 2 550 172.0 41.4 J16 85.308 28.404 2 263 小滑坡× √ 2 804 191.8 33.3 J17 85.303 28.408 1 728 小滑坡× √ 2 805 49.1 42.9 J18 85.312 28.399 3 628 小滑坡× √ 2 741 172.1 25.6 J19 85.358 28.411 351 小滑坡× √ 2 984 135.1 37.0 J20 85.364 28.382 21 275 √ √ 3 152 256.7 48.2 J21 85.365 28.382 17 310 √ √ 3 152 256.7 42.7 注：√表示可以监测到对应滑坡体，×表示无法监测到对应滑坡体。
下载: 导出CSV
表 3 聂拉木山体滑坡编目表
Table 3 Landslides Catalog of Nyalam
序号 滑坡位置 滑坡体面积/m2 高程/m 坡向/(°) 坡度/(°) 东经/(°) 北纬/(°) N01 85.997 28.090 46 565 3 457 123.3 37.6 N02 85.998 28.087 7 153 3 264 177.7 31.7 N03 85.990 28.084 13 099 3 769 166.3 60.8 N04 85.997 28.085 6 360 3 238 270 26.9 N05 85.991 28.046 38 639 3 265 274.9 57.1 N06 85.994 28.039 53 518 3 049 171.8 32.8 N07 85.981 28.032 51 028 3 157 100.9 68.8 N08 85.980 28.021 12 191 2 602 106.2 45.4 N09 85.963 27.983 309 600 2 375 157.2 60.8 N10 86.039 27.975 381 165 3 462 106.5 42.2
下载: 导出CSV
-
[1] Wikipedia. April 2015 Nepal Earthquake[OL]. https://en.wikipedia.org/wiki/April 2015 Nepal earthquake, 2016
[2] Feng G C, Li Z W, Shan X J, et al. Geodetic Model of the 2015 April 25 Mw 7.8 Gorkha Nepal Earthquake and Mw 7.3 Aftershock Estimated from InSAR and GPS Data[J]. Geophysical Journal International, 2015, 203(2):896-900 doi: 10.1093/gji/ggv335
[3] 陶舒. 汶川地震滑坡遥感信息提取及灾害危险性评价研究[D]. 北京: 首都师范大学, 2009 Tao Shu. Information Extraction of Landslides Induced by the Wenchuan Earthquake and Susceptibi-lity Assessment Using Remotely Sensed Data[D]. Beijing: Capital Normal University, 2009
[4] 许冲, 徐锡伟, 沈玲玲, 等. 2014年鲁甸MS 6.5地震触发滑坡编录及其对一些地震参数的指示[J].地震地质, 2014, 36(4):1186-1203 https://www.researchgate.net/profile/Chong_Xu/publication/287556048_Inventory_of_landslides_triggered_by_the_2014_MS65_Ludian_earthquake_and_its_implications_on_several_earthquake_parameters/links/580ea47a08ae7525273d29d3.pdf?origin=publication_list Xu Chong, Xu Xiwei, Shen Lingling, et al. Inventory of Landslides Triggered by the 2014 MS 6.5 Ludian Earthquake and Its Implications on Several Earthquake Parameters[J].Seismology and Geology, 2014, 36(4):1186-1203 https://www.researchgate.net/profile/Chong_Xu/publication/287556048_Inventory_of_landslides_triggered_by_the_2014_MS65_Ludian_earthquake_and_its_implications_on_several_earthquake_parameters/links/580ea47a08ae7525273d29d3.pdf?origin=publication_list
[5] 武雪玲, 沈少青, 牛瑞卿. GIS支持下应用PSO-SVM模型预测滑坡易发性[J].武汉大学学报·信息科学版, 2016, 41(5):665-671 http://ch.whu.edu.cn/CN/abstract/abstract5447.shtml Wu Xueling, Shen Shaoqing, Niu Ruiqing. Landslide Susceptibility Prediction Using GIS and PSO-SVM[J]. Geomatics and Information Science of Wuhan University, 2016, 41(5):665-671 http://ch.whu.edu.cn/CN/abstract/abstract5447.shtml
[6] Lacroix P, Zavala B, Berthier E, et al. Supervised Method of Landslide Inventory Using Panchromatic SPOT5 Images and Application to the Earthquake-Triggered Landslides of Pisco (Peru, 2007, Mw8.0)[J]. Remote Sensing, 2013, 5(6):2590-2616 doi: 10.3390/rs5062590
[7] Leprince S, Barbot S, Ayoub F, et al. Automatic and Precise Orthorectification, Co-registration and Sub-pixel Correlation of Satellite Images, Application to Ground Deformation Measurements[J]. IEEE Transactions on Geoscience and Remote Sen-sing, 2007, 45(6):1529-1558 doi: 10.1109/TGRS.2006.888937
[8] 董晓燕, 丁晓利, 李志伟, 等.一种新的SAR像素偏移量估计流程及其在同震形变监测中的应用[J].武汉大学学报·信息科学版, 2011, 36(7):789-792 http://ch.whu.edu.cn/CN/abstract/abstract592.shtml Dong Xiaoyan, Ding Xiaoli, Li Zhiwei, et al. A New Chain for SAR Pixel Offset Eestimation and Its Application in Co-seismic Deformation Measurement[J].Geomatics and Information Science of Wuhan University, 2011, 36(7):789-792 http://ch.whu.edu.cn/CN/abstract/abstract592.shtml
[9] Ding C, Feng G C, Li Z W, et al. Spatio-Temporal Error Sources Analysis and Accuracy Improvement in Landsat 8 Image Ground Displacement Measurements[J].Remote Sensing, 2016, 8(11):937-958 doi: 10.3390/rs8110937
[10] Binet R, Bollinger L. Horizontal Coseismic Deformation of the 2003 Bam (Iran) Earthquake Measured from SPOT-5 THR Satellite Imagery[J]. Geophysi-cal Research Letters, 2005, 320(2):287-294 doi: 10.1029/2004GL021897
[11] Leung T K P, Admas J. The 2005, Mw 7.6 Kashmir Earthquake, Rupture Kinematics from Sub-pixel Correlation of ASTER Images and Seismic Waveforms Analysis[J]. Earth and Planetary Science Letters, 2006, 249(3):514-528 https://thesis.library.caltech.edu/2257/
[12] 冯光财, 许兵, 单新建, 等.基于Landsat 8光学影像的巴基斯坦Awaran Mw7.7地震形变监测及参数反演研究[J].地球物理学报, 2015, 58(5):1634-1644 doi: 10.6038/cjg20150515 Feng Guangcai, Xu Bing, Shan Xinjian, et al. Coseismic Deformation and Source Parameters of the 24 September 2013 Awaran, Pakistan Mw7.7 Earthquake Derived from Optical Landsat 8 Satellite Images[J]. Chinese Journal Geophysics, 2015, 58(5):1634-1644 doi: 10.6038/cjg20150515
[13] Avouac J P, Meng L, Wei S, et al. Lower Edge of Locked Main Himalayan Thrust Unzipped by the 2015 Gorkha Earthquake[J].Nature Geoscience, 2015, 8(9):1-6 http://palgrave.nature.com/ngeo/journal/v8/n9/abs/ngeo2518.html
[14] Kargel J S, Leonard G J, Shugar D H, et al. Geomorphic and Geologic Controls of Geohazards Induced by Nepal's 2015 Gorkha Earthquake[J]. Science, 2016, 351(6269):1-12 http://ecite.utas.edu.au/107060
[15] Jane Q. Killer Landslides: The Lasting Legacy of Nepal's Quake[OL]. http://www.nature.com/news/killer-landslides-the-lasting-legacy-of-nepal-s-quake-1.19803, 2016
[16] Scherler D, Leprince S, Strecker M R. Glacier-surface Velocities in Alpine Terrain from Optical Satellite Imagery-Accuracy Improvement and Quality Assessment[J]. Remote Sensing of Environment, 2008, 112(10):3806-3819 doi: 10.1016/j.rse.2008.05.018
-
期刊类型引用(2)
1. 陆浩楠，叶世榕，张庆兰，夏朋飞，鄂盛龙. 利用COSMIC-1/2 IonPrf产品分析中国区域近15年电离层时空变化特征. 武汉大学学报(信息科学版). 2024(05): 765-774 . 
百度学术
2. 唐丝语，周春红，黄智，袁洪. 国际参考电离层模型输出参数峰值高度性能. 空间科学学报. 2021(06): 898-904 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 1602
- HTML全文浏览量: 184
- PDF下载量: 332
- 被引次数: 4
