近20年河湟谷地糧食作物種植結構對氣候變化的響應

陳有祿，周 強，陳 瓊，毛旭峰，劉峰貴

(青海師范大學生命與地理科學學院，青海 西寧 810008)

近百年來由于溫室效應的影響，全球平均氣溫升高了0.5℃[1,2]，特別是北半球高緯度地區氣候變暖比較明顯[3]。因而影響到區域農業生產結構和全球糧食安全，已成為全球變化研究的熱點[4]。研究表明，氣候變暖對中國農業的產量、熟制種植結構都產生一定影響[5-6]。未來氣候繼續變暖將使我國農作物的種植結構發生巨大改變，據初步估算，2030年我國種植業產量將因全球變暖而減少5%～10%，其中小麥、水稻和玉米三大農作物都可能會減產[7]，西北地區的研究也表明溫度升高對該地區農業種植結構和產量產生了重大影響[8-12]。熱量資源增加使喜溫農作物的范圍向高海拔地區擴展，在水資源條件適當的情況下對喜溫農作物的種植非常有益。河湟谷地地處青藏高原和黃土高原的過渡地帶，屬于干旱、半干旱氣候，全球氣候變化背景下該地區氣候如何變化？該地區糧食作物的種植結構對氣候變化如何響應？本文試圖通過對河湟谷地氣候變化的實際和糧食作物種植結構的變化來探討糧食作物種植結構對氣候變化的響應關系，以期為地方農業生產和規劃提供參考。

1 研究區概況

河湟谷地位于青藏高原東北部，區域范圍包括湟水河流域及黃河流域龍羊峽至寺溝峽一段，行政單元包括青海省的西寧市(含城西、城北、城中、城東、大通、湟中、湟源等4區3縣)、海東市(含平安、樂都、民和、互助、化隆、循化等2區4縣)以及海南藏族自治州的貴德縣和黃南藏族自治州的同仁縣、尖扎縣，共6區10縣，總面積約3.3×104km2。

河湟谷地系湟水河、黃河及其支流沖積形成的河谷地區[13]，地勢相對于青海其它地區海拔較低，氣候溫和，是青海省重要的農業生產基地，但是由于溫度條件的限制，河湟谷地農業生產熟制為一年一熟，單位產量較低，耕地面積有限，隨著近20年來氣候變化和生產技術進步，部分地區由原來的一年一熟逐漸演變為二年三熟[14，15]，且這種種植區域有擴大的趨勢，尤其是氣候變化使河湟谷地農業種植結構發生顯著變化。

圖1 研究區位置Fig.1 Location of the study area

2 數據來源與分析方法

2.1 數據來源

2.1.1 氣象數據 本文選取1991～2010年河湟谷地12個氣候站點氣象數據主要來自中國氣象數據共享服務網(http://data.cma.cn)。

2.1.2 糧食作物數據 1990～2012年各縣、區的糧食作物數據來自《青海省農業統計年鑒》[16]和《青海省經濟統計年鑒》[17]。

2.1.3 問卷數據 樣點主要均勻的分布在民和、互助、平安、大通、湟中、湟源、海晏一區六縣的25個村，總共發放825份問卷，收回800份問卷，有效率達97%。

2.2 分析方法

2.2.1 數據分析 實測經驗和理論推算均認為降水量在350～500 mm的地區發展旱作農業比較可靠，“有收無收在于水”是旱作農業的普遍共識，河湟谷地地區降水量在400 mm左右，是適合旱作農業的地區，同時春、夏季節降水量接近300 mm，占全年降水量的3/4左右，雨熱同季，有利于農作物生長。但與其它旱作農業區相比較，該地區屬于青藏高原與黃土高原的過渡區域，海拔在1 650 m到3 500 m之間，農作物生長不僅受到水分因素的影響，同時也受因海拔導致的溫度變化而產生的影響，其影響因素具有雙重性，部分地區溫度因素的限制可能大于水分因素。

氣溫和海拔高度影響下的積溫是農作物生長的重要物理基礎，本文選擇1991～2010年積溫變化作為本地區重要的限制性因子。河湟谷地近年氣溫上升比較顯著，秋季和冬季的增溫變化比較顯著，20世紀90年代比60年代和80年代的平均氣溫增高了0.3～0.5℃[18]。河湟谷地川水地區收獲后，≥0℃的余熱資源為1 068～1 613℃·d,占全年的39%～47%；淺山地區作物收獲后，≥0℃的余熱資源為496～622℃·d，占全年的35%～40%[14]。收獲后余熱資源的總量較大，適合于農作物的生長，如冬小麥，蔬菜等，為河湟谷地余熱資源的有效利用和種植結構的調整起到了鋪墊作用。

河湟谷地糧食作物的種植結構主要有小麥、薯類、青稞、油菜和蠶豆等經濟作物，熟制為一年一熟。近20年來河湟谷地農作物總播種面積呈現增加態勢，其中經濟作物的播種面積增加趨勢明顯，糧食作物的種植面積明顯下降，冬小麥和喜溫作物玉米的種植是一個新的發現[19]，并且近20年來種植面積不斷擴大。因此，選擇小麥、薯類、玉米3種主要糧食作物的播種面積變化來反映農作物種植結構變化。一般作物種植結構的變化可用作物各自所占總播種面積的百分比(A)或者作物各自占主要糧食作物的播種面積(糧食作物總播種面積之和)的百分比(B)進行判斷，前者百分比變化往往受經濟作物種植面積變化的影響較大，無法準確判斷糧食作物種植結構的變化是經濟作物增加導致的還是氣候變化導致的，影響分析判斷的結果；后者則是在糧食作物播種總面積的內部進行種植結構比較，總面積的變化不會影響其內部種植結構的比較結果，能夠真實反映糧食作物種植結構的變化判斷。因此，本文采用后者(B)對河湟谷地糧食作物的種植結構進行判斷。

2.2.2 快速聚類分析及數據處理 河湟谷地糧食作物的種植結構雖然只有三大類,但近20年的不同地區的組合類型相對復雜，屬于大樣本資料,因此需要將數據對象進行深度挖掘，按照其相似性歸成若干類別。聚類分析是數據深度挖掘常用的技術和方法，其目的是可以將數據對象分成若干個類，使同一個類中的對象高度相似，而不同類中的對象差異最大[20]。

快速聚類法是采用歐氏距離法進行計算，其原則是每個樣本點距離本組的矩心(本組數據平均值所構成的點)的最近距離來進行分類。聚類過程中必須使得組內的因子和預測量的方差變小，而使組間的距離增大。一般先給出一定的判別值，作為變動范圍的限制值，其公式可寫為[21]：

(1)

其中，k表示每個樣本中有k個變量；xi表示第一個樣本在第i個變量上的取值；yi表示第二個樣本在第i個變量上的取值。

(1)初定組別。

(2)組的分割。計算4組內每一個因子的標準差及在該組中預報量的凝聚直徑記為因子的標準差。

CK=maxyi-minyi(i∈nk;k=1，2，…，G)

將4組內每一個因子的標準差和每一組的極值與事先給定的閾值δ及θ進行比較，

若在3個因子中，每一個因子的標準差最大者大于δ或Ck>θ，則將該組以α值分裂為兩組。確定α值的辦法：若每一組的極值≤θ，則取第k組預報量y的中值，否則取標準差最大的那個因子的平均值作為分組界限。

(3)計算組間的距離，重復(2)-(3)步，直到組無法分割或合并為止[22]。

在實際操作中，我們將1990、1995、2003、2006、2010年各縣的小麥、薯類、玉米種植面積占3種作物總播種面積的百分比用快速聚類方法結合河湟谷地的實際情況我們把河湟谷地過去20年的種植結構分為4類，即以小麥為主、以小麥和薯類為主、以薯類為主、以薯類、小麥和玉米多元發展趨勢為主(表1)，在此過程中共有840個數據參加了分組，從而獲得了每種種植結構類型的空間分布及其多年來的變化。

表1 河湟谷地糧食作物種植結構類型及其聚類中心

3 結果分析

通過實測氣象數據分析(圖2)和800份問卷的訪談調查證明(圖3)，河湟地區近20年來降水穩定，水分條件基本能夠保證大部分區域農作物生長的需求，熱量條件的變化可能是造成河湟谷地糧食作物種植結構差異的主要原因。器測氣溫變化和農戶對氣溫變化的感知基本一致。氣象數據顯示近20年降水變化不大，一直處在350～400 mm左右，適宜于農作物的生長，但在調查中67.75%的農戶認為氣候比較干燥，只有24.88%的人認為氣候比較濕潤，與實際降水相差較大，存在差異的主要原因可能是河湟谷地降水趨于集中，降水變率增大，頻度減小，春夏季節降水集中，秋冬季降水減少。

圖2 近20年河湟谷地年均溫度和降水變化Fig.2 Nearly 20 years in Yellow River-Huangshui River Valley(YHV) annual average temperature and precipitation changes

圖3 河湟谷地農戶對氣候變化感知Fig.3 Farmers perceive the intensity of climate change in YHV

本文將1990、1995、2003、2006、2010年作為典型年份，其中1990年代表氣候變暖開始的第一個轉折點[23]，1995年代表氣候顯著變暖的年份，2003、2006、2010年分別代表河湟谷地玉米、新經濟作物等新種農作物的種植及擴大年份。然后根據河湟谷地各縣、市、區小麥、薯類、玉米3種糧食作物的播種面積占其總播種面積的百分比劃分不同的種植類型，再根據其變化情況劃分為不同的結構類型。快速聚類過程中運用5年的數據得出河湟谷地過去20年的4類種植結構類型，1990和1995年河湟谷地種植結構以小麥為主，2003年河湟谷地種植結構演變為小麥和薯類為主，2006年河湟谷地種植結構演變為薯類為主，到2010年后河湟谷地種植結構演變成以薯類、小麥和玉米等多元化農作物的種植結構。

1990年以來，河湟谷地薯類的種植面積隨著溫度的升高不斷地擴大[24]。種植結構的更大變化主要發生在2003年左右，積溫的變化可能是導致種植結構變化的主要原因，2000年以后貴德縣、樂都區、民和回族土族自治縣、尖扎縣、循化撒拉族自治縣≥10℃的積溫均已達到≥2 700℃，部分地區當時就出現了喜溫農作物玉米的種植(表2、圖5、圖6)，因為玉米全生育期所需要的≥10℃的積溫為≥2 700℃[25,26]。由表2和當地實際情況可知河湟谷地≥10℃積溫≥2 700℃的地區都有玉米的種植，進一步說明積溫的升高才使得玉米種植范圍不斷地擴大。

表2 青藏高原≥10℃有效積溫區劃與相關作物系統[26]

3.1 糧食作物播種面積變化

近20年來河湟谷地主要糧食作物的總播種面積呈現減少趨勢(圖4)，糧食作物內部種植結構的面積比重也呈現出不同的變化趨勢(圖5)。

圖4 近20年河湟谷地主要糧食作物播種面積的變化Fig.4 Changes of main grain crops planting areas in nearly 20 years in YHV

注：(a)1995年積溫等值線;(b)1995年薯類比重;(c)1995年小麥比重;(d)2003年玉米比重;(e)2010年積溫等值線;(f)2010年薯類比重;(g)2010年小麥比;(h)2010年玉米比重;(o)1995～2003年積溫變化;(p)1995～2010年積溫變化;(q)1990年;(r)1995年Note: (a)accumulated temperature contour line in 1995;(b)potato proportion in 1995;(c)wheat proportion in 1995;(d)cornproportion in 1995;(e) accumulated temperature contour line in 2010;(f)potato proportion in 2010;(g)wheat proportion in 2010;(h)cornproportion in 2010;(o)variation of accumulated temperature between 1995 to 2003;(p)variation of accumulated temperature between 1995 to 2010;(q)1990 year;(r)1995 year圖5 1995和2010年代河湟谷地10℃積溫均值變化與主要糧食作物播種面積百分比變化Fig.5 Changes of the mean accumulative temperature and the planting percentage in the 1995s and 2010s in YHV

20世紀90年代前半葉，河湟谷地糧食作物種植結構主要以小麥為主，播種面積約為123.5×104hm2，1992年開始下降，至2010年減少至52.1×104hm2，播種面積在3種農作物中的比例從80%下降到37.4%；薯類種植面積在2000年之前一直徘徊在33×104hm2左右，約占當時3種糧食作物總播種面積的19%左右，2000年以后播種面積呈現上升趨勢，至2010年播種面積增加到78.9×104hm2，約占3種糧食作物播種面積的56.7%；2001年新的變化是本地區開始出現喜溫作物玉米的種植[27]，播種面積一直呈現增加趨勢，糧食作物種植結構由原來的以小麥為主逐漸演變為以薯類、小麥和玉米多元發展趨勢的種植結構。

3.2 糧食作物種植格局的變化

20世紀90年代，河湟谷地小麥種植區域主要分布在湟水中游的西寧、大通、湟中與民和等市、縣(圖5)，種植結構比重分別為11.81%(西寧)、10.32%(大通)、14.60%(湟中)、5.36%(民和)，至2010年，該地區小麥種植結構比重降低到0.54%(西寧)、6.05%(大通)、6.06%(湟中)、3.04%(民和)，遠小于上世紀90年代。

薯類和玉米的種植結構比重則呈現大幅增加的趨勢(薯類播種面積所占的比重由1995年的18.96%增加到現在的56.60%；玉米播種面積由2003年的0.68%增加到了現在的6.03%)，其空間分布格局也呈現出顯著變化的趨勢，20 世紀90年代，薯類作物主要分布在大通河流域的門源縣，至2010年，薯類種植范圍擴大到河湟谷地的互助、民和、樂都和大通等區縣，2003年玉米種植出現在河湟河下游氣候較好地區，其種植面積和種植結構比重不斷增加(湟水河下游民和回族土族自治縣玉米的種植面積由2003年的0.64%增加到2010年的5.8%，湟水河中游互助縣2010年開始種植玉米，黃河上游的貴德縣玉米種植面積由2003年的0.002%增加到了2010年的0.09%)。樂都區、尖扎縣、循化撒拉族自治縣等從2001年出現了玉米的種植(圖5、圖6)。

注：(a)2003年;(b)2006年;(c)2010年。圖6 10℃積溫差值和不同時期河湟谷地主要糧食作物種植格局Fig.6 10 ℃ accumulated temperature difference with main grain crops structure in different periods in YHV

3.3 糧食作物種植格局同溫度變化間的關系

河湟谷地主要糧食作物的種植結構伴隨著該地區積溫的變化而變化，該區域的北部、西南和東南地區積溫變化明顯和實際農作物種植結構變化相吻合。由(圖6)可知積溫變化明顯的地區種植結構變化較大，特別是西北地區、西南地區、東南地區變化比較明顯。1990年到1995年積溫基本沒有發生變化，同期小麥和薯類的種植面積基本沒有發生變化。到了2003年，河湟谷地的東部積溫較高基本都≥2 500℃，同期河湟谷地的東部地區較1990年糧食作物種植結構的變化非常明顯，由原來的以小麥為主的種植結構變成了21世紀初的以小麥和薯類為主的種植結構和以薯類為主的種植結構，說明積溫的升高促進了該地區的種植結構比重的變化。已有相關研究證明，小麥和薯類屬于喜涼作物，但是薯類生長所需要的溫度比小麥的高[28]，因此隨著積溫的增加薯類的種植結構比重也會呈現增加的趨勢[24]。

從整個河湟谷地來看，其糧食作物的種植結構的變化和積溫的變化存在一定程度的關系，積溫的差值圖變化明顯的地區糧食作物的種植結構也有所變化，特別是河湟谷地北部、西南和東南地區積溫變化較明顯。其中，貴德縣、尖扎縣、民和縣、樂都區和門源縣的種植結構由原來的以小麥為主的種植結構變成了現在以薯類和玉米為主的種植結構。

4 討論和結論

4.1 討 論

本文在論述河湟谷地糧食作物種植結構變化和氣候變化之間的關系時，借鑒了云雅如，方修琦(2005)假定其它氣候因子不變的方法。假定近20年來河湟谷地土壤肥力、地形、海拔因素不變，考慮氣溫和降水因素的方法，相關數據表明河湟谷地近20年來降水一直保持在一個相對穩定的水平，因此重點考慮溫度因素的影響，較真實的分析出河湟谷地近20年來種植結構和種植格局與溫度的關系。

表3 河湟谷地過去20年主要糧食作物種植結構的變化結果

由于河湟谷地處于特殊的地理位置，故將實際氣候數據和人們的感知數據相結合，有助于更加明確地了解河湟谷地氣候的變化。雖然本文和云雅如，方修琦(2005)的思想基本一致，但是在文章中也存在較大的差異。在處理數據的過程中，本文總結出了青藏高原≥10℃有效積溫區劃與相關作物系統之間的關系和河湟谷地過去20年主要糧食作物種植結構的變化，為河湟谷地農業發展提供參考依據。除了氣溫對河湟谷地種植結構的影響外，政府政策因素，農業技術還有人類活動對糧食作物的種植格局都會影響糧食作物的種植結構。在此由于特殊地區的條件限制故沒有加入這些影響因素的分析。

4.2 結 論

1)過去20年，在河湟谷地糧食作物種植結構比重方面，小麥的種植結構比重呈現下降的趨勢，薯類和玉米的種植比重明顯增加，該地區的種植結構由原來的以小麥為主的種植結構變成了21世紀初的以小麥、薯類和玉米不同組合的種植結構方式。

2)在種植格局方面，薯類的種植格局向北部和東部地區移動，小麥的種植格局在整個河湟谷地都在大范圍減少，玉米的種植格局向西南和東南地區移動。

3)河湟谷地主要糧食作物的種植格局和溫度之間存在一定程度的關系。過去20年里，隨著積溫的增加，喜溫作物的種植面積大幅度增加，溫度變化越明顯的地區糧食作物的種植結構的變化也越明顯。整個河湟谷地的喜溫農作物的種植界線隨著溫度的增加向西南和東南移動并呈現增加的趨勢。