豫中地區(qū)冬小麥臨界氮稀釋曲線與氮營養(yǎng)指數(shù)模型的建立

岳松華 劉春雨 黃玉芳 葉優(yōu)良河南農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，河南鄭州 450002

?

豫中地區(qū)冬小麥臨界氮稀釋曲線與氮營養(yǎng)指數(shù)模型的建立

岳松華 劉春雨 黃玉芳 葉優(yōu)良*
河南農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，河南鄭州 450002

摘 要：為了驗(yàn)證小麥臨界氮稀釋曲線在豫中地區(qū)的適用性，選擇大穗型品種周麥16和多穗型品種豫麥49-198進(jìn)行了連續(xù)3年（2011年10月至2014年6月）的田間定位試驗(yàn)，設(shè)置5個(gè)氮水平（純氮0、120、180、240、360 kg hm-2），分別構(gòu)建了小麥地上部生物量臨界氮濃度稀釋曲線模型和氮營養(yǎng)指數(shù)模型。結(jié)果表明，地上部生物量隨著施氮量增加顯著提高，而后趨于平緩；小麥臨界氮濃度與地上部生物量符合冪函數(shù)關(guān)系，其 R2分別為 0.8203 （豫麥49-198，P＜0.01）和0.7981 （周麥16，P＜0.01），且模型在年度間具有較好的穩(wěn)定性。根據(jù)模型得到的氮營養(yǎng)指數(shù)可以較好地評(píng)價(jià)小麥植株氮營養(yǎng)狀況，且小麥氮營養(yǎng)指數(shù)和相對(duì)莖基部硝酸鹽之間具有顯著的線性相關(guān)。說明本文所構(gòu)建的小麥臨界氮濃度稀釋模型適用于本地區(qū)，且可以用來估測(cè)小麥植株氮素盈虧水平。

關(guān)鍵詞：小麥；臨界氮稀釋曲線；氮營養(yǎng)指數(shù)；莖基部硝酸鹽

本研究由國家自然科學(xué)基金項(xiàng)目（31471935）和國家公益性行業(yè)（農(nóng)業(yè)）科研專項(xiàng)（201103003）資助。

This study was supported by the National Natural Science foundation of China （31471935） and the Special Fund for Agro-scientific Research in the Public Interest （201103003）

第一作者聯(lián)系方式∶ E-mail∶ 1069342241@qq.com

URL∶ http∶//www.cnki.net/kcms/detail/11.1809.S.20160321.1056.008.html

小麥?zhǔn)俏覈蠹Z食作物之一，河南省是我國的主要糧食生產(chǎn)基地，其中豫中地區(qū)的小麥種植面積和產(chǎn)量分別占全省的 37%和45%，是河南省小麥的主產(chǎn)區(qū)［1］。氮是小麥生長(zhǎng)發(fā)育必需的營養(yǎng)元素，近年來，生產(chǎn)上盲目施肥、過度施氮，導(dǎo)致氮肥利用效率低下，不僅沒有達(dá)到預(yù)期增產(chǎn)目標(biāo)，而且造成環(huán)境污染、生產(chǎn)成本提高等問題［2-3］。研究表明，適量氮肥供應(yīng)能夠調(diào)控小麥生長(zhǎng)發(fā)育，改善其光合性能，從而提高小麥產(chǎn)量和品質(zhì)［4］。

明確小麥植株各個(gè)生育時(shí)期的臨界氮濃度是合理施用氮肥的基礎(chǔ)，臨界氮濃度指作物獲得最大生物量增長(zhǎng)所需要的最低氮濃度，即作物獲得最大生物量增長(zhǎng)所需要的最少氮素營養(yǎng)［5-6］。Greenwood等［7］的研究表明，作物地上部含氮量（N）與地上部總生物量（W）之間呈冪函數(shù)（N = aW-b）關(guān)系，并提出了C3、C4作物的臨界氮濃度與地上干物質(zhì)間的通用模型，即N = 5.17W-0.5（C3作物）和N = 4.11W-0.5（C4作物）；Lemaire等［8］利用多個(gè)試驗(yàn)結(jié)果的平均值，將上述模型修正為 N = 4.8DM-0.34（C3作物）和N = 3.6DM-0.34（C4作物），式中DM為地上部生物量。但由于供試作物有限，上述模型仍無法代表所有 C3、C4作物，隨后在玉米［5］、棉花［9］、水稻［10］、番茄［11］等作物上分別進(jìn)行了參數(shù)調(diào)試，并建立了相應(yīng)的模型。在小麥上，國內(nèi)外學(xué)者根據(jù)不同試驗(yàn)、不同品種對(duì)小麥地上部含氮量與地上部總生物量之間的關(guān)系也進(jìn)行了大量的研究，并先后構(gòu)建了小麥臨界氮濃度與地上部生物量之間的關(guān)系模型［12-14］，為小麥氮素營養(yǎng)診斷和精確施肥管理提供了理論基礎(chǔ)。

基于臨界氮濃度，氮營養(yǎng)指數(shù)（nitrogen nutrition index）被定義為實(shí)際氮濃度與臨界氮濃度的比值，可以作為評(píng)價(jià)氮營養(yǎng)狀況的一個(gè)可靠指標(biāo)［15-16］。然而，臨界氮濃度與地上部干物重的冪函數(shù)曲線模型的參數(shù)可能因試驗(yàn)地區(qū)或作物的不同而出現(xiàn)差異［7，12，17］，因此，對(duì)于特定作物來說，研究其地區(qū)適應(yīng)性及不同品種之間的差異顯得十分必要。

本研究選擇兩種類型的小麥主栽品種，通過 3年定位試驗(yàn)，以探討利用小麥臨界氮濃度稀釋曲線來估測(cè)小麥植株氮素盈虧水平和診斷小麥氮素營養(yǎng)的可行性，為小麥氮肥精確管理提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)地概況

河南省禹州市順店鎮(zhèn)康城村（34°27′ N，113°34′E）屬暖溫帶半濕潤(rùn)季風(fēng)季候，土壤類型為潮土，質(zhì)地黏壤，播前耕層土壤（0～20 cm） pH值為8.2，含全氮1.04 g kg-1、有機(jī)質(zhì)20.5 g kg-1、速效鉀142 mg kg-1、速效磷20 mg kg-1；0～30 cm土層中無機(jī)氮（Nmin）為58 kg hm-2，耕層土壤容重為1.45 g cm-3。

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

2011年10月至2014年6月連續(xù)3個(gè)小麥生長(zhǎng)季，選用多穗型品種豫麥49-198和大穗型品種周麥16，設(shè)5個(gè)施氮（尿素，含純氮46%）水平，純氮量分別為0 （N0）、120 （N120）、180 （N180）、240 （N240）和360 kg hm-2（N360），1/2在播前施入，1/2在拔節(jié)期追施；播前隨氮肥同時(shí)施入P2O5和K2O各90 kg hm-2，其中磷肥為過磷酸鈣（含P2O512%），鉀肥為氯化鉀（含K2O 60%）。小區(qū)長(zhǎng)8 m，寬6 m，重復(fù)3次，隨機(jī)區(qū)組排列。兩個(gè)小麥品種的播量均為150 kg hm-2。

于同一地塊定位試驗(yàn)，種植模式為小麥-玉米一年二熟制。2011—2012年度于 2011年 10月 17日播種，10月26日出苗，2012年6月5日收獲；2012 —2013年度于2012年10月12日播種，10月19日出苗，2013年6月5日收獲；2013—2014年度于2013 年10月13日播種，10月22日出苗，2014年6月1日收獲。為保證不同氮水平對(duì)土壤造成的差異一致，玉米季氮、磷、鉀肥總用量與小麥季相同，即也設(shè)同樣5個(gè)氮水平，但1/3氮肥在播前施入，2/3在大喇叭口期追施，而磷、鉀肥全部在播種前一次性施入。

1.3 地上部生物量和植株氮含量的測(cè)定

在小麥播后68 d （越冬期）、132～137 d （返青期）、156～165 d （拔節(jié)期）、196～198 d （開花期）、232～236 d（成熟期），取每小區(qū)1 m雙行植株樣，105℃殺青，75℃烘干后稱取重量。

從上述各生育期所取植株樣品中選擇代表性的10株（開花期和成熟期為10個(gè)單莖），105℃殺青，75℃烘干后粉碎，充分混勻，用凱氏定氮法測(cè)定小麥整株的全氮含量，利用反射儀（RQ-flex）以比色法直接測(cè)定植株莖基部的硝酸鹽含量［18］。

1.4 小麥臨界氮濃度稀釋曲線

綜合 Justes等［12］和薛曉萍等［9］提出的臨界氮濃度稀釋曲線的計(jì)算方法，主要包括4個(gè)步驟∶ （1）對(duì)比分析不同氮素水平試驗(yàn)下每次取樣的地上部干物質(zhì)重及相應(yīng)的植株氮濃度值，用方差分析對(duì)作物生長(zhǎng)受氮素營養(yǎng)限制與否的施氮水平進(jìn)行分類；（2）對(duì)施氮量不能滿足作物最大生長(zhǎng)需求的試驗(yàn)資料，其地上干物質(zhì)與氮濃度值的關(guān)系以線性曲線擬合；（3）對(duì)作物生長(zhǎng)不受氮素影響的施氮水平，用其地上部干物質(zhì)的平均值代表最大干物質(zhì)；（4）每個(gè)取樣日的理論臨界氮濃度由上述線性曲線與橫坐標(biāo)上最大干物質(zhì)垂線的相交點(diǎn)的縱坐標(biāo)決定。小麥臨界氮濃度稀釋曲線通用模型為Nc= aDM-b，式中，Nc為地上部臨界氮濃度值（%），DM 為地上最大干物質(zhì)積累量（t hm-2），a和b為參數(shù)，a為地上部干物質(zhì)為1 Mg hm-2時(shí)的氮濃度，b為控制該曲線斜率的統(tǒng)計(jì)參數(shù)。

1.5 模型的檢驗(yàn)

采用回歸估計(jì)標(biāo)準(zhǔn)誤差（RMSE）分析模擬值和真實(shí)值之間的符合度，RMSE值越小，模擬值與真實(shí)值的一致性越好，偏差越小，即模型的預(yù)測(cè)精度越高；同時(shí)，用模擬值與真實(shí)值之間1∶1直方圖來直觀顯示模型的擬合度和可靠性［19］。

式中Pi為實(shí)測(cè)值，Oi為對(duì)應(yīng)的模擬值，N為樣本量。

本文用2011—2012和2012—2013年度數(shù)據(jù)分別曲線擬合建立模型，將2013—2014年度試驗(yàn)數(shù)據(jù)帶入所建模型中，計(jì)算出2013—2014年度各個(gè)取樣日的臨界氮濃度值，為模擬值；根據(jù)臨界氮濃度稀釋曲線的計(jì)算方法，計(jì)算2013—2014年度各個(gè)取樣日的實(shí)際臨界氮濃度值，為真實(shí)值，利用公式對(duì)模型進(jìn)行檢驗(yàn)。

1.6 小麥氮營養(yǎng)指數(shù)

為了更精確地反映小麥植株含氮量是否適宜，氮營養(yǎng)指數(shù)模型被用來定量評(píng)估小麥體內(nèi)的氮素狀況［20］。NNI = Na/Nc，式中NNI為氮營養(yǎng)指數(shù)，Na為小麥地上部氮濃度的實(shí)測(cè)值，Nc為根據(jù)臨界氮濃度稀釋模型求得的臨界氮濃度值。NNI = 1表明小麥體內(nèi)氮素營養(yǎng)水平處于最佳狀態(tài)，NNI ＞ 1表明氮營養(yǎng)過剩，NNI ＜ 1表明氮營養(yǎng)不足［21］。

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

采用Microsoft Excel 2010和SPSS20.0處理試驗(yàn)數(shù)據(jù)，用 Duncan's法進(jìn)行處理間多重比較；采用OriginPro 8.1繪圖。

2 結(jié)果與分析

2.1 施氮量對(duì)小麥地上部生物量累積的影響

3年試驗(yàn)結(jié)果表明，地上部生物量均隨著施氮量的增加而增加，N180、N240和N360處理間差異不顯著，但都顯著高于 N0和 N120處理（表 1）。根據(jù)前人報(bào)道的臨界氮濃度稀釋曲線模型計(jì)算方法［9，12］，分別對(duì)比分析豫麥49-198和周麥16的地上部生物量，均滿足統(tǒng)計(jì)意義上的不等式DM0＜ DM120＜ DM180= DM240= DM360。

2.2 小麥植株臨界氮濃度稀釋曲線模型的建立

以2011—2012和2012—2013年度數(shù)據(jù)擬合模型，發(fā)現(xiàn)隨地上部生物量的增長(zhǎng)，小麥臨界氮濃度（Nc）呈逐漸下降的趨勢(shì)，兩品種的模型結(jié)構(gòu)相同，但臨界氮稀釋曲線斜率有明顯差異，因此分別對(duì)兩品種進(jìn)行臨界氮濃度稀釋曲線擬合，方程的決定系數(shù)分別為0.8203和0.7981，均達(dá)極顯著水平（圖1）。說明這 2個(gè)方程可以用來表征小麥植株氮含量與地上部生物量的關(guān)系。

兩個(gè)小麥品種臨界氮濃度模型（Nc= aDM-b）中的參數(shù)b有一定的差異，但變化范圍較??；而參數(shù)a的變化范圍較大，在3.756～4.103之間，說明在相同干物質(zhì)情況下，豫麥49-198的臨界氮濃度值高于周麥16。

2.3 臨界氮濃度稀釋曲線模型的驗(yàn)證

利用 2013—2014年度數(shù)據(jù)對(duì)上述模型驗(yàn)證可以看出，兩品種的 5個(gè)取樣期的臨界氮濃度模擬值與真實(shí)值的誤差分別為-0.32%～0.64%和-0.21%～0.60%，按公式（1）求得兩品種臨界氮濃度的 RMSE值分別為0.38和0.35 （表2）。由圖2可知，模擬值與真實(shí)值之間擬合度達(dá)到極顯著水平，二者之間基本符合y = x的1∶1線性關(guān)系，說明模擬值與真實(shí)值之間具有較高的吻合度，模擬精度較高。由此可見，臨界氮濃度稀釋曲線模擬模型在年度間具有較好的穩(wěn)定性，可作為氮營養(yǎng)診斷曲線。

2.4 小麥植株地上部氮素營養(yǎng)指數(shù)模擬模型的建立

為了驗(yàn)證本文建立的小麥臨界氮稀釋曲線的合理性，進(jìn)一步建立了小麥植株地上部氮素營養(yǎng)指數(shù)模擬模型。兩品種的氮營養(yǎng)指數(shù)隨著施氮量的增加而上升，隨小麥生育時(shí)期的推進(jìn)呈波動(dòng)狀態(tài)，且 3年變化趨勢(shì)基本一致（圖3）。N0處理由于沒有施用氮肥，小麥氮營養(yǎng)指數(shù)從返青期到成熟期逐漸下降；而其他處理的氮素營養(yǎng)指數(shù)在拔節(jié)期有明顯的下降趨勢(shì)，表明拔節(jié)期對(duì)氮素需求量比較大，追施氮肥后，氮素營養(yǎng)指數(shù)明顯升高。在小麥整個(gè)生育時(shí)期，N0、N120處理的氮素營養(yǎng)指數(shù)均低于 1，表明植株受氮素的制約；N240、N360處理的氮素營養(yǎng)指數(shù)均大于 1，表明氮素供應(yīng)充足；而 N180處理的氮素營養(yǎng)指數(shù)一直在1附近波動(dòng)，表明施氮量適宜。由此可見，氮營養(yǎng)指數(shù)是一個(gè)較好的小麥植株氮營養(yǎng)狀況診斷指標(biāo)。

2.5 氮營養(yǎng)指數(shù)與植株相對(duì)莖基部硝酸鹽之間的關(guān)系

小麥植株莖基部硝酸鹽含量與植株全氮有密切的關(guān)系，當(dāng)作物輕微缺氮時(shí)，植株全氮還未發(fā)生變化，而植株莖基部硝酸鹽已發(fā)生明顯的變化，當(dāng)供氮超過作物需求時(shí)，莖基部硝酸鹽也會(huì)比全氮有大幅度的提高，通過莖基部硝酸鹽能夠更快、更易判斷冬小麥氮素營養(yǎng)狀況［25］。為了驗(yàn)證本文所求氮營養(yǎng)指數(shù)在評(píng)價(jià)植株氮素營養(yǎng)上的準(zhǔn)確性，分別分析了2個(gè)小麥品種氮營養(yǎng)指數(shù)與相對(duì)莖基部硝酸鹽之間的關(guān)系。圖4和圖5顯示，兩品種的相對(duì)莖基部硝酸鹽均隨著氮素營養(yǎng)指數(shù)的增加而增加，并呈線性關(guān)系，除豫麥49-198越冬期相對(duì)莖基部硝酸鹽與氮素營養(yǎng)指數(shù)方程的決定系數(shù)較小、擬合度不顯著外，其他均達(dá)到極顯著水平?？梢姡獱I養(yǎng)指數(shù)與相對(duì)莖基部硝酸鹽之間關(guān)系密切，利用氮營養(yǎng)指數(shù)來評(píng)價(jià)小麥植株氮素營養(yǎng)狀況的方法是準(zhǔn)確的。

表1 不同施氮量對(duì)小麥地上部生物量累積的影響Table 1 Effect of nitrogen application rate on shoot biomass accumulation in wheat （t hm-2）

圖1 小麥地上部生物量臨界氮濃度稀釋曲線Fig. 1 Critical nitrogen concentration dilution curve in shoot biomass of wheat

表2 小麥臨界氮濃度（CNC）模擬值與真實(shí)值Table 2 Simulated and observed critical nitrogen concentrations （CNC） of wheat

圖2 小麥臨界氮濃度模擬值與真實(shí)值的關(guān)系Fig. 2 Relationship between simulated and observed critical nitrogen concentrations of wheat

圖3 不同氮肥處理對(duì)小麥不同生育時(shí)期植株氮營養(yǎng)指數(shù)（NNI）的影響Fig. 3 Effects of different nutrition fertilizer treatments on nitrogen nutrition index （NNI） in different growth stages of wheat

3 討論

為解決目前小麥生產(chǎn)上過量施氮的普遍問題，亟需加強(qiáng)作物的氮素營養(yǎng)診斷研究，建立快速有效的技術(shù)方法，而明確小麥主要生育階段的臨界氮濃度是科學(xué)診斷植株氮營養(yǎng)狀況，指導(dǎo)合理施氮的基礎(chǔ)［22］。研究表明，作物在生長(zhǎng)過程中體內(nèi)氮濃度與其生物量的關(guān)系符合冪函數(shù)關(guān)系，但作物種類和品種、種植環(huán)境等均影響模型參數(shù)，因此國內(nèi)外眾多學(xué)者先后建立和修正了多個(gè)臨界氮濃度稀釋模型。如 Justes等［12］構(gòu)建的高蛋白小麥 Nc模型（Nc= 5.35DM-0.44），薛曉萍等［9］在安陽和南京構(gòu)建的 2個(gè)棉花Nc模型（Nc= 3.387DM-0.131和Nc= 2.858DM-0.131），梁效貴等［5］構(gòu)建的華北夏玉米 Nc模型（Nc= 34.914DM-0.4134），Colnenne等［23］建立的油菜Nc模型（Nc= 4.48DM-0.25），Tei等［24］構(gòu)建的馬鈴薯Nc模型（Nc= 3.90DM-0.27），這些研究所得到的模型形式與 Greenwood等提出的模型一致，但參數(shù)不同。在小麥上，與 Justes等［12］的模型相比，Yue等［13］建立的中國華北平原小麥Nc模型（Nc= 4.15DM-0.38）和趙犇等［14］利用揚(yáng)麥16和寧麥13構(gòu)建的模型（Nc= 4.65DM-0.44和 Nc= 4.33DM-0.45），其參數(shù)相差較大；而本研究針對(duì)豫麥49-198和周麥16擬合的兩個(gè)方程（Nc= 4.103DM-0.435和Nc= 3.756DM-0.388）與Yue等［13］和趙犇等［14］的模型參數(shù)差異相對(duì)較小。

圖4 豫麥49-198氮營養(yǎng)指數(shù)（NNI）與相對(duì)莖基部硝酸鹽的關(guān)系Fig. 4 Relationship between nitrogen nutrition index （NNI） and nitrate concentration in stem base of Yumai 49-198

圖5 周麥16氮營養(yǎng)指數(shù)（NNI）與相對(duì)莖基部硝酸鹽的關(guān)系Fig. 5 Relationship between nitrogen nutrition index （NNI） and nitrate concentration in stem base of Zhoumai 16

本研究利用 3年5個(gè)氮水平的數(shù)據(jù)，構(gòu)建豫中地區(qū)兩個(gè)主要小麥品種臨界氮濃度稀釋模型，通過回歸估計(jì)標(biāo)準(zhǔn)誤差及模擬值與真實(shí)值之間 1︰1的直方圖驗(yàn)證可知，這兩個(gè)方程在不同年度間具有較好的穩(wěn)定性，可以作為豫中地區(qū)小麥氮營養(yǎng)診斷曲線。這兩個(gè)方程的參數(shù)a差異較大，也與Justes等［12］模型中的參數(shù) a有很大差異，其主要原因可能是品種因素。本研究所采用的兩個(gè)小麥品種，代表了生產(chǎn)中多穗型和大穗型兩個(gè)類型的小麥品種，是我們前期從常用10個(gè)品種對(duì)比試驗(yàn)中篩選出來的。不同類型小麥品種對(duì)氮素營養(yǎng)的反應(yīng)不同［18，25-26］，多穗型豫麥49-198的植株吸氮能力大于周麥16。據(jù)報(bào)道，臨界氮稀釋曲線的斜率（參數(shù)b）受氣候、土壤、密度等環(huán)境條件的影響較小［9，14］，而本文所構(gòu)建的2個(gè)方程其參數(shù) b有明顯差異，否是與這 2個(gè)品種的籽粒蛋白質(zhì)含量差異有關(guān)，還需進(jìn)一步研究。

Lemaire等［16］基于作物臨界氮濃度稀釋模型提出氮營養(yǎng)指數(shù)的概念，可定量動(dòng)態(tài)描述作物的氮營養(yǎng)狀況。趙犇等［14］研究表明，小麥氮營養(yǎng)指數(shù)在0.37～1.28之間，并隨施氮量的增加而升高。在本研究中，2011—2012和 2012—2013年度該指標(biāo)介于0.48～1.33之間，而2013—2014年度則在0.51～1.69之間，特別是小麥開花到成熟期間，氮營養(yǎng)指數(shù)遠(yuǎn)遠(yuǎn)高于前兩年數(shù)據(jù)，這可能與2013—2014年度氣溫較常年偏高，降水量大幅減少有關(guān)。目前還沒有關(guān)于氣溫、降水量等氣象因素對(duì)氮營養(yǎng)指數(shù)影響的報(bào)道，本課題組擬對(duì)此系統(tǒng)研究。

本研究表明，氮營養(yǎng)指數(shù)在拔節(jié)期迅速降低，說明此時(shí)植株對(duì)氮需求量較大，導(dǎo)致氮吸收不足和氮濃度降低；N180處理的氮營養(yǎng)指數(shù)在 1附近波動(dòng)，反映了小麥的適宜施氮量是180 kg hm-2，這與本課題組多年多點(diǎn)的研究結(jié)果一致［18，25-26］。由此也看出，用氮營養(yǎng)指數(shù)來評(píng)價(jià)植株氮營養(yǎng)狀況的方法是可靠的，可以用來診斷調(diào)控小麥氮素營養(yǎng)。小麥植株莖基部硝酸鹽含量與植株全氮關(guān)系密切，利用莖基部硝酸鹽含量能夠更快、更方便地判斷小麥氮素營養(yǎng)狀況［27］。本研究發(fā)現(xiàn)，氮營養(yǎng)指數(shù)與小麥相對(duì)莖基部硝酸鹽顯著相關(guān)，因此認(rèn)為氮營養(yǎng)指數(shù)模型可以用來評(píng)價(jià)小麥植株氮素盈虧水平，建議采用小麥臨界氮濃度稀釋曲線診斷小麥氮營養(yǎng)。

4 結(jié)論

建立了兩種穗型小麥品種的地上部臨界氮濃度稀釋曲線模型，該模型具有較好的穩(wěn)定性，可以說明豫中地區(qū)植株地上部生物量和氮濃度之間的關(guān)系。進(jìn)一步利用臨界氮濃度建立了小麥植株地上部氮素營養(yǎng)指數(shù)的模擬模型，可以較好地評(píng)價(jià)小麥植株氮營養(yǎng)狀況。根據(jù)本研究建立的氮素營養(yǎng)指數(shù)模型，推薦豫中地區(qū)冬小麥最佳施氮量為180 kg hm-2，在該施氮水平可獲得最高產(chǎn)量。

References

［1］ 李春喜，余海波，姜麗娜，趙艷玲，邵云，張黛靜. 耕作施肥模式對(duì)豫中冬小麥干物質(zhì)積累的影響. 河南農(nóng)業(yè)科學(xué)，2012，41（9）∶ 26-30 Li C X，Yu H B，Jiang L N，Zhao Y L，Shao Y，Zhang D J. Effect of tillage and fertilization models on dry matter accumulation of winter wheat in mid Henan province. Henan Agric Sci，2012，41（9）∶ 26-30 （in Chinese with English abstract）

［2］ 李歡歡，黃玉芳，王玲敏，張立花，葉優(yōu)良. 河南省小麥生產(chǎn)與肥料施用狀況. 中國農(nóng)學(xué)通報(bào)，2009，25（18）∶ 426-430 Li H H，Huang Y F，Wang L M，Zhang L H，Ye Y L. Wheat production and fertilizer application in Henan province. Chin Agric Sci Bull，2009，25（18）∶ 426-430 （in Chinese with English abstract）

［3］ 張福鎖，王激清，張衛(wèi)峰，崔振嶺，馬文奇，陳新平，江榮風(fēng).中國主要糧食作物肥料利用率現(xiàn)狀與提高途徑. 土壤學(xué)報(bào)，2008，4∶ 915-922 Zhang F S，Wang J Q，Zhang W F，Cui Z L，Ma W Q，Chen X P，Jiang R F. China's main food crop fertilizer utilization present situation and way to improve. Acta Pedol Sin，2008，4∶ 915-922 （in Chinese with English abstract）

［4］ 陸增根，戴廷波，姜東，荊奇，吳正貴，周培南，曹衛(wèi)星. 氮肥運(yùn)籌對(duì)弱筋小麥群體指標(biāo)與產(chǎn)量和品質(zhì)形成的影響. 作物學(xué)報(bào)，2007，33∶ 590-597 Lu Z G，Dai T B，Jiang D，Jing Q，Wu Z G，Zhou P N，Cao W X. Nitrogen management on groups of weak gluten wheat yield and quality indicators and the influence of the formation. Acta Agron Sin，2007，33∶ 590-597 （in Chinese with English abstract）

［5］ Plenet D，Lemaireet G. Relationships between dynamics of nitrogen uptake and dry matteraccumulation in maize crops∶ determination of critical N concentration. Plant & Soil，2000，216∶ 65-82

［6］ Ulrich A. Physiological bases for assessing the nutritional requirements of plants. Annu Rev Plant Physiol，1992，3∶ 207-228

［7］ Greenwood D J，Lemaire G，Grosse G，Cruz P，Draycott A，Neeteson J J. Decline in percentage N of C3and C4crops with increasing plant mass. Ann Bot，1990∶ 181-190

［8］ Lemaire G，Gastal F，Cruz P，Greenwood D J，Draycott A. Relationships between plant-N，plant mass and relative growth rate for C3and C4crops. In∶ Scaife A ed. Proceedings of the First European Society of Agronomy Congress，Paris，1990. pp 1-5

［9］ 薛曉萍，周治國，張麗娟，王以琳，郭文琦，陳兵林. 棉花花后臨界氮濃度稀釋模型的建立及在施氮量調(diào)控中的應(yīng)用. 生態(tài)學(xué)報(bào)，2006，26∶ 1781-1790Xue X P，Zhou Z G，Zhang L J，Wang Y L，Guo W Q，Chen B L. Cotton flower after critical nitrogen concentration dilution model and application in the application rate control. Acta Ecol Sin，2006，26∶ 1781-1790 （in Chinese with English abstract）

［10］ 劉小軍，Syed Tahir，Ata-Ul-Karim，陸震洲，曹衛(wèi)星，朱艷. 水稻臨界氮濃度稀釋模型及氮素營養(yǎng)診斷研究. 中國稻米，2013，19（4）∶ 149 Liu X J，Syed Tahir，Ata-Ul-Karim，Lu Z Z，Cao W X，Zhu Y. Critical nitrogen concentration diluted rice model and nitrogen nutrition diagnosis. China Rice，2013，19（4）∶ 149 （in Chinese with English abstract）

［11］ 王新，馬富裕，刁明，樊華，崔靜，賈彪，何海兵，劉其. 滴灌番茄臨界氮濃度、氮素吸收和氮營養(yǎng)指數(shù)模擬. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（18）∶ 99-108 Wang X，Ma F Y，Diao M，F(xiàn)an H，Cui J，Jia B，He H B，Liu Q. Drip irrigation tomato critical nitrogen concentration，nitrogen absorption and nitrogen nutrition index simulation. Trans CASE，2013，29（18）∶ 99-108 （in Chinese with English abstract）

［12］ Justes E，Mary B，Meynard J M，Machet J M，Thelier-Huche L. Determination of a critical nitrogen dilution curve for winter wheat crops. Ann Bot，1994，4∶ 397-407

［13］ Yue S C，Meng Q F，Zhao R F，Li F，Chen X P，Zhang F S. Critical Nitrogen dilution curve for optimizing nitrogen management of winter wheat production in the North China Plain. Agron J，2012，104∶ 523-528

［14］ 犇趙，姚霞，田永超，劉小軍，曹衛(wèi)星，朱艷. 基于臨界氮濃度的小麥地上部氮虧缺模型. 應(yīng)用生態(tài)學(xué)報(bào)，2012，23∶3141-3148 Zhao B，Yao X，Tian Y C，Liu X J，Cao W X，Zhu Y. Based on the critical nitrogen concentration of wheat aboveground nitrogen deficiency model. Chin Appl Ecol，2012，23∶ 3141-3148 （in Chinese with English abstract）

［15］ Ziadi N，Brassard M，Bélanger G，Cambouris A N，Tremblay N，Nolin M C，Claessens A，Parent L E. Critical N curve and N nutrition index for corn in eastern Canada. Agron J，2008，100∶ 271-276

［16］ Lemaire G，Jeuffroy M H，Gastal F. Diagnosis tool for plant and crop N status in vegetative stage theory and practices for crop N management. Eur J Agron，2008，28∶ 614-624

［17］ Bélanger G，Walsh J R，Richards J E，Milburn P H，Ziadi N. Critical nitrogen curve and nitrogen nutrition index for potato in eastern Canada. Am J Potato Res，2001，78∶ 355-364

［18］ 劉小寧，劉海坤，黃玉芳，葉優(yōu)良. 施氮量、土壤和植株氮濃度與小麥赤霉病的關(guān)系. 植物營養(yǎng)與肥料學(xué)報(bào)，2015，21∶306-317 Liu X N，Liu H K，Huang Y F，Ye Y L. Nitrogen application rate of soil and plant nitrogen and the relationship between wheat gibberellic disease. J Plant Nutr Fert，2015，21∶ 306-317 （in Chinese with English abstract）

［19］ 劉其. 基于過程的滴灌春小麥臨界氮需求量定量化模擬模型.石河子大學(xué)碩士學(xué)位論文，新疆石河子，2013. pp 25-31 Liu Q. Quantitative Simulation Model of Critical Nitrogen Requirement of Spring Wheat Based on Process. MS Thesis of Shihezi University，Shihezi，China，2013. pp 25-31 （in Chinese with English abstract）

［20］ Lemaire G，Gastal F，Salette J. N Uptake and Distribution in Plant Canopies∶ Diagnosis of the Nitrogen Status in Crops. Berlin∶Springer，1997

［21］ Naud C，Mkowski D，Jeuffroy M H. Is it useful to combine measurements taken during growing season with dynamic model to predict the nitrogen status of winter wheat. Eur Agron，2008，28∶ 287-300

［22］ 梁效貴，張經(jīng)廷，周麗麗，李旭輝，周順利. 華北地區(qū)夏玉米臨界氮稀釋曲線和氮營養(yǎng)指數(shù)研究. 作物學(xué)報(bào)，2013，39∶292-299 Liang X G，Zhang J T，Zhou L L，Li X H，Zhou S L. Summer maize in north China the critical nitrogen dilution curve and nitrogen nutrition index. Acta Agron Sin，2013，39∶ 292-299 （in Chinese with English abstract）

［23］ Colnenne C，Meynard J M，Reau R，Justes E，Merrien A. Determination of a critical nitrogen dilution curve for winter oilseed rape. Ann Bot，1988，81∶ 311-317

［24］ Tei F，Benincasa P，Guiducci M. Critical nitrogen concentration in processing tomato. Eur J Agron，2002，18∶ 45-55

［25］ 王桂良. 超高產(chǎn)小麥干物質(zhì)累積、氮素營養(yǎng)及診斷研究. 河南農(nóng)業(yè)大學(xué)碩士學(xué)位論文，河南鄭州，2009. pp 41-42 Wang G L. Super-High-Yield Wheat Dry Matter Accumulation and Nitrogen Nutrition Research and Diagnosis. MS Thesis of Henan Agricultural University，Zhengzhou，China，2009. pp 41-42 （in Chinese with English abstract）

［26］ 王玲敏. 氮素調(diào)控對(duì)高產(chǎn)小麥干物質(zhì)轉(zhuǎn)移、氮素吸收及土壤硝態(tài)氮累積的影響. 河南農(nóng)業(yè)大學(xué)碩士學(xué)位論文，河南鄭州，2012. pp 37-38 Wang L M. Nitrogen Regulation on High-Yield Wheat Transfer of Dry Matter，Nitrogen Absorption and the Influence of Soil NO3Accumulation. MS Thesis of Henan Agricultural University，Zhengzhou，China，2012. pp 37-38 （in Chinese with English abstract）

［27］ 賈良良，陳新平，張福鎖. 葉綠素儀與植株硝酸鹽濃度測(cè)試對(duì)冬小麥氮素營養(yǎng)診斷準(zhǔn)確性的比較研究. 華北農(nóng)學(xué)報(bào)，2007，22（6）∶ 157-160 Jia L L，Chen X P，Zhang F S. Comparative study on the accuracy of nitrogen nutrition diagnosis of winter wheat by the test of the cholophyll meter and the nitrate concentration of the plant. Acta Agric Boreali-Sin，2007，22（6）∶ 157-160 （in Chinese with English abstract）

Simulating Critical Nitrogen Dilution Curve and Modeling Nitrogen Nutrition Index in Winter Wheat in Central Henan Area

YUE Song-Hua，LIU Chun-Yu，HUANG Yu-Fang，and YE You-Liang*
College of Resources and Environment，Henan Agricultural University，Zhengzhou 450002，China

Abstract：To test the applicability of critical nitrogen dilution curve for estimating nitrogen （N） nutrient status of winter wheat in central Henan area，we carried out a site-fixed experiment from October 2011 to June 2014 using the large-spike variety Zhoumai 16 and the multi-spike variety Yumai 49-198. Five N application levels，namely pure N 0，120，180，240，and 360 kg ha-1，were designed to build the critical N concentration dilution curve of shoot biomass and N nutrition index model. The aboveground biomass increased obviously with rising N level and then tended to a mild level. The relationship between critical N concentration and shoot biomass of wheat fit the power function with R2of 0.8203 for Yumai 49-198 and 0.7981 for Zhoumai 16 （P ＜ 0.01）. These models showed good stability across years. The N nutrition index estimated by the model could be used for evaluating the N nutrition situation of wheat plants，and a significant linear correlation was observed between the N nutrition index and the nitrate concentration in wheat stem base. These results indicate that the critical N concentration dilution model is applicable in diagnosing N nutrient status of wheat plants in central Henan area.

Keywords：Wheat；Critical nitrogen dilution；Nitrogen nutrition index；Stem sap nitrate

DOI：10.3724/SP.J.1006.2016.00909

*通訊作者（

Corresponding author）∶ 葉優(yōu)良，E-mail∶ ylye2004@163.com，Tel∶ 0371-63558290

收稿日期Received（）∶ 2015-09-06；Accepted（接受日期）∶ 2016-03-14；Published online（網(wǎng)絡(luò)出版日期）∶ 2016-03-21.