摘要：利用衡水市3個氣象站1965-2012年的逐日氣象數(shù)據(jù)，采用作物需水量模型以及相關統(tǒng)計方法，分析了近48年來衡水地區(qū)冬小麥灌漿期至成熟期需水量的變化特征。結果表明，1965-2012年衡水地區(qū)冬小麥生長后期需水量呈波動下降的趨勢，趨勢變化顯著，近48年來下降了37.9 mm，人工灌水量的減少趨勢快于需水量的變化，近48年來減少了62.4 mm；需水量存在時間序列上的突變現(xiàn)象，突變發(fā)生在1976年，1977年以后變化更加明顯；在冬小麥后期需水量波動減少的過程中存在明顯的7年變化周期；冬小麥后期生長需水量的變化與日照時數(shù)呈最大正相關，與相對濕度的變化呈最大負相關。

關鍵詞：冬小麥后期生長；作物需水量模型；變化

中圖分類號：S162.8；S512.1+1 文獻標識碼：A 文章編號：0439-8114（2016）14-3546-03

DOI：10.14088/j.cnki.issn0439-8114.2016.14.003

Abstract：Based on the daily meteorological data of three stations from 1965 to 2012 in Hengshui area， the change characteristics of winter wheat’s water requirement during filling stage to mature stage in 48 years of Hengshui area was analyzed by using the model of crop water requirement and the related statistical methods. The results showed a significant trend of reduction about the water requirement during the later growth stage of winter wheat from 1965 to 2012 in Hengshui area， it had 37.9 mm dropped in the past 48 years，the reduction of artificial irrigation water was more faster， and it had 62.4 mm dropped in the past 48 years. The water requirement existed a jumping phenomenon in the time series of recent 48 years， the time of jumping phenomenon appeared in 1976， and it was more obvious after 1977. In the fluctuating downward process of winter wheat’s water requirement during the later growth stage， it presented a variational period of 7 years. Sunshine time was the most positive correlation about the change of winter wheat’s water requirement during the later growth stage， the relative humidity was the most negative correlation to it.

Key words： the later growth stage of winter wheat； the model of crop water requirement； change

氣候變暖已是不爭的事實，政府間氣候變化專門委員會（IPCC）第五次會議上指出[1]，近130年來地球表面平均溫度上升了0.85 ℃，與IPCC第四次評估報告相比[2]，地球表面溫度又升高了0.11 ℃。氣候的變暖對全球生態(tài)系統(tǒng)產生了極大影響，而農業(yè)是受氣候變化影響最直接的脆弱行業(yè)[3-8]。衡水市是農業(yè)大市，耕地面積56.7萬hm2，農業(yè)人口人均占有耕地1 667.5 m2。在衡水市的農業(yè)種植結構中，主要以冬小麥-夏玉米輪作結構為主，在其主要作物之一的冬小麥中，其后期生長的抽穗至成熟期是冬小麥子粒灌漿成型的關鍵期，這段時期內的需水量顯得尤為重要。在氣候變暖的大背景下，華北地區(qū)氣溫升高尤為明顯[9]，衡水市的氣溫隨著大環(huán)境的升高而升高，在這種情況下，衡水地區(qū)冬小麥后期管理的需水量變化目前尚無研究，本研究以作物需水量模型公式為基礎，討論了近48年來衡水市冬小麥后期需水量的變化規(guī)律，以期為該段種植時間內冬小麥的需水管理提供合理化的理論指導。

1 資料與方法

1.1 資料來源

采用河北省氣象信息中心提供的衡水市冀州、故城、饒陽3個氣象站1965-2012年逐日氣溫（含最高、最低）、降水、日照、風速、相對濕度等氣象要素，取3個氣象站氣象要素的區(qū)域平均值作為衡水地區(qū)的逐日氣象資料。

1.2 試驗方法

采用作物系數(shù)法[10]計算冬小麥后期生長需水量的模型，計算公式如下：

ETc=Kc×ET0 （1）

式中，ETc為作物需水量（mm）；Kc為作物系數(shù)；ET0為潛在蒸散量（mm）。本研究中，Kc依照作物需水量計算指南（FAO-56）的建議數(shù)值[11]取0.4。

潛在蒸散量ET0的計算公式如下：

各類參數(shù)依照世界糧農組織（FAO）的推薦值[11]進行計算。式中，Rn表示冠層表面凈輻射；T表示日平均氣溫；es表示飽和水汽壓；ea表示實際水汽壓； γ表示濕度計常數(shù)；U2表示2 m高處的日平均風速；Δ是飽和水汽壓-氣溫關系曲線在T處的切線斜率；G表示土壤熱通量，本計算中忽略不計[12]。

由于抽穗至成熟期的起止日期不固定，為了更好地進行統(tǒng)計與計算，將該段日期固定為5月5日至6月10日，在日期的計算過程中采用日序的計數(shù)方式[13]，即1月1日記為1，2月2日記為33。統(tǒng)計方法采用線性趨勢、Mann-Kendall突變檢驗法、Morlet小波分析等方法[14，15]。

2 結果與分析

2.1 冬小麥生長后期需水量變化趨勢

對冀州、故城、饒陽3個氣象站的逐日氣象要素求區(qū)域平均值，得到衡水地區(qū)的氣象要素逐日平均值，通過公式（1）和（2）計算得到近48年來冬小麥后期需水量的時間序列（圖1），由圖1可知，1965-2012年，衡水地區(qū)冬小麥抽穗至成熟期平均需水量為141.0 mm，需水量最多出現(xiàn)在1971年，需水量為184.4 mm，最少出現(xiàn)在1998年，需水量為107.5 mm，近48年的時間變化呈波動下降的趨勢，變化量為-0.79 mm/年，下降趨勢通過了0.01的顯著性檢驗，變化明顯， 48年來下降了37.9 mm。

2.2 冬小麥生長后期灌水量變化趨勢

用同期的需水量減去降水量代表冬小麥后期的灌水量，由圖2可知，冬小麥后期的人工灌水量同樣呈波動下降的變化趨勢，灌水量最少的年份是1977年，由于該年降雨量大，降水量大于需水量，該年的灌水量為-71.6 mm，本研究中該年的灌水量取值為0 mm，表明該年的冬小麥生長后期自然降水完全能夠滿足冬小麥的后期生長，不需要再進行人工灌溉；灌溉水量最多年出現(xiàn)在1968年，灌水量為174.8 mm，表明近48年里，該年在冬小麥生長后期的自然降水最少；1965-2012年衡水地區(qū)冬小麥生長后期多年平均灌水量為98.9 mm，變化趨勢為-1.3 mm/年，近48年來下降了62.4 mm，其相關系數(shù)r=0.435，通過了0.01的顯著性檢驗，表明下降趨勢明顯。

2.3 冬小麥后期生長需水量時間序列變化的突變分析

由衡水地區(qū)冬小麥后期生長需水量時間變化序列得出冬小麥后期生長的需水量呈較為明顯的下降趨勢，為了進一步分析在近48年的時間序列中，受綜合氣候因素的影響冬小麥后期生長的需水量變化是否存在突變，本研究引入Mann-Kendall突變檢驗方法，該方法是一種非參數(shù)統(tǒng)計檢驗方法，可分析時間序列資料的變化趨勢和突變點[14]。其正序曲線UF與逆序曲線UB的交叉點為突變點，通過了±1.96（0.05的顯著性檢驗），表明變化趨勢顯著[14]。由圖3可知，衡水地區(qū)冬小麥生長后期的需水量在1976年發(fā)生突變現(xiàn)象，在1977年變化趨勢通過了0.05的顯著性檢驗，變化更加明顯，但由該圖可以看出UF線在2000年以后變化又趨于平緩。

2.4 冬小麥后期生長需水量時間序列的周期性變化分析

由圖1的需水量曲線可知，冬小麥后期生長需水量的變化雖呈遞減的趨勢，但不是單調遞減，而是呈波動變化的趨勢，因此，用Morlet小波實部結合其方差分析時間序列的周期性變化。小波分析源于傅立葉分析，是對傅立葉分析方法的突破性發(fā)展，通過將時間序列分解到時間頻率域內，得出時間序列顯著的周期變化動態(tài)以及周期變化動態(tài)的時間格局[14]。由圖4可知，冬小麥后期生長需水量存在準6～8年的變化性周期，同時結合小波方差圖（圖5）可見，其變化周期表現(xiàn)為7年，即需水量在近48年的時間序列變化中，存在7年增減交替變化的周期。

2.5 冬小麥后期生長需水量與氣候因子的相關性分析

隨著氣候的變暖，冬小麥后期的需水量卻呈減少的趨勢，因此對大田環(huán)境的主要氣象因子與冬小麥后期需水量的變化進行相關性分析，用以得出影響需水量變化最大的氣象因子。利用SPSS統(tǒng)計軟件對冬小麥生長后期需水量與平均氣溫、風速、日照時數(shù)、相對濕度等氣象要素進行相關分析，由各氣象因子與需水量的相關系數(shù)（表1）可以看出，與冬小麥生長后期需水量正相關最大的氣象因子是日照時數(shù)，其次是風速，再次是平均氣溫，相關系數(shù)全部通過了0.01的顯著性檢驗；負相關最大的是相對濕度，相關系數(shù)通過了0.01的顯著性檢驗。該結論也與河北日照時數(shù)減少[16]、風速減小[17]、半干旱區(qū)相對濕度變化不明顯[18]的結論一致。

3 小結與討論

衡水地區(qū)冬小麥生長后期需水量呈波動下降的趨勢，且趨勢變化明顯，該變化趨勢也與中國大部分地區(qū)潛在蒸散量呈減少的趨勢相一致[19]。1965-2012年需水量趨勢變化量減少了37.9 mm，需水量的降低和自然降水的作用使人工灌水量在近48年的減少趨勢明顯快于需水量的變化，近48年來，人工灌水量減少了62.4 mm。

冬小麥生長后期的需水量存在時間序列上的突變現(xiàn)象，突變交點發(fā)生在1976年，在1977年以后需水量的變化更加明顯，但這種變化在2000年以后有所放緩。冬小麥生長后期的需水量在波動下降的變化中存在增-減的周期性變化，其7年的變化周期在檢測中表現(xiàn)得尤為明顯。冬小麥生長后期需水量的變化與日照時數(shù)呈最大正相關，其次是風速，與相對濕度的變化呈最大負相關。該相關性也表明了其與氣候因子變化的一致性。

冬小麥生長后期是子粒產量形成的重要生長階段，其對水分的需求就顯得尤為敏感，本研究在相關模型和數(shù)據(jù)的基礎上計算發(fā)現(xiàn)，衡水地區(qū)冬小麥生長后期需水量呈減少的趨勢，這對該地區(qū)冬小麥后期管理以及合理利用水資源有一定的指導意義。

參考文獻：

[1] IPCC. IPCC Fifth Assessment Report（AR5）[M].Berne：Switzerland Press，2013.4-5.

[2] IPCC. IPCC Forth Assessment Report（AR4）[M]. Cambridge： Cambridge University Press，2007.1-12.

[3] 高素華，郭建平.提高中國三北地區(qū)氣候資源利用率的對策研究[J].資源科學，1999，21（4）：51-54.

[4] 紀瑞鵬，張玉書，馮 銳，等.遼寧省農業(yè)氣候資源變化特征分析[J].資源科學，2007，29（2）：74-82.

[5] 葉彩華，欒慶祖，胡 寶.北京農業(yè)氣候資源變化特征及其對不同種植模式玉米各生育期的影響[J].自然資源學報，2010，25（8）：1350-1364.

[6] 劉志娟，楊曉光，王文峰，等.氣候變化背景下我國東北三省農業(yè)氣候資源變化特征[J].應用生態(tài)學報，2009，20（9）：2199-2206.

[7] 楊曉光，劉志娟，陳 阜.全球氣候變暖對中國種植制度可能影響Ⅰ.氣候變暖對中國種植制度北界和糧食產量可能影響的分析[J].中國農業(yè)科學，2010，43（2）：329-336.

[8] 李克南，楊曉光，劉志娟，等.全球氣候變化對中國種植制度可能影響分析Ⅲ.中國北方地區(qū)氣候資源變化特征及其對種植制度界限的可能影響[J].中國農業(yè)科學，2010，43（10）：2088-2097.

[9] 馬潔華，劉 園，楊曉光，等.全球氣候變化背景下華北平原氣候資源變化趨勢[J].生態(tài)學報，2010，30（14）：3818-3827.

[10] 劉 鈺，PEREIRA L S.對FAO推薦的作物系數(shù)計算方法的驗證[J].農業(yè)工程學報，2000，16（5）：26-30.

[11] FAO. Evapotranspiration from a reference surface[DB/OL].http：//www.fao.org/nr/water/docs/ReferenceManualV32.pdf，2012-09.

[12] 吳紹洪，尹云鶴，鄭 度，等.近30年中國陸地表層干濕狀況研究[J].中國科學D輯地球科學，2005，35（3）：276-283.

[13] 吳云龍，李麗平，劉勝堯，等.石家莊地區(qū)冬小麥適播期時空變化分析[J].中國生態(tài)農業(yè)學報，2015，23（6）：768-773.

[14] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術[M].北京：氣象出版社，1999.69-72.

[15] 黃嘉佑.氣象統(tǒng)計與預報方法[M].北京：氣象出版，2004.28-30.

[16] 郭艷嶺，邱新法，張素云.1965-2005年河北日照時數(shù)時空分布特征及影響因子[J].干旱氣象，2010，28（3）：297-303.

[17] 張素云，周順武，郭艷嶺，等.近47a河北風速季節(jié)變化特征及影響因子分析[J].干旱區(qū)資源與環(huán)境，2010，24（12）：93-99.

[18] 靳英華，廉士歡，周道瑋，等.全球氣候變化下的半干旱區(qū)相對濕度變化研究[J].東北師大學報（自然科學版），2009，41（4）：134-138.

[19] 高 歌，陳德亮，任國玉，等.1956-2000年中國潛在蒸散量變化趨勢[J].地理研究，2006，25（3）：378-387.