氣候變化對豫東冬小麥產區(qū)的綜合影響

張怡　史本林

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

關鍵詞：氣候變化；冬小麥；豫東；極端天氣；干旱；產量

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

收稿日期：2014

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

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

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

1材料與方法

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

2結果與分析

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

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

豫東地區(qū)各季節(jié)平均氣溫除秋季外，整體上呈現(xiàn)波動的上升趨勢，從圖3中的線性趨勢可以明顯看出，其增速為冬>秋>春>夏。各季節(jié)平均氣溫均在上世紀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%左右。夏季降水與年降水的變化趨勢相關明顯，冬、春兩季降水逐漸減少，只有秋季

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

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

2.2.1氣象產量計算

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

2.2.2氣候生產力計算氣候生產力通常是指用氣候條件為主要因子進行估算的農業(yè)生產潛力，即在各項氣象要素與環(huán)境因素最為適宜的情況下，單位耕地面積所能達到的最高產量，故常被稱為氣象產量。項目通過比較后，選用李斯（Lieth）方法來計算小麥的氣候生產力，該方法有計算便捷且能準確表達氣候影響的特點。其根據全球作物產量和年均氣溫及年降雨量之間的對應關系，進而用實際蒸散量來估算作物產量，即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為氣象產量[kg/（hm2·年）]；v為年均蒸散量（mm）；r為年均降雨量（mm）；L為年均最大蒸散量（mm），為氣溫t（℃）的三次函數(shù)，公式經驗系數(shù)為30 000。

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

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

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

3結論

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

參考文獻：

[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]丁一匯，任國玉，趙宗慈，等. 中國氣候變化的檢測及預估[J]. 沙漠與綠洲氣象，2007，1（1）：1-10.

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

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

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

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

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

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

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

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

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

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

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

[14]趙俊芳，郭建平，鄔定榮，等. 2011—2050 年黃淮海冬小麥，夏玉米氣候生產潛力評價[J]. 應用生態(tài)學報，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格局及其對氣候波動的響應[J]. 地理學報，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]. 地理信息科學，2011，13（4）：526-533.

[27]王紹武，葉瑾琳，龔道溢，等. 近百年中國年氣溫序列的建立[J]. 應用氣象學報，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]. 江蘇農業(yè)科學，2015，43（2）：340-343.endprint