淮河流域夏季旱涝前兆信号及预测

符芳兵， 薛联青， 任 磊

(河海大学 水文水资源学院， 江苏 南京 210098)

1 研究背景

全球气候变暖会导致气象灾害的增加，如干旱、洪水、极端气温等[1]，某些地区更容易出现干旱和洪涝，尤其是北半球中高纬度的环球灾害带[2-3]。旱涝影响范围广，对人类生活威胁性大，世界各地对其风险管理问题高度关注，而旱涝的识别和早期预警是其首要任务之一[4]。在我国，华东和华中地区的旱涝受灾程度分别位居第一和第二[5]，而淮河流域是我国南北气候的过渡带，在黄河夺淮的影响下，旱涝频繁[6]，易发生旱涝急转事件[7]。因此，对该流域开展旱涝发生规律的识别以及前兆信号的研究，对确保农业生产安全、减小气象灾害的影响以及保障人民经济健康发展，实现人与自然的和谐发展有重要的意义。

众多学者研究表明海温和大气环流与旱涝有密切的关系，他们已经对我国华北和南方地区旱涝的气候背景特征及前兆信号进行了分析研究。例如，彭瑞等[8]讨论了西江流域夏季严重旱涝的气候背景及前兆强信号；周晋红等[9]探讨了山西春季降水与500 hPa环流场及太平洋海温场的异常关系；吴瑶[10]通过波包传播诊断(WPD)、EEOF、相关系数等方法对西南地区连旱事件的前兆信号进行了研究；王春学等[11]采用多锥度-奇异值分解等方法，研究了四川盆地区域性暴雨的时空变化特征，并对其前兆信号进行了初步探讨。他们大部分都选取了典型旱涝年份，分析同期海温和大气环流异常变化剧烈的区域，作为旱涝的前兆信号，很少从时间尺度上考虑旱涝发生前期降水与海温、大气环流变化的相关关系，也很少建立前兆信号与旱涝关系模型，因此难以进行旱涝的预测。

本文以旱涝频发的淮河流域为研究对象，采用改进的Z指数法，识别流域夏季旱涝，划分灾害等级，分析其时空变化规律和周期特征；计算研究区域50年来降水与北太平洋海温、大气环流变化的相关系数，从时空角度找出旱涝前兆信号出现的时间和地点；以这些前兆信号作为预测因子，利用逐步回归法建立旱涝指标值与预测因子的数学模型，进行旱涝指标值的预测。

2 资料及方法

2.1 研究区域概况

淮河流域地处中国东部，介于东经111°55′～121°25′、北纬30°55′～36°36′之间，跨河南、湖北、安徽、江苏、山东五省，总面积约27×104km2。地形大体由西北向东南倾斜，流域西、南、东北部为山区和丘陵区，约占总面积的1/3，其余为平原、湖泊和洼地，约占2/3[12]。流域位于中国南北气候过渡带，属暖温带半湿润季风气候区，年平均气温在11-16℃，受特殊的地理位置和气候要素影响，旱涝频繁[7]。1949年至今，淮河流域在全流域或局部流域尺度上发生过洪涝灾害10 多次，多年发生严重干旱，给人民生命财产安全造成了巨大的损失[13]。通过对降水资料的分析对比，选取了研究区域28个气象站点，区域范围与气象站点分布见图1。

图1 淮河流域气象站点及水系分布图

2.2 数据来源

(1)气象资料来源于中国气象数据共享服务网提供的中国地面气候资料日值数据集(V3.0)，选取了淮河流域28个气象站1968-2017年夏季(6、7、8月)的逐日降水资料。为保证资料的可靠性，在进行计算之前对其进行处理，包括缺测值的补全、异常值的剔除等。

(2)海温和500 hPa位势高度数据来源于美国国家海洋和大气管理局。其中海温数据集(COBE SST2 and Sea-Ice)为全球范围1°×1°格点的海温，截取了北太平洋地区的数据。全球气压位势高度数据集(NCEP/NCAR再分析数据)精度为2.5°×2.5°，包括17个气压水平的位势高度，选取了欧亚大陆范围内的500 hPa气压对应的位势高度。

2.3 改进的旱涝评估方法

日降水量在时间上不服从正态分布，而采用服从PersonⅢ型正态分布的Z指数方法来拟合降水量，能得到较好的结果[14]。Z指数的大小可以反映旱涝的严重程度，但淮河流域范围较大，包含28个站点，考虑到旱涝的空间分布不均，同一时段各站的旱涝轻重程度不同，如果根据这些站点夏季降水平均值计算整个流域的旱涝等级，势必会影响评价结果的准确性。参考鞠笑生等[15]的研究成果，对Z指数法做出了改进。

普通的Z指数法计算公式与旱涝等级划分标准如下：

(1)

式中：Cs为偏态系数；Ji为第i年降水量标准化序列。其计算公式为：

(2)

(3)

表1 Z指数的旱涝等级划分标准

改进的Z指数法是先利用式(1)计算流域各个站点历年夏季降水的Z指数，再根据表1划分旱涝等级，然后统计出各个旱涝等级包含的站数，最后根据式(4)计算整个流域的旱涝指标，根据表2划分淮河流域的旱涝等级。

L=2n1+1.5n2+n3,G=n5+1.5n6+2n7

(4)

式中：L、G分别为雨涝指标和干旱指标;n1、n2、n3、n5、n6、n7分别为Z指数为1、2、3、5、6、7级的站数。系数2和1.5旱涝程度的权重，即重涝、重旱站数越多，则研究区域旱涝程度越严重。

表2 淮河流域旱涝等级划分

2.4 旱涝前兆信号识别与旱涝指数预测模型

根据降水与北太平洋海温、500 hPa大气压高度的相关关系，从海温场和大气环流场中寻找淮河流域夏季旱涝的前兆信号，将其作为旱涝指标的预测因子，利用逐步回归法建立淮河流域夏季旱涝指标值与预测因子的数学模型。1968-1997年为模拟期，1998-2017年为验证期。具体实施步骤如下：

(1)裁剪出北太平洋地区(120°E～80°W，10°S～70°N)1968-2017年1°×1°格点的海温数据和欧亚大陆范围(30°E～150°E，10°N～70°N)1968-2017年2.5°×2.5°格点的500 hPa气压对应的位势高度数据，这些数据为淮河流域夏季(6-8月)前的月尺度数据。筛选出流域1968-2017年夏季降水数据。

(2)根据公式(5)计算每个格点的海温或500 hPa气压位势高度与降水的Pearson相关系数，利用ArcGIS软件进行空间插值，生成相关系数分布图。

(5)

式中：变量x为降水量，mm；变量y为海温，℃，或者500 hPa气压位势高度，gpm。

(3)对相关系数进行显著性分析，挑选出通过99%、95%置信度的检验的区域作为淮河流域夏季旱涝的前兆信号关键区域。

(4)以这些前兆信号作为流域夏季旱涝指标值的预测因子，其中海温预测因子用H表示，500 hPa高度预测因子用D表示。建立旱涝指标值Y(对应于前文中的L-G值)与H、D的多元线性回归模型：

Y=β0+β1H1+β2H2+…+βmHm+α0+α1D1+

α2D2+…+αmDm+ε

(6)

式中：β0、β1、β2、…、βm和α0、α1、α2、…、αm为回归系数；ε是遵从正态分布(N，σ2)的残差。

将旱涝指标，海温和或500 hPa气压位势高度数据代入公式(6)，角标n和m分别代表年份和第m个预测因子，得到如下结构式：

(7)

求解出各项回归系数β和α以及残差ε后，对回归效果进行统计检验，回归效果符合要求时，就可以利用回归方程对流域降水旱涝指标值进行模拟预测。

3 结果与分析

3.1 旱涝的识别与特征分析

3.1.1 旱涝等级评价与发展趋势 利用上述改进的Z指数法计算了近50年来淮河流域降水的干旱指标和雨涝指标，并对其进行旱涝等级的划分(图2)。其中正常年份占总年份46%，雨涝年份占26%，干旱年份占28%，淮河流域旱涝严重。

从淮河流域夏季旱涝类型可以看出，1978、1988、1999和2014年为严重干旱年份，1991、2003、2005和2007年为严重洪涝年份。1978年，黄淮海流域发生特大旱，5月初仅山东、山西、河北三省受旱面积就达881.46×104hm2，其中山东500×104hm2，9月初，安徽、江西、江苏三省受旱面积336.66×104hm2，其中安徽农作物枯死16.66×104hm2。1988年6-7月中旬，鄂北、皖中、苏北地区发生重旱，湖北到7月底受旱农田面积达286.7×104hm2，100余万人饮水困难。1999年6-9月，长江以北地区出现罕见的大旱，武汉36万人饮水告急，安徽境内淮河干流几乎断流。2014年6-8月，黄淮等地区持续高温少雨，淮河上中游来水偏少7-8成，主汛期发生严重夏伏旱[18]。1991年5-7月梅雨期间，受暴雨和大风袭击，淮河流域冬小麦大面积倒伏，部分农田出现渍涝，仅河南、安徽、江苏等地受灾农作物就达200×104hm2。2003年，淮河流域发生了新中国成立以来仅次于1954年的大洪水，洪涝受灾面积达384.69×104hm2，成灾面积259.14×104hm2，绝收面积112.95×104hm2，受灾人口3 729.7×104人，因灾死亡29人，倒塌房屋75.7×104间，直接经济损失285.8×108元。2005年7月，淮河上游发生较大洪水，淮河支流洪汝河发生大洪水，淮北中南部及江淮北部大部分地区降雨量达200～400 mm，仅安徽省受灾人口达797.9×104人，农作物受灾面积80.8×104hm2，绝收19.0×104hm2，倒塌房屋2.14×104间，直接经济损失23.72×108元。2007年，我国气候异常，降雨分布不均，东部地区呈中间多、南北少的态势，淮河流域降雨偏多，发生流域性大洪水，沿淮河南、安徽、江苏三省，洪涝受灾面积257.70×104hm2，其中因涝农作物受灾面积236.70×104hm2，洪涝受灾面积的91.8%。

图2反映了淮河流域夏季旱涝指标及对应的旱涝等级随时间的变化。从图2中可以看出，淮河流域夏季旱指标呈下降趋势，下降速率为-0.44/10a，涝指标呈上升趋势，上升速率为0.03/10a。这也印证了全球气候变暖，降水量增多这一观点。虽然旱涝指标线性趋势变化很小，但涝灾年代间变化大，具有周期性，从80年代中期到90年代末，淮河流域夏季干旱发生频率高，90年代中期到2010年左右洪涝严重。

图2 淮河流域夏季旱涝指标及旱涝等级变化

3.1.2 旱涝的空间分布 由淮河流域旱灾频率空间分布图(图3(a))可知，淮河流域整体上干旱频率中部低，东西两边高，其中安徽北部地区干旱频率最低。以秦岭-淮河线为界，淮河以北地区干旱频率较高，17个站的干旱频率均值达25%，淮河以南地区干旱频率均值为23.8%。28个站点中有4个站干旱频率达30%以上，有11个站干旱频率达25%以上。山东的兖州地区和江苏的射阳地区干旱频率最高，达32%，河南的郑州、西华和江苏的高邮地区干旱频率最低，都在17%左右。

从淮河流域涝灾频率空间分布图(图3(b))可以看出，河南信阳地区和山东费县地区雨涝频率最高，分别为32%和30%，安徽亳州地区和江苏东台地区雨涝频率最低，分别为18%和16%。28个站中虽然只有2个站雨涝频率在30%以上，但在25%以上的有15个站。秦岭-淮河线以北地区平均雨涝频率为24%，淮河以南地区为25.3%。在河南省东南部，山东省西南部，江苏省北部形成了3个明显的雨涝高值区，而安徽地区各站的雨涝频率都比较均匀。

总体来说，淮河流域呈现出南涝北旱的特征，但旱涝频率差别不大。全流域平均雨涝频率大于干旱频率，分别为24.6%和24.3%。河南信阳地区干旱和雨涝频率都很高，说明该地区降水量年际间变化相差大，应该特别注意旱涝和用水安全。河南西华地区旱涝频率均比较低。

3.2 旱涝前期海温、大气环流信号

3.2.1 北太平洋海温关键区和关键时段 从淮河流域夏季降水与前期各月北太平洋海温的相关关系出发，从空间尺度上寻找出海温关键区，从时间尺度上寻找出海温关键时段。图4给出了淮河流域夏季降水与前期(上一年12月-本年5月)各月海温场的相关系数分布。本文原本从上一年9月份开始研究淮河流域夏季降水与海温的关系，但上一年12月之前的相关系数不显著，图略。

图3 淮河流域旱涝频率空间分布

图4 淮河流域夏季旱涝与海温场的相关系数分布(红色和绿色分别代表通过99%、95%置信度的检验)

总体来说从上一年12月到本年5月，北太平洋海温与淮河流域夏季降水相关性越来越高，相关区域也越来越大。上一年12月，赤道西南太平洋出现了显著的负相关区域，千岛群岛附近出现正相关区域；1月开始，千岛群岛附近的显著性增加，正相关区域面积也增大；2月开始，南部赤道附近的Nino4区出现显著的正相关区域；3月份，千岛群岛附近，中太平洋地区和Nino4区的显著性达到最大，相关区域面积也达到最大；4月到5月中太平洋地区的负相关区域逐渐向西南移动，到达菲律宾海附近。海温关键区域大致都分布于180°经线以西的西太平洋地区。从北部的千岛群岛附近，到中部的菲律宾海附近，到南部赤道附近的Nino4区，相关系数呈“+-+”的分布结构。在3、4和5月份这一规律更加明显，这个3个月相关程度最高，相关区域也最大，因此以这3个月作为海温关键时段。

3.2.2 500 hPa高度场关键区和关键时段 为了从大气环流方面找出影响淮河流域旱涝发生的关键区域和关键时段，与海温相同，研究了夏季降水与前期(上一年12月-本年5月)各月500 hPa高度场的相关性。图6给出了欧亚大陆范围内各月500 hPa高度场与淮河流域夏季降水的相关系数分布。

由图5可知，上一年12月份，东欧平原地区的500 hPa高度场与淮河流域夏季降水为显著的负相关关系，西西伯利亚一带出现显著的正相关区域；1月份，西亚地区出现显著的正相关区域；2月份开始，东西伯利亚地区变成了正相关区域；3月份，西北太平洋地区出现显著的负相关区域；4月份，东欧平原地区的负相关区域已变成显著的正相关区域，西亚地区的正相关区域也变成了负相关区域，而东西伯利亚一带的正相关区域依旧显著；5月份，西亚地区的负相关区域向东移动，而西亚以西的地中海地区出现显著的正相关区域，东西伯利亚一带的正相关区域显著性达到最强。从相关关系的移动方向来看，可以总结一个规律：从上一年12月份开始，东欧平原500 hPa高度场的负相关区域逐渐向东南方向移动，东西伯利亚一带的正相关区域逐渐向西移动，且越来越显著。

图5 淮河流域夏季旱涝与500 hPa高度场的相关系数分布(红色和绿色分别代表通过99%、95%置信度的检验)

对图5的分析可知影响淮河流域夏季旱涝的500 hPa高度场关键区域大致分布于东西伯利亚山地、东欧和西西伯利亚平原以及西亚地区。旱涝与前期各月500 hPa高度场都有一定的相关性，但在3-5月份相关系数分布相似，也比较稳定，与海温相同，以这3个月为500 hPa高度场关键时段。

3.3 旱涝的预测模型及效果检验

3.3.1 预测因子的选取结果 前一小节分析了影响淮河流域夏季降水的海温和大气环流方面的前兆信号，得到8个海温和500 hPa气压高度关键区域，以其作为流域夏季降水预测因子，前期关键时段(3、4和5月)各月的预测因子不同，见下表3。

表3 旱涝的预测因子

3.3.2 旱涝的逐步回归预测效果 利用选取的关键时段(3、4和5月)各月的预测因子，使用SPSS统计软件建立淮河流域夏季旱涝指标值(L-G)预测线性模型。其模型方程如下：

3月：

Y1=0.541H1-6.441H2+18.619H3+

723.265D1+485.982D2-6792.057

(8)

4月：

Y2=1.115H1-64.806H4+37.276H5-

6448.784D1+5496.505D5+4260.61

(9)

5月：

Y3=3.206H1-18.555H4+3.155H5-

232.083D1+9583.62D5-53416.032

(10)

其中海温的单位为℃，500 hPa气压高度的单位为103gpm。

根据模型方程模拟出淮河流域夏季降水旱涝指数，与实际计算值进行对比发现(图6、7)，模拟期和验证期内4、5月份的模型方程较3月份模拟效果好，能够拟合出绝大部分极值。由F检验临界值表查出，F0.01(5,30)=3.70，3月、4月、5月模型方程的F检验值分别是F1=4.2、F2=9.6、F3=7.8，均大于F0.01。所以建立的预测方程是显著的，方程式是可行的。用复相关系数检验模型方程的优劣，3个模型方程的在模拟期和验证期的复相关系数(表4)均在0.7以上，表明模型的适用性也较好。

综上所述，淮河流域夏季降水的旱涝指数可以用前期3、4和5月3个模型方程及其预测因子进行预测，而4、5月份的模型方程最为准确。

图6 模拟期淮河流域1968-1997年夏季降水

图7 验证期淮河流域1998-2017年夏季降水旱涝指数模拟效果图

表4 回归方程的效果检验

4 结论与讨论

利用改进的Z指数法对淮河流域近50年来旱涝等级进行了评估，根据逐步回归法建立了旱涝指标值的预测模型，得出以下主要结论：

(1)旱涝等级评价结果中的严重旱涝年份与淮河流域历史旱涝年份是对应的，说明评估方法较为准确。

(2)淮河流域旱涝灾害年份占总年份的54%，旱涝灾情严重，呈现出南涝北旱的特征，河南信阳地区干旱和雨涝频率都很高，河南西华地区旱涝频率均比较低。流域夏季旱指标下降速率为-0.44/10a，涝指标上升速率为0.03/10a，周期变化以1990年为界限。

(3)影响淮河流域夏季旱涝的海温关键区大致都分布于180°经线以西的西太平洋地区，从南到北分别为千岛群岛，菲律宾海到太平洋中部，Nino4区。500 hPa高度场关键区大致分布于东欧平原和西西伯利亚平原、东西伯利亚山地以及西亚地区。3、4和5月份是海温和500 hPa高度场信号的关键时段。

(4)利用逐步回归法建立的方程，具有较好的预测精度和实用性，能够对淮河流域夏季降水的旱涝指数进行预测，其预测结果将在今后的实际工作中得到进一步检验。