干旱作为一种反复出现的极端气候事件, 对农牧业、水资源以及自然和社会生态系统产生了许多负面影响^[1]。据国际灾害数据库统计, 1960—2016年全球因干旱造成的损失平均每年约2210亿美元。然而, 随着全球变暖强度的增加, 干旱的发生频率和持续时间正在呈现出明显的上升趋势^[2-3], 全球陆地的极端干旱面积也在不断扩大^[4], 其中以亚欧大陆和非洲大陆最为显著^[5-6]。大量的事实已经揭示了近几十年中国干旱发生频率、持续时间和范围不断扩大的事实^[7], 尤其是发生在我国北方地区的重度和极端干旱事件更加频繁^[8-10]。此外, 旱灾的发生区域也开始不断地向我国南部和东部的湿润、半湿润地区扩展^[11], 这种趋势可能继续在21世纪持续下去^[12]。有研究发现, 如果中国不对干旱灾害给予应有的重视和积极有效的应对, 到2030年, 仅中国东北地区, 就可能造成3500万农民的农业收入损失一半以上^[13]。

干旱指标是干旱监测的基础与核心。为了评估气候干旱的频率、持续时间、严重程度和空间范围, 已经开发了数百个干旱指数^[14]。对于同一干旱事件的特征, 不同的干旱指数因所需变量和计算方法差异, 使其结果也不尽相同。干旱指标均是建立在特定的时空范围内, 有其相应的尺度^[15], 因为在不同流域和地区, 干旱变化特征及表现会因降水、气温、土壤水文等因素的空间分异而存在较大差异, 因此, 研究干旱指数的区域适用性具有重要的科学意义^[16]。然而, 对于何种干旱指数能够充分表征干旱条件方面, 学术界目前并未达成共识。但是, 有一些干旱指数更受欢迎, 例如由McKee等^[17]提出的标准化降水指数SPI不仅计算简单, 且能较好地反映不同时间尺度和不同区域的干旱状况。然而, 在全球变暖的情况下, SPI仅考虑降水资料, 未考虑影响干旱的其他因素如温度、蒸散等, 限制了该方法进一步的适用。此外, 考虑了温度和前期天气条件对干旱的影响的Palmer干旱指数(PDSI, Palmer Drought Severity Index), 其物理机制较为明确, 非常适用于干旱对全球变暖响应的研究^[18], 但其计算繁杂, 对资料要求较高, 部分参数只能靠经验估计, 致使计算精度大大降低。加之时间尺度固定(9—12个月之间)而无法有效地应用于干旱的多时间尺度研究中^[19]。标准化降水蒸散指数(SPEI, Standardized Precipitation Evapotranspiration Index)是Vicente-Serrano^[20]为检测全球干旱分布状况, 在SPI指数的基础上提出的, 可以反映干旱的持续时间和积累, 拥有多尺度特征和能够衡量温度变化对干旱影响的双重优点, 使其成为评估、监测和评估全球变暖背景下干旱的理想指标^[21-23]。已有学者基于SPEI方法评价了中国西北^[24-25]、西南^[26]、东北^[27-28]、华中^[29]和华北^[30]等区域的干旱时空格局, 并得到了较好的验证。

呼伦贝尔草原作为亚欧大陆草原景观的一部分^[31], 拥有森林、草原、疏林沙地和沼泽湿地等多种生态系统类型, 是全球气候变化过程中的生态敏感区和脆弱区^[32-33]。同时, 该区也是我国北方地区重要的防风固沙和水源涵养生态功能区、京津上游的生态屏障区、生物多样性维持区和重要的畜牧业生产基地^[34]。然而, 干旱是生态系统初级生产力变化的重要驱动因素^[35], 尤其在干旱和半干旱地区, 干旱还可能引发一系列的环境问题^[36-37], 以及包括食物安全、环境难民、经济发展等农牧业干旱、水文干旱和社会经济干旱等社会经济问题。近年来呼伦贝尔草原区的气候变化, 尤其是干旱的频发对该区草原畜牧业和旱作农业产生了重大影响, 集中表现在草地群落结构简单化、物种多样性减少、土地沙化不断扩展、鼠虫害频发等^{[34, 38]}。因此, 对呼伦贝尔草原气象干旱发生的特征进行准确、有效的监测, 不仅可为该区农牧业干旱、水文干旱和社会经济干旱提供研究背景^[39]和干旱风险的管理依据^[40], 还可为该区农牧民生活的改善以及生态环境的保护提供决策支持。

1 研究区概况

呼伦贝尔草原位于内蒙古自治区东北部, 大兴安岭以西, 位于47°05′—53°20′N, 115°31′—123°00′E。行政范围包括海拉尔市、满洲里市、陈巴尔虎旗、新巴尔虎右旗、新巴尔虎左旗、鄂温克族自治旗, 总面积约83350 km²。其中, 陈巴尔虎旗、新巴尔虎右旗、新巴尔虎左旗、鄂温克族自治旗是著名的牧业四旗, 四旗的草原面积占到全市草原总面积的90%以上。

呼伦贝尔草原处于干旱半干旱、季风与非季风以及农区和牧区的过渡地带^[41-42], 地形为平坦、辽阔、波状起伏的高平原, 海拔在650—700 m。该区属温带大陆性季风气候, 春季干旱少雨, 大风频发；夏季温凉短促, 水热集中；秋季降温迅速, 日温差大, 秋霜早；冬季严寒漫长, 冷空气活动频繁。各区域间气候差异较小, 年均温约-1.0—1.0℃, 年均降水量250—350 mm。该区植被类型复杂多样, 草原植物资源约1000余种, 是欧亚草原的重要组成部分^[31]。该区东部由于受到大兴安岭山地的影响, 气候较湿润, 属于森林草原交错带；从大兴安岭西麓向西至呼伦湖方向, 植被类型由森林草原, 草甸草原过渡到典型草原。呼伦贝尔草原区域分布着额尔古纳河、根河、海拉尔河、伊敏河等河流, 以及零星分布的大小湖泊, 河流和湖泊沿岸一般发育各类型草甸或者沼泽^[43]。

2 数据与方法 2.1 数据

气象及灾情数据来源于国家气象科学数据共享服务平台(http://data.cma.cn/), 该数据集已经过基本的质量控制, 包括气候极值范围检查、内部一致性检查和时间一致性检查, 数据可靠性和连续性均能满足研究的需求。本文选取呼伦贝尔典型草原和草甸草原区的4个气象站点(图 1)1960—2017年的逐月降水量及月均气温：海拉尔(49.22°N, 119.75°E)、满洲里(49.57°N, 117.43°E)、新巴尔虎左旗(48.22°N, 118.27°E)和新巴尔虎右旗(48.67°N, 116.82°E)。此外, 为验证本文计算的SPEI指数在呼伦贝尔草原区的有效性, 选取了1960—1999年我国重大干旱(重旱、特大旱)灾害事件以及满洲里和鄂温克族自治旗的农业气象灾情旬值数据和2005—2016年《中国气象灾害年鉴》干旱事件发生的记录作为验证。

2.2 方法 2.2.1 标准化降水蒸散指数(SPEI)

标准化降水蒸散指数(SPEI)是降水量与潜在蒸散发量差值序列的累积概率做正态标准化后的指数。为了揭示以草原生态系统为主(天然草场面积占研究区总面积的80%以上)的区域整体干旱特征, 对差异较小的4个气象站点逐月的气温和降水取均值来计算其标准化降水蒸散指数(SPEI)。在验证该指数表征呼伦贝尔草原区干旱事件发生的有效性时, 则分别计算各站点的SPEI指数。SPEI指数是降水量与潜在蒸散发量差值序列的累积概率做正态标准化后的指数。首先, 采用Thornthwaite法计算潜在蒸散量, 然后计算逐月降水和蒸散的差值, 并建立不同时间尺度的水分盈亏累积序列。由于原始数据序列中可能存在负值, 所以, 采用3参数的Log-logistic概率分布对累积概率密度进行标准化, 最终计算得出SPEI值, 计算过程如下^[19]：

首先, 计算潜在蒸散量(PET), Vicente-Serrano推荐使用Thornthwaite方法：

(1)

(2)

式中, PET为潜在蒸散量；T_i为月平均温度；H为年热量指数；A为常数, 由H决定, A=0.49+0.179H-0.0000771H²+0.000000675H³。

其次, 计算逐月降水量与蒸散量的差值：

(3)

式中, D_i为降水量与蒸散量的差值, P_i为月降水量, PET_i为月蒸散量。

第三, 采用3个参数的Log-logistic概率分布对D_i数据序列进行正态化, 计算每个数值对应的SPEI指数：

(4)

其中, 参数α、β、γ的计算如下：

(5)

(6)

(7)

其中, Γ为阶乘函数, ω₀、ω₁、ω₂为数据序列D_i的概率加权矩：

(8)

(9)

N为参与计算的月份数。

最后对累计概率密度进行标准化：

(10)

当累计概率P≤0.5时：

(11)

(12)

式中, 常数c₀, c₁, c₂, d₁, d₂, d₃同SPI计算过程, c₀ =2.515517, c₁=0.802853, c₂=0.010328, d₁=1.432788, d₂=0.189269, d₃=0.001308。

由于该方法是基于“历史同月的累积水分亏缺量服从Log-logistic分布”这一假设, 因此, 为了验证呼伦贝尔草原累积水分亏缺量是否符合Log-logistic分布, 本文对1个、3个和12个月尺度下的累积水分亏缺量序列与Log-logistic分布进行了K-S(Kolmogorov-Smirnov)检验。结果表明, 在1个、3个和12个月尺度下, K-S检验的概率P值(分别为0.241、0.167、0.508)均大于显著性水平α=0.05, 没有足够理由拒绝零假, 即累积水分亏缺量样本来自Log-logistic分布总体。因此, 可判定该序列来自Log-logistic分布总体, 基于该分布的SPEI指数在呼伦贝尔草原干旱表征上具备数学统计理论基础。

此外, 为了反映研究区干旱过程的细节与总体特征, 分析各个站点月、季以及年际干旱演变特征(分别记作SPEI-1、SPEI-3和SPEI-12)。3个月时间尺度SPEI与草原生长季的关系密切；12个月时间尺度SPEI能较清晰的反映长期干旱变化特征, 因此, 本文以3个月为时间尺度(SPEI-3)的5月、8月、11月和次年2月作为表征春季、夏季、秋季、冬季，12个月为时间尺度的SPEI(SPEI-12)表征年际干旱, 并参照相关文献^{[19, 23]}对干旱等级进行了划分(表 1)。

2.2.2 Hurst指数和R/S分析法

Hurst指数通常是用来分析时间序列的分形特征和长期记忆过程。它是由Hurst在分析尼罗河水文数据时提出的, 经常用于分析长时间序列相关性。其中, R/S分析法(也被称为重标极差分析法)被广泛用于计算H指数。R/S原理简述如下^[44]：

设有时间序列, {ξ(t)}, t=1, 2, 3…, 对于任意正整数≥1, 有均值序列：

(1)

由此求得累计离差：

(2)

极差R定义为：

(3)

标准差S定义为：

(4)

R、S、τ满足一般关系式：

(5)

式中, c为常数；R(τ)/S(τ)为重标极差；H为Hurst指数。

H范围从0到1。当H=0.5时, 表明该时间序列为随机序列；H>0.5时表示该过程具有连续特性, 并且未来趋势与过去趋势一致, 而H < 0.5时则表示未来趋势与过去趋势相反。

2.2.3 其他方法

采用气候倾向率法研究SPEI的时空变化趋势, Mann-Kendall统计检验用于研究SPEI时间序列的突变, Morlet小波分析用于计算本研究中干旱发生的周期性。

3 结果分析 3.1 呼伦贝尔草原气候变化趋势

由图 2可知, 1960—2017年呼伦贝尔草原多年平均降水量为286.66 mm, 且以3.65 mm/10a的速率呈不显著减少趋势, 其中, 1986年和1998年分别是近58年来该区域降水量最少和最多的两个年份, 年均降水量分别为137.4 mm和567.225 mm。过去58年该区域多年平均气温为-0.06℃, 并以0.35℃/10a的速率呈上升趋势(P < 0.001), 其中, 1969年和2007年分别是过去58年来该区域年平均温度最高和最低的年份, 其值分别为-2.30℃和2.31℃。

由图 3可知, 过去58年呼伦贝尔草原春季多年平均降水量为32.68 mm, 占年平均降水量的11.40%, 并以1.75 mm/10a的速率呈不显著增多趋势, 最低值和最高值出现的年份分别是1972年(6.65 mm)和2005年(68.6 mm)的春季。过去58年该区域多年春季平均气温为1.58℃, 并以0.46℃/10a的速率呈上升趋势(P < 0.001)。其中, 1966年和1998年分别是多年春季平均气温的最低值和最高值, 分别为-1.77℃和5.08℃。同样地, 夏季多年平均降水量为202.72 mm, 占年平均降水量的70.72%, 过去58年以7.75 mm/10a的速率呈减少趋势, 最低值和最高值出现的年份分别是2016年(80.05 mm)和1998年(432.68 mm)。该区域夏季多年平均气温为19.28℃, 过去58年以0.39℃/10a的速率呈上升趋势(P < 0.001), 最低值和最高值出现的年份分别是1983年(17.44℃)和2001年(22.12℃)的夏季。

秋季多年平均降水量为43.42 mm, 占年平均降水量的15.15%, 过去58年以1.142 mm/10a的速率呈增加趋势, 最低值和最高值出现的年份分别是1986年(4.8 mm)和1998年(103.55 mm)。该区域秋季多年平均气温为0.25℃, 过去58年以0.28℃/10a的速率呈上升趋势(P=0.002), 最低值和最高值出现的年份分别是1981年(-1.68℃)和2004年(2.79℃)。冬季多年平均降水量为7.82 mm, 占年平均降水量的2.73%, 过去57年以1.33 mm/10a的速率呈增加趋势(P < 0.001), 最低值和最高值出现的年份分别是1973年(1.9 mm)和2012年(17.925 mm)。该区域冬季多年平均气温为-21.34℃, 过去57年以0.24℃/10a的速率呈不显著上升趋势, 最低值和最高值出现的年份分别是1986年(-26.59℃)和2004年(-17.35℃)。

3.2 SPEI指数的年际变化特征

图 4为呼伦贝尔草原平均SPEI指数年际变化及其Mann-Kendall (M-K)检验曲线。呼伦贝尔草原的年尺度的SPEI值在1960—2017年间呈下降趋势, 线性倾向率为-0.218/10a(P < 0.01), 表明, 近58年来呼伦贝尔草原呈干旱化趋势, 其中, 1999年和2001年分别是该地区近58年来最湿润和最干旱的两个年份, 其SPEI值分别为2.012和-2.284。UF曲线显示近58年呼伦贝尔草原平均SPEI指数经历了下降上升和下降的趋势, 其中UF值在1984—1994年间(除了1986和1987年)为正值, 表明SPEI指数在此期间总体呈上升趋势, 1984年之前以及1994年之后, UF值均小于0, 尤其是2004年之后, UF曲线超出了0.05的置信区间, 说明SPEI指数下降的趋势显著。在置信区间内, UF和UB曲线相交于1998年, 说明1998年是SPEI指数突变的开始, 也表明了呼伦贝尔草原自1998年之后干旱开始加剧。

图 5为呼伦贝尔草原季节尺度的SPEI指数年际变化及M-K突变检验。由春季SPEI指数的年际变化可知, 近58年来该地区春季呈干旱化趋势, 其SPEI指数为-0.142/10a(P=0.05), 2005年和2017年分别是最湿润和最干旱的两个年份, 其SPEI值分别为1.564和-2.077。在显著性水平0.05的临界线之间, UF、UB曲线相交于1982年, 这是呼伦贝尔草原春季干旱突变的开始。由夏季SPEI指数的年际变化可知, 1960—2017年间呼伦贝尔草原夏季呈现显著的干旱化趋势(P < 0.01)。其中, 1998年和2006年的夏季是最湿润和最干旱的两个年份, 其SPEI指数分别为2.023和-1.9。由M-K检验可知, 1998年是呼伦贝尔草原夏季干旱突变的开始。由秋季SPEI指数的年际变化可知, 呼伦贝尔草原近58年秋季呈现不显著干旱化趋势。其中, 1970和2007年的秋季是最湿润和最干旱的两个年份, 其SPEI指数分别为1.6和-2.333。UF和UB曲线交于1995年, 可知1995年是呼伦贝尔草原秋季干旱突变的开始。由冬季SPEI指数的年际变化可知, 呼伦贝尔草原近58年冬季呈现湿润化趋势, 其中, 1973年和2012年的冬季分别是最干旱和最湿润的两个年份, 其SPEI指数分别为-1.960和1.665。UF和UB曲线交于1998年, 可知1998年冬季是呼伦贝尔冬季湿润化突变的开始。

3.3 干旱事件频率分析

根据表 1中SPEI指数评估等级, 统计了近58年呼伦贝尔草原4个站点每月的SPEI值, 得到该区域不同年代干旱发生的频率和不同季节干旱发生频次(图 6)。容易看出, 1960—2017年呼伦贝尔草原不同等级干旱发生的频率均呈增加趋势, 其中, 中旱发生频率最高, 重旱次之, 极旱最少。此外, 在1960s、1970s、1980s和1990s, 中旱占同年代发生干旱的比率在74.16%—78.87%, 而在21世纪初的18年, 中旱发生频率的占比下降到55.05%, 同时, 与其他年代相比, 该时期也是中旱及以上等级干旱发生频率最高的年代。可见, 21世纪初期该区域不但干旱发生频率增加, 而且发生重旱和极旱的频率显著提高。例如, 2000—2017年极旱发生频率为12.96%, 显著高于其他年代, 而重旱发生频率为32.41%, 是发生最少年代(1970s)的4倍。

就不同季节干旱发生频次而言, 过去58年来, 夏季发生的干旱频次最多, 为123次, 其次为冬季(120次), 春季(116次)和夏季(112次), 但总体而言差异不大。从不同等级干旱在4个季节的分布可以看出, 中旱均是发生频次最多的干旱, 重旱次之, 极旱最少。此外, 冬季是发生中旱次数最多的季节(75次), 夏季是发生重旱最多的季节(45次), 秋季是发生极旱次数最多的季节(21次)。

3.4 干旱趋势分析 3.4.1 Morlet小波分析

对呼伦贝尔草原1960—2017年间的年尺度与季节尺度的SPEI值进行了小波分析, 选用了水文气象常用的Morlet连续复小波变换, 其中, 小波系数大小表示信号强弱, 等值线中心为正表示该年份偏湿(干旱化指数偏高), 为负则表示该年份偏旱。由图 7的小波系数等值线图及小波方差图可知, 在年尺度上, 明显存在着8—11年和17—19年的短周期和40—50年的长周期。其中, 在17—19年的时间尺度上出现4次明显的旱湿交替；8—11年的尺度上有明显的7次旱湿交替, 40—50年尺度上有1次明显的旱湿交替。结合小波方差图得知, 在8年、11年、18年、30年和44年时间尺度上震荡明显, 且44年、11年和8年分别为第一主周期和第二、三周期。由于当前(2017年)的小波系数图在44年长周期和11年短周期上并未闭合, 说明呼伦贝尔草原在11年短期内将持续偏旱趋势, 而在44年长周期上将持续偏湿趋势。

从图 8可以发现, 呼伦贝尔草原近58年的春季SPEI在4—6年、10—12年、20—25年和43—48年尺度上存在着明显的高低值交替现象, 且由小波方差可以看出, 22年、4年和44年分别为第一主周期和第二、三周期。在22年的主周期上, 1960—2017年, 呼伦贝尔草原春季SPEI经历了4次旱湿交替的周期, 且当前该区春季正处于旱湿交替的过渡时期, 未来短期内仍将有偏旱的趋势。从夏季SPEI的小波系数和小波方差图容易发现, 该区明显存着9年、17年和44年的震荡周期, 其中, 在44年周期上的信号最强, 为第一主周期, 9年和17年分别为第二、三周期。在44年尺度上, 该区域夏季目前处于偏湿期, 且有持续偏湿的趋势；在9年和17年尺度上, 该区域夏季正目前处于偏旱期, 且未来有持续偏旱的趋势。从秋季SPEI的小波系数和小波方差图容易发现, 在6年、15年、29年和44年的周期上存在较明显的震荡现象, 其中, 29年尺度的能量最强, 为第一主周期, 44年和15年分别为第二、三周期。在29年和44年的周期上, 目前正处于偏湿时期, 且有持续偏湿的趋势。然而, 在15年尺度上, 下一周期有转向偏干的趋势。由呼伦贝尔草原冬季SPEI的小波方差图可知, 冬季存在44年的第一主周期, 29年和17年分别为第二、三周期。在44年和29年尺度上, 目前该区冬季正处于偏干期, 在17年尺度上, 则处于偏湿期。

3.4.2 基于Hurst指数的未来干旱趋势分析

对于呼伦贝尔草原年际干旱的未来趋势运用R/S分析法进行预测(表 2), 年际SPEI的Hurst指数0.834>0.5, C_M>0, 表明时间序列前后具有持续性。未来呼伦贝尔草原的年际SPEI指数变化趋势与过去58年变化趋势一致, 即SPEI指数将继续呈现下降趋势, 干旱化的态势也将继续加剧, 这与11年尺度的小波分析结果较一致。在季节尺度上, SPEI的Hurst指数均>0.5, 且C_M>0, 表明未来四季的SPEI指数将与过去58a变化的趋势一致。然而, 不同季节变化趋势的强度不同, 其中, 冬季的SPEI的Hurst指数最大(0.920), 持续性变化最强, 表明呼伦贝尔草原未来冬季SPEI指数持续上升的可能性高于其他各季, 即未来研究区的冬季将呈现湿润的趋势, 这与小波分析的第三主周期(17年尺度)分析的结果也较一致。夏季SPEI的Hurst指数次之(0.802), 表明该区未来夏季SPEI指数将继续呈现下降趋势, 干旱化将加剧, 这与小波分析的第二、三周期(9年和17年尺度)分析的结果较一致。秋季SPEI和春季SPEI的Hurst指数也均大于0.5, 表明呼伦贝尔草原未来春季和秋季的干旱趋势将继续加剧, 但可能性小于冬季和夏季。秋季SPEI的Hurst分析结果与其小波分析的第三主周期(15年尺度)分析的结果一致, 春季SPEI与其第一主周期(22年尺度)的分析结果较一致。

3.5 SPEI指数在呼伦贝尔草原的适用性分析

SPEI的适用性分析在全国尺度^[45]、东北地区^[46]、西北地区^[23]、长江中下游地区^[47]和内蒙古地区^[48]均进行过探讨, 然而, 在以草原景观尺度的适用性的研究较少。表 3为1960—2016年来呼伦贝尔草原区典型干旱事件与同期SPEI验证结果, 可以看出, 在本区域实际发生的干旱事件与本文计算得到的同期SPEI所表征的干旱事件吻合度较高, 表明SPEI指数在呼伦贝尔草原区的干旱监测与分析中具有较好的适用性。

4 讨论

事实上, 多个研究表明了气温在干旱变化中所起的重要作用^{[20, 49]}, 而我国干旱事件持续增加主要和气温大幅升高、降水变化不显著有关^[50]。其中, 中国北方的干旱化趋势与太平洋海温的年代际异常有关, 特别是与太平洋年代际振荡(PDO, Pacific Decadal Oscillation)存在显著的位相对应关系^[51-52]。在2000年以后, PDO由暖位相开始转换为一个冷位相, 北方降水趋于增多^[51]。然而本文基于SPEI指数, 对我国北方以呼伦贝尔草原为代表的草原生态系统的干旱特征的研究结果发现, 呼伦贝尔草原在气温显著上升, 降水减少的背景下干旱化趋势显著, 这虽然与我国北方气候变化的大背景是相似的^[8], 然而2000年以后该区域SPEI指数以及所发生的干旱事件均表明, 该区域干旱化趋势更加显著。一方面, 该区域降水量总量少, 年内分配不均, 其中, 夏季降水占全年降水总量的70.72%, 可能更易出现旱情^[53]。此外, 马柱国等^[50]基于GRACE(Gravity Recovery and Climate Experiment)卫星数据反演的陆地水储量(TWS, Terrestrial Water Storage)对中国区域干旱化问题的分析结果发现中国北方降水量从2001年以后有一个增加的趋势, 但陆地水储量TWS却是减小趋势, 说明我国北方干旱化仍在加剧。

值得说明的是, Hurst指数虽然在一定程度上可以预测呼伦贝尔草原未来的干旱趋势, 但并不能预测未来的干旱趋势所持续的时间。Morlet小波虽然可以判别时间序列中所包含多时间尺度周期性的大小及这些周期在时域中的分布, 但其结果往往存在着多重时间尺度上的复杂嵌套结构, 且在不同的尺度周期中, 表现出不同的干湿振荡规律, 造成结果的不确定性。本文将R/S与Morlet小波分析相结合并考虑时间序列的变化趋势, 有利于提高对呼伦贝尔草原未来尤其是近期干旱预测的准确性。

此外, 干旱具有发生的隐蔽性、形成的复杂性以及多尺度变化特征, 目前开发的任一干旱指标均很难达到时空上的普适性条件, 因为干旱指标大都建立在特定的时间和区域范围内^[30]。因此, 这给干旱的评估和预测带来了严峻挑战。虽然呼伦贝尔草原区过去58年发生的干旱事件均通过本文使用的SPEI方法得到了一致的验证, 在很大程度上表征了该区域的干旱情况, 但反过来, 仍存在SPEI指数较低的时期, 而事实上并未发生干旱事件的现象。这也恰恰说明了干旱事件的发生不仅取决于温度、降水和蒸散等物理性质, 还取决于生态、社会和经济因素^[54-55]。此外, 本文未能揭示该区域干旱发生的物理机制以及受限于气象站点的数量而忽略了空间异质性, 具有一定的局限性。因此, 未来可通过使用多种观测数据(包括物理的、社会的和经济的)和评估方法(如模型的、指标的)为特定区域和生态系统开发综合的干旱指数, 并分析不同维度对干旱发生的相对贡献, 以便进一步了解干旱发生的机制。

5 结论

1960—2017年, 呼伦贝尔草原的气温呈现显著上升趋势, 而降水呈不显著减少趋势。然而, 季节性的温度和降水变化趋势差异较大。多年来夏季降水呈减少趋势, 其他季节均呈不同速率的增加趋势, 其中以冬季增加的趋势最显著, 但冬季的气温并未发生显著变化。春、夏和秋季的年际温度变化均呈显著上升趋势, 其中春季气温上升速率最快, 达到0.46℃/10a, 夏季次之(0.39℃/10a)。

对不同年代干旱发生频次及季节分布的分析发现, 中旱发生频率在不同年代以及不同季节的分布均是最高的, 重旱次之, 极旱最少。与其他年代相比, 21世纪初(2000—2017), 干旱发生的总频率以及重旱和极旱的频率均显著高于其他年代。虽然不同季节发生干旱的总频次差异不大, 但不同等级干旱的季节分配差异较大, 其中, 冬季、夏季和秋季分别是发生中旱、重旱和极旱次数最多的季节。

通过Morlet小波和R/S分析对呼伦贝尔草原未来干旱趋势的分析发现, 该区域在11年尺度上呈现旱湿交替的周期性的可能性较大, 且未来几年有持续偏旱趋势。然而, 不同季节的周期性和趋势具有差异性, 其中, 春季在22年尺度上呈现旱湿交替的可能性较大, 未来有持续偏旱的趋势。夏季在9年尺度上呈现旱湿交替的可能性较大, 未来有持续偏旱的趋势。秋季在15年尺度上呈现旱湿交替的可能性较大, 未来有持续偏旱的趋势。冬季在17年尺度上呈现旱湿交替的可能性较大, 未来有持续偏湿的趋势。此外, SPEI指数与同期干旱事件发生的吻合度较高, 表明该指数在呼伦贝尔草原区的干旱监测与分析中具有较好的适用性。