魯永新 ，鄒 萍 ，張中平，武 勇

（1.楚雄州氣象局，云南 楚雄 675000；2.楚雄州農(nóng)業(yè)技術(shù)推廣研究所，云南 楚雄 675000）

小麥?zhǔn)窃颇鲜〉闹饕杭Z食作物之一，常年種植46.67萬hm2左右，最高年達(dá)77.33萬hm2，主要種植區(qū)分布在云南的紅河、曲靖、楚雄、大理、臨滄、保山、玉溪、麗江、昆明等地，秋季種植，春季收獲，全生育期時段正好與云南季風(fēng)氣候的干季相對應(yīng)[1-2]。因受云南的地理、氣候、經(jīng)濟(jì)等因素的影響，絕大部分山區(qū)、半山區(qū)水利基礎(chǔ)設(shè)施建設(shè)滯后或不能滿足干季農(nóng)業(yè)生產(chǎn)用水的基本需求，小麥生長期內(nèi)基本無水利灌溉條件，主要在自然氣候條件下生長。隨著城鎮(zhèn)化的發(fā)展，云南耕地垂直空間還將向山區(qū)、半山區(qū)坡地轉(zhuǎn)移[3]，也就是說隨著時間的推移小麥種植將向廣大山區(qū)、半山區(qū)推移，小麥產(chǎn)業(yè)發(fā)展?jié)摿薮?，干季氣候因素的對小麥生產(chǎn)的影響權(quán)重也將會隨之增加。因此，氣候因素影響將是云南小麥產(chǎn)量和質(zhì)量的關(guān)鍵性和決定性因子。

楚雄是最具云南季風(fēng)氣候影響區(qū)典型的半干旱區(qū)的代表區(qū)域之一。全區(qū)域內(nèi)多年平均降雨量670～1 100 mm，其中干季（11月至次年4月）降雨80～130mm。氣候因素、農(nóng)業(yè)水利基礎(chǔ)建設(shè)及農(nóng)業(yè)灌溉條件在云南全省的廣大山區(qū)、半山區(qū)具有很好的代表性。以楚雄州小麥種植區(qū)為代表進(jìn)行了氣候生態(tài)因素年際間的動態(tài)變化對小麥種植影響的綜合評價試驗研究，為綜合評價云南半干旱季風(fēng)氣候區(qū)小麥?zhǔn)粘傻哪觌H變化提供依據(jù)。

1 材料和方法

1.1 研究試驗資料來源

氣候資料來源于楚雄州姚安、楚雄、大姚、祿豐、武定、雙柏等10個縣市1961～2011年50 a間的地面氣象觀測資料；小麥產(chǎn)量（單產(chǎn)）資料1985年以后來源于楚雄州統(tǒng)計局“楚雄州大小春糧食產(chǎn)量”統(tǒng)計資料，1985年以前來源于楚雄州農(nóng)業(yè)技術(shù)推廣研究所“國家現(xiàn)代小麥產(chǎn)業(yè)技術(shù)體系楚雄州綜合試驗站”歷史統(tǒng)計資料。

1.2 小麥種植區(qū)氣候生態(tài)因子選取

云南絕大多數(shù)小麥種植區(qū)一般在10中旬至11月上旬播種，12月中旬至次年1月上旬分蘗結(jié)束，次年1月至2月拔節(jié)孕穗，3月抽穗開花，4月至5月上旬灌漿成熟，5月中旬收獲前收獲結(jié)束，全生育期約165～200 d[1]。全生育期間的降雨量、氣溫、日照時數(shù)的時間和空間分布是影響云南地麥生長發(fā)育的關(guān)鍵因子[4]。因此，以楚雄州小麥種植區(qū)為研究典型對象，選擇了6個氣象影響因子作為氣候指標(biāo)[1-5]，分別為x1（11月至次年4月總?cè)照諘r數(shù)）；x2（10月降雨量）；x3（12月中旬次年2月時段總降雨量）；x4（3月下旬至4月上旬初時段平均氣溫）；x5（3月中旬后期至4月上旬末時段平均氣溫）；x6（3月上旬至4月上旬時段總降雨量）。

1.3 綜合評價方法和原理

利用多元統(tǒng)計學(xué)的主成分分析方法建立評價模型，根據(jù)評價模型得出的綜合評價值來判別動態(tài)氣候生態(tài)條件影響下的小麥產(chǎn)量。主成分分析原理是利用降維的思想，把眾多變量轉(zhuǎn)化為少數(shù)幾個綜合指標(biāo)，這幾個綜合指標(biāo)為原來的線性組合；綜合指標(biāo)保留了原始變量的主要信息，彼此之間又不相關(guān)，能使復(fù)雜的問題簡單化，使原因子的信息集中到幾個主要因子上進(jìn)行分析[6-8]。主成分分析方法的步驟是：（1）對原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，以消除量綱不同的影響；（2）求無量綱后數(shù)據(jù)的相關(guān)矩陣R；（3）求R的特征值、特征向量及貢獻(xiàn)率；（4）根據(jù)貢獻(xiàn)率確定主成分的個數(shù)；（5）以線性組合中權(quán)重較大的幾個影響因子作綜合意義解釋，求主成分值，并計算綜合得分值，對研究對象進(jìn)行評價。

1.4 評價樣本的確定

根據(jù)小麥種植研究對象區(qū)（楚雄）秋播春夏收小麥50年單產(chǎn)變化作為評價體系樣本，以小麥生長期內(nèi)降雨量、氣溫、日照資源的時段內(nèi)分布及年際間的動態(tài)變化，考察變量影響的內(nèi)在關(guān)系，將小麥種植代表地區(qū)1961～2011年連續(xù)50個小麥種植年的6個氣候影響指標(biāo)因子值組成一個P＝6、N＝50的矩陣Xpn；前m個主成分組成的矩陣Ymn反映了“集中”原矩陣Xpn中的因子信息[8]。

1.5 數(shù)據(jù)處理結(jié)果及解釋

利用SPSS軟件，對矩陣X6×50進(jìn)行主成分分析。用方差最大法進(jìn)行正交旋轉(zhuǎn)，使成分負(fù)載向0和1兩極分化。計算主成分分值采用回歸法。

2 結(jié)果與分析

2.1 試驗數(shù)據(jù)分析

經(jīng)數(shù)據(jù)分析選出了3個主成分（m=3），處理結(jié)果見表1，樣品相關(guān)矩陣的特征向量如表2所示。

表1 X pn的協(xié)方差陣特征值及方差和累積方差貢獻(xiàn)率

表2 樣本相關(guān)矩陣的特征向量（V矩陣）

2.2 主成分分析結(jié)果表達(dá)式

主成分分析結(jié)果表達(dá)式如下：

2.3 主成分的意義及相關(guān)解釋

表1中的荷載率絕對值的大小表明了各因子與主成分的相關(guān)程度，荷載率越大的因子，與其主成分的相關(guān)性越好，越能顯現(xiàn)出該主成分的影響效應(yīng)。一般當(dāng)|荷載率|≥0.6時，認(rèn)為該要素與主成分的相關(guān)程度足夠大[9-10]。在主成分Y6×50中的前3個向量的累積方差已達(dá)83.094%（見表1），因此只需計算3個主成分。

從表2中的荷載率可看出：第1主成分中，荷載率較大的變量依次是X4，X5，X6，這3個因子反映了小麥生長發(fā)育中后期熱量與水分。其中X4代表小麥揚(yáng)花期平均氣溫，X5代表小麥灌漿乳熟期的平均氣溫，指地麥生育期熱量因素的貢獻(xiàn)較大，X6代表小麥抽穗揚(yáng)花期的降雨量，從變量因子荷載率符號可以看出抽穗揚(yáng)花至乳熟期低溫陰雨也會產(chǎn)生負(fù)影響，與實(shí)際相符；第2個主成分中的荷載率較大的變量依次是X3，X1，這2個因子反映了小麥分蘗至拔節(jié)期降雨量和出苗至乳熟期光照強(qiáng)度，即水熱條件，其中，X3的因子荷載率為正，與X1的因子荷載率符號相反，小麥分蘗至拔節(jié)期降雨量（地麥的水分敏感期降雨量）呈正影響效應(yīng)，對小麥后期產(chǎn)量形成至關(guān)重要，該時段光照過強(qiáng)可能與高溫有關(guān)，會存在前期徒長或后期逼熟現(xiàn)象，對小麥生長及產(chǎn)量形成產(chǎn)生負(fù)效應(yīng)，與地麥生產(chǎn)實(shí)踐十分吻合；第3個主成分中荷載率較大的變量依次是X2和X9，這2個因子主要反映了小麥生長發(fā)育期的水分供給及影響，即水分供給狀況，其中X2為小麥播種前后的時段降雨量，因子荷載率為正，說明此時段降雨量對小麥出苗產(chǎn)生正效應(yīng)，X6為小麥抽穗揚(yáng)花期陰雨天氣，因子荷載率為負(fù)，說明該時段陰天宜少，降雨量適中為好，若陰天多會在絕大多數(shù)山區(qū)發(fā)生低溫，對產(chǎn)量形成產(chǎn)生負(fù)效應(yīng)，與生產(chǎn)實(shí)際相吻合。

2.4 建立綜合評價模型

以各主成分貢獻(xiàn)率為權(quán)重（數(shù)），構(gòu)建綜合主成分模型，即綜合評價模型如下：

圖1 1962～2011年楚雄小麥種植區(qū)綜合評價得分變化

根據(jù)綜合評價模型計算年際間得分值，評價50 a來小麥年際間產(chǎn)量動態(tài)變化。模型計算出的綜合評價結(jié)論中，幾個氣候因素影響突出年的綜合評價結(jié)果與實(shí)際非常吻合（圖1），1974年，1986年，1989年2004年和2010年是楚雄州歷史上罕見的干旱，或干旱低溫冷害影響特別突出的幾年，這幾年的小麥產(chǎn)量減產(chǎn)幅度均大于3成；模型計算結(jié)論也驗證了氣候因素影響突出的幾個小麥豐產(chǎn)年，分別是 1983、1987、1993、2001、2002、2006、2007、2008、2011 年。

3 小結(jié)

綜合評價模型計算得到的楚雄州小麥主要種植區(qū)1962～2011年50 a間的綜合評價得分值增減變化與50 a間小麥實(shí)際產(chǎn)量的動態(tài)增減變化具有很好的對應(yīng)關(guān)系，評價結(jié)果及氣候生態(tài)因子對小麥生產(chǎn)的影響效應(yīng)與多年實(shí)踐相符。研究試驗中采用主成分分析方法以楚雄州氣候因素的年際間變化對小麥產(chǎn)量的動態(tài)影響作用綜合評價結(jié)果，排除了人為因素干擾，避免了通常只考慮較少影響因子指標(biāo)進(jìn)行評價而造成的評價試驗結(jié)果片面性。綜合評價模型研究結(jié)果可操作性強(qiáng)，容易推廣應(yīng)用于工作實(shí)踐中。

[1] 魯永新，楊永生.楚雄州農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與氣候[M].云南民族出版社，2004.77-79.

[2] 王 宇.云南山地氣候[M].昆明：云南科技出版社，2006.281-292.

[3] Tong Shao Yu,Lu Yong Xin.A Study on vertical spatial transfer of cultivated land in Yunnan Province from 1986 to 2008[J].Agricultural science&technology,2011,12（1）:127-132.

[4] 鄒 萍，武 勇.氣候條件對旱地小麥生產(chǎn)影響及增產(chǎn)措施[C].楚雄州科協(xié)論文集[A].2010.10-12.

[5] 馮秀藻，陶炳炎.農(nóng)業(yè)氣象學(xué)原理[M].北京，氣象出版社，1991.311-322.

[6] 袁志發(fā)，周靜芋.多元統(tǒng)計分析[M].北京：科學(xué)出版社，2003.188-206.

[7] 魯永新，張映翠，王恩超，等.楚雄州烤煙種植的生態(tài)區(qū)劃方法研究[J].中國農(nóng)業(yè)氣象，2009，30（S1）：93-96.

[8] 施 能.氣象科研與預(yù)報中的多元分析方法[M].北京：氣象出版社，2001.75-100.

[9] 魏占雄，謝衛(wèi)東.三江源興海地區(qū)西北針茅產(chǎn)量影響因子的主成分分析[J].安徽農(nóng)業(yè)科學(xué)，2011,39（17）：10513-10515.

[10] 呂 新，張 偉，胡昌浩，等.玉米種植區(qū)氣候生態(tài)因素優(yōu)勢綜合評價研究[J].干旱區(qū)研究，2005，（3）：387-390.