張怡　史本林

摘要：分析了豫東冬小麥產(chǎn)區(qū)近50年的水、溫要素的變化趨勢，通過數(shù)學(xué)統(tǒng)計與模型模擬農(nóng)業(yè)氣候資源的變化和極端氣候事件發(fā)生的規(guī)律，預(yù)測未來農(nóng)業(yè)氣候資源變化趨勢和極端天氣事件出現(xiàn)的規(guī)律，探討氣候變化和干旱災(zāi)害對豫東地區(qū)冬小麥產(chǎn)量的影響，結(jié)果表明，豫東冬小麥產(chǎn)區(qū)氣候趨于“暖濕化”，對冬小麥產(chǎn)量的提升較為有利；冬小麥生長中后期的光照、溫度和降水適宜度與冬小麥產(chǎn)量的相關(guān)性顯著；豫東冬小麥產(chǎn)區(qū)干旱受災(zāi)面積呈減小趨勢，但成災(zāi)面積呈增大趨勢。

關(guān)鍵詞：氣候變化；冬小麥；豫東；極端天氣；干旱；產(chǎn)量

中圖分類號： S162.5文獻標(biāo)志碼： A文章編號：1002-1302（2015）02-0336-04

收稿日期：2014

基金項目：國家自然科學(xué)基金（編號：41140019）；河南省科技發(fā)展計劃[編號：商（20093004）、商（20103012）]。

作者簡介：張怡（1981—），男，河南柘城人，碩士，工程師，主要從事氣象預(yù)報預(yù)測與災(zāi)害防護的研究。Tel：（0370）3235767；E-mail：zhangyi8733873@sina.com。近50年中國平均氣溫變暖趨勢達到每50年0.6～1.1 ℃，升溫速率可達每10年0.25 ℃，大于北半球平均的升溫速率。氣候變化下，農(nóng)業(yè)氣象災(zāi)害發(fā)生頻次增加，極端氣候事件（洪水、干旱、極端高溫和低溫冷害等）對未來農(nóng)業(yè)生產(chǎn)的影響更大，使農(nóng)業(yè)種植結(jié)構(gòu)和布局改變，局部地區(qū)的農(nóng)業(yè)氣象災(zāi)害事件加劇，農(nóng)業(yè)生產(chǎn)的不穩(wěn)定性增加，就中國平均狀況而言，氣候變化將對土地生產(chǎn)潛力產(chǎn)生不利影響，如不采取應(yīng)對措施，氣候變化將對農(nóng)業(yè)產(chǎn)量、農(nóng)業(yè)經(jīng)濟造成損害。目前，氣候變化引起的冬小麥生育期極端低溫和干熱風(fēng)等極端氣候事件頻發(fā)，已成為限制冬小麥穩(wěn)產(chǎn)增產(chǎn)的主要影響因素。同時，主要氣象因素如溫度和降水的區(qū)域差異較為突出，不同區(qū)域社會經(jīng)濟環(huán)境的不同導(dǎo)致區(qū)域應(yīng)對能力的不同也會使氣候變化對農(nóng)業(yè)的影響產(chǎn)生差異[1-15]。

隨著我國糧食消費需求的不斷增長、耕地面積逐漸減少、水資源嚴(yán)重短缺，氣候因素變化對糧食生產(chǎn)的約束日益凸現(xiàn)，我國糧食安全面臨著嚴(yán)峻挑戰(zhàn)。冬小麥?zhǔn)俏覈谝淮蠹Z食作物，在保障我國糧食安全中的地位非常突出。黃淮海地區(qū)是我國冬小麥重要產(chǎn)區(qū)，該區(qū)冬小麥播種面積和產(chǎn)量約占全國的40%。受大陸性季風(fēng)氣候的影響，在冬小麥生育期間，干旱、霜凍、低溫冷害、 抽穗開花期陰雨、干熱風(fēng)等多種氣象災(zāi)害發(fā)生頻繁，嚴(yán)重制約冬小麥的高產(chǎn)穩(wěn)產(chǎn)[14-20]。因此我們探討了豫東地區(qū)氣候變化趨勢及其對冬小麥產(chǎn)量的影響、又特別對影響最為嚴(yán)重的干旱災(zāi)害進行了特征分析，以期能夠更為系統(tǒng)地了解此種變化對豫東冬小麥的綜合影響。

1材料與方法

本研究以豫東地區(qū)商丘市所轄8個國家氣象觀測站（商丘、永城、夏邑、睢縣、寧陵、柘城、民權(quán)、虞城、）1961—2010年共計50年逐日、旬、月、年的溫度、降水等觀測數(shù)據(jù)為信息源，并進行了質(zhì)量控制與訂正，保證資料的同一性和延續(xù)性。分別以四季的中間月份代表相應(yīng)季節(jié)，以逐年同期距平法分析季節(jié)氣象要素指標(biāo)。以各站點1971—2010年共30年氣溫與降水量資料，應(yīng)用線性回歸和F檢驗法分別進行趨勢與階段性分析和顯著性檢驗。使用數(shù)學(xué)統(tǒng)計與Thornthwaite Memoriae 模型，在未來氣候預(yù)測結(jié)果的基礎(chǔ)上，分析了冬小麥核心產(chǎn)區(qū)氣候變化對其產(chǎn)量的影響。氣象資料為國家氣象信息中心收集整理的中國地面氣候資料年值數(shù)據(jù)集，干旱數(shù)據(jù)來源于農(nóng)業(yè)部種植業(yè)司“災(zāi)情數(shù)據(jù)庫”和河南統(tǒng)計年鑒，其中l(wèi)967、1968、1969年3年數(shù)據(jù)缺失[21-25]。

2結(jié)果與分析

2.1豫東地區(qū)氣候變化

通過分析（圖1、圖2）可知，豫東地區(qū)20 年氣候變化基本特征為溫度顯著升高，降雨量逐漸增加，日照時數(shù)逐漸減少。歷年季節(jié)降雨量、日照時數(shù)分析結(jié)果表明，夏、秋兩季降雨量增加幅度較冬季大；夏、秋兩季日照時數(shù)降幅明顯，春季日照時數(shù)降幅較弱。

豫東地區(qū)各季節(jié)平均氣溫除秋季外，整體上呈現(xiàn)波動的上升趨勢，從圖3中的線性趨勢可以明顯看出，其增速為冬>秋>春>夏。各季節(jié)平均氣溫均在上世紀(jì)90年代初出現(xiàn)一次增幅變緩甚至為負的波動，隨后增幅變快，春、秋尤為

顯著，這種變化在年平均溫度變化曲線中也有明顯體現(xiàn)。

近20年豫東地區(qū)年降水量整體上呈波動式增加（圖2）。自1991年起5年降水量有一明顯波動，呈現(xiàn)為先減少后增加，在1995年達到最大，隨后驟然下降，并在1997年達到20年的最低值，隨后以2～3年為周期震蕩爬升，其中2003 年達到20年的最高值，隨后一直在低位波動前行緩慢上升。通過分季節(jié)統(tǒng)計可知豫東地區(qū)降水量主要集中在夏季，為全年降水量的66%以上，其次為春和秋兩季分別約占14%左右，冬季降水量最少，僅為全年降水的5%左右。夏季降水與年降水的變化趨勢相關(guān)明顯，冬、春兩季降水逐漸減少，只有秋季

降水量為緩慢增加。通過對近20年降水量增加分析，發(fā)現(xiàn)夏季對年降水量增加的貢獻最大（圖4），說明夏季降水量變化對年降水量的變化影響最為明顯。

2.2豫東地區(qū)氣候變化對冬小麥產(chǎn)量的影響

2.2.1氣象產(chǎn)量計算

式中：Y為糧食單產(chǎn)；Yt 是趨勢產(chǎn)量，為以原始產(chǎn)量資料為基礎(chǔ)進行模擬的關(guān)于時間趨勢的函數(shù)；Yw 是氣象產(chǎn)量，主要受短周期變化因子（氣象因素）影響；e是減產(chǎn)分量，主要受病蟲害、社會因素等隨機因素影響，一般因隨機性太高且比重較小，計算中常被忽略[18-21]；故簡化為：Y=Yt+Yw ；則氣象產(chǎn)量為：Yw=Y-Yt。

2.2.2氣候生產(chǎn)力計算氣候生產(chǎn)力通常是指用氣候條件為主要因子進行估算的農(nóng)業(yè)生產(chǎn)潛力，即在各項氣象要素與環(huán)境因素最為適宜的情況下，單位耕地面積所能達到的最高產(chǎn)量，故常被稱為氣象產(chǎn)量。項目通過比較后，選用李斯（Lieth）方法來計算小麥的氣候生產(chǎn)力，該方法有計算便捷且能準(zhǔn)確表達氣候影響的特點。其根據(jù)全球作物產(chǎn)量和年均氣溫及年降雨量之間的對應(yīng)關(guān)系，進而用實際蒸散量來估算作物產(chǎn)量，即Thornthwaite Memoriae 模型[16-19]：endprint

Pv=30 000[1-e-0.000 969 5（v-20）]。

其中：v=1.05r/1+（1.05r/L）2，L=300+25t+0.05t3。

式中：Pv為氣象產(chǎn)量[kg/（hm2·年）]；v為年均蒸散量（mm）；r為年均降雨量（mm）；L為年均最大蒸散量（mm），為氣溫t（℃）的三次函數(shù)，公式經(jīng)驗系數(shù)為30 000。

近20年來豫東地區(qū)冬小麥產(chǎn)量呈波動上升趨勢，一元回歸分析可知單產(chǎn)增速為198 kg/年以上，最大值為2008年的 6 851 kg/hm2，最低值為1994年的3 035 kg/hm2，從此也能看出單產(chǎn)的增速非常顯著（圖5）；通過分析可知上世紀(jì)90年代，氣象產(chǎn)量與氣候生產(chǎn)力起初步調(diào)較為一致，均為下降，隨后分歧明顯，步調(diào)相反，且變化劇烈；2000年以后，兩者變化均趨于平穩(wěn)，且變化趨勢逐漸吻合（圖6、圖7）。針對氣溫與降水同時變化，以氣溫正負1～3 ℃和降水正負10%～30%的條件兩兩結(jié)合，假設(shè)出49種情況，進而計算冬小麥的氣候生產(chǎn)力，其結(jié)果表明“冷干型”氣候?qū)Χ←溕a(chǎn)最為不利；“暖干型”氣候加劇了水分的不足，造成生產(chǎn)力下降；“冷濕型”氣候，導(dǎo)致冬小麥產(chǎn)量增幅不明顯；“暖濕型”氣候，有利于冬小麥產(chǎn)量的提高；而依據(jù)相關(guān)研究結(jié)果對中國大陸未來氣候變化進行預(yù)測，豫東地區(qū)未來80～90年內(nèi)氣候可能向“暖濕型”變化，冬小麥氣候生產(chǎn)力也有較大幅度提高，對冬小麥產(chǎn)量的提升較為有利[20-24]。

2.3豫東地區(qū)干旱災(zāi)害特征分析

豫東地區(qū)干旱數(shù)據(jù)（受災(zāi)及成災(zāi)面積等）來源于農(nóng)業(yè)部 “災(zāi)情數(shù)據(jù)庫”與河南統(tǒng)計年鑒，其中1967、1968、1969年3年數(shù)據(jù)缺失[21-25]，氣象資料來自國家氣象信息中心，選取該地區(qū)4個具有代表性的商丘、開封等站點自1951年起60年的數(shù)據(jù)。將干旱及氣象數(shù)據(jù)進行相關(guān)統(tǒng)計分析，結(jié)果表明，該區(qū)域農(nóng)業(yè)干旱災(zāi)害呈現(xiàn)明顯的周期性波動變化，輕重干旱災(zāi)情年際交替發(fā)生，干旱受災(zāi)種植面積的年際波動周期為4～14 年，干旱成災(zāi)種植面積的年際波動周期為4～17 年；整體出現(xiàn)干旱受災(zāi)面積呈減小趨勢但成災(zāi)面積呈增大趨勢的變化（圖8）；作物的干旱受災(zāi)率和成災(zāi)率呈顯著相關(guān)，兩者隨時間變化出現(xiàn)的周期性波動特征較為一致（圖9）；作物干旱受災(zāi)率異常指數(shù)與成災(zāi)率異常指數(shù)呈周期性負相關(guān)變化，即負正交替出現(xiàn)，變化周期為3～10年（圖10）；作物干旱災(zāi)害強度指數(shù)即成災(zāi)面積與受災(zāi)面積的比例，呈波動上升的趨勢，周期為3～6 年（圖11）。說明該區(qū)域的資源環(huán)境、氣候變化和社會經(jīng)濟條件均對其干旱災(zāi)害的形成有一定的影響[19-29]。

3結(jié)論

本研究基于多年的氣象觀測數(shù)據(jù)和冬小麥實際產(chǎn)量數(shù)據(jù)，應(yīng)用數(shù)學(xué)統(tǒng)計分析與模型模擬方法，重點研究了豫東地區(qū)的氣候變化趨勢，及該變化趨勢對已發(fā)生和將發(fā)生的冬小麥產(chǎn)量的影響，同時對種植區(qū)歷年的干旱災(zāi)害的特征進行了分析，結(jié)果表明，（1）豫東地區(qū)氣候變暖主要發(fā)生于20世紀(jì)80年代以后，春、冬兩季氣候變暖較為顯著，且氣候趨于“暖濕化”發(fā)展；（2）氣溫、降水量、蒸發(fā)量與極端溫度為影響冬小麥的主要氣候因子，“暖濕型”氣候有利于冬小麥生產(chǎn)力的提高，“冷干型”氣候?qū)Χ←溕a(chǎn)最為不利；未來幾十年內(nèi)氣候可能向“暖濕型”變化，對研究區(qū)糧食作物產(chǎn)量的提升較為有利；（3）近50 年降水是限制冬小麥生長發(fā)育的主要因子，氣候因子匹配效果變差對冬小麥的生長不利；在冬小麥生長的中后期光照、溫度和降水的適宜度與冬小麥產(chǎn)量的相關(guān)性顯著；（4）豫東地區(qū)農(nóng)業(yè)干旱災(zāi)害呈周期性波動，干旱受災(zāi)面積呈減小趨勢但成災(zāi)面積呈增大趨勢；資源環(huán)境、氣候變化和社會經(jīng)濟條件均對干旱災(zāi)害的形成有一定的影響。

參考文獻：

[1]Piao S L，Ciais P，Huang Y，et al. The impacts of climate change on water resources and agriculture in China[J]. Nature，2010，467（2）：43-51.

[2]丁一匯，任國玉，趙宗慈，等. 中國氣候變化的檢測及預(yù)估[J]. 沙漠與綠洲氣象，2007，1（1）：1-10.

[3]武永利，盧淑賢. 近45 年山西省氣候生產(chǎn)潛力時空變化特征分析[J]. 生態(tài)環(huán)境學(xué)報，2009，18（2）：567-571.

[4]姚玉壁，李耀輝，王毅榮，等. 黃土高原氣候與氣候生產(chǎn)力對全球氣候變化的響應(yīng)[J]. 干旱地區(qū)農(nóng)業(yè)研究，2005，23（2）：202-208.

[5]高素華，潘亞茹，郭建平.氣候變化對植物氣候生產(chǎn)力的影響[J]. 氣象，1994，20（1）：30-33.

[6]黨安榮，閻守邕，吳宏歧，等. 基于GIS 的中國土地生產(chǎn)潛力研究[J]. 生態(tài)學(xué)報，2000，20（6）：910-915.

[7]謝俊奇，蔡玉梅，鄭振源，等. 基于改進的農(nóng)業(yè)生態(tài)區(qū)法的中國耕地糧食生產(chǎn)潛力評價[J]. 中國土地科學(xué)，2004，18（4）：31-37.

[8]封志明，楊艷昭，張晶，等. 從柵格到縣域：中國糧食生產(chǎn)的資源潛力區(qū)域差異分析[J]. 自然資源學(xué)報，2007，22（5）：747-755.

[9]劉德祥，董安祥，薛萬孝，等. 氣候變暖對甘肅農(nóng)業(yè)的影響[J]. 地理科學(xué)進展，2005，22（2）：49-58.

[10]張強，鄧振鏞，趙映東，等. 全球氣候變化對我國西北地區(qū)農(nóng)業(yè)的影響[J]. 生態(tài)學(xué)報，2008，11（3）：1210-1218.

[11]蘇占勝，陳曉光，黃峰，等. 寧夏農(nóng)牧交錯區(qū)（鹽池）草地生產(chǎn)力對氣候變化的響應(yīng)[J]. 中國沙漠，2007，7（3）：430-435.

[12]郭建平，高素華，潘亞茹. 東北地區(qū)農(nóng)業(yè)氣候生產(chǎn)潛力及其開發(fā)利用對策[J]. 氣象，1995，21（2）：3-9.endprint

[13]羅永忠，成自勇，郭小芹. 近40年甘肅省氣候生產(chǎn)潛力時空變化特征[J]. 生態(tài)學(xué)報，2011，31（1）：221-229.

[14]趙俊芳，郭建平，鄔定榮，等. 2011—2050 年黃淮海冬小麥，夏玉米氣候生產(chǎn)潛力評價[J]. 應(yīng)用生態(tài)學(xué)報，2011，22（12）：3189-3195.

[15]Laux P，Jackel G，Tingem R M，et al. Impact of climate change on agricultural productivity under rainfed conditions in Cameroon：Amethod to improve attainable crop yields by planting date adaptations[J]. Agricultural and Forest Meteorology，2010，150：1258-1271.

[16]Condon A G. Breeding for high water-use efficiency[J]. Journal of Experimental Botany，2003，54（1）：20.

[17]莫興國，劉蘇峽，林忠輝，等. 華北平原蒸散和GPP格局及其對氣候波動的響應(yīng)[J]. 地理學(xué)報，2011，66（5）：589-598.

[18]Lobell D B，Sibley A，Ivan O M J. Extreme heat effects on wheat senescence in India[J]. Nature Climate Change，2012，2（3）：186-189.

[19]Pongratz J，Lobell D B，Cao L，et al. Crop yields in a geoengineered climate[J]. Nature Climate Change，2012，2（2）：101-105.

[20]余衛(wèi)東，黃玉超，康紹鈞. 商丘市極端天氣氣候事件變化趨勢[J]. 氣象科技，2005（33）：127-130.

[21]IPCC. Climate change 2007：Synthesis report.Contribution of working groups Ⅰ，Ⅱ and Ⅲ to the fourth assessment report of the intergovernmental panel on climate change[R].Genena：IPCC，2007.

[22]Root T L，Price J T，Hall K R，et al. Fingerprints of global warming on wild animals and plants[J]. Nature，2003，421（6918）：57-60.

[23]Thomas C D，Cameron A，Green R E，et al. Extinction risk from climate change[J]. Nature，2004，427（6970）：145-148.

[24]丁一匯，任國玉，石廣玉. 氣候變化國家評估報告（Ⅰ）：中國氣候變化的歷史和未來趨勢[J].氣候變化研究進展，2006，2（1）：3-8.

[25]任國玉，初子瑩，周雅清，等. 中國氣溫變化研究最新進展[J]. 氣候與環(huán)境研究，2005，10（4）：701-716.

[26]范澤孟，岳天祥，陳傳法，等. 中國氣溫與降水的時空變化趨勢分析[J]. 地理信息科學(xué)，2011，13（4）：526-533.

[27]王紹武，葉瑾琳，龔道溢，等. 近百年中國年氣溫序列的建立[J]. 應(yīng)用氣象學(xué)報，1998，9（4）：392-401.

[28]郭志梅，廖啟龍，李雄. 中國北方地區(qū)近50 年來氣溫變化趨勢及其突變性[J]. 干旱區(qū)地理，2005，28（2）：176-182.

[29]FAO/IIASA/ISRIC/ISS-CAS/JRC. Harmonized world soil database（version 1.1）[R].Rome：FAO and IIASA，2009.趙孟良，王麗慧，孫雪梅，等. 干旱脅迫下菊芋可溶性碳水化合物的積累及分配規(guī)律[J]. 江蘇農(nóng)業(yè)科學(xué)，2015，43（2）：340-343.endprint