干热风是我国北方小麦扬花灌浆期间出现的一种高温、低湿并伴有一定风力的灾害性天气, 主要分为高温低湿型、雨后青枯型、旱风型3种类型^[1]。高温低湿型在小麦扬花灌浆过程中都可能发生, 多发生在开花后20 d左右至蜡熟期, 气温突升, 空气湿度骤降, 并伴有较大风速。雨后青枯型发生在小麦成熟前10 d内, 其特征是雨后猛晴, 气温骤升, 湿度剧降, 造成小麦热枯逼熟。旱风型主要发生在新疆和西北黄土高原的多风地区, 干旱年份出现较多, 风速大、湿度低, 并与一定高温配合, 使叶片失水卷缩甚至撕裂破碎。干热风发生时, 小麦叶片叶绿素含量和光合速率下降, 受精结实率降低, 灌浆速度下降、灌浆期缩短, 从而导致减产, 危害轻的年份减产10%以下, 严重年份减产10%—20%甚至更多^[2]。黄淮海地区是我国冬小麦主产区, 也是干热风危害最重、影响范围最广的区域, 该地区干热风主要有高温低湿和雨后青枯两种类型^[1]。近10a(2006—2015年), 该区域冬小麦种植面积占全国冬小麦面积的71%, 产量则占全国的80%^[3], 该区域小麦产量的多少直接影响全国粮食总产水平, 分析干热风灾害发生规律与时空演变特征, 对冬小麦安全生产、防灾减损和提质增效具有重要意义。

我国有关干热风的研究开始于20世纪50年代后期, 20世纪80年代初, 北方13省(市)小麦干热风科研协作组针对干热风伤害机理、气象指标、气候区划、预报方法、防御技术等进行了联合攻关, 为后来开展相关研究奠定了基础^{[1, 4-9]}。其后人们针对省级或地市级尺度的干热风发生演变规律、风险区划、预测防御等开展了较多研究^[10-26], 针对省级以上尺度的研究相对较少^[27-30], 多数研究表明黄淮海地区干热风发生频率总体呈减少趋势。干热风影响定量评估技术仍处于探索阶段, 刘静等^[31]建立了干热风灾损评估模型, 用于评估宁夏引黄灌区春小麦干热风影响; 朱玉洁等^[32]探讨了基于作物模型的干热风灾损提取方法; 赵俊芳等^[33]对黄淮海地区冬小麦干热风灾损分离提取方法进行了探讨; 李颖等^[34]对基于遥感的大面积干热风灾害监测评估进行了尝试, 这些研究对干热风灾害的科学防御提供了技术支撑。

作为一种农业气象灾害, 干热风与小麦生育期紧密关联, 其影响也并非单一类型干热风作用的结果, 已有研究通常只考虑单一类型干热风^{[11-23, 29-30]}, 或将高温低湿型和青枯型分开分析^[25-27], 难以全面评价干热风总体特征, 研究时段也多以自然旬月划定^[11-19], 虽然能反映其发生规律, 但与小麦实际发育期吻合度有所欠缺, 可能造成统计样本缺失, 或部分扬花灌浆期以外的干扰样本影响分析结果。本文基于公认的干热风灾害指标, 以及气象和冬小麦生育期观测资料, 对黄淮海冬小麦主产区高温低湿和雨后青枯两种类型干热风的总体发生规律进行研究, 以期更好地揭示黄淮海地区干热风灾害的时空演变特征, 为区域大范围干热风监测预警及防御提供参考。

1 资料与方法

研究区域为黄淮海冬小麦主产区, 包括河北中南部、北京、天津、山东, 以及河南、安徽和江苏3省淮河以北地区, 基于区域内65个农业气象观测站的资料进行干热风时空特征分析(图 1)。所用气象资料为1961—2015年共55 a的逐日最高气温、14:00相对湿度、14:00风速、降水量数据, 作物资料为1981—2015年的冬小麦发育期数据, 以上资料均来自国家气象信息中心。

通常采用日最高气温、14:00相对湿度和14:00风速组合确定干热风日, 用不同等级干热风日的天数组合确定干热风过程。高温低湿型干热风日指标为, 轻度：日最高气温≥32℃、14:00相对湿度≤30%、14:00风速≥2 m/s; 重度：日最高气温≥35℃、14:00相对湿度≤25%、14:00风速≥3 m/s^[1]。雨后青枯型干热风日指标为, 小麦成熟前10 d内, 有1次小到中雨, 雨后猛晴, 3 d内有1 d日最高气温≥30℃, 14:00相对湿度≤40%, 14:00风速≥3 m/s^[1]。根据小麦千粒重降低标准^[1], 将青枯型干热风日记为重干热风日, 进行干热风日数和过程的综合分析。干热风过程的判定方法见表 1^[1-2]。

干热风危害发生在小麦生育后期, 籽粒形成期至成熟前的30 d内, 主要受害时期是灌浆中后期^[1], 本研究分析时段为成熟前的30 d, 而不局限于自然旬月时间。1961—1980年无发育期资料, 分析时用1981—2015年发育期计算各站平均发育日期, 以发育期常年值确定干热风分析时段。

利用Excel、VB.net程序进行数据处理与计算, 其中干热风日数和干热风过程变化的气候倾向率采用线性回归方法确定^[35-36], 空间分布图利用ArcMap反距离权重插值模块制作。

2 结果与分析 2.1 干热风日时空特征 2.1.1 干热风日数空间分布

干热风年平均日数是反映干热风强度和危害轻重的基本特征值, 与小麦灌浆速率下降关系密切, 通常干热风日数愈多, 危害愈重^[1]。从图 2可以看出, 河北大部、山东中北部、河南北部年平均干热风日数在3 d以上, 其中河北中南部的石家庄、保定、沧州西部、衡水、邢台等地区年平均干热风日数达5.0—6.5 d, 是黄淮海冬麦区出现干热风最多的地区。从轻、重干热风日数分布来看, 京津冀、山东北部、河南东北部平均轻干热风日数为2.0—3.6 d, 河北中南部、鲁西北重干热风日数则有2.0—3.1 d。

图 3、图 4分别为黄淮海地区干热风年平均日数频率与保证率空间分布。黄淮海地区大部年平均出现1—4 d干热风的频率最高, 其中河南大部、苏皖北部、山东中西部大部和半岛北部频率超过50%;河北中南部出现1—4 d的频率为40%—50%, 有30%—40%的年份出现5—8 d干热风, 另有10%—20%的年份可能出现9 d以上的干热风天气, 为干热风天气多发区(图 3)。京津冀、山东中部和北部、河南中部和北部年干热风日数80%保证率的日数为1—4 d, 在河北中南部的保定、石家庄东部、沧州西部、衡水西部、邢台等地有50%保证率出现5—8 d干热风(图 4)。

2.1.2 干热风时间分布

了解干热风发生初日与集中出现时期, 对于提前做好干热风防御具有重要意义。根据各站点1961—2015逐年干热风发生初日, 统计计算干热风平均初日、最早初日并绘制空间分布图。受到冬小麦进入灌浆期的早晚、内陆沿海差异、地形地势等因素的影响, 干热风平均初日总体从南向北、从西向东逐渐推迟, 豫西、豫南出现最早(5月16日—20日), 山东半岛东部和冀东北最晚(5月31日—6月4日), 其他地区一般在5月21日—30日(图 5)。干热风最早初日空间分布特征与平均初日较为相似, 豫西、豫南最早(4月下旬), 山东半岛东部和冀东北最晚(5月16日以后), 其他地区一般在5月1日—15日(图 6)。

统计1961—2015年冬小麦成熟前30 d内逐候(生育进程5 d)干热风出现的日数及比例(表 2), 可以看出随着时间推移, 干热风轻日数、重日数和总日数比例均呈逐候增多趋势, 第4—6候干热风轻日、重日和总日数发生比例分别为73.5%、93.1%和81.1%, 尤其重干热风日, 有86.9%出现在成熟前10 d内(第5—6候)。干热风日数随时间变化的特征, 与春末夏初北方雨季来临前多晴少雨、气温迅速回升的天气背景密切相关, 也说明灌浆中后期是干热风危害集中期和防御关键期, 成熟前10 d内尤其注意防范雨后青枯型干热风引起的早衰逼熟。

2.1.3 干热风日数时空变化

将一年内所有站点的轻、重干热风日数及干热风总日数分别累计并求其平均值, 绘制时间序列图(图 7)。可以看出, 1961—2015年期间, 黄淮海地区干热风日数总体呈减少趋势, 其中重干热风日数下降趋势更明显, 轻干热风则无明显变化趋势。从年代际对比来看, 20世纪60年代平均干热风日数最多, 影响最重, 90年代平均日数最少, 影响最轻, 其中有4 a区域平均干热风日数超过6 d, 分别是1962年(6.2 d)、1965年(6.8 d)、1967年(6.1 d)、2001年(8.1 d), 2001年为1961—2015年干热风日数最多的年份。

取各站点1961—2015年逐年干热风日数, 分站点计算其气候倾向率并绘制成图, 用以表征各地干热风日数的增减趋势(图 8)。黄淮海冬麦区大部干热风日数呈减少趋势, 其中山东西部和北部、河北东南部、河南东北部等地减少趋势更为明显; 豫西南、京津和冀东北等地呈增加趋势, 其中天津增加趋势较为明显。干热风日数的时空变化特征说明该地区干热风威胁总体在减轻, 其中在多发区、危害集中区减轻趋势更为明显。

2.2 干热风过程时空特征 2.2.1 过程次数空间分布

干热风年平均过程次数也是干热风强度和危害程度的重要特征值, 是干热风日数、过程极端最高气温、最低相对湿度, 以及持续危害高温的综合体现。通常一个地区年干热风过程次数越多, 过程强度越强, 对小麦造成的危害也越重。干热风年平均过程次数的空间分布趋势(图 9), 与日数分布总体是一致的, 京津冀大部、山东北部、豫北部分地区年平均干热风过程次数超过1.5次, 其中河北中部年均过程次数为2.0—2.4次, 为干热风过程次数最多的地区, 豫西南和豫东南、苏皖北部、鲁东南和半岛地区年平均过程次数少于1次。

2.2.2 干热风过程时空变化

从黄淮海地区年平均干热风过程次数时间序列图可以看出(图 10), 1961—2015年期间, 黄淮海地区干热风过程次数总体呈减少趋势, 其中重过程次数下降趋势明显, 轻过程次数的下降趋势不显著。干热风年平均过程次数与年平均日数变化趋势的一致性, 表明干热风过程次数与日数之间有极为密切的相关性和同步性, 日数越多, 过程次数越多。干热风过程次数变化趋势的空间分布与日数变化趋势的空间分布也较为一致, 大部地区呈减少趋势, 其中鲁西北、冀东南、豫东北等地减少趋势更为明显, 京津、豫西南等地则呈弱增多趋势(图 11)。

2.2.3 干热风过程持续日数

干热风过程强弱还反映在过程持续日数上, 持续日数越多, 过程越强。一般来说, 持续2 d的干热风可使局部麦田受害, 持续3—4 d可使1/3—1/2的地块受害, 持续5—9 d小麦基本都会受害, 处于乳熟-蜡熟的小麦就该收获了, 若持续10 d以上不论处在哪个灌浆阶段都能普遍干枯提前收获^[7]。统计1961—2015年不同持续日数的过程占干热风过程总数的比例, 结果表明干热风过程持续日数集中在1—4 d, 共占过程总数的95.3%, 持续5 d以上的过程所占比例不足5%。其中, 单个重干热风日构成的轻过程所占比例最高, 为37.5%;其次为持续2 d的过程, 占36.9%;持续3d、4d的过程分别占15%、5.9%。干热风过程持续日数最高值出现在1992年的河北饶阳, 共持续9 d。

3 结论与讨论

1961—2015年, 就空间平均分布而言, 河北中南部、河南北部、山东北部和西部等地年平均干热风日数、干热风过程均为黄淮海冬麦区高值区, 其中地处太行山东部的河北中南部由于强烈的下沉增温和焚风效应, 以及干旱、盐碱土壤等综合影响, 形成黄淮海地区的干热风重发区。

从时间分布来看, 受小麦生育进程、内陆沿海差异、地形地势等因素影响, 干热风平均初日和最早初日总体从南向北、从西向东逐渐推迟。在春末夏初北方雨季到来前多晴少雨、气温迅速回升的天气背景下, 随着小麦灌浆进程, 干热风发生日数呈逐渐增多趋势, 灌浆中后期干热风日数多、程度重, 是干热风危害的集中期和防御关键期。

从多年变化趋势来看, 黄淮海地区干热风日数、过程次数总体均呈减少趋势, 其中重干热风日数、重过程次数减少趋势更为明显; 从地区差异来看, 冀东南、鲁西北、豫东北干热风日数和过程次数减少趋势更为明显, 表明黄淮海冬麦区干热风威胁总体在减轻。温、湿、风要素满足指标的空间格局是干热风随时间和空间分布的基础, 在黄淮海地区冬小麦生长后期, 温、湿、风要素达到干热风标准的日数均呈减少趋势^[29], 这对当地小麦高产、稳产有利。而在干热风3个致灾因子中, “热”是主导因子, “干”居第二位, “风”则起辅助作用^{[29, 37]}, 其中温度条件相对容易达标^[27], 而在气候变暖背景下, 极端天气气候事件频发, 部分年份仍存在发生较重干热风的可能。如2001年区域平均干热风日数达8.1 d, 成为1961—2015年干热风日数最多的年份, 因此对干热风的防御仍不能掉以轻心。

本文研究区域针对黄淮海冬小麦主产区, 所指的干热风不是气象上泛指的干热风天气, 而是指小麦灌浆成熟期间发生的干热风。由于各地冬小麦发育期存在差异, 本研究对干热风分析时段的选取与冬小麦发育期相吻合, 可以较好地反映不同地区干热风发生的客观规律, 避免了整个区域按相同时段分析带来的偏差。

干热风对小麦灌浆和产量的影响, 与干热风出现的时段、强度、频率, 以及田间小气候和小麦抗干热风能力密切相关。现有干热风指标主要是基于20世纪80年代初的研究成果, 随着小麦品种更替与改良、耕作栽培技术进步, 特别是灌溉条件的改善和“一喷三防”技术的推广, 黄淮海地区冬小麦抗干热风能力大大提高, 有时虽出现较大范围干热风天气, 但对产量并未造成大的影响。此外, 干热风危害是高温低湿瞬时强度和持续时间综合作用的结果, 随着气象自动观测技术的普及, 现有观测手段已能满足逐时甚至分钟观测的需求, 为研究高温低湿的累积效应提供了基础。因此, 有必要结合当今农业生产条件、气象自动观测技术和气候资源特点对干热风指标做适当修订, 特别是考虑土壤墒情、田间小气候以及小麦生理参数等的综合影响, 以使监测评价结果与当今农业生产实际更相适应。