2001～2021 年海南岛植被覆盖时空动态变化与驱动力分析

张阳阳，刘铁冬，2，龚文峰，2，段璇瑜，孙雨欣，邱梓轩，2

（1. 海南大学林学院，海南 海口 570228；2. 海口市智慧林业重点实验室，海南 海口 570228）

植被覆盖度（fractional vegetation cover，FVC）是描述生态系统特征的基础数据［1-2］，是反映区域植被生长的重要指标［3］，是研究区域植被的重要参数。探究长时序的植被覆盖变化对研究中观尺度的生态格局变化具有重要意义；对区域的生态保护与修复具有重要的参考意义，探究区域植被覆盖时空动态变化的驱动力，对区域的生态文明建设、城市和谐发展具有重要意义［4-6］。

近年来，许多学者从不同的角度出发对不同地区的长时序植被覆盖变化特征进行了研究，并取得了一定的成果。刘可等［7］利用1981～2012 年的GIMMS NDVI3G数据研究了中国范围内各生态系统植被活动对气候变化的响应规律，指出近30 a 来我国NDVI 总体呈非稳定性上升；刑愿等［8］利用MODIS NDVI 数据分析了贵州省2001～2018 年不同土地覆盖类型的NDVI 的变化特征，指出在植被生长季，林地的NDVI 值最大。在植被覆盖度变化驱动力的研究中，多以自然因素的研究为主，如江萍等［9］使用SPOT/VEGETATION NDVI 数据分析了新疆1998～2018年植被NDVI时空变化特征和驱动力，指出水分是影响新疆植被生长的主导因素；张圳等［10］利用MODIS NDVI 数据分析了太湖2000～2015 年水生植被NDVI 的时空特征，指出不同植被类型NDVI变化的主导气候因素不同；Zhang 等［11］利用MODIS NDVI 数据分析了2000～2018 年三江源地区植被覆盖的时空变化特征，指出气候变化是影响该地区植被变化的主要因素。目前针对海南岛植被覆盖的研究较少，对影响植被覆盖因素的探究也相对较少。罗红霞等［12］通过研究指出气温对海南岛植被覆盖的影响强于降水，但未深入分析人为活动因素对植被覆盖的影响。

海南岛是我国热带季雨林和热带雨林的分布地之一，是热带雨林主要区和典型代表区，2021 年，海南热带雨林国家公园被列入中国首批国家公园名单。随着海南热带雨林国家公园的成立，海南岛也将成为我国热带雨林生态保护和修复的重点区域。因此，本研究探究海南岛植被覆盖的时空变化特征及驱动因素，分析海南岛植被覆盖未来演变趋势，研究人类活动与植被覆盖变化的相关性，以期能为海南岛生态保护与修复工作提供理论依据和参考。

1 材料与方法

1.1 研究区概况

海南岛位于中国南部（N 18°10′～20°10′，E 108°37′～111°05′），中间高，四周低，以中部山区为核心，向四周逐级下降，陆地面积3.44 万 km2，处于热带，属于热带季风性气候，年日照时数1 793～2 590 h，年平均气温22.5～25.6 ℃，各地的年降雨量923～2 459 mm。

1.2 数据来源与处理

1.2.1 遥感数据 本文所使用的数据来自美国国家航天局（national aeronautics and apace administration，NASA）发布的MOD13Q1 植被指数产品，时间分辨率为16 d，空间分辨率为250 m，时间序列为2001～2021 年。 使用MRT（MODIS reprojection tool）工具进行影像数据的镶嵌和地图投影转换，利用ArcGIS10.6结合海南岛的矢量地图数据得到海南岛区域的NDVI数据来计算海南岛的植被覆盖度，用最大值合成法（maximum value composite，MVC）获取海南岛2001～2021 年每年的植被覆盖度［13］。

图1 海南岛行政区划图Figure 1 Administrative zoning map of Hainan Island

1.2.2 相关统计数据 人类活动数据来自于2002～2021年《海南统计年鉴》。实际获取的数据范围为2001～2020年，选取国内生产总值（亿元）、人均国内生产总值（万元/人）、年末常住人口（万人）、城镇化率（%）、接待游客总人数（万人）和累计造林面积（万hm²）作为研究用的人类活动数据。

1.3 研究方法

1.3.1 植被覆盖度计算 植被覆盖度（FVC）与NDVI存在极显著的线性相关关系，采用像元二分模型可将NDVI数据转换为植被覆盖度［14-15］，计算公式如下：

式中：NDVIveg为纯植被像元的NDVI值；NDVIsoil为纯裸地像元的NDVI值。根据频率统计表，选取累计频率为0.5%的NDVI值作为NDVIsoil，选取累计频率为99.5%的NDVI值作为NDVIveg。

1.3.2 Theil-Sen Median 斜率估计与Mann-Kendall 显著性检验法 Theil-Sen Median 斜率估计是常用的非参数统计研究变化趋势的方法，适用于检验长时间序列的植被覆盖度变化趋势［16-17］，计算公式如下：

式中：FVCj和FVCi分别为第j年和第i年植被覆盖度；Median为中值函数。若Slope＜0，表示在该时间段内植被覆盖度呈下降趋势；若Slope＞0，表示在该时间段内植被覆盖度呈上升趋势。Slope的绝对值越大，植被覆盖度下降或上升的趋势越明显。

Mann-Kendall显著性检验法是常用的非参数趋势检验的工具之一，可用于长时间序列下的数据变化趋势检验［18-19］。本次研究所使用的是修正后的Mann-Kendall显著性检验法，排除了数据间的自相关性，提高了检验结果的准确性［20］。通过计算P值进行显著性检验，计算公式如下：

对于任意待检序列Xt（t=1，2，…，n），可定义统计量S：

式中：xi和xj为时间序列相应数据；sgn（xj-xi）为符号函数。

当n≥0时，统计量S近似服从正态分布，其期望为0，方差为：

式中：ρs(i)为观测值秩的自相关函数。

构造标准化的检验统计量Z：

P值：

若P＜0.01为极显著，0.01≤P＜0.05为显著，若P≥0.05 为不显著。结合Slope值与P值，将植被覆盖度变化趋势分为：极显著性减少、显著性减少、不显著性减少、不显著性增加、显著性增加、极显著性增加。

1.3.3 变异系数 变异系数（CV）能反映数据的相对离散程度［21-22］，通过计算植被覆盖度的变异系数，可以研究植被覆盖度变化的稳定性情况，计算公式如下：

式中：σ为植被覆盖度的标准差为植被覆盖度的平均值。CV值越大，植被覆盖度年际变化越大；反之，植被覆盖度年际变化越小。根据海南岛植被覆盖度稳定性的实际情况，可将稳定性等级分为：低波动变化、相对较低波动变化、中度波动变化、相对较高波动变化、高波动变化。

1.3.4 Hurst 指数 Hurst 指数是定量描述时间序列信息长期依赖性的主要方法之一，广泛应用于气象、水文等领域［23-26］。本文利用Hurst指数定量描述海南岛植被覆盖度长时间序列下的持续性或反持续性，研究海南岛植被覆盖度未来的演变趋势，计算公式如下：

存在时间序列｛ζ（t）｝，t=1，2，…，对于任意正整数τ≥1，定义其均值序列为：

累计离差序列：

极差序列：

标准差序列：

如果关系R（τ）/S（τ）ατH成立，则说明时间序列｛ζ（t）｝，t=1，2，…存在Hurst 现象，式中H值即为Hurst 指数，取值范围为（0，1），当0＜H＜0.5 时，该时间序列未来的变化趋势与过去的变化趋势呈负相关；H=0.5时，该时间序列未来的变化趋势与现在的变化趋势无关；0.5＜H＜1时，该时间序列未来的变化趋势与过去的变化趋势呈正相关［27］。结合Slope值和H值预测海南岛植被覆盖度未来的演变趋势，可将植被覆盖度未来的演变趋势分为：强持续性增加、强持续性减少、弱持续性增加、弱持续性减少、弱反持续性增加、弱反持续性减少、强反持续性增加、强反持续性减少。

在MATLAB R2018b 中，利用海南岛近21 a 每年的植被覆盖度的栅格数据，根据公式，分别计算Slope值、P值和H值，并在ArcGIS10.6 中进行统计分析。

1.3.5 相关分析和主成分分析 在R 3.6.1软件中用Pearson相关性分析的方法，确定海南岛植被覆盖度和各人类活动因素之间的相关性，并通过主成分分析（PCA），进一步确定主要的影响因素。

2 结果与分析

2.1 植被覆盖度动态变化特征

2.1.1 空间分布特征 海南岛的植被覆盖度整体较高，多年平均植被覆盖度为0.87。其中，年均植被覆盖度大于0.8 的区域占总面积的87.77%，集中分布在海南岛的中部山区，大部分地区属于热带雨林国家公园范围内，属于海南岛的自然保护区域，人为因素影响小，植被的生长状况较好；年均植被覆盖度小于0.4的区域占总面积的0.09%，分布于沿海的鱼塘、沙滩区域和正在进行开发建设的区域；年均植被覆盖度介于0.4～0.8的区域主要位于城市区域和城乡交界处，占总面积的12.14%（图2）。

图2 海南岛2001～2021年植被覆盖度年均值空间分布图Figure 2 Spatial distribution of annual mean FVC in Hainan Island from 2001 to 2021

2.1.2 时序变化特征 海南岛植被覆盖度2001年为0.852，2021年为0.881，近21 a植被覆盖度整体呈上升趋势，上升的速率为0.001 6/a（图3），增长率为3.40%。2006 年增长速率最快，海南加大生态环境建设力度，积极建设自然保护区，加大造林的力度，全省造林面积3.12 万 hm²，森林覆盖率达55.6%。2002～2005、2009～2010、2013～2014、2017～2019年4 个时段植被覆盖度呈下降趋势，2005 年的植被覆盖度为近21 a来的最低值。其中，2004～2005年，出现了全岛大范围的春冬连旱现象，全年降水明显减少［28］，造成植被覆盖度严重减少；2009～2010年减少速率最快，可能是因为国际旅游岛的建设促进了旅游业和房地产业的快速发展，建筑用地面积的快速增长，导致植被覆盖度减少速率加快。

图3 海南岛2001～2021年植被覆盖度年际变化趋势Figure 3 Interannual trend of FVC in Hainan island from 2001 to 2021

2.2 植被覆盖度动态变化规律

海南岛2001～2021年植被覆盖度主要呈增加趋势，增加区域占总面积的81.15%，主要增加区域集中分布在儋州、澄迈中部和临高南部（图4）。其中，极显著性增加和显著性增加的区域占总面积的47.54%，主要分布在儋州、临高、澄迈、白沙、五指山、琼中西部、乐东北部和昌江南部，这可能是因为2002年“退耕还林”政策实施以来，该地区的林地面积显著增加，使该地区的植被覆盖度亦随之增加（图5）。

图4 海南岛2001～2021年植被覆盖度年际变化趋势空间分布图Figure 4 Spatial distribution of interannual trend of FVC in Hainan Island from 2001 to 2021

图5 海南岛2001～2021年植被覆盖度显著性变化空间分布图Figure 5 Spatial distribution of significant change of FVC in Hainan island from 2001 to 2021

减少区域占总面积的18.85%，重点分布在各市县的沿海区域和城市建设区域。其中，显著性减少和极显著性减少的区域共占总面积的5.09%，主要分布在海口北部、陵水南部、文昌西部沿海区域、乐东南部和西部、三亚南部和中部，这些区域目前多为城市开发区和建成区，受人为干扰影响较大，植被覆盖度各类变化程度所占的比重如表1所示。

表1 植被覆盖度变化显著性统计Table 1 Statistical table of significance of FVC change

2.3 植被覆盖度稳定性分析

海南岛2001～2021 年植被覆盖度变异系数（CV）的平均值为0.04，整体波动幅度较低。低波动变化主要集中分布在中部山区，占总面积的48.00%；相对较低波动变化分布在城市和乡镇的开发边界外，占总面积的43.40%；中度波动变化、高度波动变化和相对较高波动变化的区域分布在人为活动频繁且开发较大的沿海区域、城镇开发区域和城市建设区域（图6）。由此推断，海南岛的植被覆盖受人为因素的影响可能大于自然因素的影响。植被覆盖度稳定性各类所占的比重如表2所示。

表2 植被覆盖度稳定性统计Table 2 Statistical table of FVC stability

图6 海南岛2001～2021年植被覆盖度稳定性空间分布图Figure 6 Spatial distribution of FVC stability in Hainan island from 2001 to 2021

2.4 植被覆盖度演变趋势预测

海南岛2001～2021 年植被覆盖度Hurst指数的平均值为0.47，以弱反持续性为主，占总面积的65.20%（图7）。未来植被覆盖度以减少为主，反持续性增加和持续性减少的面积占总面积的60.82%，其中，弱反持续性增加占总面积的53.56%，该区域绝大部分的植被覆盖度较高；持续性增加和反持续性减少的面积占总面积的39.18%，以弱反持续性减少和弱持续性增加为主，在空间上与弱反持续性增加交错分布（图8）。植被覆盖度各类变化趋势持续性所占的比重如表3所示。

表3 植被覆盖度变化趋势持续性统计表Table 3 Statistical table of FVC change trend persistence

图7 海南岛植被覆盖度Hurst指数空间分布Figure 7 Spatial distribution of FVC Hurst index in Hainan Island

图8 海南岛植被覆盖度未来变化趋势预测Figure 8 Prediction of future FVC change trend in Hainan Island

总体来说，海南岛的植被覆盖度在未来有60.82%的区域呈减少的趋势，这是人为活动与气候变化共同影响的结果，人类活动范围的增加和全球极端天气频繁爆发都将导致植被覆盖度的减少。因此，若不通过人为活动进行积极引导，加强对生态环境的保护与修复，海南岛整体植被覆盖度水平将面临降低的风险。

2.5 人类活动对植被覆盖度影响分析

通过相关性分析，各人类活动因素与海南岛植被覆盖度的相关性强弱顺序依次为：人均国内生产总值＞国内生产总值＞累计造林面积＞年末常住人口＞城镇化＞接待游客总人数（图9～10）。

图9 植被覆盖度与人类活动因素相关性分析Figure 9 Correlation analysis between FVC and human activity factors

图10 人类活动因素主成分分析碎石图Figure 10 Human activity factors principal component analysis gravel map

人均国内生产总值与植被覆盖度的相关系数为0.904，呈极显著正相关（P＜0.01），随着人均国内生产总值的提高，人民的物质生活水平随之提高，对生态文明建设的力度也随之加大。累计造林面积与植被覆盖度相关系数为0.887，呈极显著正相关（P＜0.01）。城镇化与植被覆盖度的相关系数分别为0.847，呈极显著正相关（P＜0.01）。随着城镇化的提高，人们也越来越重视城市绿地的建设，城市绿地系统建设工作也取得重大成效，2020 年海南省城市绿地面积共223.28 km2。主成分分析表明，人均国内生产总值和累计造林面积为主要的驱动因素，二者的累计贡献率为99.3%。

3 讨论

海南岛近21 a 的植被覆盖度呈上升趋势，这与罗红霞等［29］和李伟光等［30］的研究结果一致。本文回归分析得到的植被覆盖上升的速率低于前人的研究结果，可能是因为本文研究所用的NDVI 的时间序列更长，在研究过程中将NDVI 数据转化为植被覆盖度，减少了异常值的影响，更能体现植被覆盖度变化的规律。

海南岛2001～2021年植被覆盖度的变化趋势以增加为主。近年来海南省坚持生态立省，积极进行植树造林和退耕还林工程，2001～2021年，全省累计人工造林50.55万 hm²。在城市建设过程中，人们的生态保护意识不断加强，对城市绿地系统的建设也更加重视，目前海南省城市建成区的绿化覆盖率达到了41.10%。同时，植被覆盖度呈减少趋势的面积占总面积的18.85%，大部分面积集中在沿海区域，人类的开发利用是造成沿海区域植被覆盖度减少的重要原因之一。研究表明，近年来，海南岛海岸线利用程度综合指数大幅度提高［31］。未来植被覆盖的变化以减少为主要趋势，人为活动的干扰是植被覆盖度减少的原因之一，自然灾害、气候变化等在一定程度上也减少了区域的植被覆盖度。

4 结论

1） 海南岛近21 a植被覆盖度总体呈上升趋势，植被覆盖度变化趋势具有明显的空间异质性，全岛植被覆盖度主要呈增加趋势，以不显著性增加为主。海南岛西北部增加趋势较大，以极显著性增加为主；东北部和西南部区域增加不明显，且城市建设区域和部分沿海区域呈负增加；中部和东部地区植被覆盖度增加不明显，仍保持在较高的水平。

2） 海南岛2001～2021年的植被覆盖度最稳定的区域主要集中分布在中部和东部，西北部地区植被覆盖度变化波动中等，东北部和西南部的城市建设区、沿海区域植被覆盖度的变化波动较大。Hurst分析表明，海南岛60.82%面积的植被覆盖度有减少的趋势，以弱反持续性增加为主，主要分布在海口南部和东部、澄迈南部、儋州、白沙、昌江东南部、琼中、屯昌和定安。

3） 人均国内生产总值与植被覆盖度的相关性最强，累计造林面积和人均国内生产总值是影响植被覆盖度变化的主要人类活动因素。