高萌萌 ，刘 琼 ，王 轶 ，李小磊 ，石 鹏

（1.中国地质环境监测院，北京 100081；2.内蒙古地质环境监测院，内蒙古 呼和浩特 010020）

植被是生态系统的重要组成部分，是全球陆地生态系统变化和生态环境演化的重要指示因子[1-3]，对水平衡、地球表面能量交换和生物地球化学循环至关重要[4-5]。植被生长既取决于植被类型本身的特性，也受光热条件、水分、土壤的影响。光热条件与气温密切相关[6]。水分的影响较为复杂，由于地下水——土壤——植被——大气间的水分迁移是连续过程[7]，植被生长同时受地下水、土壤水和大气降水等因子的影响。在干旱半干旱区，以灌木、草地、森林及农田等多种形式存在的多样化植被系统有助于保持土壤及水分，监测其长期动态变化，并研究植被演化与气候、水分等相关关系，对于预测土地退化风险意义重大。

目前从区域尺度研究干旱半干旱区植被演化及其影响因子相关性分析的主要手段是遥感技术。Fabricante等[8]研究了巴塔哥尼亚北部草原植被指数与降水的关系，认为干旱的生态系统中，NDVI对降水之间存在滞后响应；Khan等[9]研究了巴基斯坦及其邻近南亚国家植被生长对当前气候变化的相应，认为NDVI值变化取决于土地覆盖类型及降雨量。国内相关研究主要集中在柴达木盆地、鄂尔多斯高原、黑河流域等西北干旱半干旱地区。杨泽元等[10]研究了陕北风沙滩地区植被生长与地下水位埋深的关系，确定了风沙滩地区生态安全水位埋深为1.5～5 m；张二勇等[11]定量研究了鄂尔多斯盆地内蒙古能源基地生态植被与地下水关系，分析了低湿植被、低湿植被与沙地植被过渡区的水位埋深阈值分别为3，4.8 m；金晓媚等[12-13]分析了海流兔河流域和柴达木盆地植被与地下水位、气象因子、地形地貌之间的关系，认为不同植被类型对地下水依赖性大；王旭升等[14]研究了鄂尔多斯高原植被盖度与浅层地下水的依赖关系，认为植被盖度同时受气候条件和地下水的影响；赵捷等[15]研究了黑河上中游流域长序列植被盖度数据，认为植被受气温的影响大于降水，而对降水的响应更为迅速；李元春等[16]在甘南和川西北地区基于植被NDVI值分析了草地植被变化趋势及驱动因素，认为温度和降水对研究区草地NDVI值的变化主要呈正向驱动，且NDVI值与温度呈正相关的面积大于降水。

西辽河平原位于内蒙古东部平原，属中温带半干旱区的农牧交错带。平原西部广泛分布风积沙地、沙丘，是我国半干旱区典型的由沙漠向草原过渡的荒漠化草原区。目前有关西辽河平原的植被演化及其与环境因子的定量关系研究较少，影响因子考虑不够全面。本文基于2000—2019年16天合成的MODIS NDVI数据，分析了该研究区近20年来植被生长变化趋势，并从影响植被生长的水热条件出发，开展NDVI与气象因素（降水、气温）、土壤湿度、地下水埋深等因子的定量关系研究，为该流域生态环境保护治理及水资源合理开发利用提供技术支撑。

1 研究区概况

西辽河平原地势平坦开阔，由西南向东北逐渐倾斜，主要地貌单元为河谷冲击平原和风积沙丘。北部主要为草原牧区，中部及南侧边缘为农区，沙地主要分布在西南部，面积约5.76×104km2。该地区多年平均降水量375 mm，多年平均蒸发量1 900 mm，降水主要集中在6—9月。区内天然植被主要有沙蒿、沙柳、苔草、芨芨草、柠条等，受水热条件影响，生长季为6—9月。主要过境河流有西辽河及新开河，1999年断流至今。自20世纪60年代起，地下水开采致使沙地周边的地下水水位下降，植被生长受到影响，荒漠化程度加重。近年来经过全面综合治理，该区域内沙地面积逐渐变小，生态环境持续向好发展。

2 研究数据及方法

2.1 数据及处理

（1）遥感数据及预处理

归一化差值植被指数（NDVI，normalized difference vegetation index）是目前广泛应用于反映植被生长状况的指标[17-20]，NDVI值越大，说明植被生长状况越好，NDVI值越小，说明植被生长状况越差。NDVI数据可直接从MOD13Q1遥感数据中提取，数据来源于美国NASA网（https://ladsweb.modaps.eosdis.nasa.gov/），空间分辨率250 m，时间分辨率16 d，数据时间为2000—2019年。覆盖研究区需要4幅，2000—2019年1—12月共1 236幅。在MRT软件中提取NDVI后进行投影转换，为了消除云、大气及太阳高度角的影响，采用最大值合成法对每个月的NDVI数据进行处理，即求取当月各栅格的最大值。

（2）地下水水位埋深和气象数据

西辽河平原地下水埋深数据采用2019年国家级地下水监测工程监测数据，监测频率为1次/h，研究区内及周边地下水监测井共125个。气象数据来源于中国气象数据共享网，西辽河平原内及周边气象站站点共9个，获取2000—2019年月均降水量、温度等相关指标。监测井和气象站点位置见图1。

图1 研究区气象站点和地下水监测站点分布图Fig.1 Distribution of meteorological stations and groundwater monitoring stations in the study area

（3）土壤湿度数据

GLDAS（the Global Land Data Assimilation System）是美国航天局NASA提供的全球尺度的陆面同化数据集，包含了全球降雨量、水分蒸发量、地表径流、地下径流、土壤湿度、地表积雪的分布以及温度和热流分布等数据。本文选取GLDAS NOAH模型中土壤湿度月解数据，分辨率为0.25 °×0.25 °。本文选用0～10 cm表层土壤湿度数据，单位为kg/m2。

（3）土地利用数据

土地利用数据来源于中国科学院地理科学与资源研究所资源环境科学与数据中心，本文选取2000年和2018年的土地利用类型数据，该数据是基于Landsat 8遥感影像，通过人工目视解译生成，数据类型是1 km栅格。通过重分类，将研究区土地利用数据分为10 类，分别是耕地、草地、建设用地、沙地、沼泽地、林地、水域、滩地、盐碱地和裸土地（图2）。

图2 2000年和2018年研究区土地利用类型分布图Fig.2 Distribution of land use types in the study area in 2000 and 2018

2.2 研究方法

（1）一元线性回归趋势法

用一元线性回归趋势分析方法分析研究区的植被变化趋势，即对2000—2019年的生长季（6—9月份）平均NDVI值求斜率，模拟每一个栅格的变化趋势：

式中：n——样本的数量，n=20；

i——年序号，i的值取1，2， ···, 20；

NDVIi—— 第i年的NDVI值；

Slope——2000—2019年每年平均NDVI的线性回归率。

计算得到的斜率值可以反映在20年的时间序列中，研究区平均NDVI值的变化趋势，每个像元点的slope值均代表该点20年的变化趋势。若slope＞0，说明NDVI增加，植被呈增长趋势，植被状况变好；若slope=0，说明NDVI不变，植被基本无变化；若slope＜0，说明NDVI减小，植被呈减小趋势，植被状况变差。

(2) Mann-Kendall突变检测

Mann-Kendall法是一种非参数统计检验方法，不受少数异常值干扰，且不需要样本遵从一定的分布，被广泛应用于突变检测[21-23]。本文利用该方法检测NDVI变化过程中的突变点，检测方法为：

根据西辽河平原年均NDVI时间序列（x1,x2,···xn），其秩序列为：

定义统计量：

式中：E(Sk)——Sk的均值；

Var(Sk)——Sk的方差。

UFk计算结果组成UF曲线，同样以此方法引用到反序列中，计算得到UBk， 使得UBk=-UFk（k=1，2，···,n）。若UFk和UBk的值大于0，表明序列呈上升趋势，反之呈下降趋势。当曲线超出置信度区间时，表明上升或下降趋势显著。如果UFk和UBk两条曲线出现交点，且交点在临界线之间，那么该点即为突变点。

（3）相关性分析法

对西辽河平原植被NDVI与气象因子（降水、温度、湿度）的相关性进行分析，相关性系数rxy为：

式中：——因子样本值的平均值。

rxy是因子x和y之间的相关系数，表示两因子相关性的统计指标，该指标的取值区间为[-1,1]，且rxy＞0为正相关，rxy的 绝对值越大，表明x和y的相 关性越密切。

3 结果与讨论

3.1 植被时空变化特征

3.1.1 植被随时间变化趋势

分析2000—2019年每年生长季的全域NDVI均值表明，西辽河平原植被NDVI值整体表现出波动升高的趋势，该时段内NDVI值变化大致分为3个阶段：2000—2005年NDVI值呈较平稳的上升趋势，2006—2012年NDVI值上下起伏变化较大，整体处于上升趋势，2013—2019年NDVI值呈平稳的上升趋势（图3）。2007年植被NDVI值达到极小值0.41，2019年植被NDVI值达到最大值0.56。对2000—2019年年均NDVI值进行Mann-Kendall突变检验，UF值始终大于0，表明20年来NDVI值呈上升趋势。2012年UF曲线超出置信区间，说明2012年之前由于震荡较大，上升趋势不显著，2012年之后呈显著上升趋势。UF和UB统计量相交于2014年，但由于交点不在置信区间范围内，因此NDVI值不存在显著突变点（图4）。

图3 研究区2000—2019年NDVI年际变化趋势图Fig.3 Interannual variation trend of NDVI in the study area from 2000 to 2019

图4 研究区NDVI值突变检验曲线图Fig.4 Mutation test curve of NDVI in the study area

3.1.2 植被空间分布特征及变化趋势

研究区2000—2019年生长季的NDVI值求平均值，得到西辽河平原20年来平均NDVI的空间分布图，见图5(a)。西辽河平原NDVI在空间上呈现“东高西低”的分布特征，结合土地利用类型，NDVI值较大

图5 研究区2000—2019年平均NDVI值和植被变化趋势分布图Fig.5 Spatial distribution maps of NDVI and variation trends of vegetation in the study area during 2000—2019

分析研究区2000—2019年生长季平均NDVI值一元线性回归趋势（式1），得到2000—2019年NDVI值变化趋势分布图，见图5(b)，研究区内大部分地区植被呈增长趋势，面积占比约92.5%，小部分地区植被呈减少趋势，面积占比约7.5%。结合用地类型分析可知，草地是植被减少的主力区域，说明草地退化较为明显。

3.2 气候对植被的影响

本文对研究区20年来月均气象数据（降水量、气温）和相应的全域NDVI月均值进行相关性分析，见图7(a)(b)。结果表明：研究区NDVI值与降水量呈正相关关系，相关系数r=0.86，NDVI值随降水量的增加而升高。而NDVI值与气温的相关关系具有分段特点，当气温在10℃以下，NDVI值趋于水平无明显变化，在0.2上下浮动；当气温大于10℃时，随着温度的升高NDVI值呈增大趋势。通过NDVI值和气温的线性拟合，二者呈正相关关系，相关系数r=0.78。的区域分布在平原中部的西辽河、新开河沿岸地带以及平原东南部地区，主要地类是耕地。NDVI值较小的区域主要分布在研究区西部及北部的科尔沁沙地，主要地类是沙地。不同地类的NDVI值由大到小依次为耕地＞林地＞沼泽地＞滩地＞草地＞盐碱地＞沙地，且随季节变化趋势一致，夏季＞秋季＞春季＞冬季（图6）。

图6 研究区不同土地利用类型NDVI值随季节变化趋势图Fig.6 Seasonal variation trend of NDVI of different land use types in the study area from 2000 to 2019

3.3 土壤湿度对植被的影响

植被生长状况与土壤水分条件有关，土壤湿度在陆地表面的水文过程中起着至关重要的作用，表层土壤湿度湿化还是干化与植被的绿化还是恶化之间存在相关关系[24]。分析研究区内2000—2019年月均NDVI值及对应月份0～10cm土壤湿度的相关性，见图7(c)。结果显示：NDVI值随表层土壤湿度的增加而增大，存在较强正相关关系，r=0.81。

图7 研究区NDVI值与降水量、气温、表层土壤湿度的相关性分析Fig.7 Correlation of NDVI with precipitation, air temperature and surface soil moisture in the study area

3.4 地下水埋深对植被的影响

本文利用研究区及周边125个地下水潜水监测井的2019年6月监测数据，与对应的NDVI值进行相关性分析。研究区处在农牧交错带，耕地和牧草地主要靠抽地下水进行灌溉，因此生长季植被指数最高的耕地，地下水位埋深反而最大。为了消除人为抽水因素影响，本文选择研究区内的天然植被区，研究植被与地下水埋深的关系。

6月天然植被区地下水埋深为0～17.45 m，将地下水埋深以0.2 m为步长，统计每个步长对应的NDVI平均值，生成NDVI值随地下水埋深变化趋势图（图8）。地下水埋深为0～4 m时，NDVI值随埋深的增加而增大；地下水埋深4～10 m时，NDVI值趋于稳定，在0.4上下浮动；地下水埋深大于10 m时，NDVI值随埋深的增加而急剧减小。NDVI值在地下水埋深3 m左右达到最大值0.43。

图8 研究区天然植被NDVI值和地下水埋深关系图Fig.8 Relationship between NDVI of natural vegetation and groundwater level depth in the study area

西辽河平原降水主要集中在植被生长季，因此降水对该地区植被生长有很重要的影响。本文在讨论生长季植被与地下水埋深相关性时，没有排除降水对植被生长的影响。降水和地下水对植被的双重影响机制还需要进一步研究。

3.5 人类活动对植被的影响

人类活动如土地开垦、植树造林、城市扩张等，通过改变土地利用类型，对地表植被分布和变化趋势产生明显的影响。

本文统计了研究区2000年和2018年不同土地利用类型的面积及变化（表1）。结果显示，耕地、林地面积显著增加，分别增加了1 873，1 312 km2，草地和水域面积显著减小，分别减小了4 187，583 km2。结合近20年来NDVI值的变化趋势，耕地和林地面积的显著增加是NDVI值整体呈增加趋势的主要原因之一，说明耕地开垦及植树造林在一定程度上改善了当地的生态环境。但是由于耕地面积增加，会导致灌溉需水量增加，地下水位降低，草地和水域面积呈减少趋势，不利于形成生态环境改善的良性循环。

表1 研究区2000和2018年土地利用类型面积统计Table 1 Area of different land use types in the study area in 2000 and 2018 /km2

4 结论

（1）研究区20年来全域NDVI均值整体呈增加趋势，表明植被生长变好，最高值0.56，最低值0.41。在空间上呈现“东高西低”的分布特征，不同用地类型NDVI值由大到小依次为：耕地＞林地＞沼泽地＞滩地＞草地＞盐碱地＞沙地。

（2）研究区内大部分地区NDVI值呈增长趋势，面积为5.33×104km2，占比92.50%，小部分地区NDVI值呈减少趋势，面积为0.43×104km2，占比7.50%。

（3）研究区NDVI值与降水量、气温、土壤湿度关系密切，呈正相关关系，相关系数分别为0.86，0.78，0.81。降水对植被生长影响最大。

（4）地下水埋深为0～4 m时，NDVI值随埋深的增加而增大；地下水埋深4～10 m时，NDVI值趋于稳定，在0.4上下浮动；地下水埋深大于10 m时，NDVI值随埋深的增加而急剧减小。NDVI值在地下水埋深为3 m左右达到最大值0.43。

（5）人类活动如土地开垦、植树造林是近20年来NDVI值呈增加趋势的主要原因之一，在一定程度上改善了当地生态环境。