干旱是中国主要的气象灾害之一, 近年来, 在全球变暖的背景下^[1], 全国遭受干旱的地区在扩大, 干旱造成的损失大幅增加^[2-3], 且未来几十年旱情可能呈加重趋势^[4-5]。长江中下游流域人口密集、经济发达, 是我国重要的粮、油、棉生产基地, 该区虽水资源丰富, 但年内降水时间分布不均且年际变化较大, 区域内季节性干旱时有发生, 造成的影响不容忽视^[6-8], 如2011年春旱给该区农业、人畜饮水、江河及湖泊水位、渔业生产、水运等带来了严重影响^[9]。而且相关研究表明, “骤发性干旱”在我国南方地区发生频率更高^[4], 南方干旱会对生态系统^{[8, 10-11]}产生严重的影响。

干旱指数是研究干旱的重要手段, 黄晚华等利用标准化降水指数(SPI)^[12]和降水距平百分率^[13]对中国南方季节性干旱进行了分析, 王文等^[6]用Palmer指数(PDSI)对长江中下游地区干旱特征进行了研究。SPEI综合考虑了气温和降水对干旱的影响, 弥补了降水距平百分率和SPI未考虑气温对干旱影响的缺点, 同时与PDSI相比, 具有计算简单、多时间尺度、多空间比较等优点^[14-16], 在我国湿润地区有很好的适用性^[17-18]。马彬等^[19]基于SPEI对中国东部季风区的研究表明, 气候出现暖干化, 局部地区干旱出现极端化。李亮^[20]基于综合气象干旱指数和SPEI对长江中下游地区季节性干旱特征进行了分析, 但仅仅是基于代表站点。张余庆等^[15]基于SPEI分析了赣江流域旱涝演变的周期性特征。赵林等^[21]用SPEI分析了湖北省年尺度干旱站次比、强度和频率。以上研究多在大区域或省域尺度开展, 流域尺度干旱变化特征的研究对于流域水资源管理和生态平衡具有重要意义^[22], 但是基于SPEI对长江中下游流域不同时间尺度干旱站次比、强度、频率的综合分析相对较少。考虑到干旱对长江中下游流域自然生态环境及社会经济影响的严重性, 本文利用SPEI指数, 结合干旱站次比、强度及频率, 分析长江中下游流域1961—2015年年尺度和各季节干旱时空变化特征, 并对干旱和流域气温、降水变化及ENSO的关系进行了研究, 可以进一步丰富人们对该区域干旱发生规律的认识, 以期为流域水资源规划和防灾减灾提供科学依据。

1 资料来源与研究方法 1.1 研究区概况

长江是中国最长的河流, 上游和中游的分界点为湖北宜昌, 中游和下游的分界点为江西湖口。长江中下游流域由中游干流区(主要分布于湖北省)、汉江流域(主要分布于陕西省、河南省及湖北省)、洞庭湖流域(主要分布于贵州省、湖南省和湖北省)、鄱阳湖流域(主要分布于江西省)、长江下游干流区(主要分布于安徽省和江苏省)及太湖流域(主要分布于江苏省、上海市及浙江省)组成^[23](图 1)。长江中下游流域多为海拔较低的丘陵和平原, 以亚热带季风气候为主, 东亚季风活动明显^[7]。

1.2 数据来源

气象数据来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn), 选用长江中下游流域1961—2015年资料序列较长的129个气象站点的逐日降水和气温数据, 经过进一步计算得到月降水量和月平均气温。该数据经过了严格的质量控制, 其中有6个站点存在数据缺测, 缺测时间主要为1967和1968年的部分月份, 用相邻站点线性回归方法对缺测数据进行插补, 插补后的数据经过极值检验和时间一致性检验。湖北省的五峰站迁站海拔差异较大, 故未使用该站数据, 气象站点分布情况如图 1所示。多变量ENSO指数(MEI)数据来源于http://www.esrl.noaa.gov/。Niño3.4区海洋表面温度距平(SSTA)和南方涛动指数(SOI)数据来源于http://www.cdc.noaa.gov/。

1.3 研究方法 1.3.1 标准化降水蒸散指数

标准化降水蒸散指数(SPEI)是Vicente-Serrano等在标准化降水指数(SPI)的基础上, 考虑水分亏缺和累积效应两个因素, 用降水量和潜在蒸散量的差值偏离平均状态的程度来表征干旱, 其中潜在蒸散量的算法主要为Thornthwaite法和Penman-Monteith法, 因为基于两种算法的SPEI值在长江中下游流域差异较小^[24], 且Thornthwaite法具有计算简便的优点, 本文采用该方法计算潜在蒸散量, SPEI具体计算方法参照文献^[25]。参考相关研究^[26], 根据累积概率分布函数将SPEI值分为9个旱涝等级：极涝(SPEI≥2.0)、重涝(1.5≤SPEI＜2.0)、中涝(1.0≤SPEI＜1.5)、轻涝(0.5≤SPEI＜1.0)、正常(-0.5＜SPEI＜0.5)、轻旱(-1＜SPEI≤-0.5)、中旱(-1.5＜SPEI≤-1.0)、重旱(-2.0＜SPEI≤-1.5)、极旱(SPEI≤-2.0)。3个月时间尺度SPEI能反映季节尺度的旱涝情况, 与农业旱涝关系密切；12个月时间尺度SPEI能较清晰的反映长期旱涝变化特征, 并对河流径流量、下层土壤含水量以及水库储水量有较好的反映^[15]。本文选用3、12个月时间尺度的SPEI分析研究区季节和年尺度干旱时空演变特征。春、夏、秋、冬季干旱分别由3个月时间尺度5月、8月、11月和次年2月的SPEI值表示, 年尺度干旱由12个月时间尺度12月份的SPEI值表示。

1.3.2 干旱站次比(drought stations proportion, P_j)

用某一区域内干旱发生站数占全部站数的比例来评价干旱影响范围的大小^[12]。

式中, M是研究区总气象站数；m为发生干旱的站数；下标j为不同年份的代号。

1.3.3 干旱强度(drought severity, S_ij)

用于评价某时段区域干旱严重程度。

式中, j表示年份, i表示不同站号, m为发生干旱的站数, |SPEI_i|表示j年i站发生干旱时SPEI的绝对值。

1.3.4 干旱频率(drought frequency, F_i)

用于评价某站在某时段干旱发生的频繁程度^[12]。

式中, 下标i为不同站号, N为某站有气象资料的总年数；n为该站发生干旱的年数。

由于长江中下游流域水资源较为丰富, 轻旱造成的影响相对较小, 本文中的干旱站次比、强度和频率均按中旱(含中旱及以上)等级统计。采用线性倾向估计法分析SPEI及干旱站次比与强度的年际变化趋势^[27]，显著性检验为F检验。采用反距离加权插值法对SPEI变化趋势和干旱频率进行空间化处理。用偏相关分析法研究SPEI和气温、降水的关系^[28], 用相关分析对干旱和ENSO各表征因子的关系进行研究^[28]。

2 结果分析 2.1 干旱时间变化特征 2.1.1 SPEI变化趋势

由表 1可知, 1961—2015年长江中下游流域年尺度SPEI呈下降趋势, 倾向率为-0.020/10a。春季和秋季SPEI也均呈下降趋势, 即干旱化趋势, 倾向率分别为-0.121/10a和-0.093/10a, 春季干旱化趋势显著(P < 0.05)。夏季和冬季SPEI呈上升趋势, 即湿润化趋势, 倾向率分别为0.077/10a和0.049/10a。

2.1.2 干旱站次比和干旱强度

为了进一步说明干旱变化趋势, 本文对1961—2015年长江中下游流域干旱站次比和强度进行了分析。年尺度干旱站次比和强度均呈增加趋势, 倾向率分别为1.065%/10a和2.15/10a, 表明年尺度干旱呈加重趋势。年干旱站次比在50%以上的年份有1966、1978和2011, 1978年最高, 达68.22%。年干旱强度最大的年份为1978年, 达150.64, 其次为2011(118.77)和1966(100.74)(图 2)。

春旱站次比为极显著(P < 0.01)的增加趋势, 倾向率为3.969%/10a, 春旱强度为显著(P < 0.05)的增加趋势, 倾向率为8.785/10a。春旱站次比在50%以上的年份有2000、2001、2007和2011年, 2011年最高, 达87.60%。春旱强度最大的年份为2011年, 达235.61, 其次为2007(125.04)和2001(105.29)(图 2)。

夏旱站次比和强度均呈现减小趋势, 倾向率分别为-2.145%/10a和-4.479/10a, 表明夏旱呈减弱趋势。夏旱站次比在50%以上的年份为1978年, 达61.24%。夏旱强度最大的年份为1978年, 达137.12, 其次为1972(104.21)和2013(90.57)(图 2)。

秋旱站次比和强度均呈增加趋势, 倾向率分别为2.934%/10a和6.004/10, 表明秋旱呈加重趋势。秋旱站次比在50%以上的年份有1979、1998和2007, 1998和2007年均达到63.57%。秋旱强度最大的年份为1998年, 达149.64, 其次为2007(121.76)和1979(104.37)(图 2)。

冬旱站次比和强度均呈下降趋势, 倾向率分别为-1.365%/10a和-1.939/10a, 表明冬旱呈减轻趋势。冬旱站次比在50%以上的年份有1962、1967和1998, 1998年最大, 达91.47%。冬旱强度最大的年份为1998年, 达230.54, 其次为1962(158.71)和1967(111.47)。年和各季节干旱站次比和强度的变化趋势与干湿变化趋势相对应, 由干旱站次比与强度确定的重旱年与相关文献的记录相吻合^{[3, 7, 20]}(图 2)。

2.2 干旱空间变化特征 2.2.1 SPEI变化趋势

图 3是1961—2015年长江中下游流域年尺度及各季节SPEI变化趋势的空间分布图。1961—2015年, 长江中下游流域年尺度SPEI呈上升、不变和下降趋势的站点分别占总站点数的44.96%、1.55%和53.49%。区域西部的洞庭湖流域、中游干流区及汉江流域SPEI主要为下降趋势, 倾向率集中在-0.1/10a—0之间；东部的鄱阳湖流域、下游干流区及太湖流域SPEI主要为上升趋势, 倾向率集中在0—0.1/10a之间。

96.90%的站点春季SPEI呈下降趋势, 表明春季长江中下游流域整体以干旱化趋势为主, SPEI下降幅度空间差异明显, 在鄱阳湖流域东南部较小, 向外围逐渐增加。26.36%的站点春季SPEI呈显著(P < 0.05)下降趋势, 主要位于太湖流域、下游干流区东部、洞庭湖流域中西部及汉江流域, SPEI变化倾向率多小于-0.15/10a。

82.95%的站点夏季SPEI呈上升趋势, 表明夏季以湿润化趋势为主。SPEI上升幅度空间差异明显, 上升幅度最高的区域在太湖流域, 在0.15/10a以上, 其次为下游干流区、鄱阳湖流域东部及汉江流域中部, 在0.1/10a—0.15/10a之间。10.08%的站点夏季SPEI显著(P < 0.05)上升, 在太湖流域分布集中。

92.25%的站点秋季SPEI呈下降趋势, 表明秋季以干旱化趋势为主。SPEI下降幅度在鄱阳湖流域中部和北部及洞庭湖流域北部较低, 向南北两侧逐渐升高, 在汉江流域东部、下游干流区东部及太湖流域下降幅度较大, 低于-0.15/10a。13.95%的站点秋季SPEI呈显著(P < 0.05)下降趋势, 主要分布于太湖流域、下游干流区东部及汉江流域东部。

分别有72.87%和27.13%的站点冬季SPEI呈上升和下降趋势, 表明以湿润化趋势为主。呈下降趋势的站点集中分布于洞庭湖流域中西部及汉江流域, 其他区域以上升趋势为主。11.63%的站点SPEI上升趋势显著, 主要分布于太湖流域和下游干流区, 上升幅度多大于0.15/10a。值得注意的是, 太湖流域大部分站点在春季和秋季呈显著(P < 0.05)的干旱化趋势, 而在夏季和冬季呈显著(P < 0.05)的湿润化趋势。洞庭湖流域中西部和汉江流域在春、秋、冬和年尺度均以干旱化趋势为主。

2.2.2 干旱频率变化特征

1961—2015年和各时间段干旱频率如图 4所示^{[2, 21]}。

1961—2015年, 长江中下游流域年尺度、春、夏、秋、冬干旱频率分别为17.56%、16.27%、17.24%、17.56%和16.53%, 年尺度干旱和秋旱频率较高, 其次为夏旱。年尺度干旱频率在18%以上的区域主要在汉江流域、洞庭湖流域西部和东部、鄱阳湖流域西部和东部、中游干流区中部。春旱频率在18%以上的区域主要为汉江流域中部和洞庭湖流域西部。夏旱频率在18%以上的区域主要分布在汉江流域、中游干流区、洞庭湖流域南部和鄱阳湖流域南部。秋旱频率在18%以上的区域主要为洞庭湖流域西北和东部、鄱阳湖流域中北部和下游干流区。冬旱发生频率在18%以上的区域主要为汉江流域东部、中游干流区东部、鄱阳湖流域西北部、下游干流区及太湖流域南部。

不同年代干旱频率差异显著。年尺度干旱频率年代际变化为减-减-增-增, 在21世纪初达到最高, 区域平均干旱频率为24.40%。春旱频率年代际变化为减-增-增-增, 在21世纪初达到最高, 区域平均干旱频率为27.05%。夏旱频率在1960s最高, 年代际变化为减-减-减-增, 在1990s达到最低, 21世纪初有所增加。秋旱频率年代际变化为增-减-增-减, 在1990s达到最高, 为27.11%, 在21世纪初下降为24.22%。冬旱频率在1960s最高, 年代际变化为减-增-增-减。年尺度、春季和夏季干旱频率从1990s到21世纪初均表现为增加趋势, 秋旱频率虽未增加, 但仍高于1960s、1970s和1980s。

2.3 干旱变化影响因素分析 2.3.1 SPEI和气温、降水的关系

由于各季节SPEI变化趋势具有较好的空间一致性, 从区域角度对各季节SPEI变化和气温、降水的关系进行分析。由表 2可知, 1961—2015年长江中下游流域各季节气温均呈增加趋势, 除了夏季, 均通过了0.01的置信度检验, 春季气温上升最快, 为0.27℃/10a。夏季和冬季降水量均呈增加趋势, 倾向率分别为16.26 mm/10a和5.40 mm/10a, 夏季降水增加显著(P < 0.05)；春季和秋季降水量呈减少趋势, 倾向率分别为-7.62mm/10a、-5.08mm/10a。整体上SPEI和气温呈极显著(P < 0.01)的负相关, 和降水呈极显著(P < 0.01)的正相关。春季和秋季降水量的减少和气温的上升共同导致春、秋季呈干旱化趋势；夏季气温变化不明显、但降水量显著增加, 使夏季呈湿润化趋势；冬季气温和降水均呈增加趋势, 但冬季气温对SPEI的影响较小, 相关性仅为-0.60, 所以冬季呈湿润化趋势。

2.3.2 干旱和ENSO的关系

厄尔尼诺/南方涛动(ENSO)与中国各地的干旱有密切的联系^[29], 其特征值常用赤道太平洋中东部海洋表面温度距平值(SSTA)、南方涛动指数(SOI)和多变量ENSO指数(MEI)表示。当发生厄尔尼诺事件(ENSO暖事件)时, SSTA和MEI为正值, SOI为负值；发生拉尼娜事件(ENSO冷事件)时, SSTA和MEI为负值, SOI为正值^[30]。

对干旱和ENSO各特征值的相关性进行统计(表 3), 表明冬季ENSO对次年春旱影响显著, 相对于MEI和SSTA, 冬季SOI和次年春季SPEI、干旱站次比及强度相关性更高, 相关性极显著(P＜0.01)。冬季SOI和次年春季SPEI呈负相关, 与干旱站次比和强度呈正相关, 表明冬季SOI值越大, 越易发生干旱, 即拉尼娜事件在冬季发生时, 次年春季更易发生干旱。1961—2015年, 共有14年冬季发生拉尼娜事件^[31], 次年春季平均干旱频率、强度和站次比分别为24.35%、51.54和24.36%；而冬季未发生拉尼娜事件时, 次年春季平均干旱频率、强度和站次比分别为13.37%、24.50和13.50%, 由此进一步说明拉尼娜事件在冬季发生时对次年春旱的影响。春季MEI及SSTA与同年夏季SPEI相关性显著(P＜0.05), 但是和干旱站次比及强度并不存在显著的关系。

3 结论与讨论 3.1 主要结论

(1) 区域尺度上, 1961—2015年, 长江中下游流域年尺度、春季和秋季均呈干旱化趋势, 春季干旱化趋势显著；夏季和冬季均呈湿润化趋势。年尺度、春季和秋季干旱站次比及强度均呈增加趋势, 春旱站次比与强度增加趋势显著；夏季和冬季干旱站次比和强度均呈下降趋势。

(2) 空间变化上, 对于年尺度, 汉江流域、中游干流区及洞庭湖流域以干旱化趋势为主；鄱阳湖流域、下游干流区和太湖流域以湿润化趋势为主。对于季节尺度, 春季和秋季分别有96.90%和92.25%的站点呈干旱化趋势；夏季和冬季分别有82.95%和72.87%的站点呈湿润化趋势。

(3) 干旱频率时空分布差异显著, 年尺度和春季干旱频率在21世纪初均达到最高；年尺度、春季和夏季干旱频率从20世纪90年代到21世纪初均呈增加趋势。

(4) 春、秋季干旱化趋势与降水量的减少及气温的上升相关, 夏、冬季降水量的增加使得夏、冬季呈湿润化趋势。相对于MEI和SSTA, 冬季SOI和次年春旱相关性更强, 相关性极显著, 冬季发生拉尼娜事件时, 次年春季更易发生干旱。

3.2 讨论

本文综合分析SPEI、干旱站次比和强度的变化特征, 表明长江中下游流域年尺度、春季和秋季干旱加重；夏季和冬季干旱减弱。黄晚华等^[12]基于SPI对中国南方干旱的研究以及王文举等^[32]基于SPEI对湖北省的研究同样发现, 春旱和秋旱加重, 夏旱和冬旱减轻。不同学者分别对中国^[2]和中国东部季风区^[19]的研究发现, 近20年干旱事件增加；有学者^{[21, 32]}对湖北省的研究表明, 2000年以后干旱加重；本研究发现, 相对于之前的年代, 长江中下游流域年尺度和春季干旱频率在21世纪初均达到最高, 在干旱频率较高和干旱化趋势显著的区域应注重旱灾的防御。本文分析的气象干旱是其他类型干旱研究及干旱风险评估和区划的基础, 干旱变化趋势对于农业、生态、社会经济等的影响仍值得探讨^{[8, 33-35]}。

气温在干旱变化中起着重要的作用^{[2, 16]}, Chen等^[5]发现, 近20年中国干旱事件持续而显著的增加主要和气温大幅升高、降水却没有较大变化有关。本研究表明, 长江中下游流域春、秋季干旱化趋势受气温上升和降水减少共同影响。当冬季发生拉尼娜事件时, 次年春季更易发生干旱, 这是由于拉尼娜事件会使西太平洋副热带高压势力减弱, 暖湿气流无法深入长江中下游流域, 从而造成干旱^{[9, 36]}。本文仅从气温、降水及ENSO的角度对干旱变化的影响因素进行了分析, 机理方面有待进一步研究。