旱地春小麦产量对不同生育阶段温度变化的响应模拟

奥鹏网院作业

旱地春小麦产量对不同生育阶段温度变化的响应模拟
王钧1, 2，李广1，闫丽娟3，刘强2，聂志刚2

（1甘肃农业大学林学院，兰州730070；2甘肃农业大学信息科学技术学院，兰州 730070；3甘肃农业大学农学院，兰州 730070）

摘要：【目的】定量分析陇中黄土丘陵沟壑区旱地春小麦不同生育阶段温度变化对小麦产量的影响，为陇中黄土丘陵沟壑区制定合理种植模式提供理论指导和决策依据。【方法】基于灰色关联度分析不同生育阶段温度与旱地春小麦产量间的关联性，确定影响旱地春小麦产量的关键生育阶段。应用APSIM（Agricultural Production System Simulator）模型模拟不同生育阶段温度变化条件下春小麦的产量，通过二次多项回归分析、单因素边际效应分析方法研究旱地春小麦产量对不同生育阶段温度变化的响应机制。【结果】（1）APSIM模型对模拟陇中黄土丘陵沟壑区不同播种期旱地春小麦的产量和生育期具有较高的适用性，模拟产量和生育期的均方根误差（RMSE）平均值分别为39.95 kg·hm-1和2.78 d，归一化均方根误差（NRMSE）平均值分别为1.55%和1.87%，模型有效性指数（ME）平均值分别为0.73和0.83；（2）陇中黄土丘陵沟壑区旱地春小麦全生育期平均温度总体呈上升趋势，但不同生育阶段的温度增幅各不相同，播种—出苗、出苗—分蘖、分蘖—拔节、拔节—孕穗、孕穗—开花、开花—灌浆和灌浆—成熟阶段的平均气温倾向率分别为0.44℃·（10a）-1、0.34℃·（10a）-1、0.17℃·（10a）-1、0.41℃·（10a）-1、0.49℃·（10a）-1、0.52℃·（10a）-1和0.35℃·（10a）-1；（3）不同生育阶段温度变化对旱地春小麦产量的影响大小依次为灌浆—成熟、开花—灌浆、播种—出苗、孕穗—开花、拔节—孕穗、出苗—分蘖、休闲和分蘖—拔节阶段；（4）在其他生育阶段温度不变的条件下，播种—出苗阶段温度每提高0.5℃，春小麦产量平均增加0.45%；孕穗—开花阶段温度每提高0.5℃，春小麦产量平均减少0.34%。开花—灌浆阶段温度每提高0.5℃，春小麦产量平均减少0.65%；灌浆—成熟阶段温度每增加0.5℃，春小麦产量平均减少1.09%。【结论】APSIM模型对陇中黄土丘陵沟壑区不同播种期旱地春小麦的产量和主要生育期有较好的模拟效果。研究区旱地春小麦全生育期平均温度总体呈上升趋势，但不同生育阶段温度增幅不同。不同生育阶段温度变化对春小麦生长和发育影响不同。播种—出苗阶段增温倾向于增产，其他阶段增温倾向于减产。

关键词：春小麦产量；不同生育时期；温度变化；APSIM模型；灰色关联度分析

0 引言
【研究意义】全球气候变化是21世纪人类面临最为严重的环境问题之一。全球气候变化主要表现于温度上升、降水格局改变以及极端气候事件增加等，IPCC预测21世纪全球温度升幅可能超过1.5—2.5℃ [1-2]。我国陇中黄土丘陵沟壑区是全球增温最为显著的地区之一，也是对气候变化响应最为脆弱地区之一[3]。小麦作为陇中黄土丘陵沟壑区的主要粮食作物，近年来气候变暖对其生产影响较为严重，尤其对过度依赖自然条件的旱作小麦[4-5]。气候变暖对小麦生长各发育阶段的影响并不完全一致，同时各阶段气候变化对产量及产量构成要素的影响也存在着差异。因此，阐明和量化不同生育阶段温度变化对陇中黄土丘陵沟壑区旱地春小麦产量的影响程度具有重要的研究意义。【前人研究进展】目前，国内外大量研究表明，温度变化对小麦生长的影响十分复杂，对作物不同生长发育阶段产生不同的影响[6-11]。高美玲等[12]利用Meta回归分析了田间增温对小麦产量及生育期持续时间的影响；董莉霞等[13]利用APSIM模型研究了逐日最低温度和最高温度变化对陇中旱地春小麦产量的影响；任新庄等[14]利用APSIM模型分析了陇中旱地春小麦产量对降水与温度变化的响应；张凯等[15-16]利用大田开放式红外增温装置，模拟增温对半干旱雨养区春小麦物质生产与分配的影响，以及抽穗—乳熟阶段增温对陇中黄土高原春小麦产量的影响。PRADHAN等[17]研究了出苗到开花期遭遇高温对春小麦产量的影响。【本研究切入点】陇中黄土丘陵沟壑区旱地春小麦不同生育阶段的温度变化趋势差异显著，且不同生育阶段温度变化对旱地春小麦产量影响大小也各不相同。研究不同生育阶段温度变化对旱地春小麦生长过程的影响有助于进一步揭示温度变化对该地区旱地春小麦生产影响的机理[18]。【拟解决的关键问题】本研究利用APSIM模型模拟陇中黄土丘陵沟壑区不同播期旱地春小麦的产量和生育期，以验证模型对模拟温度变化条件下旱地小麦生长的适用性，基于灰色关联度分析不同生育阶段温度变化与旱地春小麦产量的关联性，从而确定影响小麦产量的主要生育阶段，通过APSIM（Agricultural Production System Simulator）模型模拟分析旱地春小麦4个主要生育阶段11种温度变化水平对应的春小麦产量，并运用二次多项式逐步回归分析、单因素边际效应分析等研究方法，定量分析不同生育阶段温度变化对春小麦产量的影响，从而为陇中黄土丘陵沟壑区合理种植模式的制定提供理论指导和决策依据。

1 材料与方法
1.1 研究区概况
试验地位于甘肃省定西市安定区安家坡乡，该地区属于典型黄土高原丘陵沟壑区，平均海拔高度为2 000 m，多年平均降雨量为395 mm，年均蒸发量为1 540 mm，年平均气温为6.5℃，年均≥0℃积温为2 935.1℃，年均≥10℃积温为2 240.2 ℃，年均太阳辐射为141.6 kJ·cm-2，为典型的半干旱雨养农业区。

1.2 试验设计
为验证不同温度条件下APSIM模型模拟陇中黄土丘陵沟壑区旱地春小麦生育期和产量的适用性，在2016—2018年进行大田试验，试验共设3个播期，分别为早播（ESW，每年3月3日播种）、正常播种（NSW，每年3月18日播种），晚播（LSW，每年3月31日播种），小麦品种为‘定西35号’。试验小区面积为20 m×4 m，采用传统耕作方式，播种深度为30 cm，每个小区重复3次，采用随机区组排列。每个小区各施纯N 105 kg·hm-2，P2O5 105 kg·hm-2。

1.3 APSIM模型
APSIM是澳大利亚农业生产系统研究组（APSRU）开发的农业生产系统模拟模型，该模型具有较强的地域适用性和农业机理基础，对作物种植制度、作物生理生态机理等具有较好的模拟能力。APSIM模型采用模块化管理模式，主要由土壤参数、作物参数和管理参数3个模块组成，能够准确预测不同气候、品种、土壤和管理因素下作物的产量形成、土壤温度、土壤水分等指标的变化过程[19]；在不同气候带、不同作物品种、多种土壤和管理措施条件下，准确地模拟农业系统中各组分的变化情况[20]；同时模型还可以模拟轮作体系下不同耕作措施的土壤温度、土壤含水量、生育进程、降水量、土壤氮素等之间的关系，模拟各种气候背景组合条件下的作物产量和产量构成因素，为全球气候变化背景下的决策、耕作措施的合理制定提供可靠的科学依据[21]。

APSIM模型所需研究区气象数据来源于甘肃省气象局提供的定西市安定区（1970—2018年）历史气象资料。作物品种参数参考李广等[22]关于APSIM在陇中黄土丘陵沟壑区适用性研究中的实测数据（表1）。土壤参数为研究区田间的实测数据（表2）。APSIM模型的管理参数根据研究区的实际管理措施进行设定。

1.4 生育期预测模型
本研究利用APIM模型中的调整积温（TT’）来计算旱地春小麦的生育期，当调整积温达到该生育期阶段的目标积温后，小麦生长就进入了下一个新的生育期。其中，调整积温为某一生育期开始到结束时每日积温在受到光周期、春化和其他环境因子影响后的总和。调整积温计算公式如下：

pagenumber_ebook=56,pagenumber_book=906
式中，△TT为日积温，采用三基点温度进行计算，三基点温度分别为0℃、26℃和34℃；fD和fV分别为光周期因子和春化因子；fW、fP和fN分别为土壤水分、氮和磷含量对小麦生育期的限制因素，其值均为1。

其中，播种至发芽阶段小麦的生长速度主要受温度和土壤水分的影响；发芽至出苗阶段小麦的生长速度主要受播种深度、温度的影响；出苗至分蘖阶段小麦的生长速度主要受温度和春化因子fV（公式2）的影响；分蘖至孕穗阶段小麦的生长速度主要受光周期因子fD（公式3）、温度和小麦品种特性的影响；抽穗至开花阶段小麦的生长速度主要受温度的影响；开花至成熟阶段小麦的生长速度受温度和小麦品种特性（灌浆因子）的影响。

表1 “定西35号”小麦生长发育参数
Table 1 Growth parameters of wheat cultivar of “Dingxi 35”

pagenumber_ebook=56,pagenumber_book=906
表2 研究区土壤属性参数
Table 2 Soil property parameters in the experiment site

pagenumber_ebook=57,pagenumber_book=907
pagenumber_ebook=57,pagenumber_book=907
式中，RV为春化敏感因子；RP为光周期敏感因子；LP为日长，由日序、当地纬度（35°58′）和晨昏蒙影（-6°）根据标准的天文公式计算而来；V为总春化指数，其值为日春化指数（△V）减去反春化指数（△Vd）（公式4）。其中，当日最高温低于30℃，且日最低温低于15℃时，日春化指数计算公式如公式5所示；当日最高温大于30℃，且总春化指数小于10时，将发生反春化，反春化指数计算公式如公式6所示。

pagenumber_ebook=57,pagenumber_book=907
1.5 APSIM产量形成模型
APSIM模型模拟的小麦产量主要由小麦籽粒个数、灌浆速率、以及灌浆前干物质积累转化来决定，小麦产量计算如公式7所示[23]。

pagenumber_ebook=57,pagenumber_book=907
pagenumber_ebook=57,pagenumber_book=907
式中，Ywheat为小麦产量（kg·hm-2），Pstem为每公顷茎数，Dg为每茎籽粒重（g），Dgm为每茎最大籽粒重（g），Sgm为最大籽粒大小（g），Qmeal为谷物粉干物质重（g），Ng为单位茎籽粒数，Rp为潜在灌浆速率（(g·grain-1·d-1），hT为日平均温度影响灌浆速率因子，Tmax日最高温度（℃），Tmin日最低温度（℃）,fN,grain为氮素限制灌浆因子，Rg为单位茎籽粒数(grain·g-1·stem-1)，Qstem为开花时茎干物质重（g），Q为开花时小麦干物质重（g），I为辐射截获，RUE为辐射利用效率，fs为胁迫因子，fc为碳因子，Qleaf为开花时叶干物质重（g），Qhead开花时小麦果实干物质重（g），fT,photo为温度限制性因子，fN,photo为氮因子，hN,poten为潜在灌浆速率（g·grain-1·d-1），hN,min最小潜在灌浆速率（g·grain-1·d-1），hN,grain为氮限制灌浆因子（APSIM模型中默认为1），CN为茎或叶的氮浓度，CN,crit和CN,min为茎或叶的临界和最小氮浓度。公式（7—15）表明，开花期前，温度是APSIM模型小麦茎干物质形成的限制性因子，由茎干物质重影响小麦籽粒数；进入灌浆期后，温度会影响小麦的灌浆速率，从而影响小麦的籽粒重量。

1.6 灰色关联度分析
本文将研究区1970—2018年各年旱地春小麦不同生育阶段（休闲、播种—出苗、出苗—分蘖、分蘖—拔节、拔节—孕穗、孕穗—开花、开花—灌浆、灌浆—成熟）的年平均温度作为灰色系统的一个因素，应用灰色关联法对不同生育阶段平均温度与旱地春小麦模拟产量的关联度进行综合分析。

1.6.1 确定参考数列和比较数列 以1970—2018年研究区旱地春小麦各年产量作为参考数列：X0={x0(k)|,k=1,2,…,n}，以1970—2018年各年旱地春小麦不同生育阶段的平均温度作为比较数列：Xi={xi(k)|,k=1,2,…,n}，其中i=1,2,3,…,m，m为生育阶段个数。

1.6.2 数据归一化处理 利用公式（16）将原始数据归一化，使数据映射到[0, 1]的范围内。

pagenumber_ebook=58,pagenumber_book=908
式中，min为数据的最小值，max为数据的最大值。

1.6.3 计算灰色关联系数 计算参考数列X0与比较数列Xi对应点的绝对差值，即△i(k)=|x0(k)-xi(k)|。ξi(k)为参考数列X0和比较数列Xi在第k个评价指标上的关联系数（公式17）。

pagenumber_ebook=58,pagenumber_book=908
式中，minimink|x0(k)-xi(k)|和maximaxk|x0(k)-xi(k)|分别为关联系数的二级最小差和最大差，λ为分辨系数，λ∈[0,1]，本文取λ=0.5。

1.6.4 计算灰色关联度 计算灰色关联度（公式18），在求出灰色关联度后，按灰色关联度大小排序，该值越大表明参考数列和比较数列的关联度越高，比较数列越接近于参考数列。

pagenumber_ebook=58,pagenumber_book=908
式中，n为年份总数，Ri为灰色关联度。

1.7 数据处理
本研究通过2016—2018年3个播期的大田试验结果，对不同温度条件下APSIM模型模拟旱地春小麦生育期和产量的适用性进行验证，并通过灰色关联度分析确定旱地春小麦休闲、播种—出苗、出苗—分蘖、分蘖—拔节、拔节—孕穗、孕穗—开花、开花—灌浆、灌浆—成熟阶段温度变化与旱地春小麦产量的关联度。基于灰色关联度分析结果，根据IPCC未来全球平均温度升幅2℃的控制目标，在±2.5℃的范围内将4个与小麦产量关联度最大生育阶段的逐日温度降低或者增加0.5℃作为1个试验处理的气象数据，基于APSIM模型利用调整后的气象数据模拟1970—2018年不同气象条件下的小麦产量，利用二次多项回归分析、单因素边际效应分析方法定量分析旱地春小麦产量对不同生育阶段温度变化的响应。

2 结果
2.1 APSIM模型的适应性分析
2.1.1 APSIM模型模拟旱地春小麦生育期的适应性分析 出苗期、孕穗期、开花期和灌浆期是旱地春小麦生长过程中的关键生育期，因此本研究利用2016—2018年研究区3种不同播期条件下旱地春小麦的出苗期、孕穗期、开花期和灌浆期的实测数据对生育期的模拟结果进行对比分析。在3种不同播期条件下，旱地春小麦的4个主要生育期的模拟值与实测值均匀的分布在1∶1线两侧，表现出良好的一致性（图1）。

pagenumber_ebook=58,pagenumber_book=908
图1 春小麦4个主要生育期模拟值与实测值之间的关系
Fig.1 Relationship of simulated and observed values of four major phenology dates of spring wheat

为了进一步验证3种不同播期条件下旱地春小麦出苗期、孕穗期、开花期和灌浆期的模拟值与实测值的拟合程度，本研究选取均方根误差（RMSE）、归一化均方根误差（NRMSE）和模型有效性指数（ME）3种统计指标对模拟结果进行分析（表3）。结果显示，APSIM模型对3种不同播期的旱地春小麦4个主要生育期的模拟效果均较好，出苗期的RMSE为2.71—3.46 d，NRMSE为2.53%—3.71%，ME为0.68—0.83；孕穗期的RMSE为2.52—3.11 d，NRMSE为1.55%—1.81%，ME为0.84—0.90；开花期的RMSE为2.52—3.00 d，NRMSE为1.41%—1.62%，ME为0.80—0.84；灌浆期的RMSE为2.00—2.71 d，NRMSE为1.09%—1.41%，ME为0.83—0.92，相比之下模型对灌浆期的预测精度相对较高。

表3 春小麦4个主要生育期模拟结果检验
Table 3 Test results of simulation on four major phenology dates of spring wheat

pagenumber_ebook=59,pagenumber_book=909
2.1.2 APSIM模型模拟旱地春小麦产量的适应性分析 运用APSIM模型模拟2016—2018年3种不同播期的小麦产量，根据实测产量数据和模拟产量数据对模型模拟结果的有效性进行检验。结果显示，APSIM模型对早播旱地春小麦产量模拟结果的RMSE为45.08 kg·hm-1，NRMSE为1.59%，ME为0.71；对正常播旱地春小麦产量模拟结果的RMSE为39.15 kg·hm-1，NRMSE为1.54%，ME为0.73；对晚播旱地春小麦产量模拟结果的RMSE为35.61 kg·hm-1，NRMSE为1.53%，ME为0.74。研究结果表明，APSIM模拟对不同播期旱地春小麦产量的模拟具有较高的准确性（表4）。

表4 模拟产量与实测产量的相关性
Table 4 Relationships of simulated and observed values of yield

pagenumber_ebook=59,pagenumber_book=909
2.2 春小麦不同生育阶段平均气温的年际变化分析
从1970—2018年陇中黄土丘陵沟壑区旱地春小麦不同生育阶段平均温度的年际变化趋势，可见播种—出苗、出苗—分蘖、分蘖—拔节、拔节—孕穗、孕穗—开花、开花—灌浆和灌浆—成熟阶段的平均气温倾向率分别为0.44、0.34、0.17、0.41、0.49、0.52和0.35℃·(10a)-1（图2）。研究结果表明，陇中黄土丘陵沟壑区旱地春小麦各生育阶段平均温度总体呈上升趋势，但不同生育阶段的温度变化幅度各不相同，各生育阶段温度增幅大小依次为开花—灌浆、孕穗—开花、播种—出苗、拔节—孕穗、灌浆—成熟、出苗—分蘖和分蘖—拔节。

2.3 旱地春小麦产量对不同生育阶段温度的响应
不同生育阶段温度对春小麦生长发育和产量的影响程度不一，为了研究不同生育阶段温度对陇中黄土丘陵沟壑区旱地春小麦生产的影响程度，本文以旱地春小麦休闲、播种—出苗、出苗—分蘖、分蘖—拔节、拔节—孕穗、孕穗—开花、开花—灌浆和灌浆—成熟阶段的平均气温和小麦模拟产量为研究对象，建立旱地春小麦不同生育阶段平均温度和产量的灰色关联矩阵。对不同生育阶段平均温度与旱地春小麦产量数据做灰色关联度分析。研究结果表明，灌浆—成熟阶段平均气温是影响陇中黄土丘陵沟壑区旱地春小麦产量的主要因素，其次依次为开花—灌浆、播种—出苗、孕穗—开花、拔节—孕穗、出苗—分蘖、休闲和分蘖—拔节阶段（表5）。

pagenumber_ebook=60,pagenumber_book=910
图2 1970—2018年研究区不同生育阶段平均温度年际变化
Fig.2 Interannual variations of average temperatures at different growth stages from 1970-2018

表5 旱地春小麦产量与不同生育阶段内平均温度的灰色关联系数
Table 5 Grey correlation coefficients between the average temperatures at different growth stages and spring wheat yield

pagenumber_ebook=60,pagenumber_book=910
2.4 不同生育阶段温度变化对春小麦产量的影响
将1970—2018年旱地春小麦播种—出苗、孕穗—开花、开花—灌浆和灌浆—成熟阶段的逐日温度在±2.5℃的范围内降低或者增加0.5℃，并基于该气象数据利用APSIM模型模拟调整后气象条件下的小麦产量（表6）。可以看出，不同生育阶段温度变化对旱地春小麦产量的影响各不相同，播种—出苗阶段增温对春小麦产量呈正效应，温度每增加0.5℃，小麦产量平均提高0.45%。孕穗—开花、开花—灌浆和灌浆—成熟阶段增温均与春小麦产量呈负效应，其中，孕穗—开花阶段温度每增加0.5℃，小麦产量平均降低0.34%；开花—灌浆阶段温度每增加0.5℃，小麦产量平均降低0.65%；灌浆—成熟阶段温度每增加0.5℃，小麦产量平均降低1.09%。研究结果表明，孕穗—开花、开花—灌浆和灌浆—成熟阶段增温会使研究区旱地春小麦减产，其中灌浆—成熟阶段增温对小麦产量的影响最大。同时由于播种时研究区温度较低，增温可以提高春小麦的出苗率，从而提高小麦产量。

表6 不同生育阶段温度变化对春小麦产量的影响
Table 6 Influences of temperature changes during different growing stages on spring wheat yield

pagenumber_ebook=61,pagenumber_book=911
同列不同小写字母表示处理间在0.05水平差异显著 Different small letters in the same column indicate significant differences among treatments at 0.05 level

2.5 春小麦产量和不同生育阶段温度变化的回归分析
以播种—出苗阶段温度变化量X1、孕穗—开花阶段温度变化量X2、开花—灌浆阶段温度变化量X3和灌浆—成熟阶段温度变化量X4为自变量，以APSIM模拟的春小麦产量Y为因变量，建立二次多项式逐步回归方程（公式19）。该二次多项式逐步回归方程的相关系数R为0.98，F=28.97，大于F0.05，说明二次多项式逐步回归方程达到了显著水平；决定系数R2为0.96，说明二次回归方程可以较好拟合春小麦产量变化过程与不同生育阶段温度变化的对应关系。

pagenumber_ebook=61,pagenumber_book=911
由于二次回归方程自变量均为温度，在量纲相同的情况下，偏回归系数可以反映某一生育阶段温度变化对小麦产量产生的效应。其中，播种—出苗阶段温度变化量 X1的偏回归系数为25.67，为正效应；孕穗—开花阶段温度变化量 X2的偏回归系数为-17.84，为负效应；开花—灌浆阶段温度变化量 X3的偏回归系数为-30.15，为负效应；灌浆—成熟阶段温度变化量 X4的偏回归系数为-45.67，为负效应。研究结果表明，研究区不同生育阶段温度变化对小麦产量的影响程度大小依次为X4＞X3＞X1＞X2，且随温度升高，负效应大于正效应。

为了进一步分析单生育阶段温度变化对的小麦产量的影响，对二次多项式逐步回归方程进行降维，分析其他生育阶段温度变化固定在零水平时的产量效应（公式20—23）。其中，一次项系数为直接贡献率，正负号表示贡献方向；二次项系数为间接贡献率，该系数为正时，表示叠加正效应；该系数为负时，表示多项式到一定程度（临界点）会出现报酬递减。从回归结果（公式20—23，图3）可以看出，播种—出苗阶段和孕穗—开花阶段温度变化与小麦模拟产量效应为直线，而开花—灌浆和灌浆—成熟阶段温度变化与小麦产量效应为曲线。结果表明，播种—出苗阶段温度升高对春小麦产量的影响呈正效应，而孕穗—开花、开花—灌浆和灌浆—成熟阶段温度升高均对春小麦产量的影响呈负效应。这主要是由于播种前期研究区温度较低，增温可以提高旱地春小麦的出苗率，从而提高粮食产量。而孕穗—开花、开花—灌浆和灌浆—成熟阶段升温加快了旱地小麦的生长发育进程，导致小麦生育期缩短，不利于粮食产量的提高。

pagenumber_ebook=62,pagenumber_book=912
pagenumber_ebook=62,pagenumber_book=912
边际产量能够分析不同生育阶段温度变化对小麦产量的影响速率，春小麦产量的单因素边际效应可以通过二次多项式逐步回归方程降维，并对其求一阶偏导得到：播种—出苗阶段，dy/dx=25.67；孕穗—开花阶段，dy/dx=-17.84；开花—灌浆阶段，dy/dx=-30.15-13.64X3；灌浆—成熟阶段，dy/dx=-45.76-15.5X4。

pagenumber_ebook=62,pagenumber_book=912
图3 不同生育阶段温度变化对旱地春小麦产量影响的单因素效应
Fig.3 Single-factor effect of temperature change of different growth stages on spring wheat yield in dryland

pagenumber_ebook=62,pagenumber_book=912
图4 不同生育阶段温度变化对旱地春小麦产量影响的单因素边际效应
Fig.4 The single-factor marginal effect of temperature change of different growth stages on spring wheat yield

单因素边际效应分析反映了各种因素边际产量效应随着温度波动的变化情况，播种—出苗和孕穗—开花阶段温度的边际效应为常数，开花—灌浆和灌浆—成熟阶段温度边际效应随温度增加而呈递减趋势（图4）。结果表明，播种—出苗和孕穗—开花阶段温度变化对春小麦产量的影响不显著；而开花—灌浆和灌浆—成熟阶段温度变化对春小麦产量的影响表现为随温度升高，春小麦边际产量呈递减趋势。与播种—出苗、孕穗—开花阶段温度变化相比，开花—灌浆和灌浆—成熟阶段温度变化对春小麦产量的影响程度更大，对春小麦产量形成起主导作用。因此，温度升高引起春小麦减产，其主要原因是由于开花—灌浆、灌浆—成熟阶段温度升高而引起的。

不同生育阶段合理的温度可以有效提高旱地春小麦的生长发育和粮食产量，通过二次多项式逐步回归方程的交互项系数可以进一步分析不同生育期阶段温度变化的协同效应。其中，多元回归方程（公式19）X1X2项的交互系数为7.13，表明播种—出苗和孕穗—开花阶段增温之间呈正效应。而X1X3项的交互系数为-16.89，表明播种—出苗和开花—灌浆阶段增温之间呈负效应。通过二次多项式逐步回归方程的各项系数表明，影响旱地春小麦产量的因素不仅取决于播种—出苗、孕穗—开花、开花—灌浆和灌浆—成熟阶段某个单独生育阶段的温度变化，还取决于他们之间的交互效应，只有5者同时满足条件，才能使旱地春小麦的产量达到最大。

3 讨论
温度是植物生长发育的重要气象因子，温度升高会导致作物产量及其构成发生变化[24-26]。温度变化对小麦各生长发育阶段的影响不完全一致，各阶段温度变化对产量及产量构成要素的影响也存在着差异[27]。侯学会等[28]研究认为，黄土高原冬小麦生产主要受生育期内月平均温度高低的影响。肖国举等[29]研究认为，孕穗期增温将导致宁夏引黄灌区小麦光合速率明显下降，从而影响干物质累积和千粒重，最终影响当地小麦生产。张秀云等[30]研究认为，黄土高原半干旱区5月下旬至6月上旬是气温影响春小麦生长发育的关键时期，气候变暖及气温偏高会导致春小麦产量下降。本研究结果表明，不同生育阶段温度对旱地春小麦产量的影响大小依次为：灌浆—成熟阶段＞开花—灌浆阶段＞播种—出苗阶段＞孕穗—开花阶段＞拔节—孕穗阶段＞出苗—分蘖阶段＞休闲阶段＞分蘖—拔节阶段，且孕穗—开花、开花—灌浆和灌浆—成熟阶段增温对小麦产量呈负效应，而播种—出苗阶段增温对小麦产量呈正效应，这与赵鸿等[31]研究认为开花—乳熟阶段的温度增加和产量的相关性最大，出苗—拔节、开花—成熟阶段的温度增加以及拔节—孕穗阶段的温度降低，是造成小麦产量增加的直接原因的研究结果不一致，这可能主要与2个研究区的气候条件不同有关。

有研究表明，温度每升高1℃，小麦产量将会减少4%—7%。相反，也有研究表明温度每升高1.0—1.5℃，小麦产量会提高15%—20%[32-33]。任新庄等[14]研究表明温度每升高 0.5℃，陇中旱地春小麦产量最大降幅为4.92%，最小降幅为1.94%，平均降幅为3.24%。张凯等[16]研究表明，在灌浆后期到乳熟期间，温度每增加1—2℃，春小麦的穗粒数和穗粒质量显著减少，产量较对照降低25.4%—45.5%。普宗朝等[34]研究表明，苗期平均气温增加对小麦产量的影响呈正效应，影响系数为0.5—32.0 kg·hm-2·℃-1，而开花期—成熟期温度增加对小麦产量影响呈负效应，影响系数为-22.9—-8.0 kg·hm-2·℃-1。而本研究表明，播种—出苗阶段温度每提高0.5℃，春小麦产量平均增加0.45%；孕穗—开花阶段温度每提高0.5℃，春小麦产量平均减少0.34%。开花—灌浆阶段温度每提高0.5℃，春小麦产量平均减少0.65%；灌浆—成熟阶段温度每增加0.5℃，春小麦产量平均减少1.09%。研究结果与前人相似，只是在产量随温度变化的幅度上有一定的差异，差异主要是由研究区的特殊性和试验设计的不同造成的。

本研究利用APSIM模型模拟陇中黄土丘陵沟壑区旱地春小麦4个主要生育阶段（播种—出苗阶段、孕穗—开花阶段、开花—灌浆阶段、灌浆—成熟阶段）温度变化对旱地春小麦产量的影响，但模拟过程并未考虑其他生育期温度变化、小麦品种、降水、CO2浓度变化以及管理措施等因素对研究结果的影响，造成研究结果存在不确定因素。此外，气候变化会导致病虫害和极端天气的增加，同时病虫害和极端天气也是影响小麦产量的重要因素，在后续研究中可以考虑增加病虫害、极端天气等边界条件，从而进一步提高模型的模拟精度。

4 结论
本文以陇中黄土丘陵沟壑区为研究区域，以旱地春小麦为研究对象，结合APSIM模型和大田试验，通过灰色关联度分析、二次多项回归分析、单因素边际效应分析等研究方法定量分析了旱地春小麦产量对不同生育阶段温度变化的响应机制。研究表明，APSIM模型对模拟陇中黄土丘陵沟壑区不同播种期旱地春小麦的产量和主要生育期具有较好的模拟效果。研究区旱地春小麦全生育期平均温度总体呈上升趋势，但不同生育阶段的温度增幅各不相同。不同生育阶段温度变化对研究区旱地春小麦产量影响大小各不相同，其中，灌浆—成熟、开花—灌浆、播种—出苗、孕穗—开花阶段温度变化是影响陇中黄土丘陵沟壑区旱地春小麦产量的关键生育阶段。

References

[1] ZHAO H Y, GONG L J, QU H H, ZHU H X, LI X F, ZHAO F.The climate change variations in the northern greater khingan mountains during the past centuries.Journal of Geographical Sciences, 2016,26(5): 585-602.

[2] 秦大河, STOCKEER T.IPCC第五次评估报告第一工作组报告的亮点结论.气候变化研究进展, 2014, 10(1):1-6.QIN D H, STOCKER T.Highlight of the IPCC working group I Fifth assessment report.Climate Change Research, 2014, 10(1):1-6.(in Chinese)

[3] 张凯, 王润元, 冯起, 王鹤龄, 赵鸿, 赵福年, 阳伏林, 雷俊.模拟增温和降水变化对半干旱区春小麦生长及产量的影响.农业工程学报, 2015, 31(1): 161-170.ZHANG K, WANG R Y, FENG Q, WANG H L, ZHAO H, ZHAO F N, YANG F L, LEI J.Effect of simulated warming and precipitation change on growth characteristics and grain yield of spring wheat in semi-arid area.Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(1): 161-170.(in Chinese)

[4] KHEIRI M, SOUFIZADEH S, GHAFFARI A, AGHAALIKHANI M,ESKANDARI A.Association between temperature and precipitation with dryland wheat yield in northwest of Iran.Climatic Change, 2017,141(4): 703-717.

[5] 裴雪霞, 党建友, 张定一, 武雪萍, 赵娟.近54年来晋南气候变化及其对旱地小麦产量的影响.麦类作物学报, 2016, 36(11):1502-1509.PEI X X, DANG J Y, ZHANG D Y, WU X P, ZHAO J.Climate change during nearly 54 years in south of shanxi and its effect on wheat yield in dry land.Journal of Triticeae Crops, 2016, 36(11):1502-1509.(in Chinese)

[6] 孙正国.模拟增温对小麦生长及土壤的调控效应.核农学报, 2016,30(10): 2041-2048.SUN Z G.Control response of temperature enhancement on growth and soil health of wheat.Journal of Nuclear Agricultural Sciences,2016, 30(10): 2041-2048.(in Chinese)

[7] 房世波, 谭凯炎, 任三学.夜间增温对冬小麦生长和产量影响的实验研究.中国农业科学, 2010, 43(15): 3251-3258.FANG S B, TAN K Y, REN S X.Winter wheat yields decline with spring higher night temperature by controlled experiments.Scientia Agricultura Sinica, 2010, 43(15): 3251-3258.(in Chinese)

[8] 张凯, 王润元, 王鹤龄, 赵鸿, 赵福年, 阳伏林, 雷俊.田间增温对半干旱区春小麦生长发育和产量的影响.应用生态学报, 2015,26(9): 2681-2688.ZHANG K, WANG R Y, WANG H L, ZHAO H, ZHAO F N, YANG F L, LEI J.Effects of increasing field temperature on growth,development and yield of spring wheat in semi-arid area. Chinese Journal of Applied Ecology, 2015, 26(9): 2681-2688.(in Chinese)

[9] 田云录, 陈金, 邓艾兴, 郑建初, 张卫建.开放式增温下非对称性增温对冬小麦生长特征及产量构成的影响.应用生态学报, 2011,22(3): 681-686.(in Chinese)TIAN Y L, CHEN J, DENG A X, ZHENG J C, ZHANG W J.Effects of asymmetric warming on the growth characteristics and yield components of winter wheat under free air temperature increased.Chinese Journal of Applied Ecology, 2011, 22(3): 681-686.(in Chinese)

[10] DALLA M A, EITZINGER J, KERSEBAUM K C, TODOROVIC M,ALTOBELLI F.Assessment and monitoring of crop water use and productivity in response to climate change.The Journal of Agricultural Science, 2018, 156: 575-576.

[11] 尚艳, 赵鸿, 柴守玺.气候变化与品种更新对黄土高原半干旱雨养农业区冬小麦的影响.干旱地区农业研究, 2017, 35(5):66-72.SHANG Y, ZHAO H, CHAI S X.Effects of climate change and cultivars change on winter wheat in Semi-arid region of Loess Plateau in northeast China.Agricultural Research in the Arid Areas, 2017,35(5): 66-72.(in Chinese)

[12] 高美玲, 张旭博, 孙志刚, 孙楠, 李仕冀, 高永华, 张崇玉.中国不同气候区小麦产量及发育期持续时间对田间增温的响应.中国农业科学, 2018, 51(2): 386-400.GAO M L, ZHANG X B, SUN Z G, SUN N, LI S J, GAO Y H,ZHANG C Y.Wheat yield and growing period in response to field warming in different climatic zones in China.Scientia Agricultura Sinica, 2018, 51(2): 386-400.(in Chinese)

[13] 董莉霞, 李广, 刘强, 燕振刚, 罗珠珠.旱地春小麦产量对逐日最低温度和最高温度变化响应的模拟与分析.中国生态农业学报,2013 , 21(8): 1016-1022.DONG L X, LI G, LIU Q, YAN Z G, LUO Z Z.Simulation analysis of spring wheat grain yield response to mean daily minimum and maximum temperature in drylands.Chinese Journal of Eco-Agriculture,2013, 21(8): 1016-1022.(in Chinese)

[14] 任新庄, 闫丽娟, 李广, 聂志刚, 王钧, 罗永忠.陇中旱地春小麦产量对降水与温度变化的响应模拟.干旱地区农业研究, 2018,36(3): 125-129.REN X Z, YAN L J, LI G, NIE Z G, WANG J, LUO Y Z.Simulation of the effects of precipitation and temperature change on spring wheat yield in dryland of central Gansu.Agricultural Research in the Arid Areas, 2018, 36(3): 125-129.(in Chinese)

[15] 张凯, 王润元, 王鹤龄, 赵鸿, 齐月, 赵福年, 雷俊.模拟增温对半干旱雨养区春小麦物质生产与分配的影响.农业工程学报, 2016,32(16): 223-232.ZHANG K, WANG R Y, WANG H L, ZHAO H, QI Y, ZHAO F N,LEI J.Effect of simulated warming on dry matter production and distribution of rainfed spring wheat in semi-arid area.Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(16):223-232.(in Chinese)

[16] 张凯, 王润元, 王鹤龄, 赵鸿, 赵福年, 阳伏林, 陈斐, 齐月, 雷俊.温度升高和降水减少对半干旱区春小麦生长发育及产量的协同影响.中国生态农业学报, 2019, 27(3): 413-421.ZHANG K, WANG R Y, WANG H L, ZHAO H, ZHAO F N, YANG F L, CHEN F, QI Y, LEI J.Influence of climate warming and rainfall reduction on semi-arid wheat production.Chinese Journal of Eco-Agriculture, 2019, 27(3): 413-421.(in Chinese)

[17] PRADHAN G P, PRASAD P V V, FRITZ A K, KIRKHAM M B.Effects of drought and high temperature stress on synthetic hexaploidy wheat.Functional Plant Biology, 2012, 39(3): 190-198.

[18] WCISLO R, MILLER S S, DZWINEL W.Particle automata model in simulation of Fusarium graminearum pathogen expansion.Journal of Theoretical Biology, 2016, 389: 110-122.

[19] JIN Z N, RISHI P, JOHN S, ZHUANG Q L.Crop model and satellite imagery-based recommendation tool for variable rate N fertilizer application for the US corn system.Precision Agriculture, 2017, 18:779-800.

[20] FAN D L, YANG F Y, PAN Z H, SU X Y.Development of an improved model to evaluate vulnerability in spring wheat under climate change in Inner Mongolia.Sustainability, 2018, 10(12):4581-4596.

[21] PALMER J, THORBURN P J, BIGGS J S, DOMINATI E J,PROBERT M E, MEIER E A, HUTH N I, DODD M, SNOW V,LARSEN J R, PARTON W J.Nitrogen cycling from increased soil organic carbon contributes both positively and negatively to ecosystem services in wheat agro-ecosystems.Frontiers In Plant Science, 2017, 8: 1-12.

[22] 李广, 黄高宝, WILLIAM B, 陈文.APSIM 模型在黄土丘陵沟壑区不同耕作措施中的适用性.生态学报, 2009, 29: 2655-2663.LI G, HUANG G B, WILLIAM B, CHEN W.Adaptation research of APSIM model under different tillage systems in the Loess hill-gullied region.Acta Ecology Sinica, 2009, 29: 2655-2663.(in Chinese)

[23] ZHENG B Y, CHENU K, DOHERTY A, CHAPMAN S.The APSIMWheat Module (7.5 R3008) Documentation.Toowoomba:APSRU, 2014: 20-25.

[24] 杨扬, 罗贤, 李荣平, 王莹, 温日红, 杨桂娟, 焦敏, 张琪, 张晓月.气象要素对植物物候影响及其驱动机制研究进展.气象与环境学报, 2016, 32(5): 154-159.YANG Y, LUO X, LI R P, WANG Y, WEN R H, YANG G J, JIAO M, ZHANG Q, ZHANG X Y.A review of the effect of meteorological factors on plant phenology and its driving mechanisms.Journal of Meteorology and Environment, 2016, 32(5):154-159.(in Chinese)

[25] 霍治国, 尚莹, 邬定荣, 吴立, 范雨娴, 王培娟, 杨建莹, 王纯枝.中国小麦干热风灾害研究进展.应用气象学报, 2019, 30(2):129-138.HUO Z G, SHANG Y, WU D R, WU L, FAN Y X, WANG P J, YANG J Y, WANG C Z.Review on disaster of hot dry wind for wheat in China. Journal of Applied Meteorological Science, 2019, 30(2):129-138.(in Chinese)

[26] 李阔, 熊伟, 潘婕, 林而达, 李迎春, 韩雪.未来升温1.5oC与2.0oC背景下中国玉米产量变化趋势评估.中国农业气象, 2018,39(12): 765-777.LI K, XIONG W, PAN J, LIN E D, LI Y X, HAN X.Trend evaluation on changes of maize yield in China under global warming by 1.5℃and 2.0℃.Chinese Journal of Agrometeorology, 2018, 39(12):765-777.(in Chinese)

[27] 孟自力, 陈昆, 闫向泉, 朱倩, 倪雪峰, 朱伟.气象因子变化及其对小麦生产的影响.安徽农业科学, 2018, 46(7): 27-29.MENG Z L, CHEN K, YAN X Q, ZHU Q, NI X F, ZHU W.Research status of change of meteorological factors and its influence on wheat production.Journal of Anhui Agricultural Sciences, 2018, 46(7):27-29.(in Chinese)

[28] 侯学会, 隋学艳, 姚慧敏, 梁守真, 王猛.中国北方麦区冬小麦物候期对气候变化的响应.麦类作物学报, 2019, 39(2): 202-209.HOU X H, SUI X Y, YAO H M, LIANG S Z, WANG M.Response of winter wheat phenology to climate change in northern China.Journal of Triticeae Crops, 2019, 39(2): 202-209.(in Chinese)

[29] 肖国举, 张强, 张峰举, 罗成科, 王润元.增温对宁夏引黄灌区春小麦生产的影响.生态学报, 2011, 31(21): 6588-6592.XIAO G J, ZHANG Q, ZHANG F J, LUO C K, WANG R Y.The impact of rising temperature on spring wheat production in the Yellow River irrigation region of Ningxia.Acta Ecologica Sinica, 2011,31(21): 6588-6592.(in Chinese)

[30] 张秀云, 王鹤龄, 雷俊.气候暖干化对半干旱区春小麦产量形成的影响.生态环境学报, 2015, 24(4): 569-574.ZHANG X Y, WANG H L, LEI J.Impacts of climate warming and drying on spring wheat yield in a semi-arid region.Ecology and Environmental Sciences, 2015, 24(4): 569-574.(in Chinese)

[31] 赵鸿, 何春雨, 李凤民, 王润元, 杨启国, 邓振镛, 王鹤龄.气候变暖对高寒阴湿地区春小麦生长发育和产量的影响.生态学杂志,2008, 27(12): 2111-2117.ZHAO H, HE Y C, LI F M, WANG R Y, YANG Q G, DENG Z Y,WANG H L.Effects of climate warming on spring wheat growth and yield in high-altitude cold and dankness region.Chinese Journal of Ecology, 2008, 27(12): 2111-2117.(in Chinese)

[32] HATFIELD J L, BOOTE K J, KIMBALL B A, ZISKA L H,IZAURRALDE R C, ORT D, THOMSON A M, WOLFE D.Climate impacts on agriculture: Implications for crop production.Agronomy Journal, 2011, 103(2): 351-370.

[33] TACK J, BARKLEY A, NALLEY L L.Effect of warming temperatures on US wheat yields.Proceedings of the National Academy of Sciences of the USA, 2015, 112(22): 6931-6936.

[34] 普宗朝, 张山清, 徐文修, 李景林, 王荣晓.气候变化对伊犁河谷冬小麦产量的影响.中国农学通报, 2014, 30(15): 173-182.PU Z C, ZHANG S Q, XU W X, LI J L, WANG R X.Impact of climate change to winter wheat yield in Yili river valley of Xinjiang.Chinese Agricultural Science Bulletin, 2014, 30(15): 173-182.(in Chinese)

Simulation of Spring Wheat Yield Response to Temperature Changes of Different Growth Stages in Drylands

WANG Jun1,2, LI Guang1, YAN LiJuan3, LIU Qiang2, NIE ZhiGang2
(1College of Forestry, Gansu Agricultural University, Lanzhou 730070; 2College of Information Science and Technology, Gansu Agricultural University, Lanzhou 730070; 3Agronomy College, Gansu Agriculture University, Lanzhou 730070)

Abstract: 【Objective】This research quantitatively analyzed the effects of temperature changes in different growth stages on the yield of spring wheat in the Loess Plateau Gully Region of Central Gansu, so as to provide theoretical guidance and decisionmaking basis for making reasonable cropping pattern in the Loess Plateau Gully Region of Central Gansu 【Method】Specifically, it investigated the relationship between temperature changes in different growth stages and the yield of spring wheat in dryland based on grey relational analysis, so as to identify the key growth stages of spring wheat.The yield of spring wheat was simulated using the APSIM (Agricultural Production System Simulator) model under different temperature conditions in different growth stages.The mechanism of yield response to temperature changes in different growth stages was studied through quadratic polynomial regression analysis and single factor marginal effect analysis.【Result】(1) The results showed that the APSIM model had good applicability in simulating yield and growth stages of spring wheat in the Loess Plateau Gully Region of Central Gansu.The values of the average root mean square error (RMSE), normalized root mean square error (NRMSE), and model effectiveness index (ME) of the yield of the simulation model were 39.95 kg·hm-1, 1.55% and 0.73, respectively.The values of the RMSE, NRMSE, and ME of the simulated model phenology dates were 2.78 d, 1.87% and 0.83, respectively.(2)The average temperature increased steadily during all growth stages of spring wheat in the Loess Plateau Gully Region of Central Gansu, and the average temperatures of the different growth stages showed an increasing trend.The increase of average annual temperature were 0.44℃·(10a)-1, 0.34℃·(10a)-1, 0.17℃·(10a)-1, 0.41℃·(10a)-1, 0.49℃·(10a)-1, 0.52℃·(10a)-1 and 0.35℃·(10a) -1 in sowing-emergence, emergence-tillering, tillering-jointing, jointing-booting, booting-flowering, flowering-grain filling, grain filling-maturity, respectively.(3) The importance of temperature changes in different growth stages on the yield of spring wheat were ranked as follows: grain filling-maturity, flowering-grain filling, sowing-emergence, booting- flowering, jointing-booting, emergencetillering, fallow, and the tillering-jointing stage.(4) Under constant temperatures at the other growth stages, a temperature increase of 0.5℃during the sowing-emergence stage could increase the wheat yield by 0.45%; A temperature increase of 0.5℃ during the booting-flowering stage could reduce the wheat yield by 0.34%; a temperature increase of 0.5℃ during the flowering-filling stage could reduce the wheat yield by 0.65%; A temperature increase of 0.5℃ during the filling-maturity stage could reduce the wheat yield by 1.09%.When temperature changes during the sowing-emergence and booting-flowering stages were the same, the effects of the temperature changes on the yield of spring wheat were similar.【Conclusion】APSIM model has achiedved good simulation results on the yield and growth stages of spring wheat of different sowing dates in the Loess Plateau Gully Region of Central Gansu.The trend of temperature increase obviously existed in the study area during the whole growth stages of spring wheat, but temperature increment was different in different growth stages.The effects of temperature changes in different growth stages on the yield of spring wheat were different.Warming was beneficial to increasing yield during sowing to emergence, but warming led to yield reduction during other growth stages.

Key words: spring wheat yield; different growth stages; temperature change; APSIM; grey relation analysis

doi: 10.3864/j.issn.0578-1752.2020.05.004

开放科学（资源服务）标识码（OSID）：

pagenumber_ebook=54,pagenumber_book=904
收稿日期：2019-08-02；

接受日期：2019-09-12

基金项目：甘肃农业大学学科建设基金（GAU-XKJS2018258）、国家自然科学基金（31660348）、甘肃省重点研发计划（18YF1NA070）、甘肃省高等学校协同创新团队项目（2018C-16）、甘肃省财政专项（GSCZZ-20160909）、甘肃农业大学青年导师基金项目（GAU-QNDS-201702）、甘肃农业大学信息科学技术学院发展基金项目（GAU-XKFZJJ-2020-12）、甘肃省自然科学基金项目（18JR3RA169）

联系方式：王钧，E-mail：julianwong82@163.com。

通信作者李广，E-mail：lig@gsau.edu.cn

（责任编辑 杨鑫浩）