陈有禄，周 强，陈 琼，毛旭峰，刘峰贵

(青海师范大学生命与地理科学学院，青海 西宁 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年里，随着积温的增加，喜温作物的种植面积大幅度增加，温度变化越明显的地区粮食作物的种植结构的变化也越明显。整个河湟谷地的喜温农作物的种植界线随着温度的增加向西南和东南移动并呈现增加的趋势。