基于GGE雙標(biāo)圖對(duì)苜蓿品種豐產(chǎn)性和穩(wěn)定性綜合評(píng)價(jià)

陳彩錦，張尚沛，師尚禮，曾燕霞，吳 娟，尚繼紅，高 婷，張 蓉

(1.甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院，甘肅 蘭州 730070；2.寧夏農(nóng)林科學(xué)院固原分院，寧夏 固原 756000；3.寧夏農(nóng)林科學(xué)院動(dòng)物科學(xué)研究所，寧夏 銀川 756002；4.寧夏農(nóng)林科學(xué)院植物保護(hù)研究所，寧夏 銀川 756002)

紫花苜蓿(MedicagosativaL.，以下簡稱苜蓿)營養(yǎng)價(jià)值高、適口性好，具有“牧草之王”的美譽(yù)[1]，在世界范圍內(nèi)作為栽培作物或飼草廣泛種植。近幾年，隨著草食畜牧業(yè)的快速發(fā)展，寧夏南部山區(qū)(以下簡稱寧南山區(qū))優(yōu)質(zhì)苜蓿的需求量逐年攀升，種植面積也在大范圍擴(kuò)增。尤其在“十三五”農(nóng)業(yè)產(chǎn)業(yè)發(fā)展及草畜產(chǎn)業(yè)成為自治區(qū)主導(dǎo)產(chǎn)業(yè)之后，寧南山區(qū)由于獨(dú)特的生境成為寧夏苜蓿產(chǎn)業(yè)發(fā)展的主戰(zhàn)場[2]。但是，苜蓿生產(chǎn)在該區(qū)域長期以來存在品種退化、混雜、產(chǎn)量低、抗逆性及品質(zhì)差等突出問題[3]。因此，要實(shí)現(xiàn)苜蓿的豐產(chǎn)及穩(wěn)產(chǎn)，使其在畜牧業(yè)中充分發(fā)揮作用，首先必須選育或篩選出適宜不同區(qū)域的品種[4]。而牧草新品種選育或篩選需要對(duì)品種開展區(qū)域試驗(yàn)與分析，使得該品種在最適宜的區(qū)域獲得高產(chǎn)和穩(wěn)產(chǎn)，以提高牧草的綜合生產(chǎn)力[5]。品種區(qū)域試驗(yàn)分析的中心任務(wù)是正確把握品種與環(huán)境互作的模式[6]，GGE雙標(biāo)圖(GGE biplot)是分析多因素相互作用的一種方法，常用于分析基因與環(huán)境互作關(guān)系、不同處理?xiàng)l件下各指標(biāo)間的關(guān)系等[7]。它以圖解的方式，借助輔助線有效反映出何品種在何種環(huán)境下表現(xiàn)最好，有利于清晰鑒別出高產(chǎn)、穩(wěn)產(chǎn)的品種及有效分辨力的環(huán)境[8]。GGE雙標(biāo)圖在作物上應(yīng)用較多，其中柴繼寬[4]、李四龍[8]、王磊[9]、尚毅[10]、王曉斌[11]、張海娟[12]、嚴(yán)威凱[13]、許乃銀[14]、張勇[15]等人利用GGE雙標(biāo)圖分別對(duì)燕麥(AvenasativaL.)、花生(ArachishypogaeaLinn.)、水稻(OryzasativaL.)、甘藍(lán)型油菜(BrassicanapusL.)、馬鈴薯(SolanumtuberosumL.)、大麥(HordeumvulgareL.)、大豆(Glycinemax(Linn.)Merr.)、棉花(Gossypiumspp)、小麥(TriticumaestivumL.)等作物的品種、種質(zhì)資源、雜交組合等開展了分析研究，確定出了具有適應(yīng)性、豐產(chǎn)性、穩(wěn)定性的品種、資源、雜交材料及經(jīng)濟(jì)效益最高的種植密度，篩選出了具有代表性、鑒別力的區(qū)域等。GGE雙標(biāo)圖在牧草方面應(yīng)用主要集中在青貯玉米、飼用小黑麥等方面，其中廖正巧[16]利用GGE雙標(biāo)圖分析了青貯玉米19份材料穩(wěn)定性和豐產(chǎn)性，綜合篩選出了3份高產(chǎn)、穩(wěn)產(chǎn)的品種；劉晶等人[17]利用GGE雙標(biāo)圖研究了4個(gè)小黑麥品種基因型與環(huán)境互作效應(yīng)的穩(wěn)定性，篩選出了草產(chǎn)量最高、穩(wěn)定、營養(yǎng)品質(zhì)最優(yōu)的小黑麥新品系及其適宜種植的區(qū)域?？傊?，基因型主效應(yīng)(G)加基因型與環(huán)境互作(G×E)效應(yīng)的GGE雙標(biāo)圖模型與變異系數(shù)法、適定性參數(shù)法、回歸分析法、高穩(wěn)系數(shù)法、加性主效應(yīng)乘積交互作用(additive main effects and multiplicative interaction，AMMI)模型相比較，不僅在分析基因型與環(huán)境互作(G×E)方面具有分析結(jié)果精確的優(yōu)點(diǎn)，而且能鑒定篩選出高產(chǎn)、穩(wěn)產(chǎn)、適應(yīng)性廣的理想品種，并對(duì)試點(diǎn)辨別力和代表性進(jìn)行篩選，對(duì)品種合理布局做出評(píng)價(jià)[18-19]。因此，本研究應(yīng)用此方法對(duì)引進(jìn)的6個(gè)紫花苜蓿品種在寧南山區(qū)3個(gè)不同區(qū)域的適應(yīng)性、穩(wěn)定性，以及試點(diǎn)的代表性、區(qū)分力進(jìn)行分析研究，以期篩選出性狀優(yōu)良、適應(yīng)性強(qiáng)、穩(wěn)定性和豐產(chǎn)性高的品種及品種區(qū)域試驗(yàn)的理想環(huán)境[8]，為區(qū)域苜蓿種植提供技術(shù)依據(jù)及數(shù)據(jù)參考，也為當(dāng)?shù)氐牟菽翗I(yè)發(fā)展提供品種基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)地基本概況

試驗(yàn)地位于寧夏固原市原州區(qū)頭營鎮(zhèn)徐河村，原州區(qū)彭堡鎮(zhèn)彭堡三隊(duì)，以及寧夏固原市隆德縣觀莊鄉(xiāng)前莊村寧夏農(nóng)林科學(xué)院試驗(yàn)基地，3個(gè)不同地區(qū)具體地理環(huán)境狀況如表1所示。

表1 試驗(yàn)地地理環(huán)境狀況

1.2 試驗(yàn)設(shè)計(jì)

本試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì)，留苗密度60 萬株·hm-2，試驗(yàn)6個(gè)品種，每個(gè)品種設(shè)置3次重復(fù)，共18個(gè)小區(qū)，小區(qū)面積15 m2，小區(qū)行長5 m，每個(gè)小區(qū)13行，行間距25 cm，區(qū)間距80 cm，重復(fù)間距100 cm，小區(qū)四周設(shè)有保護(hù)行，人工開溝拉線條播。2015年進(jìn)行試驗(yàn)的建植及養(yǎng)護(hù)管理，其中E1，E2，E3試驗(yàn)點(diǎn)土壤基礎(chǔ)理化性質(zhì)分別為：E1土壤類型為黃綿土，有機(jī)質(zhì)10.95 g·kg-1，全鹽含量為1.24 g·kg-1，速效氮含量為50.10 mg·kg-1，速效磷含量為10.27 mg·kg-1，速效鉀含量為243.17 mg·kg-1，耕作層pH 7.00，前茬為蕎麥；E2土壤類型為黑壚土，有機(jī)質(zhì)含量為40.60 g·kg-1，全鹽含量為0.69 g·kg-1，速效氮含量為98.00 mg·kg-1，速效磷含量為52.20 mg·kg-1，速效鉀含量為265.00 mg·kg-1，耕作層pH 7.80，前茬為冬小麥；E3土壤類型為淺黑壚土，有機(jī)質(zhì)含量為42.10 g·kg-1，全鹽含量為0.44 g·kg-1，速效氮含量為104.60 mg·kg-1，速效磷含量為73.00 mg·kg-1，速效鉀含量為318.06 mg·kg-1，耕作層pH 7.80，前茬為馬鈴薯。試驗(yàn)地土壤肥力和施肥水平等同于大田生產(chǎn)，且肥力均勻。播前秋季翻耕1次，翻深20～30 cm，第二年春季解凍后旋耕耙耱各1次，基施磷酸二銨225 kg·hm-2之后播種，在整個(gè)試驗(yàn)期間，均無灌溉，但由于出苗期間降雨導(dǎo)致板結(jié)，人工破板1次，并對(duì)試驗(yàn)的斷行進(jìn)行補(bǔ)苗1次，其余管理也等同于大田生產(chǎn)。2016—2017年進(jìn)行相關(guān)產(chǎn)量指標(biāo)測定，測定方法：每個(gè)品種3次重復(fù)，每重復(fù)選3個(gè)樣方，共9個(gè)樣方，樣方在初花期刈割前選擇1 m2進(jìn)行刈割，留茬高度0.06～0.08 m，稱鮮重，然后自然風(fēng)干并稱干重。

1.3 參試品種

參試品種共6個(gè)，品種具體名稱、代碼及來源見表2。

表2 苜蓿品種及來源

1.4 數(shù)據(jù)統(tǒng)計(jì)與分析

本試驗(yàn)的相關(guān)數(shù)據(jù)的初步處理采用Excel軟件；采用GGE雙標(biāo)圖對(duì)各品種和環(huán)境的最佳適應(yīng)區(qū)域，品種的豐產(chǎn)性、穩(wěn)定性，以及試點(diǎn)的代表性及區(qū)分力進(jìn)行分析；聯(lián)合方差分析采用DPS 16.5高級(jí)版統(tǒng)計(jì)分析軟件進(jìn)行統(tǒng)計(jì)分析。

1.5 分析方法

GGE雙標(biāo)圖最初由Gabriel提出，它有2部分構(gòu)成，分別為雙標(biāo)圖和GGE[20-22]，是利用雙標(biāo)圖表示基因(G)與基因環(huán)境互作(G×E)兩個(gè)因素之間的關(guān)系(G+GE)，雙標(biāo)圖屬于一種2維矩陣圖，通過對(duì)含有品種環(huán)境的2維矩陣進(jìn)行特征值分解(singular value decomposition,SVD)，將一個(gè)含有品種環(huán)境的區(qū)域試驗(yàn)數(shù)據(jù)集進(jìn)行主成分分析，其中解釋變異最多的主成分叫第一主成分(PC1)，解釋變異處于第二位的叫第二主成分(PC2)[23]。

其中在不同區(qū)域試驗(yàn)中，產(chǎn)量一般由以下幾個(gè)組分組成：

Yge=μ+αg+βe+φge+εge

其中μ為總體平均值，αg為基因型主效應(yīng)，βe為環(huán)境主效應(yīng)，φge為基因型與環(huán)境互作效應(yīng)，εge為總體誤差，Yge是基因型g在環(huán)境e中的性狀值[4]。假定地點(diǎn)都同等重要，則Yge減去各試點(diǎn)均值，就剩基因型主效αg和基因型與環(huán)境互作效應(yīng)φge，合稱GGE[4]。GGE雙標(biāo)圖不僅在同一個(gè)生態(tài)區(qū)域內(nèi)，同時(shí)顯示各試驗(yàn)環(huán)境對(duì)品種的區(qū)分力和對(duì)目標(biāo)環(huán)境的代表性，顯示各個(gè)品種的高產(chǎn)穩(wěn)定性，還在品種區(qū)域試驗(yàn)中可以把環(huán)境直觀劃分為若干品種生態(tài)區(qū)，并揭示各生態(tài)區(qū)內(nèi)最適應(yīng)的品種[21]。

2 結(jié)果與分析

2.1 產(chǎn)量表現(xiàn)

由表3可知，在3個(gè)不同的試驗(yàn)區(qū)，2016—2017年2年度苜蓿品種產(chǎn)量平均值最高的試驗(yàn)區(qū)是E2，為14 816.7 kg·hm-2，E3試驗(yàn)區(qū)最低，為10 314.0 kg·hm-2，且前者比后者高43.7%；在參試品種中，2年度平均產(chǎn)量最高的品種是‘甘農(nóng)7號(hào)’，最低的品種是‘阿爾岡金’，前者比后者高44.9%。

表3 寧南山區(qū)2016—2017年苜蓿品種產(chǎn)量平均值

2.2 各品種的變異特性

對(duì)表4不同試驗(yàn)區(qū)苜蓿品種產(chǎn)量結(jié)果聯(lián)合方差分析可知：基因型與環(huán)境互作(G×E)、環(huán)境(E)、基因型(G)內(nèi)部之間均存在顯著的差異(P<0.05)，且產(chǎn)量引起的基因(G)變異、環(huán)境(E)變異、基因型與環(huán)境互作(G×E)變異的平方和分別占總處理的平方和的12.3%，16.3%，19.4%，其中基因型與環(huán)境互作(G×E)變異是3種變異中最主要的變異。這表明，基因型與環(huán)境互作(G×E)、環(huán)境(E)、基因型(G)對(duì)苜蓿品種產(chǎn)量影響顯著，其中的基因型與環(huán)境(G×E)互作是引起苜蓿各品種產(chǎn)量之間差異的最主要原因，環(huán)境(E)因素次之，基因型(G)因素最小。

表4 品種試驗(yàn)方差分析

2.3 各品種與環(huán)境的最佳適應(yīng)區(qū)域

GGE雙標(biāo)圖是把一些“角頂”品種的標(biāo)志點(diǎn)用直線連接起來形成一個(gè)多邊形，把其他品種都囊括在內(nèi)，從原點(diǎn)起作各邊的垂線，把整個(gè)疊圖分為若干扇形區(qū)域。每個(gè)試驗(yàn)地點(diǎn)自然落于某個(gè)區(qū)內(nèi)，每個(gè)區(qū)內(nèi)的“角頂”品種就是這個(gè)區(qū)內(nèi)在每一試驗(yàn)地點(diǎn)上表現(xiàn)最好的品種，位于多邊形內(nèi)部，靠近原點(diǎn)的品種是接近平均產(chǎn)量且對(duì)環(huán)境變化不敏感的品種。如圖1所示，GGE雙標(biāo)圖被劃分為3個(gè)區(qū)，設(shè)以E1，E2試點(diǎn)落于的扇形區(qū)為第1區(qū)，順時(shí)針旋轉(zhuǎn)依次為第2區(qū)、第3區(qū)。E1，E2落于第1扇形區(qū)，其中產(chǎn)量最高的是G3，表現(xiàn)最好，E3落于第一扇形和第二扇形的交界處，且在第二扇形區(qū)產(chǎn)量最高、表現(xiàn)最好的品種是G2，G5品種沒有地點(diǎn)可以落到，說明此品種在所有試驗(yàn)地點(diǎn)表現(xiàn)均不是最好的。由于以上3個(gè)試驗(yàn)地點(diǎn)在2個(gè)扇形區(qū)，且其中1個(gè)試驗(yàn)地點(diǎn)落在2個(gè)生態(tài)區(qū)域的交接處，另外2個(gè)落于1個(gè)生態(tài)區(qū)，由此可以將生態(tài)區(qū)劃分為2個(gè)，以E1，E2為代表的1個(gè)，以E3為代表的1個(gè)。品種的產(chǎn)量可能跟環(huán)境溫度、降雨量等存在一定的相關(guān)性。因此，G3品種可能耐旱性較好，在寧南半干旱區(qū)產(chǎn)量表現(xiàn)較好，G2品種可能抗寒性強(qiáng)，返青早，從而在寧南半陰濕冷涼地區(qū)產(chǎn)量表現(xiàn)較好。

圖1 GGE雙標(biāo)圖對(duì)苜蓿品種的適應(yīng)性分析

2.4 品種的豐產(chǎn)性與穩(wěn)定性

在特定品種生態(tài)區(qū)內(nèi)，理想的品種既是高產(chǎn)的又是穩(wěn)產(chǎn)的。在GGE雙標(biāo)圖中，X軸(第一主成分得分-PC1)代表了品種的平均產(chǎn)量，Y軸(第二主成分得分-PC2)代表了品種的穩(wěn)定性(即偏離平均產(chǎn)量的傾向)，此值越接近于0，穩(wěn)定性越好[6]。如圖2所示，參試苜蓿品種平均產(chǎn)量效應(yīng)為91.5%，品種的穩(wěn)定性(PC2)效應(yīng)為7.7%，且這兩個(gè)效應(yīng)綜合了大部分的變異信息(99.2%)。因此，根據(jù)此結(jié)果推斷以上的分析具有一定的可靠性。圖中圓圈代表的是平均環(huán)境，經(jīng)過原點(diǎn)和平均環(huán)境的這條帶箭頭的直線為平均環(huán)境軸(AEC)，在平均環(huán)境軸上沿著箭頭的方向產(chǎn)量逐漸增高。且在AEC縱軸左邊的品種產(chǎn)量低于平均值，在縱軸右邊的品種產(chǎn)量高于平均值。因此產(chǎn)量從高到低依次是G3>G6>G4>G1>G2>G5，且G3，G6，G4產(chǎn)量高于平均值。AEC縱軸代表了品種的穩(wěn)定性，AEC橫軸上垂線的長短顯示品種穩(wěn)定性的大小[22]。圖2顯示，穩(wěn)定性最好的品種是G5，G4，G6，但G5品種產(chǎn)量低于平均值，是穩(wěn)定低產(chǎn)的品種，因而是表現(xiàn)最差品種；G3品種產(chǎn)量最高，但不穩(wěn)定，屬于高產(chǎn)不穩(wěn)定品種。因此綜合高產(chǎn)及穩(wěn)定性，表現(xiàn)最好的是G4，G6品種。

圖2 GGE雙標(biāo)圖對(duì)苜蓿品種的穩(wěn)定性和豐產(chǎn)性分析

2.5 區(qū)分力與代表性

理想的試驗(yàn)環(huán)境應(yīng)當(dāng)具備對(duì)目標(biāo)生態(tài)區(qū)的代表性和參試品種較強(qiáng)區(qū)分力。反映在圖示中，理想的試點(diǎn)應(yīng)具有較大的第一主成分值和接近于0的第二主成分值，這樣的試點(diǎn)能有效鑒別出高產(chǎn)且穩(wěn)定的品種[24]。圖3中的小圓圈代表的是“平均環(huán)境”，它的位置取決于各試驗(yàn)環(huán)境坐標(biāo)的平均值，代表目標(biāo)環(huán)境，帶箭頭且通過原點(diǎn)和平均環(huán)境的稱為平均環(huán)境向量或平均環(huán)境軸(AEA)，各試驗(yàn)點(diǎn)向量的長度是區(qū)分能力的度量，向量越長，對(duì)品種的區(qū)分力越強(qiáng)。試驗(yàn)點(diǎn)向量與平均環(huán)境向量的角度是其對(duì)目標(biāo)環(huán)境的代表性的度量，角度越小代表性越強(qiáng)，角度越大代表性越弱[21]，如果兩個(gè)之間的夾角為鈍角，則這個(gè)地方不適合作為試驗(yàn)點(diǎn)，為此，平均環(huán)境軸箭頭所指方向是對(duì)試驗(yàn)點(diǎn)區(qū)分力和代表性的綜合評(píng)價(jià)。在圖3中，3個(gè)試驗(yàn)地點(diǎn)與平均環(huán)境軸之間的夾角都明顯小于90°，其中E2

圖3 用GGE雙標(biāo)圖分析苜蓿各品種的代表性和區(qū)分力

3 討論

產(chǎn)量是飼草作物生產(chǎn)潛力評(píng)價(jià)最有效的方法和飼草生產(chǎn)中的重要指標(biāo)[25]。在本研究中，2年度苜蓿平均產(chǎn)量最高的地點(diǎn)是E2，比E1，E3地點(diǎn)高12.5%，30.38%，最高的品種是‘甘農(nóng)7號(hào)’，比產(chǎn)量最低品種‘阿爾岡金’高31.0%。以上不同品種在不同區(qū)域內(nèi)的產(chǎn)量表現(xiàn)存在一定的差異，主要是因?yàn)槠贩N與環(huán)境之間存在著一個(gè)適應(yīng)性問題，不同的品種控制其產(chǎn)量的基因不同，外界因素的影響在不同的試點(diǎn)有選擇的進(jìn)行表達(dá)，這與前人的研究結(jié)果相符[26]。同一個(gè)品種在不同區(qū)域的表現(xiàn)也存在差異，在本研究中，同一品種在3個(gè)試點(diǎn)產(chǎn)量最大相差21.6%～51.0%，這表明一個(gè)品種產(chǎn)量表現(xiàn)除自身的基因(G)決定之外，環(huán)境(E)、基因型與環(huán)境的交互作用(G×E)也起關(guān)鍵作用。本研究中基因型(G)、環(huán)境(E)、基因型與環(huán)境的交互作用(G×E)對(duì)產(chǎn)量的影響均表現(xiàn)為顯著，且基因型與環(huán)境的交互作用(G×E)影響最大，這與孫萬斌等人[27]的研究結(jié)果相符。因此，在種植品種的選擇中，需根據(jù)品種的特性因地制宜選擇具體的環(huán)境，只有將品種與環(huán)境耦合起來，才能最大限度的發(fā)揮地域優(yōu)勢及品種的生產(chǎn)潛力。

精確的試點(diǎn)對(duì)于品種評(píng)價(jià)至關(guān)重要。一個(gè)合適的試點(diǎn)決定了一個(gè)品種的適應(yīng)性、推廣價(jià)值和利用價(jià)值。本研究中利用GGE雙標(biāo)圖模型將3個(gè)試點(diǎn)劃分為2個(gè)不同的生態(tài)區(qū)域，其中E1，E2為1個(gè)，E3為1個(gè)。這與當(dāng)前的具體實(shí)際環(huán)境存在一致性，其中E1，E2位于寧南半干旱區(qū)域，E3位于寧南半陰濕冷涼區(qū)域。E1，E2試點(diǎn)同一品種的產(chǎn)量也存在較大的差異，這可能是因?yàn)橐陨?個(gè)區(qū)域降雨量雖然相同，但其余氣候因子、土壤因子及生物因子等還存在差異，或者是因?yàn)橥黄贩N受2個(gè)地區(qū)相同環(huán)境的影響程度不同，從而產(chǎn)量存在較大差異，以上可能原因還有待于進(jìn)一步研究確定。試點(diǎn)的選擇既要考慮對(duì)各產(chǎn)區(qū)環(huán)境的代表性，也要考慮對(duì)品種差異的區(qū)分能力[4]。利用GGE雙標(biāo)圖對(duì)本研究中3個(gè)試點(diǎn)的代表性和區(qū)分力進(jìn)行分析可知，E2的代表性最好，但區(qū)分力最差，E1的代表性處于3個(gè)試點(diǎn)的第2位，但區(qū)分力最好，E3代表性最差，但區(qū)分力處于3個(gè)試點(diǎn)的第2位。E2試點(diǎn)因?qū)ζ贩N沒有區(qū)分力，此試點(diǎn)作用不大，可以淘汰，E3有一定的區(qū)分能力但代表性不強(qiáng)，可以用于淘汰G3等不穩(wěn)定品種，但不能用于優(yōu)良品種的選擇，因此3個(gè)試點(diǎn)中最理想的試點(diǎn)是E1。以上試點(diǎn)出現(xiàn)代表性很好，但對(duì)品種沒有區(qū)分力的原因可能是包括環(huán)境因素在內(nèi)其它因素引起的，因此對(duì)于品種對(duì)試點(diǎn)的適應(yīng)需進(jìn)行多年多點(diǎn)試驗(yàn)研究。

豐產(chǎn)性和穩(wěn)定性也是決定一個(gè)品種是否有推廣利用價(jià)值的重要指標(biāo)。在本研究中，利用GGE雙標(biāo)圖模型，篩選出了飼草高產(chǎn)和穩(wěn)產(chǎn)的品種G4和G6，也篩選出了在E1試點(diǎn)高產(chǎn)的品種G3，同時(shí)也得到了一個(gè)結(jié)果，豐產(chǎn)性和穩(wěn)定性之間存在一定差異，本研究中最穩(wěn)定的品種G5產(chǎn)量最低，是穩(wěn)定低產(chǎn)的品種。因此，在寧南山區(qū)大范圍推廣應(yīng)用的品種首選G4和G6。但在局部半干旱區(qū)域，G3品種也是一個(gè)很好的選擇?？傊?，為使品種的區(qū)域布局合理化，發(fā)揮增產(chǎn)潛力，應(yīng)對(duì)品種最佳適應(yīng)區(qū)域進(jìn)行分析[28]。

以上僅為苜蓿品種建植后2年的數(shù)據(jù)，結(jié)果在當(dāng)?shù)鼐哂幸欢ǖ膮⒖純r(jià)值，但需進(jìn)一步試驗(yàn)驗(yàn)證。同時(shí)品種的適應(yīng)性還需綜合分析品種生育期、農(nóng)藝性狀、抗逆性、抗病性等方面，環(huán)境因素還需分析各試驗(yàn)點(diǎn)土壤理化性質(zhì)等方面，這些研究將在后續(xù)進(jìn)一步開展。

4 結(jié)論

本研究結(jié)果表明，基因型(G)、環(huán)境(E)、基因型與環(huán)境互作(G×E)3者引起的變異對(duì)苜蓿各品種產(chǎn)量的影響顯著，但基因與環(huán)境互作效應(yīng)是引起變異的最主要原因。同時(shí)利用GGE雙標(biāo)圖模型分析可知，在寧南山區(qū)高產(chǎn)穩(wěn)產(chǎn)的苜蓿品種是‘三得利’和‘苜蓿王’；試點(diǎn)頭營鎮(zhèn)徐河村和彭堡鎮(zhèn)彭堡三隊(duì)生態(tài)區(qū)相似，與觀莊鄉(xiāng)前莊村不同，最理想試點(diǎn)為頭營鎮(zhèn)徐河村。