孫建軍，丁超明，申培林，郭新海，郭玉璽，王德新

（河南省農(nóng)業(yè)科學院，鄭州 450002）

氮素是植物生長發(fā)育所必需的營養(yǎng)元素，合理有效地施用氮肥在提高作物產(chǎn)量的同時會改善其品質(zhì)［1］。21 世紀以來，人口、資源、環(huán)境的矛盾日益突出，中國農(nóng)業(yè)生態(tài)環(huán)境面臨多方面的嚴峻挑戰(zhàn)，施肥對環(huán)境的影響受到越來越多的關注［2］。主要農(nóng)作物以增加施氮量來提高產(chǎn)量，但氮肥過多施用不僅增產(chǎn)效果不顯著，還會造成環(huán)境污染，給社會帶來不利影響。在施氮量較低條件下，產(chǎn)量和氮肥利用率較高的品種具有較高的干物質(zhì)積累量、莖蘗成穗率、葉面積指數(shù)、粒葉比及群體生長速率［3］。在小麥栽培中，適宜濃度的氮肥能夠促進小麥農(nóng)藝性狀，提高小麥的品質(zhì)和產(chǎn)量，施肥水平過低或過高均不利于小麥農(nóng)藝性狀特別是產(chǎn)量相關性狀的表達。

學者們對GGE 雙標圖在作物區(qū)試數(shù)據(jù)分析中的優(yōu)勢進行了分析，并提出在用于分析區(qū)試數(shù)據(jù)的同時，還可對任何2 位數(shù)據(jù)進行分析［4-6］。GGE 雙標圖與其他圖解法相比，具有獨特的優(yōu)勢，它是綜合品種總體效應（G）和品種×環(huán)境互作（GE）的數(shù)學模型［7］。GGE 雙標圖在鮮食糯玉米［8］和普通玉米［9］、水稻［10］、花生［11］、冬油菜［12］等作物的區(qū)試和不同類二維數(shù)據(jù)的分析中得以應用。本研究利用方差分析、多重比較和GGE 雙標圖分析不同施氮量對小麥農(nóng)藝性狀表達的影響，為小麥的合理施氮提供參考。

1 材料與方法

1.1 試驗材料

供試品種為半冬性、中熟品種周麥16，來源于河南省周口市農(nóng)業(yè)科學院。

1.2 試驗設計

試驗于2020—2021 年在河南省農(nóng)業(yè)科學院現(xiàn)代化研究基地進行，前茬作物為玉米（同一品種，非小區(qū)試驗），地勢平坦、肥力均勻、排灌方便，測定試驗地土壤有機質(zhì)10.9 mg∕kg，堿解氮41.2 mg∕kg，速效磷27.52mg∕kg，速效鉀97.35mg∕kg，基肥在播種前施入。肥力水平設置5 個處理（表1），后期追肥保持一致，組間隨機排列，設3 次重復，小區(qū)面積7.2 m2，6 行區(qū)，播種量一致。

表1 不同氮肥施用量處理

1.3 田間管理

生育期內(nèi)根據(jù)土壤干旱情況適當灌水，各小區(qū)灌水量一致，選擇合理的化學藥劑噴殺病蟲害。

1.4 性狀調(diào)查

在小麥整個生育期內(nèi)，按照《植物新品種特異性、一致性和穩(wěn)定性測試指南普通小麥小麥DUS 測試指南GB∕T 19557.2—2017》進行觀測，取20 株進行性狀調(diào)查。成熟后選擇典型株40株進行取樣、考種。

1.5 數(shù)據(jù)處理

采用Excel 2016 整理試驗數(shù)據(jù)，采用SAS 軟件進行數(shù)據(jù)方差分析和多重比較，用AI軟件處理圖片，利用R 語言的GGEBiplot GUI 包制作GGE 雙標圖。

2 結果與分析

2.1 不同年份、施氮量對小麥農(nóng)藝性狀的影響

不同施氮量、年份與小麥株高、穗長度、旗葉寬、旗葉長、不育小穗率、千粒重和穗粒重的差異顯著性見表2。由表2 可知，小麥所有測試性狀在不同施氮量處理間差異均達到極顯著水平（P＜0.01），而年份間差異不顯著，年份和處理互作中僅有穗長度差異達到顯著差異水平（P＜0.05）。

表2 不同年份、處理的農(nóng)藝性狀的相關系數(shù)

根據(jù)方差分析結果，由于2 年間各農(nóng)藝性狀差異不顯著，只對2 年農(nóng)藝性狀數(shù)據(jù)處理間進行多重比較，結果見表3。由表3 可知，株高、穗長度隨肥力水平的增加呈先增加后降低的趨勢，N2 處理株高最高，為80.83 cm，N3 處理穗長度最長，為9.72 cm；旗葉長和旗葉寬均隨施氮量先增后減，在N4 處理時均為最大；不育小穗率則隨施氮量的增加而降低，N1處理量最大，為16.34%；千粒重在N4 和N5 處理間呈顯著差異。

表3 不同處理的小麥農(nóng)藝性狀

2.2 試驗處理與農(nóng)藝性狀GGE 雙標圖分析

GGE 雙標圖法可直觀地判斷不同處理下表現(xiàn)最好的指標以及各指標的相互關系。將原始數(shù)據(jù)減去各指標的均值后形成數(shù)據(jù)集，使數(shù)據(jù)中只含有處理主效應G 和處理與各指標互作效應GE，對GGE 做單值分解，以第一主成分和第二主成分為代表，按照第一和第二主成分值將所有處理和指標繪于二維圖上即形成GGE 雙標圖。對采用的各指標數(shù)據(jù)進行標準化處理后，可認為這些指標在各處理中具有同等重要的作用。把處理的標志點（N1、N2 和N5）用直線連接起來，形成一個把所有處理都囊括在內(nèi)的三角形，從原點（0，0）起作各邊的垂線，把整個雙標圖分為3 個扇區(qū)，不同指標落于各個區(qū)內(nèi)，每個區(qū)內(nèi)的頂角處理即為該扇區(qū)內(nèi)所有指標的最優(yōu)處理。試驗處理與農(nóng)藝性狀GGE 雙標圖分析見圖1。

由圖1a可知，第一主成分的效應為PC1=70.91%，第二主成分的效應為PC2=23.11%，GGE 雙標圖上展示了各處理真實信息的94.02%。周麥16 的7 個農(nóng)藝性狀分別落于3 個扇區(qū)中。N1 區(qū)內(nèi)不育小穗率指標值最高。N2 區(qū)內(nèi)株高指標值最高。其中N3、N4 處理包含在N5 區(qū)域內(nèi)，穗長度、旗葉寬、旗葉長、千粒重、穗粒重指標值最高。

圖1b 中單箭頭的橫軸代表各處理在所有性狀下的影響程度，軸上的小圓圈代表平均值，小圓圈箭頭方向表示影響程度，越向右影響程度越高。從圖1b 可知，N4 和N5 在箭頭的右測，均高于均值，其中N5 影響程度最高，隨后分別為N4、N3、N2、N1?？v軸代表各處理與各性狀相互作用的傾向性，處理越偏離平均軸越不穩(wěn)定，而且具有負作用。其中N2處理最不穩(wěn)定，N3、N4 處理影響作用則較穩(wěn)定。

圖1c 中顯示性狀的區(qū)分力和代表性。GGE 雙標圖中農(nóng)藝性狀向量間的夾角可以判斷農(nóng)藝性狀間的相關性。從坐標原點到農(nóng)藝性狀標志點的連線作為農(nóng)藝性狀向量，以某一個農(nóng)藝性狀向量為起始，順時針方向旋轉，其他農(nóng)藝性狀向量與該線的夾角余弦值即為二者農(nóng)藝性狀間的相關系數(shù)。農(nóng)藝性狀向量與平均軸的角度是其對目標農(nóng)藝性狀代表性的度量，代表性與角度呈負相關，代表性隨著角度變小而變強，隨著角度變大而變?nèi)?。向量的長度代表其區(qū)分力，向量長度越長，表明農(nóng)藝性狀對各處理的表現(xiàn)具有較強的區(qū)分能力。由圖1c 可知，穗長度、旗葉寬、旗葉長、和穗粒重與平均環(huán)境軸的夾角較小，而且向量長度相對較長，表明對各處理差異的區(qū)分力和目標農(nóng)藝性狀的代表性都較好，同時旗葉寬和穗粒重間存在較大的相關性。不育小穗率呈負相關。千粒重和株高具有較長的向量長度，表明對各施肥處理具有較高的區(qū)分力，而其夾角過大表明其代表性較弱。

圖1d 表示理想氮肥施用量的選擇，以最佳氮肥量為圓心做多層同心圓，根據(jù)與最佳氮肥量的接近程度，GGE 雙標圖對各處理的效應優(yōu)劣進行排序。越靠近同心圓中心，則表示該處理結果越準確，反之亦然。從圖1d可知，N4靠近最小同心圓，說明其是綜合表現(xiàn)最好；N1、N2 距離同心圓較遠，表現(xiàn)相對較差。

圖1 試驗處理與農(nóng)藝性狀GGE 雙標圖分析

3 小結與討論

氮素是小麥生長發(fā)育所必需的大量營養(yǎng)元素。調(diào)控小麥生長發(fā)育的多個過程，合理運用氮肥可以協(xié)調(diào)小麥產(chǎn)量與品質(zhì)的關系［13］?；輹喳惖龋?4］研究發(fā)現(xiàn)，與不施氮相比，長期施氮小麥平均增產(chǎn)67.1%，生物量提高52.0%，收獲指數(shù)提高9.5%，穗數(shù)和穗粒數(shù)分別提高32.5% 和40.0%，千粒重下降7.1%。宮明波等［15］研究發(fā)現(xiàn)，隨施氮量的增加，葉綠素呈增長趨勢，但葉面積、地上干物質(zhì)量、穗數(shù)、穗粒數(shù)及產(chǎn)量均呈先增后減的趨勢；蛋白質(zhì)及濕面筋含量隨施氮量增加而提高，氮肥過量施用會產(chǎn)生鹽害，抑制小麥生長，導致減產(chǎn)。

GGE 雙標圖在國內(nèi)外作物區(qū)試試驗中應用廣泛［16-18］，它可以將不同環(huán)境和因素及因素之間的關系以突變的形式直觀表現(xiàn)出來。本研究發(fā)現(xiàn)，處理間各農(nóng)藝性狀差異均達到顯著性水平，而年份間差異均不顯著，小麥穗長度、旗葉寬、旗葉長、千粒重和穗粒重隨著氮肥量增加，整體均呈增長趨勢，但后期增長趨勢不明顯，其中千粒重在N4 和N5 處理間差異顯著。株高隨施氮量的增加呈先增加后降低趨勢，在N2 處理時達到最大。不育小穗率則隨施氮量的增加而降低，年份和處理互作中，僅穗長度達到顯著差異；利用GGE 雙標圖分析方法綜合評價出施肥量在N3 和N4 處理下，各個性狀的表現(xiàn)最為理想。影響小麥農(nóng)藝性狀標的因素不僅有氮素，還有種植密度、播期、氣候等，這些方面需要做進一步研究。