生态学报  2023, Vol. 43 Issue (1): 419-431

文章信息

苏日罕, 郭恩亮, 王永芳, 银山, 顾锡羚, 康尧, 包山虎
SU Rihan, GUO Enliang, WANG Yongfang, YIN Shan, GU Xiling, KANG Yao, BAO Shanhu
1982-2020年内蒙古地区极端气候变化及其对植被的影响
Extreme climate changes in the Inner Mongolia and their impacts on vegetation dynamics during 1982-2020
生态学报. 2023, 43(1): 419-431
Acta Ecologica Sinica. 2023, 43(1): 419-431
http://dx.doi.org/10.5846/stxb202112073469

文章历史

收稿日期: 2021-12-07
网络出版日期: 2022-08-29
1982-2020年内蒙古地区极端气候变化及其对植被的影响
苏日罕1,2 , 郭恩亮1,2 , 王永芳1,2 , 银山1,3 , 顾锡羚1,2 , 康尧1,2 , 包山虎1,3     
1. 内蒙古师范大学地理科学学院, 呼和浩特 010022;
2. 内蒙古自治区蒙古高原灾害与生态安全重点实验室, 呼和浩特 010022;
3. 内蒙古自治区遥感与地理信息系统重点实验室, 呼和浩特 010022
摘要: 内蒙古地处生态环境脆弱区,对气候变化尤为敏感。在全球气候变暖背景下,探究极端气候变化及其影响显得尤为重要。基于内蒙古地区115个气象站点1982-2020年的逐日气象数据,从强度、持续时间、频率3个维度出发计算了18个极端气候指数,在综合分析极端气候的时空变化特征的基础上,运用地理探测器和皮尔逊相关分析方法,定量评估极端气候对该区植被的影响。结果表明:(1)极端暖指数均呈增加趋势,说明1982-2020年期间内蒙古地区极端偏暖现象增多。(2)持续干旱日数与持续湿润日数呈减少趋势,说明39年来内蒙古地区连续性无降水天数和降水天数均减少。(3)极端气候指数与归一化植被指数(NDVI)的相关关系表现出明显的空间异质性,表明内蒙古不同区域NDVI对各极端气候指数的响应程度不同。(4)因子探测器结果表明极端降水指数相对于极端气温指数来说,对内蒙古植被生长变化的影响较大。研究结果可为内蒙古地区防灾减灾与生态修复工程提供一定的科学依据。
关键词: 极端气候    归一化植被指数(NDVI)    时空变化    地理探测器    内蒙古    
Extreme climate changes in the Inner Mongolia and their impacts on vegetation dynamics during 1982-2020
SU Rihan1,2 , GUO Enliang1,2 , WANG Yongfang1,2 , YIN Shan1,3 , GU Xiling1,2 , KANG Yao1,2 , BAO Shanhu1,3     
1. College of Geography Science, Inner Mongolia Normal University, Hohhot 010022, China;
2. Key Laboratory of Disaster and Ecological Security on Mongolian Plateau, Hohhot 010022, China;
3. Key Laboratory of Remote Sensing and Geographic Information System of Inner Mongolia Autonomous Region, Hohhot 010022, China
Abstract: In recent years, a growing number of reports and studies have shown that extreme weather events will become more frequent, which will cause serious impacts and losses to society, economy and human life. In the context of global warming, it is especially significant to explore the characteristics of extreme climate change and its impact on ecosystem. Vegetation plays a key role in the terrestrial ecosystem. The Normalized Difference Vegetation Index (NDVI) is the best indicator for monitoring the growth state and spatial distribution of terrestrial vegetation, and is often used to reflect the response of vegetation dynamics to climate change in the region. The Inner Mongolia is located in a fragile ecological environment and is particularly sensitive to climate change. Therefore, it is necessary to explore the dynamic change of extreme climate in Inner Mongolia and its impact on vegetation, which can provide a scientific basis for disaster mitigation and ecological restoration projects in the Inner Mongolia. Based on the daily meteorological data of 115 meteorological stations in the Inner Mongolia from 1982 to 2020, this paper selected and calculated 18 extreme climate indices from three dimensions of intensity, duration, and frequency, and comprehensively analyzed the temporal and spatial variation characteristics of extreme climate. On this basis, the geographical detectors and Pearson correlation analysis methods were used to quantitatively evaluate the impact of extreme climate on vegetation in the region. The results showed that: (1) the extreme warm temperature indices all showed an increasing trend, indicating that extreme warm phenomenon in the Inner Mongolia has increased during the period 1982-2020. (2) Both consecutive dry days and consecutive wet days showed a decreasing trend, indicating that the continuous drought and precipitation events in the Inner Mongolia have decreased in the past 39 years. (3) The correlation between the extreme climate indices and the NDVI showed obviously spatial heterogeneity, which is shown as follows: the NDVI and extreme warm indices are positively correlated in the east, but negatively correlated in the central region; the NDVI was significantly positively correlated with extreme precipitation intensity and frequency indices in the central Inner Mongolia. (4) The factor detector results showed that the extreme precipitation index had a greater impact on the vegetation growth and changes than the extreme temperature index in the Inner Mongolia. In general, the paper can provide the theoretical evidence for disaster mitigation and ecological environment conservation in the Inner Mongolia.
Key Words: extreme climate    NDVI    spatio-temporal change    geographical detector    the Inner Mongolia    

极端气候事件是指特定时段内某类气候要素统计量值显著偏离平均值、且超出其观测统计量值区间上下限的阈值的事件[1]。与平均态的气候变化相比, 极端气候事件由于其突发性和破坏性强, 会给社会、经济和人类生活造成严重的影响和损失。IPCC第六次评估报告中指出, 在全球变暖进一步加剧的背景下, 热浪、强降水、干旱等极端气候事件发生的频次和强度将增加[2]。一系列研究又指出欧洲、蒙古高原、美国、伊朗等地区均出现极端暖事件逐渐增多, 极端冷事件逐渐减少的现象[38]。因此, 极端气候的研究已成为当今社会的热点话题。

联合国政府间气候变化专门委员会第六次评估报告(IPCC6)中指出,气候变化监测和指标专家组(ETCCDI)定义了27个代用气候指数, 主要集中在对极端气候的描述上[9]。目前, 国内外不同研究小组都一致认同代用极端气候指数对极端气候变化监测的重要性[10]。除了极端气候事件监测外, 其对陆地生态系统的影响研究也是热点问题之一[11]。植被连接着土壤、水圈和大气圈[12], 其呼吸作用、光合作用和植被的蒸散发等对生态系统能量收支、水热平衡和生物循环有着重要作用[13]。归一化植被指数(NDVI)是监测陆地植被生长状态及空间分布变化的最佳指示因子[14], 常用于反映区域内植被动态对气候变化的响应[15]。例如, 高滢等对陕西省极端气候指数与NDVI进行相关分析, 发现NDVI对极端气候的响应具有空间差异性, 比如NDVI与日最低气温极大值(TNx)在陕西省北部呈负相关, 在南部呈正相关[16]。张晋霞等发现, 气温日较差(DTR)和冷夜日数(Tn10p)是影响新疆大部分地区植被NDVI的主导指数[17]

内蒙古位于中国北部边疆, 其植被类型多样, 生态环境脆弱, 对气候变化尤为敏感。极端气候事件会对生态系统土壤水分和养分状况、植物生长发育和生理特性、群落结构、生产力和碳循环产生影响[18]。目前, 有关内蒙古地区极端气候变化的研究已有不少, 但是从极端气候的强度、持续时间、频率等多种维度深入探究其演变规律及对植被生长影响的研究较少。因此本文从极端气候的强度、持续时间和频率角度出发, 选取18个极端气候指数对内蒙古地区极端气候的时空演变规律进行研究, 在此基础上, 利用皮尔逊相关系数法和地理探测器模型, 分别从线性和空间异质性角度入手, 探究极端气候对植被动态变化的影响, 研究结果对该区生态环境的保护与筑牢区域生态安全具有重要意义。

1 材料与方法 1.1 研究区概况

内蒙古自治区(97°12′—126°04′E, 37°24′—53°23′N)位于我国北部边疆地区, 总面积为118.3万km2。内蒙古地形复杂、地貌单元多样, 平均海拔1000m, 地形地势由西向东、从南到北逐渐下降(图 1)。内蒙古位于北半球中高纬度内陆区, 属于西北干旱、半干旱气候向东南沿海湿润、半湿润季风气候的过渡带[19]。年均温为-4—9℃, 具有夏季短暂而炎热、冬季寒冷而漫长、春季风大少雨、秋季气温剧降, 昼夜温差大, 日照时间充足等特点。降水量少且不均匀, 全区降水量平均为375mm, 由东北向西南递减。

图 1 研究区概况 Fig. 1 Overview of the study area
1.2 数据来源与处理 1.2.1 NDVI数据

本研究使用的NDVI数据来源于国家地球系统科学数据中心提供的NOAA CDR NDVI数据集(http://www.geodata.cn)。选取的数据集时间跨度为1982—2020年, 空间分辨率为5km, 时间分辨率为1d。最大值合成法(MVC) 可以从多时间数据中选择最高NDVI值来代表某像元值, 并具有消除云、大气与太阳高度角等干扰的优点[20]。因此, 本研究利用MVC对NDVI数据集进行合成处理, 构建了1982—2020年内蒙古地区的NDVI年最大值来监测该区植被动态及其对极端气候的响应。

1.2.2 气象数据

气象数据来源于内蒙古气象局115个站点的逐日最高气温、日最低气温和日降水量数据, 若有缺失的站点数据, 利用历史同时期多年平均值进行了插补处理。使用Rclimdex1.0对原始数据进行质量分析后计算了1982—2020年18个极端气候指数(包括10个极端气温指数和8个极端降水指数), 再用反距离权重法进行了空间插值, 并将选取的极端气候指数从强度、持续时间和频率3个维度进行了归类[21](表 1)。

表 1 极端气候指数及定义 Table 1 Definitions of extreme climate indices
类别Category 代码ID 名称Name 定义Definition 单位Unit
极端气温强度指数 TNn 日最低气温极小值 每月内日最低气温的最小值
Extreme temperature intensity TXn 日最高气温极小值 每月内日最高气温的最小值
indices TNx 日最低气温极大值 每月内日最低气温的最大值
TXx 日最高气温极大值 每月内日最高气温的最大值
DTR 气温日较差 年内日最高气温与最低气温的差值
极端气温持续时间指数
Extreme temperature duration index
GSL 生长季长度 日平均气温第一次连续6天以上大于5 ℃至日平均气温第一次(6月1日后)连续6天小于5 ℃的日数 d
极端气温频率指数 Tx90p 暖昼日数 日最高气温>90%分位值的日数 d
Extreme temperature frequency Tn90p 暖夜日数 日最低气温>90%分位值的日数 d
indices Tx10p 冷昼日数 日最高气温<10%分位值的日数 d
Tn10p 冷夜日数 日最低气温<10%分位值的日数 d
极端降水强度指数 Rx1day 1日最大降水量 每月最大1日降水量 mm
Extreme precipitation intensity Rx5day 5日最大降水量 每月连续5日最大降水量 mm
indices PRCPTOT 雨日降水总量 雨日(日降水量≥1mm)降水总量 mm
极端降水持续时间指数 CDD 持续干燥日数 日降水量<1mm的最长连续日数 d
Extreme precipitation
durationindices
CWD 持续湿润日数 日降水量≥1mm的最长连续日数 d
极端降水频率指数 R10 中雨日数 每年日降水量≥10mm的总日数 d
Extreme precipitation R20 大雨日数 每年日降水量≥20mm的总日数 d
frequency indices R25 暴雨日数 每年日降水量≥25mm的总日数 d
1.3 研究方法 1.3.1 森氏斜率(Theil-Sen)和曼-肯德尔(Mann-Kendall)显著性检验

Theil-Sen斜率是一种非参数估计方法, 常用于分析某一要素在长时间序列中的变化趋势[22]。相对于简单的线性回归分析, 该方法可以减少异常值的影响[23]。因此, 本文用此方法计算了极端气候指数与NDVI的空间变化趋势。计算公式如下:

(1)

式中, median表示中位数函数, β表示某极端气候指数或NDVI的变化趋势, xixj表示不同年份的序列数据, titj为时间序列。当β>0时, 表示某极端气候指数或NDVI呈上升趋势;β<0时, 则表示存在下降趋势。

Mann-Kendall检验(简称MK检验)是一种非参数统计检验方法, 由于其不需要样本遵循一定的分布, 也不受少数异常值的干扰, 在水文、气象等要素时间序列趋势变化中广泛应用[16]。计算公式如下:

假设时间序列X1, X2, …, Xn。定义检验统计量S:

(2)
(3)

式中, xixj为时间序列数据, n为数据个数。当n>10时, 统计量S为正态分布, 标准化检验统计量Z通过下式计算:

(4)
(5)

式中, 当|Z|>1.65、1.96和2.58时, 表示变化趋势分别通过了置信度90%、95%和99%的显著性检验, 说明存在显著变化趋势。

1.3.2 相关分析

本文采用Pearson相关分析方法分析了极端气候指数与NDVI之间的相关关系。计算公式如下:

(6)

式中, Rxy为变量xy的相关系数, xi为第i年的NDVI值, X为NDVI多年平均值, yi为第i年极端气候指数值, Y为极端气候指数多年平均值。相关系数Rxy的取值范围为[-1, 1], |R|越大则表示变量之间的相关性越强。

1.3.3 地理探测器

地理探测器是探测空间分异性, 以及揭示其背后驱动力的一组统计学方法[24], 包括因子探测器、交互探测器、风险探测器和生态探测器。本文采用因子探测器测度内蒙古地区18个极端气候指数对植被NDVI变化的解释力。计算公式如下:

(7)

式中, Q为极端气候指数对NDVI变化趋势的解释力, Q⊂[0, 1], q值越大表示解释力越强;h(1, 2, …, L)为变量分类或分区;NhN分别为层h和区域单元数;σh2σ2分别是h层和区域内Y值的方差。

2 结果与分析 2.1 内蒙古地区极端气候指数变化特征 2.1.1 极端气候指数时间变化

从极端气温指数的年际变化(图 2)发现, 极端气温强度指数日最低气温极小值(TNn)、日最高气温极小值(TXn)分别以0.192℃/10a、0.074℃/10a的速率减小, 而日最低气温极大值(TNx)和日最高气温极大值(TXx)分别以0.521℃/10a、0.561℃/10a的速率增加。但是日气温极大值的增加速率明显大于日气温极小值的减小速率, 可能这种不对称的变化是导致气温日较差(DTR)缓慢减小(0.018℃/10a)的原因。生长季长度(GSL)以3.559d/10a的速率增加, 说明近39年内蒙古地区极端气温持续时间呈增加趋势。暖昼日数(Tx90p)和暖夜日数(Tn90p)分别以1.977d/10a、2.285d/10a的速率增加, 而冷昼日数(Tx10p)、冷夜日数(Tn10p)分别以1.095d/10a、1.512d/10a的速率减少, 说明内蒙古地区出现暖昼日数和暖夜日数的频率增加, 而出现冷昼日数和冷夜日数的频率减少。从昼夜变化中看, 内蒙古夜间变化速率比白天大, 可以进一步说明, 夜间变暖可能导致了内蒙古地区气温日较差呈减少趋势。

图 2 极端气温指数时间变化趋势 Fig. 2 Temporal variation trend of extreme temperature indices TNn: 日最低气温极小值Monthly minimum value of daily minimum temperature; TXn: 日最高气温极小值Monthly minimum value of daily maximum temperature; TNx: 日最低气温极大值Monthly maximum value of daily minimum temperature; TXx: 日最高气温极大值Monthly maximum value of daily maximum temperature; DTR: 气温日较差Diurnal temperature range; GSL: 生长季长度length of growing season; Tx90p: 暖昼日数Warm days; Tn90p: 暖夜日数Warm nights; Tx10p: 冷昼日数Cold days; Tn10p: 冷夜日数Cold nights

从极端降水指数的年际变化(图 3)看出, 1日最大降水量(Rx1day)以0.238mm/10a的速率增加, 5日最大降水量(Rx5day)以0.663mm/10a的速率减少, 说明内蒙古地区极端降水呈单日降水强度增加, 连续降水强度减少趋势。雨日降水总量(PRCPTOT)则以1.317mm/10a的速率增加。在1982—2020年期间持续干燥日数(CDD)以2.229d/10a的速率减少, 并且持续湿润日数(CWD)也以0.09d/10a的速率减少。总的来说, 39年来内蒙古地区连续性无降水天数和降水天数均在减少。中雨日数(R10)、大雨日数(R20)和暴雨日数(R25)呈缓慢增加趋势, 增加速率分别为0.008d/10a、0.003d/10a、0.001d/10a。由此可见, 1982—2020年内蒙古地区极端降水的频率呈增加趋势, 其中雨日数增加最明显。除此之外, 发现Rx1day、Rx5day、PRCPTOT、R10、R20和R25等极端降水指数最大值都出现在1998年, 可能与1998年气候异常导致的我国大部分地区暴雨偏多有关[2526]

图 3 极端降水指数时间变化趋势 Fig. 3 Temporal variation trend of extreme precipitation indices Rx1day: 1日最大降水量Monthly maximum one-day precipitation; Rx5day: 5日最大降水量Monthly maximum five-day precipitation; PRCPTOT: 雨日降水总量Total amount of precipitation on rainy days; R10:中雨日数Number of days per year with precipitation amount≥10mm; R20:大雨日数Number of days per year with precipitation amount≥20mm; R25:暴雨日数Number of days per year with precipitation amount≥25mm; CDD: 持续干燥日数Consecutive dry days; CWD: 持续湿润日数Consecutive wet days
2.1.2 极端气候指数空间变化

从极端气温指数的空间变化(图 4)发现, TNn呈上升和下降的面积分别占总面积的35.5%和64.5%, 其中仅5个和8个站点通过了显著性检验(P<0.05)。TXn呈显著上升和显著下降的站点仅有2个。TNx呈增加的面积占99.5%, 其中45个站点通过显著性检验(P<0.01)。TXx呈增加的面积占99.6%, 31个站点呈极显著上升趋势(P<0.01)。气温日较差(DTR)的空间变化趋势具有明显差异, 呈上升和下降的面积各占37.9%和62.1%。其中有22个站点呈极显著(P<0.01)上升趋势, 主要分布于赤峰市和鄂尔多斯市, 有29个站点呈极显著的下降趋势, 主要分布于阿拉善盟和巴彦淖尔市和兴安盟的部分地区。生长季长度(GSL)呈增加的面积占97.8%, 共有22个站点通过了显著性检验(P<0.01), 零星分布于全区, 说明39年来内蒙古地区极端气温持续时间普遍增加。Tx90p和Tn90p在99%地区呈增加的趋势且分别有97和91个站点通过了显著性检验(P<0.01)。相反, Tx10p和Tn10p呈减少的面积各占99%和97%, 且分别有52和68个站点呈极显著的减少趋势(P<0.01)。整体而言, 1982—2020年期间内蒙古地区极端气温呈上升趋势。

图 4 极端气温指数空间变化趋势 Fig. 4 Spatial variation trend of extreme temperature indices

从极端降水指数的空间变化(图 5)看出, Rx1day呈上升和下降的面积分别占48%和52%, Rx5day呈上升和下降的面积分别占41.7%和58.3%, PRCPTOT呈上升和下降的面积分别占36.3%和63.7%。其中Rx1day、Rx5day与PRCPTOT呈显著上升(P<0.05)趋势的站点分布于西部的阿拉善盟和巴彦淖尔市, 而显著下降的站点分布于东部的赤峰市和通辽市。说明39年来内蒙古西部地区的极端降水强度在增加, 而中部和东部地区的极端降水强度在减少。CDD呈减少的面积占86%, 呈增加的面积占14%, 但是变化并不显著。CWD呈减少的面积占69%, 主要分布于中部和东部地区, CWD呈增加的面积占31%, 主要分布于西部地区。说明39年来内蒙古西部干旱区的连续性降水天数呈增加趋势, 而中部和东部地区的连续性降水天数呈减少趋势。R10、R20、R25呈增加的面积各占61%、65%、57%, 主要分布在西部地区, 呈减少的面积各占39%、35%、43%, 主要分布在中部地区。综上所述, 1982—2020年内蒙古西部干旱区极端降水的强度、持续湿润时间、频率均在增加, 中部典型草原区极端降水的强度、持续湿润时间、频率均在减少, 表明西部干旱区的干旱程度有所减轻。

图 5 极端降水指数空间变化趋势 Fig. 5 Spatial variation trend of extreme precipitation indices
2.2 内蒙古地区NDVI变化特征

空间分布上, 39年来内蒙古地区NDVI呈自东北向西南逐渐降低的空间格局(图 6)。NDVI高值主要分布于东部地区, 低值主要分布于西部地区。主要是因为干湿度地带性引起了内蒙古植被类型从东到西依次为森林、草甸草原、典型草原、荒漠草原、戈壁荒漠植被类型分布导致的空间差异。

图 6 内蒙古NDVI空间分布及变化趋势图 Fig. 6 Distribution and spatial variation trend of NDVI NDVI: 归一化植被指数Normalized Difference Vegetation Index

在空间变化趋势上(图 6), 根据Sen斜率和MK显著性检验结果划分出了六个等级: 极显著减少(P<0.01)、显著减少(P<0.05)、不显著减少、不显著增加、显著增加(P<0.05)和极显著增加(P<0.01)。结果显示内蒙古地区NDVI呈增加的面积占总面积的76%, 其中呈极显著增加趋势的占64%, 主要分布于内蒙古东部的呼伦贝尔市、通辽市、赤峰市和西南部乌兰察布市、呼和浩特市、鄂尔多斯市。NDVI呈减少的面积占总面积的24%, 其中0.04%的地区通过了显著性检验, 减少趋势并不明显。NDVI的变化趋势具有明显空间差异, 一方面是由于不同地区水热条件和植被类型不同, 从而导致气候变化对植被NDVI产生不同的影响, 引起植被变化的空间差异。另一方面也说明“三北防护林”等生态建设工程取得了显著成效[27]

2.3 极端气候指数对NDVI变化的影响力

从因子探测器结果发现(图 7), 内蒙古地区不同极端气候指数对植被NDVI变化趋势的因子解释力(Q值)排序为: R25>Rx1day>R20>R10>Rx5day> PRCPTOT>GSL>TXn>TNx>CWD>CDD>TXx>DTR>Tx90p>TNn>Tn10p>Tn90p>Tx10p, 且均通过了显著性检验。其中, R25、Rx1day、R20、R10、Rx5day、PRCPTOT等极端降水指数的因子解释力均在0.3以上, 说明内蒙古地区植被变化对极端降水的强度和频率最为敏感。Tx90p、Tn90p、Tx10p、Tn10p、TNn等极端气温指数的因子解释力较小, 均在0.1以下, 说明内蒙古地区植被变化对极端气温的敏感性不强。

图 7 内蒙古地区不同极端气候指数对植被NDVI变化的Q值 Fig. 7 Q valve of extreme climate indices to NDVI changes in Inner Mongolia Q Valve: 因子解释力The contribution rate of extreme climate indices to NDVI
2.4 内蒙古地区极端气候与NDVI的相关性 2.4.1 极端气温与NDVI的相关性

图 8表示内蒙古地区年NDVI对各极端气候指数的响应程度, 通过的显著性检验均在α=0.05水平。结果显示, NDVI与TNn、TXn的相关性普遍不高, 全区98%的地区相关系数在-0.4-0.4之间。TNx、TXx在内蒙古中部和西部地区与NDVI呈显著负相关, 可能是由于该地区下垫面为典型草原和荒漠草原区, 加之气候干旱且水分条件较差, 因此日气温的升高一定程度上加快了植被的蒸腾和土壤水分蒸发, 从而对植被的生长产生负面影响[28];在内蒙古东部地区与NDVI呈显著正相关, 有研究指出大兴安岭北部森林区植被气候生产潜力与气温呈正相关且气温是影响该区植被气候生产潜力的主导因子[29], 由此可以确定, 日气温的升高促进了该区植被生长。NDVI与Tx90p、Tn90p呈显著正相关的地区主要东部地区, 呈负相关的地区主要分布于中西部地区, 但是相关性并不显著。与之相反, Tx10p、Tn10p与NDVI在东部地区呈显著负相关, 在中西部地区呈显著正相关。可能是对于东部地区森林和草甸草原区来说, 暖昼日数和暖夜日数的增加会使土壤温度增加且保持平衡, 从而减轻低温带来的伤害, 增强了植被的碳固定能力、初级生产力及呼吸作用, 为植被创造了良好的生长条件[30];而对于中部典型草原和荒漠草原地区来说冷昼日数和冷夜日数的增加利于土壤保持水分、减少蒸散量, 从而利于该区植被的生长。DTR与NDVI呈显著负相关的地区主要分布于中部地区, 与NDVI呈显著正相关的地区主要分布于东部地区。GSL与NDVI呈显著正相关的地区主要分布于中部和东部地区, 呈负相关的地区未通过显著性检验。

图 8 极端气温指数与年NDVI的相关性 Fig. 8 Correlation between extreme temperature indices and annual NDVI
2.4.2 极端降水与NDVI的相关性

图 9可以看出, NDVI与Rx1day、Rx5day、PRCPTOT、R10、R20、R25在内蒙古中部和西部地区均呈显著正相关(P<0.05), 可能是由于该地区水资源相对匮乏, 极端降水的增加在一定程度上减缓了旱情, 从而促进了植被生长。该结果与上文中计算的因子探测器的结果较为一致。NDVI与Rx1day、Rx5day、PRCPTOT、R10、R20、R25呈负相关的地区主要分布于东部地区, 但未通过显著性检验, 相关性总体较弱。NDVI对极端降水强度和频率的响应具有空间规律性, 具体表现为: NDVI与Rx1day、Rx5day、PRCPTOT、R10、R20、R25在内蒙古中部典型草原区均呈正相关。NDVI与CWD呈显著正相关的地区主要分布于中部地区, 呈负相关的地区主要分布于东部森林区。森林区土壤水分较充足, 因此极端降水持续时间过长可能会导致森林植被根系因水淹而出现无氧呼吸, 影响植被的养分吸收效率, 进而影响植被的正常发育[18]。因此, CWD与NDVI在东部森林区呈显著负相关。此外, 图中显示NDVI与CDD相关性不高。

图 9 极端降水指数与年NDVI的相关性 Fig. 9 Correlation between extreme precipitation indices and annual NDVI
3 讨论

全球正经历以气候变暖为主要特征的变化。本文对内蒙古地区1982—2020年极端气候变化进行研究发现, 内蒙古地区正经历极端偏暖现象, 这与李春兰、雅茹等的研究结果也是一致的[3132]。Alexander等在全球极端气候变化研究中指出, 全球70%的陆地呈现出冷夜日数持续减少而暖夜日数持续增加趋势越发显著[33]。因此, 内蒙古地区符合全球气候变化趋势。在昼夜变化中发现, 内蒙古夜间变化速率比白天大, 故夜间温度的升高可能进一步加速了气候变暖。在极端降水变化中发现, 西部地区极端降水的强度、持续湿润时间、频率均在增加, 说明内蒙古西部干旱区的干旱程度有所减轻, 与前人的研究结果较为一致[34]

对比极端气温和极端降水对NDVI的影响发现, 影响内蒙古地区植被的主要生长因素是降水, 且对极端降水的强度和频率的响应最为敏感, 并呈现出明显的空间特征, 具体表现为: NDVI与极端降水在内蒙古中部典型草原区呈正相关。罗琳玲在内蒙古植被长势与气候驱动分区研究中的发现在森林生态区、草原生态区和荒漠草原生态区3个分区中, 内蒙古草原生态区受降水的驱动力为最强, 这与本研究结果一致[35]

本文也存在一些不足之处, 影响内蒙古植被变化的因素不仅限于气候, 人为因素同样起着重要作用。内蒙古东部地区极端降水在减少, NDVI却呈现显著增加趋势, 这离不开内蒙古东部地区“三北防护林”、“退耕还林(草)”等生态工程的实施[3637]。因此, 在今后研究中持续强化极端气候与人为因素对植被的影响机制研究。由于研究区东西跨度大, 水热条件和下垫面因子也存在较大差异, 尤其是植被类型多样且遵循有显著的空间分布规律, 因此, 极端气候变化对内蒙古不同植被类型区NDVI的影响也是在以后研究的重点。

4 结论

本文利用NDVI数据, 以及内蒙古115个气象站点的逐日最高气温、日最低气温和日降水量数据, 计算出18个极端气候指数, 并采用趋势分析、MK检验、相关性分析、地理探测器等方法, 探究内蒙古地区1982—2020年极端气候的强度、持续时间、频率时空变化特征及其对植被的影响, 得到以下主要结论:

(1) 在时间变化趋势上, TNx、TXx、Tx90p、Tn90p和GSL等极端气温指数均呈增加趋势, 说明近39年内蒙古地区极端高温的强度、持续时间和频率均在增加。极端降水持续时间指数CDD、CWD呈减少趋势, 说明内蒙古地区连续性无降水天数和降水天数均减少。

(2) 在空间变化趋势上, 发现, GSL呈增加的面积占总面积的97.8%, 说明, 内蒙古地区极端气温持续时间普遍增加。从极端降水指数的空间变化发现, 39年来内蒙古西部地区的极端降水强度和频率在增加, 而中部和东部地区的极端降水强度和频率在减少。

(3) 极端气候指数与NDVI的相关性具有明显的空间差异。具体表现为: NDVI与TNx、TXx、Tx90p、Tn90p等极端暖指数在东部多呈正相关, 而在中部地区多呈负相关;NDVI与Rx1day、Rx5day、PRCPTOT、R10、R20、R25等极端降水指数在内蒙古中部地区均呈显著(P<0.05)正相关。

(4) R25、Rx1day、R20、R10、Rx5day、PRCPTOT等极端降水指数的因子解释力均在0.3以上, 说明极端降水的强度和频率对内蒙古地区植被变化有较大影响。

参考文献
[1]
韩丹丹, 穆兴民, 高鹏, 赵广举, 孙文义, 田鹏. 黄土高原地区植被变化及其对极端气候的响应. 水土保持通报, 2020, 40(2): 247-254. DOI:10.13961/j.cnki.stbctb.2020.02.036
[2]
周波涛, 钱进. IPCC AR6报告解读: 极端天气气候事件变化. 气候变化研究进展, 2021, 17(6): 713-718.
[3]
Gruza G, Rankova E, Razuvaev V, Bulygina O. Indicators of climate change for the Russian federation. Weather and Climate Extremes. Dordrecht: Springer Netherlands, 1999, 42(1): 219-242.
[4]
Choi G, Collins D, Ren G Y, Trewin B, Baldi M, Fukuda Y, Afzaal M, Pianmana T, Gomboluudev P, Huong P T T, Lias N, Kwon W T, Boo K O, Cha Y M, Zhou Y Q. Changes in means and extreme events of temperature and precipitation in the Asia-Pacific Network region, 1955-2007. International Journal of Climatology, 2009, 29(13): 1906-1925. DOI:10.1002/joc.1979
[5]
Fathian F, Ghadami M, Haghighi P, Amini M, Naderi S, Ghaedi Z. Assessment of changes in climate extremes of temperature and precipitation over Iran. Theoretical and Applied Climatology, 2020, 141(3-4): 1119-1133. DOI:10.1007/s00704-020-03269-2
[6]
Klein Tank A M G, Konnen G P. Trends in indices of daily temperature and precipitation extremes in Europe, 1946-99. Journal of Climate, 2003, 16(22): 3665-3680. DOI:10.1175/1520-0442(2003)016<3665:TIIODT>2.0.CO;2
[7]
Wang L, Yao Z J, Jiang L G, Wang R, Wu S S, Liu Z F. Changes in climate extremes and catastrophic events in the Mongolian plateau from 1951 to 2012. Journal of Applied Meteorology and Climatology, 2016, 55(5): 1169-1182. DOI:10.1175/JAMC-D-14-0282.1
[8]
Aguilar E, Peterson T C, Obando P R, Frutos R, Retana J A, Solera M, Soley J, García I G, Araujo R M, Santos A R, Valle V E, Brunet M, Aguilar L, Álvarez L, Bautista M, Castañón C, Herrera L, Ruano E, Sinay J J, Sánchez E, Oviedo G I H, Obed F, Salgado J E, Vázquez J L, Baca M, Gutiérrez M, Centella C, Espinosa J, Martínez D, Olmedo B, Espinoza C E O, Núñez R, Haylock M, Benavides H, Mayorga R. Changes in precipitation and temperature extremes in Central America and northern South America, 1961-2003. Journal of Geophysical Research Atmospheres, 2005, 110(D23): D23107. DOI:10.1029/2005JD006119
[9]
封国林, 龚志强, 支蓉. 气候变化检测与诊断技术的若干新进展. 气象学报, 2008, 66(6): 892-905. DOI:10.3321/j.issn:0577-6619.2008.06.005
[10]
胡宜昌, 董文杰, 何勇. 21世纪初极端天气气候事件研究进展. 地球科学进展, 2007, 22(10): 1066-1075. DOI:10.3321/j.issn:1001-8166.2007.10.011
[11]
王晓利, 侯西勇. 1982-2014年中国沿海地区归一化植被指数(NDVI)变化及其对极端气候的响应. 地理研究, 2019, 38(4): 807-821.
[12]
杜佳梦, 包刚, 佟斯琴, 黄晓君, 温都日娜, 美丽, 包玉海. 1982-2015年蒙古国植被覆盖变化及其与气候变化和人类活动的关系. 草业学报, 2021, 30(2): 1-13.
[13]
罗琳玲, 易桂花, 张廷斌, 李景吉, 王海军, 别小娟, 王娜, 杜清江. 内蒙古植被长势变化及其对气候变化的响应. 测绘科学, 2021, 46(5): 126-134.
[14]
吴霞, 张永宏, 王长军, 樊丽琴, 李磊. 基于MODIS的宁夏植被覆盖时空变化特征分析. 江苏农业科学, 2021, 49(2): 204-208.
[15]
耿庆玲, 陈晓青, 赫晓慧, 田智慧. 中国不同植被类型NDVI对气候变化和人类活动的响应. 生态学报, 2022, 49(09): 1-12.
[16]
高滢, 孙虎, 徐崟尧, 张世芳. 陕西省植被覆盖时空变化及其对极端气候的响应. 生态学报, 2022, 42(03): 1-12.
[17]
张晋霞, 徐长春, 杨秋萍. 2001-2017年新疆NDVI变化及其对极端气候的响应. 水土保持通报, 2020, 40(5): 250-256, 275, 341.
[18]
张彬, 朱建军, 刘华民, 潘庆民. 极端降水和极端干旱事件对草原生态系统的影响. 植物生态学报, 2014, 38(9): 1008-1018.
[19]
穆少杰, 李建龙, 陈奕兆, 刚成诚, 周伟, 居为民. 2001-2010年内蒙古植被覆盖度时空变化特征. 地理学报, 2012, 67(9): 1255-1268.
[20]
龙鑫, 李静, 柳钦火. 植被指数合成算法综述. 遥感技术与应用, 2013, 28(6): 969-977.
[21]
李春兰. 蒙古高原多时空尺度极端气候变化特征及其影响研究[D]. 上海: 华东师范大学, 2019.
[22]
Sen P K. Estimates of the regression coefficient based on Kendall's tau. Journal of the American Statistical Association, 1968, 63(324): 1379-1389.
[23]
李明, 葛晨昊, 邓宇莹, 王贵文, 柴旭荣. 黄土高原气象干旱和农业干旱特征及其相互关系研究. 地理科学, 2020, 40(12): 2105-2114.
[24]
王劲峰, 徐成东. 地理探测器: 原理与展望. 地理学报, 2017, 72(1): 116-134.
[25]
孙继昌, 刘金平, 周国良. 1998年暴雨洪水. 水文, 1998, 18(S1): 97-103.
[26]
林婧婧, 张强. 中国气候态变化特征及其对气候变化分析的影响. 高原气象, 2015, 34(6): 1593-1600.
[27]
陈赛赛, 孙艳玲, 杨艳丽, 王中良. 三北防护林工程区植被景观格局变化分析. 干旱区资源与环境, 2015, 29(12): 85-90.
[28]
许旭, 李晓兵, 梁涵玮, 黄玲梅. 内蒙古温带草原区植被盖度变化及其与气象因子的关系. 生态学报, 2010, 30(14): 3733-3743.
[29]
赵慧颖, 田宝星, 宫丽娟, 曲辉辉, 季生太, 李秀芬, 张兴林. 近308年来大兴安岭北部森林植被气候生产潜力及其对气候变化的响应. 生态学报, 2017, 37(6): 1900-1911.
[30]
李峰, 李琪, 薛红喜, 吴东丽, 王云龙. 温度对克氏针茅草原生态系统生长季碳通量的影响. 农业环境科学学报, 2012, 31(7): 1453-1459.
[31]
雅茹, 丽娜, 银山, 包玉海. 1960-2015年内蒙古极端气候事件的时空变化特征. 水土保持研究, 2020, 27(3): 106-112.
[32]
Li C L, Wang J, Hu R C, Yin S, Bao Y H, Ayal D Y. Relationship between vegetation change and extreme climate indices on the Inner Mongolia Plateau, China, from 1982 to 2013. Ecological Indicators, 2018, 89: 101-109.
[33]
Alexander L V, Zhang X, Peterson T C, Caesar J, Gleason B, Klein Tank A M G, Haylock M, Collins D, Trewin B, Rahimzadeh F, Tagipour A, Rupa Kumar K, Revadekar J, Griffiths G, Vincent L, Stephenson D B, Burn J, Aguilar E, Brunet M, Taylor M, New M, Zhai P, Rusticucci M, Vazquez-Aguirre J L. Global observed changes in daily climate extremes of temperature and precipitation. Journal of Geophysical Research Atmospheres, 2006, 111(D5): D05109.
[34]
Li C L, Leal Filho W, Wang J, Fudjumdjum H, Fedoruk M, Hu R C, Yin S, Bao Y H, Yu S, Hunt J. An analysis of precipitation extremes in the Inner Mongolian plateau: spatial-temporal patterns, causes, and implications. Atmosphere, 2018, 9(8): 322.
[35]
罗琳玲. 内蒙古不同生态区植被长势对气候变化的响应[D]. 成都: 成都理工大学, 2020.
[36]
陈赛赛, 孙艳玲, 杨艳丽, 王中良. 三北防护林工程区植被景观格局变化分析. 干旱区资源与环境, 2015, 29(12): 85-90.
[37]
Li S D, Liu M C. The development process, current situation and prospects of the conversion of farmland to forests and grasses project in China. Journal of Resources and Ecology, 2022, 13(1): 120-128.