引言
寒露风害即“秋季低温寒害”，是南方双季稻地区。寒露风在安徽省称为秋分寒，是安徽省危害水稻的三大低温冷害之一。寒露风引起的低温主要影响双季晚稻的抽穗扬花，严重时花而不实或抽不出穗，极易导致空壳、瘪粒而减产。
长期以来，中国不少学者对寒露风的时空分布特征、寒露风对晚稻的影响与预防措施、寒露风预测等进行了研究。在寒露风的时空分布特征方面，刘文英等分析了1959—2008年江西省寒露风发生情况和演变趋势以及寒露风出现日数与空壳率的关系；黄晚华等系统地研究了1961—2008年湖南省寒露风的空间分布特征和出现规律，并应用GIS和定量评估技术给出了湖南省寒露风综合气象风险区划；刘丽英等统计分析了广东寒露风时空分布特征；李艳兰等和戴剑波等分别探讨了广西省和桂北的寒露风气候变化特征。在寒露风预测方面，中央气象局研究所早在20 世纪 70 年代就提出了长江中下游寒露风出现早晚趋势的长期预报指标；娄秀荣建立了长江中下游寒露风灰色系统预测模型；安晓明用因子搭配法制作桂西北寒露风长期预报；居晶琳建立湖南桃源县寒露风灰色预测模型；欧阳兆云建立广西德保县寒露风初日的模糊数学预测模型；肖富明对影响广西寒露风的多种因子进行普查、分析和挑选后，建立了广西寒露风气候趋势预测概念模型。辛学飞等研究给出湘西北寒露风偏早年份的主要环流特征。可见2000年以后，由于气候变暖，寒露风灾害明显减弱，对寒露风预测的研究工作逐渐减少。
总的来看，有关寒露风时空分布和气候趋势预测的研究主要集中在江西、湖南、两广地区，而针对安徽省的相关研究十分缺乏，系统的研究文献几乎没有。
因此，本研究旨在理解安徽省寒露风初日的时空分布规律及其与9月平均气温的联系，构建适用于安徽省的寒露风初日趋势预测模型，以提高安徽省寒露风气候预测专项服务水平，以期为政府和双季晚稻的生产部门提供决策依据，对预防和减少农业气象灾害具有重要的现实意义。
1 资料与方法
1.1 资料来源
考虑到在安徽寒露风对双季晚稻危害严重，安徽双季晚稻主要分布在沿江和江南，因此本研究范围为安徽省沿江和江南地区。使用的数据为1961—2013年安徽省沿江江南除天柱山和黄山2个高山站外的37个气象观测站的逐日平均气温数据集，数据来源于安徽省气候中心，经安徽省气象信息中心审核和整编，包括对资料进行要素允许值范围检查、极值检查、气候学界限值检查、时间一致性检查、内部一致性检查和空间一致性检查的质量控制。寒露风初日气候趋势预测模型所用的数据来源于国家气候中心的126项环流指数。
1.2 研究方法
寒露风指标采用安徽农业气象业务服务中使用的指标，即秋季日平均气温连续3天<20℃，该标准与气象行业标准《寒露风等级》中的中度干冷型寒露风等级和气象行业标准《水稻冷害评估技术规范》中的寒露风轻度灾害致灾指标一致。因此，寒露风初日指秋季日平均气温连续3天<20℃的初日。在寒露风气候趋势预测中，主要考虑寒露风初日，而寒露风初日的预测重点是寒露风偏早的研究和预报。
本研究主要分析方法包括：利用线性回归分析方法，探索寒露风初日的长期变化趋势；利用线性相关分析方法，研究寒露风初日和9月气温之间的联系及其与126项环流指数之间的联系。利用多元逐步回归分析方法，筛选有用的环流指数并建立寒露风初日和环流指数之间的线性回归预测模型。
2 结果分析
2.1 寒露风初日的空间分布特征
首次寒露风天气过程出现的早晚可以用寒露风初日出现的日期来表示。为计算方便，本研究以9月1日为日序1，依次顺推10月1日的日序为31。首先，给出不同气候平均态（1961—1990年和1981—2010年）下寒露风初日日序平均值的空间分布（见图1），可以看到，1961—1990年沿江东部寒露风初日平均值出现在9 月 26—27 日，中部在 9 月 27—28 日，而沿江西部在 9月29日—10月2日。江南寒露风初日平均值相对沿江要早4~7天，出现在9月22—25日，其中太平和宁国寒露风初日平均值出现的最早，在9月22—23日。
1981—2010 年与 1961—1990 年气候平均态下寒露风初日平均值的空间分布型基本一致，两者的空间相关系数达到0.83，但出现日期明显推迟2~3天。寒露风初日平均值在沿江东部和中部出现于9月30日—10 月 1 日，在沿江西部出现在 10 月 3 日，在江南出现在9 月 25—29 日，其中太平和旌德出现的最早，出现在24—25 日。【图1】
寒露风初日平均值的空间分布型与9月平均气温的空间分布型（图略）十分一致。在1961—1990年和1981—2010 年气候平均态上，寒露风初日和 9 月平均气温之间的空间相关系数分别达到0.91和0.92，均通过0.001信度的显着性检验。9月平均气温的高低与寒露风初日平均值的早晚关系密切，沿江平均气温高于江南，尤其沿江西部最高，对应着寒露风初日平均值沿江晚于江南，且沿江西部最晚。江南太平、旌德、宁国和泾县温度最低，对应着寒露风初日平均值出现的最早。1981—2010年相比1961—2010年9月平均气温明显升高，对应着寒露风初日平均值明显推迟。
2.2 寒露风初日的时间变化特征
从1961—2013年寒露风初日变化趋势的空间分布（见图2）来看，寒露风初日日序只有沿江西部呈现-0.4~0.8 d/10 a 的下降趋势，大部分地区为上升趋势。【图2】
表明沿江江南大部地区寒露风初日的出现均为推迟趋势，而沿江西部寒露风初日的出现有提前趋势。沿江东部寒露风初日日序为显着的上升趋势，变化趋势达到1~1.7 d/10 a，通过0.05信度的显着性检验，当涂的变化趋势最大，达到1.7 d/10 a。江南寒露风初日日序为明显的上升趋势，变化趋势为0.4~1 d/10 a，但未通过0.05信度的显着性检验，其中江南东部和南部上升趋势多在0.8~1 d/10 a，趋势大于江南北部。表明寒露风初日推迟趋势在沿江东部最为显着，其次是江南东部和南部，再次是江南北部。
这与9月平均气温变化趋势密切相关，两者的空间相关系数达到0.82。9月平均气温变暖趋势对应着寒露风初日日序的上升趋势，也就是寒露风初日的推迟趋势。沿江东部气温显着变暖趋势对应着寒露风初日显着推迟趋势，而沿江西部气温变化趋势微弱，对应着寒露风初日微弱提前趋势。江南气温变暖趋势南部大于北部，对应寒露风初日推迟趋势南部较北部明显。
2.3 区域寒露风初日的时间变化特征
考虑到寒露风出现早晚的区域一致性较强，气候趋势预测业务中用区域寒露风初日代表沿江江南寒露风出现的早晚。借鉴气象行业标准《寒露风等级》中区域寒露风的定义，定义沿江江南区域内出现寒露风的站点达到该区域总站数35%以上的寒露风初日为区域寒露风初日。1961—2013年区域寒露风初日与各站寒露风初日的相关关系在 0.42~0.93，多数在 0.6~0.8，均通过 0.001 信度的显着性检验，表明区域寒露风可以代表各站的时间变化特征。从1961—2013年区域寒露风初日出现早晚的年际变化（见图3）可看出，区域寒露风初日最早出现在9月10日（1981年、1984年和2006年），最晚出现在10月14日（1998年），早晚相差1个多月（34天），多年平均值在9月26日，寒露风初日出现在9月的年份有36年，占总年数的68%。由于安徽省晚稻抽穗扬花主要发生在9月中下旬，也就意味着安徽双季晚稻3年中有2年会受到寒露风的影响。1961—2013年区域寒露风初日日序呈现0.6 d/10 a的上升趋势，通过0.001信度的显着性检验，表明寒露风初日存在显着的推迟趋势。按年代际统计来看，区域寒露风初日在20世纪80年代和70年代出现的较早，而在 20世纪90年代和21世纪前10年出现的较晚。值得注意的是21世纪以来，相较于20世纪90年代，寒露风初日明显提前，2006年寒露风初日异常偏早位居历史并列第1，2010—2013年连续4年寒露风初日早于气候平均值，表明在气候变暖的背景下，区域寒露风初日整体趋于推迟，但仍会出现异常偏早的情况。【图3】
2.4 区域寒露风初日的气候趋势预测
在气候预测业务中，8月开展寒露风出现早晚的气候趋势预测专项服务。实际业务中通常将寒露风趋势的预测转化为9月平均气温趋势的预测，9月平均气温偏低时，寒露风初日偏早，反之亦然。寒露风初日与9 月平均气温之间的密切联系可从上文得到证实。该方法简单省事，但预测准确率受局限，1961—2013年寒露风初日和9月平均气温的距平符号有39年一致，仍有14年不一致。另外，在气候变暖背景下，气候预测模式过多考虑这一背景，一直预测气温偏高趋势。
如国家气候中心业务使用的月动力延伸集合模式在2005—2013 年均预测 9 月气温偏高，虽有 5 年趋势预测正确，但并不能准确预测气温偏低的年份，而气温偏低的年份往往正是寒露风出现早、危害大的年份。鉴于此，本研究构建基于前期物理因子的区域寒露风初日趋势预测模型。
寒露风初日气候趋势预测模型的建立过程主要有3 个步骤：首先，从前期（上一年 10 月至当年 7 月）126项环流指数中普查可能的预测因子，需满足2个条件：（1）前期环流指数与寒露风初日存在显着的高相关关系，通过0.05信度的显着性检验；（2）前期环流指数与寒露风初日的距平符号一致率达到70%及以上。将满足条件的环流指数作为可能的预测因子；然后，利用逐步回归法，并结合各个可能预测因子之间的独立性检验进行预测因子的筛选；最后，利用最小二乘法用筛选的预测因子拟合寒露风初日的回归预测模型。
考虑到126项环流指数从1980年开始才完整，模型验证时段选为 1981—2010 年，气候平均态使用1981—2010 年。采用交叉检验的办法去验证寒露风初日趋势预测模型的预测效果，具体为逐一去除1981—2010 年中的任一年用于建模，将去除年份作为预测年，由此可获得 1981—2010年任一年的预测结果。2011—2013 年作为独立检验时段，具体是利用1981 年至预测的前一年用来建模，然后预测当年的寒露风初日趋势。由图4可知，寒露风初日趋势实况和预测的相关系数为0.89，通过0.001信度的显着性检验，明显高与 9 月平均气温距平的相关关系(0.57)。
1981—2010 年预测模型交叉检验阶段只有 1988 年和1992 年寒露风初日趋势预测错误，其余 28 年均正确，预测正确率为93.3%，远高于1981—2010年9月平均气温与寒露风初日的距平符号一致率(64%)，预测正确率相对提高46%。2011—2013年独立检验时段均预测寒露风初日偏早，预测趋势均正确，也好于9月平均气温与寒露风初日的距平符号一致率(67%)，两者在2011 年和 2012 年距平符号一致，在 2013 年 9 月平均气温正距平而寒露风初日偏早（负距平）。表明寒露风初日趋势预测模型的预测准确率要明显优于直接利用9月平均气温距平符号的趋势预测方法。【图4】
3 结论
（1）1961—1990年气候平均态下，寒露风初日沿江东部出现在9月26—27日，中部9月27—28日，西部9 月 29 日—10 月 2 日，江南比沿江出现的要早 4~7 天，出现在9月22—25日。1981—2010年与1961—1990年寒露风初日气候平均态的空间分布型基本一致，但出现日期明显推迟2~3天。
（2）在气候变暖背景下，1961—2013年沿江江南大部地区寒露风初日呈现推迟趋势，而沿江西部有提前趋势。寒露风初日推迟趋势中沿江东部最为显着，其次是江南东部和南部，再次是江南北部。
（3）寒露风初日与9月平均气温不论是气候平均态和多年变化趋势的空间分布型，还是两者在1961—2013 年的时间演变特征均十分一致，为利用 9 月平均气温趋势预测寒露风初日趋势提供科学依据。
（4）在气候变暖背景下，1961—2013年区域寒露风初日呈现显着的推迟趋势，但其存在年代际差异，21世纪以来明显提前，仍会出现异常偏早的情况。
（5）构建基于前期物理因子的寒露风初日趋势预测模型，该预测模型不论在1981—2010年交叉检验时段，还是在2011—2013年独立检验时段，对于寒露风初日趋势的预测正确率明显优于目前业务上直接利用9 月平均气温距平符号的趋势预测方法。
4 讨论
本研究构建了安徽省寒露风初日趋势预测模型，该模型预测准确率明显优于以往的业务预测方法，具有优异的业务应用前景，有助于提高寒露风气候预测专项服务水平，对于预防和减轻寒露风引起的农业气象灾害具有重要意义。
在气候变暖背景下，寒露风呈现推迟减弱的趋势，生产中对防御寒露风灾害有所轻视。但值得注意的是，21 世纪以来寒露风初日明显提前，尤其 2010—2013 年连续 4 年偏早，这一点必须引起各有关部门和粮农的高度重视。
寒露风专项气候预测业务虽然在寒露风初日趋势方面已取得良好的服务效果，然而目前气候预测业务水平与服务需求仍存在巨大差距。在政府决策部门和双季晚稻生产部门十分关注的寒露风引起的低温程度、发生时段、气象灾害风险等级等预测服务需求方面，气候预测水平十分有限，预测业务仍然空白。有待从寒露风气象灾害角度加强对寒露风专项预测业务技术的研发。
参考文献 [1] 陈利东,黄永森,匡昭敏.寒露风对玉林 2011 年晚稻产量的影响调查分析[J].中国农学通报,2013,29(18):109-113. [2] 安徽省气象局资料室.安徽气候[M].安徽:安徽科学技术出版社,1983:28-30. [3] 刘凤舞,刘志澄,矫梅燕.安徽省防汛抗旱气象手册[M].北京:气象出版社,1998:24-26. [4] 于波,鲍文中,吴必文,等.安徽省农业气象业务服务手册[M].北京:气象出版社,2013:154-157. [5] 刘文英,张显真,简海燕.江西近 50 年寒露风演变趋势及其对双季晚稻的影响[J].气象与减灾研究,2009,32(4):67-70.