随着粤港澳大湾区国家战略的提出及“一核一带一区”区域发展新格局的构建, 广东省的生态问题得到广泛关注。在粤港澳大湾区国家战略中, 广东省承担着大湾区生态保护屏障功能, 其森林覆盖率达59.08%^[1], 植被作为区域重要的生态屏障, 对区域的生态安全和可持续发展起着不可替代的作用。植被覆盖度(Fractional Vegetation Coverage, FVC)作为重要的生态气候参数^[2], 区域地表植被覆盖度的提取及变化分析, 对揭示生态系统环境变化、植被恢复和重建布局具有重要意义^[3]。近年来, 遥感数据因其区域性及时效性等优势日渐成为区域植被动态监测的重要数据^[4], 由遥感数据反演的植被指数(Ⅵ)可以定量表达植被的生长状况, 已被广泛用于植被动态的参数化表达^[5]。归一化差值植被指数(NDVI)和增强植被指数(EVI)是目前反映植被生长状况常用的两个植被指数^[6-7]。其中, NDVI与植被分布密度呈显著线性相关, 是研究植被覆盖度的重要参数^[8], 也是应用最为普遍的植被定量监测指标^[9]。何全军等^[10]选取RVI、NDVI和EVI构建广东省植被指数序列并进行相关植被指数的应用分析, 结果表明RVI和NDVI的相关性高且结果稳定, 是广东省植被指数序列构建的最佳选择, 而EVI具有数据异常现象。

目前, 基于遥感信息估算植被覆盖度可归纳为经验模型法和植被指数转换法。像元二分法是利用NDVI对植被覆盖度进行提取的一种方法, 很大程度上削弱了大气、土壤背景和植被类型的影响, 普适性较好^[11-12], 是当前应用最为广泛的植被指数法。前人基于此已取得了大量成果: 穆少杰等^[13]探索了2001—2010年内蒙古地区植被覆盖度的空间格局和变化规律；李林叶等^[14]基于MODIS-NDVI数据, 利用像元二分模型从区域和像元尺度分析了呼伦贝尔草原2000—2016年植被覆盖度时空变化和驱动机制；齐亚霄等^[15]利用天山北坡2001—2015年MODIS-NDVI数据, 结合像元二分法进行植被覆盖度提取, 并对其空间分布特征、随时间及地形的动态变化和面积加权重心进行分析。

近年来, 广东省作为沿海开放省份, 随着城镇化和工业化进程快速推进, 土地利用类型发生较大变化, 对广东省植被覆盖的空间格局也产生了深远的影响。地学信息图谱分析法通过图谱单元来记录时空复合信息^[16], 常用于土地利用变化的研究^[17-18], 然而, 将土地利用变化图谱与植被覆盖变化相结合的研究鲜有报道。“一核一带一区”区域发展新格局要求全省区域协调发展, 因此研究广东省植被覆盖变化及其与土地利用变化的相互驱动关系对于地区植被覆盖保持、土地利用结构优化、协调生态环境保护和经济发展具有重要意义。基于此, 本文利用GIS技术, 结合NDVI数据, 采用像元二分模型对2000—2020年间广东省植被覆盖时空变化特征进行分析, 并结合土地利用变化图谱特征, 探究土地利用与植被覆盖之间的相互影响驱动关系, 以期为广东省生态安全格局构建及可持续发展决策提供理论支持和科学依据。

1 研究区概况及数据处理 1.1 研究区概况

广东省地处中国大陆的南端, 位于109°45′—117°20′E和20°09′—25°31′N之间, 现辖21个地市, 地势北高南低, 地貌以山地、丘陵占优势^{[1, 19]}(图 1)。广东属于东亚季风区, 具有热带、亚热带季风和海洋等气候特点, 常年植被覆盖较为丰富, 地带性植被以北热带季雨林、南亚热带季风常绿阔叶林、中亚热带典型常绿阔叶林和沿海的红树林等为主^[20]。作为改革开放的排头兵、先行地、实验区, 广东省承担着先行先试、敢为人先, 率先建设生态文明、实现人与自然和谐共生现代化的历史使命^[21]。

1.2 数据来源及处理

本文NDVI数据来源于美国国家航天局NASA提供的MOD13Q1级产品(http://lpdaac.usgs.gov), 为更准确反映植被覆盖状况, 选取植被生长季(4—10月份)作为研究时段。时间范围为2000(4月—10月)—2020(4月—10月), 数据格式为EOS-HDF, 空间分辨率和时间分辨率分别为250 m和16 d。NDVI数据的预处理方法如下: 首先使用MRT(MODIS Reprojection Tools)工具进行格式、投影转换和波段提取, 然后利用广东省行政边界进行研究区提取, 最后采用最大值合成法^{[10, 22-24]}排除云和大气等干扰得到年NDVI影像。

此外, 文中使用的广东省土地利用类型数据来源于中国科学院地理科学与资源研究所资源环境科学与数据中心(http://www.resdc.cn), 空间分辨率1 km, 选取2000年和2020年两期土地利用数据。本文旨在揭示广东省植被覆盖规律和与土地利用之间的响应关系, 植被覆盖与土地利用数据之间分辨率的差异不会对研究结果产生影响, 因此未对两者的分辨率进行统一。

2 研究方法 2.1 植被覆盖度计算

植被覆盖度和NDVI之间存在极显著的线性相关关系, 通过建立二者之间的转换关系, 可直接提取植被覆盖度信息。像元二分模型^[12-14]的原理: 假设遥感传感器所观测到的每个像元的信息S可表达为由绿色植被部分所贡献的信息S_v与由裸土部分所贡献的信息两部分组成, 将S分解为S_v和S_s两部分:

(1)

假设有植被覆盖面积比例即为植被覆盖度C_i, 则裸土面积比例为1-C_i。如果全由植被所覆盖的纯像元所得的遥感信息为S_veg, 则S_v=C_i×S_veg, 同理得S_s=(1-C_i)×S_soil, 则:

(2)

(3)

NDVI是反映地表植被生长状态的重要指标因子, 利用植被覆盖度和NDVI之间极显著的线性相关关系, 在像元二分模型的基础上对植被覆盖度进行提取^[15]。基于上式(2)、(3)反演得到21年植被覆盖度的时间序列如下:

(4)

式中, NDVI_V为植被覆盖部分的NDVI值, NDVI_S为土壤部分的NDVI值, C_i为植被覆盖度。NDVI_V和NDVI_S的取值是像元二分模型应用的关键, 由于大气云层、地表湿度和光照等条件的影响, NDVI_S不再是一个接近于0的定值, 其变化范围通常在-0.1—0.2之间^[13]。对于纯植被像元来说, 植被类型及其构成、植被的空间分布和植被生长的季相变化都会造成NDVI_V值的时空变异。结合实际情况, 根据整幅影像上NDVI的灰度分布, 选取累计频率为0.5%的NDVI值为NDVI_S, 频率为99.5%的NDVI值为NDVI_V, 最后得到2000—2020年的植被覆盖度分布图。

根据水利部2008年颁布的《土壤侵蚀分类分级标准》中植被覆盖度分级标准和相关研究^{[3, 15]}, 将植被覆盖度划分为5个等级: ＜30%(低覆盖度)、30%—45%(中低覆盖度)、45%—60%(中等覆盖度)、60%—75%(中高覆盖度)、＞75%(高覆盖度)。

2.2 植被覆盖度年际变化趋势 2.2.1 趋势分析

趋势性分析是以单个像元的时间变化趋势来反映整体空间变化规律, 基于像元的一元线性回归分析可模拟区域时空格局变化趋势。利用ENVI的Layer Stacking工具对21年植被覆盖度波段进行合成, 然后利用趋势分析功能对21年植被覆盖度进行一元线性回归, 生成的多波段影像包含斜率、截距和R²。其计算公式为^[25]:

(5)

式中, θ_slope为单个像元植被覆盖度回归方程的斜率；C_i为第i年的植被覆盖度；n为检测时段的年数；变量i从1到n, 为研究时段年份的序号。回归方程的斜率(θ_slope)代表植被覆盖度各栅格单元的变化趋势, 当θ_slope＞0时, 表明该栅格单元所代表的植被覆盖度随时间变化而增大, 反之当θ_slope＜0时, 表明该栅格单元所代表的植被覆盖度随时间变化呈下降趋势。

2.2.2 显著性检验

采用F检验对变化趋势进行显著性检验, P＜0.01为极显著变化, 0.01＜P＜0.05为显著变化, P＞0.05为不显著变化^{[26, 27]}。统计量计算公式如下:

(6)

式中, 称为误差平方和, 称为回归平方和, y_i为第i年植被覆盖度的像元值, 为其回归值, y为研究时段植被覆盖度的平均值, n为研究年数, 即n=21。根据检验结果将变化趋势分为6个等级: 极显著退化(θ_slope＜0, P＜0.01)；显著退化(θ_slope＜0, 0.01＜P＜0.05)；不显著退化(θ_slope＜0, P＞0.05)；不显著改善(θ_slope＞0, P＞0.05)；显著改善(θ_slope＞0, 0.01＜P＜0.05)；极显著改善(θ_slope＞0, P＜0.01)。

2.3 土地利用变化图谱分析 2.3.1 建立土地利用变化图谱

目前主要是通过对不同时期的土地利用数据进行空间叠加, 以反映一定时期内的土地利用格局变化, 而图谱就是“空间与过程”研究的时空复合体。因此, 借鉴已有的图谱融合方法^[28-29], 对前后两期土地利用类型的图谱单元代码按照地图代数叠加运算而实现图谱融合, 从而记录图谱单元的演化过程。具体操作为:

(7)

式中, C为表征研究时段内土地利用方式演化特征的图谱单元栅格图；A为前一时期即2000年广东省土地利用类型的图谱单元栅格属性值；B为后一时期即2020年广东省土地利用类型的图谱单元栅格属性值。根据土地利用现状遥感监测数据, 本文基于耕地, 林地, 草地, 水域, 城乡、工矿、居民用地和未利用地6个一级土地利用类型对2000年、2010年和2020年三期图谱单元进行融合, 由此得到广东省2000—2020年的土地利用变化图谱和2000—2010—2020年土地利用变化模式图谱特征。

2.3.2 土地利用变化图谱特征分析

为了更明显表现土地利用类型图谱单元的变化特征, 引入土地利用变化的变化比率和空间分离度指标。变化比率表示某种土地利用变化类型的图谱单元面积占研究区内所有土地利用变化类型图谱单元总面积的比率；空间分离度描述图谱单元在空间分布上的分散程度, 分离度越大, 图谱单元在空间上分布越离散。计算公式如下^[30]:

(8)

(9)

式中, C_ij、F_ij分别表示变化比率和空间分离度；N_ij、A_ij分别表示初期(t)第i种土地利用类型转变为末期(t+Δt)第j种土地利用类型的图谱单元数和面积；n是土地利用类型的个数。

3 结果与分析 3.1 植被覆盖空间分异及年际变化 3.1.1 植被覆盖度的空间分布特征

利用广东省2000—2020年的植被覆盖度图, 计算得到多年平均植被覆盖度图(图 2)。可以看出, 广东省85%以上的区域处于中高植被覆盖度及以上, 中等植被覆盖度、中低植被覆盖度和低植被覆盖度分别占研究区总面积的6.27%、3.66%和4.37%。从空间分布来看, 广东省植被覆盖空间集聚性显著, 高植被覆盖区主要分布在粤北地区(梅州、清远、河源、韶关、云浮5市), 因粤北地区以丘陵、山地为主, 是原生型亚热带常绿阔叶林、天然针叶林的集中分布区, 植被覆盖度远高于全省其他地区；粤西地区(湛江、茂名、阳江3市)是全省热量资源最丰富的区域, 植被种类丰富, 植被覆盖度仅次于粤北地区；粤东地区(潮州、揭阳、汕头、汕尾4市)是广东省人口密度较大的地区, 频繁的人类活动使得区域原生植被保存受到影响, 植被覆盖度低于粤北、粤西地区^[31]；而低植被覆盖区则主要集中在珠三角地区, 这与广东省的地势特征、城市化进程以及经济发展水平的区域差异密切相关^[32]。

3.1.2 植被覆盖度的年际变化特征

对研究区2000—2020年4—10月的平均植被覆盖度进行统计分析发现(图 3), 2000—2020年广东省平均植被覆盖度在0.7和0.76之间变化, 整体以0.23%/a的变化率呈现波动增长趋势。对研究区不同等级的植被覆盖度进行统计分析(图 4), 结果发现: 2000—2020年间, 研究区的高植被覆盖度面积大致呈增长态势, 中高植被覆盖度面积呈减少趋势, 中等、中低和低植被覆盖度面积呈波动状态。其中: (1)高植被覆盖度所占比例在46%—67%, 2000年占比最低, 2020年达最高；(2)中高植被覆盖度所占比例在18%—36%, 2000年占比最高, 2020年达最低；(3)中等、中低和低植被覆盖度所占比例在14%—18%之间波动, 其中2000年为最高, 占17.83%；2016年达最低, 为13.93%。总体来看, 高植被覆盖度所占比例呈逐渐增加的趋势, 这与广东省2000—2018年NDVI的时空变化趋势一致^[1]。

由于年际变化是所有像元进行平均的结果, 为深入分析植被覆盖变化趋势的空间差异, 对2000—2020年植被覆盖度进行逐像元线性拟合。结果发现: 植被覆盖极显著退化的区域主要集中分布在珠三角、粤东沿海和粤西沿海等地, 增幅较大的区域则零散分布在各个地市。F检验表明, 绝大部分像元植被覆盖度不存在显著的增减趋势(图 5), 呈极显著和显著改善的像元占比36.92%, 呈极显著和显著退化的像元占比5.84%(表 1)。

为了进一步比较各地市植被覆盖度的年际变化差异, 对广东省、珠三角、粤东、粤西和粤北地区植被覆盖显著变化情况进行统计(表 2)。结果表明: 广东省极显著或显著改善面积大于退化面积, 整体呈改善状态；从各地区来看, 极显著和显著改善区主要集中在粤北地区, 极显著和显著退化区则主要集中在珠三角地区；从面积差来看, 粤北地区极显著变化面积差最大, 占广东省的44.9%, 显著变化面积差占比52.8%。这表明广东省近21年的植被恢复和生态建设取得了一定成果, 主要集中分布在粤北地区。

3.2 土地利用变化图谱特征 3.2.1 图谱单元的数量变化特征

2000—2020年, 广东省土地利用变化的总面积达57341 km², 存在30类土地利用变化图谱(表 3、图 6)。从变化比率看, 耕地和林地的相互转化(编码12和21)占土地利用变化总量的45.6%, 是该时段广东省最主要的土地利用变化类型。被占用与补充比率相差最大的为城乡、工矿、居民用地, 被占用比率(7.83%)小于补充比率(16.49%), 相差8.7%, 呈现以城乡、工矿、居民用地扩张为主要特征；其次为耕地: 耕地累计被占用比率为(38.84%)大于补充比率(33.33%), 相差5.5%, 这说明耕地转换过程被占用大于补充, 耕地数量净减少；林地被占用规模也较明显: 林地被占用比率(35.76%)大于补充比率(33.1%), 相差2.7%；草地、水域和未利用地被占用和补充比率差较小, 整体变化不大。广东省整体呈现耕地和林地被占用, 城乡、工矿、居民用地补充的土地利用变化特征。

从图谱单元的变化比率和空间分离度来看(图 6), 耕地与林地(编码12、21)、耕地与城乡、工矿、居民用地(编码15、51)和林地与草地(编码23、32)互换的分离度最低且相近；未利用地与林地, 草地, 城乡、工矿、居民用地, 水域之间的相互转化较少, 变化比率均小于0.1%, 而空间分离度较大, 说明相关变化图谱单元在区域上占比较小, 分布较分散。

基于表 3的分析, 用补充比率减去被占用比率得到补充和占用比率差, 对广东省21个地市耕地和林地, 草地, 水域, 城乡、工矿、居民用地的比率差逐个进行分析(图 7), 发现: (1)各地市耕地和林地面积补充和占用比率差均小于0, 城乡、工矿、居民用地比率差均大于0, 草地和水域的比率差在0上下浮动, 说明广东省整体耕地和林地整体面积减少, 城乡、工矿、居民用地增加, 草地和水域面积变化不大, 与之前结论相符；(2)除去肇庆、珠海、佛山等地, 耕地和林地与城乡、工矿、居民用地变化曲线基本以x轴为轴对称, 说明被占用的耕地和林地面积用于补充城乡、工矿、居民用地；(3)相比其他地区, 珠三角地区占用大量的耕地和林地用于发展城乡、工矿、居民用地, 以广州、东莞最为显著。总的来说, 广东省呈现占用耕地和林地转换成城乡、工矿、居民用地的显著特征, 以珠三角地区最为显著。

3.2.2 耕地和林地变化模式的图谱特征

为了进一步研究耕地和林地利用变化模式, 通过融合2000、2010和2020年3期土地利用数据, 构建耕地和林地利用变化模式图谱。由表 4可知, 全省耕地利用变化模式以后期变化和前期变化为主, 反复变化次之, 持续变化所占比例较小；林地利用变化模式则有所不同, 以后期变化为主要变化类型, 反复变化次之, 前期变化和持续变化所占比例较小。

(1) 前期变化型。前期变化型模式指土地仅在2000—2010年间出现转移的类型。在耕地的前期变化模式中, “耕地→林地→林地(编码122)”(75.1%)和“耕地→城乡、工矿、居民用地→城乡、工矿、居民用地(编码155)”(16.5%)占本阶段转移总面积的91.6%, 其中以“耕地→林地→林地(编码122)”为主要转移类型, 面积达7269.8 km²；在林地的前期变化模式中, “林地→耕地→耕地(编码211)”(70.8%)为主要转移类型, 面积约为3534.8 km², 远小于“耕地→林地→林地(编码122)”, 其次为“林地→城乡、工矿、居民用地→城乡、工矿、居民用地(编码255)”(12.8%)。

(2) 后期变化型。后期变化型模式指土地仅在2010—2020年间出现转移的类型。耕地的后期变化型模式中“耕地→耕地→林地(编码112)”(52%)和“耕地→耕地→城乡、工矿、居民用地(编码115)”(29%)占本阶段转移总面积的81%, 林地则以“林地→林地→耕地(编码221)”(63.8%)和“林地→林地→草地(编码223)”(19.7%)为主要转移类型。

(3) 反复变化型和持续变化型。反复变化型是指土地仅在2010年发生转移, 2000和2020年土地利用类型保持一致；持续变化型指土地在2000—2010年间持续发生转移。耕地和林地的反复变化模式分别以“耕地→林地→耕地(编码121)”(61.4%)和“林地→耕地→林地(编码212)”(67.3%)为主要转移类型；持续变化型模式则分别以“耕地→林地→城乡、工矿、居民用地(编码125)”(23.4%)和“林地→耕地→城乡、工矿、居民用地(编码215)”(24.1%)为主要转移类型。

结合耕地和林地的四种土地利用类型变化模式, 可见城乡、工矿、居民用地在耕地和林地的占用和补充过程中起了重要的作用, 反映出城乡建设过程中对耕地和林地的大量占用。

3.3 植被动态与土地利用变化响应

结合2000—2020年广东省植被动态变化及土地利用变化特征来看(图 5、图 8), 改善是最主要的植被动态变化形式, 占比高达78%。植被改善区以耕地和林地、林地和草地的相互转化为主, 说明土地利用类型之间的相互转化会使得植被覆盖发生变化。土地利用变化图谱中变化最为显著的是“耕地→林地(编码12)”图谱单元类型, 分布相对分散；其次图谱变化比较明显的是“耕地→城乡、工矿、居民用地(编码15)”、“林地→城乡、工矿、居民用地(编码25)”、“耕地→水域(编码14)”、“水域→耕地(编码41)”和“水域→城乡、工矿、居民用地(编码45)”, 主要集中分布在珠三角、粤东沿海和粤西沿海等地, 与植被覆盖显著和极显著退化区的分布相似(图 5), 说明向城乡、工矿、居民用地转化是造成植被覆盖退化的主要原因。

通过2000年、2010年和2020年不同土地利用类型的植被覆盖度变化分析可知(图 9)不同土地利用类型面积占比发生明显变化, 耕地面积先减后增, 整体呈减少趋势, 这与耕地的前期变化型模式保持一致, 主要表现为耕地转化为林地和城乡、工矿、居民用地；林地面积先增后减, 总面积有所减少, 这与耕地和林地前后期的互相转化有关；草地和水域面积变化不大；城乡、工矿、居民用地面积一直保持增长趋势, 总面积明显增加。各土地利用类型的植被覆盖度存在明显差异, 其中耕地、林地和草地以中高植被覆盖度及以上为主。2000—2020年间, 耕地、林地和草地高植被覆盖占比持续增加, 中高和中等植被覆盖占比持续减少, 其中林地以高覆盖为主, 高植被覆盖面积显著增加；水域和城乡用地都以低覆盖为主。结合植被覆盖度明显的上升趋势分析, 各土地利用类型对植被覆盖度的贡献各异, 以林地的贡献最大。

综合图 5、图 8和图 9可以发现: (1)城市化进程导致耕地、林地等向城乡、工矿、居民用地转化, 使得植被覆盖发生退化；(2)在生态城市指引下带来的林地面积的恢复令城市绿化度增加, 进而改善植被覆盖。

4 讨论与结论

作为衡量地区生态的重要因子, 由于不同研究采取的置信度水平存在差异可能造成结果的差异, 但基本都是利用NDVI和植被覆盖度的正相关关系进行提取^[15]。过去已有大量研究运用不同数据源从不同时间和空间尺度对植被覆盖动态变化进行了研究^{[1, 33]}, 得出我国和广东省总体植被覆盖逐渐绿化的趋势。

4.1 讨论

近年来, 广东深入落实自然资源“两统一”职责, 推动新时代广东生态保护修复取得良好成效^[21], 植被覆盖整体呈现波动增长趋势。围绕广东省“三线一单”生态环境分区管控方案和“一核一带一区”区域发展新格局, 结合本文研究结果, 对珠三角核心区、沿海经济带和北部生态发展区实施差异性生态保护修复策略。

(1) 珠三角核心区对标国际一流湾区, 强化创新驱动和绿色引领, 实施更严格的生态环境保护要求。珠三角地区植被覆盖较低, 呈现“总体改善、部分退化”的显著变化趋势。植被退化区主要表现为大量的耕地和林地向城乡、工况、居民用地转化, 说明想要推进珠三角地区生态保护修复工作, 必须重点加强城市空间与绿地系统、生态空间的统筹协调, 以“逆城镇化”的发展趋势, 推动产业高端化发展和促进生态环境质量协调化提升。生态修复重点区域主要集中在珠三角核心区, 实施生态建设和生态修复, 在保护生态空间的基础上合理开展生产建设活动^[36]。

(2) 沿海经济带主要包含珠三角沿海7市和粤东、粤西地区。粤东、粤西地区植被覆盖较高, 沿海少部分区域呈现植被退化趋势。沿海经济带东西两翼地区生态环境保护与经济发展介于珠三角核心区和北部生态发展区之间, 发展不平衡不充分的问题较为突出, 远未达到珠三角沿海地区发展水平。对于沿海经济带而言, 坚持新型城镇化的建设道路, 在推进环境治理项目、应用绿色科技的基础上, 充分利用珠三角中心的辐射作用, 发挥地区特色和资源禀赋, 最大化减少产业同质化, 打造生态环境与经济社会协调发展区, 实现生态价值提升和产业高质量发展的匹配。

(3) 北部生态发展区即粤北地区, 其大部分区域为优先保护单元, 以维护生态系统功能为主, 禁止或限制大规模、高强度的工业和城镇建设。北部生态发展区, 作为广东省的重点生态屏障, 以中高和高植被覆盖为主, 整体呈现为植被改善。随着生态环境分区管控方案的落实, 北部生态发展区将保持植被改善的发展态势, 打造广东生态安全屏障。然而, 如何在高生态保护水平的前提下实现经济高质量发展, 是北部生态发展区面临的重大课题。立足北部生态发展区的生态功能定位, 坚持生态优先、绿色发展, 推进产业生态化和生态产业化的同时, 健全生态补偿体制机制, 实现乡村振兴。

习近平总书记的重要论述表明, 新时代我们要协调好新型城镇化、逆城镇化和乡村振兴的关系^[37], 以差异性生态保护修复策略和产业化、高端化和绿色化的发展思路, 实现广东省生态环境和经济水平的协同发展。

4.2 结论

本文基于2000—2020年MODIS-NDVI数据分析了广东省近21年植被覆盖变化特征, 并结合土地利用变化图谱, 探究土地利用与植被覆盖度之间的相互影响驱动关系, 得出以下结论:

(1) 广东省多年平均植被覆盖度介于0.7—0.76, 85%以上的区域处于中高植被覆盖度及以上, 植被空间集聚性特征显著, 高植被覆盖区主要分布在粤北地区, 低植被覆盖区则主要集中在珠三角区域, 这与广东省的地势特征、土地利用类型、城市化进程以及经济发展水平的区域差异密切相关。广东省植被覆盖度整体以0.23%/a的变化率呈现波动增长趋势, 高植被覆盖区面积逐年增加, 全省植被呈正增长的区域面积占比达78.18%；近5.84%的区域发生显著或极显著退化, 主要集中于佛山、中山、广州、东莞等珠三角核心城市。随着广东省沿海地区低端产业向中部、北部转移以及经济不发达地区城市化、工业化进程的加速, 导致广东省中等植被覆盖度区域向中部的肇庆、清远、惠州、云浮扩散, 较高、高植被覆盖区域逐渐向北部和东部方向转移, 其中高植被覆盖区域退缩到经济较为落后但森林资源丰富的粤北部, 如韶关、河源、梅州一带^[33]。

(2) 广东省的土地利用类型变化主要呈现以下3个特征: ①耕地和林地的相互转化是2000—2020年广东省主要的土地利用变化类型；②珠三角、粤东和粤西沿海地区的土地利用变化类型以占用耕地和林地向城乡、工矿、居民用地转化为主；③城乡、工矿、居民用地在耕地和林地的占用和补充过程中起了重要的作用, 反映出城乡建设过程中对耕地和林地的大量占用。

(3) 国内外相关研究表明, 人为活动主导下的土地利用变化与自然环境因子的时空分异构成了区域植被覆盖变化的重要驱动力来源^[34]。土地利用变化对于植被覆盖变化的影响主要源于土地转移过程中短时限的人为干扰和转移完成后较长时限的用地功能变化, 往往表现为群落物种组成、生物生产能力、生态系统类型等方面的综合变化^[35], 并在较大尺度上表现出明显的时空差异特征。研究结果也表明, 研究区植被覆盖度的变化与土地利用变化存在响应联系, 土地利用变化是研究区植被覆盖动态变化的重要影响因素之一。各土地利用类型对植被覆盖度的贡献各异, 以林地的贡献最大；向城乡、工矿、居民用地转化的图谱与植被显著和极显著退化区高度重合, 说明耕地和林地向城乡、工矿、居民用地转化是造成植被覆盖退化的主要原因。