玉米作为世界第一大粮食作物, 种植面积广泛、生产潜力巨大, 其在粮食安全及国民经济发展中占据重要地位^[1-2]。IPCC第五次评估报告指出, 预计21世纪全球平均气温增幅可能超过1.5—2℃, 以全球气候变暖为主要特征的气候要素变化已对世界环境和经济发展造成了巨大威胁^[3-4]。气候变化条件下, 农业气候资源也相应发生变化^[5-7]。由于农业生产活动对气候条件的依赖程度很高, 任何程度的气候变化都将给农业生产带来显著影响, 特别是作物生长季内水热条件的变化^[8]。玉米生产受气候变化的影响十分显著, 因此, 如何应对气候变化对农作物带来的影响成为当前研究的热点话题。

近年来, 国内外学者关于气候变化对粮食生产的影响进行了广泛研究, 其方法主要是利用不同作物生长模型以及统计学方法研究不同时空尺度、作物类型生育期的温度、降水、日照、蒸散量等气候变量的时空特征及其与农业生产的相关性评价^[9-12]。如何奇瑾等^[13]利用ArcGIS空间分析对我国春玉米潜在种植分布区进行研究, 并进行了玉米种植区划；余会康等^[14]研究表明增温效应对东北地区水稻生育过程产生重要影响, 水稻冷害发生频率逐年减少；薛亚荣等^[15]认为降水是影响塔城地区棉花灌溉需水量最大的气象因子；Golovko等^[16]研究了在相对寒冷潮湿的土壤和气候条件下培育早熟玉米杂交种的可能性和前景；Labell等^[17]认为极端高温对非洲玉米生产影响显著, 日最高温度大于30℃的积温每增加1℃, 将导致玉米减产1.7%；同时, Otkin和Wu等^[18-19]研究认为, 干旱对玉米种植产生巨大影响。综合来看, 以往研究缺乏对作物各生育期内气温和降水时空特征的协同研究, 致使难以明确气候变化背景下的区域水热协同变异过程对作物生育的影响；另外, 对年代际间春玉米水热变化趋势及差异水平的研究鲜有发现, 而这在指导地区农业可持续与稳定发展、中长期规划与种植结构调整中具有潜在价值。

辽宁省作为国家重点粮食产区之一, 在《全国新增1000亿斤粮食生产能力规划(2009—2020年)》中, 被定位为粮食生产核心区^[20]。辽宁省玉米种植面积较广阔, 在辽宁省14个地级市中, 12个以玉米为第一种植作物^[21-22]。近年来, 随着全球气候变化的不断发展, 东北地区气候总体存在以“气温显著增高, 降水、日照时数减少”为主要特点的暖干化趋势, 农业生产活动受到影响^[23]。选取辽宁省1970—2018年近49年春玉米生育期内19个气象站点逐日气象数据, 划分为1970—1979、1980—1989、1990—1999、2000—2009、2010—2018年5个研究区间, 从年代际尺度上分析气候变化背景下辽宁省玉米生育期内气温和降水气候资源的时空变化特征及差异水平, 旨在为优化其农业生产布局, 调整作物种植结构提供理论与决策支持。

1 研究区域概况

辽宁省位于我国东北地区的东南部, 38°43′—43°26′N, 118°53′—125°46′E之间, 南濒黄海、渤海, 总面积为14.8万km²(图 1)。辽宁省是我国重要的粮食主产区之一, 现代农业是辽宁省发展规划的重要部分, 其中玉米种植占据重要地位。辽宁省属温带大陆性季风气候, 且受海洋性气候影响显著；年平均气温7—11℃, 年降水量为600—1100 mm, 时空分布不均, 东部多西部少, 且多集中7—8月。辽宁省气候自东南向西北由温带的半湿润向半干旱过渡, 具有日照丰富、雨热同季、积温较高、降水较少的特征^[24]。近年来, 气候变暖下东北地区因纬度偏高, 增暖明显, 降水量减少, 热量指数和水分亏缺指数随年代呈波动上升趋势, 农业生产受到较大影响^[25]。

2 资料与方法 2.1 数据来源

数据来源于中国气象数据网, 主要包括辽宁省近49年(1970—2018年)春玉米全育期内19个标准气象站点的逐日气象数据(图 2), 主要包括日最高气温, 日平均气温以及日降水量和各个站点的经纬度等地理特征数据。数据分析基于Excel 2010和SPSS软件进行数据的统计分析, 采用ANOVA进行方差分析, 显著性水平P=0.05。空间分析采用ArcGIS 10.0软件的克里金插值方法(Kriging)对气象数据进行插值, 生成空间栅格数据, 按掩膜提取后得到相应的空间分布图。

2.2 研究方法 2.2.1 辽宁省春玉米生育期划分及作物系数K_C

根据实际生产实践经验并结合辽宁省水利电力厅编制的《辽宁省抗旱资料汇编》及相关资料^[26-28], 玉米生产经历播种、出苗、拔节、抽雄、灌浆、乳熟、成熟等主要发育时期, 将辽宁省春玉米全育期划分为5个阶段(表 1), 具体划分时间为播种期、苗期、拔节期、抽穗期、成熟期, 共历时163 d。根据李晶等^[29]对已有对辽宁玉米作物系数的研究, 结合靳春香^[30]对辽宁省春玉米作物系数的总结, 故本文选取K_C(表 2)。

2.2.2 热量指标

气温是影响玉米生长的重要条件之一, 不仅影响作物的发育速度, 同时也影响作物的产量。玉米为喜温作物, 气温的变化会直接影响玉米的整个生长过程, 玉米生长的最低温度为10—12℃, 适宜温度为25—31℃^[31-33]。积温学说或者有效积温学说(≥10℃)是评价作物生长期内热量的有效依据, 目前学者多使用平均气温作为判断热量的主要指标。然而在作物实际生长过程中, 过高或过低的温度都会对作物生长产生不利影响^[34]。为更好地分析玉米生育期内的热量条件, 本文将采用度日指标^[35-37]。

生长度日(Growing degree days, GDD)是分析作物生长期内的热量指标, 文中是指日平均气温≥10℃的有效积温。高温度日(Heat degree days, HDD)是分析作物极端高温的热量指标, 文中是指日平均气温≥31℃的积累温度。本文利用度日指标衡量春玉米全育期及各生育期的温度累积程度, 分析辽宁省春玉米不同生育期受温度累积的可能程度及变化趋势。

(1)

(2)

式中, T_ave代表日平均气温, T_base代表玉米全育期下限温度, 文中取值为10℃；T_max代表日最高气温, T_opt代表玉米生育期的上限温度, 文中取值为31℃。

2.2.3 气候倾向率

气候倾向率是指气候要素多年变化, x_i表示样本量为n的某气候变量, t_i表示x_i对应的时间, 建立x_i与t_i的一元回归方程。

(3)

式中, a指回归常数, b指回归系数, 通常以b的10倍作为各气候要素的气候倾斜率。利用气候倾向率分析1970—2018年5个不同年代际间辽宁省春玉米全育期及各个生育期的GDD、HDD和降水量的变化趋势。

2.2.4 干燥度指数

干燥度指数(Aridity index)是评价一个地区干湿程度的重要指标, 通常以某地区水分收支与热量平衡来计算^[2]。计算公式如下：

(4)

式中, Aridity代表玉米在生育期内的干燥度, P代表玉米在生育期内的相应累积降水量；ET₀代表玉米在生育期内的相应作物蒸散量之和。如果Atidity小于0, 表明该时段里降水量大于蒸散量, 干旱不会发生, 因此可以用0代替所有负值；如果Aridity大于0, 表明该时段里干旱可能发生, 并且Aridity越大, 干旱越严重。

农作物蒸散量通常指农作物在生长期间蒸发蒸腾所消耗的水资源量, 通常采用FAO(联合国粮食及农业组织)推荐的Penman-Monteith模型。

(5)

式中, R_n为净辐射MJ m^-2 d ^-1；G为热通量密度MJ m^-2 d ^-1；T(℃)为日平均温度；u₂为2m高处风速(m/s)；e_s为饱和空气水汽压(kPa)；e_a为空气水汽压(kPa)；Δ为饱和水汽压-温度曲线的梯度(kPa/℃)；γ为湿度计常数；G_n、G_d为固定常数, 在估算逐日潜在蒸散发量时取值分别为900和0.34。

由于不同地区地貌地形条件的不同, 根据不同地区作物系数K_C对蒸散量ET₀进行修订, 得到作物需水量ETc。因此, 采用修订后的干燥度指数公式为：

(6)

(7)

3 结果与分析 3.1 不同年代辽宁省春玉米全育期生长度日GDD时空特征分析

辽宁省不同年代际间春玉米各个生育阶段GDD平均值及气候倾向率见表 3。总体来看, 1970—2018年辽宁省春玉米全育期GDD以44.23℃ · d/10a的速率上升, 不同年代际春玉米生长度日(GDD)在苗期、成熟期以及全育期存在显著差异, 且显著性水平通过0.05检验；表明不同年代际温度在苗期、成熟期以及全育期变化幅度存在明显差异。此外, 各个生育阶段内春玉米GDD呈现递增趋势, 在播种期增长速率最高, 气候倾向率均呈正向变化。其中, 2010—2018年间春玉米GDD在各个生育阶段增长趋势明显大于其他阶段, 将意味着春玉米播种期的提前, 生长期的延长^[38]。同时, 近年来GDD增幅明显, 2000—2018年苗期春玉米GDD明显大于1970—1979年(P < 0.05), 1970—1979年成熟期春玉米GDD明显低于1990—2018年, 2000—2009年GDD值显著高于其他年际(P < 0.05), 2000—2018年全育期内春玉米GDD明显大于1970—1989年(P < 0.05), 表明春玉米全育期生长热量总体在持续增加。

辽宁省春玉米1970—2018年不同年代际全育期GDD空间分布如图 3所示。近49年, 辽宁省春玉米全育期GDD范围在1317.23—1980.53℃ · d之间；GDD空间分布上由西北向东南沿海递增, 局部地区显著增长。GDD高值区范围主要分布在辽宁省西部和中部地区, 低值范围主要分布在辽宁省东部地区, 与王培娟^[39]对东北地区春玉米生育期热量研究的结果一致, 同时与整个东北春玉米冷害危险性指数呈自西南向东北方向递增的空间分布格局相呼应^[40]。由图 3可以看出, 1970—1979年GDD差值不大, 1980—1999年西北部地区GDD增加明显, 2000—2018年辽宁省GDD整体增加, 高值区主要向朝阳、鞍山及营口等地区集中, 其主要原因是朝阳、鞍山及营口等地区温度增大于其他地区, 为0.5—1℃, 其他地区平均增幅多为0—0.5℃^[41]；低值区主要出现在清原、宽甸及周围地区, 温度的区域分布特点可能受地形及海洋影响^[42]。总体上, GDD呈现上升趋势, 空间分布上由西北向东南扩增, 同时地区间差异逐渐减小, 辽宁省春玉米全育期的热量增加有利于春玉米发育及产量增加。

3.2 不同年代辽宁省春玉米全育期高温度日HDD时空特征分析

辽宁省不同年代际间春玉米各个生育阶段HDD平均值及气候倾向率见表 4。总体来看, 1970—2018年辽宁省春玉米全育期HDD以4.77 ℃ · d/10a的速率上升, 在2010—2018年变化速率最为明显, 且不同年代际间春玉米高温度日(HDD)在相同生育阶段内未达到显著差异水平, 表明HDD在不同年代际不同生育期虽有差异, 但差异显著性水平未达到0.05。整体上春玉米HDD在拔节期的均值明显大于其他生育阶段, 在2000—2009年表现尤为显著, 其次为成熟期。同时, 从各个生育阶段来, 近49年春玉米HDD增长期主要在抽穗期和拔节期, 在抽穗期HDD的增长最为明显。据历史资料显示, 2000—2009年, 辽宁省干旱发生频繁, 2000, 2001年、2007年、2009年均出现了较大旱情, 可能与气温升高、降水减少有关。

辽宁省春玉米1970—2018年不同年代际全育期高温度日HDD空间分布如图 4所示。近49年来, HDD在不同年代际的范围变化幅度明显, 在空间上整体表现出不断增大的趋势, 辽西北地区增长显著；在空间变化上由西北地区向东南沿海增大, 其主要原因可能是西部、中部地区高温日数普遍于沿海和东部山区^[43]。1970—1999年间HDD增长趋势较2000年后的增长趋势低；2000—2018年HDD增长趋势明显上升, 其中2010—2018年间HDD显著高于其他阶段；高值区主要在朝阳和叶柏寿地区集中, 可能与其降水少、蒸发大、植被少等因素有关^[24], 其次在辽宁中部地区, 主要是鞍山一带, 与其地区气温增加幅度显著有关；整体来看, 以朝阳为中心的辽西北地区高值区变化范围明显, 80年代末开始突变的气候变暖持续影响至今, 朝阳地区气温整体呈现波动上升趋势且最高和最低增温的强弱在时间上表现出不一致性^[44]。HDD的增长不仅表明春玉米全育期内热量的增加, 也意味着蒸发加剧, 作物生育期水分亏缺指数增加, 高温干旱的风险将大大增强^[2]。总体上, 近49年来HDD呈现增长趋势, 在空间上与生长度日GDD分布趋势大体一致。

3.3 不同年代辽宁省春玉米全育期降水量时空特征分析

辽宁省不同年代际间春玉米各个生育阶段降水量平均值及气候倾向率见表 5。1970—2018年辽宁省春玉米全育期降水量以4.93 mm/10 a的速率下降, 年际变化大, 波动性强；且均值间差异没有达到显著水平, 表明降水量在不同年代际不同生育期差异显著性水平未达到0.05。其中, 2000—2009年降水量明显低于其他年代, 减少时期主要在抽穗期和成熟期, 春玉米全育期内降水量在2010—2018年变化速率最为明显；从不同玉米生育期阶段来, 近49年春玉米生育期降水量增长主要表现在播种期和苗期, 而春玉米在拔节期、抽穗期以及成熟期的降水量存在不同程度的降低, 随着降水减少、气温升高将对作物生物量积累和产量形成都有不利影响。

辽宁省春玉米1970—2018年不同年代际全育期降水量空间分布如图 5所示。近49年, 辽宁省春玉米全育期降水量在不同年代际间变率大, 空间上呈东南—西北递减趋势, 其中中部地区年降水量变化较大, 而西部和东部变化较小^[45], 主要是487.56—628.35 mm降水量的空间分布变化较大。辽宁省春玉米降水量在346.75—909.94 mm, 高值区主要分布在宽甸和丹东等东南地区, 低值区主要分布叶柏寿和朝阳等西北地区, 主要与辽宁省地形有关, 西北地区受海洋性气候影响较小, 气温高, 蒸发旺盛^{[40, 46]}。不同年代际间辽宁省春玉米生育期年均降水量空间分布明显变化, 地区间存在明显差异, 与韩丛利对辽宁省洪涝灾害时间变化特征分析中提到的降水在时空上分布差异的研究结果一致^[47]；1990—1999年辽宁省春玉米全育期内整体降水量高于其他年代际；而2000—2009年降水量显著减少。总体而言, 近49年辽宁省春玉米生育期降水量空间分布差异不明显, 但与辽宁省春玉米生育期生长度日GDD、高温度日HDD的空间分布呈相反趋势。

4 讨论

以温度升高和降水减少为主要标志的辽宁省气候暖干化现象对当地农业产生了重要影响, 高温和干旱是限制农作物生产的重要因子, 其危害的严重度很大程度上取决于温度^[48-49]。由于气温变化是导致作物生育期内活动积温变化的直接原因, 因此本文的热量条件重点考虑了温度, 并未考虑日照辐射的影响。本文通过对辽宁省春玉米全育期GDD、HDD以及降水量在年代际上的时空变化特征及差异水平研究发现, 辽宁省春玉米全育期GDD、HDD呈不断上升的趋势, 虽然未来春玉米生产能够获得良好的热量资源, 但叶柏寿、朝阳及鞍山等地极端高温、干旱的风险极有可能大大增加, 气温升高会加剧土壤蒸发, 引起土壤干旱, 从而严重影响玉米产量^[50]。此外, 辽宁省春玉米全育期内降水量呈现降低趋势, 且西北部降水量较少, 未来春玉米干旱风险加大。

4.1 不同年代际辽宁省春玉米干燥度指数时空特征分析

农业作为受干旱影响最显著的部门之一, 由于气温升高或者降水减少, 作物不能从土壤或者空气中获得生长发育足够的水分, 造成农作物生长发育不良, 导致作物减产或者品质下降的一种灾害现象^[51]。前人对辽宁省玉米干旱研究颇多, 但多集中在以Penman-Monteith、Z指数、标准化降水指数(SPI)、帕默尔干旱指数(PDSI)、综合气象干旱指数(CI)等方法为基础上的干旱风险研究^{[27, 52]}, 通过重点研究气温和降水组合条件下辽宁省春玉米生产的气候条件, 虽然春玉米全育期的热量资源增加, 但降水量在减少, 表明辽宁省春玉米全育期水热时空变化给玉米生产带来了极端高温、干旱等风险变化。在此背景下引入干燥度指数(Aridity index), 对气候变化下辽宁省春玉米在1970—2018年全育期发生干旱的风险性进行研究。

利用修订后干燥度指数公式分析1970—2018年5个年代际辽宁省春玉米全育期干旱发生的时空变化特征, 干燥度指数时空变化分布如图 6所示。近49年来, 辽宁省春玉米全育期在不同年代际干燥度指数空间分布不同, 总体上由东南地区向西北地区增大, 其中, 东部地区整体变化不明显, 西部及中部地区在不同年代际分布不同；此外, 辽西北地区干燥度指数最为显著, 高值中心出现在朝阳、叶柏寿等地及周边地区。其中, 1970—1979年、2000—2009年干燥度指数空间分布特征基本一致；1980—1989年、2010—2018年干燥度指数空间分布特征基本一致；1990—1999年干燥度指数整体较低。据史料记载, 1994、1997、2001以及2014年辽宁中部地区发生不同程度干旱事件, 2000年辽宁省出现大面积干旱事件, 2007年遭受了严重夏旱, 2009年辽西北发生特大干旱事件^[53-54]。结合GDD、HDD以及降水量的时空变化特征, 1970—1989年间辽宁省春玉米全育期干燥度指数变化不大, 高值区主要集中在西北地区。其中, 1980—1989年瓦房店以南地区变化明显, 大连及周边地区变化显著。在此阶段, 西北地区GDD、HDD增高, 降水量整体上变化不显著。1990—1999年整体干燥度指数为近49年最低。在此阶段, 西北地区GDD增大, HDD由西北地区向东北增加, 降水量整体由东南向西北增加, 为近49年最多, 因此干燥度指数较低。2000—2018年间干燥度指数较1990—1999年明显增加, 仍以辽西北地区最为显著。在此阶段, GDD以及HDD总体由西北向东南方向继续增加, 2000—2009年降水量明显减少, 2010—2018年降水量较上个阶段有所增加, 干燥度指数明显增加。

综上所述, 辽宁省春玉米全育期干燥度指数在不同年代际存在差异, 但GDD、HDD的持续增加并未使得干旱风险范围发生明显变化；同时, 降水量大范围的减少, 干燥度指数会产生一定变化, 结果表明GDD、HDD的升高以及降水量的减少可能带来干旱风险；总体而言, 辽宁省西北部地春玉米发生干旱风险性最高。

4.2 干燥度指数与GDD、HDD以及降水量的聚合分析

散点图可以清晰展示数据的分布和聚合情况。1970—2018年辽宁省春玉米全育期干燥度指数和GDD、HDD以及降水量的散点分布图如图 7所示。由图 7可知, 干燥度指数和GDD没有明显聚合分布, GDD的升高不能引起干燥度指数的明显变化, 结果表明辽宁省春玉米全育期GDD的增长引发干旱风险的可能性较低。干燥度指数和HDD呈现一定聚合分布, 随着HDD的增大, 干燥度指数随之增大, 结果表明辽宁省春玉米全育期HDD的增高会造成一定程度的干旱事件。干燥度指数和降水量的聚合分布程度较高, 降水量的减少会引起干燥度指数的明显变化, 结果表明辽宁省春玉米全育期降水量的减少引发干旱的风险较大。

综上所述, 辽宁省春玉米全育期GDD的增加引发干旱的风险性较低, 而HDD的增长以及降水量的减少对干旱发生的影响较大。因此, 辽宁省应通过优化作物布局, 调整种植结构, 培育栽种技术以及改进管理模式等措施应对干旱风险, 例如, 调整播种期, 培育新品种, 改进灌溉技术等。而未来辽宁省春玉米如何应对气候变化, 降低干旱风险的对策与措施还需要进一步深入研究。

5 结论

(1) 近49年春玉米全育期GDD以44.23 ℃ · d/10a的速率上升, 范围在1317.29—1980.53 ℃ · d之间波动, 同时苗期、成熟期及全育期GDD在不同年代际间存在显著差异；在空间分布上, 1979—1999年间地区差异不大；从2000年开始, 地区间差异逐渐缩小, 高值区向朝阳、鞍山等地集中, 低值区以清原、宽甸地区为中心。良好的热量条件的增加有利于玉米产量和品质的提升。

(2) 近49年春玉米全育期HDD以4.77 ℃ · d/10a的速率上升, 2000—2018年增长趋势明显, 增长期主要集中在抽穗期；在空间分布上, 高值区主要集中在叶柏寿、朝阳、鞍山等地, 整体上以朝阳为中心的辽西北地区高值区变化范围明显。HDD的增加将不利于春玉米生长, 辽西北地区干旱风险发生可能性增加。

(3) 近49年辽宁省春玉米全育期降水量整体以4.93 mm/10a的减少趋势, 各个年代际间变化明显, 拔节期、抽穗期及成熟期降水年际变化大, 波动性强；空间上, 降水量由东南沿海向西北方向递减, 其中中部地区年降水量变化较大, 不同年代际间地区差异明显, 1990—1999年降水量分布变化较大。

(4) 辽宁省春玉米全育期干燥度指数在不同年代际存在差异。从空间分布上, GDD、HDD的对干旱风险范围变化不明显, 但降水量大范围的减少, 干燥度指数会产生一定变化；从干燥度指数与GDD、HDD以及降水量的聚合分析看, 辽宁省春玉米生育期GDD的升高对干旱风险影响较小, 而HDD的升高和降水量的减少引发干旱的风险性较高。

(5) 虽辽宁省春玉米全育期GDD、HDD近49年来呈现不断上升的趋势, 高值区多集中在朝阳、叶柏寿等西部地区；但降水量呈现降低趋势, 且西北部降水量较少, 未来西部地区春玉米干旱风险加大。下一步有必要针对农业气候变化及其对农业影响等开展分区研究, 在相对干旱地区, 可以通过扩大旱田水浇地面积防御干旱威胁, 减轻干旱影响；未来应重点加强对辽宁省春玉米生产气象条件的监测, 同时提高玉米生产用水利用率和利用效率, 降低干旱带来的危害。