2000—2020年中国植被时空变化格局及其对气候作用的滞后响应分析

doi:10.6046/zrzyyg.2024381

2000—2020年中国植被时空变化格局及其对气候作用的滞后响应分析

许颖^,¹, 张润泽², 郭兵^,²

1.河南财经政法大学城乡规划学院，郑州 454000

2.山东理工大学建筑工程与空间信息学院，淄博 255000

Spatiotemporal variations of 2000—2020 vegetation in China and their time-lag responses to climate

XU Ying^,¹, ZHANG Runze², GUO Bing^,²

1. Henan University of Economics and Law，School of Urban and Rural Planning，Zhengzhou 454000，China

2. School of Civil Engineering and Geomatics，Shandong University of Technology，Zibo 255000，China

通讯作者: 郭兵（1987-），男，博士，教授，主要研究方向为全球变化与土地退化。Email：guobingjl@163.com。

责任编辑: 李瑜

收稿日期: 2024-11-27 修回日期: 2025-02-12

基金资助:

河南省科技攻关项目“融合航空影像和激光雷达点云豫西黄河流域切沟监测及体积预测模型研究”(242102321113)
河南省高等学校重点科研项目“融合机载Lidar点云和航空影像的豫西黄河流域切沟监测及体积估算模型优化”(24B420001)
国家自然科学基金项目“水-土-植耦合新视角下黄河三角洲盐渍化演变机理的跨尺度不变性与差异性研究”(42471329)

Received: 2024-11-27 Revised: 2025-02-12

作者简介 About authors

许颖（1986-），女，博士，讲师，主要研究方向为遥感数据处理与应用。Email：y_xu_lidar@huel.edu.cn。

摘要

全球气候变化背景下，中国植被生态系统发生了深刻变化。然而，中国不同生态子区、不同历史时期植被生态演变机制如何，有何差异，亟须深入探索。该研究基于归一化植被指数（normalized differential vegetation index，NDVI）数据，采用重心模型、时滞分析、地理探测器和偏相关分析等方法对中国6大生态子区进行时空演变格局及驱动机制的分析。结果表明：①2000—2020年间，中国大陆地区植被覆盖度呈现出自东向西逐渐递减的趋势；②6大生态子区植被覆盖度均呈现增加趋势，其中中南地区增幅（变化斜率）最高，为0.003 9，华东地区增幅（变化斜率）最低，为0.002；③2000—2010年间，植被覆盖等级升高的面积占比为92%，2010—2020年间植被覆盖等级升高的面积占比为71%；④在全国范围内，栽培植被、灌丛与降水滞后性普遍为前1—3月，栽培植被、针叶林与气温滞后性普遍为当前月，阔叶林与气温滞后性普遍为前1—2月，不同地区不同植被类型滞后时间存在异质性；⑤东北、华北、西北、西南地区植被演变主导影响因子均为降水，华东地区植被演变主导影响因子为土地利用和国内生产总值（gross domestic product，GDP），而中南地区植被演变主导影响因子为降水、土地利用。研究成果可为不同生态区植被的修复和保护提供重要的数据支撑。

关键词： NDVI; 植被覆盖度; 时滞效应; 重心; 地理探测器; 中国

Abstract

Amid global changes，China’s vegetation ecosystem has undergone profound transformations. However，there is an urgent need to thoroughly explore the mechanisms underlying the ecological evolution of vegetation in different ecological subregions and historical periods，as well as their differences. Therefore，based on normalized difference vegetation index （NDVI） data，this study investigated the spatiotemporal evolution of vegetation across six major ecological subregions in China and its driving mechanisms using methods such as the gravity center model，lag analysis，geographical detectors，and partial correlation analysis. The results indicate that from 2000 to 2020，mainland China witnessed a decreasing trend in vegetation coverage from east to west. Vegetation coverage increased in all six ecological subregions，with the highest increase （slope of change） observed in the south-central part of China （0.003 9） and the lowest in eastern China （0.002）. From 2000 to 2010，regions with increased vegetation coverage accounted for 92%，and this proportion dropped to 71% from 2010 to 2020. Heterogeneous lag times were observed across different vegetation types in varying regions. Specifically，cultivated vegetation and shrubland generally exhibited a 1 to 3-month lag in response to precipitation；cultivated vegetation and coniferous forests presented a lag limited to the current month in relation to temperature，and broadleaf forests generally displayed a 1 to 2-month lag in response to temperature. Precipitation is identified as the dominant factor driving vegetation changes in North China and the northeastern，northwestern，and southwestern parts of China. In eastern China，land use and gross domestic product （GDP） represent the primary driving force behind vegetation change. In the south-central part of China，both precipitation and land use serve as dominant factors. The results of this study can provide significant data support for vegetation restoration and protection in different ecological regions.

Keywords： normalized difference vegetation index （NDVI）; vegetation coverage; time-lag effect; gravity center; geographical detector; China

PDF (10872KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

许颖, 张润泽, 郭兵. 2000—2020年中国植被时空变化格局及其对气候作用的滞后响应分析[J]. 自然资源遥感, 2025, 37(5): 254-266 doi:10.6046/zrzyyg.2024381

XU Ying, ZHANG Runze, GUO Bing. Spatiotemporal variations of 2000—2020 vegetation in China and their time-lag responses to climate[J]. Remote Sensing for Land & Resources, 2025, 37(5): 254-266 doi:10.6046/zrzyyg.2024381

0 引言

植被作为地表生境状况的表征，对于流域生态系统的生态平衡、水文、地貌过程具有重要的调节作用^[1]，植被对环境变化具有极高的敏感性，被视为全球变化的“指示器”^[2]。植被也是连接大气圈、水圈和土壤的天然纽带，在水土保持、维护气候和生态稳定性方面起到至关重要的作用^[3]。陆地植被的生长不仅受到气候变化的影响，也受到人类活动的影响。近年来，伴随着全球变暖、城市扩张和人口增长等情况的凸显，对植被也造成一定程度的影响^[4]。深入分析植被时空变化特征及其影响因子，对于实现“双碳”目标，维护生态平衡具有重要意义。

近年来，众多学者基于遥感影像数据对不同地区植被覆盖情况及驱动机制进行了相关分析。鲁佳琪等^[5]基于降水蒸散发指数（standardize precipitation evaporation index，SPEI）数据集、陆地储水量数据对黄河流域内蒙古段气象干旱与水文干旱对植被基于归一化植被指数（normalized differential vegetation index，NDVI）变化的影响进行研究，发现气象干旱和水文干旱分别集中在1月、10月和1月，大部分地区的植被短期内受水文干旱的影响较大；袁江龙等^[6]基于NDVI数据探究了2000—2020年昆仑山植被覆盖度的时空变化规律，并采用地理探测器对NDVI变化的主导因子进行探究，发现该地区NDVI在时间上呈波动增加趋势，空间上呈现“西高东地、北高南低”的分布格局；胡仁杰等^[7]基于比值植被指数（ratio vegetation index，RVI）探讨了NDVI的饱和性，并评估了NDVI饱和性对青藏高原地区高寒草甸地上生物量时空动态变化分析的影响，发现全年6月—9月初（全年第161—257天）饱和性呈现先自东南向西北延伸，后自西北向东南消退的变化趋势；申子傲等^[8]基于NDVI，利用slope趋势分析和变异系数对2000—2020年河西内陆河流域植被覆盖变化进行分析，并采用地理探测器对影响植被NDVI变化的自然因子和人为因子进行分析，发现该地区植被覆盖总体呈现波动上升的趋势，自然因素对植被影响占据主导地位；Mason等^[9]基于 MODIS NDVI 数据对澳大利亚植被动态的研究表明，研究区植被出现先减少后增加的趋势；Muradyan等^[10]对亚美尼亚山区的植被变化特征和生长滞后性进行了研究，发现 NDVI 最高的时段为植被生长季，且植被变化与气候的变化响应存在滞后性。然而，在全球变化背景下，中国植被生态系统发生了深刻变化，并且不同生态子区的植被变化影响因子存在显著差异，亟需深入探索^[11]。同时，不同植被类型、不同子区相同植被NDVI对气温与降水时滞效应存在显著差异，目前尚不清晰，仍需深入探索^[12]。

本研究基于MODIS NDVI时序数据集，采用重心模型、时滞效应分析法、偏相关分析及地理探测器，在分区视角下探讨了中国不同生态子区植被NDVI时空演变格局，并对影响植被NDVI变化的时滞时间、主导驱动因子进行分析，为区域植被生态系统的保护和修复提供重要的决策支持。

1 研究区概况及数据源

1.1 研究区概况

如图1所示，中国位于亚洲东部，太平洋西岸，国土面积约960万km² （图1），海岸线达1.8万km。中国地势西高东低，呈现3级阶梯状分布，高原、山岭、平原、丘陵、盆地5种基本地形在中国均有分布。由于国土面积广阔、经纬度跨度大（南北跨越纬度50°，东西跨越经度60°），使气温、降水组合形式多样，并构成复杂多样的气候环境，其中季风气候最为显著，这也使得中国经常发生气象灾害和地质灾害。由于南北纬度跨越大，使得气温南北空间异质性明显，在秦岭-淮河线以南，主要是亚热带与热带地区，大兴安岭北端及两侧分布于寒温带，中温带介于暖温带与寒带之间，面积广阔。受东南季风的影响，降水分布呈现出自东南向西北递减的趋势。东南沿海地区会最先得到季风带来的水汽，形成丰富的降水，而西南地区深居内陆，使得水汽较少，降水量少。植被类型丰富，主要有草原、热带雨林、常绿阔叶林、落叶阔叶林以及针叶林，其中阔叶林广泛分布于温带、暖温带地区，针叶林主要分布于寒温带等地区。土壤类型复杂多样，分布随地理位置、地形高度变化而呈有规律更替的现象，既表现为水平地带性，又与垂直地带性、地域性等有关。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 研究区概况（审图号：GS京（2025）1936号）

Fig.1 Overview of the study area

1.2 数据源及预处理

本研究所使用的植被NDVI数据、植被类型数据均来源于国家青藏高原数据中心（https://data.tpdc.ac.cn/zh-hans/data/），其中NDVI数据集原始格式均为.tif，空间分辨率为250 m，时间分辨率为月尺度，所需时间跨度为2000—2020年，植被类型数据为矢量格式，利用ArcGIS 10.7的面转栅格工具将其转为.tif格式。气候站点数据来源于中国气象数据网（http://data.cma.cn/），年降水、年均温、年积温以及年日照数据等均由该数据网获取，利用ArcGIS 10.7的Geostatistical Analyst工具进行插值获取1 km格网数据。土地利用数据、社会经济数据（人口密度、国内生产总值（gross domestic product，GDP））以及6大分区数据来源于中国科学院资源环境科学与数据中心（http://www.resdc.cn/DOI/），时间分辨率为1 a，其中土地利用数据空间分辨率为30 m，社会经济数据空间分辨率为1 km，高程、坡度数据通过SRTM 90 m中提取得到，数据来源于地理空间数据云（http://www.Gscloud.cn）。上述数据集均采用ArcGIS10.7进行数据异常值的剔除、裁剪、重采样（1 km）及重投影（正弦投影（sinusoidal projection，SIN）→克拉索夫斯基投影）。

1.3 研究方法

1.3.1 相关性分析

相关系数能够表征不同要素之间的相关程度^[13]，利用基于像元的皮尔逊相关系数分析和探讨植被NDVI与典型气候因子的相关性及其显著水平。其皮尔逊相关系数、偏相关系数计算公式如下：

（1）r_x_，_y=$\frac{\sum _{i=1}^{n}（{x}_{i}-\stackrel{-}{x}）（{y}_{i}-\stackrel{-}{y}）}{\sqrt{\sum _{i=1}^{n}（{x}_{i}-\stackrel{-}{x}{）}^{2}\sum _{i=1}^{n}（{y}_{i}-\stackrel{-}{y}{）}^{2}}}$，

（2）${r}_{x{y}_{1}，{y}_{2}}$=$\frac{{r}_{x{y}_{1}}-{r}_{x{y}_{2}}{r}_{{y}_{1}{y}_{2}}}{\sqrt{1-{r}_{x{y}_{2}}^{2}}\sqrt{1-{r}_{{y}_{1}{y}_{2}}^{2}}}$，

式中：r_x_，_y为x变量与y变量的相关系数（${r}_{x{y}_{1}}$，${r}_{x{y}_{2}}$，${r}_{{y}_{1}{y}_{2}}$与之相类似），其中正值表示正相关，否则反之；x_i与y_i分别为2个变量在第i年的数值；${r}_{x{y}_{1}，{y}_{2}}$为控制变量y₂的情况下x与y₁的偏相关系数。

当控制多个变量时，有高阶偏相关系数，计算公式如下：

（3）${r}_{xy，{y}_{1}{y}_{2}\dots {y}_{g}}$=$\frac{{r}_{xy，{y}_{1}{y}_{2}\dots {y}_{g-1}}-{r}_{x{y}_{g}，{y}_{1}{y}_{2}\dots {y}_{g-1}}{r}_{y{y}_{g}，{y}_{1}{y}_{2}\dots {y}_{g-1}}}{\sqrt{1-{r}_{x{y}_{g}，{y}_{1}{y}_{2}\dots {y}_{g-1}}^{2}}\sqrt{1-{r}_{y{y}_{g}，{y}_{1}{y}_{2}\dots {y}_{g-1}}^{2}}}$。

偏相关系数采用t显著性检验（Student’s t-test），显著性检验公式如下：

（4）t=$\frac{r}{\sqrt{1-{r}^{2}}}\sqrt{n-m-1}$，

式中：r为偏相关系数；n为样本数目；m为变量数目。

本研究采用二阶偏相关系数，通过查表获取在相应自由度n-m-1下p=0.05以及p=0.01时的t值，并与已计算出的t值进行比较，以进行p<0.05（显著相关）以及p<0.01（极显著相关）置信水平下对植被NDVI与不同因子之间的线性关系性显著性检验。

1.3.2 趋势分析

使用一元线性回归分析中国2000—2020年NDVI的变化趋势^[14]。计算公式为：

（5）θ_slope=$\frac{n\times \sum _{i=1}^{n}i\times NDV{I}_{i}-\sum _{i=1}^{n}i\sum _{i=1}^{n}NDV{I}_{i}}{n\times \sum _{i=1}^{n}{i}^{2}-{\left(\sum _{i=1}^{n}i\right)}^{2}}$，

式中：i为年份；n为时序年数；NDVI_i为第i年NDVI；θ_slope为趋势线斜率，其值为正或为负时分别表示植被NDVI随时间变化呈改善或退化趋势。若θ_slope>0，代表植被NDVI有所增加；若θ_slope<0，则表示植被NDVI呈减少趋势。为检验样本的总体均数是否存在显著差异，变化趋势的显著性检验采用F检验（Fisher’s test）。计算公式为：

（6）F=U×$\frac{n-2}{Q}$，

（7）U=$\sum _{i=1}^{n}（{\hat{y}}_{i}-\overline{y}{）}^{2}$，

（8）Q=$\sum _{i=1}^{n}（{y}_{i}-\hat{{y}_{i}}{）}^{2}$，

式中：U为回归平方和；Q为残差平方和；y_i为第i年的NDVI；${\hat{y}}_{i}$为NDVI回归值；$\overline{y}$为NDVI均值；n为年数。

如表1所示，依据slope及显著性水平（p<0.05，显著；p<0.01，极显著；相应F值由查找表获取），将植被NDVI增长情况划分为5级。

表1 显著性分级

Tab.1 Significance classification

变化类型	F检验显著性水平	NDVI变化斜率
明显增加	p≤0.01	k≥0.001
略微增加	0.01<p≤0.05	k≥0.001
基本不变	p>0.05	-0.001<k<0.001
略微减少	0.01<p≤0.05	k≤-0.001
明显减少	p≤0.01	k≤-0.001

新窗口打开| 下载CSV

1.3.3 地理探测器

地理探测器模型能够揭示某因子对植被NDVI的解释力大小^[15]。其相关程度用q值度量，q的取值范围为[0,1]，其值越大，表明该因子对NDVI影响力越大，反之则越小。表达式为：

（9）q=1-$\frac{\sum _{h=1}^{L}{N}_{h}{{\sigma }_{h}}^{2}}{N{\sigma }^{2}}$=1-$\frac{SSW}{SST}$，

（10）SSW=$\sum _{h=1}^{L}$N_hσ_h² ，

（11）SSW=SST=Nσ² ，

式中：h=1，…，L；L为变量Y或因子X的分类；N_h和N分别为层h和全区内单元数；σ_h²为类h的Y值的方差；σ²为全区的Y值的方差；SSW为层内方差之和；SST为全区总方差。

2 研究结果

2.1 2000—2020年间中国陆地植被空间分布格局分析

采用年最大值合成法合成逐年NDVI，对逐年NDVI 2000—2020年间共21 a数据取均值（图2），依据相关学者划分方法^[16-18]，将其划分为5等级，分别为[0,0.2）（低植被覆盖）,[0.2,0.4）（中低植被覆盖）,[0.4,0.6）（中植被覆盖）,[0.6,0.8）（中高植被覆盖）,[0.8,1]（高植被覆盖）。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 不同等级NDVI空间分布（审图号：GS京（2025）1936号）

Fig.2 Spatial distribution of different NDVI levels

如图2所示，中高植被覆盖占比面积最大，为28.3%，主要分布在中国东北、东南沿海等地区，其中东北地区境内有大兴安岭，广袤林区分布广袤，黑土肥沃，能为植物提供充足养分，植被覆盖率较高。中国低植被覆盖区面积占比次之，为26.2%，主要集中在塔克拉玛干沙漠、青藏高原西部等地区。青藏高原地势较高，高山高原气候寒冷干燥，沙漠景观干旱少雨，均不适于植被生长，致使植被覆盖率最低。高植被覆盖区面积占比为25%，主要分布在中国南方地区，该地区气候为季风气候，夏季湿热多雨，植被长势较好，适于常绿树种的生长。整体分析，自东向西，植被覆盖率逐渐降低，西方深居内陆，气候干燥，不适于植被生长，自北向南，呈现出由高植被覆盖变为高植被覆盖与中高植被覆盖混合分布的特点，南方如四川盆地等地区，四面环山，地势较周围低，导致植被覆盖较周围低，为中高植被覆盖。

2.2 2000—2020年间中国陆地植被不同时间尺度下时序变化分析

为分析植被覆盖时序变化特征，将整个中国分为6大分区，分别为东北地区、华北地区、华东地区、西北地区、西南地区和中南地区，并分别统计各大分区在2000—2020 年间逐年植被NDVI 均值，并采用折线图进行分析。如图3所示。在2000—2020年间，华东地区植被NDVI拟合值增幅为0.002，在6大地区中增幅最低，整体呈现出波动增加的趋势，拟合优度R²为0.617 1，增势较为显著。年均最大值出现在2017年，为0.784，年均最低值出现在2002年，为0.725。在2015—2016年间增幅最高，为0.015 41，在2010—2011年间降幅最高，为0.015。中南地区植被NDVI拟合值增幅为0.003 9，拟合值增幅最高，整体呈现出波动增加的趋势，拟合优度R²为0.848 4，增势较为显著。最大值出现在2017年，为0.817，最低值出现在2001年，为0.724。在2002—2003年期间增幅最高，为0.028 1；在2010—2011年间，降幅最高，为-0.017 8。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 不同区域NDVI年际变化

Fig.3 Interannual variation of NDVI in different regions

综上分析，不同区域植被覆盖等级可分为3大部分，其中东北、华东、中南地区植被覆盖最高，华北、西南地区次之，西北地区植被覆盖度最低。其中东北地区有优质黑色土壤，华东、中南地区气候适宜，适于植被生长，西北地区深居内陆，气候大陆性显著，降水稀少，有广阔的沙漠，不适于植被生长。

2.3 不同等级面积变化

为分析不同时段植被覆盖等级变化的面积变化情况，将研究时段划分为2大历史时期，分别为2000—2010年与2010—2020年，并采用转移矩阵对不同时段内面积变化进行统计，如图4所示。图4（a）中，在2000—2010年间，植被覆盖等级发生变化的区域中，由中高植被覆盖转为高植被覆盖面积占比最高，为54%，主要分布在南方地区，如安徽省、湖南省、贵州省等地区，在中国东北部内蒙古高原与大兴安岭、东北平原交界处也有分布；其次为由中植被覆盖转为中高植被覆盖，面积占比约为14%，主要集中在二三级阶梯交界处。整体分析可知，植被覆盖等级升高的面积为92%，植被覆盖等级降低的面积占比为8%，表明在2000—2010年间植被覆盖等级呈现出升高的趋势。图4（b）中，在2010—2020年间，植被覆盖等级发生变化的区域中，由中高植被覆盖转为高植被覆盖面积占比最高，为43%，主要分布在中国南方地区，如贵州省、云南省等地区，在中国东北部地区也有少量分布；其次为中植被覆盖转为中高植被覆盖，面积占比约为11%，主要分布在内蒙古高原与黄土高原交界处。整体分析可知，植被覆盖等级升高的面积占比为71%，植被覆盖等级较低的面积占比为29%，表明在2011—2020年间植被覆盖等级呈现出升高的趋势。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 不同时期不同等级NDVI转移变化（审图号：GS京（2025）1936号）

Fig.4 Changes in NDVI levels across different periods

综上所述，在2000—2020年间，中国大陆地区植被覆盖整体呈现由低等级向高等级转移的趋势。

2.4 滞后分析

由于植被的生长受到气温与降水的影响，而气温与降水反映到植被的生长需要一定的时间，即植被对气温与降水的表现具有一定的滞后性。同时不同地理位置、不同植被类型对滞后性均有一定影响。如图5所示，通过植被NDVI与不同时期的气温与降水的相关系数以表征滞后性的大小。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 不同植被类型滞后性空间分布（审图号：GS京（2025）1936号）

Fig.5 Spatial distribution of lag effects for different vegetation types

如图5（a）所示，植被NDVI与降水滞后性普遍为前1—2月与前1—3月。其中前1—2月主要分布在二三级阶梯交界处，华东平原、青藏高原地区，前1—3月主要分布在东北平原、长江中下游平原等地区。在武夷山脉等地区有少量滞后前二月与前三月分布。如图5（b）所示，植被NDVI与气温滞后性普遍为当前月、前1—3月、前1—2月。其中当前月主要分布在华东平原地区，前1—3月主要分布在二三级阶梯交界处、中国南方沿海一带等地区，前1—2月主要分布在辽东丘陵、青藏高原、东南丘陵等地区。

不同植被类型分析，栽培植被、灌丛与降水滞后性普遍为前1—3月，草甸、草丛与降水滞后性普遍为前1—2月，栽培植被、针叶林与气温滞后性普遍为当前月，阔叶林与气温滞后性普遍为前1—2月。不同地区不同植被类型滞后时间存在异质性。

对不同地区不同植被类型滞后性进一步分析，结果见图6。如图6（a）所示，东北地区针叶林植被受前1—2月的降水影响最大，相关系数为0.82，栽培植被等植被类型受前1—3月降水影响最大。如图6（b）所示，东北地区不同植被类型均受前1—2月气温影响最大，当前月气温对针阔林影响次之，相关系数为0.92。如图6（c）所示，华东地区栽培植被、针叶林和草丛受当前月降水影响最大，相关系数分别为0.73，0.72和0.82，阔叶林受前1—3月降水影响最大，相关系数为0.60，其次为前一月，相关系数为0.59。如图6（d）所示，华东地区灌丛与阔叶林受前二月气温影响最大，相关系数分别为0.42和0.34，其他植被类型均受前1—2月气温影响最大。如图6（e）所示，西北地区不同植被类型均受前1—2月降水影响最大，与草原、阔叶林和灌丛的相关系数分别为0.78，0.74和0.76，前1—3月降水对所有植被类型的影响次之。如图6（f）所示，西北地区不同植被类型均受前1—2月气温影响最大，其中与沼泽、阔叶林和针叶林的相关系数分别为0.96，0.93和0.89。如图6（g）所示，华北地区阔叶林、草甸和草原均受前一月降水影响最大，相关系数分别为0.85，0.86和0.74，沼泽、灌丛和草丛受前1—2月降水影响最大，相关系数分别为0.92，0.86和0.83。如图6（h）所示，华北地区针叶林、阔叶林和草原均受前1—3月气温影响最大，相关系数分别为0.82，0.77和0.81。如图6（i）所示，西南地区草甸和高山植被受前1—2月降水影响最大，相关系数分别为0.83和0.91。如图6（j）所示，西南地区灌丛受前1—3月气温影响最大，相关系数为0.67，针叶林和草丛受前二月气温影响最大，相关系数分别为0.47和0.55。如图6（k）所示，中南地区不同植被类型均受前1—3月降水影响最大，其中与阔叶林、针叶林和栽培植被的相关系数分别为0.56，0.54和0.53。如图6（l）所示，中南地区沼泽和草丛受前1—3月气温影响最大，相关系数分别为0.79和0.76，其他植被类型受前1—2月气温影响最大，其中与阔叶林、针叶林的相关系数分别为0.83和0.80。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 不同分区滞后性

Fig.6 Lag effects in different regions

2.5 不同历史时期、不同子区交互主导因子分析

不同因子之间的交互作用对植被NDVI的影响具有一定的影响，且不同地区交互作用不同，采用地理探测器对影响植被NDVI变化的交互因子进一步分析。

如图7所示，东北地区和华北地区大部分因子之间的交互作用为非线性增强，其余地区大部分因子之间的交互作用为双因子增强。在东北地区、华东地区、西北地区和西南地区土地利用在多个因子的组合中均表现出较高的解释力，分别为土地利用∩降水（2010年，q=0.87）、降水∩土地利用（2020年，q=0.640）、降水∩土地利用（2000年，q=0.814）、降水∩土地利用（2010年，q=0.854），表明土地利用对相应地区植被NDVI的变化影响较大，在大陆东部沿海等地区，随着人类活动水平的提高，城镇化水平逐渐提高，城市扩张程度加剧，西北内陆等地区沙漠广阔，土地利用程度低，降水稀少，均在一定程度上影响植被NDVI的变化。在华北地区，降水对植被NDVI的影响贯穿始终，西南地区，降水在2010年和2020年对植被NDVI的影响较强，该地区位于第三级阶梯，海拔高，降水少，对植被的生长有一定的影响。在中南地区，交互主导因子由气候因子主导逐渐变为人类活动主导（2000年、2010年：积温∩日照，q=0.870，0.876；2020年：GDP∩土地利用，q=0.428），表明人类活动对中南地区植被生长的影响逐渐增强。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 交互因子

Fig.7 Interaction factors

3 讨论

3.1 2000—2020年间中国陆地植被时空变化的原因

中国陆地植被覆盖度自东向西逐渐降低。高植被覆盖区主要分布在第三级阶梯，在中国南方地区，气候以季风气候为主，夏季高温多雨，气候湿润，且南方海拔地势低，来自太平洋的暖湿气流能够深入中国南方地区，充分满足植物的生长需要，促使植被覆盖度较高^[19]。在中国东北地区也有高植被覆盖的分布，主要位于东北平原等地区，该地区冬季寒冷漫长，有利于有机质积累，夏季温暖多雨，有利于有机质的分解，从而形成肥沃的黑土，为植被生长提供了充足的养分条件，适合种植大量农作物^[20]。此外，东北地区三面环山，开口向南，有利于海洋水汽进入，降水较多，同时该地区纬度高，蒸发量小，有冻土层存在，地表水较丰富，气候较为湿润，形成冷湿的环境，易于大量森林的形成^[21]。低植被覆盖区主要分布在中国第一级阶梯，该地区有广阔的塔克拉玛干沙漠，地势较高的青藏高原，沙漠地区昼夜温差大，降水稀少，砂质土壤，普通作物难以固定根部，青藏高原海拔高，气温低，来自太平洋以及印度洋的水汽难以抵达，均不适于植被生长^[22]。

华北、华东、西北地区中，NDVI变化趋势的震荡也较为明显，且6大区域在2008—2009年间NDVI均呈现下降趋势，依据薛港等^[23]和宁晓春等^[24]的研究，华北、华东、西北地区降水、气温波动变化大，降水、气温的波动变化影响植被的生长，进而导致植被NDVI的变化。2008—2009年间由于发生极端气候（极端高温、极端降水、极端强对流等天气）的原因^[25]，影响植被的生长，进而导致2008—2009均出现植被NDVI下降的现象。

2000—2020年间，中国不同地区植被覆盖度均呈现出不同程度增加的趋势，这可能与国家的宏观调控有关。早在2010年，国家发布关于全国林地保护利用规划纲要的批复，要求做好林地保护^[26]。其中中南地区增幅最高，为0.003 9，这可能与当地积极响应国家政策，退耕还林还草有关。

3.2 中国植被变化主导影响因素及其变化原因

在2000年影响东北地区植被变化主导因子为降水，东北地区东南部地形以山脉为主，地势高而气温低，降水主要以雪的形式，不利于植被生长^[27]。而在2010年和2020年，影响因子均为土地利用，这可能与随着生活水平的提高，城镇化水平逐渐升高有关。城镇化对与植被的生长具有两面性，合理的城镇化可以提高植被覆盖程度，如在城市中实行合理的保护林地等政策^[28]；不合理的城镇化如大量开垦荒地，提高城镇化水平，可能会对原有林地造成一定程度的破坏，降低植被覆盖^[29]。在华北地区，影响植被变化的主导因子均为降水，该地区位于半湿润半干旱地区，植被生长对气候变化极为敏感，该地区生长季降水量相对稀少，水分对该地区的植被生长非常重要，使得降水对该地区植被生长起主导作用^[30]。在2000年、2010年和2020年，影响华东地区植被变化的主导因子分别为土地利用、GDP和土地利用，均与人类活动有关。华东地区作为我国经济最发达的地区之一，人类活动对植被变化具有重要影响，人类对植被变化的影响是双面的，合理的退耕还林有利于植被的恢复，增加植被覆盖，而人类的乱砍乱伐，不合理的开垦荒地，提高城镇化水平，均会降低植被覆盖^[31]。西北、西南地区植被变化的主导因子均为降水，其中，西北地区水资源短缺、生态环境脆弱，干旱缺水仍是制约该地区经济和社会发展的重要问题之一。水资源作为植被生长的重要条件之一，少量的水资源会大幅度限制植被的生长^[32]。青藏高原地区地势海拔高，降水较少，降水的增加会对植被的生长起促进作用；而在云南等地区，当处于植被生长季时，降水量充沛，过高的降水反而会抑制植被的生长^[33]。在2000年、2010年和2020年，影响中南地区植被变化的主导因子分别为降水、降水和土地利用，由气候因子逐渐变为人类活动因子，这可能与人类生活水平的提高，城镇化水平不断提高，对植被生长变化的影响逐渐变强有关^[34]。

4 结论

本研究基于MODIS NDVI时序数据集，采用重心模型、时滞效应分析法、偏相关分析及地理探测器，在分区视角下探讨了中国不同生态子区植被NDVI时空演变格局，并对影响植被NDVI变化的时滞时间、主导驱动因子进行分析。主要结论如下：

1）中国大陆地区植被覆盖度呈现出自东向西逐渐递减的趋势。

2）6大分区2000—2020年间植被覆盖度均呈现出增加的趋势，其中中南地区增幅最高，为0.003 9，华东地区增幅最低，为0.002。

3）在2000—2010年间，植被覆盖等级升高的面积占比为92%，2010—2020年间植被覆盖等级升高的面积占比为71%，表明2000—2020年间中国大陆地区植被覆盖度逐渐增加。

4）在全国范围内，栽培植被、灌丛与降水滞后性普遍为前1—3月，栽培植被与气温滞后性普遍为当前月，阔叶林与气温滞后性普遍为前1—2月。不同地区不同植被类型滞后时间存在异质性。

5）对东北、华北、西北和西南地区植被覆盖影响最强的主导因子均为降水，对华东地区植被覆盖影响最强的因子为土地利用和GDP，对中南地区植被覆盖影响最强的因子为降水和土地利用。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

颜明, 贺莉, 王随继, 等.

基于NDVI的1982—2012年黄河流域多时间尺度植被覆盖变化

[J]. 中国水土保持科学, 2018, 16(3):86-94.

[本文引用: 2]

Yan

, He

, Wang

S J

, et al.

Changing trends of NDVI in the Yellow River Basin from 1982 to 2012 at different temporal scales

[J]. Science of Soil and Water Conservation, 2018, 16(3):86-94.

[本文引用: 2]

[2]

王伟, 阿里木·赛买提, 吉力力·阿不都外力.

基于地理探测器模型的中亚NDVI时空变化特征及其驱动因子分析

[J]. 国土资源遥感, 2019, 31(4):32-40.doi:10.6046/gtzyyg.2019.04.05.