留言板

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

姓名
邮箱
手机号码
标题
留言内容
验证码

基于断点检测的延安市富县植被变化研究

季烨云 刘鹏举 高影 任怡 余涛 牛思圆

季烨云, 刘鹏举, 高影, 任怡, 余涛, 牛思圆. 基于断点检测的延安市富县植被变化研究[J]. 浙江农林大学学报, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
引用本文: 季烨云, 刘鹏举, 高影, 任怡, 余涛, 牛思圆. 基于断点检测的延安市富县植被变化研究[J]. 浙江农林大学学报, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
JI Yeyun, LIU Pengju, GAO Ying, REN Yi, YU Tao, NIU Siyuan. Vegetation change based on break-point detection in Fu County, Yan’an City[J]. Journal of Zhejiang A&F University, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
Citation: JI Yeyun, LIU Pengju, GAO Ying, REN Yi, YU Tao, NIU Siyuan. Vegetation change based on break-point detection in Fu County, Yan’an City[J]. Journal of Zhejiang A&F University, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457

基于断点检测的延安市富县植被变化研究

doi: 10.11833/j.issn.2095-0756.20210457
基金项目: 国家重点研发计划项目(2017YFD0600906);“三北”防护林体系管理与服务平台(林防合字2019-17)
详细信息
    作者简介: 季烨云(ORCID: 0000-0001-7711-2084),从事地理信息系统技术与应用研究。E-mail: 15850781932@163.com
    通信作者: 刘鹏举(ORCID: 0000-0002-5502-6334),副研究员,从事林业地理信息系统应用与开发研究。E-mail: liupengju@ifrit.ac.cn
  • 中图分类号: TP79

Vegetation change based on break-point detection in Fu County, Yan’an City

  • 摘要:   目的  多项生态工程实施以来,陕北黄土高原地区环境明显改善,探究该区域植被的时空变化,为生态环境建设与差异化决策提供依据。  方法  以陕西省延安市富县为研究对象,利用2000—2019年生长季的Landsat TM/ETM+/OLI影像,运用DBEST算法对归一化植被指数(NDVI)时间序列进行断点检测,并提出多阶段趋势分析方法分析富县植被变化时空特征。  结果  富县植被的断点检测结果表明:66.47%的区域断点NDVI变化幅度小于0.2,主要分布在富县西部和东北地区,植被相对稳定,未发生剧烈变化;33.53%的区域断点NDVI变化幅度大于0.2,其中断点数量多于4个的区域仅占5.88%,集中分布于道路、河流沿线,变化频繁,与人为活动有关。当前阶段趋势分析表明:断点NDVI变化幅度大于0.2的区域,改善的面积占总面积的24.57%,退化面积仅占2.12%;开始时间均在2014年之前,时间分布跨度大,空间异质性强,揭示了富县植被变化的多样性及复杂性;富县有林地、疏林地、灌木林地、未成林造林地4种主要林地均趋于改善,改善面积占林地面积的96.23%。  结论  2000—2019年,富县植被整体呈现改善趋势,生态建设工程取得良好的效果;植被变化的时空特征存在差异,后续决策需因地制宜。图6表1参25
  • 图  1  2019年富县基础地理信息示意图

    Figure  1  Basic geographic information of Fu County in 2019

    图  2  2019年富县林地类型与分布示意图

    Figure  2  Forest land types and distribution of Fu County in 2019

    图  3  富县归一化植被指数变化幅度及其面积占比      

    Figure  3  Change magnitude and area percent of NDVI in Fu County      

    图  4  富县不同NDVI变化幅度小于0.2与大于0.2的断点数目及其面积占比

    Figure  4  Number of breakpoints and area percent of NDVI (<0.2 and >0.2) in Fu County

    图  5  2000—2019年富县NDVI时间序列总体趋势类别及其面积占比

    Figure  5  Overall trend level and area percent of NDVI time series in Fu County from 2000 to 2019

    图  6  2000—2019年富县NDVI时间序列当前阶段趋势及其面积占比

    Figure  6  Current stage trend and area percent of NDVI time series in Fu County from 2000 to 2019

    表  1  富县林地类型趋势情况统计

    Table  1.   Statistics on trend of forest land types in Fu County

    林地类型改善/%基本不变/%退化/%
    有林地96.980.862.16
    疏林地94.321.913.77
    灌木林地94.411.813.78
    未成林造林地92.971.915.12
    其他用地86.662.5010.84
    下载: 导出CSV
  • [1] HEIMANN M, REICHSTEIN M. Terrestrial ecosystem carbon dynamics and climate feedbacks [J]. Nature, 2008, 451(7176): 289 − 292.
    [2] 徐涵秋. 城市遥感生态指数的创建及其应用[J]. 生态学报, 2013, 33(24): 7853 − 7862.

    XU Hanqiu. A remote sensing urban ecological index and its application [J]. Acta Ecol Sin, 2013, 33(24): 7853 − 7862.
    [3] 许文鑫, 周玉科, 梁娟珠. 基于变化点的青藏高原植被时空动态变化研究[J]. 遥感技术与应用, 2019, 34(3): 667 − 676.

    XU Wenxin, ZHOU Yuke, LIANG Juanzhu. A breakpoints based spatio-temporal analysis of Tibetan Plateau vegetation [J]. Remote Sensing Technol Appl, 2019, 34(3): 667 − 676.
    [4] CAO Qian, WU Jianguo, YU Deyang, et al. The biophysical effects of the vegetation restoration program on regional climate metrics in the Loess Plateau, China [J]. Agric For Meteorol, 2019, 268: 169 − 180.
    [5] HANSEN M C, POTAPOV P V, MOORE R, et al. High-resolution global maps of 21st-century forest cover change [J]. Science, 2013, 342(6160): 850 − 853.
    [6] 孙雷刚, 刘剑锋, 徐全洪. 河北坝上地区植被覆盖变化遥感时空分析[J]. 国土资源遥感, 2014, 26(1): 167 − 172.

    SUN Leigang, LIU Jianfeng, XU Quanhong. Remote Sensing based temporal and spatial analysis of vegetation cover changes in Bashang area of Hebei Province [J]. Remote Sensing Land Resour, 2014, 26(1): 167 − 172.
    [7] 刘咏梅, 马黎, 黄昌, 等. 基于MODIS-Landsat时空融合的陕北黄土高原植被覆盖变化研究[J]. 西北大学学报(自然科学版), 2019, 49(1): 62 − 70.

    LIU Yongmei, MA Li, HUANG Chang, et al. Study on the change of vegetation coverage of Loess Plateau in northern Shaanxi Province based on MODIS-Landsat fusion data [J]. J Northwest Univ Nat Sci Ed, 2019, 49(1): 62 − 70.
    [8] QAMER F M, SHEHZAD K, ABBAS S, et al. Mapping deforestation and forest degradation patterns in western Himalaya, Pakistan [J]. Remote Sensing, 2016, 8(385): 1 − 17.
    [9] 沈文娟, 李明诗, 黄成全. 长时间序列多源遥感数据的森林干扰监测算法研究进展[J]. 遥感学报, 2018, 22(6): 1005 − 1022.

    SHEN Wenjuan, LI Mingshi, HUANG Chengquan. Review of remote sensing algorithms for monitoring forest disturbance from time series and multi-source data fusion [J]. J Remote Sensing, 2018, 22(6): 1005 − 1022.
    [10] 刘宝柱, 方秀琴, 何祺胜, 等. 基于MODIS数据和BFAST方法的植被变化监测[J]. 国土资源遥感, 2016, 28(3): 146 − 153.

    LIU Baozhu, FANG Xiuqin, HE Qisheng, et al. Monitoring the changes of vegetation based on MODIS data and BFAST methods [J]. Remote Sensing Land Resour, 2016, 28(3): 146 − 153.
    [11] FORKEL M, CARVALHAIS N, VERBESSELT J, et al. Trend change detection in NDVI time series: effects of inter-annual variability and methodology [J]. Remote Sensing, 2013, 5(5): 2113 − 2144.
    [12] STOW D, PETERSEN A, HOPE A, et al. Greenness trends of Arctic tundra vegetation in the 1990s: comparison of two NDVI data sets from NOAA AVHRR systems [J]. Int J Remote Sensing, 2007, 28(21): 4807 − 4822.
    [13] 李斌. 青藏高原植被时空分布规律及其影响因素研究[D]. 北京: 中国地质大学, 2016.

    LI Bin. Study of Vegetation’s Spatial and Temporal Distribution and Influencing Factors on the Tibetan Plateau[D]. Beijing: China University of Geosciences, 2016.
    [14] 海月, 杨广斌, 李若男, 等. 基于时序突变检测的植被空间变化特征识别方法研究——以海河北部山区为例[J]. 生态学报, 2020, 40(24): 1 − 10.

    HAI Yue, YANG Guangbin, LI Ruonan, et al. Recognition of vegetation spatial variation based on time-series mutation detection: a case study of the mountainous area of Northern Haihe River Basin [J]. Acta Ecol Sin, 2020, 40(24): 1 − 10.
    [15] KENNEDY R E, YANG Zhiqiang, COHEN W B. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr-Temporal segmentation algorithms [J]. Remote Sensing Environ, 2010, 114(12): 2897 − 2910.
    [16] VERBESSELT J, HYNDMAN R, NEWNHAM G, et al. Detecting trend and seasonal changes in satellite image time series [J]. Remote Sensing Environ, 2010, 114(1): 106 − 115.
    [17] JAMALI S, JONSSON P, EKLUNDH L, et al. Detecting changes in vegetation trends using time series segmentation [J]. Remote Sensing Environ, 2015, 156: 182 − 195.
    [18] 刘吉彤, 安如. 基于DBEST的2003—2012年济南市植被变化研究[J]. 甘肃科学学报, 2018, 30(4): 17 − 23.

    LIU Jitong, AN Ru. Study on vegetation change in Ji’nan during 2003−2012 based on DBEST [J]. J Gansu Sci, 2018, 30(4): 17 − 23.
    [19] 王恩鲁, 汪小钦, 陈芸芝. 时间序列植被覆盖度断点检测方法研究[J]. 地球信息科学学报, 2017, 19(10): 1355 − 1363.

    WANG Enlu, WANG Xiaoqin, CHEN Yunzhi. The breakpoints detection method using time series of vegetation fractional coverage [J]. J Geo-Inf Sci, 2017, 19(10): 1355 − 1363.
    [20] SHEN Xiaoji, AN Ru, FENG Li, et al. Vegetation changes in the Three-River Headwaters Region of the Tibetan Plateau of China [J]. Ecol Indic, 2018, 93: 804 − 812.
    [21] 王颖洁, 刘良云, 王志慧. 基于时序Landsat数据的三江平原植被地表类型变化遥感探测研究[J]. 遥感技术与应用, 2015, 30(5): 959 − 968.

    WANG Yingjie, LIU Liangyun, WANG Zhihui. Land covermapping based on Landsat time-series stacks in Sanjiang Plain [J]. Remote Sensing Technol Appl, 2015, 30(5): 959 − 968.
    [22] 贾铎, 王藏姣, 牟守国, 等. 基于NDVI时间序列轨迹的草原露天矿区植被时空动态特征[J]. 应用生态学报, 2017, 28(6): 1808 − 1816.

    JIA Duo, WANG Cangjiao, MU Shouguo, et al. Vegetation spatial and temporal dynamic characteristics based on NDVI time series trajectories in grassland opencast coal mining [J]. Chin J Appl Ecol, 2017, 28(6): 1808 − 1816.
    [23] 陈昀琳. 基于Landsat和MODIS NDVI时序数据的青海湖流域植被覆盖度提取及其变化分析[D]. 北京: 中国地质大学, 2019.

    CHEN Yunlin. Fractional Vegetation Cover in Qinghai Lake Watershed Extraction and Change Analysis based on Time-series NDVI Data of Landsat and MODIS[D]. Beijing: China University of Geosciences, 2019.
    [24] 李杰, 张军, 刘陈立, 等. 基于MODIS-NDVI的中老缅交界区近16年植被覆盖时空变化特征[J]. 林业科学, 2019, 55(8): 9 − 18.

    LI Jie, ZHANG Jun, LIU Chenli, et al. Spatiotemporal variation of vegetation coverage in recent 16 years in the border region of China, Laos, and Myanmar based on MODIS-NDVI [J]. Sci Silv Sin, 2019, 55(8): 9 − 18.
    [25] 宋怡, 马明国. 基于SPOT VEGETATION数据的中国西北植被覆盖变化分析[J]. 中国沙漠, 2007, 27(1): 89 − 94.

    SONG Yi, MA Mingguo. Study on vegetation cover change in Northwest China based on SPOT VEGETATION data [J]. J Desert Res, 2007, 27(1): 89 − 94.
  • [1] 胡瑞祺, 方卫军, 王义平.  天目山不同海拔大蛾类种群多样性及时间动态 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20230469
    [2] 李思源, 叶真妮, 毛勇伟, 陈玉玲, 曾纳.  4种遥感数据时空融合模型生成高分辨率归一化植被指数的对比分析 . 浙江农林大学学报, 2023, 40(2): 427-435. doi: 10.11833/j.issn.2095-0756.20220381
    [3] 徐悦, 吴筱萌, 陈诗燕, 程鸿浩, 张淋, 周夏芝, 邹运鼎, 毕守东.  茶园蜘蛛与广翅蜡蝉的时间和空间关系 . 浙江农林大学学报, 2022, 39(5): 1067-1079. doi: 10.11833/j.issn.2095-0756.20210734
    [4] 杨绍钦, 王翔, 许澄, 商天其.  基于MODIS时间序列数据的竹林地上生物量估算 . 浙江农林大学学报, 2022, 39(4): 734-741. doi: 10.11833/j.issn.2095-0756.20210431
    [5] 庞天虹, 钱婕妤, 付建新, 顾翠花, 张超.  桂花己糖激酶基因家族成员的序列及表达分析 . 浙江农林大学学报, 2021, 38(2): 225-234. doi: 10.11833/j.issn.2095-0756.20200370
    [6] 魏玮, 李俊敏, 孙丽英, 戎均康, 周伟.  小麦黄花叶病抗性鉴定及抗性亲本简单重复序列 . 浙江农林大学学报, 2016, 33(1): 71-79. doi: 10.11833/j.issn.2095-0756.2016.01.010
    [7] 王琳, 景元书, 张悦.  基于MODIS的长江中下游地区植被净第一性生产力时空变化规律 . 浙江农林大学学报, 2015, 32(6): 829-836. doi: 10.11833/j.issn.2095-0756.2015.06.002
    [8] 李亚丹, 杜华强, 周国模, 谷成燕, 徐小军, 孙少波, 高国龙.  雷竹叶绿素与高光谱植被指数关系及其反演模型 . 浙江农林大学学报, 2015, 32(3): 335-345. doi: 10.11833/j.issn.2095-0756.2015.03.002
    [9] 刘鑫, 张金池, 谢德晋, 庄家尧, 邵永昌, 张水锋.  南京近郊麻栎林冠层蒸腾时间变化及降温效应 . 浙江农林大学学报, 2015, 32(4): 529-536. doi: 10.11833/j.issn.2095-0756.2015.04.006
    [10] 伍周玲, 梁东军, 郭明, 王晓萌, 范文翔.  新型固定化乙酰胆碱酯酶传感器的制备及氨基甲酸酯农残检测 . 浙江农林大学学报, 2014, 31(3): 457-464. doi: 10.11833/j.issn.2095-0756.2014.03.020
    [11] 尹建新, 祁亨年, 冯海林, 杜晓晨.  一种基于混合纹理特征的木板材表面缺陷检测方法 . 浙江农林大学学报, 2011, 28(6): 937-942. doi: 10.11833/j.issn.2095-0756.2011.06.017
    [12] 桂仁意, 金爱武, 高培军, 曹福亮.  银杏核糖体DNA内转录间隔序列初步分析 . 浙江农林大学学报, 2007, 24(1): 17-21.
    [13] 邝粉良, 刘宁, 仓决卓玛, 李建川, 杨乐, 李福秀.  藏北黑颈鹤繁殖前期的昼间活动时间分配 . 浙江农林大学学报, 2007, 24(6): 686-691.
    [14] 鲍滨福, 余文军, 姜志宏, 俞友明.  热压温度和养护时间对快速固化水泥刨花板性能的影响 . 浙江农林大学学报, 2006, 23(2): 129-134.
    [15] 鲁庆彬, 王小明, 丁由中, 李新庆.  扬子鳄发情期昼间行为时间分配及活动节律 . 浙江农林大学学报, 2006, 23(6): 608-613.
    [16] 周莉花, 郝日明, 赵宏波.  蜡梅花粉活力检测方法筛选及保存时间观察 . 浙江农林大学学报, 2006, 23(3): 270-274.
    [17] 钱俊, 金永明, 沈哲红, 俞友明, 楼永生.  电场强度与作用时间对马尾松木材含水率梯度的影响 . 浙江农林大学学报, 2005, 22(2): 193-197.
    [18] 金永明, 鲍滨福, 陈庆虎, 余学军, 俞友明, 叶良明, .  陈放时间和灰木比对快速固化水泥刨花板性能的影响 . 浙江农林大学学报, 2003, 20(3): 232-235.
    [19] 吴章文.  亚热带森林旅游区夏季舒适温度的持续时间 . 浙江农林大学学报, 2003, 20(4): 380-384.
    [20] 陈波.  陆地植被净第一性生产力对全球气候变化响应研究的进展 . 浙江农林大学学报, 2001, 18(4): 445-449.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20210457

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2022/3/540

图(6) / 表(1)
计量
  • 文章访问数:  490
  • HTML全文浏览量:  140
  • PDF下载量:  36
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-06-29
  • 修回日期:  2021-12-23
  • 录用日期:  2021-12-31
  • 网络出版日期:  2022-05-23
  • 刊出日期:  2022-05-23

基于断点检测的延安市富县植被变化研究

doi: 10.11833/j.issn.2095-0756.20210457
    基金项目:  国家重点研发计划项目(2017YFD0600906);“三北”防护林体系管理与服务平台(林防合字2019-17)
    作者简介:

    季烨云(ORCID: 0000-0001-7711-2084),从事地理信息系统技术与应用研究。E-mail: 15850781932@163.com

    通信作者: 刘鹏举(ORCID: 0000-0002-5502-6334),副研究员,从事林业地理信息系统应用与开发研究。E-mail: liupengju@ifrit.ac.cn
  • 中图分类号: TP79

摘要:   目的  多项生态工程实施以来,陕北黄土高原地区环境明显改善,探究该区域植被的时空变化,为生态环境建设与差异化决策提供依据。  方法  以陕西省延安市富县为研究对象,利用2000—2019年生长季的Landsat TM/ETM+/OLI影像,运用DBEST算法对归一化植被指数(NDVI)时间序列进行断点检测,并提出多阶段趋势分析方法分析富县植被变化时空特征。  结果  富县植被的断点检测结果表明:66.47%的区域断点NDVI变化幅度小于0.2,主要分布在富县西部和东北地区,植被相对稳定,未发生剧烈变化;33.53%的区域断点NDVI变化幅度大于0.2,其中断点数量多于4个的区域仅占5.88%,集中分布于道路、河流沿线,变化频繁,与人为活动有关。当前阶段趋势分析表明:断点NDVI变化幅度大于0.2的区域,改善的面积占总面积的24.57%,退化面积仅占2.12%;开始时间均在2014年之前,时间分布跨度大,空间异质性强,揭示了富县植被变化的多样性及复杂性;富县有林地、疏林地、灌木林地、未成林造林地4种主要林地均趋于改善,改善面积占林地面积的96.23%。  结论  2000—2019年,富县植被整体呈现改善趋势,生态建设工程取得良好的效果;植被变化的时空特征存在差异,后续决策需因地制宜。图6表1参25

English Abstract

季烨云, 刘鹏举, 高影, 任怡, 余涛, 牛思圆. 基于断点检测的延安市富县植被变化研究[J]. 浙江农林大学学报, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
引用本文: 季烨云, 刘鹏举, 高影, 任怡, 余涛, 牛思圆. 基于断点检测的延安市富县植被变化研究[J]. 浙江农林大学学报, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
JI Yeyun, LIU Pengju, GAO Ying, REN Yi, YU Tao, NIU Siyuan. Vegetation change based on break-point detection in Fu County, Yan’an City[J]. Journal of Zhejiang A&F University, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
Citation: JI Yeyun, LIU Pengju, GAO Ying, REN Yi, YU Tao, NIU Siyuan. Vegetation change based on break-point detection in Fu County, Yan’an City[J]. Journal of Zhejiang A&F University, 2022, 39(3): 540-546. doi: 10.11833/j.issn.2095-0756.20210457
  • 植被是生态系统的重要组成部分,对全球气候变化和生态系统的结构功能都有显著影响[1-3]。研究地表植被覆盖及其长期变化不仅可以揭示地表空间变化规律,而且有利于分析区域生态系统及环境变化情况[4],对因地制宜开展植被的保护管理工作具有重要意义。黄土高原位于黄河中游地区,是世界上水土流失最为严重的区域之一。陕北是黄土高原的中心部分,水土流失严重、生态环境脆弱,但自三北防护林、天然林资源保护、退耕还林还草等工程实施以来,该区域生态环境得到了明显改善。因此有必要开展工程实施效果阶段性评估,便于为今后该区域的林业生态建设和保护提供理论依据。

    时间序列的遥感数据可以动态监测大面积的植被生长状况[5-6]、发展变化趋势特征[78]、扰动情况[9-10]等。现有的大多数研究中,对归一化植被指数(normalized difference vegetation index,NDVI)时间序列的分析多针对空间和时间总体变化展开[11],如STOW等[12]提出用一元线性回归方法来分析每个栅格像元多年的变化趋势,以此来计算植被的绿度变化率(greenness rate of change,GRC)。当植被遭遇砍伐、病虫害、火灾等自然灾害或人类活动影响时,植被绿度变化趋势可能会出现突然的转变。然而,一元线性回归方法只能表现植被变化的整体趋势,不能监测到以上趋势发生突然变化的情况,具有很大的局限性[13-14]。因此,学者们提出了专用于处理遥感时序数据的断点检测方法,如LandTrendr(Landsat based detection of Trends in Disturbance and Recovery)方法[15]、BFAST(breaks for additive seasonal and trend)方法[16]等。由于上述方法受限于卫星传感器、数据类型和数据序列长度等方面的不足,JAMALI等[17]提出了DBEST(detecting breakpoints and estimating segments in trend)算法,该算法不仅对卫星传感器类型和数据序列长度没有限制,而且支持周期性和非周期性植被指数时间序列,检测效果较好且具有普遍的适用性[18-20]

    本研究以延安市富县为例,利用2000—2019年Landsat影像构建NDVI时间序列,运用DBEST算法识别NDVI时间序列的断点,分析断点的时空分布,并研究植被变化多阶段趋势特征,为生态修复工程规划和长期生态环境保护决策提供依据。

    • 富县位于陕西省北部,延安市南部,35°44′06″~36°23′23″N,108°29′30″~109°42′54″E,总面积为4 182 km2,属于渭北黄土高原丘陵沟壑地带。该县属大陆性暖温带季风气候,四季分明,年平均气温为7~9 ℃,年平均降水量为500~600 mm。区域内分布着丰富的森林资源,乔木树种主要有栎类Quercus spp.、桦木Betula spp.、刺槐Robinia pseudoacacia等;灌木主要有胡颓子Elaeagnus pungens、马桑Coriaria nepalensis、白刺花Sophora davidii等。居民地、河流、铁路和公路分布概况如图1所示。

      图  1  2019年富县基础地理信息示意图

      Figure 1.  Basic geographic information of Fu County in 2019

    • Landsat TM/ETM+/OLI多光谱遥感影像来源于美国地质调查局提供的地表反射率(surface reflectance,SR)产品,即L2级别产品(https://earthexplorer.usgs.gov/),空间分辨率为30 m,时间分辨率为16 d,光谱范围包括可见光、近红外、短波红外和热红外波段[21]。由于2003年之后的Landsat 7 ETM+影像出现数据条带丢失的问题,因此需先利用ENVI软件中的Landsat_gapfill插件去除此类影像中的条带[22]。然后利用Fmask云检测算法对所有影像进行云、云影识别提取[23]。最后利用红色和近红外波段计算NDVI,并基于最大值合成法[24](maximum value composite,MVC)合成了2000—2019年的NDVI时间序列数据。

      本研究还收集了研究区2019年的林地类型,数据来源于国家森林资源智慧管理平台中的“森林资源一张图”数据,用于分析不同林地的植被变化情况。各类型面积如图2所示:有林地2 664.77 km2,疏林地82.80 km2,灌木林地554.95 km2,未成林造林地253.85 km2,其他用地625.63 km2,包括人造地表、水体、耕地等。

      图  2  2019年富县林地类型与分布示意图

      Figure 2.  Forest land types and distribution of Fu County in 2019

    • DBEST算法主要是由时间序列分解和趋势分割两部分组成[17]。时间序列分解主要是指植被指数时间序列数据的预处理过程,而趋势分割是实现断点检测的过程。由于本研究对生长季NDVI数据采用最大值生成年度NDVI数据,因此时间序列分解部分未进行以局部移动平均为基础的季节趋势分解(seasonal-trend decomposition procedure based on Loess,STL),仅利用水平位移第1阈值$ {\theta }_{1} $、水平位移第2阈值$ {\theta }_{2} $、持续时长$ \varphi $共3个参数检测候选断点。趋势分割(即断点检测)包括3个过程:①找到时间序列的转点(即峰/谷点、距离阈值点和水平位移点);②通过下一个转点减去当前转点处的NDVI值来计算局部趋势变化(trend local change,TLC);③所有转点按TLC降序重新排列,利用贝叶斯信息准则(bayesian information criterion,BIC)确定断点数量,若BIC准则得到的断点数量少于转点数量,则最终的断点数为BIC准则得到的结果。

      运用DBEST算法检测断点,设置参数如下:水平位移第1阈值($ {\theta }_{1}) $为0.1,即当时序NDVI数据变化超过0.1时认为该点是候选水平位移点;水平位移第2阈值$ ({\theta }_{2} $)为0.2,即当候选水平位移点前后最小步长范围内的子序列均值变化超过0.2认为候选水平位移点符合第2阈值;持续时长($ \varphi $)为2 a,在该时段内只能有1个水平位移点,以检测NDVI短时间内的平均变化;以上3个参数可确定时间序列中的水平位移点。变化级别($ \beta $)、泛化百分比($ \delta $)、断点数目(m)共3个参数可根据需求三选一达到不同的目的。由于富县植被覆盖率较高,NDVI位于−0.2~0.2之间的变化在原始影像上的目视解译不太明显,因此本研究将$ \beta $设为0.2,可保留$ \beta $大于0.2的断点,以便于检测较为显著的变化。算法输出变化开始时间(即断点位置)、变化结束时间、变化幅度和变化持续时长。

    • 本研究利用一元线性回归方法来拟合每个像元年内NDVI最大值(INDV(max))的年际变化趋势[24]

      $$ {S_{{\rm{NDVI}}}} = \frac{{n \displaystyle\sum\limits_{i = 1}^n {\left( {iI_{{\rm{NDV}}\left( {{\rm{max}}} \right)}^i} \right)} - \displaystyle\sum\limits_{i = 1}^n i \sum\limits_{i = 1}^n {I_{{\rm{NDV}}\left( {{\rm{max}}} \right)}^i} }}{{n \displaystyle\sum\limits_{i = 1}^n {{i^2}} - {{\left( {\displaystyle\sum\limits_{i = 1}^n i } \right)}^2}}} 。 $$

      其中:i为年序号;n为总年数;INDV(max)为年内NDVI最大值;$ {I}_{\mathrm{N}\mathrm{D}\mathrm{V}\left(\mathrm{m}\mathrm{a}\mathrm{x}\right)}^{i} $为第i年的INDV(max)$ {S}_{\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}} $为NDVI年内最大值的年际变化趋势。将基于一元线性回归方法得到的富县2000—2019年NDVI时间序列的趋势定义为总体趋势。依据宋怡等[25]提出的分类标准将趋势划分为3种类型:退化(${S}_{\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}}$≤−0.001 0)、基本不变($-0.001\;0 < {S}_{\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}}$≤0.001 0)、改善(${S}_{\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}} > 0.001\,0$)。

    • ① 对NDVI时间序列进行断点检测,分析断点及NDVI变化幅度的时空分布。② 针对未检测到断点或NDVI变化幅度小于0.2的区域,研究其总体趋势;针对断点NDVI变化幅度大于0.2的区域,对其进行多阶段趋势分析。多阶段趋势分析过程如下:以NDVI时间序列的所有断点及对应的变化结束时间对其进行分段,分别进行拟合,从而得到趋势序列。断点及对应的变化结束时间生成的拟合段定义为变化段,除了变化段以外的拟合段定义为趋势段。当前阶段趋势定义为最末一段趋势段。趋势类型按2.2.1中标准分类。对于多趋势段区域,研究其当前阶段趋势,以分析植被的现状。对于持续时长小于5 a的趋势段,由于时间太短,不进行趋势分析。③ 根据林地类型分析林地总体趋势和当前阶段趋势的变化情况。

    • 应用DBEST算法对富县NDVI时间序列数据进行断点检测,时间序列中第1断点NDVI变化幅度的空间分布与不同幅度面积占比如图3所示。发生正向变化的区域占全县面积的79.34%;负向变化的区域占20.66%,多位于道路两侧及居民点附近(图1)。NDVI变化幅度为−0.2~0.2的区域占总面积66.47%,主要位于西部及东北地区,植被发生的变化较小;变化幅度绝对值超过0.2的区域占总面积的33.53%,位于富县的东部及河流道路周围区域,属于NDVI剧烈变化的区域。

      图  3  富县归一化植被指数变化幅度及其面积占比      

      Figure 3.  Change magnitude and area percent of NDVI in Fu County      

      图4显示:2000—2019年间富县NDVI变化幅度大于0.2的区域断点个数的空间分布及各类型的面积占比。其中,66.47%的区域断点NDVI变化幅度小于0.2,主要分布在富县西部和东北地区,说明NDVI变化幅度小,植被相对稳定,因此不作进一步分析。断点NDVI变化幅度大于0.2的区域占总面积的33.53%,其中断点数量为1或2个的区域占总面积的21.59%,位于东部地区;断点数量多于4个的区域,仅占5.88%,集中分布于道路、河流沿线,植被变化频繁,与人为活动有关。

      图  4  富县不同NDVI变化幅度小于0.2与大于0.2的断点数目及其面积占比

      Figure 4.  Number of breakpoints and area percent of NDVI (<0.2 and >0.2) in Fu County

    • 富县NDVI总体趋势有很强的空间差异。对趋势值进行分类,得到不同类别的分布及面积(图5)。20 a间富县超过97%的区域植被呈改善趋势;退化区域占比不足2%,植被改善的面积远大于退化的面积。退化区域集中分布在富县东部区域、道路及河流沿线,在西北部散落分布,其中东部部分人造地表退化面积较大;基本不变的区域面积最小,不足1%,分布较零散。

      图  5  2000—2019年富县NDVI时间序列总体趋势类别及其面积占比

      Figure 5.  Overall trend level and area percent of NDVI time series in Fu County from 2000 to 2019

    • 基于3.1断点检测结果,对2000—2019年富县NDVI时间序列进行多阶段趋势分析。由于66.47%的区域无断点,因此该区域只分析总体趋势;33.53%的区域具有1个及以上断点,对该区域时间序列进行分段,分析多阶段趋势中的当前阶段趋势。当前阶段趋势统计结果及分布(图6A)表明:27.57%的区域当前阶段趋势持续时长大于5 a,植被当前阶段改善的面积远大于退化的面积,改善面积占24.54%,分布于东部区域及道路河流沿线的周边区域;退化面积仅占2.12%,基本不变面积不足退化面积的一半,分布较零散。当前阶段趋势的开始时间(图6B)均在2014年之前,时间分布跨度大;开始于2002年的区域面积最大,占12.76%,且均为改善趋势,表明经历一次变化以后,区域植被一直呈现改善状态;其次是开始于2014年的区域,占3.68%,主要位于道路、河流沿线;其余年份所占面积均不足2%。

      图  6  2000—2019年富县NDVI时间序列当前阶段趋势及其面积占比

      Figure 6.  Current stage trend and area percent of NDVI time series in Fu County from 2000 to 2019

    • 图2所示:富县植被退化的区域大部分分布在西北、西南、东北以及中部的一些地区,如今中部地区以及一些道路沿线已经变为永久性的建设用地;其他大部分均已改善,大多变为有林地和疏林地。2019年富县林地类型中,有林地占63.72%,疏林地占1.98%,灌木林地占13.27%,未成林造林地占6.07%,4种林地占富县总面积的85.04%。如表1所示:有林地、疏林地、灌木林地、未成林造林地的改善面积占其总面积的比例分别为96.98%、94.32%、94.41%、92.97%,其他用地改善和基本不变的比例分别为86.66%和2.50%,退化的比例达10.84%,由于其他用地包括建设用地,说明大部分退化与城市和道路建设相关。整体来说,富县生态建设工程取得了良好的效果,大部分林地处于改善的趋势,特别是有林地和灌木林地。

      表 1  富县林地类型趋势情况统计

      Table 1.  Statistics on trend of forest land types in Fu County

      林地类型改善/%基本不变/%退化/%
      有林地96.980.862.16
      疏林地94.321.913.77
      灌木林地94.411.813.78
      未成林造林地92.971.915.12
      其他用地86.662.5010.84
    • 本研究基于断点检测的多阶段趋势分析方法,研究了陕西省富县近20 a植被的时空变化特征,得到以下结论:①第1断点检测结果表明:富县发生正向变化的区域面积大,占79.34%。检测结果中NDVI变化幅度小于0.2的区域占66.47%,主要分布在富县西部和东北地区,植被相对稳定;大于0.2的区域占33.53%,其中断点数量为1或2个的区域比例为21.59%,位于东部地区;断点数量多于4个的区域,仅占5.88%,集中分布于道路、河流沿线,变化频繁,与人为活动有关。②植被的多阶段趋势表明:富县共有27.57%的区域当前阶段趋势持续时长大于5 a,当前阶段趋势为改善的面积占24.54%,分布于东部区域及道路河流沿线的周边区域;为退化的面积占2.12%,基本不变面积仅占0.91%,分布均较零散。当前阶段趋势的开始时间均在2014年之前,时间分布跨度大,空间异质性强,揭示了富县植被变化的多样性及复杂性。③近年来富县有林地、疏林地、灌木林地、未成林造林地均是改善情况明显,其中有林地和灌木林地中的改善比例最大。林地中呈改善态势的面积占全县总面积的81.83%,富县生态建设工程取得了良好的效果。同时对于检测到退化区域的林地,需要进一步分析,采取有效措施提高质量。

      1999年以来,富县积极响应国家号召,实施退耕还林、荒山荒地造林、封山育林等措施。这些生态工程建设促进了植被的恢复,然而在不同地区的效果是存在差异的。本研究获得的NDVI的变化直接反映了植被恢复的效果,但DBEST算法参数设置对分析结果的影响、趋势变化与影响因素的关系尚不清楚,在后续研究中,需从多尺度应用DBEST算法,探究其应用技术,并综合考虑气温、降水、人为因素等对趋势变化的影响。

参考文献 (25)

目录

    /

    返回文章
    返回