自然资源遥感, 2025, 37(1): 221-231 doi: 10.6046/zrzyyg.2023226

技术应用

2000—2020年山西省植被生长状态对气象干旱的响应

赵富,, 王莉,, 马禹昂, 蒋瑞霞

河南理工大学测绘与国土信息工程学院,焦作 454003

Responses of vegetation growth conditions to meteorological drought in Shanxi Province from 2000 to 2020

ZHAO Fu,, WANG Li,, MA Yuang, JIANG Ruixia

School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454003, China

通讯作者: 王 莉(1981-),女,博士,副教授,主要从事地理信息系统研究。Email:wangli29@hpu.edu.cn。

责任编辑: 李瑜

收稿日期: 2023-07-24   修回日期: 2023-10-10  

基金资助: 河南省科技攻关项目“基于GIS的大学生体质健康评价与预警技术研究”(192102310025)
河南理工大学国家社科基金预研项目“我国城乡青少年体质健康时空差异及体制提升研究”(CSKY2023-04)

Received: 2023-07-24   Revised: 2023-10-10  

作者简介 About authors

赵 富(2000-),男,硕士研究生,主要从事气象学与地理信息技术及应用研究。Email: zhaofu_2000@163.com。

摘要

基于2000—2020年标准化降水指数(standardized precipitation index, SPI)与植被状态指数(vegetation condition index, VCI),采用变分模态分解(variational mode decomposition, VMD)、Mann-Kendall趋势分析、Pearson相关系数等方法,分析了山西省不同植被类型下气象干旱变化趋势,量化了植被生长状态对气象干旱的响应时间。结果表明: ①20世纪以来,山西省整体气象干旱逐步缓解,但季尺度上呈现由春至冬旱情加重区不断扩大的态势; ②不同植被类型下气象干旱均得到缓解,且随着时间尺度增长,态势越发显著,而季尺度上由春至冬缓解态势渐弱,加重态势渐强; ③植被生长状态受气象干旱影响显著,年尺度上主要以正相关为主,季尺度上高相关性区域由春至冬逐渐收缩至西北部和东北部,而响应时间呈现出春冬长于秋夏的特征; ④不同植被类型下植被生长状态对气象干旱的响应均在夏季最为迅速,且耕地是对气象干旱最敏感的地类。

关键词： SPI; VCI; 气象干旱; 响应关系; 山西省

Abstract

Based on the standardized precipitation index (SPI) and vegetation condition index (VCI) from 2000 to 2020, this study analyzed the trends in meteorological drought across different vegetation types in Shanxi Province using methods such as variational mode decomposition (VMD), Mann-Kendall trend analysis, and Pearson correlation coefficient. Accordingly, this study quantified the response time of vegetation growth conditions to meteorological drought. The results indicate that from the beginning of the 2000s, the overall meteorological drought in Shanxi Province has gradually eased. However, on a seasonal scale, areas with increasingly aggravated drought continuously expand from spring to winter. Meteorological drought has alleviated across various vegetation types, with the alleviation becoming increasingly significant with an increase in the time scale. In contrast, on a seasonal scale, the drought relief gradually weakens from spring to winter, during which drought aggravation progressively strengthens. Vegetation growth conditions are significantly influenced by meteorological drought. On the annual scale, there is a predominantly positive correlation between both. On the seasonal scale, areas with a strong correlation between both gradually contract from spring to winter, when such areas are dominated by the northwestern and northeastern areas of the province. Additionally, the response time of vegetation to drought is longer in spring and winter compared to autumn and summer. Across different vegetation types, the responses of vegetation growth conditions to meteorological drought prove the most rapid during the summer, and cultivated lands are identified as the most sensitive land type to meteorological drought.

Keywords： SPI; VCI; meteorological drought; response relationship; Shanxi Province

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

本文引用格式

赵富, 王莉, 马禹昂, 蒋瑞霞. 2000—2020年山西省植被生长状态对气象干旱的响应[J]. 自然资源遥感, 2025, 37(1): 221-231 doi:10.6046/zrzyyg.2023226

ZHAO Fu, WANG Li, MA Yuang, JIANG Ruixia. Responses of vegetation growth conditions to meteorological drought in Shanxi Province from 2000 to 2020[J]. Remote Sensing for Land & Resources, 2025, 37(1): 221-231 doi:10.6046/zrzyyg.2023226

0 引言

伴随着全球变暖和气候异常,世界上越来越多的地区长期降水匮乏,促使干旱成为频发性自然灾害,虽然其演变过程缓慢,但影响范围广、危害程度深,对生态环境、农业生产及社会经济都造成了严重威胁[1-2]。根据现有研究,干旱常分为气象干旱、农业干旱、水文干旱和社会经济干旱[3],其中气象干旱常发生在缺乏降水或降水量长期显著低于平均水平的地区,通常是其他类型干旱发生的主要驱动力[4],也是造成植被生长状态恶化的主要原因。

气象干旱常用基于气象站点的干旱监测指数表征,如标准化降水指数(standardized preciptation index,SPI)、帕默尔干旱指数(Palmer drought severity index, PDSI)、标准化降水蒸散指数(standardized precipitation evapotranspiration index, SPEI)、综合气象干旱指数(composite index of meteorological drought, CI)[5]等,其中SPI理论完善、易于计算且特征明显,已被广泛应用于评估不同时间尺度的气象干旱[4,6]。而植被生长状态常采用基于遥感的干旱监测指数表征,如归一化植被指数(normalized difference vegetation index, NDVI)、植被状态指数(vegetation condition index,VCI)、温度状况指数(temperature condition index, TCI)等,其中NDVI是衡量植被覆盖率、反映植被长势的重要指标,但对于不同植被类型的解释能力有限,易受云、阴影和大气干扰[7-8],故在此基础上提出VCI用于植被健康状况评估和水分监测,目前已有诸多学者研究证明了VCI在黄土高原及西北地区存在良好的适用性。李家誉等[9]以VCI作为植被状况指标,得到了黄土高原植被状况时空变化特征并分析了其与气象干旱的相关性; 李明等[10]采用SPI和VCI分别表征气象干旱和农业干旱,进而探究二者的时滞关系,证明了VCI能较好反映黄土高原地区的干旱情况; 张华等[11]通过对NDVI和VCI所表征的植被生长状况与SPEI的相关性进行对比分析,证明了VCI相较于NDVI适用性更好。

对不同指数之间的响应关系的研究常以2种思路展开。其一,以统计面雨量的方式,按以点代面的形式,分析2个指数时间序列的相关性,进而量化响应时间、拟合响应关系。刘永佳等[6]以时序SPI和标准化径流指数(standardized runoff index, SRI)为基础,研究无定河、窟野河和沁河流域水文干旱对气象干旱的响应时间,得到无定河、窟野河和沁河流域气象干旱向水文干旱的年传播时间为5~7个月; Dai等[12]以SPI和标准化土壤水分指数(standardized soil moisture index, SSMI)表征气象干旱和农业干旱,通过选择不同时间尺度的SPI和1个月时间尺度的SSMI之间的最大条件概率值,确定从气象干旱到农业干旱的传播时间,结果表明夏季和秋季干旱传播时间短于春季和冬季,说明夏季和秋季干旱传播速度快于春季和冬季。其二,通过逐项元计算2个指数影像集的相关性,直观分析响应时间在空间上存在的差异。靖娟利等[7]采用趋势分析和相关性分析,以SPEI和NDVI影像集为基础,研究植被变化及其对气象干旱的响应特征,结果表明滇黔桂地区NDVI对短时间尺度SPEI响应比较敏感; 吕振涛等[13]借助NDVI和SPEI,研究蒙古国植被对气象干旱的敏感性,结果表明大部分区域植被生长受干旱影响较强,且相关性显著区域主要在中部地区,各植被指数与干旱指数呈较明显的正相关关系,且草地是最敏感的地类。本研究将2者融合,在完成空间相关性分析的基础上,得到可用来量化响应关系的平均响应时间,同时完成不同植被类型下植被生长状态对气象干旱响应时间的比较。

山西省位于黄土高原东翼,远离海洋,再加上东西有太行、吕梁两山耸立,阻隔水汽,故降水稀少,干旱频发。全省林草丰富,且以夏、秋为主要粮食产季,产量易受干旱影响。多数学者以黄土高原为研究区,聚焦于整体的干旱时空特征分析或干旱响应关系研究[9-10,14-15],本研究选取黄土高原中属典型半干旱半湿润区的山西省为研究区,借助Google Earth Engine(GEE)平台,基于气象站点数据和遥感监测数据,将变分模态分解(variational model decomposition,VMD)与Mann-Kendall相结合,分析山西省2000—2020年不同植被类型下气象干旱变化趋势,同时采用点、面结合的思路,研究其植被生长状态对气象干旱的响应机制,量化响应时间,为山西省应对气象变化、实现生态建设提供理论支撑。

1 研究区概况及数据源

1.1 研究区概况

山西省地形错综复杂、河谷纵横,西部是以吕梁山为主干的黄土高原,东部是以太行山为主脉形成的块状山地,汾河、桑干河贯穿山谷之中,盆地、平原沿河由东北至西南分散,形成 “两山夹一川”的地形格局,属我国生态环境脆弱地带。全省的植被类型多样,除主要城市建设区以外,大部分区域被林地、草地和耕地所覆盖,其中草地主要分布在地形起伏较小的丘陵和山地,占总面积的61.32%(图1)。同时,受温带大陆性季风气候和地理因素影响,山西省冬夏气温悬殊,省内降水分布不均,降水季节性显著。全省各地年降水量介于358~621 mm之间,且主要集中在夏季6—8月,约占全年降水量的70%,其余季节降雨稀少,总体水资源贫乏,旱情严重[16-17]。

图1

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

图1   山西省地形、植被类型及气象站点分布

Fig.1   Distribution of terrain, vegetation types and meteorological stations in Shanxi Province


1.2 数据源及其预处理

1.2.1 气象数据

本研究所采用的山西省2000—2020年逐日降水数据来自国家气象信息中心整理的国家级地面气象站基本气象要素日值数据集(V3.0)(http://data.cma.cn/),将监测时序不完整的气象站点剔除,选取24个气象站点反映山西省降水情况(图1(a))。

1.2.2 MODIS数据

本研究所采用的MODIS数据来自GEE平台(https://code.earthengine.google.com/),分别获取了植被指数产品MOD13A1(空间分辨率为500 m,时间分辨率为16 d)和土地覆盖类型产品MCD12Q1(空间分辨率为500 m,时间分辨率为1 a)。对于MOD13A1产品,为减少云雾和霾的影响,采用最大值合成法(maximum value composite, MVC)将NDVI原始序列合成为月尺度序列[18]。对于MCD12Q1产品,为了避免地类变化造成的影响,采用众数合成法进行多年合成,再将得到的土地覆盖类型数据基于国际地圈生物圈计划(International Geosphere Biosphere Programme, IGBP)中全球植被分类方案重新划分为4类: 耕地(包括天然植被的作物和镶嵌农田)、草地(包括乔木草原、稀疏草原和草原)、林地(包括树栖针叶林、树栖阔叶林、落叶针叶林、落叶阔叶林、混交林、封闭灌木林和稀疏灌木林)、无植被区(包括城市和建筑区、冰、雪、裸露土地或稀疏植被覆盖、湿地和水域)[19-20]。

1.2.3 其他数据

本研究所采用的全国行政区划数据来自中国科学院资源环境科学与数据中心(https://www.resdc.cn/), 山西省30 m分辨率DEM数据来自地理空间数据云(http://www.gscloud.cn/)。

2 研究方法

2.1 SPI指数

SPI是由McKee等[21]于1993年提出的气象干旱指数,其计算所需参数少且易于获取,2009年被世界气象组织推荐为主要的气象干旱指数,用于旱情监测与跟踪[22]。一般而言,SPI计算步骤大致为: 首先,计算某一时期的累积降水量,构建时间序列; 其次,使用概率分布估计时间序列的累积概率; 最后,将累积概率转换为标准正态分布[23-24]。具体的计算公式如下[25-26]:

SPI=Gt-(u3t+u2)t+u1[(l3t+l2)t+l1]t+1，
t=ln1Y(x)2，

式中: x为降水量; G为正负系数,当Y(x)>0.5时,G=1,当Y(x)≤0.5时,G=-1; u1,u2,u3,l1,l2,l3为常数,分别为: u1=2.515 517,u2=0.802 853,u3=0.010 328,l1=1.432 788,l2=0.189 269,l3=0.001 308; Y(x)2为与Γ函数相关的降水量分布概率,由Γ函数概率密度积分公式求得,即

Y(x)=2βγΓ(γ)∫0xxγ-1e-x/βdx，
Γ(γ)=∫0∞xγ-1e-xdx，

式中γ和β分别为Γ分布函数的形状和尺度参数。

基于气象站点降水数据计算得到的SPI具有多时间尺度特性,记为SPI1,SPI2,…,SPI24。其中SPI3可用来分析气象干旱的季节特性[27],SPI12则可表征气象干旱的年际变化[9],适用于本研究中植被生长状态对气象干旱的响应与时滞关系探究。

2.2 VCI指数

VCI是通过对某年某月的NDVI与研究时段内所有年份同月的NDVI最大和最小值归一化得到的植被指数,单幅影像可分析地区内不同植被的相对生长状况,多幅影像则可分析植被的生长变化趋势。具体的计算公式如下[11,28-29]:

VCIi=NDVIi-NDVIminNDVImax-NDVImin×100，

式中: VCIi为某年第i月VCI; NDVIi为某年第i月NDVI; NDVImax和NDVImin分别为研究序列中第i月的NDVI最大值和最小值。

2.3 趋势分析

Mann-Kendall趋势分析是一种非参数检验方法,具有不受数据概率分布形式影响,且不被异常值干扰的优点,常用于气象和水文时间序列的趋势检验[4,24,26]。本研究取|z|>1.64,1.96和2.58分别代表原假设通过90%,95%和99%的显著性检验。同时,基于最小二乘法的线性拟合,建立线性回归模型,计算斜率,表征时间序列趋势。

2.4 VMD分解

VMD分解是一种自适应、完全非递归的模态变分和信号处理方法[30]。其可根据实际情况灵活确定分解模态数,随后匹配每种模态的最佳中心频率和有限带宽,得到原始信号的最佳分解方式,具有更坚实的数学基础[31]。VMD对于构建和求解变分问题,具体过程为: 首先构造带约束条件的变分问题,约束条件为分解后的分量和与原信号相等,表达式为[30⇓-32]:

min{uk},{ωk}∑k‖∂tδ(t)+jπt·uk(t)e-jωkt‖22，
s.t.∑k=1Kuk=f(t)，

式中: K为分解模态数; {uk}和{ωk}分别为k个模态分量和中心频率; e-jωkt为混合算子; ‖‖22为频率梯度2范数的平方; δ(t)为狄拉克分布; f(t)为输入信号。

接着引入拉格朗日乘法算子与二次惩罚因子将约束变分问题转换为非约束变分问题,最后采用交替方向乘子法即可解决上述问题。表达式为:

L({uk},{ωk},λ(t))=α∑k‖∂tδ(t)+jπt·ukte-jωkt‖22+‖f(t)-∑kukt‖22+λ(t),f(t)-∑kukt，

式中: α为二次惩罚因子,在高斯噪声存在的情况下保证信号重构精度; λ(t)为拉格朗日乘法算子,使约束条件保持严格性。

在实现VMD分解时,需提前设置惩罚因子、噪声容限、收敛容错和分解模态数,进行信号初始化和常量判断。其中,惩罚因子α影响分解的精度,太小导致分解模态精度不高,但太大又面临着分解时间增多,甚至可能造成死循环。对于带有明显频率成分的信号,选用较小的惩罚因子能更好地捕捉信号的频率特征,而对于嘈杂或频率特征不明显的信号,可用较大的惩罚因子来限制模态的带宽,以提高噪声抑制能力; 噪声容限tau主要控制对噪声的容忍度,过大的噪声容限可能会导致信号损失,对于高噪声信号,可以选用较大的噪声容限以容忍噪声成分,对于低噪声信号,可以选用较小的噪声容限以更严格地限制噪声; 收敛容错tol控制误差大小,较小的收敛容错会提高分解精度,但也会导致耗费更长的计算时间。此三者会对分解结果的好坏产生间接影响。此外,信号初始化init设置为1进行均匀初始化; 常量判断DC根据是否需要考虑直流成分而定,如需分析信号的直流成分,将其设置为1,否则为0。根据先验知识并结合多次实验结果,本研究取α=2 000,tau=0.024 4,tol=1E-6,init=1,DC=0。

除以上参数外,由于分解模态数K直接决定分解结果是否正确,故研究确定合适的分解模态数常作为相关研究的重点[32]。而本研究采用中心频率法确定分解模态数,大致过程为:

1)将K=2,3,…,n依次代入VMD,计算各分解结果的高频模态的中心频率,记为ω2,ω3,…,ωn。

2)若其中|ωt-ωt+1|≤0.001⋂|ωt+1-ωt+2|≤0.001,表示中心频率出现平缓趋势,K=t即为最优解。

2.5 相关性分析

Pearson相关系数常用于研究不同干旱指数间的相关性和一致性。其计算公式为[33]:

R=∑i=121∑j=112[(xij-x-)(yij-y-)]∑i=121∑j=112(xij-x-)2∑i=121∑j=112(yij-y-)2，

式中: R为相关系数; xij,yij分别为第i年第j月的VCI值和SPI值。本研究基于GEE平台,利用最大值合成的月NDVI求得VCI,再采用均值法合成年尺度VCI,同时选取插值而来的SPI12,逐像元分析年尺度的相关性。对于月和季尺度的相关性分析,逐像元计算各月多年VCI和相应的1~24个月时间尺度SPI之间的相关系数,得到288幅相关系数图,通过最大值合成各月、季(春季为3—5月,夏季为6—8月,秋季为9—11月,冬季为12—次年2月)最大相关系数图,其对应的SPI时间尺度即为滞后时间[6,7,13]。

3 结果与分析

3.1 气象干旱变化趋势

3.1.1 年与月尺度气象干旱变化趋势

通过经验贝叶斯克里金法对各站点SPI的Mann-Kendall趋势检验结果和线性拟合斜率进行插值,分析山西省年、月尺度SPI空间变化趋势(图2)。其中,年SPI以临汾、榆社和平定为界,呈现出北升南降的总趋势,东北部太行山、五台山附近以8%/a的速率上升,东南部安泽、长治等地以-0.98%/a的速率下降; SPI1全省总体呈上升趋势,以五台山为中心上升最为明显,速率为0.1%/a,南部以隰县、介休和榆社为界,南部上升趋势不明显,速率趋近于0.01%/a; SPI3趋势变化特征与SPI1类似,区别在于五台山附近上升速率为0.3%/a; SPI12大部分地区呈上升趋势,五台山附近以0.8%/a的速率上升,东南部以安泽、襄垣两地为中心辐射阳城,形成下降三角区。

图2

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

图2   年与月尺度SPI变化趋势

Fig.2   Variation trend of SPI at the interannual and monthly scales


3.1.2 季尺度气象干旱变化趋势

分析山西省季尺度SPI变化趋势可知(图3),春季SPI全省98%的地区呈上升趋势,除朔州、大同交界处和朔州西北部外,其余地区旱情均得到缓解; 夏季SPI全省北部呈上升趋势,南部以安泽、襄垣、阳城和长治为中心形成闭环,夏旱呈现加重趋势; 秋季SPI以兴县、平定为界,以北呈上升趋势,以南呈下降趋势; 冬季全省97%的地区旱情加重,只有灵丘和朔州2地SPI呈上升趋势。综合分析,山西省季尺度SPI呈上升趋势的区域,随着季节更替,一直处于由南至北方向上的动态移动,同时该区域的面积也在不断减小。

图3

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

图3   季尺度SPI变化趋势

Fig.3   Variation trend of SPI at the seasonal scale


3.1.3 不同植被类型下气象干旱变化趋势

为了探究不同植被类型下SPI的变化趋势,采用均值法合成面雨量,计算出不同植被类型下多时间尺度SPI,再分析趋势(表1)。以草地春季SPI序列为例,采用VMD分解为9个模态分量(图4),记为IMF1,…,IMF9,其中IMF1为趋势分量,反映了原始序列的总体趋势,同时结合趋势检验结果中草地春季SPI斜率为0.049,呈上升趋势,与IMF1一致。

表1   不同植被类型下年、月与季尺度SPI变化趋势

Tab.1  Variation trend of SPI at different scales under different vegetation types

土地类型类别年SPI1SPI3SPI12春季夏季秋季冬季
耕地z①0.6040.6091.693*2.420**1.4190.635-0.544-0.974
s0.0440.0010.0020.0030.0460.033-0.025-0.031
草地z2.477**1.0752.818***6.555***1.3591.1780.0300.616
s0.0820.0010.0030.0060.0490.0480.002-0.021
林地z0.9670.4071.4392.461**1.5710.574-0.574-0.974
s0.0340.0000.0010.0020.0500.033-0.025-0.033
无植被区z1.2980.7532.141**3.661***1.4801.238-0.513-0.876
s0.0510.0010.0020.0030.0400.050-0.027-0.035

① z表示显著性检验结果,s表示线性拟合斜率,*表示通过90%显著性检验,**表示通过95%显著性检验,***表示通过99%显著性检验。

新窗口打开| 下载CSV


图4-1

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

图4-1   草地春季SPI的VMD分解

Fig.4-1   VMD decomposition of grassland SPI in spring


图4-2

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

图4-2   草地春季SPI的VMD分解

Fig.4-2    VMD decomposition of grassland SPI in spring


3.2 植被生长状态对气象干旱的响应

3.2.1 年尺度植被生长状态与气象干旱的相关性

分析年尺度VCI与SPI12的相关系数及显著性检验结果可看出(图5),二者主要以正相关为主,占全省面积的94%,负相关区域主要集中在安泽和阳城附近,极显著正相关和显著正相关主要分布在以隰县、介休、太谷和平定为界的山西北部,占正相关区域面积的25%,以南主要为不显著正相关和不显著负相关。

图5

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

图5   年尺度VCI与SPI12的相关系数及显著性检验空间分布

Fig.5   Spatial distribution of correlation coefficient and significance test between interannual VCI and SPI12


3.2.2 季尺度植被生长状态对气象干旱的响应

通过对每月VCI和对应月份多时间尺度SPI求解相关系数,反映植被状态对不同降水累积时间下的气象干旱响应过程。将每月的相关系数图按季节划分合成最大相关系数图(图6),并得到相应的响应时间(图7)。

图6

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

图6   季尺度VCI与不同时间尺度SPI的最大相关系数空间分布

Fig.6   Spatial distribution of maximum correlation coefficient between seasonal VCI and SPI at different time scales


图7

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

图7   季尺度VCI对SPI的响应时间

Fig.7   Response time of seasonal VCI to SPI


总体而言,各季VCI与SPI均以正相关为主,并且R>0.55的区域,随着季节更替,由春至冬收缩至西北和东北部。其中,春季低相关性区域(R<0.45)主要集中在太原及山西中南部的无植被区,占全省总面积的4%,高相关性区域(R>0.75)占全省总面积的26%,在西北部兴县、五寨,及中部隰县、介休、榆社集中分布; 夏季低相关性区域范围扩大至安泽、阳城,高相关性区域由西北部转移至东北部的右玉、大同,由中部向北转移至离石; 秋季以隰县、榆社和平定为界,北部主要为R>0.55的区域,其中高相关性区域占11%,低相关性区域占R<0.45区域的56%; 冬季R>0.55的区域收缩至山西西北吕梁山和东北太行山地区,低相关性区域占全省面积的40%。

VCI对SPI的响应时间存在空间差异,且随季节变化明显。春季响应时间以长期为主,响应时间小于12个月的区域主要集中在东北部的天镇、灵丘和五台山地区,仅占全省面积的18%,其余地区由北至南响应时间从13~24个月逐渐延长; 夏季响应时间仍以长期为主,响应时间大于12个月的区域占61%,响应时间小于12个月的区域主要沿兴县、五寨、朔州、右玉、大同和天镇一线分布,南部盐湖、垣曲、阳城和长治一线响应时间相较春季也明显缩短; 秋季响应时间呈长短期离散分布,其中响应时间为7个月、18个月和11个月的区域分别占32%,11%和9%; 冬季响应时间主要为9~18个月,响应时间超18个月的区域主要分布在太原及南部平原等地,占全省面积的17%。

3.2.3 月尺度植被生长状态对气象干旱的响应

分析全省各月VCI与1~24个月时间尺度SPI的相关系数均值(图8),可以看出整个山西省VCI对SPI的响应时间随季节变化明显,R>0.4主要出现在4—8月,其余各月相关系数随着SPI时间尺度增加也不断增大。从中选取各月响应时间均值作为对应季节的响应时间(如春季响应时间,计算3,4,5月最大相关系数对应的SPI时间尺度的均值),春季为13.7个月,夏季为10个月,秋季为14.7个月,冬季为22.7个月,说明夏季气象干旱对植被生长状态的影响反应最快。

图8

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

图8   各月VCI与1~24个月时间尺度SPI的相关系数均值

Fig.8   The average value of the correlation coefficient between VCI and SPI with a time scale of 1 to 24 months


3.2.4 不同植被类型下植被生长状态对气象干旱的响应

为进一步细化VCI对SPI响应的空间差异,统计分析了不同植被类型下生长季VCI与SPI的最大相关系数及对应时间尺度(表2)。结果表明,耕地、草地、林地和无植被区的最大相关系数分别出现在7月、7月、5月和6月,其对应的响应时间分别为6个月、6个月、3个月和5个月。而耕地、草地、林地和无植被区最短响应时间分别出现在5月、6月、5月和6月。4种地类在春夏两季的响应时间相较秋冬要短,且在冬季响应时间都最长。同时,耕地在各月的响应时间都要短于草地和林地,说明耕地对气象干旱更敏感。

表2   不同植被类型下生长季最大相关系数及对应时间尺度

Tab.2  Maximum correlation coefficients and corresponding time scales of growing seasons under different vegetation types

土地
类型		4月5月6月7月8月9月10月
R响应时
间/月		R响应时
间/月		R响应时
间/月		R响应时
间/月		R响应时
间/月		R响应时
间/月		R响应时
间/月
耕地0.3880.4120.4750.4960.4470.37160.249
草地0.49150.47160.4650.5060.49190.46160.439
林地0.33240.4930.3750.25180.2870.24160.3712
无植被区0.3270.2980.4150.3860.27100.21160.1610

新窗口打开| 下载CSV


4 讨论与结论

4.1 讨论

4.1.1 气象干旱变化特征

山西省“两山夹一川”的地理格局,直接造成山区多雨、平原少雨的降雨分布,再加之受东亚季风减弱影响[34],省内东南部降水异常减少,北部无明显变化,气象干旱年际变化也呈现出南重北轻的趋势。这与李明等[10]的研究结果一致,一致性也体现在春、夏和秋3季,但在冬季存在差异,可能是由于研究时间跨度不同造成的。此外,本研究用于表征气象干旱的SPI主要受降水影响,在晏利斌[35]、侯青青等[36]和安彬等[37]的研究中山西省降水变化趋势与本研究的结果具有空间一致性也可得出这一结论。

同时,土壤含水量也能有效地反映旱情,根据薛少博等[38]研究,山西省在2002—2020年间水储量明显减少,且由西北至东南减少趋势更为明显,同时,在东南部0~200 cm土层土壤水分均流失严重,这与气象干旱年际变化相吻合。

4.1.2 植被生长状态对气象干旱的响应

基于VCI与SPI探究植被生长状态对气象干旱的响应,各月相关系数空间分布以正相关为主,这与李明等[10]的研究结果相似,但相关程度不同,可能是由于研究区域大小不同,气候变化存在空间差异性。同时,考虑到SPEI相较于SPI将蒸散发纳入对干旱的响应,VCI则是在NDVI基础上求得的植被指数,故二者相关性与SPI-VCI存在相似性,这在史尚渝等[39]的研究中得到验证。

此外,探究植被生长状态对气象干旱的响应,应综合考虑多种因素影响。其一,相关指数选取。SPI作为仅依靠降水数据得到的干旱指数,在气象干旱表征中存在一定局限性,同时在探究响应机制时也会存在不足,所以对于气象干旱指数应选取综合了气温、地形、海拔、降水和蒸散发等因子的干旱指数来表征。VCI主要通过植被覆盖度和植被水分来量化植被生长状态,而山西省以夏、秋为主要粮食产季,冬季农作物种植较少且常有积雪覆盖,春季农作物长势不明显,故农作物面积与实际耕地面积常常不符,在借助VCI反映农作物收割前后的耕地旱情时存在偏差,可能出现VCI值比实际偏小的情况。其二,影响因素的考量。植被生长状态不仅受气象因素影响,同时还与土壤状况、土地类型、人类活动有密切关联。自1999年以来,山西省积极响应国家退耕还林还草政策,植被覆盖率大幅增长[40],大面积区域植被得到改善。但政策施行初期伴随着植被蒸腾作用增强和气候异常,植被返青所需的土壤水含量不足,植被生长受限[41]。植被的生长发育是受到多种因素共同作用的结果。综上所述,选取更合理的评价指数,综合考虑各种影响因子,将是未来气象干旱评估重要的发展方向。

4.1.3 对于不同植被类型的考虑

针对不同植被类型研究气象干旱变化趋势,分析植被生长状态与气象干旱的响应关系,为了避免长时间跨度地类变化造成的影响,本研究采用众数合成法对土地覆盖类型进行多年合成,对于此问题其他学者也有不同的解决办法,如尉毓姣等[42]将研究时段划分为4期土地覆盖数据,分时段提取不同类型植被的NDVI和SPEI; 吴林霖等[43]使用2005年和2010年数据分别提取2000—2008年、2009—2013年间的典型植被类型; 李家誉等[9]统计了黄土高原地区1990年和2015年耕地、林地和草地3种土地利用类型未发生改变的格点。不同的方法有不同的可行性,但如何更合理地将随时间变化的各种因素纳入对气象干旱研究的考虑范围,得出更科学的结论是目前亟待解决的问题。

4.2 结论

研究气象干旱的变化趋势,分析气象干旱与植被生长状态的关系,有利于更好地了解干旱的演变规律及其对社会经济造成的影响。本研究的主要结论如下:

1)2000—2020年山西省年与12个月尺度气象干旱均呈现北部缓解、东南部加重的态势,而在1个月和3个月尺度则整体缓解,不同尺度下的相同之处在于缓解态势最明显的区域均在五台山附近; 山西省气象干旱存在明显的季节性差异,呈现由春至冬旱情加重区不断扩大的态势。不同植被类型(耕地、草地、林地和无植被区)下气象干旱均得到缓解,且随着时间尺度增长,态势越发显著,而季尺度上由春至冬缓解态势渐弱,加重态势渐强。

2)山西省年尺度植被生长状态与气象干旱主要以正相关为主,且极显著与显著正相关区域向五台山附近集中,表明植被受气象干旱影响显著,且阻隔季风与湿气的山地植被对降水更为敏感; 季尺度上高相关性区域由春至冬逐渐收缩至西北部和东北部,表明春季植被生长状态受降水限制最强,但在季节更替中,该限制逐渐减弱。而植被生长状态对气象干旱的响应时间存在空间差异,且随季节变化,以长期响应为主的区域不断缩减,短期响应的区域逐步从东北部向西北部再向全省蔓延。同时经过量化分析各月VCI与1—24个月时间尺度SPI的相关系数均值,发现各月响应时间存在明显季节性差异,其中夏季气象干旱对植被生长状态的影响反应最快。对于不同植被类型,最大相关系数集中出现在5—7月,也说明夏季植被生长状态对气象干旱响应更快,同时耕地在各月的响应时间都要短于草地和林地,说明耕地对气象干旱更敏感。

参考文献

Varela V, Vlachogiannis D, Sfetsos A, et al.

Projection of forest fire danger due to climate change in the French Mediterranean region

[J]. Sustainability, 2019, 11(16):4284.

[本文引用: 1]

Zhou Z Q, Shi H Y, Fu Q, et al.

Assessing spatiotemporal characteristics of drought and its effects on climate-induced yield of maize in Northeast China

[J]. Journal of Hydrology, 2020,588:125097.

[本文引用: 1]

Heim R R Jr.

A review of twentieth-century drought indices used in the United States

[J]. Bulletin of the American Meteorological Society, 2002, 83(8):1149-1166.

[本文引用: 1]

Li J, Wu C, Xia C A, et al.

Assessing the responses of hydrological drought to meteorological drought in the Huai River basin,China

[J]. Theoretical and Applied Climatology, 2021, 144(3):1043-1057.

[本文引用: 3]

吴燕锋, 巴特尔·巴克, 李维, 等.

基于综合气象干旱指数的1961—2012年阿勒泰地区干旱时空演变特征

[J]. 应用生态学报, 2015, 26(2):512-520.

[本文引用: 1]

基于阿勒泰地区7个气象站1961&mdash;2012年的逐日气温、降水等气象资料，运用趋势分析法、M-K突变检验法、小波分析法并结合ArcGIS软件，探究了阿勒泰地区干旱的时空演变特征.结果表明： 综合气象干旱指数能较好地反映阿勒泰地区干旱状况.研究期间，阿勒泰地区虽然年和四季的干旱频次和不同等级干旱覆盖范围都呈减少趋势，但是旱情仍较为严重，年和季节的干旱发生频繁，且多年存在重旱和特旱发生在整个阿勒泰地区；秋季和冬季的干旱频次分别在1997和1983年发生减少突变，夏季的干旱频次首先在1962年发生增加突变，随后在1991年发生减少突变，年均和春季的干旱频次无突变发生；年和四季的干旱频次都存在明显的周期性.干旱频次和不同干旱等级在空间上的分布表明，东部清河县旱情较为严重，中部阿勒泰市、富蕴县、布尔津县、福海县旱情次之，西部地区的哈巴河和吉木乃县的旱情较轻.&nbsp;

Wu Y F, Bake B, Li W, et al.

Spatio-temporal variation of drought condition during 1961 to 2012 based on composite index of meteorological drought in Altay region,China

[J]. Chinese Journal of Applied Ecology,2015, 26(2):512-520.

[本文引用: 1]

刘永佳, 黄生志, 方伟, 等.

不同季节气象干旱向水文干旱的传播及其动态变化

[J]. 水利学报, 2021, 52(1):93-102.

[本文引用: 3]

Liu Y J, Huang S Z, Fang W, et al.

Propagation and dynamic change of meteorological drought to hydrological drought in different seasons

[J]. Journal of Hydraulic Engineering, 2021, 52(1):93-102.

[本文引用: 3]

靖娟利, 王永锋, 和彩霞.

滇黔桂地区NDVI变化及其对SPEI的响应特征

[J]. 长江流域资源与环境, 2022, 31(8):1763-1775.

[本文引用: 3]

Jing J L, Wang Y F, He C X.

Spatial and temporal variations of NDVI and its response to SPEI in Yunnan-Guizhou-Guangxi region

[J]. Resources and Environment in the Yangtze Basin, 2022, 31(8):1763-1775.

[本文引用: 3]

张更喜, 粟晓玲, 郝丽娜, 等.

基于NDVI和scPDSI研究1982—2015年中国植被对干旱的响应

[J]. 农业工程学报, 2019, 35(20):145-151.

[本文引用: 1]

Zhang G X, Su X L, Hao L N, et al.

Response of vegetation to drought based on NDVI and scPDSI data sets from 1982 to 2015 across China

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(20):145-151.

[本文引用: 1]

李家誉, 佘敦先, 张利平, 等.

黄土高原植被变化对气象干旱多尺度响应特征与机制

[J]. 水土保持学报, 2022, 36(6):280-289.

[本文引用: 4]

Li J Y, She D X, Zhang L P, et al.

Multi-scale response characteristics and mechanism of vegetation to meteorological drought on the Loess Plateau

[J]. Journal of Soil and Water Conservation, 2022, 36(6):280-289.

[本文引用: 4]

李明, 葛晨昊, 邓宇莹, 等.

黄土高原气象干旱和农业干旱特征及其相互关系研究

[J]. 地理科学, 2020, 40(12):2105-2114.

DOI:10.13249/j.cnki.sgs.2020.12.017      [本文引用: 4]

选择标准化降水指数（SPI）和植被状态指数（VCI）分别作为评价黄土高原气象干旱和农业干旱的指标，使用干旱频率和Sen斜率分析了黄土高原地区干旱的分布特征与变化趋势，并探讨了气象干旱与农业干旱的相关性。结果表明：① 黄土高原西部干旱频率总体上高于东部。气象干旱和农业干旱变化趋势在空间上表现有所不同，黄土高原西部、北部气象干旱呈不显著减缓趋势，东部和南部呈不显著加重趋势，但绝大部分地区的农业干旱呈减缓趋势，尤其是400 mm等降水量一线两侧区域。② 季节上，黄土高原夏季和秋季气象干旱频率较高，春季和冬季气象干旱频率相对较低。黄土高原农业干旱频率春季最高，夏季其次，VCI对农业干旱实时监测的适用性更强。③ 不同季节，农业干旱滞后气象干旱的时间长短不同，冬季滞后约2个月，春季滞后约1个月，夏季和秋季滞后少于1个月。黄土高原一熟制种植区的SPI-12值与VCI值具有较好的正相关性。研究结果可以为黄土高原的干旱监测和预警、干旱区划以及干旱灾害风险评估提供科学依据。

Li M, Ge C H, Deng Y Y, et al.

Meteorological and agricultural drought characteristics and their relationship across the Loess Plateau

[J]. Scientia Geographica Sinica, 2020, 40(12):2105-2114.

DOI:10.13249/j.cnki.sgs.2020.12.017      [本文引用: 4]

In recent years, the combined effects of climate change and soil and water conservation measures have led to more complex drought conditions in the Loess Plateau, which has caused huge agricultural and socioeconomic losses. In this paper, standardized precipitation index (SPI) and vegetation condition index (VCI) are selected as indicators to evaluate the meteorological and agricultural droughts on the Loess Plateau, respectively. The distribution characteristics and trends of droughts in the Loess Plateau are analyzed using drought frequency and Sen slope, and the correlation between meteorological drought and agricultural drought is also discussed. The results show that: 1) The drought frequency in the western Loess Plateau is generally higher than that in the east. The trends of meteorological droughts are different from that of agricultural droughts in space. The meteorological droughts in the western and northern Loess Plateau show no significant slowing down trends, and the meteorological droughts in the eastern and southern Loess Plateau show no significant increasing trends. However, the agricultural droughts in most areas are slowing down, especially in the areas on both sides of the line with a rainfall of 400 mm. 2) The frequency of meteorological drought in summer and autumn is higher than that in spring and winter. The frequency of agricultural drought is the highest in spring and the second in summer. VCI is more suitable for real-time monitoring of agricultural drought. 3) In different seasons, the lag time of agricultural drought to meteorological drought is different. The lag time is about 2 months in winter, 1 month in spring, and less than 1 month in summer and autumn. There is a good positive correlation between SPI-12 and VCI in one-harvest regions of the Loess Plateau. The results can provide scientific basis for drought monitoring and early warning, drought regionalization and drought disaster risk assessment in the Loess Plateau.

张华, 徐存刚, 王浩.

2001—2018年西北地区植被变化对气象干旱的响应

[J]. 地理科学, 2020, 40(6):1029-1038.

DOI:10.13249/j.cnki.sgs.2020.06.019      [本文引用: 2]

基于2001?2018年逐月的MODIS NDVI数据，以归一化植被指数（Normalized Difference Vegetation Index, NDVI）和植被状态指数（Vegetation Condition Index, VCI）作为植被生长状况指标，结合2001?2018年的月降水和月均温数据计算的标准化降水蒸散指数值，分析西北地区植被状况和气象干旱指数的变化趋势及其空间分布特征，以及多时间尺度下植被对气象干旱的响应。结果表明：2001?2018年西北地区植被的生长状况整体呈好转趋势，但空间分布上差异明显，东部植被改善状况高于中西部地区。近18 a西北地区5种不同时间尺度标准化降水蒸散指数（Standardized Precipitation Evapotranspiration Index，SPEI）均值整体上均呈增加趋势，表明干旱程度降低；空间上，干旱化趋势整体上表现为中西部高，东部低。植被生长状况在大部分区域均与SPEI呈现不同程度的正相关，总体表现为，西北地区东部植被对气象干旱的响应程度最高，西部次之，中部最低；不同植被类型中，草地对SPEI-12的响应最强，耕地次之，而林地的响应最弱；各植被类型在生长季的多数月份中对SPEI-3和SPEI-12的响应普遍较高。

Zhang H, Xu C G, Wang H.

Response of vegetation change to meteorological drought in Northwest China from 2001 to 2018

[J]. Scientia Geographica Sinica, 2020, 40(6):1029-1038.

DOI:10.13249/j.cnki.sgs.2020.06.019      [本文引用: 2]

Northwest China is one of the regions with severe drought in China. Based on the monthly MODIS NDVI data from 2001 to 2018, the normal difference vegetation index NDVI and the vegetation condition index VCI were used as indicators of vegetation growth status, combined with the SPEI values calculated from the monthly precipitation and monthly mean temperature data from 2001 to 2018. Using the univariate linear regression trend analysis method and Pearson correlation coefficient to analyze the variation trend and spatial distribution characteristics of vegetation condition and meteorological drought in Northwest China, as well as the response of vegetation to meteorological drought in multiple time scales. The results showed that the growth of vegetation in Northwest China showed a trend of improvement in the period from 2001 to 2018, but the spatial distribution was significantly different. The improvement of vegetation in the eastern region was higher than that in the central and western regions. The average SPEI of the five different time scales in Northwest China all showed an overall increase trend in the past 18 years, indicating that the degree of drought decreased. In space, the overall trend of aridification was higher in the central and western regions and lower in the east. The vegetation growth status was positively correlated with SPEI in most areas. The overall performance showed that the vegetation in the eastern part of the Northwest China had the highest response to meteorological drought, followed by the western part and the middle part. Among the different vegetation types, the sensitivity of grassland to SPEI-12 was the strongest, followed by arable land, while the response of forest land was slower; the response of each vegetation type to SPEI-3 and SPEI-12 was generally higher in most months of the growing season.

Dai M, Huang S, Huang Q, et al.

Propagation characteristics and mechanism from meteorological to agricultural drought in various seasons

[J]. Journal of Hydrology, 2022,610:127897.

[本文引用: 1]

吕振涛, 李生宇, 彭中敏, 等.

蒙古国植被对干旱响应的敏感性研究

[J]. 地理研究, 2021, 40(11):3016-3028.

DOI:10.11821/dlyj020210105      [本文引用: 2]

本文分析了蒙古国不同区域、不同土地覆被类型区的植被生长状况对干旱响应的敏感性特征,并探讨了成因。研究表明：① 2001&#x02014;2019年,蒙古国虽然发生了3次较为严重的干旱事件,但整体上干旱程度呈轻微降低趋势,SPEI出现轻微上升,总体增速为0.001%/a;植被生长状况出现了好转,植被指数也呈上升趋势,总体增速为0.15%/a。② 蒙古国大部分区域植被生长受干旱影响较强,各植被指数与干旱指数呈较明显的正相关关系,在NDVI与SPEI-12的相关性分析中,正相关面积占比达76.36%;而在杭爱山脉和肯特山脉以北降水量丰沛的区域和阿尔泰山脉以南荒漠化严重和植被极为稀少区域,植被生长受干旱影响较弱,植被指数与SPEI相关性较小。③ 不同类型植被对干旱响应的敏感性也有差异,其中草地敏感性最强（0.22）,而森林敏感性最弱（-0.04）。④ 干旱持续时间长短,对植被生长影响有较大差异,蒙古国大部分区域的植被对延续3个月至半年的干旱最为敏感,而年内以7月份干旱对植被影响最强。⑤ 植被对干旱响应的敏感性实质是区域水分平衡的植被影响,对同种植被而言多年平均气温越高或降水量越少,植被对干旱响应的敏感性越强,反之则越弱。本文的研究结果对蒙古国因地制宜开展荒漠化防治具有重要参考价值。

Lyu Z T, Li S Y, Peng Z M, et al.

Sensitivity of vegetation responses to drought in Mongolia

[J]. Geographical Research, 2021, 40(11):3016-3028.

[本文引用: 2]

Zhao A, Yu Q, Feng L, et al.

Evaluating the cumulative and time-lag effects of drought on grassland vegetation:A case study in the Chinese Loess Plateau

[J]. Journal of Environmental Management, 2020,261:110214.

[本文引用: 1]

Han Z, Huang Q, Huang S, et al.

Spatial-temporal dynamics of agricultural drought in the Loess Plateau under a changing environment:Characteristics and potential influencing factors

[J]. Agricultural Water Management, 2021,244:106540.

[本文引用: 1]

刘立文, 段永红, 徐立帅, 等.

山西省农业干旱时空变化特征

[J]. 灌溉排水学报, 2020, 39(2):114-121.

[本文引用: 1]

Liu L W, Duan Y H, Xu L S, et al.

Characteristics of drought spatio-temporal transformation in Shanxi Province

[J]. Journal of Irrigation and Drainage, 2020, 39(2):114-121.

[本文引用: 1]

袁瑞强, 龙西亭, 王鹏, 等.

山西省降水量时空变化及预测

[J]. 自然资源学报, 2015, 30(4):651-663.

[本文引用: 1]

Yuan R Q, Long X T, Wang P, et al.

Tempo-spatial variation and forecast of precipitation in Shanxi Province

[J]. Journal of Natural Resources, 2015, 30(4):651-663.

DOI:10.11849/zrzyxb.2015.04.011      [本文引用: 1]

Shanxi Province, where there are widely distributed mountains and karst systems, is the area of surface runoff yield and groundwater recharge for the North China Plain (NCP). Knowledge of tempo-spatial variations and forecast of precipitation in Shanxi could improve the understanding of droughts in the NCP. In this study, spatial distribution of precipitation was clarified based on time series clustering. Features of temporal variance of precipitation were revealed using breaking point test and trend analysis methods. A wavelet- support vector machine model was established to predict the precipitation in the future. As results showed, six sub-regions with different precipitation variation characteristics were identified, which are the North High-latitude Area, the Taihang Alpine Area, the North Mountains and Plateau Area, the Middle Basin Area, the Middle-West Mountains and Plateau Area, and the South Basin Area. More precipitation presented in the Taihang Alpine Area and the South Basin Area, while lessprecipitation occurred in the North High-latitude Area. Precipitation in basins was less than in mountains and plateau area on the same latitude. From 1965 to 1991, annual precipitation decreased significantly in the South, the East, the Middle and the North sequentially. The average rate during 1957 to 2012 is -1.5 mm/a. As a result, the area of Shanxi Province entered an arid period in whole since 1990s. Accordingly, runoff declined and groundwater recharge reduced. Decrease in precipitation of Shanxi in recent decades was evidently one of the driving factors of droughts in the NCP. The study region will get out of the less precipitation period before the year of 2020 with an average gradient of 16 mm/a, which would relieve droughts in the NCP.

Holben B N.

Characteristics of maximum-value composite images from temporal AVHRR data

[J]. International Journal of Remote Sensing, 1986, 7(11):1417-1434.

[本文引用: 1]

Cao S, Zhang L, He Y, et al.

Effects and contributions of meteorological drought on agricultural drought under different climatic zones and vegetation types in Northwest China

[J]. Science of the Total Environment, 2022,821:153270.

[本文引用: 1]

刘远, 周买春.

3种IGBP分类系统的土地覆盖数据在韩江流域的对比分析

[J]. 遥感技术与应用, 2017, 32(3):575-584.

DOI:10.11873/j.issn.1004-0323.2017.3.0575      [本文引用: 1]

土地覆盖是陆面过程模型、生态水文模型的重要输入，对模拟结果影响较大。针对基于IGBP分类系统的土地覆盖数据IGBP DISCover、MODIS和MICLCover，以韩江流域为研究区域，分析了3种数据土地覆盖空间分布，通过人工判读样地Google卫星影像对其精度进行了评价。结果表明：IGBP DISCover土地覆盖类型少，空间变异性小，而MODIS和MICLCover类型多，空间变异性大；从大类看，3种数据土地覆盖空间分布较吻合;3种数据林地划分存在明显差异，其中MICLCover林地成分较好地反映了流域实际情况；MODIS和MICLCover耕地面积与流域实际相符，IGBP DISCover耕地面积明显偏大；3种数据建设用地、裸地、湿地和水域的面积较实际小，说明1 km的空间分辨率过粗，不能细致反映流域中面积较小或较分散的土地覆盖，尤其是IGBP DISCover数据。Google卫星影像样地检验表明，MICLCover精度明显高于IGBP DISCover和MODIS，IGBP DISCover的主要错误是将许多其他土地覆盖划分为耕地，而MODIS的主要错误是将许多其他土地覆盖划分为草地。

Liu Y, Zhou M C.

Comparative analysis on three land cover datasets based on IGBP classification system over Hanjiang River basin

[J]. Remote Sensing Technology and Application, 2017, 32(3):575-584.

[本文引用: 1]

McKee T B, Doesken N J, Kleist J.

The relationship of drought frequency and duration to time scales

[C]// Proceedings of the 8th Conference on Applied Climatology. American Meteorological Society,1993:179-183.

[本文引用: 1]

Svoboda M, Fuchs B A. Handbook of drought indicators and indices[M]. Switzerland: World Meteorological Organization Geneva,2016:13-14.

[本文引用: 1]

Ma L, Huang Q, Huang S, et al.

Propagation dynamics and causes of hydrological drought in response to meteorological drought at seasonal timescales

[J]. Hydrology Research, 2022, 53(1):193-205.

[本文引用: 1]

Zhou Z, Liu S, Ding Y, et al.

Assessing the responses of vegetation to meteorological drought and its influencing factors with partial wavelet coherence analysis

[J]. Journal of Environmental Management, 2022,311:114879.

[本文引用: 2]

Han Z, Huang S, Huang Q, et al.

Propagation dynamics from meteorological to groundwater drought and their possible influence factors

[J]. Journal of Hydrology, 2019,578:124102.

[本文引用: 1]

吕纯月. 基于SPI指数的中国夏季干旱区域性特征及环流异常研究[D]. 南京: 南京信息工程大学, 2021.

[本文引用: 2]

Lyu C Y. Study on regional characteristics and circulation anomalies of summer drought in China based on SPI index[D]. Nanjing: Nanjing University of Information Science & Technology, 2021.

[本文引用: 2]

马景钊, 郝璐.

基于SPI和SPEI指数的锡林郭勒草原干旱时空变化特征

[J]. 草业科学, 2021, 38(12):2327-2339.

[本文引用: 1]

Ma J Z, Hao L.

Temporal and spatial variation of drought in the Xilingol grassland based on the standardized precipitation index and standardized precipitation evapotranspiration index

[J]. Pratacultural Science, 2021, 38(12):2327-2339.

[本文引用: 1]

Kogan F N.

Application of vegetation index and brightness temperature for drought detection

[J]. Advances in Space Research, 1995, 15(11):91-100.

[本文引用: 1]

Li W, Wang Y, Yang J.

Cloudy region drought index (CRDI) based on long-time-series cloud optical thickness (COT) and vegetation conditions index (VCI):A case study in Guangdong,south eastern China

[J]. Remote Sensing, 2020, 12(21):3641.

[本文引用: 1]

Dragomiretskiy K, Zosso D.

Variational mode decomposition

[J]. IEEE Transactions on Signal Processing, 2014, 62(3):531-544.

[本文引用: 2]

李文武, 石强, 王凯, 等.

基于变分模态分解和深度门控网络的径流预测

[J]. 水力发电学报, 2020, 39(3):34-44.

[本文引用: 2]

Li W W, Shi Q, Wang K, et al.

Runoff prediction based on variational mode decomposition and deep gated network

[J]. Journal of Hydroelectric Engineering, 2020, 39(3):34-44.

[本文引用: 2]

吴文轩, 王志坚, 张纪平, 等.

基于峭度的VMD分解中k值的确定方法研究

[J]. 机械传动, 2018, 42(8):153-157.

DOI:10.16578/j.issn.1004.2539.2018.08.030      [本文引用: 2]

在变分模态分解（Variational mode decomposition，VMD）中，因为分解层数k值属于自定义变量，所以在取值时分解结果会随着k值的变化而得到不同的结果，k值的取值直接影响着结果的准确性，k值取得过大或者过小都会对结果造成影响。基于以上问题，提出了一种利用峭度确定k值的方法。选取k值为2~n的整数，计算当k值为2~n时，每一个k值相关系数最大分量的峭度，并绘制峭度的变化曲线，若在该曲线内峭度没有峰值且单调递增，则继续计算当k为n+1时相关系数最大分量的峭度值，重复以上步骤，以峭度最大作为优化的标准，当峭度最大的时候，k为最佳值。用实际故障信号验证了方法的可行性，为VMD方法的研究提出了一种新的思路。

Wu W X, Wang Z J, Zhang J P, et al.

Research of the method of determining k value in VMD based on kurtosis

[J]. Journal of Mechanical Transmission, 2018, 42(8):153-157.

[本文引用: 2]

邢愿. 基于不同时间尺度的贵州省农业干旱对气象干旱的响应机制[D]. 贵阳: 贵州师范大学, 2021.

[本文引用: 1]

Xing Y. Response mechanism of agricultural drought to meteorological drought based on different time scales in Guizhou Province[D]. Guiyang: Guizhou Normal University, 2021.

[本文引用: 1]

丁一汇, 司东, 柳艳菊, 等.

论东亚夏季风的特征、驱动力与年代际变化

[J]. 大气科学, 2018, 42(3):533-558.

[本文引用: 1]

Ding Y H, Si D, Liu Y J, et al.

On the characteristics,driving forces and inter-decadal variability of the East Asian summer monsoon

[J]. Chinese Journal of Atmospheric Sciences, 2018, 42(3):533-558.

[本文引用: 1]

晏利斌.

1961—2014年黄土高原气温和降水变化趋势

[J]. 地球环境学报, 2015, 6(5):276-282.

[本文引用: 1]

Yan L B.

Characteristics of temperature and precipitation on the Loess Plateau from 1961 to 2014

[J]. Journal of Earth Environment, 2015, 6(5):276-282.

[本文引用: 1]

侯青青, 裴婷婷, 陈英, 等.

1986—2019年黄土高原干旱变化特征及趋势

[J]. 应用生态学报, 2021, 32(2):649-660.

DOI:10.13287/j.1001-9332.202102.012      [本文引用: 1]

干旱作为极端气候事件之一,其频率和强度的变化影响到区域水资源,而干旱半干旱区植物生长的主要限制因素是水分,因此,研究黄土高原干旱时空特征及未来变化趋势对当地的生态环境具有重要意义。本研究基于1986—2019年降水和温度逐月格点数据,计算标准化降水蒸散发指数(SPEI)和干旱发生频率,并运用Mann-Kendall检验和Sen斜率估计方法,探讨了黄土高原的年、季尺度干旱时空分布及变化特征,最后利用NAR神经网络结合Hurst指数对黄土高原未来干旱趋势进行预测。结果表明: 研究期间,黄土高原总体呈现干旱化趋势,且年际和季节尺度的干旱发生频率在空间上差异较大。其中,年际、春季和冬季以黄土高原东南部和西部干旱发生频率最高,夏季和秋季以西北部干旱发生频率最高。夏季以中度干旱发生频率最高,年际及其他季节以轻度干旱发生频率最高。黄土高原春、夏季呈现干旱化趋势,秋、冬季研究区大部分区域干旱趋势减轻。黄土高原年际、春季、夏季的SPEI值在未来一段时间内仍处于下降趋势,即干旱化趋势加重,且夏季的Hurst指数最大,持续性变化最强,未来持续干旱的可能性高于其他季节。

Hou Q Q, Pei T T, Chen Y, et al.

Variations of drought and its trend in the Loess Plateau from 1986 to 2019

[J]. Chinese Journal of Applied Ecology, 2021, 32(2):649-660.

[本文引用: 1]

安彬, 肖薇薇, 朱妮, 等.

近60 a黄土高原地区降水集中度与集中期时空变化特征

[J]. 干旱区研究, 2022, 39(5):1333-1344.

DOI:10.13866/j.azr.2022.05.01      [本文引用: 1]

全球气候变暖背景下,厘清黄土高原地区降水集中度（PCD）和集中期（PCP）时空变化规律,对该区水土保持及水资源利用具有重要意义。基于1960&#x02014;2019年黄土高原地区55个气象站点逐日降水观测资料,利用趋势分析、空间插值、相关分析等方法,分析了近60 a及退耕还林（草）工程前后黄土高原PCD与PCP时空变化特征。结果表明：（1） 1960&#x02014;2019年黄土高原年均降水量呈减少趋势,年均PCD逐渐减弱、PCP逐渐提前。（2） 黄土高原地区年均降水量和PCD自东南向西北递减,PCP自东向西逐渐递减,但差异不大。PCD变化趋势自东北向西南呈&#x0201c;下降-上升-下降&#x0201d;相间分布,PCP表现出东部推迟、西部提前的态势。（3） 退耕还林（草）工程后,黄土高原年均降水量增加、PCD减弱、PCP推迟。降水量自东向西呈&#x0201c;偏少-偏多&#x0201d;相间分布,其变化趋势以偏多为主;PCD自东北向西南呈&#x0201c;偏低-偏高-偏低&#x0201d;分布规律,其变化趋势以偏低为主;PCP以偏高为主,偏高幅度呈南多北少、东多西少态势,其趋势变化也以偏多为主。（4） 1960&#x02014;2019年黄土高原地区PCD、PCP与年降水量皆以正相关为主。

An B, Xiao W W, Zhu N, et al.

Temporal and spatial variations of precipitation concentration degree and precipitation concentration period on the Loess Plateau from 1960 to 2019

[J]. Arid Zone Research, 2022, 39(5):1333-1344.

DOI:10.13866/j.azr.2022.05.01      [本文引用: 1]

With global warming as a foundation, it is significantly important that the law of precipitation concentration degree (PCD) and precipitation concentration period (PCP) changes for the Loess Plateau are clarified for the conservation of regional soil and water and utilization of water resources. Based on the daily precipitation observation data compiled from 55 meteorological stations on the Loess Plateau from 1960 to 2019, the temporal and spatial variation characteristics of PCD and PCP on the Loess Plateau in the last 60 years, along with before and after the project of returning farmland to forest (grassland), were analyzed using trend analysis, spatial interpolation, and correlation analysis. The average annual precipitation, PCD, and PCI on the Loess Plateau had a decreasing (downward) trend from 1960 to 2019 according to the results. In addition, the PCD gradually weakened, and the PCI shifted to an earlier date. The spatial distribution characteristics of annual precipitation, PCD and PCI on the Loess Plateau were complex. Annual precipitation and the PCD decreased from southeast to northwest, while PCP gradually decreased from east to west. However, the difference was not significant. The variation trend of the PCP had an alternating distribution law of “decline rise decline” from northeast to southwest. Overall, the PCP had a trend of delaying in the east and advancing in the west. After the project of returning farmland to forest (grassland) was completed, the average annual precipitation increased, the PCD weakened, and the PCP was postponed. The precipitation on the Loess Plateau is distributed alternately from east to west, and its change trend is mainly more. The distribution law of PCD is “low high low” from northeast to southwest, and its variation trend is mainly low. The PCP is mainly on the high side, and the range of the high side is more in the south than the north and more in the east than the west, and the trend change is mainly on the high side. The PCD, PCI, and annual precipitation in the Loess Plateau are mainly positively correlated.

薛少博, 李鹏, 于坤霞, 等.

2002—2020年黄土高原土壤水变化及其相关性分析

[J]. 水土保持学报, 2021, 35(5):221-226.

[本文引用: 1]

Xue S B, Li P, Yu K X, et al.

Soil water change and its correlation analysis in the Loess Plateau from 2002 to 2020

[J]. Journal of Soil and Water Conservation, 2021, 35(5):221-226.

[本文引用: 1]

史尚渝, 王飞, 金凯, 等.

黄土高原地区植被指数对干旱变化的响应

[J]. 干旱气象, 2020, 38(1):1-13.

[本文引用: 1]

Shi S Y, Wang F, Jin K, et al.

Response of vegetation index to meteorological drought over loess plateau

[J]. Journal of Arid Meteorology, 2020, 38(1):1-13.

[本文引用: 1]

王一, 郝利娜, 许强, 等.

2001—2019年黄土高原植被覆盖度时空演化特征及地理因子解析

[J]. 生态学报, 2023, 43(6):2397-2407.

[本文引用: 1]

Wang Y, Hao L N, Xu Q, et al.

Spatio-temporal variations of vegetation coverage and its geographical factors analysis on the Loess Plateau from 2001 to 2019

[J]. Acta Ecologica Sinica, 2023, 43(6):2397-2407.

[本文引用: 1]

付建新.

山西黄河流域不同土地利用类型NDVI时空变化及其对气温、降水的响应

[J]. 水土保持研究, 2023, 30(3):364-372.

[本文引用: 1]

Fu J X.

Temporal and spatial changes of NDVI of different land covers and their responses to temperature and precipitation in the Yellow River basin of Shanxi

[J]. Research of Soil and Water Conservation, 2023, 30(3):364-372.

[本文引用: 1]

尉毓姣, 朱琳, 曹鑫宇, 等.

基于转移函数分析的蒙东地区不同类型植被变化对干旱的响应

[J]. 地球科学, 2023, 48(9):3539-3551.

[本文引用: 1]

Wei Y J, Zhu L, Cao X Y, et al.

Response of different types of vegetation changes to drought in eastern inner Mongolia based on transfer function analysis

[J]. Earth Science, 2023, 48(9):3539-3551.

[本文引用: 1]

吴林霖, 王思远, 马元旭, 等.

中亚地区植被对气候变化的响应机制初探

[J]. 遥感学报, 2022, 26(11):2248-2267.

[本文引用: 1]

Wu L L, Wang S Y, Ma Y X, et al.

Response of vegetation to climate change in Central Asia with remote sensing and meteorological data

[J]. National Remote Sensing Bulletin, 2022, 26(11):2248-2267.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编：100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发