海南岛植被覆盖变化驱动因子及相对作用评价
1.
2.
Relative roles of climate changes and human activities in vegetation variables in Hainan Island
1.
2.
责任编辑: 陈 理
收稿日期: 2019-01-4 修回日期: 2019-04-1 网络出版日期: 2020-03-15
基金资助: |
|
Received: 2019-01-4 Revised: 2019-04-1 Online: 2020-03-15
作者简介 About authors
罗红霞(1985-),女,硕士,助理研究员,主要从事生态环境遥感监测研究。Email:120081008@163.com。 。
为了解海南岛植被覆盖变化影响因子,利用MODIS归一化植被指数(normalized difference vegetation index,NDVI)数据,采用残差分析法、趋势分析法和相对作用分析法对海南岛2001—2015年间植被覆盖变化特征进行研究,定量分析了气候变化和人类活动对海南岛植被覆盖变化的贡献。研究表明,2001—2015年间,海南岛植被总体呈现增加趋势,其增长速率为0.024/10 a,植被覆盖增加区域占全岛的77.77%,植被退化区域所占比例仅为22.23%。在植被改善区域,气候变化和人类活动的相对作用分别为31.04%和68.96%,其中相对作用大于50%的区域(即以人类活动为主)占整个植被改善区域的80.79%; 植被退化区域,气候变化和人类活动的相对作用分别为35.03%和64.97%,其中相对作用大于50%的区域(即以人类活动为主)占整个植被退化区域的75.59%。整体来看,气候变化对海南岛植被生长有促进作用,人类活动对海南岛地区植被的建设作用大于破坏作用。
关键词:
Using MODIS normal difference vegetation index (NDVI) data from 2001 to 2015 and climate data acquired from 7 meteorological stations in Hainan Island, the authors analyzed vegetation variability under the influence of climate variations and human activities based on residual analysis method, unary trend curve regression model and relative role analysis method in Hainan Island. The results showed that NDVI observed by remote sensing (MODIS NDVI) increased over the past 15 years with a rate of 0.024 per ten years in Hainan Island. The proportions of increasing vegetation and decreasing vegetation areas were 77.77% and 22.23%, respectively. As identified by remote sensing observations in Hainan Island, the relative roles of climate changes and anthropogenic activities in vegetation increase areas were 31.04% and 68.96%, while the roles of climate variations and human activities in vegetation decrease areas were 35.03% and 64.97%, respectively, indicating that human activities played a major role in vegetation changes. The area of increased vegetation mainly influenced by human activities (relative role >50%) accounted for 80.79% of the whole increase vegetation area, which was associated with the large-scale rubber plantation in these areas. In contrast, the area of reduced vegetation mainly induced by human activities (relative role > 50%) accounted for 75.59% of the reduced vegetation regions, which may be induced by the increased urbanization and urban expansion in coastal regions. On the whole, climate changes could promote the vegetation growth, whereas human activities played a greater role in the vegetation increasing than in the vegetation decreasing.
Keywords:
本文引用格式
罗红霞, 戴声佩, 李茂芬, 李玉萍, 郑倩, 胡盈盈.
LUO Hongxia, DAI Shengpei, LI Maofen, LI Yuping, ZHENG Qian, HU Yingying.
0 引言
植被变化受气候变化和人类活动的共同作用。遥感技术的快速发展极大地促进了大尺度植被动态变化的研究进程,国内外许多学者开展了植被变化的驱动机制研究[6,7,8,9]。针对海南岛植被变化的研究,章明等[10]基于SPOT VEGETATION数据分析了2001—2010年间海南岛植被覆盖时空变化特性; 李伟光等[11]利用MODIS增强植被指数产品分析了海南岛热带雨林、橡胶林及农田植被覆盖变化特征,指出降水是橡胶林的主要影响因子,而雨林则受温度影响; 林培松等[12]研究发现人类活动是导致海南岛西部沙漠化的主要原因。总体来看,对海南岛植被变化驱动因子的研究,针对气候变化分析较多,而人类活动对植被变化的影响研究侧重于定性描述,定量化研究较少。基于此,本文采用残差分析法、趋势分析法和相对作用分析法开展海南岛植被覆盖变化特征分析,分离气候变化和人类活动的贡献,定量评价两者对海南岛植被覆盖变化的相对作用大小,为海南岛生态环境相关政策制定提供科学支撑。
1 研究区概况及数据源
1.1 研究区概况
海南岛位于中国最南端,地理坐标为N18°10'~20°10',E108°37'~111°05',北隔琼州海峡与广东省雷州半岛相望,是我国仅次于台湾岛的第二大岛。海南岛植被茂密,森林覆盖率超过60%; 岛上耕地主要分布于四周沿海平原区域,中部和南部地区则主要以林地和园地为主,其中在儋州市、琼中黎族苗族自治县、白沙黎族自治县和澄迈县等地多为橡胶林[13]。海南岛属于热带季风气候区,年平均气温在22~26 ℃之间,全岛年均降水量在1 600 mm以上,降水空间分布呈现东多西少的特征,中部和东部年均降水量为2 000~2 400 mm,西南部沿海年均降水量为1 000~1 200 mm; 由于热带风暴或台风的影响,70%~90%的降水发生在5—10月之间[14]。研究区土地利用类型及数字高程模型(digital elevation model,DEM)如图1所示。
图1
1.2 数据来源
本文使用的遥感数据为2001—2015年的MOD13Q1 NDVI数据(空间分辨率为250 m,时间分辨率为16 d),行列号为h28v06和h28v07,数据来源于美国国家航空航天局; 已进行了消除水汽、云汽、重气溶胶影响等处理。为进一步消除异常值的影响,基于ArcGIS10.0软件采用最大值合成法(maximum value composites,MVC)合成年最大NDVI数据,最终形成2001—2015年间海南岛逐年NDVI数据集。
表1 各气候指标与NDVI的相关系数
Tab.1
气候指标 | 相关系数 |
---|---|
年均气温 | -0.253*① |
年均降水量 | 0.667** |
干燥度指数 | -0.196 |
生物温暖指数 | -0.253 |
生物干湿度指数 | 0.559* |
①: **和*分别表示相关系数通过了P<0.05和P<0.01显著性检验。
图2
图2
主要气候指标与NDVI的空间相关性
Fig.2
Correlation coefficients between main climate factors and NDVI
表1结果表明,气候指标中年均气温、干燥度指数和生物温暖指数与NDVI呈负相关,其中年均气温与NDVI的相关系数在0.01置信水平上达到显著水平; 年均降水量和生物干湿度指数与NDVI呈正相关,相关系数分别为0.667和0.559,分别通过了P<0.05和P<0.01的显著性检验。因此,本文采用年均降水量和生物干湿度指数2个指标逐栅格像元计算NDVI预测值。
土地利用类型数据来源于中国生态系统与生态功能区划专题数据库(
2 研究方法
2.1 残差分析法
利用年最大NDVI和5个气候指标中相关性最大的2个指标做回归分析,逐栅格建立气候指标与NDVI之间的二元一次回归关系模型,得到海南岛2001—2015年间逐栅格像元的NDVI预测值。具体公式为
式中: NDVIp为NDVI预测值; β0,β1和β2分别为回归方程的待定系数; x1和x2分别为与NDVI相关性较大的2个气候指标变量,即年均降水量和生物干湿度指数。
式中: NDVIt为NDVI遥感观测值; ε为NDVI观测值与预测值的残差,当ε>0时,表示人类活动对植被生长有利; 当ε<0时,表示人类活动不利于植被生长; 当ε≈0时,表示人类活动对植被变化不造成影响。应用残差分析法分离气候变化对海南岛植被覆盖变化的影响,可以准确识别人类活动在植被覆盖变化过程中起到的作用。
2.2 趋势分析法
趋势分析法可以模拟每个栅格像元的变化趋势,通过逐像元时间变化特征反映整个研究区域空间变化规律以及区域时空格局演变。具体公式为
2.3 相对作用分析法
表2 不同情景下气候变化和人类活动在植被覆盖变化过程中的相对作用分析法①
Tab.2
植被 变化 | 情景 | Slope (NDVIC) | Slope (NDVIH) | 气候变化的相对作用/% | 人类活动的相对作用/% | 说明 | |
---|---|---|---|---|---|---|---|
植被增加区 | 1 | >0 | >0 | 气候变化和人类活动共同作用导致植被覆盖增加 | |||
2 | >0 | <0 | 100 | 0 | 气候变化导致植被覆盖增加 | ||
3 | <0 | >0 | 0 | 100 | 人类活动导致植被覆盖增加 | ||
植被减少区 | 1 | <0 | <0 | 气候变化和人类活动共同作用导致植被覆盖减少 | |||
2 | <0 | >0 | 100 | 0 | 气候变化导致植被覆盖减少 | ||
3 | >0 | <0 | 0 | 100 | 人类活动导致植被覆盖减少 |
①: Slope(NDVIC)为NDVI预测值变化趋势; Slope(NDVIH)为NDVI残差值变化趋势; CNDVI为2015年通过气候指标计算的NDVI预测值与2001年NDVI预测值的差值; HNDVI为2015年NDVI残差值与2001年NDVI残差值的差值。
3 结果与分析
3.1 海南岛植被覆盖时空变化特征
2001—2015年间,海南岛植被覆盖总体呈现增长态势,增速为0.024/10 a。从空间分布来看,植被覆盖增加区域(Slope>0)占全岛的77.77%,植被退化区域(Slope<0)仅占全岛的22.23%(图3和表3)。从图3可以看出,减少的区域主要分布在海口市北部、文昌市西南部、陵水黎族自治县东部和三亚市南部沿海区域,东方市、琼中黎族苗族自治县、乐东黎族自治县、琼海市和万宁市也有部分分布。增加的区域集中分布于儋州市、澄迈县、屯昌县、临高县南部和文昌市北部等区域。从表3中可以看出,各市县中儋州市植被增加区域所占比例最大(8.04%),而植被减少区域面积最大的是乐东黎族自治县和三亚市,所占比例分别为2.67%和2.40%。
图3
图3
海南岛2001—2015年植被覆盖变化趋势(气候变化和人类活动共同作用影响)
Fig.3
Spatial distribution of land cover change trend from 2001 to 2015 of Hainan Island
表3 海南岛各市县植被增加和减少区域所占比例
Tab.3
地区 | 整体变化影响 | 气候变化影响 | 人类活动影响 | |||
---|---|---|---|---|---|---|
植被增加区域 | 植被减少区域 | 植被增加区域 | 植被减少区域 | 植被增加区域 | 植被减少区域 | |
海口市 | 4.92 | 1.73 | 4.43 | 2.22 | 4.70 | 1.94 |
文昌市 | 5.19 | 1.67 | 4.90 | 1.91 | 4.94 | 1.87 |
琼海市 | 4.25 | 1.00 | 2.94 | 2.31 | 4.26 | 0.98 |
万宁市 | 4.43 | 1.30 | 3.15 | 2.58 | 4.40 | 1.34 |
陵水黎族自治县 | 2.38 | 0.91 | 1.98 | 1.29 | 2.31 | 0.97 |
三亚市 | 3.05 | 2.40 | 3.09 | 2.35 | 2.76 | 2.68 |
乐东黎族自治县 | 5.28 | 2.67 | 5.49 | 2.47 | 4.55 | 3.40 |
东方市 | 3.99 | 1.99 | 3.56 | 2.43 | 3.92 | 2.07 |
昌江黎族自治县 | 3.25 | 1.15 | 2.50 | 1.90 | 3.22 | 1.18 |
白沙黎族自治县 | 5.71 | 0.85 | 4.55 | 2.03 | 5.69 | 0.89 |
儋州市 | 8.04 | 0.89 | 7.53 | 1.38 | 7.65 | 1.26 |
临高县 | 3.42 | 0.32 | 3.08 | 0.66 | 3.36 | 0.38 |
澄迈县 | 5.42 | 0.82 | 4.85 | 1.40 | 5.28 | 0.97 |
定安县 | 3.22 | 0.48 | 1.85 | 1.86 | 3.27 | 0.44 |
琼中黎族苗族自治县 | 6.58 | 1.78 | 5.15 | 3.24 | 6.89 | 1.50 |
屯昌县 | 2.97 | 0.82 | 2.74 | 1.05 | 2.89 | 0.91 |
五指山市 | 2.98 | 0.51 | 2.09 | 1.42 | 2.98 | 0.53 |
保亭黎族苗族自治县 | 2.76 | 0.86 | 2.37 | 1.25 | 2.71 | 0.91 |
3.2 气候变化和人类活动对植被覆盖变化影响的空间特征
海南岛NDVI预测值时间序列趋势分析表明,NDVI预测值增速为0.004/10 a,该值远远小于NDVI观测值的变化趋势(0.024/10 a)。从图4(a)统计分析得出,海南岛地区气候变化导致的NDVI预测值变化呈增加趋势的面积占总面积的66.25%,NDVI预测值变化呈减少趋势的面积占33.75%,植被变化呈增加趋势的面积大于呈减少趋势的面积。总体来看气候变化对海南岛植被生长有促进作用。其中,NDVI预测值呈增加趋势的区域主要分布于海口市东部、琼中黎族苗族自治县中部、文昌市北部以及儋州市、临高县、澄迈县的大部分区域; NDVI预测值呈减少趋势的区域主要分布于乐东黎族自治县、三亚市、东方市西北部地区、琼海市西南部地区以及昌江黎族自治县西北部地区。海南岛NDVI残差趋势分析显示,NDVI残差序列变化趋势为0.020/10 a,该变化趋势与NDVI观测值的变化趋势(0.024/10 a)较为相近。15 a间的海南岛NDVI残差变化呈增加趋势的面积占总面积的75.41%,呈减少趋势的面积占总面积的24.59%(图4(b)); 总体来看,人类活动对海南岛地区植被的建设作用大于破坏作用。NDVI残差变化呈增加趋势的区域主要分布在儋州市中北部、临高县、澄迈县和定安县,这些区域人类活动对植被覆盖增加起到了积极的作用。NDVI残差变化呈减少趋势的区域主要分布在海口市北部、乐东黎族自治县中北部、三亚市、保亭黎族苗族自治县中南部、东方市东部、琼海市西南部地区以及琼中黎族苗族自治县、屯昌县、陵水黎族自治县、万宁市也有零星分布。
图4
图4
海南岛2001—2015年气候变化和人类活动对植被覆盖变化影响
Fig.4
Spatial distribution of the effects of climate changes and human activities from 2001 to 2015 of Hainan Island
3.3 气候变化和人类活动对植被覆盖变化的相对作用
通过以上分析,可以看出海南岛植被变化有些区域受气候变化大于人类活动,而在有些区域受人类活动影响大于气候变化,为了详细区分海南岛植被覆盖变化中气候变化和人类活动的作用,采用2.3节提及的相对作用分析法进行研究。
植被覆盖改善区域,气候变化和人类活动的相对作用空间分布如图5所示。气候变化的相对作用为31.04%,人类活动的相对作用为68.96%,表明2001—2015年间人类活动是导致海南岛植被改善的主要因素。统计分析显示,相对作用大于50%的区域(即以人类活动为主)占整个植被退化区域的80.79%,而相对作用小于50%的区域(以气候变化为主)占19.21%。空间分布上看,以人类活动影响为主的区域主要分布在儋州市、澄迈县、定安县、白沙黎族自治县和五指山市等地区。这可能与近年来政府出于发展经济和生态保护双重考虑,大规模发展橡胶种植有关[5]; 以儋州市和澄迈县为例,海南统计年鉴数据显示[23,24],2001年2个地区橡胶种植面积分别为18 200 hm2和5 905 hm2,2015年分别增加到89 101 hm2和53 622 hm2,种植面积分别增加了4倍和9倍。气候变化影响为主的区域分布于琼中黎族苗族自治县、屯昌县和保亭黎族苗族自治县等中部山区,可能因为中部山区森林茂密,降水多且雨季长,是水源涵养功能最重要的区域,2008年中部山区总产水量为34.29亿m3,占全省总用水量的77.88%[25],良好的气候条件对植被生长十分有利。
图5
图5
海南岛2001—2015年植被覆盖改善区气候变化和人类活动的相对作用
Fig.5
Spatial distribution of relative roles of climate changes and human activities for increasing vegetation areas in Hainan Island form 2001 to 2015
植被退化区域,气候变化和人类活动的相对作用空间分布如图6所示。气候变化的相对作用为35.03%,人类活动的相对作用为64.97%,表明2001—2015年间人类活动是导致海南岛植被退化的主要因素,这一结果与黄宝荣等[26]的研究一致。统计分析显示,相对作用大于50%的区域(即以人类活动为主)占整个植被退化区域的75.59%,而相对作用小于50%的区域(以气候变化为主)占24.41%。从空间分布上来看,人类活动影响为主的区域主要分布在海口市北部沿海、文昌市、陵水黎族自治县和万宁市东部沿海区域、三亚市、乐东黎族自治县和东方市大部分地区。这些区域人类活动对植被覆盖产生了负面作用,如海口市国际旅游岛建设、城镇化快速发展,造成城市用地过度扩张,建成区面积从2001年的34 km2增加到2015年的152.4 k
图6
图6
海南岛2001—2015年植被退化区气候变化和人类活动的相对作用
Fig.6
Spatial distribution of relative roles of climate changes and human activities for decreasing vegetation areas in Hainan Island form 2001 to 2015
3 结论
基于2001—2015年的MODIS NDVI数据,采用残差分析法、趋势分析法和相对作用分析法研究了海南岛不同地区气候变化和人类活动对海南岛植被覆盖变化的相对作用影响,得出如下结论:
1)2001—2015年间,海南岛植被覆盖总体呈现增长态势,增速为0.024/10 a。植被覆盖增加区域占全岛的77.77%,植被退化区域仅占全岛的22.23%。各市县中儋州市植被增加区域所占比例最大。
2)植被改善区域,气候变化的相对作用为31.04%,人类活动的相对作用为68.96%,相对作用大于50%的区域(即以人类活动为主)占整个植被改善区域的80.79%,而相对作用小于50%的区域(以气候变化为主)占19.21%。植被退化区域,气候变化的相对作用为35.03%,人类活动的相对作用为64.97%,相对作用大于50%的区域(即以人类活动为主)占整个植被退化区域的75.59%,而相对作用小于50%的区域(以气候变化为主)占24.41%。
3)采用回归模型预测NDVI值,选取的气候指标不同、回归模型的不同都会造成预测NDVI的差异,如何依据研究区域典型植被特点有针对性选取气侯指标,探讨不同数学模型精确模拟NDVI值变化,仍需进行大量的深入分析研究。
参考文献
Vegetation establishment on the southern Israeli coastal sand dunes between the years 1965 and 1999
[J].
Fractional vegetation cover estimation in arid and semi-arid environments using HJ-1 satellite hyperspectral data
[J].
2000—2010年陕北地区植被NDVI时空变化及其与区域气候的关系
[J].
Spatio-temporal change of vegetation NDVI and its relations with regional climate in northern Shaanxi Province in 2000—2010
[J].
海南省生态旅游发展的SWOT分析
[J].
Analysis of eco-tourism development based on SWOT on Hainan Island
[J].
海南岛主要森林类型时空动态及关键驱动因子
[J].
The dymics of spatial and temporal changes to forested land and key factors driving change on Hainan Island
[J].
Global correlation analysis for NDVI and climatic variables and NDVI trends:1982—1990
[J].
Drivers of greening trend across vertically distributed biomes in temperate arid Asia
[J].
2000—2014秦巴山区植被覆盖时空变化特征及其归因
[J].
Spatial temporal variation of vegetation coverage in Qinling-Daba Mountains in relation to environmental factors
[J].
2001—2014年黄土高原植被覆盖时空变化特征及其归因
[J].
Spatiotemporal analyses and associated driving forces of vegetation coverage change in the Loess Plateau
[J].
基于SPOT VEGETATION数据的海南岛年际植被变化研究
[J].
Study of annual vegetation variability in Hainan Island based on SPOT VEGETATION data
[J].
海南岛典型植被区EVI特征及其对气象因子的响应
[J].
Typical vegetation area EVI characteristics and responses to meteorological factors in Hainan Island
[J].
近20 a来海南岛西部土地沙漠化与气候变化关联度研究
[J].
Correlativity between land desertification and climate variability in west of Hainan Island during past nearly 20 years
[J].
基于PALSAR雷达数据与多时相TM/ETM+影像的海南土地利用分类研究
[J].
Land utilization mapping in Hainan Island by using ALOS PALSAR and multi-temporal Landsat TM/ETM+ imagery
[J].
Tropical birds are declining in the Hainan Island of China
[J].
The relative roles of climate variations and human activities in vegetation change in North China
[J].
KIRA指标的拓展及其在中国植被与气候关系研究中的应用
[J].
Development of KIRA’s indices and its application to vegetation climate interaction study of China
[J].
黄土高原植被覆盖变化与气候和人类活动的关系
[J].
Vegetation cover,climate and human activities on the Loess Plateau
[J].
三江源地区人类活动对植被覆盖的影响
[J].
Influence of human activities on variation of vegetation cover in the Three-River Source Region
[J].
基于遥感的植被年际变化及其与气候关系研究进展
[J].
Advance in the inter-annual variability of vegetation and its relation to climate based on remote sensing
[J].
近10年陕甘宁黄土高原区植被覆盖时空变化特征
[J].
The spatial-temporal changes of vegetation restoration on Loess Plateau in Shaanxi-Gansu-Ningxia region
[J].
气候变化和人类活动在鄂尔多斯地区沙漠化过程中的相对作用研究
[J].
Assessing the relative role of climate change and human activities in sandy desertification of Ordos region,China
[J].
华北地区植被覆盖变化及其影响因子的相对作用分析
[J].
Analysis of the vegetation cover change and the relative role of its influencing factors in North China
[J].
海南岛中部山区生态系统水源涵养功能研究
[J].
Study on the water conservation function of the ecosystem in the central mountainous area on Hainan Island
[J].
海南岛生态环境脆弱性评价
[J].
Assessment of eco-environmental vulnerability of Hainan Island,China
[J].
城市用地扩张影响因素及其极限规模测算——以海口市为例
[J].
The influence factors of urban land expansion and calculation of its utmost scale:A case study of Haikou City
[J].
三亚市土地利用/覆被变化及其驱动机制研究
[J].
Land use and land cover change and its driving factors in Sanya
[J].
RS与GIS支持下近50 a海南岛西部土地沙漠化时空演变过程研究
[J].
Spatio-temporal changes of sandy desertification in recent 50 years in western Hainan Island based on RS and GIS
[J].
雨影区对海南岛西部土地沙漠化影响的风洞实验研究
[J].
Study of influence of rain shadow region on sandy desertification in west of Hainan Island based on wind tunnel experiments
[J].
海南西部干旱的地理边缘解析
[J].
Analysis of drought in western Hainan from the perspective of geographical fringes
[J].
基于SVM的湖泊咸度等级遥感信息提取方法——以内蒙古巴丹吉林沙漠为例
[J].
Extraction of remote sensing information for lake salinity level based on SVM:A case from Badain Jaran desert in Inner Mongolia
[J].
/
〈 |
|
〉 |
