冬小麦属喜凉作物，灌浆阶段天气凉爽有利于提高千粒重^[1]，但在我国冬小麦主产区小麦灌浆成熟期常会遭遇高温天气，高温危害是小麦取得高产的主要障碍之一。随着全球气候变暖，一些地区小麦灌浆期的高温胁迫发生频率增加^[2]，极端温度事件将给农作物产量造成更大损失^[3]，我国未来极端高温事件也将明显增加^[4]，年平均高温日数呈上升趋势^{[5, 6]}，冬小麦生育后期高温危害势必对高产稳产带来更大影响，威胁我国冬小麦产量安全。关于高温胁迫对小麦的光合、生理、产量及籽粒品质的影响国内外已开展了较多的试验^{[7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]}，一些试验研究表明灌浆期高温胁迫致使冬小麦产量降低^{[15, 16, 17]}，在不同灌浆阶段高温胁迫的作用不同^{[18, 19]}。然而，针对灌浆期高温胁迫程度与冬小麦籽粒产量的定量关系研究还较少，究竟灌浆期短期高温对冬小麦产量的影响遵从何种规律，气候异常时出现极端高温是否将导致冬小麦产量陡降，以及高温胁迫对冬小麦灌浆与产量的作用机制都有待深入研究。华北地区是我国冬小麦主产区之一，本研究试图通过田间模拟高温试验探索冬小麦灌浆中后期短期高温对华北地区冬小麦产量的影响规律及其作用机制，以期为评估未来极端高温事件对该地区冬小麦产量安全威胁提供科学依据。在冬小麦灌浆中后期，穗数和穗粒数已经基本确定，粒重大小决定了最终籽粒产量的高低，因此，本研究通过粒重的变化来反映产量对高温胁迫的响应。

本研究的试验于2012—2013年在中国气象科学研究院河北固城生态与农业气象试验站进行，两年试验的冬小麦品种均为济麦22，属半冬性品种。试验在大田中进行，各小区播种密度和田间管理措施完全一致，增温处理在冬小麦灌浆中后期(开花11 d后至成熟^[19])进行，处理前冬小麦长势基本相同，在自然状态下采用红外辐射器增温方法对小麦冠层进行辐射增温。设两组增温处理，一组为在冬小麦灌浆中后期选择6个晴热天进行增温处理(TI 6 d)，另一组为在冬小麦灌浆中后期选择9个晴热天进行增温处理(TI 9 d)，由于天气条件的波动，6 d处理和9 d处理不一定连续。每组设3个增温幅度(T1,T2,T3)，通过采用不同功率红外辐射器及其安装高度产生不同增温幅度，辐射器由3根红外灯管并排组成，3个增温幅度的红外辐射器功率和安装高度分别为3×500 W,距冠层顶部0.6 m；3×600 W，距冠层顶部0.4 m；3×700 W,距冠层顶部0.4 m。另外设置1组对照(CK)。每个幅度的增温处理和对照均设3个小区重复，每个小区的面积为1.5 m×2.0 m，红外辐射器水平悬挂在小区中央冬小麦冠层上方，辐射器照射区域0.8 m×1.2 m，辐射器灯管的反射罩宽度2.5 cm，其对阳光的遮挡作用可以忽略。增温处理只在晴热天气的午间进行，时间为10:00—15:00，在每个处理和对照安装1个带辐射罩的温度传感器，传感器安装在照射区域中央的冬小麦冠层高度处，采集器每10 min采集1次数据。

根据距试验地40 km的涿州农气站两个生长季各月温度距平资料(表1)，2011—2012冬小麦生长季温度较常年偏高，而2012—2013冬小麦生长季温度偏低，2012—2013季从冬小麦苗期开始气温多低于上一季，特别是起身期以后气温持续低于上年。两个生长季未出现明显的气象和病虫灾害情况。

小麦开花后，在各处理小区的中央区域选定高度和大小相近的40穗，做出标记，开花5 d后在对照小区开始取样观测籽粒灌浆速率，每隔4 d取样1次，增温处理开始后，各小区同时观测。每次在每小区取5穗，剥出籽粒在80 ℃恒温下烘干后称重，计算平均粒重，小麦成熟后考种得到最终粒重。同期的百叶箱气温、风速和空气湿度值来自本站的自动气象站观测资料。

冠层气温小时值为小时内6次读数的平均，午间冠层气温平均值为11:00—15:00冠层气温小时值的平均。高温有效度时值为冬小麦灌浆中后期(2012年5月20日—6月8日，2013年5月26日—6月18日)冠层气温小时值中高于30 ℃部分的和。为了消除年份间除灌浆中后期气温以外其它条件对粒重的影响，对两年的粒重数据进行了标准化处理，由分年取得的粒重与午间冠层平均气温的回归关系式或粒重与高温有效度时值的回归关系式计算多年平均午间冠层气温或高温有效度时值为零所对应的粒重值，然后各处理的粒重数据与它的比值即为标准化后的粒重相对值。采用Excel 软件进行数据分析处理和作图，用 MATLAB 软件作粒重增长曲线模拟和分析，用SPSS统计软件进行显著性分析。

2012年开花14 d后进行9 d增温处理(T3)与对照的千粒重变化和2013年开花16 d后进行9 d增温处理(T2)与对照的千粒重变化见图1。图1直观表明，增温后，千粒重增速较对照明显减小，同时，增温处理的籽粒增长过程比对照提前终止，增温使灌浆持续时间比对照缩短约3—5 d，从观察到的冠层表征显示，增温处理没有灼伤叶片，但使叶片和麦穗提前变黄。应用Logistic模型模拟图1中的4条灌浆进程曲线后求得，增温处理使2012年最大灌浆速率从2.41 mg d^-1 粒^-1下降到2.00 mg d^-1 粒^-1；使2013年最大灌浆速率从2.09 mg d^-1 粒^-1下降到1.70 mg d^-1 粒^-1。由此可见，灌浆期高温胁迫造成粒重下降的机制是降低灌浆速率和缩短灌浆持续时间的双重影响。

图 1 高温对灌浆进程及灌浆速率的影响 Fig.1 Impacts of high temperature on process and rate of grain filling

图选项

表2给出了2012年和2013年冬小麦灌浆中后期短期(6 d和9 d)增温处理的冠层气温及其粒重数据。首先可见，由于两个生长季气候条件不同，两年间对照的温度和粒重均相差很大；同时从表中可看出，各组试验中，冬小麦最终粒重都随增温处理期间午间平均冠层气温的升高而减小。两年不同的背景温度条件使试验取得了更宽温度范围内粒重与冠层气温的关系数据。

经回归分析，分别得出各组试验中粒重与冠层气温的相关关系如下：

点击浏览公式

点击浏览公式

点击浏览公式

点击浏览公式

式中,T为增温处理期间的午间(11:00—15:00)平均冠层气温(℃)；G为粒重(mg/粒)；* 表示通过P≤0.10的显著性检验，** 表示通过P≤0.05显著性检验。

分析表明，在不同年份、不同高温处理持续时间下，粒重与高温处理期间的平均午间冠层气温间都呈显著的负线性相关关系。

前面的分析中所用温度资料均为冬小麦群体冠层气温，为了与常规气象资料联系起来，本文利用冬小麦灌浆中后期每天11:00—15:00对照小区冠层小时平均气温和同时的百叶箱气温资料，分别分析得到冠层气温与百叶箱气温(式5)，百叶箱午间(11:00—15:00)平均气温与日平均气温(式6)和日最高气温(式7)间的关系如下：

点击浏览公式

点击浏览公式

点击浏览公式

式中,T_a为百叶箱气温； T_ca为冠层气温； T_md 为百叶箱午间(11:00—15:00)平均气温； T_av为百叶箱日平均气温； T_M为百叶箱日最高气温； ***表示通过P≤0.01的显著性检验。

结果表明，在冬小麦灌浆中后期，气温较低时，冠层气温与百叶箱气温接近，气温较高时，冠层气温则高于百叶箱气温，天气越炎热，两者相差越大。其原因可能在于，在冬小麦灌浆中后期，土壤湿度较低，午后的冠层温度一般高于空气温度^[20]，到后期，随着根系吸水能力下降，蒸腾作用减小，午间光照越强，冠层温度越高，冠气温差越大，冠层气温也就越高于百叶箱气温。根据观测资料，在灌浆中期以前，冠层气温一般接近或低于百叶箱气温，这可能与前期土壤湿度较高及冠层蒸腾能力较强有关。因此，式5只适应于冬小麦灌浆中后期午间时段。同时分析得到，午间平均气温与日平均气温和日最高气温间均存在非常密切的线性相关关系(式6、式7)。

根据涿州农气站历史资料分析得出，冬小麦乳熟至成熟期的多年(1991—2010)日平均气温为23.6 ℃，该温度值可以视作为当地近20年冬小麦灌浆中后期的平均温度条件，由上述关系式换算成灌浆中后期午间平均冠层气温为30 ℃。

因为两年试验冬小麦生长背景气象条件不同，造成不能将两年的粒重与冠层气温关系作统一分析。为了消除年份间除灌浆中后期气温以外其它条件对粒重的影响，分年对粒重数据进行了标准化处理，当地灌浆中后期午间冠层气温的多年平均值取30 ℃。图2给出了经标准化处理后两年试验中9d增温和6d增温的粒重相对值与平均冠层气温的关系，图2显示，经标准化后，两年试验数据遵从一致的相关关系，粒重相对值随冠层气温升高呈线性下降，回归关系均通过极显著检验(P≤0.01)，表明两年试验得到的高温减产规律相同。

图 2 标准化后粒重相对值与灌浆中后期平均冠层气温的关系, Fig.2 Relationship between normalized kernel weight and average canopy air temperature at mid to late grain filling stage

图选项

根据图2中两组回归关系可计算出，相对于多年平均冠层气温(30 ℃)，9d高温处理期间午间平均冠层气温每升高1 ℃，粒重降低3.4%，6d高温处理平均冠层气温每升高1 ℃，粒重降低约2.6%，这表明，高温强度相同时，持续时间较长导致更大的减产率。以当地灌浆中后期多年平均气温下的粒重作为正常值，依据通常气象灾害分级对应的减产百分率^[21]，以减产5%—10%为轻度受灾，10%—20%为中度受灾，而20%—30%为重度受灾，分别可以得出，在冬小麦灌浆中后期，如果出现9d平均午间冠层气温达到31.4、32.9 ℃和35.9 ℃，或6d平均午间冠层气温达到 32 ℃、33.8 ℃和37.6 ℃，则冬小麦将分别受到轻、中、重度高温减产危害，据此由式5—式7计算出相对应的日平均最高气温分别为31、32.1、34.4 ℃和31.5、32.8、35.6 ℃。

前面用平均冠层气温作为特征量分析了短期高温处理对粒重的影响，它们间很好的相关关系一定程度上与增温处理均选择了晴热天气有关，如果用灌浆中后期全阶段的平均气温来表示高温程度，则可能因期间出现多个低温日而使其不能真实反映期间高温天气的热害作用。为了体现灌浆中后期全阶段所有不利于冬小麦籽粒灌浆的高温对粒重的累积影响，本文用灌浆中后期高温积热(累计高温有效度时，℃·h)来作为全阶段的高温热害特征量，即冬小麦开花11d后至成熟期间所有小时气温值高于某一阈值的部分的累加。本文以试验地区冬小麦灌浆中后期多年平均气温对应的午间冠层气温值(30 ℃)作为阈值，计算了两年各处理的高温有效度时值。为了消除两年间冬小麦生育期其它条件对粒重的影响，分年对粒重数据进行标准化处理后得到粒重相对值。由两年各处理灌浆中后期累计高温有效度时和粒重相对值数据得到两者关系如图3，统计结果显示，粒重与反映灌浆中后期高温积热的有效度时值之间存在很好的线性相关关系，统计结果通过P≤0.05的显著性检验。分析结果表明，高温有效度时值较好地代表了灌浆中后期高温的强度及其持续时间对粒重的双重影响。按照2.2.3 中的方法，相应地可以计算出试验地区冬小麦灌浆中后期轻、中、重高温热害对应的全阶段累计高温有效度时值分别为144、215(℃·h)和358(℃·h)。

图 3 标准化后粒重相对值与灌浆中后期高温有效度时的关系 Fig.3 Relationship between normalized kernel weight and effective accumulated heat at mid to late grain filling stage

图选项

本研究的两年田间试验期间，分别只有1d的风速和湿度条件达到干热风的标准(限于篇幅，资料略)，因此，试验结果表现的主要是高温胁迫的影响效应。试验结果表明，灌浆中后期短期高温可以导致产量的显著降低，粒重随高温处理期间的平均冠层气温而线性下降，这与曲曼丽早期的田间观测结果规律一致^[15]，本试验中，灌浆中后期冠层气温差异是导致粒重差异的唯一因子，所以两者的相关关系更为密切。因试验处理增温幅度有限，其结果未能反映更高温度下的产量响应，但据此可以认为，在一定的温度变幅内，短期高温不至造成产量陡降。

作物受高温伤害的程度随高温强度越强或持续时间越长而越严重^[22]，本研究的分析结果显示，高温胁迫引起的冬小麦籽粒减产率也随高温胁迫强度增加、胁迫时间延长而加大，可见，冬小麦遭受高温胁迫危害后的减产幅度取决于高温强度及其持续时间两个因素。本文提出的高温有效度时有机结合了高温强度和胁迫时间的综合影响，全面反映了灌浆中后期高温危害程度，随着自动气象站的普及，由小时平均温度数据计算高温有效度时值切实可行。

试验表明，高温胁迫对冬小麦灌浆的影响表现在既抑制灌浆速率，同时缩短灌浆持续时间。在高温胁迫下，光合器官的功能受到伤害以及同化物向库器官的分配受到阻碍可能是高温胁迫致使灌浆速率下降的原因^{[9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]}，而高温胁迫使小麦灌浆期缩短则与后期持续高温导致植株体内SOD、POD酶的活性下降，加速根系与地上部组织衰老有关^{[24, 25]}。为了应对冬小麦生育后期高温危害，在生产中，除选用耐热性强的品种外，采取一定的栽培管理措施也可以减轻热害，根据研究，提高拔节期追施氮肥的比例可以延缓后期植株衰老,降低灌浆中后期高温胁迫对粒重、 蛋白质质量和产量的不利影响^[26]。

小麦灌浆期高温热害是一种偶发天气现象，其发生的时期、强度和持续天数因年因地而异。控制试验表明，高温胁迫对小麦产量的影响因高温处理的灌浆阶段、高温强度及其持续时间而不同^{[18, 19, 27]}，同时，冬小麦对高温的响应还与品种的耐热性密切相关^[10]，因此，要全面揭示冬小麦后期高温胁迫对产量的影响规律需要开展大量的试验研究。此外，本文中给出的试验地区高温热害的部分温度阈值来自于两年试验，也需要更多观测数据的验证。

致谢： 感谢中国气象科学研究院固城生态与农业气象试验站工作人员对本试验的大力协助，感谢梁宏博士在MATLAB软件应用上的帮助。