中国生态农业学报(中英文)  2020, Vol. 28 Issue (3): 337-348  DOI: 10.13930/j.cnki.cjea.190723
0

引用本文 

陈怀亮, 李树岩. 气候变暖背景下河南省夏玉米花期高温灾害风险预估[J]. 中国生态农业学报(中英文), 2020, 28(3): 337-348. DOI: 10.13930/j.cnki.cjea.190723
CHEN H L, LI S Y. Prediction of high temperature disaster risks during summer maize flowering under future climate warming background in Henan Province[J]. Chinese Journal of Eco-Agriculture, 2020, 28(3): 337-348. DOI: 10.13930/j.cnki.cjea.190723

基金项目

国家重点研发计划课题(2018YFD0300704,2017YFD0300304)和中国气象局气候变化专项(CCSF201408)资助

通信作者

李树岩, 主要从事农业气象灾害评估及作物模型应用研究。E-mail:lsy_126com@126.com

作者简介

及其

文章历史

收稿日期:2019-10-12
接受日期:2019-12-13
气候变暖背景下河南省夏玉米花期高温灾害风险预估*
陈怀亮1,2, 李树岩1,3     
1. 中国气象局/河南省农业气象保障与应用技术重点开放实验室 郑州 450003;
2. 河南省气象局 郑州 450003;
3. 河南省气象科学研究所 郑州 450003
摘要:为预估未来气候变暖背景下夏玉米花期高温灾害风险,根据河南省19个农业气象观测站夏玉米抽雄期常年观测资料和未来RCPs(representative concentration pathways)气候变化情景数据,构建夏玉米花期高温风险评价指标,开展河南省夏玉米花期高温灾害时空特征及风险演变分析。其中RCPs气候情景数据包括基准气候条件(1951-2005年,RCP-rf)和未来(2006-2050年)RCP 4.5(中)、RCP 8.5(高)两种浓度路径数据。以抽雄普遍期及之后7 d确定为夏玉米花期,并内插匹配气候情景格点数据。以花期最高气温≥ 32℃和≥ 35℃作为轻度和重度高温灾害发生阈值,根据轻、重度夏玉米花期高温发生频率和高温积害,建立风险评价指标并分级。结果表明,RCP-rf情景下全省夏玉米花期高温发生频率在20.5%~81.0%(≥ 32℃)和3.9%~51.9%(≥ 35℃)。与基准条件相比,≥ 32℃高温发生频率增加9.1%(RCP 4.5)和11.0%(RCP 8.5),≥ 35℃高温发生频率增加8.7%(RCP 4.5)和8.3%(RCP 8.5)。RCP-rf情景下全省夏玉米花期高温积害在48.5~200.9℃·d(≥ 32℃)和9.8~138.5℃·d(≥ 35℃)。与基准条件相比,≥ 32℃高温积害增加25.4℃·d(RCP 4.5)和25.6℃·d(RCP 8.5),≥ 35℃高温积害增加25.8℃·d(RCP 4.5)和31.4℃·d(RCP 8.5)。由综合风险分析可知,RCP-rf情景下夏玉米花期高温灾害高值风险区主要分布在新乡、郑州、许昌、漯河、周口及其以东以北的地区(商丘除外),约占夏玉米主栽区面积的30.1%;RCP 4.5情景下高值风险区扩大至洛阳和南阳以东的大部分地区,约占夏玉米主栽区面积的63.4%;RCP 8.5情景下高值风险区面积进一步向西扩大,约占夏玉米主栽区面积的占76.3%。
关键词RCP气候情景    气候变化    夏玉米    花期    高温风险    河南省    
Prediction of high temperature disaster risks during summer maize flowering under future climate warming background in Henan Province*
CHEN Huailiang1,2, LI Shuyan1,3     
1. China Meteorological Administration/Henan Key Laboratory of Agrometeorological Support and Applied Technique, Zhengzhou 450003, China;
2. Henan Meteorological Service, Zhengzhou 450003, China;
3. Henan Institute of Meteorological Sciences, Zhengzhou 450003, China
Abstract: To predict the risk of high temperature disasters during the summer maize flowering period under future climate change scenarios, a high temperature risk index was established based on the long-term observational data of the summer maize booting period from 19 agro-meteorological observational stations in Henan Province and future climate change scenarios (representative concentration pathways, RCPs). Based on the high temperature risk index, the spatio-temporal characteristics and risk of high temperature disasters during the summer maize flowering period were analyzed in Henan Province. The climate change scenarios (RCPs) included the baseline climate scenario (RCP-rf, 1951-2005) and two future climate change scenarios with moderate emission (RCP 4.5) and high emission (RCP 8.5) during 2006-2050. The summer maize flowering period at each site was defined as the period from the booting stage to 7 days after booting, and then, it was interpolated into the whole province to match the gridded date of climate scenarios. The maximum temperature of higher than 32℃ and 35℃ during flowering were used as the threshold for occurrence of light and heavy high temperature disasters. A comprehensive risk index was successfully established, and the risk level was divided by integrating high temperature frequency and accumulation of temperature higher than 32℃ and 35℃ during the summer maize flowering period. The results showed that under the RCP-rf scenario, the high temperature frequency during the summer maize flowering period was 20.5%-81.0% (≥ 32℃) and 3.9%-51.9% (≥ 35℃). Compared with the baseline, the frequency of high temperature days of more than 32℃ during the summer maize flowering period increased by 9.1% (RCP 4.5) and 11.0% (RCP 8.5), and that of higher than 35℃ increased by 8.7% (RCP 4.5) and 8.3% (RCP 8.5). Under the RCP-rf scenario, the accumulation of temperature higher than 32℃ during the summer maize flowering period ranged 48.5-200.9℃·d, and that of higher than 35℃ ranged 9.8-138.5℃·d. Compared with the baseline, the accumulation of temperature higher than 32℃ during the summer maize flowering period increased by 25.4℃·d (RCP 4.5) and 25.6℃·d (RCP 8.5); and that of higher than 35℃ increased by 25.8℃·d (RCP 4.5) and 31.4℃·d (RCP 8.5). According to the comprehensive risk analysis of high temperature during the summer maize flowering period, the high-risk zones under RCP-rf scenario were distributed in the eastern and northern areas, which adjacent to Xinxiang, Zhengzhou, Xuchang, Luohe, and Zhoukou, except Shangqiu, accounting for approximately 30.1% of the main summer maize planting area. The high-risk areas covered most of the areas in the east of Luoyang and Nanyang, accounting for 63.4% of the main area of summer maize planting under the RCP 4.5 scenario and 76.3% of the main area of summer maize planting under the RCP 8.5 scenario.
Keywords: RCP scenarios    Climate change    Summer maize    Flowing period    High temperature risk    Henan Province    

玉米(Zea mays)是我国第一大粮食作物, 在国民经济发展中具有举足轻重的地位, 河南是夏玉米种植大省, 播种面积331.7万hm2, 总产1 752.1万t[1]。随着玉米产量水平的提高, 穗粒数对产量的决定作用愈加显著[2]。花粉活力是影响穗粒数的重要原因之一, 虽然玉米雄穗的耐热和散粉特性存在着基因型间差异, 但高温是造成玉米花粉活力减退的重要原因[3-4]。异常高温影响花粉结构和功能, 造成花粉的数量和活力下降或花期不相遇, 最终导致籽粒败育结实率下降和减产[5-7]。赵龙飞等[8]研究表明在开花期对玉米品种‘浚单20’进行高温处理会使穗粒数下降20.4%~22.0%, 百粒重下降8.8%~10.5%。玉米抽雄到吐丝期适宜温度为25~28 ℃, 但河南夏玉米抽雄散粉、吐丝期一般处于7月下旬到8月初的高温期, 特别是在全球气候变暖的大背景下, 夏玉米花期常遭遇高温干旱等极端天气, 作物产量急剧下降[9-10], 且目前生产上还缺乏有效的防御或缓解措施, 严重威胁玉米的安全生产[11]。花期高温已成为影响河南省夏玉米生长的主要气象灾害之一。

在高温灾害指标的研究方面, 气象上通常将35 ℃作为高温绝对阈值, 但农业生产中即使出现相对高温也可使作物发生生理障碍而减产[12]。针对夏玉米花期高温的危害特征和影响机理前人做了大量深入研究[13-20], 形成了一些高温判定指标。如陈朝辉等[17]研究发现, 玉米处于38 ℃的极端高温下3 d, 就会停止散粉, 授粉率与温度呈负相关。徐翠莲等[18]研究揭示, 温度和湿度对花粉数量与寿命影响最大, 散粉期温度≥28 ℃花粉量会明显减少; 若温度高于35 ℃甚至38℃时, 花粉活力很快丧失。王海梅[19]研究表明大于32 ℃高温胁迫将对河套灌区玉米生理指标和产量构成产生影响。徐美玲[20]指出玉米雌花活力与温度相关, 气温越高寿命越短, 气温达34.3~37.8 ℃时, 花丝寿命只有72 h。但这些生理机制方面的高温热害指标, 较难应用在生产实践中。李德等[21]根据极端最高气温、日最高气温≥35 ℃的日数和日最高气温高于35 ℃的积害量及其期间的平均最小相对湿度构建了淮北平原夏玉米高温热害综合气候指标。王丽君[22]以日最高温大于32 ℃的高温累积度日和日最高温大于32 ℃的天数作为评价玉米生育期内极端高温风险的两个指标。这些指标从气象监测数据和灾情调查资料出发, 更便于应用。

在高温风险评估和区划方面针对水稻(Oryza sativa)开展的研究较多, 对玉米高温风险研究相对较少, 主要包括刘哲等[23-24]利用日高温时长概率分布函数计算了黄淮海各县区玉米花期高温热害风险概率和空间分布, 并利用MODIS数据分析了2011—2014年黄淮海夏玉米高温风险。尹小刚等[25]根据≥30 ℃的积温和日数分析东北地区高温对玉米生产的影响。在已有研究中结合气候情景数据, 探讨未来气候变化条件下花期高温灾害风险演变特征的研究甚少。当前气候变暖已是不争的事实, 未来气候增温还会持续[26-27], 夏玉米抵御高温热害, 保障安全生产的形势将更加严峻。

本研究重点在于应用未来气候变化数据分析夏玉米花期高温风险演变。首先, 借鉴前人研究成果提取夏玉米花期≥32 ℃或≥35 ℃高温日数、发生频率和高温积害, 构建夏玉米花期高温风险评估指标。然后引入RCP气候情景数据, 分析未来不同排放情景下夏玉米花期高温灾害的时空分布特征, 预估未来夏玉米花期高温风险。研究结果对于调整玉米生产系统与抗逆栽培, 制定气候变化适应对策具有重要意义。

1 材料与方法 1.1 研究区域概况

河南省位于黄河中下游(110°~117°E, 31.5°~36.5°N), 北亚热带向暖温带过渡地区, 地势西高东低, 主要作物系统为冬小麦-夏玉米轮作。夏玉米6月上旬播种, 9月下旬收获, 夏玉米生长季太阳总辐射1 900~2 400 MJ·m-2, 降水量400~600 mm, ≥10 ℃积温2 600~3 100 ℃·d, 适宜玉米生长。研究区域为河南省夏玉米主要种植区, 其中信阳地区主要种植水稻, 不作为研究区。

1.2 数据来源 1.2.1 RCP气候情景数据

研究未来气候变化对农业的影响需要结合气候情景数据, IPCC第5次评估报告(AR5)中采用了融入政策因素的代表性浓度路径(representative concentration pathways, RCPs)情景预估未来气候变化趋势[28]。包括历史气候模拟(historical)RCP-rf数据和RCP 4.5、RCP 8.5两种排放情景下未来气候变化预估数据, 气候模式模拟的1961—2005年气温与观测值的偏差除青藏高原外大部分地区在±1 ℃之间[28]。模式水平分辨率为50 km。时间尺度上基准气候条件(RCP-rf)为1951—2005年, RCP 4.5和RCP 8.5为2006—2050年。RCP气候情景数据覆盖河南全省区域共165个格点, 其空间分布如图 1所示。

图 1 河南省RCP情景数据示意图 Fig. 1 Sketch map of climate data in Henan Province under RCP scenarios
1.2.2 夏玉米发育期数据

来源于河南省19个夏玉米农业气象观测站, 其空间分布如图 2所示。观测年代选取记录完整的1988—2017年共30年。

图 2 河南省夏玉米主产区和农业气象站分布示意图 Fig. 2 Map of summer maize production area and distribution of agro-meteorological stations in Henan Province 阴影部分为信阳地区, 以种植水稻为主, 不作为研究区域。The shaded area is Xinyang City where rice is the main crop, so is not included in the study area.
1.3 夏玉米花期确定

以抽雄普遍期作为花期的开始日期, 抽雄普遍期向后推迟7 d为花期的结束日期。根据1988—2017年的农业气象观测资料, 计算夏玉米花期开始和结束日期的多年平均值, 即花期常年值。根据距离相近原则, 将花期起止日期常年值推算到河南省区域内的RCP气候格点数据上。夏玉米抽雄期日序(花期开始日期)的常年值如图 3所示。

图 3 河南省夏玉米抽雄期日序常年值 Fig. 3 Average Julian day series of booting of summer maize in Henan Province 信阳地区以种植水稻为主, 不作为研究区域。Xinyang City is not included in the study area because where rice is the main crop.
1.4 夏玉米花期高温热害指标

降志兵等[5]研究表明温度达到32 ℃且持续60 min, 玉米小花受精率和总结实率显著低于对照; 当温度高于35 ℃时小花受精率和总结实率也显著低于32 ℃处理。因此, 将≥32 ℃和≥35 ℃作为两个不同程度高温热害的致灾阈值, 即≥32 ℃为轻度受害, ≥35 ℃为重度受害。综合考虑夏玉米花期高温热害发生的频率和强度, 结合前人研究成果[21-22, 29], 确定≥32 ℃和≥35 ℃的“高温日数”和“高温积害”作为评价花期高温影响的指标。

1) 高温日数:在1.2中确定的夏玉米花期7 d内, 当日最高气温≥32 ℃时作为轻度高温日, 日最高气温≥35 ℃时作为一个重度高温日, 每年花期轻、重高温灾害发生的总日数为轻、重度高温日数。

2) 高温积害:以花期最高气温≥32 ℃或≥35 ℃的高温累积值, 来表征高温灾害发生的严重程度, 单位为℃·d。计算公式如下:

${\rm{T}}{{\rm{H}}_i} = \sum\nolimits_{i = 1}^n {{T_{hi}}} $ (1)
${T_{hi}} = \left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} 0\\ {{T_{{\rm{max}}}}} \end{array}}&{\begin{array}{*{20}{c}} {{\rm{(}}{T_{{\rm{max}}}} < {T_0}{\rm{)}}}\\ {{\rm{(}}{T_{{\rm{max}}}} \ge {T_0}{\rm{)}}} \end{array}} \end{array}} \right. $ (2)

式中: THi表示≥32 ℃或≥35 ℃的高温积害值, 即THi可表示为TH32或TH35; Thi为当日积害值; Tmax为逐日最高气温; T0为夏玉米花期最高温的临界值, 即32 ℃或35 ℃。

1.5 夏玉米花期高温灾害发生频率

统计多年夏玉米花期高温日数除以花期总日数, 可计算高温灾害发生的频率, 公式如下:

${P_i} = {n_i}/{N_i}\times100\% $ (3)

式中: Pi为花期≥32 ℃或≥35 ℃高温的发生频率, 即Pi可表示为P32P35; ni为夏玉米花期内≥32 ℃或≥35 ℃高温日数; Ni为夏玉米花期总日数。

1.6 夏玉米花期高温风险综合指数及风险等级划分

气候风险指数可用概率乘以强度来表示, 夏玉米花期高温风险综合指数计算方法如下:

$I = {\omega _1} \times {P_{32}} \times {\rm{T}}{{\rm{H}}_{32}} + {\omega _2} \times {P_{35}} \times {\rm{T}}{{\rm{H}}_{35}} $ (4)

式中: I为夏玉米花期高温风险综合指数, ω1ω2分别为不同程度高温影响的权重系数, P32P35分别为≥32 ℃或≥35 ℃的高温发生频率, TH32和TH35分别为对应等级的高温积害值。

ω1ω2确定方法:降志兵等[5]研究指出高温处理1 h后, ≥32 ℃或≥35 ℃处理受精结实率均与对照呈显著差异, 较对照分别降低49%和65%, 因此权重系数取值为二者损失率的相对比例, 计算方法如下:

${\omega _{\rm{1}}} = \frac{{49{\rm{\% }}}}{{49{\rm{\% }} + 65{\rm{\% }}}} = 0.43 $ (5)
${\omega _2} = \frac{{6{\rm{5\% }}}}{{49{\rm{\% }} + 65{\rm{\% }}}} = 0.57 $ (6)

以RCP-rf基准条件为参考, 将夏玉米花期高温发生风险指数划分为轻、中、重3级。分级的阈值计算公式如下:

${I_a}{\rm{ = }}{a_i} \times {I_{{\rm{max}}}} $ (7)

式中:Ia为分级的阈值; Imax为全区域最高高温风险综合指数; ai为分级系数, 其中轻度与中度风险分级系数取0.4, 中度与重度风险分级系数取0.7。根据阈值大小, 将计算得到的分级阈值直接应用到RCP 4.5和RCP 8.5情景下进行轻、中、重风险等级划分。

1.7 数据处理分析方法

采用Microsoft Excel进行数据处理, 利用Suffer软件选择克里金插值方法作图。

2 结果与分析 2.1 夏玉米花期高温日数及灾害发生频率 2.1.1 年际变化

夏玉米花期≥32 ℃高温日数变化如图 4a1-a3所示。RCP-rf情景下≥32 ℃高温日数为0.4~6.8 d, 多年平均4.2 d, 呈显著上升趋势(P < 0.05); RCP 4.5情景下多年平均为4.7 d, 变化趋势不显著; RCP 8.5情景下增温幅度更大, ≥32 ℃高温日数平均为4.8 d, 呈显著上升趋势(P < 0.05)。夏玉米花期≥35 ℃高温日数变化如图 4b1-b3所示, RCP-rf情景下≥35 ℃高温日数多年平均为2.0 d, 也呈显著上升趋势(P < 0.05); RCP 4.5情景下多年平均为2.7 d, 变化趋势均不显著; RCP 8.5情景下多年平均为2.8 d, 呈显著上升趋势(P < 0.05)。

图 4 RCP情景下河南省夏玉米花期≥32 ℃和≥35 ℃的高温日数 Fig. 4 Changes of days of temperature higher than 32 ℃ and 35 ℃ during summer maize flowering in Henan Province under RCP scenarios

夏玉米花期≥32 ℃和≥35 ℃高温发生频率与高温日数的变化趋势一致(图略), RCP-rf情景下夏玉米花期高温≥32 ℃和≥35 ℃发生频率多年平均为61.3%和28.8%; RCP 4.5情景下夏玉米花期高温发生频率上升趋势不显著, 多年平均为69.6%(≥32 ℃)和38.3%(≥35 ℃); RCP 8.5情景下夏玉米花期高温发生频率呈显著上升趋势(P < 0.05), 多年平均为70.4%(≥32 ℃)和40.1%(≥35 ℃)。≥35 ℃高温发生频率的波动性增加, 未来情景下变异系数分别为15.4%(RCP 4.5)和13.9%(RCP 8.5), 均高于基准条件的12.1%。

2.1.2 频率空间分布

夏玉米花期≥32 ℃高温日数及灾害发生频率空间分布如图 5a1-a3所示。RCP-rf情景下夏玉米花期≥32 ℃高温日数全省为1.7~5.7 d, 高温发生频率全省为20.5%~81.0%;除三门峡、洛阳和南阳的西部外, 全省其他大部分地区≥32 ℃高温发生频率均在70%以上。RCP 4.5情景下夏玉米花期≥32 ℃高温日数为2.4~6.3 d, 高温发生频率全省为30.6%~89.9%; ≥32 ℃高温发生频率的高值区主要分布在新乡、郑州、平顶山及南阳一线以东的大部分地区, 频率大于80%。RCP 8.5情景下夏玉米花期≥32 ℃高温发生频率全省为36.1%~87.3%, 高温日数2.8~6.1 d, 发生频率大于80%的高值区分布与RCP 4.5情景相似。与基准条件相比, 未来排放情景下夏玉米花期≥32 ℃高温发生日数分别增加0.6 d (RCP 4.5)和0.5 d (RCP 8.5), 发生频率增加9.1% (RCP 4.5)和11.0% (RCP 8.5)。

图 5 RCP情景河南省夏玉米花期≥32 ℃和≥35 ℃的高温日数及发生频率空间分布 Fig. 5 Spatial distribution of days and occurrence frequency of temperature higher than 32 ℃ and 35 ℃ during summer maize flowering stage in Henan Province under RCP scenarios 信阳地区以种植水稻为主, 不作为研究区域。Xinyang City is not included in the study area because where rice is the main crop.

夏玉米花期≥35 ℃高温日数及灾害发生频率空间分布如图 5b1-b3所示。RCP-rf情景下夏玉米花期≥35 ℃高温日数为0.3~3.6 d, 高温发生频率全省为3.9%~51.9%;其中焦作、洛阳、平顶山和驻马店一线以东的大部分地区≥35 ℃高温发生频率均在40%以上。RCP 4.5情景下夏玉米花期≥35 ℃高温发生频率全省为8.8%~59.7%, 高温日数为0.8~4.2 d, ≥35 ℃高温发生频率的高值区主要分布在新乡、郑州、平顶山及驻马店一线以东的大部分地区, 高温发生频率在50%以上。RCP 8.5情景下夏玉米花期≥35 ℃高温发生频率全省为12.7%~56.3%, 高温日数为0.9~3.9 d, 发生频率大于50%高值区较RCP 4.5情景分布范围更广。与基准条件相比, 未来排放情景下夏玉米花期≥35 ℃高温发生日数分别增加0.6 d (RCP 4.5)和0.7 d (RCP 8.5), 发生频率增加8.7% (RCP 4.5)和8.3% (RCP 8.5)。

2.2 夏玉米花期高温积害 2.2.1 年际变化

夏玉米花期≥32 ℃高温积害年际变化如图 6a1-a3所示。RCP-rf情景下≥32 ℃高温积害多年平均为151.3 ℃·d, 呈显著上升趋势(P < 0.05); RCP 4.5情景下变化趋势不显著, 多年平均为174.9 ℃·d; RCP 8.5情景下夏玉米花期≥32 ℃高温积害也呈显著上升趋势(P < 0.05), 多年平均为177.9 ℃·d。夏玉米花期≥35 ℃高温积害年际变化如图 6b1-b3所示, RCP-rf情景下≥35 ℃高温积害多年平均为75.5 ℃·d, 呈显著上升趋势(P < 0.05); RCP 4.5情景下变化趋势不显著, 多年平均为99.5 ℃·d; RCP 8.5情景下夏玉米花期≥35 ℃高温积害也呈显著上升趋势(P < 0.05), 多年平均为106.1 ℃·d。

图 6 RCP情景河南省夏玉米花期≥32 ℃和≥35 ℃的高温积害变化 Fig. 6 Change of accumulation of temperature higher than 32 ℃ and 35 ℃ during summer maize flowering in Henan Province under RCP scenarios

不同情景下≥32 ℃高温积害变异系数差别较小, 但≥35 ℃高温积害变异系数在未来情景下为48.9%(RCP 4.5)和46.1%(RCP 8.5), 均高于基准条件的38.6%, 表明重度高温灾害发生的波动性更大。

2.2.2 空间分布

夏玉米花期≥32 ℃高温积害空间分布如图 7a1-a3所示。RCP-rf情景下全省≥32 ℃高温积害为48.5~200.9 ℃·d, 其中焦作、洛阳、南阳一线以东的大部分区域为高温积害高值区, 在180 ℃·d以上, 约占全省夏玉米主栽区面积的62%。RCP 4.5情景下全省≥32 ℃高温积害为73.4~231.3 ℃·d, 大于180 ℃·d的高值区主要分布在焦作、济源、洛阳和南阳一线以东的大部分地区, 分布面积约占全省夏玉米主栽区面积的82%。RCP 8.5情景下全省≥32 ℃高温积害为7.3~223.8 ℃·d, 分布形式与RCP 4.5情景相似, 积害在180 ℃·d以上的范围约占全省夏玉米主栽区面积的84%。与基准条件相比, 未来排放情景下夏玉米花期≥32 ℃高温积害分别增加25.4 ℃·d (RCP 4.5)和25.6 ℃·d (RCP 8.5)。

图 7 RCP情景河南省夏玉米花期≥32 ℃和≥35 ℃的高温积害空间分布 Fig. 7 Spatial distribution of accumulation of temperature higher than 32 ℃ and 35 ℃ during summer maize flowering in Henan Province under RCP scenarios 信阳地区以种植水稻为主, 不作为研究区域。Xinyang City is not included in the study area because where rice is the main crop.

夏玉米花期≥35 ℃高温积害空间分布如图 7b1-b3所示。RCP-rf情景下全省≥35 ℃高温积害为9.8~138.5 ℃·d, 安阳、濮阳、鹤壁、新乡、郑州、开封、许昌和周口的大部地区高温积害120℃·d以上, 约占全省夏玉米主栽区面积的25%。RCP 4.5情景下全省≥35 ℃高温积害为22.5~160.3 ℃·d, 大于120 ℃·d的积害高值区主要分布在焦作、济源、洛阳及南阳一线以东的大部分地区, 面积较基准条件明显增加, 约占夏玉米主栽区面积的60%。RCP 8.5情景下全省≥35 ℃高温积害为32.7~154.9 ℃·d, 大于120℃·d积害分布范围更广, 约占夏玉米主栽区面积的68%。与基准条件相比, 未来排放情景下夏玉米花期≥35 ℃高温积害分别增加25.8 ℃·d (RCP 4.5)和31.4 ℃·d (RCP 8.5)。

2.3 夏玉米花期高温综合风险分析

综合夏玉米花期高温发生频率和积害强度计算的综合风险指数空间分布如图 8所示。河南省夏玉米高温热害综合风险指数总体上呈东高西低的分布形式。RCP-rf情景下高值风险区主要分布在新乡、郑州、许昌、漯河、周口及其以东以北的地区(商丘除外), 约占夏玉米主栽区面积的30.1%;低值风险区主要分布在豫西三门峡、洛阳西部和南阳的西部。RCP 4.5情景下, 济源、洛阳和南阳一线以东地区均为高风险区; 低值风险区主要分布在豫西三门峡、

图 8 RCP情景河南省夏玉米花期高温综合风险分布 Fig. 8 Distribution of high temperature comprehensive risk during summer maize flowering in Henan Province under RCP scenarios 信阳地区以种植水稻为主, 不作为研究区域。Xinyang City is not included in the study area because where rice is the main crop.

洛阳西部和南阳的西北部, 面积较基准条件下明显减小。RCP 8.5情景下, 增温幅度更大, 高值风险区分布范围较RCP 4.5情景相比向西进一步扩大, 低值风险区面积进一步缩小。未来情景下高值风险区分布面积占夏玉米主栽区面积的63.4%(RCP 4.5)和76.3%(RCP 8.5), 较基准条件分别增加33.3% (RCP 4.5)和46.2% (RCP 8.5), 未来RCP情景下夏玉米花期高温灾害风险明显增加。

3 讨论

许多气候变化影响研究指出, 温度升高对作物生产的影响十分显著[30-34]。但大多关注的是平均温度对作物生长发育和产量的影响。然而作物产量与温度之间是非线性关系, 小于阈值产量随温度升高而增加, 高于阈值产量受灾下降, 而极端高温的发生将导致更大的产量损失[35-36]。降志兵等[5]研究表明极端高温38 ℃以上处理1 h, 玉米开花授粉将造成致命的损失, 受精结实率仅2%。未来气候变化情景下气候不稳定性增加, 经统计基准条件下, 夏玉米花期≥38 ℃极端高温发生频率为10.6%, 未来气候情景下增加到19.2%(RCP 4.5)和20.3%(RCP 8.5), 极端高温发生概率增加将面临更严重的绝收风险。本文仅重点分析了≥32 ℃和≥35 ℃两个致灾阈值的发生频率和时空分布特征, 而对产量损失更大的极端高温的影响还有待做进一步详细分析。

受气候情景数据时间分辨率的制约, 未来高温灾害风险的预估是建立在逐日气象数据基础上的。随着高温灾害研究的不断深入, 不同高温持续时间对玉米花粉活力、产量影响的差异也很显著[5], 已细化到以小时为时间尺度的研究层面, 对高温灾害影响的评估提出更精细化的气象数据需求。

本文对夏玉米花期的判定是在假设品种熟性不变的基础上完成的, 随着气候变暖, 生产上可通过改种偏晚熟的品种以充分利用生长季热量资源提高产量。夏玉米花期亦会发生改变, 有可能避过发生极端高温的高峰时段, 降低高温灾害影响。前人研究也指出播期调整是作物调节发育期趋利避害的较为有效措施[37], 通过改变播期调整花期, 躲避高温影响。但花期调整的策略和范围尚不明确, 以未来RCP气候情景下河南省夏玉米主要生长季平均最高气温变化为例(图 9), RCP-rf基准条件下, 向前向后调整花期均有降低风险的可能。RCP 4.5条件下向后调整花期可降低高温风险。RCP 8.5情景下, 向前向后调整播期均不利于躲避高温, 尤其是向前调整, 花期高温风险增加明显。但这仅就最高气温变化的“平均态”而言, 不同地区气候变化的差异性较大, 需要开展更深入的研究, 根据实际情况制定适应对策。

图 9 RCP情景下花期调整对夏玉米躲避高温风险的效果 Fig. 9 Effect of adjusting flowering date on avoiding high temperature risk under RCP scenarios

高温往往伴随着干旱, 研究表明高温风险较高的地区干旱发生频繁, 且高温干旱的协同作用大于二者的简单相加[38]。高温加速了土壤蒸发和作物蒸腾, 使干旱加剧, 玉米受灾减产风险增加。加强水肥调控是防御高温干旱的综合有效措施, 。合理的灌溉可以降低田间温度, 灌溉后玉米能得到充足的水分, 促进蒸腾作用, 降低冠层温度, 增施有机肥、施用微量元素锌肥和后期补充钾肥可以增加玉米耐热性[24], 从而有效减少高温灾害影响。

4 结论

RCP-rf情景下河南省夏玉米花期≥32 ℃高温日数为1.7~5.7 d, 发生频率为20.5%~81.0%。与RCP-rf相比, 未来RCP情景下夏玉米花期≥32 ℃高温发生日数增加0.6 d(RCP 4.5)和0.5 d(RCP 8.5), 发生频率增加9.1%(RCP 4.5)和11.0%(RCP 8.5)。RCP-rf情景下夏玉米花期≥35 ℃高温日数在0.3~3.6 d范围, 发生频率全省在3.9%~51.9%范围。与RCP-rf相比, 未来RCP情景下夏玉米花期≥35 ℃高温发生日数增加0.6 d(RCP 4.5)和0.7d(RCP 8.5), 发生频率增加8.7%(RCP 4.5)和8.3%(RCP 8.5)。RCP-rf情景下全省夏玉米花期≥32 ℃高温积害为48.5~200.9 ℃·d, ≥35 ℃高温积害为9.8~138.5 ℃·d。与RCP-rf相比, 未来RCP情景下≥32 ℃高温积害增加25.4

℃·d(RCP 4.5)和25.6 ℃·d(RCP 8.5); ≥35 ℃高温积害增加25.8 ℃·d(RCP 4.5)和31.4 ℃·d(RCP 8.5)。

夏玉米花期高温综合风险分布可知, RCP-rf情景下高值风险区主要分布在新乡、郑州、许昌、漯河、周口及其以东以北的地区(商丘除外), 约占夏玉米主栽区面积的30.1%;未来情景下高值风险区分布面积扩大至洛阳和南阳以东的大部分地区, 约占夏玉米主栽区面积的63.4%(RCP 4.5)和76.3%(RCP 4.5), 较基准条件分别增加33.3%(RCP 4.5)和46.2%(RCP 4.5), 未来气候情景下夏玉米花期高温灾害风险增加。

致谢: 感谢国家气候中心提供的利用区域气候模式所进行的中国区域未来气候变化模拟结果。

参考文献
[1]
国家统计局河南调查总队. 河南调查年鉴-2017[M]. 北京: 中国统计出版社, 2017.
Survey Office of the National Bureau of Statistics in Henan. Henan Survey Yearbook[M]. Beijing: China Statistics Press, 2017.
[2]
李金才, 崔彦宏, 董海荣, 等. 玉米花丝生长发育及受精结实能力研究进展[J]. 河北农业大学学报, 2002, 25(1): 86-89.
LI J C, CUI Y H, DONG H R, et al. The advances on the growth and development and the fertilization capability of maize (Zea mays L.) silk[J]. Journal of Agricultural University of Hebei, 2002, 25(1): 86-89.
[3]
赵龙飞, 李潮海, 刘天学, 等. 玉米花期高温响应的基因型差异及其生理机制[J]. 作物学报, 2012, 38(5): 857-864.
ZHAO L F, LI C H, LIU T X, et al. Genotypic responses and physiological mechanisms of maize (Zea mays L.) to high temperature stress during flowering[J]. Acta Agronomica Sinica, 2012, 38(5): 857-864.
[4]
陶志强, 陈源泉, 隋鹏, 等. 华北春玉米高温胁迫影响机理及其技术应对探讨[J]. 中国农业大学学报, 2013, 18(4): 20-27.
TAO Z Q, CHEN Y Q, SUI P, et al. Effects of high temperature stress on spring maize and its technologic solutions in North China Plain[J]. Journal of China Agricultural University, 2013, 18(4): 20-27.
[5]
降志兵, 陶洪斌, 吴拓, 等. 高温对玉米花粉活力的影响[J]. 中国农业大学学报, 2016, 21(3): 25-29.
JIANG Z B, TAO H B, WU T, et al. Effects of high temperature on maize pollen viability[J]. Journal of China Agricultural University, 2016, 21(3): 25-29.
[6]
杨国虎. 玉米花粉花丝耐热性研究进展[J]. 种子, 2005, 24(2): 47-51.
YANG G H. Study on the progress of heat tolerance in maize pollen and silk[J]. Seed, 2005, 24(2): 47-51. DOI:10.3969/j.issn.1001-4705.2005.02.013
[7]
张保仁, 董树亭, 胡昌浩, 等. 高温对玉米籽粒淀粉合成及产量的影响[J]. 作物学报, 2007, 33(1): 38-42.
ZHANG B R, DONG S T, HU C H, et al. Effect of high air temperature during different growth stage on starch synthesis in grain and yield in maize (Zea mays L.)[J]. Acta Agronomica Sinica, 2007, 33(1): 38-42. DOI:10.3321/j.issn:0496-3490.2007.01.007
[8]
赵龙飞, 李潮海, 刘天学, 等. 花期前后高温对不同基因型玉米光合特性及产量和品质的影响[J]. 中国农业科学, 2012, 45(23): 4947-4958.
ZHAO L F, LI C H, LIU T X, et al. Effect of high temperature during flowering on photosynthetic characteristics and grain yield and quality of different genotypes of maize (Zea Mays L.)[J]. Scientia Agricultura Sinica, 2012, 45(23): 4947-4958. DOI:10.3864/j.issn.0578-1752.2012.23.023
[9]
LIU Z J, YANG X G, HUBBARD K G, et al. Maize potential yields and yield gaps in the changing climate of northeast China[J]. Global Change Biology, 2012, 18(11): 3441-3454. DOI:10.1111/j.1365-2486.2012.02774.x
[10]
TEIXEIRA E I, FISCHER G, VAN VELTHUIZEN H, et al. Global hot-spots of heat stress on agricultural crops due to climate change[J]. Agricultural and Forest Meteorology, 2013, 170: 206-215. DOI:10.1016/j.agrformet.2011.09.002
[11]
REZAEI E E, WEBBER H, GAISER T, et al. Heat stress in cereals:Mechanisms and modelling[J]. European Journal of Agronomy, 2015, 64: 98-113. DOI:10.1016/j.eja.2014.10.003
[12]
姚凤梅, 张佳华. 1981-2000年水稻生长季相对极端高温事件及其气候风险的变化[J]. 自然灾害学报, 2009, 18(4): 37-42.
YAO F M, ZHANG J H. Change of relative extreme high temperature events and climate risk in rice growing period in China from 1981 to 2000[J]. Journal of Natural Disasters, 2009, 18(4): 37-42. DOI:10.3969/j.issn.1004-4574.2009.04.007
[13]
LIZASO J I, RUIZ-RAMOS M, RODRÍGUEZ L, et al. Impact of high temperatures in maize:Phenology and yield components[J]. Field Crops Research, 2018, 216: 129-140. DOI:10.1016/j.fcr.2017.11.013
[14]
TAO Z Q, CHEN Y Q, LI C, et al. The causes and impacts for heat stress in spring maize during grain filling in the North China Plain-A review[J]. Journal of Integrative Agriculture, 2016, 15(12): 2677-2687. DOI:10.1016/S2095-3119(16)61409-0
[15]
赵丽晓, 雷鸣, 王璞, 等. 花期高温对玉米子粒发育和产量的影响[J]. 作物杂志, 2014(4): 6-9.
ZHAO L X, LEI M, WANG P, et al. Effects of high temperature stress during flowering on maize kernel development and grain yield[J]. Crops, 2014(4): 6-9.
[16]
于康珂, 刘源, 李亚明, 等. 玉米花期耐高温品种的筛选与综合评价[J]. 玉米科学, 2016, 24(2): 62-71.
YU K K, LIU Y, LI Y M, et al. Screening and comprehensive evaluation of heat-tolerance of maize hybrids in flowering stage[J]. Journal of Maize Sciences, 2016, 24(2): 62-71.
[17]
陈朝辉, 王安乐, 王娇娟, 等. 高温对玉米生产的危害及防御措施[J]. 作物杂志, 2008(4): 90-92.
CHEN Z H, WANG A L, WANG J J, et al. Influence of high temperature on growth and development of maize[J]. Crops, 2008(4): 90-92. DOI:10.3969/j.issn.1001-7283.2008.04.029
[18]
徐翠莲, 黄晓书, 王彩霞, 等. 玉米花粉的采集、干燥和贮藏的初步研究[J]. 河南农业科学, 1996(8): 8-10.
XU C L, HUANG X S, WANG C X, et al. The study on collection, dryness and storage of corn pollen[J]. Journal of Henan Agricultural Sciences, 1996(8): 8-10.
[19]
王海梅. 高温胁迫对河套灌区玉米生理指标及产量构成要素的影响[J]. 干旱气象, 2015, 33(1): 59-62.
WANG H M. Influence of high temperature stress on physiological indexes and yield components of maize in Hetao irrigation district[J]. Journal of Arid Meteorology, 2015, 33(1): 59-62.
[20]
徐美玲. 温度对玉米花丝生活力的影响[J]. 浙江农业科学, 2002(3): 120-122.
XU M L. Impacts of temperature on vitality of pistils in maize[J]. Journal of Zhejiang Agricultural Sciences, 2002(3): 120-122. DOI:10.3969/j.issn.0528-9017.2002.03.007
[21]
李德, 孙义, 孙有丰. 淮北平原夏玉米花期高温热害综合气候指标研究[J]. 中国生态农业学报, 2015, 23(8): 1035-1044.
LI D, SUN Y, SUN Y F. Use of integrated climatic index to determine high temperature damage to summer maize at florescence in the Huaibei Plain[J]. Chinese Journal of Eco-Agriculture, 2015, 23(8): 1035-1044.
[22]
王丽君.黄淮海平原夏玉米季干旱、高温的发生特征及对产量的影响[D].北京: 中国农业大学, 2018
WANG L J. Spatio-temporal characteristics of drought, heat and its effect on yield for summer maize in Huang-Huai-Hai Plain, China[D]. Beijing: China Agriculture University, 2018
[23]
刘哲, 乔红兴, 赵祖亮, 等. 黄淮海夏播玉米花期高温热害空间分布规律研究[J]. 农业机械学报, 2015, 46(7): 272-279.
LIU Z, QIAO H X, ZHAO Z L, et al. Spatial distribution of high temperature stress at corn flowering stage in Huang-Huai-Hai Plain of China[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 272-279.
[24]
刘哲, 汪雪滢, 刘帝佑, 等. 基于MODIS数据的黄淮海夏玉米高温风险空间分布[J]. 农业工程学报, 2018, 34(9): 175-181.
LIU Z, WANG X Y, LIU D Y, et al. Spatial distribution of high temperature risk on summer maize in Huang-Huai-Hai Plain based on MODIS data[J]. Transactions of the CSAE, 2018, 34(9): 175-181.
[25]
尹小刚, 王猛, 孔箐锌, 等. 东北地区高温对玉米生产的影响及对策[J]. 应用生态学报, 2015, 26(1): 186-198.
YIN X G, WANG M, KONG Q X, et al. Impacts of high temperature on maize production and adaptation measures in Northeast China[J]. Chinese Journal of Applied Ecology, 2015, 26(1): 186-198.
[26]
李克让, 陈育峰. 中国全球气候变化影响研究方法的进展[J]. 地理研究, 1999, 18(2): 214-219.
LI K R, CHEN Y F. The progress in methodologies to study impacts of global climate change in China[J]. Geographical Research, 1999, 18(2): 214-219. DOI:10.3321/j.issn:1000-0585.1999.02.015
[27]
赵宗慈. 为IPCC第五次评估报告提供的全球气候模式预估[J]. 气候变化研究进展, 2009, 5(4): 241-243.
ZHAO Z C. Experiments of global climate models proposed for IPCC AR5[J]. Advances in Climate Change Research, 2009, 5(4): 241-243. DOI:10.3969/j.issn.1673-1719.2009.04.010
[28]
王绍武, 罗勇, 赵宗慈, 等. 新一代温室气体排放情景[J]. 气候变化研究进展, 2012, 8(4): 305-307.
WANG S W, LUO Y, ZHAO Z C, et al. New generation of scenarios of greenhouse gas emission[J]. Progressus Inquisitiones de Mutatione Climatis, 2012, 8(4): 305-307. DOI:10.3969/j.issn.1673-1719.2012.04.011
[29]
张倩, 赵艳霞, 王春乙. 我国主要农业气象灾害指标研究进展[J]. 自然灾害学报, 2010, 19(6): 40-54.
ZHANG Q, ZHAO Y X, WANG C Y. Advances in research on major agro-meteorological disaster index in China[J]. Journal of Natural Disasters, 2010, 19(6): 40-54.
[30]
郭建平. 气候变化对中国农业生产的影响研究进展[J]. 应用气象学报, 2015, 26(1): 1-11.
GUO J P. Advances in impacts of climate change on agricultural production in China[J]. Journal of Applied Meteorological Science, 2015, 26(1): 1-11.
[31]
LIU Y, WANG E L, YANG X G, et al. Contributions of climatic and crop varietal changes to crop production in the North China Plain, since 1980s[J]. Global Change Biology, 2010, 16(8): 2287-2299.
[32]
马玉平, 孙琳丽, 俄有浩, 等. 预测未来40年气候变化对我国玉米产量的影响[J]. 应用生态学报, 2015, 26(1): 224-232.
MA Y P, SUN L L, E Y H, et al. Predicting the impact of climate change in the next 40 years on the yield of maize in China[J]. Chinese Journal of Applied Ecology, 2015, 26(1): 224-232.
[33]
LOBELL D B, SCHLENKER W, COSTA-ROBERTS J. Climate trends and global crop production since 1980[J]. Science, 2011, 333(6042): 616-620. DOI:10.1126/science.1204531
[34]
TAO F L YOKOZAWA M, LIU J Y, et al. Climate-crop yield relationships at provincial scales in China and the impacts of recent climate trends[J]. Climate Research, 2008, 38(1): 83-94.
[35]
SCHLENKER W, ROBERTS M J. Nonlinear temperature effects indicate severe damages to U.S. crop yields under climate change[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(37): 15594-15598. DOI:10.1073/pnas.0906865106
[36]
LOBELL D B, BÄNZIGER M, MAGOROKOSHO C, et al. Nonlinear heat effects on African maize as evidenced by historical yield trials[J]. Nature Climate Change, 2011, 1(1): 42-45. DOI:10.1038/nclimate1043
[37]
安盼盼, 明博, 董朋飞, 等. 黄淮南部玉米产量对气候生态条件的响应[J]. 作物学报, 2018, 44(3): 442-453.
AN P P, MING B, DONG P F, et al. Response of maize (Zea mays L.) yield to climatic ecological condition on the South Yellow-Huaihe-Haihe Rivers Plain[J]. Acta Agronomica Sinica, 2018, 44(3): 442-453.
[38]
赵霞, 穆心愿, 马智艳, 等. 不同玉米杂交种对花期高温、干旱复合胁迫的响应[J]. 河南农业科学, 2017, 46(8): 32-37.
ZHAO X, MU X Y, MA Z Y, et al. Response of different maize hybrids to high temperature and drought stresses at flowering stage[J]. Journal of Henan Agricultural Sciences, 2017, 46(8): 32-37.