高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

基于多时相Landsat数据的汶川地震后王家岩滑坡植被恢复估计

, 

上一篇

下一篇

叶席, 匡鸿海. 基于多时相Landsat数据的汶川地震后王家岩滑坡植被恢复估计[J]. 西南大学学报（自然科学版）, 2022, 44(6): 150-159. doi: 10.13718/j.cnki.xdzk.2022.06.016
引用本文: 叶席, 匡鸿海. 基于多时相Landsat数据的汶川地震后王家岩滑坡植被恢复估计[J]. 西南大学学报（自然科学版）, 2022, 44(6): 150-159. doi: 10.13718/j.cnki.xdzk.2022.06.016
shu
YE Xi, KUANG Honghai. Estimation of Vegetation Restoration in Wangjiayan Landslide after Wenchuan Earthquake Based on Multi-temporal Landsat Data[J]. Journal of Southwest University Natural Science Edition, 2022, 44(6): 150-159. doi: 10.13718/j.cnki.xdzk.2022.06.016
Citation: YE Xi, KUANG Honghai. Estimation of Vegetation Restoration in Wangjiayan Landslide after Wenchuan Earthquake Based on Multi-temporal Landsat Data[J]. Journal of Southwest University Natural Science Edition, 2022, 44(6): 150-159. doi: 10.13718/j.cnki.xdzk.2022.06.016
shu

基于多时相Landsat数据的汶川地震后王家岩滑坡植被恢复估计

  • 西南大学 地理科学学院，重庆 400715

  • 基金项目: 国家自然科学基金项目(41877369)
详细信息
    作者简介:

    叶席，硕士研究生，主要从事GIS技术分析与应用的研究 .

    通讯作者: 匡鸿海，副教授

  • 中图分类号: P315

Estimation of Vegetation Restoration in Wangjiayan Landslide after Wenchuan Earthquake Based on Multi-temporal Landsat Data

  • School of Geographical Sciences, Southwest University, Chongqing 400715, China

  • 摘要: 为研究汶川地震后发生严重同震滑坡的北川县王家岩受损植被恢复情况，利用遥感技术对王家岩震后地区进行长期植被动态监测，基于研究区2007-2020年多时相Landsat遥感影像，采用归一化植被指数(NDVI)和植被覆盖恢复率(VRR)来定量评价植被生长和恢复状态. 结果表明：研究区内滑坡区域受损植被已逐渐恢复，其NDVI由地震发生后的0.385上升至2020年8月的0.624；经过12年的演替，滑坡区域植被覆盖恢复率经季节性调整后接近91.8%；根据VRR回归模型预估，滑坡区域内的植被需要26年左右才能恢复到地震前的水平.
    Abstract: In order to study the restoration of damaged vegetation in Wangjiayan, Beichuan County, where serious co-seismic landslides occurred after the Wenchuan earthquake, remote sensing technology was used to monitor the long-term vegetation dynamics in the post-earthquake area of this region. Based on the multi-temporal Landsat remote sensing images of the study area from 2007 to 2020, the normalized difference vegetation index (NDVI) and vegetation coverage recovery rate (VRR) were used to quantitatively evaluate the vegetation growth and recovery status. The results indicated that the vegetation restoration in the studied landslide area showed an upward trend, and NDVI increased from 0.385 after the earthquake to 0.624 in August 2020. After 12 years of succession, the vegetation coverage recovery rate in landslide area was close to 91.8% after seasonal adjustment. According to the VRR regression model, it is estimated that the vegetation restoration in the landslide area needs to take about 26 years to reach the level before the earthquake.
  • 加载中

  • 图 1  研究区地理位置

    下载: 全尺寸图片 幻灯片

    图 2  研究区假彩色卫星图像

    下载: 全尺寸图片 幻灯片

    图 3  研究区滑坡提取结果

    下载: 全尺寸图片 幻灯片

    图 4  研究区NDVI提取结果

    下载: 全尺寸图片 幻灯片

    图 5  研究区(a)和滑坡区(b)VRR回归模型

    下载: 全尺寸图片 幻灯片

    表 1  研究区、非滑坡区以及滑坡区NDVI

    成像时间 Landsat 距地震时间/年 NDVI
    研究区 非滑坡区域 滑坡区域
    2007/5/6 5 -1.02 0.654 0.646 0.732
    2008/7/11 5 0.16 0.553 0.572 0.385
    2009/6/12 5 1.08 0.582 0.601 0.424
    2010/3/27 5 1.88 0.369 0.378 0.281
    2011/5/17 5 3.01 0.563 0.573 0.481
    2013/5/22 8 5.03 0.632 0.636 0.522
    2014/7/12 8 6.16 0.631 0.641 0.563
    2015/10/19 8 7.43 0.505 0.518 0.477
    2016/4/28 8 7.96 0.613 0.628 0.591
    2017/5/1 8 8.97 0.627 0.639 0.550
    2018/6/5 8 10.06 0.671 0.676 0.607
    2019/8/11 8 11.16 0.668 0.678 0.618
    2020/8/13 8 12.17 0.671 0.675 0.624
    下载: 导出CSV

    表 2  研究区和滑坡区季节性调整前后的VRR

    成像时间 NDVI VRR NDVI
    偏移值    		 调整后NDVI 调整后VRR
    研究区 滑坡区 研究区 滑坡区 研究区 滑坡区 研究区 滑坡区
    2007/5/6 0.654 0.732 0.031 0.685 0.764
    2008/7/11 0.553 0.385 0 0 -0.008 0.545 0.377 0 0
    2009/6/12 0.582 0.424 0.278 0.303 -0.009 0.574 0.415 0.349 0.351
    2010/3/27 0.369 0.281 0.051 0.221 0.087 0.455 0.368 0.417 0.446
    2011/5/17 0.563 0.481 0.357 0.431 -0.026 0.537 0.455 0.525 0.533
    2013/5/22 0.632 0.522 0.440 0.494 0.001 0.633 0.523 0.643 0.614
    2014/7/12 0.631 0.563 0.729 0.406 -0.013 0.618 0.550 0.632 0.606
    2015/10/19 0.505 0.477 0.772 0.431 0.086 0.591 0.562 0.605 0.675
    2016/4/28 0.613 0.591 0.427 0.521 -0.011 0.602 0.580 0.702 0.747
    2017/5/1 0.627 0.550 0.546 0.602 0.059 0.685 0.609 0.749 0.705
    2018/6/5 0.671 0.607 0.662 0.667 0.023 0.694 0.630 0.842 0.812
    2019/8/11 0.668 0.618 0.537 0.593 0.015 0.683 0.633 0.830 0.886
    2020/8/13 0.671 0.624 0.618 0.729 0.019 0.689 0.642 0.871 0.918
    下载: 导出CSV
  • [1] CHIGIRA M, WANG W N, FURUYA T, et al. Geological Causes and Geomorphological Precursors of the Tsaoling Landslide Triggered by the 1999 Chi-Chi Earthquake, Taiwan[J]. Engineering Geology, 2003, 68(3-4): 259-273. doi: 10.1016/S0013-7952(02)00232-6
    [2] YANG M D, SU T C, HSU C H, et al. Mapping of the 26 December 2004 Tsunami Disaster by Using FORMOSAT-2 Images[J]. International Journal of Remote Sensing, 2007, 28(13-14): 3071-3091. doi: 10.1080/01431160601094500
    [3] 李小军, 曲国胜, 张晓东. 2005年巴基斯坦北部7.8级地震灾害调查与分析[J]. 震灾防御技术, 2007, 2(4): 354-362. doi: 10.3969/j.issn.1673-5722.2007.04.004
    [4] 许冲, 徐锡伟, 沈玲玲, 等. 2014年鲁甸MS 6.5地震触发滑坡编录及其对一些地震参数的指示[J]. 地震地质, 2014, 36(4): 1186-1203. doi: 10.3969/j.issn.0253-4967.2014.04.020
    [5] 陈晓利, 单新建, 张凌, 等. 地震诱发滑坡的快速评估方法研究: 以2017年MS 7.0级九寨沟地震为例[J]. 地学前缘, 2019, 26(2): 312-320. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201902029.htm
    [6] CUI X, GIBBES C, SOUTHWORTH J, et al. Using Remote Sensing to Quantify Vegetation Change and Ecological Resilience in a Semi-Arid System[J]. Land, 2013, 2(2): 108-130. doi: 10.3390/land2020108
    [7] 吴宁, 卢涛, 罗鹏, 等. 地震对山地生态系统的影响——以5.12汶川大地震为例[J]. 生态学报, 2008, 28(12): 5810-5819. doi: 10.3321/j.issn:1000-0933.2008.12.004
    [8] 王飞龙, 陈明, 王晓迪, 等. 汶川地震震中地区崩滑体的植被恢复与长期效应[J]. 水土保持通报, 2020, 40(3): 175-179. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-STTB202003025.htm
    [9] 李明威, 唐川, 陈明, 等. 汶川震区北川县泥石流流域崩滑体时空演变特征[J]. 水文地质工程地质, 2020, 47(3): 182-190. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG202003023.htm
    [10] 李京忠, 曹明明, 邱海军, 等. 汶川地震区灾后植被恢复时空过程及特征——以都江堰龙溪河流域为例[J]. 应用生态学报, 2016, 27(11): 3479-3486. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201611010.htm
    [11] 田颖颖. 2015年尼泊尔地震震区滑坡的继发性和震后演化研究[D]. 北京: 中国地震局地质研究所, 2020.
    [12] 赵旦, 张淼, 于名召, 等. 汶川地震灾后农田和森林植被恢复遥感监测[J]. 遥感学报, 2014, 18(4): 958-970. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201404016.htm
    [13] doi: https://ejournal.sinica.edu.tw/bbas/content/2002/1/bot431-08.pdf YANG C, CHEN J, PENG L, et al. Chichi Earthquake Caused landslide: Grey Prediction Model for Pioneer Vegetation Recovery Monitored by Satellite Images. Botanical Bulletin of Academia Sinica, 2002, 43: 69-75.
    [14] CHUANG C W, LIN C Y, CHIEN C H, et al. Application of Markov-Chain Model for Vegetation Restoration Assessment at Landslide Areas Caused by a Catastrophic Earthquake in Central Taiwan[J]. Ecological Modelling, 2011, 222(3): 835-845. doi: 10.1016/j.ecolmodel.2010.11.007
    [15] VERDONEN M, BERNER L T, FORBES B C, et al. Periglacial Vegetation Dynamics in Arctic Russia: Decadal Analysis of Tundra Regeneration on Landslides with Time Series Satellite Imagery[J]. Environmental Research Letters, 2020, 15(10): 105020. doi: 10.1088/1748-9326/abb500
    [16] YANG W T, QI W W, ZHOU J X. Decreased Post-Seismic Landslides Linked to Vegetation Recovery after the 2008 Wenchuan Earthquake[J]. Ecological Indicators, 2018, 89: 438-444. doi: 10.1016/j.ecolind.2017.12.006
    [17] 徐锡伟, 闻学泽, 叶建青, 等. 汶川MS 8.0地震地表破裂带及其发震构造[J]. 地震地质, 2008, 30(3): 597-629. doi: 10.3969/j.issn.0253-4967.2008.03.003
    [18] DAI F C, XU C, YAO X, et al. Spatial Distribution of LandslidesTriggered by the 2008 ms 8. 0 Wenchuan Earthquake, China[J]. Journal of Asian Earth Sciences, 2011, 40(4): 883-895. doi: 10.1016/j.jseaes.2010.04.010
    [19] 王根龙, 张军慧, 刘红帅. 汶川地震北川县城地质灾害调查与初步分析[J]. 中国地质灾害与防治学报, 2009, 20(3): 47-51. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH200903012.htm
    [20] 唐川, 齐信, 丁军, 等. 汶川地震高烈度区暴雨滑坡活动的遥感动态分析[J]. 地球科学, 2010, 35(2): 317-323. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201002017.htm
    [21] DAI Z L, HUANG Y, CHENG H L, et al. 3D Numerical Modeling Using Smoothed Particle Hydrodynamics of Flow-Like Landslide Propagation Triggered by the 2008 Wenchuan Earthquake[J]. Engineering Geology, 2014, 180: 21-33. doi: 10.1016/j.enggeo.2014.03.018
    [22] 姜高珍, 韩冰, 高应波, 等. Landsat系列卫星对地观测40年回顾及LDCM前瞻[J]. 遥感学报, 2013, 17(5): 1033-1048. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB201305001.htm
    [23] HORNING N. Review: Advanced Remote Sensing: Terrestrial Information Extraction and Applications[J]. The Quarterly Review of Biology, 2014, 89(2): 199. doi: 10.1086/676116
    [24] 李喆, 胡蝶, 赵登忠, 等. 宽波段遥感植被指数研究进展综述[J]. 长江科学院院报, 2015, 32(1): 125-130. doi: 10.3969/j.issn.1001-5485.2015.01.026
    [25] 彭文甫, 王广杰, 周介铭, 等. 基于多时相Landsat5/8影像的岷江汶川-都江堰段植被覆盖动态监测[J]. 生态学报, 2016, 36(7): 1975-1988. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201607021.htm
    [26] DEERING D W. Rangeland Reflectance Characteristics measured by Aircraft and Spacecraft Sensors[M]. Texas A&M University, 1978.
    [27] 田琰琰. 巴基斯坦克什米尔地震滑坡区植被动态变化特征[D]. 北京: 北京林业大学, 2019.
    [28] 虎振兴, 徐泓, 汪超亮, 等. 基于时间序列遥感影像的滑坡检测方法[J]. 航天返回与遥感, 2018, 39(2): 104-114. doi: 10.3969/j.issn.1009-8518.2018.02.012
    [29] YANG M D, CHEN S C, TSAI H P. A Long-Term Vegetation Recovery Estimation for Mt. Jou-Jou Using Multi-Date SPOT 1, 2, and 4 Images[J]. Remote Sensing, 2017, 9(9): 893. doi: 10.3390/rs9090893
    [30] 李登科, 范建忠, 权文婷. 陕西省植被退化及其驱动因素分析[J]. 生态学杂志, 2015, 34(10): 2907-2913. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-STXZ201510031.htm
    [31] SITCH S, SMITH B, PRENTICE I C, et al. Evaluation of Ecosystem Dynamics, Plant Geography and Terrestrial Carbon Cycling in the LPJ Dynamic Global Vegetation Model[J]. Global Change Biology, 2003, 9(2): 161-185. doi: 10.1046/j.1365-2486.2003.00569.x
    [32] CARLSON T N, RIPLEY D A. On the Relation between NDVI, Fractional Vegetation Cover, and Leaf Area Index[J]. Remote Sensing of Environment, 1997, 62(3): 241-252. doi: 10.1016/S0034-4257(97)00104-1
    [33] 廖清飞, 张鑫, 马全, 等. 青海省东部农业区植被覆盖时空演变遥感监测与分析[J]. 生态学报, 2014, 34(20): 5936-5943. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201420028.htm
    [34] doi: https://www.sciencedirect.com/science/article/pii/S030438000300557X LIN C Y, LO H M, CHOU W C, et al. Vegetation Recovery Assessment at the Jou-Jou Mountain Landslide Area Caused by the 921 Earthquake in Central Taiwan[J]. Ecological Modelling, 2004, 176(1/2): 75-81.
    [35] doi: https://www.sciencedirect.com/science/article/pii/S0169555X01000873 DAI F C, LEE C F. Landslide Characteristics and Slope Instability Modeling Using GIS, Lantau Island, Hong Kong[J]. Geomorphology, 2002, 42(3/4): 213-228.
  • 加载中
图( 5) 表( 2)
计量
  • 文章访问数:  751
  • HTML全文浏览数:  751
  • PDF下载数:  163
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-05-17
  • 刊出日期:  2022-06-20

基于多时相Landsat数据的汶川地震后王家岩滑坡植被恢复估计

    通讯作者: 匡鸿海，副教授
    作者简介: 叶席，硕士研究生，主要从事GIS技术分析与应用的研究
  • 西南大学 地理科学学院，重庆 400715
  • 收稿日期:  2021-05-17
基金项目:  国家自然科学基金项目(41877369)

摘要: 为研究汶川地震后发生严重同震滑坡的北川县王家岩受损植被恢复情况，利用遥感技术对王家岩震后地区进行长期植被动态监测，基于研究区2007-2020年多时相Landsat遥感影像，采用归一化植被指数(NDVI)和植被覆盖恢复率(VRR)来定量评价植被生长和恢复状态. 结果表明：研究区内滑坡区域受损植被已逐渐恢复，其NDVI由地震发生后的0.385上升至2020年8月的0.624；经过12年的演替，滑坡区域植被覆盖恢复率经季节性调整后接近91.8%；根据VRR回归模型预估，滑坡区域内的植被需要26年左右才能恢复到地震前的水平.

English Abstract

    全文HTML

  • 开放科学（资源服务）标志码（OSID）：

  • 近年来，全球因地震引发的地质灾害频有发生[1-5]，这些灾害对植被造成了巨大破坏，对当地自然环境和生态系统造成一定程度的影响. 地震灾区自然环境和生态系统的恢复是长期且复杂的过程，地表植被的生长状况是评价这一过程的极为重要的指标[6]，因此对地震灾区植被生长情况的监测不仅能够为该地区生态系统恢复的评估提供依据，还能预测其未来的发展趋势. 遥感数据因其观测范围广、时效性强等特点在植被的动态监测和定量评估过程中体现出了巨大的优势，是植被覆盖恢复动态监测与评价最有效的数据来源[7].

    学者们利用卫星影像和航空数据对地震灾区植被恢复做了大量研究，王飞龙等[8]利用多期Landsat影像数据分析了汶川地震震中附近震后崩滑体上的植被动态恢复变化，发现该地区经过近9年的恢复，植被覆盖度恢复到0.74，与震前相比差值为0.08，根据拟合模型预计，2022年植被覆盖度能恢复到震前水平；李明威等[9]利用不同来源的遥感影像分析了北川县泥石流流域崩滑体时空演变特征，流域内植被覆盖度在2008年“9.24”泥石流灾害后呈稳定恢复，到2016年研究区植被覆盖度已恢复至较高水平；李京忠等[10]利用中分辨率成像光谱仪—植被指数(MODIS-NDVI)对都江堰龙溪河流域植被恢复进行了定量评估，认为植被覆盖度对地震损害的响应存在滞后现象；田颖颖[11]利用多其次滑坡编目数据、MODIS-NDVI数据、无人机航拍数据等分析了2015年尼泊尔地震后地震区植被演化特征，结果发现震区植被在地震发生后以减少为主，并在震后前两年明显增加后逐渐减少，截至2019年，震区内的植被和松散斜坡物质仍没有达到稳定状态；赵旦等[12]利用机载高分影像对汶川地震影响较为严重区域的农田和森林进行监测，发现震后5年农田的恢复较低，森林恢复情况总体较好；Yang等[13]基于SPOT数据利用灰色预测模型预测台湾1999年“9.21”地震受损植被恢复情况，发现先锋植被恢复时间约为2年，相对实际情况约有2个月的滞后期；Chuang等[14]基于SPOT数据利用马尔科夫链模型对台湾“9.21”地震受损植被恢复状况做出了评估，结果表明草可以作为滑坡区植被恢复的参考指标和滑坡区植被恢复的先导物种进行单独提取；Verdonen等[15]基于QuickBird和WorldView数据利用NDVI来监测俄罗斯北部滑坡冻土的植被状况，结果显示，山体滑坡对苔原植被产生显著而持久的影响，随着气候持续变暖，热喀斯特滑坡及其对植被的影响可能会在西伯利亚西北部和其他北极地区变得越来越普遍；Yang等[16]利用多年滑坡调查数据和MODIS-NDVI数据对汶川地震震中附近受损植被恢复情况进行了研究，结果显示，震后植被恢复可能促进了震后滑坡活动的减少，汶川地震对区域地震后滑坡频率的影响可能在地震后20年内消失.

    以上这些研究在地震后植被恢复监测和分析方面做出了重要贡献. 然而，通过时间序列遥感数据分析受损植被恢复情况并对其恢复到地震前水平做出预估的相关研究目前相对较少. 2008年5月12日，四川省汶川县发生了MS 8.0地震，震中位于四川省汶川县映秀镇(31.021°N，103.367°E)，震源深度约14 km[17]. 汶川地震导致约56 000处同震滑坡发生，总面积约为811 km2[18]，灾害造成相关地区植被损坏严重. 汶川地震所引发的北川县王家岩的同震滑坡，导致北川县大量的房屋被掩埋和人员伤亡[19]，自然环境也受到重创，因此对该地区的植被恢复情况进行研究，可以为地震后自然环境和生态系统恢复、区域规划发展等提供借鉴.

1.   研究区与研究方法

    1.1.   研究区概况

  • 研究区域为四川省北川县曲山镇王家岩及其周边区域(图 1)，面积约9.67 km2. 该区域属于亚热带季风性湿润气候，四季分明，气候温和，区域内雨量充沛，年均降水量约为1 280 mm，降水集中于6~9月；研究区海拔为630~1 610 m，地貌类型以山地、河谷为主，该区域属龙门山前山与后山交界地带，紧邻汶川地震的主发震断裂映秀—北川断裂的上盘，属于较为典型的地质构造不稳定区域[20]，主要由寒武系砂岩、页岩和片岩组成[21].

    • 1.2.   数据来源及数据预处理

  • 由于Landsat影像数据具有时间序列长、易获取、分辨率较高等特点，本研究采用Landsat影像为数据源. 由于Landsat 7的ETM+传感器的扫描线校正仪(SLC)于2003年5月31日发生故障，导致每景ETM+影像存在约22%的数据扫描空隙，产生条带现象，使得该数据的数据质量与数据应用都受到严重的制约和影响[22]，因此本研究选取成像时间为2007/05/06，2008/07/11，2009/06/12，2010/03/27，2011/05/07获取的轨道号为129/38的5期Landsat 5 TM影像，以及成像时间为2013/05/22，2014/07/12，2015/10/19，2016/04/28，2017/05/01，2018/06/05，2019/08/11，2020/08/13获取的轨道号为129/38的8期Landsat 8 OLI影像为本研究的数据源(数据来源于https://earthexplorer.usgs.gov/). 利用ENVI 5.3软件对13个不同时相的遥感影像进行辐射定标、FLAASH大气校正以及目标区域掩膜提取等预处理. 经预处理后的假彩色影像如图 2所示.

    • 1.3.   归一化植被指数(NDVI)计算

  • 归一化植被指数(NDVI)是反映植被生长状况的重要指标因子，是目前应用最为广泛的指标因子之一[23]，常用于对区域植被生长状况的定量评估. NDVI能减少地形地貌的干扰，克服红光波段反射率特别小(接近0)的情况，能够将绿色植被信息与土壤背景等信息进行区分，NDVI值介于-1~1之间，其中0和负值代表无植被覆盖的表面[24]. NDVI数值越大代表植被的覆盖状况越好，植被的生物量越多[25]. NDVI通常被定义为近红外波段与可见光红光波段反射率之差与反射率之和的比值[26]，其表达式为：

    其中：NIR为近红外波段的反射率；R为可见光红光波段的反射率.

    • 1.4.   滑坡区域提取

  • 由于地震前(2007/05/06)后(2008/07/11)2个时相遥感影像所得的NDVI数据存在一定的差异，本研究将研究区内地震前后时相的NDVI差值数据划分为多个小区域，再分别调整每个小区域的阈值，提取出研究区内震后的滑坡分布范围[27]. 利用多时相进行变化检测时容易将与滑坡光谱变化特征相似的区域识别为滑坡[28]，因此同时对研究区震后的Google Earth高分辨率影像进行滑坡的目视解译，再将二者的结果相结合，得到最终的滑坡区域.

    • 1.5.   季节性调整

  • 为保证研究精度，本研究应选取云量在30%以下的遥感影像，在同震滑坡所处时期的研究区内多期遥感影像的积云均较多，为满足研究精度，需要选取各年份不同时段的遥感影像. 本研究选取的遥感影像获取时段处于3-10月之间. 虽然所选影像避开了头年11月至次年2月天气较寒冷因而植被生长条件较差的时段，但由于研究区内的地貌和气候条件使得该区域内植物生长的季候差异明显，因此需要对所选取的遥感影像进行季节性调整处理. 本研究采用Yang等[29]提出的通过简单的季节性调整来适应多时相图像的季节性变化的方法来修正选取的遥感影像之间的植被季候差异. 该方法假设自然季节的影响在遥感影像上的空间分布是均匀的，以仅受季节性影响而未受同震滑坡影响的非滑坡区作为参考区域，从每幅遥感影像中提取非滑坡区域所有像元的NDVI并取其平均值作为背景值，计算每幅遥感影像背景值之间的差异，并以此作为季节性调整的偏移值[30]. 利用所得的季节性调整偏移值对滑坡区域和研究区域的NDVI进行调整后，得到调整后的NDVI，进而对每幅遥感影像进行调整后的植被覆盖恢复率(VRR)计算.

    • 1.6.   VRR计算

  • 遥感数据是区域范围内植被覆盖度提取的最有效数据源[31]. 植被覆盖度(FVC)是反映植被生态功能的动态变量，可用于监测环境的变化和评价受损植被的恢复状况[29]. NDVI被广泛应用于植被覆盖度的计算，利用NDVI与光谱混合分析(NDVI-SMA)相结合的方法提取植被覆盖度，其原理是假设给定像元由植被与非植被覆盖的地表两部分所构成，其光谱信息只由这两个组分线性合成，它们各自的面积在像元中所占的比率即为各因子的权重，由它们的相对比例加权得到每个像元，其中植被所占像元的百分比即为该像元的FVC[32]. 植被覆盖度表达式为：

    式中：NDVIveg为全植被覆盖像元的NDVI值；NDVIsoil为裸土像元的NDVI值；FVC的取值范围为0~1. 因受噪声、植被分布情况以及邻近地物辐射等因素的影响，植被覆盖度置信度的取值主要由研究区域遥感影像的实际情况来决定，在没有大量实测数据作参考的情况下，以累积频率为5%和95%作为置信区间[33]. 本研究选择累积频率为5%的NDVI作为裸土像元的NDVI值，累积频率为95%的NDVI作为全植被覆盖像元的NDVI值来进行植被覆盖度的计算.

    植被覆盖恢复率(VRR)可以定义为植被覆盖度受损后恢复的速率，它是根据植被受损前后的差异来监测和评估植被的恢复情况，以植被覆盖度为基础来进行计算[34]. 植被覆盖恢复率的表达式为：

    式中：FVC0为同震滑坡前的植被覆盖度；FVC1为同震滑坡后的植被覆盖度；FVCi为同震滑坡后i(i=1，2，…，12)年植被覆盖度.

2.   结果

    2.1.   滑坡提取结果

  • 首先根据地震前后时相NDVI差值的阈值提取出滑坡分布的范围，同时结合震后Google Earth的高分辨率影像对研究区内滑坡进行目视解译的成果，得到研究区内滑坡区域(图 3)，其中滑坡区域面积约为0.94 km2，非滑坡区域面积约为8.73 km2.

    • 2.2.   NDVI计算结果

  • 本研究中的NDVI均通过逐像素计算得到，研究区NDVI提取结果如图 4所示. 根据滑坡提取结果，本研究将研究区域划分为非滑坡区域以及滑坡区域，并分别计算出研究区域、非滑坡区域以及滑坡区域的NDVI均值，其结果如表 1所示. 滑坡区域的NDVI在同震滑坡发生后立即急剧下降，并在随后的年份出现波动，但是NDVI总体呈现上升趋势，对比图 4，可以清晰看到滑坡区域植被恢复显著. 在地震发生前，滑坡区域的植被覆盖良好，平均NDVI为0.732；由于研究区内建筑物、河流、坑塘、水坝等的存在，使得非滑坡区域的平均NDVI低于滑坡区域，非滑坡区域的平均NDVI为0.646；整个研究区平均NDVI为0.654，表明震前该区域植被条件较好，但地震以及滑坡发生后，滑坡区域的平均NDVI急剧下降至0.385. 12年后的2020年8月13日，研究区、非滑坡区域、滑坡区域的平均NDVI分别上升到0.671，0.675，0.624. 由于选取的遥感影像时段处于3月至10月之间，存在一定的物候差异，从2009年6月12日到2010年3月27日、2014年7月12日到2015年10月19日，研究区内的NDVI水平均呈现较大辐度的下降. 因此通过一定的季节性调整方法对数据进行修正以降低植被的物候差异是十分必要的.

3.   讨论与分析

    3.1.   VRR分析

  • 本研究利用NDVI和经过季节性调整后的NDVI计算出FVC，进而得到VRR和调整后的VRR，结果如表 2所示. 为了评估整个研究区和滑坡区域的VRR变化情况，对研究区(图 5a)和滑坡区域(图 5b)构建VRR回归模型，并对其进行趋势分析.

    地震以及滑坡发生之后，研究区及滑坡区的VRR总体呈现上升趋势，在震后早期，受损植被的恢复速率相对较快，然后逐渐趋于平缓. 从整个研究区来看，2010年3月27日、2014年7月12日以及2015年10月19日其VRR存在明显异常，其中2010年3月27日的VRR相对偏低，其值为0.051，而2014年7月12日和2015年10月19日的VRR相对偏高，其值分别为0.729和0.772. 出现异常情况可能是由季节因素导致的植物物候差异、降水条件的改变，以及人类活动的影响等所致. 从滑坡区域来看，震后VRR的波动范围较小，总体呈现逐步上升的趋势. 决定系数R2和p值分别用来衡量季节调整VRR回归模型的拟合程度以及显著性水平，从研究区和滑坡区域的VRR回归模型的p值来看，VRR回归模型是显著的. 从决定系数R2来看，滑坡区域VRR回归模型R2值为0.863，研究区VRR回归模型R2值为0.688，滑坡区域的VRR回归模型优于整个研究区的VRR回归模型.

    • 3.2.   调整后的VRR分析

  • 震后研究区及滑坡区域经季节性调整后的VRR值总体呈逐渐上升趋势，研究区及滑坡区域调整后的VRR的波动范围均较小，受损植被处于逐步恢复的状态. 由研究区和滑坡区域调整后的VRR回归模型的p值可知，调整后的VRR回归模型是显著的. 对于整个研究区(图 5a)，调整后的VRR回归模型R2值为0.961，较原VRR模型的R2值(0.688)有显著提高，提高了0.273. 对于滑坡区域(图 5b)，调整后的VRR回归模型R2值为0.957，较原VRR回归模型R2值增加了0.094，存在明显提升. 滑坡区域调整后VRR回归模型R2值为0.961，而研究区调整后VRR回归模型R2值为0.863，从决定系数R2来看，研究区调整后的VRR回归模型优于滑坡区. 此外，就整个研究区而言，2010年3月27日观测到相对较低的VRR值为0.051，而2014年7月12日和2015年10月19日观测到相对较高的VRR分别为0.729和0.772. 经季节性调整后，2010年3月27日相对较低的VRR向上调整为0.416，2014年7月12日和2015年10月19日相对较高的2个VRR分别向下调整为0.632和0.604.

    在对研究区NDVI水平进行季节性调整后，对比调整前后所得的VRR数据，可以发现无论是对于整个研究区还是滑坡区域，经季节性调整后，VRR的数值波动范围均有所减小，修正了VRR的大幅异常波动情况. 同时VRR模型在R2和p值上均优于未作季节性调整的原始VRR模型，说明对NDVI进行季节性修正有助于后续的VRR计算以及未来植被恢复情况的预测，季节性调整能较好地估计植被恢复状况. 根据滑坡区域调整后的VRR模型，本研究预测该滑坡区域的受损植被大约需要26年才能完全恢复.

    • 3.3.   本研究方法对其他区域的适用性

  • 本研究采用多时相Landsat遥感影像对研究区域植被的时空变化情况进行了监测和分析：首先对研究区2007-2020年的NDVI进行了分析，再以植被覆盖度为基础建立VRR模型，以及经季节性调整后的VRR模型，并进行相应的分析. 研究区内受损植被震后早期的恢复速率相对较快，之后逐渐趋于平缓. 根据滑坡区域经季节性调整后的VRR模型数据，本研究预测该滑坡区域的受损植被大约需要26年才能完全恢复. 本研究选取的研究区域为典型的山地环境同震滑坡区域，并且研究区内植被随季节变化的物候差异明显. 本研究所采用的研究方法可以为山地环境同震滑坡引起的植被在时空范围上的变化情况以及受损植被完全恢复的合理预估提供参考，可应用于其他类似区域的相应研究.

4.   结论

  • 本研究利用多时相Landsat遥感影像对汶川地震后北川县王家岩地区受到同震滑坡影响的植被的时空变化进行了动态监测与分析，并基于2007-2020年的13景Landsat遥感影像，对研究区的NDVI进行了分析，建立了相应的VRR回归模型，同时为消除季节变化对植被的影响，对研究数据进行了季节性调整. 结果显示，在地震后12年的植被演替过程中，滑坡区域的NDVI总体水平呈上升趋势，但与地震发生前的植被条件相比还有一定差距；对于整个研究区和滑坡区域，季节性调整后的VRR的相关系数高于原始VRR系数，表明在植物生长季候明显的区域，利用多时相遥感影像进行长期的植被观测时需要保持相同的季节或对不同季节进行季节性调整. 本研究根据滑坡区域调整后的VRR模型，估计该滑坡区域的受损植被大约需要26年才能完全恢复. 本研究对于同震滑坡后山区植被恢复的长期跟踪监测与未来趋势预测具有一定参考意义，可应用于其他类似区域的研究.

    地震后同震滑坡体上的植被恢复是一个长期且复杂的过程，目前对于引起这类变化的控制因素的研究成果相对较少，本研究也仅基于长时间序列NDVI数据来统计分析震后滑坡体植被恢复特征及变化趋势，进而评估震后滑坡体植被恢复完成的时间. 例如，坡度是导致滑坡以及影响之后植被恢复的重要因素之一[35]，但本研究在分析同震滑坡后植被恢复的情况时，未结合坡度、海拔、坡向、降水、水系分布以及人类活动等影响植被生长和恢复的因子进行综合评估，在今后的研究中，可以利用最新的卫星图像和不同系列的传感器，在考虑季节变化和地形因子、降水等影响的基础上，建立可靠的植被恢复回归模型，对植被恢复进行持续监测.

参考文献 (35)

目录

版权所有：西南大学学报（自然科学版）

地址：重庆市北碚区西南大学 电话：023-68254576, 68252540, 68366165,68366917 电子邮箱：lkxb@swu.edu.cn

本系统由北京仁和汇智信息技术有限公司开发

/

返回文章
返回