中國北部冬麥區小麥生育期對生育階段積溫變化的響應*

馬倩倩，賀 勇，張夢婷，張 聰，許吟隆

?

中國北部冬麥區小麥生育期對生育階段積溫變化的響應*

馬倩倩，賀 勇**，張夢婷，張 聰，許吟隆**

（中國農業科學院農業環境與可持續發展研究所，北京 100081）

利用1993-2013年中國北部冬麥區19個農業氣象觀測站的冬小麥生育期及氣象資料，研究了冬小麥各生育階段積溫和生育期的時空變異特征。通過皮爾遜相關性分析等方法，探究冬小麥各生育階段積溫變化對生育期的影響。結果表明：（1）播種-出苗、返青-拔節階段≥0℃積溫和越冬期負積溫的值呈東高西低的空間分布，拔節-抽穗、抽穗-乳熟、乳熟-成熟和播種-成熟階段≥0℃積溫為東南低西北高，而出苗-越冬開始階段≥0℃積溫則呈東南高西北低的分布；拔節-抽穗和乳熟-成熟階段≥0℃積溫均在21%的站點上顯著減少，返青-拔節、抽穗-乳熟和播種-成熟階段≥0℃積溫及越冬期負積溫分別在26%、37%、21%和42%的站點上顯著增加，而播種-出苗和出苗-越冬開始階段≥0℃積溫的變化較小；（2）播種和出苗期呈東部晚西部早的空間分布，抽穗、乳熟和成熟期則相反；越冬開始期呈東南晚西北早的分布，返青期則相反；拔節早的站點主要位于麥區東部。播種、出苗、返青、拔節、乳熟和成熟期分別在21%、16%、37%、26%、42%和21%站點顯著推遲且多位于麥區東部，而越冬開始期和抽穗期僅在5%站點變化顯著；（3）相關分析表明，各生育階段≥0℃積溫（或越冬期負積溫）與多個生育期的相關性顯著，生育階段積溫的變化可能直接或間接影響了冬小麥的生長發育。越冬期負積溫與返青、拔節、抽穗、乳熟和成熟期相關性最大，且與冬后多個生育期呈現一致的時空變異特征，其時空變異性可能是造成冬小麥冬后生育期在時空上存在差異的原因。

冬小麥；生育期；積溫需求；越冬期負積溫；北部冬麥區

氣候的任何波動和變化都會影響作物的生長[1]，明顯地反映在作物生育期的變化上[2-4]。Oteros等[4]研究指出，在氣候變化的影響下，西班牙大部分谷物的生育期提前。Tao等[5]分析了中國水稻、小麥、玉米生育期的變化發現，溫度顯著升高使作物生育期改變。已有研究表明，氣候變暖在冬、春季尤為明顯[6-7]，而冬小麥主要生長在10月-翌年6月，受氣候變化的影響較大[4]。多數研究發現，氣候變暖導致冬小麥冬前生育期推遲，冬后生育期提前[8-11]。Hu等[8]研究發現美國大平原6個站點冬小麥的抽穗或開花期均提前。Xiao等[12]研究表明，中國冬小麥播種期推遲，抽穗和成熟期提前。但也有研究指出，冬小麥播種和出苗期呈提前趨勢，返青期變化較復雜，抽穗和成熟期推遲[3, 13-14]。可見，在氣候變化的影響下，冬小麥生育期的變化規律非常復雜，且具有較強的區域特征。

冬小麥生育期的變化受品種、熱量和水分等因素的影響，而熱量因素尤為重要[2, 10, 15-16]。Xiao等[10, 16]研究表明，氣候變暖是冬小麥生育期變化的主要原因。積溫為熱量資源的指標之一，顯著影響了冬小麥的生育期。已有研究表明，冬后各生育期與冬季積溫呈顯著負相關[17-18]，春季積溫增加可導致抽穗或開花期提前[8]。冬小麥完成一個生育階段需要一定的熱量累積，生育階段內的積溫可以反映其熱量需求特性[2]。 Tao等[2, 14]研究指出，在氣候變化的影響下，冬小麥從播種到幾乎每個生育期內的積溫多呈增加趨勢，不同生育階段內也發生了改變，但目前關于積溫對冬小麥生育期的影響研究還未細化到生育階段。因此，有必要細致地分析各生育階段內積溫對冬小麥生育期的影響。

北部冬麥區為中國冬小麥最北部種植區，冬小麥生長期內光熱資源較為適宜，但氣候波動對冬小麥的生產影響較大[19]。已有研究表明，氣候變化顯著影響了該地區冬小麥的生長發育[14, 20]。因此，明確該地區冬小麥生育期的演變規律及影響要素，能更好地適應氣候變化，提高冬小麥的產量。

基于上述研究，本文利用中國北部冬麥區19個農業氣象站1993-2013年的冬小麥生育期數據，研究冬小麥各生育階段積溫和生育期的時空變異特征，并進行相關性分析，旨在為該區冬小麥的品種選育和種植決策提供科學依據。

1 材料與方法

1.1 研究區域

依據趙廣才[19]對中國小麥區劃的研究結果，選取北部冬麥區中除遼寧省外的區域為研究對象，橫跨河北、天津、北京、山西、陜西和甘肅6省市，為一條狹長地帶（圖1）。研究區內冬小麥占糧食作物種植面積的30%～40%，全年≥10℃的積溫為3500℃·d左右，最冷月平均氣溫為?10.7～?4.1℃，冬小麥生育期太陽輻射總量為276～293kJ·cm?2，播種-成熟期＞0℃積溫為2200℃·d左右，干旱、嚴寒是影響冬小麥生產的主要問題[19]。

圖1 研究區域及站點分布

1.2 數據來源

綜合考慮數據的完整性與連續性，選取研究區域內19個農業氣象站點21a（1993-2013年）的冬小麥生育期資料，包括播種、出苗、越冬開始、返青、拔節、抽穗、乳熟和成熟期，數據來源于中國氣象數據網的《中國農作物生長發育和農田土壤濕度旬值數據集》。氣象資料來源于中國氣象數據網的《中國地面氣候資料日值數據集》，僅使用平均氣溫（℃）資料。對于個別缺測的平均氣溫數據，利用五日滑動平均法插補[21]。因選用的氣候資料數據集并不完全具有所選農業氣象站的氣象資料，故采取就近匹配原則，即對于沒有氣象資料的農業氣象站，選用最近的氣象站（位于相鄰縣或市，氣候條件接近）作為替代站與其對應，見表1。

1.3 研究方法

1.3.1 核密度估計

核密度估計是一種能從樣本估計總體的非參數密度估計方法。一方面，應用核密度估計不需要引入對數據分布的先驗假設，只從樣本本身出發獲取數據分布特征，可以用來估計任意形狀密度估計的方法，另一方面，相比于直方圖等其它密度估計方法，通過核密度估計得到的概率密度連續性更好，并且不會依賴于選取的區間長度[22-23]。生育期的核密度測算依據研究站點1993-2013年的生育日期數據集來計算其集聚特征，核密度估算值越大，則生育期出現在該日期的概率相對越大。核密度的基本表達式為

式中，f(x)為獨立分布的數據x1，x2，…，xn估計得到的概率密度函數；K為核函數；h為帶寬。

表1 農業氣象站和對應的氣象站

1.3.2 標準差和變異系數

標準差（σ）可以反映組內個體間的離散程度，變異系數（cv）可以衡量波動程度，變異系數大，說明要素的穩定性差，波動性大。設數據系列為xi，則計算式為

1.3.3 變化趨勢分析

式中，a為回歸常數，以b的10倍作為氣候傾向率。采用Student的t-test 檢驗法對各要素變化趨勢進行顯著性檢驗，P＜0.05為顯著，P＜0.01為極顯著。

1.4 數據處理

以1月1日為1，將各生育日期轉換為日序。根據實測生育期和逐日平均氣溫計算各站逐年各生育階段積溫，其中，播種-出苗、出苗-越冬開始、返青-拔節、拔節-抽穗、抽穗-乳熟、乳熟-成熟和播種-成熟階段計算≥0℃積溫（以下簡稱為積溫），而越冬期計算該階段的負積溫（下同）。各站點的生育階段積溫（或生育期）為該站多年平均值，分析各站點的生育階段積溫數據集的變異系數（絕對值），以反映生育階段積溫的空間變異特征。由于生育期數據的缺失，使得分析1993-2013年冬小麥不同生育階段積溫與生育期相關性的樣本量在374～395，因其樣本量較大，進行相關性大小的比較時，假定其樣本量相同。采用R語言進行統計分析和相關性檢驗，使用ArcGIS和R語言進行圖形制作。

2 結果與分析

2.1 北部冬麥區小麥各生育階段積溫的時空變異特征

2.1.1 空間變異

比較變異系數可知，播種-成熟和抽穗-乳熟階段積溫的變異系數較小，而越冬期負積溫的變異系數最大（25.2%），表明播種-成熟和抽穗-乳熟階段積溫的空間差異較小，而越冬期負積溫的差異則較大，其它生育階段積溫居中（圖2）。在空間上，播種-出苗階段積溫的區域平均為122.9℃·d，高于平均值的站點多集中在麥區的東北部，低于平均值的站點多位于山西、陜西和甘肅境內，整體呈東高西低的空間分布；出苗-越冬開始階段積溫的區域平均為360.4℃·d，高于平均值的站點多位于河北的昌黎、黃驊、霸州和定州、北京通縣、山西長治、陜西洛川和甘肅西峰鎮一帶，呈東南高西北低的分布；越冬期負積溫的區域平均值為?258.8℃·d，高于平均值的站點多位于麥區東部，呈東高西低的分布；返青-拔節階段積溫的區域平均值為376.1℃·d，低于平均值的站點多位于麥區南部邊界的甘肅西峰、陜西洛川和山西長治一帶，呈東高西低的分布；拔節-抽穗階段積溫的區域平均值為334.6℃·d，高于平均值的站點多位于麥區西北部，呈東南低西北高的分布；抽穗-乳熟階段積溫的區域平均值為490.4℃·d，空間分布與拔節-抽穗階段一致；乳熟-成熟階段積溫的區域平均值為329.6℃·d，空間分布與拔節-抽穗階段一致；播種-成熟階段積溫的區域平均值為2088.6℃·d，低于平均值的站點多位于甘肅的平涼和涇川、山西長治、天津靜海、河北的河間和昌黎一帶，整體呈東南低西北高的分布。

圖2 1993-2013年研究區19個農氣站冬小麥各生育階段的積溫

注：cv為19個站點生育階段積溫多年平均值數據集的變異系數；越冬期為負積溫。下同。

Note: cv is coefficient of variation of the data set that mean of accumulated temperature during growth stages at the 19 stations. Negative accumulated temperature was calculated during the overwintering period. The same as below.

2.1.2 變化趨勢

從冬小麥生育階段積溫傾向率的區域平均值來看，越冬期負積溫和返青-拔節階段積溫顯著增加，抽穗-乳熟階段積溫呈極顯著增加，其余階段積溫變化則不顯著（圖3）。從各站點來看，冬小麥播種-出苗階段，11%的站點積溫變化顯著，均位于麥區東部；出苗-越冬開始積溫呈顯著減少和增加的站點均占11%，顯著減少的站點位于麥區的中、東部，顯著增加的站點均位于麥區的西部；42%的站點越冬期負積溫顯著增加，即越冬期變冷，且均位于麥區東部；返青?拔節階段，26%的站點積溫變化顯著，均呈增加趨勢且位于麥區東南部；拔節-抽穗階段，16%的站點積溫顯著增加且均位于麥區西部，21%的站點顯著減少且多位于麥區東部；抽穗-乳熟階段，47%的站點積溫變化顯著，37%的站點為顯著增加，且多位于麥區東部；乳熟-成熟階段，21%的站點積溫顯著減少且均位于麥區東部，11%站點顯著增加且均位于麥區中部的山西省；播種-成熟階段，21%的站點積溫顯著增加，11%站點顯著減少，顯著增加的站點均位于麥區東部，呈顯著減少的站點均位于麥區東北部。

圖3 1993-2013年研究區19個站點冬小麥各生育階段積溫的傾向率

注：*、**分別表示相關系數通過0.05、0.01水平的顯著性檢驗。下同。

Note:*is P＜0.05，**is P＜0.01.The same as below。

2.2 北部冬麥區小麥生育期的時空變異特征

2.2.1 基本特征

利用所有站點和年份生育期實際觀測資料繪制各生育期的核密度估計圖，并計算平均值及其標準差，結果見圖4。由圖可見，研究區內冬小麥各生育期出現日期均符合正態分布規律。其中，播種期大概率出現在10月初，平均日序為275（以1月1日為1，播期在10月2日）；出苗期大概率出現在10月中旬，平均日序為284（10月11日）；越冬開始大概率出現在12月初，平均日序為336（12月2日）；返青期大概率出現在3月上旬，平均日序為67（3月8日）；拔節期大概率出現在4月中旬，平均日序為108（4月18日）；抽穗期大概率出現在5月上旬，平均日序為129（5月9日）；乳熟期大概率出現在6月初，平均日序為154（6月3日）；成熟期大概率出現在6月中旬，平均日序為168（6月17日）。由于冬小麥生育期數據系列包括19個站點21a的資料，因此，標準差的數值大小包含了空間和時間兩方面的變化。從各生育期的標準差來看，越冬開始和返青期的標準差較大，分別為8.64d和9.48d，而拔節、抽穗期的標準差較小，分別為6.25d和6.07d。表明1993-2013年北部冬麥區冬小麥越冬開始和返青期的時空變異較大，拔節、抽穗期的時空變異較小，其它生育期變異程度居中。

圖4 研究區19個站點1993-2013年冬小麥實測生育期的核密度估計

注：σ為每個生育期19個站點21a實測數據系列的標準差。

Note: σ is standard deviation of each measured phenology data（including the data of 21 years at 19 stations）.

2.2.2 空間變異

計算各研究站生育日期的多年平均值和19個站點生育期多年平均值數據集的標準差，分析生育期的空間變異特征。由圖5可見，各站點播種期日序的多年平均值在263～285，出苗期在272～293，越冬開始在329～341，返青期在53～81，拔節期在103～117，抽穗期在122～141，乳熟期在146～163，成熟期在157～180。在空間上，播種和出苗期早于區域平均日期的站點主要在甘肅西峰、陜西洛川和山西長治地區，整體呈西早東晚的空間分布；越冬開始晚于區域平均日期的站點主要在甘肅西峰、陜西洛川、山西介休、河北河間、天津靜海一帶，整體呈東南晚西北早的空間分布；返青期早于區域平均日期的站點主要在甘肅平涼、山西汾陽、河北的定州和黃驊、天津寶坻地區，整體呈東南早西北晚的空間分布；拔節期早于區域平均日期的站點主要集中在麥區東部；而抽穗、乳熟和成熟期早于區域平均日期的站點主要位于麥區東部的河北、北京和天津境內，整體呈東早西晚的空間分布。比較標準差的大小可知，返青和成熟期的標準差較大，分別為6.52d和6.70d，而越冬開始和拔節期的標準差較小，分別為3.73d和3.19d。表明返青和成熟期的空間差異較大，越冬開始和拔節期則較小，其它生育期居中。

2.2.3 變化趨勢

從各生育期的變化趨勢來看（圖6），線性變化傾向率的區域平均值均為正值，且播種期顯著推遲（u=1.7d·10a?1, P<0.05）。各研究站點中，播種期，21%的站點顯著推遲，呈推遲趨勢的站點位于麥區中部的山西境內和東部臨海地區；出苗期，16%的站點顯著推遲，11%的站點呈顯著提前的趨勢，變化顯著的站點均位于麥區東南部；越冬開始期，僅5%的站點變化顯著，只在麥區東北部的北京通縣顯著推遲；返青、乳熟和成熟期分別有37%、42%、21%的站點變化顯著，均呈推遲趨勢且位于麥區東部；拔節期，26%的站點顯著推遲，其中11%的站點顯著提前，顯著推遲的站點均位于麥區東部，顯著提前的站點均位于西部；僅有5%的站點的抽穗期變化顯著，在麥區西南部的甘肅西峰顯著提前，多數站點的變化幅度較小。

圖5 1993-2013年研究區19個站點實測冬小麥的生育期（年日序）

注：σ為19個站點生育期多年平均值數據集的標準差。

Note: σ is standard deviation of the data set that mean date of winter wheat phenology at the 19 stations.

圖6 1993-2013年研究區19個站點冬小麥各生育期的線性傾向率

2.3 冬小麥生育期與生育階段積溫的相關性分析

由表2的相關性分析可知，播種-出苗階段積溫與出苗期呈極顯著正相關，而與乳熟和成熟期呈顯著負相關，與其它生育期的相關性不顯著，表明出苗、乳熟和成熟期受到了播種-出苗階段積溫變化的影響，播種-出苗階段積溫增加使出苗期推遲，乳熟和成熟期提前；出苗-越冬開始階段積溫與拔節、乳熟和成熟期呈極顯著負相關，與抽穗期呈顯著負相關，與返青期的相關性則不顯著，表明出苗-越冬開始階段積溫越大，拔節、抽穗、乳熟和成熟期將提前，而對返青期影響較小；越冬期負積溫與返青、拔節、抽穗、乳熟和成熟期均呈極顯著負相關，表明越冬期負積溫減少，小麥冬后生育期提前；返青-拔節階段積溫與拔節期呈極顯著正相關，與抽穗、乳熟和成熟期呈極顯著負相關，表明返青-拔節階段積溫增加使拔節期推遲，抽穗、乳熟和成熟期提前；拔節-抽穗階段積溫與抽穗、乳熟期呈極顯著正相關，與成熟期呈顯著正相關，表明拔節-抽穗階段積溫增加，則抽穗、乳熟和成熟期均推遲；抽穗-乳熟階段積溫與乳熟期呈極顯著正相關，與成熟期相關性不顯著，表明抽穗-乳熟階段積溫增加使乳熟期推遲，而成熟期受該階段熱量需求特性的影響較小；乳熟-成熟階段積溫與成熟期呈極顯著正相關，說明乳熟-成熟階段積溫增加導致成熟期推遲；播種-成熟階段積溫與成熟期相關性不顯著。就各生育期而言，出苗期受播種-出苗階段積溫的影響顯著，返青、拔節、抽穗、乳熟和成熟期均與越冬期負積溫的相關性最大，表明生育階段積溫的變化對冬小麥生育期產生了不同程度的影響，但冬后生育期受越冬期負積溫的影響程度更大。

分析各站點冬小麥越冬期負積溫與冬后各生育期的相關性（圖7）可知，冬小麥越冬期負積溫與返青、拔節、抽穗、乳熟和成熟期分別在79%、79%、47%、26%、42%的站點上呈顯著負相關，與返青和拔節期顯著相關的站點較多。在空間上，越冬期負積溫與抽穗、乳熟和成熟期相關性顯著的站點多位于麥區的東部，在西部站點上的相關性多不顯著，而與返青和拔節期的相關性在區域上的差異較小，普遍呈顯著負相關關系。

表2 冬小麥各生育階段積溫與生育期的相關性

Note: ATSE, ATEO, ATRJ, ATJH, ATHM, ATMMand ATSMrepresent ≥0℃ accumulated temperature during sowing to emergence, emergence to overwintering, regreening to jointing, jointing to heading, heading to milk-ripe, milk-ripe to maturity and sowing to maturity periods, respectively.

圖7 1993-2013年冬小麥越冬期負積溫與冬后各生育期的相關性

3 結論與討論

冬小麥各生育階段內積溫的變化影響了其生長發育進程，研究區內21%的站點拔節-抽穗、乳熟-成熟階段積溫均顯著減少，返青-拔節、抽穗-乳熟和播種-成熟階段積溫及越冬期負積溫分別在26%、37%、21%和42%的站點顯著增加，而播種-出苗和出苗-越冬階段積溫的變化幅度較小，且各生育階段積溫的變化多存在東西部區域差異。表明冬小麥的積溫需求結構發生改變，且在不同生育階段和區域的變化存在差異。就區域平均而言，麥區冬小麥各生育期在1993-2013年的傾向率均為正值，且播種期顯著推遲，與楊建瑩等[13]對華北地區和高輝明等[14]對北部冬小麥的研究結果一致。相關性分析表明，播種-出苗階段積溫與出苗、乳熟和成熟期，出苗-越冬開始階段積溫與拔節、抽穗、乳熟和成熟期的相關性均顯著；返青-拔節、拔節-抽穗和乳熟-成熟階段積溫及越冬期負積溫與相應階段之后生育期的相關性均顯著，且多呈極顯著相關；抽穗-乳熟階段積溫與乳熟期呈極顯著相關。表明生育階段積溫的變化直接影響下個生育期，同時也可能間接影響了之后多個生育期的發生。

冬小麥在播種-出苗階段對積溫的需求影響出苗日期。在各站點平均值的空間分布上，播種-出苗階段積溫在麥區東部較高，西部則較低；東部出苗晚，西部則較早。相關性分析表明，播種-出苗階段積溫與出苗呈顯著正相關關系。但在1993-2013年，冬小麥播種-出苗階段積溫和出苗期分別在11%和26%站點上變化顯著，播種-出苗階段積溫和出苗期在時間上的變化存在差異。表明播種-出苗階段積溫需求的改變將影響冬小麥的出苗，但播期、播種深度、溫度、土壤水分狀況等因素對出苗期仍有較大的影響[4, 24-25]。

對生育期與生育階段積溫進行相關性分析可知，返青、拔節、抽穗、乳熟和成熟期與越冬期負積溫的相關性最大，且均呈極顯著負相關。在空間分布上，麥區東部站點越冬期負積溫的多年平均值多大于區域平均，拔節、抽穗、乳熟和成熟期早于區域平均日期的站點主要位于麥區東部的河北、北京和天津境內，整體均存在東西部差異。在時間尺度上，越冬期負積溫在42%的站點上顯著增加，且均位于麥區東部，西部站點變化不顯著；麥區東部多個站點冬小麥的返青、乳熟和成熟期呈顯著推遲趨勢，西部站點的變化則多不顯著；拔節期顯著推遲的站點位于麥區東部，顯著提前的站點位于麥區西部。可以看出，越冬期負積溫與冬后多個生育期呈現一致的時空變異特征。已有研究指出，越冬期負積溫影響冬小麥的生理代謝和形態建成，負積溫越多，冬小麥越易發生凍害，根系吸水慢，發芽遲，影響返青和拔節成穗[26-29]。綜上所述，越冬期負積溫影響冬小麥冬后的生長發育，其時空變異性可能是造成冬后生育期在時空上存在差異的原因。

越冬期負積溫與返青、拔節、抽穗、乳熟和成熟期的相關性在空間上差異較大。越冬期負積溫與抽穗、乳熟和成熟期顯著相關的站點多位于麥區東部，在西部站點上的相關性則多不顯著，與返青和拔節期則普遍呈顯著負相關關系。因此，麥區東部冬小麥越冬期負積溫普遍顯著增加可能是造成返青、拔節、乳熟和成熟期推遲的原因。此外，為適應氣候變暖而使用冬性減弱、熱量需求增加、開花與成熟較晚的品種[2, 30]，可能也導致了冬后生育期的推遲。但東部站點的抽穗期變化均不顯著，可能是由于越冬期負積溫增加使抽穗期推遲，而返青-拔節階段積溫增加和拔節-抽穗階段積溫減少使抽穗期提前共同作用的結果。西部多數站點的越冬期負積溫與冬后生育期相關性不顯著，一方面可能是西部站點的越冬期負積溫變化幅度較小導致，另一方面也說明這些相關性不顯著站點的抽穗、乳熟和成熟期變化的主要因素不是越冬期負積溫。

冬小麥各生育階段積溫與冬后多個生育期的相關性顯著，但相比于越冬期負積溫，其相關性均較低。這可能是由于在各生育階段內，冬小麥對積溫的需求影響其生長發育，同時，生育期在很大程度上可能也受生長前期的熱量條件、光照、土壤水分狀況及施肥等因素的影響[31]，特別是冬小麥在返青-拔節階段受光周期效應的影響尤為顯著[25]，拔節-乳熟階段又是冬小麥的需水關鍵期[32]。冬小麥在抽穗-乳熟階段主要進行灌漿，此階段的積溫與乳熟期呈極顯著正相關。并且，東部多數站點抽穗-乳熟階段積溫顯著增加，乳熟期也多顯著推遲。表明抽穗-乳熟階段的熱量需求增加使乳熟期推遲，乳熟較抽穗期推遲幅度大使抽穗-乳熟階段延長，延長了灌漿時間。Fernández等[33]研究指出，灌漿期的長度與千粒重呈正相關關系。可知，抽穗-乳熟階段積溫的增加有利于提高千粒重，增加產量。因此，選育抽穗-乳熟階段積溫需求高的品種可能是適應氣候變化的有效措施[2, 14]。

由于本研究收集的冬小麥生育期觀測站點及其資料的時間序列有限，分析結果可能與實際有差異。此外，冬小麥生育期的變化受光、溫、水等多種氣候要素的影響，而本文只進行了生育階段積溫的影響分析，在未來應進一步分析生育期對多種要素變化及要素間互作效應的響應，并探究冬小麥不同生育期階段積溫的變化對其品質和產量的影響。

[1]Wheeler T,Von Braun J.Climate change impacts on global food security[J].Science,2013,341(6145):508-513.

[2]Tao F,Zhang S,Zhang Z.Spatiotemporal changes of wheat phenology in China under the effects of temperature,day length and cultivar thermal characteristics[J].European Journal of Agronomy,2012,43:201-212.

[3]Estrella N,Sparks T H,Menzel A.Trends and temperature response in the phenology of crops in Germany[J].Global Change Biology,2007,13(8):1737-1747.

[4]Oteros J,García-Mozo H,Botey R,et al.Variations in cereal crop phenology in Spain over the last twenty-six years (1986-2012)[J].Climatic Change,2015,130(4):545-558.

[5]Tao F L,Yokozawa M,Xu Y L,et al.Climate changes and trends in phenology and yields of field crops in China,1981- 2000[J].Agricultural and Forest Meteorology, 2006,138(1): 82-92.

[6]IPCC.Climate change 2013:the physical science basis[M]. Cambridge,U K:Cambridge University Press,2013.

[7]丁一匯,任國玉,石廣玉,等.氣候變化國家評估報告I:中國氣候變化的歷史和未來趨勢[J].氣候變化研究進展,2006,2(1): 3-8,50.

Ding Y H,Ren G Y,Shi G Y,et al.National assessment report of climate change I:climate change in China and its future trend[J].Advance in Climate Change Research,2006,2(1): 3-8,50.(in Chinese)

[8]Hu Q,Weiss A,Feng S,et al.Earlier winter wheat heading dates and warmer spring in the US Great Plains[J].Agricultural and Forest Meteorology,2005,135(1):284-290.

[9]Croitoru A E,Holobaca I H,Lazar C,et al.Air temperature trend and the impact on winter wheat phenology in Romania[J]. Climatic Change,2012,111(2):393-410.

[10]Xiao D P,Tao F L,Liu Y J,et al.Observed changes in winter wheat phenology in the North China Plain for 1981- 2009[J]. International Jurnal of Bometeorology,2013,57(2): 275-285.

[11]姬興杰,朱業玉,劉曉迎,等.氣候變化對北方冬麥區冬小麥生育期的影響[J].中國農業氣象,2011,32(4):576-581.

Ji X J,Zhu Y Y,Liu X Y,et al.Impacts of climate change on the winter wheat growth stages in North China[J].Chinese Journal of Agrometeorology,2011,32(4):576-581.(in Chinese)

[12]Xiao D,Moiwo J P,Tao F,et al.Spatiotemporal variability of winter wheat phenology in response to weather and climate variability in China[J].Mitigation and Adaptation Strategies for Global Change,2015,20(7):1191-1202.

[13]楊建瑩,梅旭榮,劉勤,等.氣候變化背景下華北地區冬小麥生育期的變化特征[J].植物生態學報,2011,35(6):623-631.

Yang J Y,Mei X R,Liu Q,et al.Variations of winter wheat growth stages under climate changes in Northern China[J]. Chinese Journal of Plant Ecology,2011,35(6): 623-631.(in Chinese)

[14]高輝明,張正斌,徐萍,等.2001-2009年中國北部冬小麥生育期和產量變化[J].中國農業科學,2013,46(11):2201-2210.

Gao H M,Zhang Z B,Xu P,et al.Changes of winter wheat growth period and yield in Northern China from 2001-2009[J].Scientia Agricultura Sinica,2013,46(11):2201- 2210.(in Chinese)

[15]Mc Master G S,White J W,Hunt L,et al.Simulating the influence of vernalization,photoperiod and optimum temperature on wheat developmental rates[J].Annals of Botany,2008,102(4):561-569.

[16]Li K N,Yang X G,Tian H Q,et al.Effects of changing climate and cultivar on the phenology and yield of winter wheat in the North China Plain[J].International Journal of Biometeor- ology,2016,60(1):21-32.

[17]樊曉春,郭江勇,楊小利.西峰黃土高原冬季積溫變化對作物發育期的影響[J].中國農業氣象,2007,28(3):318-321.

Fan X C,Guo J Y,Yang X L.Influence of changes of accumulation temperature in winter oncrop growing periods in Loess Plateau of Xifeng[J].Chinese Journal of Agrometeor- ology,2007,28(3):318-321.(in Chinese)

[18]閆小珍,張隨賢,苗國柱,等.小麥發育期對冬季積溫變化的響應分析[J].氣象與環境科學,2009,32(4):37-39.

Yan X Z,Zhang S X,Miao G Z,et al.Analysis in response to changes of accumulated temperature in winter on growing period of winter wheat[J].Meteorological and Environmental Sciences,2009,32(4):37-39.(in Chinese)

[19]趙廣才.中國小麥種植區劃研究(一)[J].麥類作物學報,2010,30(5):886-895.

Zhao G C.Study on chinese wheat planting regionalization I[J].Journal of Triticeae Crops,2010,30(5):886-895.(in Chinese)

[20]王斌,顧蘊倩,劉雪,等.中國冬小麥種植區光熱資源及其配比的時空演變特征分析[J].中國農業科學,2012,45(2):228-238.

Wang B,Gu Y Q,Liu X,et al.Analysis of the temporal changes of photo-thermal resources in winter wheat growing regions in China[J].Scientia Agicultura Sinica,2012,45(2): 228-238.(in Chinese)

[21]楊帆,劉布春,劉園,等.氣候變化對東北玉米干旱指數保險純費率厘定的影響[J].中國農業氣象,2015,36(3):346-355.

Yang F,Liu B C,Liu Y,et al.Impact of climate change on pure premium rating of drought index insurance for maize in Northeast China[J].Chinese Journal of Agrometeorology, 2015,36(3):346-355.(in Chinese)

[22]張桂銘,朱阿興,楊勝天,等.基于核密度估計的動物生境適宜度制圖方法[J].生態學報,2013,33(23):7590-7600.

Zhang G M,Zhu A X,Yang S T,et al.Mapping wildlife habitat suitability using kernel density estimation[J].Acta Ecologica Sinica,2013,33(23):7590-7600.(in Chinese)

[23]謝國偉,錢雪忠,周世兵.基于非參數核密度估計的密度峰值聚類算法[J].計算機應用研究,2018,35(10):1-6.

Xie G W,Qian X Z,Zhou S B.Density peak clustering algorithm based on non-parametric kernel density estimation[J]. Application Research of Computers,2018, 35(10):1-6.(in Chinese)

[24]Menzel A,von Vopelius J,Estrella N,et al.Farmers’ annual activities are not tracking the speed of climate change[J]. Climate Research,2006,32(3):201-207.

[25]嚴美春,曹衛星,羅衛紅,等.小麥發育過程及生育期機理模型的研究I.建模的基本設想與模型的描述[J].應用生態學報,2000,11(3):355-359.

Yan M C,Cao W X,Luo W H,et al.Validation and evaluation of a mechanistic model of phasic and phenological development of wheat[J].Scientia Agricultura Sinica,2000, 11(3):355-359.(in Chinese)

[26]宋吉仁.冬小麥越冬期生理指標的比較及蔗糖代謝酶基因的表達[D].吉林:東北農業大學,2014.

Song J R.Comparison of some physiological traits in over-wintering stage and genes expression in sugar metabolism in winter wheat[D].Jilin:Northeast Agricultural University,2014.(in Chinese)

[27]莫志鴻,霍治國,葉彩華,等.北京地區冬小麥越冬凍害的時空分布與氣候風險區劃[J].生態學雜志,2013,32(12): 3197-3206.

Mo Z H,Huo Z G,Ye C H,et al.Spatial-temporal distribution and climatic risk zonation of freeze injury of winter wheat during overwintering stage in Beijing region[J].Chinese Journal of Ecology,2013,32(12):3197-3206.(in Chinese)

[28]王慧芳,顧曉鶴,董瑩瑩,等.冬小麥凍害災情及長勢恢復的變化向量分析[J].農業工程學報,2011,27(11):141-150.

Wang H F,Gu X H,Dong Y Y,et al.Monitoring freeze injury and growth recovery of winter wheat based on change vector analysis[J].Transactions of the CSAE,2011,27(11):141-150. (in Chinese)

[29]張愛芝,高志強,衛麗萍,等.北移冬小麥生長發育與氣象條件關系的研究[J].中國農業氣象,2004,25(1):27-29.

Zhang A Z,Gao Z Q,Wei L P,et al.Studies on relationship between growth and development of northward winter wheat and meteorological conditions[J].Chinese Journal of Agrometeorology,2004,25(1):27-29.(in Chinese)

[30]He L,Asseng S,Zhao G,et al.Impacts of recent climate warming,cultivar changes,and crop management on winter wheat phenology across the Loess Plateau of China[J]. Agricultural and Forest Meteorology,2015, 200(1):35-43.

[31]俄有浩,霍治國,馬玉平,等.中國北方春小麥生育期變化的區域差異性與氣候適應性[J].生態學報,2013,33(19): 6295-6302.

E Y H,Huo Z G,Ma Y P,et al.The regional diversity of changes in growing duration of spring wheat and its correlation with climatic adaptation in Northern China[J]. Acta Ecologica Sinica,2013,33(19):6295-6302.(in Chinese)

[32]信乃詮,崔讀昌,高亮之,等.中國農業氣象學[M].北京:中國農業出版社,1999.

Xin N Q,Cui D C,Gao L Z,et al.China agricultural meteorology[M].Beijing:China Agriculture Press,1999.

[33]Fernández C,Tewolde H,Elledge R,et al.Early-heading winter wheat varieties secure higher yields in southwest texas.BL-RC18[J].Bottom Lines:An Overview of Research and Extension Projects Texas A&M University Agricultural Research and Extension Center The Texas Agricultural Experiment Station The Texas A&M University System, Uvalde,TX,1998.

Responses of Winter Wheat Phenology to Accumulated Temperature during Growing Periods in Northern China Wheat Belt

MA Qian-qian，HE Yong，ZHANG Meng-ting，ZHANG Cong，XU Yin-long

（Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081, China）

Based on data collected from 19 agro-meteorological stations located in the winter wheat belt of Northern China from 1993 to 2013, this study analyzed temporal and spatial variability of winter wheat phenology and accumulated temperature in each growing period. Pearson correlation analysis and other methods were used to analyze the effects of accumulated temperature in each growing period on winter wheat phenology. The results showed that,（1）the spatial distribution of ≥0℃accumulated temperature during sowing to emergence and regreening to jointing periods along with the negative accumulated temperature value during the overwintering period (NATop) increased from east to west, and ≥0℃ accumulated temperature during jointing to heading, heading to milk-ripe, milk-ripe to maturity and sowing to maturity periods increased from southeast to northwest, while the distribution during emergency to start of overwintering period was contrary. The ≥0℃accumulated temperature for jointing to heading and milk-ripe to maturity periods decreased significantly at 21% of the investigated stations. The value of NATop, ≥0℃accumulated temperature during regreening to jointing, heading to milk-ripe and sowing to maturity periods increased significantly at 42%, 26%, 37% and 21% of the investigated stations, respectively. The variations of ≥0℃ accumulated temperature during sowing to emergency and emergency to the start of overwintering periods were much smaller. (2) Sowing and emergency date in the east was later than the counterpart in the west, of which variation trend was contrary to heading, milk-ripe and maturity date. Start of overwintering date in the southeast was later than the counterpart in the northwest, of which variation trend was contrary to regreening date. The stations of earlier jointing date were mainly located in the east. Sowing, emergency, regreening, jointing, milk-ripe and maturity date delayed significantly at 21%, 16%, 37%, 26%, 42% and 21% of the investigated stations, respectively. And most of these stations were in the east of the study area. Overwintering and heading date changed significantly only at 5% of the investigated stations. (3) Correlation analysis showed that the correlation between ≥0℃ accumulated temperature (or NATop) and multiple development stages was significant, which indicated that the growth and development of winter wheat might be directly or indirectly influenced by the accumulated temperature during growing periods. The NATop had the highest correlation with regreening, jointing, heading, milk-ripe and maturity date, and showed consistent spatiotemporal variation characteristics with multiple post-winter development stages. Spatiotemporal variability of the NATop might be the factor that caused spatiotemporal variations of winter wheat post-winter phenology.

Winter wheat; Phenology; Accumulated temperature requirement; Negative accumulated temperature during the overwintering period; Northern China winter wheat belt

10.3969/j.issn.1000-6362.2018.04.003

馬倩倩,賀勇,張夢婷,等.中國北部冬麥區小麥生育期對生育階段積溫變化的響應[J].中國農業氣象,2018,39(4):233-244

2017-08-22

。E-mail: heyong01@caas.cn; xuyinlong@caas.cn

國家科技支撐計劃項目(2013BAC09B04)；國家基金委青年科學基金項目(41501118)

馬倩倩（1992-），女，碩士生，研究方向為氣候資源與氣候變化。E-mail: maer12366@163.com