不同水分條件下葡萄臨界氮稀釋曲線模型的建立及氮素營養(yǎng)診斷*

李佳帥，楊再強,2**，李永秀,2，江夢圓，趙和麗，韋婷婷，張旭然

不同水分條件下葡萄臨界氮稀釋曲線模型的建立及氮素營養(yǎng)診斷*

李佳帥1，楊再強1,2**，李永秀1,2，江夢圓1，趙和麗1，韋婷婷1，張旭然1

（1．南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，南京 210044；2．江蘇省農(nóng)業(yè)氣象重點實驗室，南京 210044)

以1a生葡萄植株“紅提”為試材，在Venlo型試驗溫室內(nèi)進(jìn)行土壤水分和施氮量雙因素區(qū)組試驗。試驗設(shè)置4個灌水水平，分別為正常灌溉量W1（田間最大持水量的70%～80%）、輕度水分脅迫W2（60%～70%）、中度水分脅迫W3（50%～60%）和重度水分脅迫W4（30%～40%）；設(shè)置4個施氮水平，分別為1.5倍推薦施氮量（N1，25.5g×plant?1）、正常推薦施氮量（N2，17g×plant?1）、0.5倍推薦施氮量（N3，8.5g×plant?1）和不施用氮肥（N4，0g×plant?1）。每10d觀測一次植株體內(nèi)氮濃度和植株地上部生物量，利用不同水分條件下葡萄植株在一定生長時期內(nèi)所獲最大生物量時對應(yīng)的最小氮濃度值即臨界氮濃度（Nc）構(gòu)建葡萄臨界氮濃度稀釋曲線模型，并在此基礎(chǔ)上建立氮素吸收模型（Nupt）和氮素營養(yǎng)指數(shù)模型（NNI），對不同水分條件下葡萄氮營養(yǎng)狀況進(jìn)行定量診斷。結(jié)果表明：設(shè)施葡萄植株臨界氮濃度與地上部生物量存在冪函數(shù)關(guān)系，隨著灌水量的增加，葡萄植株臨界氮濃度值增大，氮素吸收量及地上部生物量也呈增加趨勢；在W1、W2水分條件下，葡萄植株生物量隨施氮量增加而增加，而W3和W4處理葡萄生物量隨施氮量增加呈先增后降的趨勢；在相同水分條件下，氮濃度隨施氮量增加而增加，隨葡萄生長進(jìn)程而降低；利用Nupt和NNI模型可對植株體內(nèi)氮營養(yǎng)元素虧缺與否進(jìn)行有效診斷。

葡萄；臨界氮濃度；氮素營養(yǎng)指數(shù)；氮素吸收；干物質(zhì)積累

葡萄（）屬多年生落葉藤本植物，在全世界果品生產(chǎn)中，其總產(chǎn)量及栽培面積均居首位。設(shè)施葡萄栽培尤以其經(jīng)濟(jì)效益高，發(fā)展前景廣闊，隨著市場需求逐漸加大，其栽培面積也逐漸增加。截至2014年，中國設(shè)施葡萄栽培面積達(dá)13.3萬hm2，居世界第一[1]。江蘇省葡萄種植面積近4萬hm2，以鮮食葡萄為主[2]。目前，施用過量氮肥以獲取更高產(chǎn)量的做法不僅浪費資源、增加成本，還會造成對土壤環(huán)境的破壞，給持續(xù)生產(chǎn)帶來嚴(yán)重威脅。因此，判斷氮肥施用量是否合理尤為重要，診斷作物體內(nèi)氮素盈虧狀況的基本方法之一，是確定作物的臨界氮濃度值[3]，臨界氮濃度為在一定生長時期內(nèi)獲得最大生物量時的最小氮濃度[4]，其對實時了解作物氮營養(yǎng)狀況，提高作物品質(zhì)及產(chǎn)量，避免肥料浪費具有重要意義。

目前，許多國內(nèi)外學(xué)者已對部分作物的氮素營養(yǎng)做出了合理診斷。Lemaire等[5]于1984年首次提出牧草的臨界氮濃度稀釋曲線模型，表明牧草在不受氮素限制時，地上部干物質(zhì)量與氮濃度存在冪函數(shù)關(guān)系，土壤含氮量成為作物干物質(zhì)量累積的主要限制因素之一。隨后Justes等[6]在1994年建立了冬小麥臨界氮濃度稀釋曲線模型Nc=5.35×DM?0.442，Sheehy等[7]對熱帶水稻的研究結(jié)果為Nc=5.35×DMM?0.5，隨后該模型在冬油菜[8]、玉米[9]、向日葵[10]、番茄[11]、包心菜[12]、馬鈴薯[13]、高粱[14]、甜椒[15]等作物上均有應(yīng)用。這些研究所得到的模型形式一致，但參數(shù)大多不同。國內(nèi)也有許多學(xué)者對作物臨界氮濃度稀釋曲線進(jìn)行研究，薛曉萍等[16]構(gòu)建了不同地區(qū)棉花的臨界、最低和最高氮濃度稀釋曲線模型，并得出該作物的最佳施氮量；梁效貴等[17?18]構(gòu)建了不同地區(qū)夏玉米的臨界氮稀釋曲線模型并在此基礎(chǔ)上建立了氮營養(yǎng)指數(shù)模型，氮營養(yǎng)指數(shù)可準(zhǔn)確診斷作物體內(nèi)氮素狀況。向友珍等[19]對甜椒的研究發(fā)現(xiàn)，隨著土壤水分提高，植株的氮素吸收量和地上部干物質(zhì)量呈先增大后減少的變化趨勢，各水分條件下最佳施氮量均介于150～225kg·hm?2。楊慧等[20]構(gòu)建了番茄植株的臨界和邊界氮濃度稀釋曲線模型，并發(fā)現(xiàn)在相同灌溉水平下，植株氮濃度隨施氮量增加而增加，隨生長進(jìn)程而降低，隨土壤水分增加植株體內(nèi)氮濃度整體提高。

前人研究表明，臨界氮濃度稀釋曲線模型、氮營養(yǎng)指數(shù)模型和氮素吸收模型可診斷植株氮營養(yǎng)狀況，明確植株需氮量。植株體內(nèi)氮濃度受水分供應(yīng)影響，適宜水分會促進(jìn)干物質(zhì)累積，當(dāng)?shù)?yīng)不足時，會稀釋植株體內(nèi)氮濃度。而前人在構(gòu)建臨界氮濃度稀釋曲線模型時僅考慮氮素因素，關(guān)于水分氮素雙因素對該模型的影響鮮有報道。因此，本研究擬通過水氮耦合的全面試驗，建立設(shè)施葡萄植株臨界氮濃度稀釋曲線模型、氮素營養(yǎng)指數(shù)模型和氮素吸收模型，并檢驗其可靠性，探究利用此模型來指導(dǎo)設(shè)施葡萄氮素營養(yǎng)診斷，以期為設(shè)施葡萄的氮素合理利用提供理論依據(jù)和技術(shù)支持。

1 材料與方法

1.1 試驗地概況

試驗于2018年3?8月在地處亞熱帶季風(fēng)氣候區(qū)的南京信息工程大學(xué)農(nóng)試站Venlo型試驗溫室（32.20°N，118.71°E）內(nèi)進(jìn)行。該地海拔28m，年平均氣溫為15.6℃。試驗期間正值春夏季，溫室內(nèi)整體日平均氣溫呈上升趨勢，平均空氣相對濕度為74.5%，平均日光合有效輻射累計值為3.22MJ·m?2。溫室頂高5.0m，肩高4.5m，寬9.6m，長30.0m，南北走向，采用自動天窗和側(cè)通風(fēng)口通風(fēng)。供試土壤為中壤土，土壤肥力均一，pH為7.4，有機質(zhì)含量為18.32g·kg?1，全氮0.86g·kg?1，全磷0.75g·kg?1，土壤體積含水量32.45%。

1.2 試驗設(shè)計

以1a生“紅提”苗木為供試葡萄品種進(jìn)行盆栽試驗，盆的上口直徑41.5cm，底部直徑26.5cm，深31cm。設(shè)置水和氮兩個因素，每個因素均四個水平，四個水分梯度[21]分別：正常水分W1（田間最大持水量的70%～80%）、輕度脅迫W2（田間最大持水量的60%～70%）、中度脅迫W3（田間最大持水量的50%～60%）、重度脅迫W4（田間最大持水量的30%～40%）。通過EM50（Decagon Device，USA）采集土壤深度25cm處水分?jǐn)?shù)據(jù)，獲取每小時的土壤水分?jǐn)?shù)據(jù)。查閱文獻(xiàn)并計算出葡萄推薦施氮量，葡萄生長初期的推薦施氮量[22?23]為純氮17g×plant?1，因此，設(shè)置四個氮素指標(biāo)分別為：1.5倍推薦施肥（N1）、正常推薦施肥（N2）、0.5倍推薦施肥（N3）、不施用氮肥（N4）。試驗共16個處理，設(shè)置如表1，每個處理重復(fù)6次，其中3次的數(shù)據(jù)均值用于建模，另外3次的數(shù)據(jù)均值用于檢驗。采集葡萄植株鮮樣均控制在9：00?11：00進(jìn)行，試驗期間共取樣9次。當(dāng)日樣品于液氮速凍，并收集于?20℃冰箱中備用，采樣全部結(jié)束后，統(tǒng)一測定相關(guān)指標(biāo)。試驗期間，除控制變量外，其它管理措施正常。

表1 葡萄植株水氮處理設(shè)置

1.3 項目測定

地上部生物量：處理后每10d每個處理選6株，試驗期間共采樣9次。將葡萄植株莖和葉用水洗凈擦干后，在105℃殺青30min，再于85℃恒溫烘干至恒重，用精度為0.001g電子天平分別測定植株地上部干物質(zhì)量。

植株氮濃度：將烘干樣品粉碎后用H2SO4?H2O2消化，最后用凱氏法蒸餾測定全氮濃度[24]。

各器官氮累積量[19]=器官氮濃度×器官干物質(zhì)量 （1）

植株地上部氮累積量[19]=莖氮累積量+葉氮累積量 （2）

1.4 模型描述

1.4.1 臨界氮濃度稀釋曲線模型

臨界氮濃度（NC）指植物在一定生長時期內(nèi)獲得最大生物量時對應(yīng)的最小氮濃度[4]。在植株生長發(fā)育過程中，若地上部氮濃度處于臨界值以下，或增加施氮量時，植株干物質(zhì)量顯著增加，說明作物的生長將受制于氮元素；反之，處于臨界值以上，或增加施氮量植株地上部干物質(zhì)量并無明顯變化，說明作物的生長不受制于氮元素，因此，地上部氮濃度值與臨界氮濃度相同時的施氮量最適宜。

構(gòu)建基于植株生物量的臨界氮濃度稀釋曲線模型，首先需確定臨界氮濃度值，構(gòu)建臨界氮稀釋曲線模型的步驟[6，16]如下：

（1）通過方差分析將不同時段觀測的植株地上部干物質(zhì)量和對應(yīng)地上部氮濃度數(shù)據(jù)，按植株生長發(fā)育受氮元素限制與否進(jìn)行分類；

（2）將植株生長受限于氮素的地上部干物質(zhì)量（橫坐標(biāo)）與氮濃度（縱坐標(biāo)）間的關(guān)系進(jìn)行曲線擬合；

（3）當(dāng)植株生長不受限于氮素影響時，干物質(zhì)量的最大值用其地上部干物質(zhì)量平均值來代表；

（4）每個理論臨界氮濃度值，即在最大干物質(zhì)量處作橫坐標(biāo)的垂線，垂線與上述擬合曲線交點對應(yīng)的縱坐標(biāo)值；

（5）曲線擬合：利用回歸擬合方程法以臨界氮濃度值Nc（g·kg?1）為因變量，以植株地上部最大干物質(zhì)量DMmax（t·hm?2）為自變量，建立臨界氮濃度稀釋曲線模型，即

式中，a為植株地上部干物質(zhì)量為1t·hm?2時的臨界氮濃度，b為決定臨界氮濃度稀釋曲線斜率的參數(shù)[25]。

1.4.2 氮素營養(yǎng)指數(shù)模型

為更精確地反映設(shè)施葡萄植株氮營養(yǎng)情況，采用氮營養(yǎng)指數(shù)模型來定量評估葡萄植株體內(nèi)氮素情況[5]。即

式中，Nt為葡萄植株氮濃度實測值（g·kg?1），Nc為模擬值（g·kg?1）。NNI=1表明植株體內(nèi)氮營養(yǎng)最適，NNI>1表明氮營養(yǎng)過多，NNI<1表明氮營養(yǎng)缺乏。

1.4.3 氮素吸收模型

植株氮吸收量（Nupt，kg·hm?2）與地上部最大干物質(zhì)量（DMmax，t·hm?2）之間的關(guān)系[16]為

將式（1）代入式（3），得到葡萄植株臨界氮吸收模型為

1.5 模型的檢驗

采用均方根誤差（RMSE）分析模擬值與實測值的擬合度，RMSE越小，說明誤差越小，模擬值擬合效果越好，即模型的模擬精度越高；同時，用模擬值與實測值之間1:1圖來顯示模型的精度[26]。

式中，Pi為實測值，Oi為模擬值，N為樣本數(shù)量。

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

采用Microsoft Excel 2010進(jìn)行數(shù)據(jù)處理和作圖，采用SPSS20.0軟件進(jìn)行數(shù)據(jù)的方差分析。

2 結(jié)果與分析

2.1 不區(qū)分水分條件下氮濃度稀釋曲線

所有處理每10d觀測數(shù)據(jù)顯示，試驗過程中，設(shè)施葡萄植株地上部生物量的變化范圍在0.10～11.29t·hm?2，氮濃度的變化范圍在0.52～8.61g·kg?1，不同階段兩者關(guān)系見圖1。由圖1a可見，當(dāng)?shù)厣喜可锪啃∮?t·hm?2（取樣數(shù)n=42）時，生物量和氮濃度的數(shù)據(jù)波動較大，二者無顯著關(guān)系。當(dāng)?shù)厣喜可锪看笥?t·hm?2時（圖1b，n=102），氮濃度隨生物量增加呈下降趨勢即存在氮稀釋現(xiàn)象，二者的關(guān)系可表示為Nc=3.66DW?0.46，方程的決定系數(shù)為0.0583，說明當(dāng)?shù)厣喜可锪看笥?t·hm?2時，氮濃度與地上部生物量存在冪函數(shù)關(guān)系。

2.2 不同水分條件下臨界氮稀釋曲線模型的構(gòu)建與驗證

圖1 葡萄植株氮濃度與地上部生物量的關(guān)系

注：圖a表示地上部生物量小于1t·hm?2（n=42），圖b表示地上部生物量大于1t·hm?2（n=102）。

Note: Fig.a indicates that the aboveground biomass is less than 1t·ha?1(n=42). Fig. b indicates that the aboveground biomass is more than 1t·ha?1(n=102).

圖2 葡萄植株地上部生物量小于1t·hm?2（1）和大于1t·hm?2（2）時不同水分條件下氮濃度與地上部生物量的關(guān)系

注：W1、W2、W3和W4分別表示田間最大持水量的70%～80%、60%～70%、50%～60%和30%～40%水分處理。下同。

Note : W1, W2, W3 and W4 is 70%?80%, 60%?70%, 50%?60% and 30%?40% of maximum water capacity in the field, respectively. The same as below.

2.2.1 模型構(gòu)建

當(dāng)生物量小于1.00t·hm?2時，設(shè)施葡萄植株地上部氮濃度值隨生物量增加無明顯下降趨勢，因此，該階段的氮濃度值不進(jìn)行模擬。當(dāng)?shù)厣喜可锪看笥?.00t·hm?2時，由于各試驗地土壤濕度不同，故將數(shù)據(jù)分為4組：正常灌溉W1（田間最大持水量的70%～80%）、輕度脅迫W2（60%～70%）、中度脅迫W3（50%～60%）、重度脅迫W4（30%～40%）。根據(jù)式（1）對植株氮濃度及對應(yīng)的地上部生物量數(shù)據(jù)進(jìn)行擬合，得到4個水分條件下設(shè)施葡萄植株的臨界氮濃度稀釋曲線模型，并分別選取最大、最小氮濃度實測值模擬得出氮濃度稀釋邊界模型（圖2），且由式（3）和式（4）得出臨界氮吸收模型，兩個模型參數(shù)見表2。分析實測數(shù)據(jù)得出，不同水氮條件下葡萄植株地上部干物質(zhì)積累過程存在顯著差異，總體上表現(xiàn)為：土壤濕度越小，植株干物質(zhì)積累越少。正常供水情況下（W1）植株地上部干物質(zhì)累積整體水平較高（0.14～11.29t·hm?2）；隨著土壤濕度的降低，W2和W3處理中植株干物質(zhì)量分別降至0.12～10.92t·hm?2、0.10～10.42t·hm?2水平；在嚴(yán)重缺水情況下（W4），植株干物質(zhì)量僅0.12～7.62t·hm?2，顯著低于其它灌水處理（P<0.05），干物質(zhì)積累量小于1.00t·hm?2的情況也逐漸增多。同時，土壤含水量對葡萄植株氮濃度的影響也較大，土壤中水分越少植株吸收的氮也越少，且各水分處理間差異顯著。正常水分條件下（W1），地上部植株體內(nèi)氮濃度整體水平較高（1.23～8.61g·kg?1）；隨著土壤濕度的降低，W2、W3處理中植株氮濃度分別降至1.16～7.90g·kg?1、0.52～6.98g·kg?1水平；在嚴(yán)重缺水情況下（W4），植株氮濃度僅0.59～3.50g·kg?1，顯著低于其它灌水處理（P<0.05）。同時，在正常供水情況下，隨著地上部生物量的增加，植株氮濃度有降低的趨勢，但在相對缺水條件下，這種降低趨勢逐漸減弱，植株氮濃度均處于較低水平。說明在水分供應(yīng)充足時植株對氮素營養(yǎng)吸收能力強，體內(nèi)氮含量相對較高，在干旱條件下，設(shè)施葡萄植株受水分制約導(dǎo)致氮吸收能力降低，從而植株氮濃度下降。

由表2可以看出，不同水分條件下設(shè)施葡萄植株臨界、最低和最高氮濃度稀釋曲線參數(shù)a不同，且隨著土壤含水量的減少各參數(shù)均呈減小趨勢，說明土壤水分適宜可提高葡萄植株臨界氮濃度，增強植株對氮的吸收能力?？赡苁怯捎谠谌彼畻l件下，葡萄植株長勢減弱，導(dǎo)致對氮的容納能力減弱，因此，植株體內(nèi)氮濃度下降。其次，參數(shù)b在不同水分條件下也不同，其隨著水分脅迫加重呈先減小后增加再減少的趨勢，說明輕度水分脅迫下葡萄植株氮濃度隨地上部生物量增加平緩降低，葡萄具有一定的耐旱性能，而到中度水分脅迫時，由于水分虧缺使得氮濃度隨植株增長迅速下降，重度水分脅迫條件下b值低是由于其植株初始氮濃度較低，降低幅度也因此較小。

2.2.2 模型驗證

采用另外3組重復(fù)的數(shù)據(jù)（n=80）對臨界氮濃度稀釋曲線經(jīng)驗?zāi)Ｐ瓦M(jìn)行驗證，如圖3所示，W1（正常水分條件）、W2（輕度水分脅迫）、W3（中度水分脅迫）和W4（重度水分脅迫）4種水分條件下，模型RMSE分別為0.273、0.263、0.277、0.283，R2分別為0.893、0.828、0.923、0.958，綜合RMSE和R2可以看出，模擬效果整體較好，說明該模型精度較高，可用于進(jìn)一步氮素營養(yǎng)診斷。

表2 不同水分條件下葡萄臨界氮濃度稀釋模型和氮素吸收模型參數(shù)值

注：Nmax和Nmin分別為模擬出的最高和最低氮濃度稀釋曲線模型（氮濃度稀釋邊界模型），Nuptmax和Nuptmin分別為最高和最低氮素吸收模型。

Note: Nmaxand Nminare the maximum and minimum nitrogen concentration dilution curve models by simulated (nitrogen concentration dilution boundary model). Nuptmaxand Nuptminare the highest and lowest nitrogen absorption models.

圖3 不同水分條件下葡萄臨界氮濃度稀釋模型模擬值與觀測值的關(guān)系

注：短線表示均方誤。下同。

Note : The short bar is mean square error.The same as below.

2.3 不同水分條件下臨界氮稀釋曲線模型的應(yīng)用

2.3.1 氮營養(yǎng)狀況診斷

圖4表示不同水分條件下設(shè)施葡萄植株在不同施氮水平處理下的營養(yǎng)指數(shù)（NNI）動態(tài)變化。由圖可知，植株的NNI隨施氮量增加而上升，且隨生長時間延長大體呈減小趨勢，表明處理后期植株對氮的需求量較大，W1（正常水分條件）、W2（輕度水分脅迫）、W3（中度水分脅迫）和W4（重度水分脅迫）水平下NNI變化區(qū)間分別為0.39～1.17、0.40～1.23、0.22～1.27和0.41～1.39，說明NNI隨土壤水分減少大致呈增加趨勢，在W1條件下N1（1.5倍施氮量）處理NNI>1，表明葡萄植株氮營養(yǎng)過剩；N2（正常施氮量）處理下前80d NNI>1，之后NNI<1，表明N2處理在前期氮營養(yǎng)充足，后期需氮量增加，導(dǎo)致營養(yǎng)不足；N3（0.5倍施氮量）和N4（無施氮）處理均為NNI<1，說明氮肥不足抑制了葡萄植株干物質(zhì)量的累積，增施氮肥后NNI明顯升高。在W2、W3條件下，NNI值呈現(xiàn)的變化規(guī)律大體與W1一致，且略大于W1條件下，在W4條件下，N1和N2處理時NNI>1，而N3和N4處理時 NNI<1，表明N1和N2處理下氮營養(yǎng)過剩，而N3和N4處理氮營養(yǎng)不足。因此經(jīng)綜合考慮得設(shè)施葡萄植株的較適施氮量應(yīng)在N1?N3之間，且更接近于N2。

2.3.2 氮吸收的營養(yǎng)診斷

由圖5可見，設(shè)施葡萄植株氮吸收量隨著施氮量和地上部生物量的增加呈增長趨勢，W1（正常水分條件）下植株氮吸收量整體較高，其后依次是W2（輕度水分脅迫）、W3（中度水分脅迫）和W4（重度水分脅迫），說明水分不足導(dǎo)致葡萄植株對氮的吸收量降低。由圖5a發(fā)現(xiàn)，臨界氮吸收量曲線（Nupt）處于N1（1.5倍施氮量）和N2（正常施氮量）曲線之間，這與NNI的診斷一致，說明灌溉量適宜時有助于葡萄植株對氮營養(yǎng)的吸收利用；由圖5b可得，相同灌溉條件下，N2處理植株氮吸收量低于N1處理，而地上部生物量高于N1處理，說明施氮量適宜時有助于葡萄植株地上部干物質(zhì)的累積，施氮量過多則適得其反，N3（0.5倍施氮量）和N4（無施氮）低于臨界氮吸收量，說明施氮不足導(dǎo)致作物基本的氮營養(yǎng)需求無法得到滿足，因此地上部干物質(zhì)量累積較慢，對葡萄植株生長產(chǎn)生抑制作用。

圖4 不同水氮條件下葡萄植株營養(yǎng)指數(shù)（NNI）的動態(tài)變化

圖5 不同水氮條件下葡萄地上部生物量與氮吸收量的關(guān)系

3 結(jié)論與討論

3.1 討論

為避免實際生產(chǎn)中過度施氮的現(xiàn)象，明確設(shè)施葡萄植株在不同生長時期所需臨界氮濃度，以合理指導(dǎo)施氮量，本研究構(gòu)建了葡萄植株在不同灌溉條件下的臨界、最高和最低氮濃度稀釋曲線模型，并對不同水氮條件下葡萄植株的氮營養(yǎng)狀況進(jìn)行了診斷，發(fā)現(xiàn)施氮量對該模型基本形式無顯著影響，但對其統(tǒng)計學(xué)參數(shù)值影響較大。所構(gòu)建的模型參數(shù)不同于前人所建立的作物臨界氮濃度稀釋曲線模型[5]，主要原因[4,20]一是前人構(gòu)建的模型是基于無氮素限制條件下得到的；二是由于作物品種、類型及內(nèi)在生理生態(tài)機制的不同；三是不同的栽培方式也會使參數(shù)產(chǎn)生差異；四是由于作物生長的氣候條件不同。任一模型應(yīng)用于特定作物時，模型參數(shù)都需重新設(shè)定。本研究表明，參數(shù)a隨土壤水分含量增加而增加，說明在適當(dāng)?shù)乃謼l件下葡萄植株對氮營養(yǎng)的吸收能力增強，有利于干物質(zhì)的積累。與楊慧等[20?21]構(gòu)建的番茄植株臨界氮濃度模型的結(jié)論一致。本研究中參數(shù)a值在正常水分條件（W1）下較大，葡萄植株地上部干物質(zhì)累積量較高；該值在重度水分脅迫（W4）下較小，是因為水分過度虧缺限制了葡萄植株對氮肥的吸收利用，從而抑制了葡萄植株的生長。本研究所構(gòu)建的在不同水氮條件下的臨界氮濃度稀釋曲線模型相比于僅由氮素單因素所構(gòu)建的模型更具普適性，但此模型僅在同一地點，由同一品種的試驗條件下構(gòu)建，因此還待進(jìn)一步完善。

NNI可實時診斷葡萄植株的氮肥營養(yǎng)狀況，量化植株受氮肥制約的程度[27]，是進(jìn)行合理施氮的重要依據(jù)。前人研究表明，NNI能較好地反映油菜、玉米、高粱、番茄、包心菜、馬鈴薯、向日葵、甜椒等植株的氮營養(yǎng)狀況，基于NNI來診斷植株生長進(jìn)程中的氮素營養(yǎng)盈虧，從而及時采取措施，調(diào)整水肥管理方法，提高作物地上部干物質(zhì)的累積速率，為精確施氮提供理論依據(jù)與指導(dǎo)。本研究對各水氮條件下葡萄植株不同時期NNI和葡萄植株地上部氮吸收進(jìn)行分析，發(fā)現(xiàn)葡萄植株NNI隨土壤水分的提高而降低，隨施氮量增加而上升，隨生長時間延長而減小的趨勢，葡萄植株氮素吸收量隨著施氮量和地上部干物質(zhì)量的增加呈增加的變化趨勢，本研究分析初步認(rèn)為，葡萄的較適施氮量介于17～25.5g×plant?1，作物干物質(zhì)量的有效積累需要吸收適宜的氮肥，則該施氮量有助于作物對氮營養(yǎng)的利用。同時，葡萄植株的地上部生物量和氮素吸收量隨著土壤水分的增加而增加，其中W1（70%～80%）適宜施氮量為N1（25.5g×plant?1）、W2（60%～70%）適宜施氮量為N2（17g×plant?1），兩者均可促進(jìn)設(shè)施葡萄植株對氮素營養(yǎng)的吸收利用，使地上部生物量增加最終分別達(dá)11.29t·hm?2和10.92t·hm?2。基于成本考慮，在W2的灌水條件下施用N2更適合應(yīng)用于實踐。本研究還認(rèn)為在一定土壤水分條件下，高氮和低氮均會使葡萄植株的生物量減少，趙薇等[28]以番茄為試驗對象也得出了一致結(jié)論。葡萄植株地上部生物量隨施氮量增加呈先增后降的變化趨勢，說明葡萄植株對氮的容納有一定限度，高于臨界氮吸收量后增施氮肥不僅不能同時提高作物干物質(zhì)量，還會抑制植株對氮的吸收利用。

3.2 結(jié)論

（1）以1a生葡萄品種“紅提”為供試材料，設(shè)置4個灌溉水平：W1（70%～80%）、W2（60%～70%）、W3（50%～60%）、W4（30%～40%），4個氮素水平（25.5、17、8.5和0g×plant?1），構(gòu)建了不同灌溉條件下設(shè)施葡萄植株的臨界、最大和最小氮濃度稀釋曲線模型，結(jié)果表明，不同時期的植株最大生物量與氮濃度之間存在冪函數(shù)關(guān)系，由模型參數(shù)得，植株的氮素吸收量和地上部生物量隨灌溉水平提高而增加。

（2）氮營養(yǎng)指數(shù)模型與氮素吸收模型診斷各水分條件下設(shè)施葡萄植株氮營養(yǎng)結(jié)果一致：W1、W2及W3灌溉水平下，每株葡萄植株較適施氮量在17～25.5g，W4條件下在8.5～17g。

（3）分析氮素營養(yǎng)對植株生物量的影響效應(yīng)認(rèn)為，水分條件為W1和W2時較利于葡萄對氮營養(yǎng)元素的吸收利用，地上部生物量較高。在適宜施氮量情況下，灌溉水平為W2時可獲得僅低于W1灌溉水平下約3.3%的生物量，大大減少了水和肥料的成本?？紤]到溫室無降水的供應(yīng)，灌水量W2（60%～70%）施氮量為17g×plant?1左右為較適水氮配比。

[1] 穆維松,高陽,王秀娟,等.我國設(shè)施葡萄生產(chǎn)的成本收益比較研究[J].中外葡萄與葡萄酒,2014,(3):20-24.

Mu W S,Gao Y,Wang X J,et al.Comparative cost-benefit analysis on grape protected cultivation[J].Sino-Overseas Grapevine & Wine,2014,(3):20-24.(in Chinese)

[2] 吉沐祥,李國平,芮東明,等.江蘇省鮮食葡萄病蟲害綠色防控技術(shù)規(guī)程[J].江蘇農(nóng)業(yè)科學(xué),2015,43(7):107-109.

Ji M X,Li G P,Rui D M,et al.Technical specification for green control of diseases and insect pests of fresh grapevine in Jiangsu province[J].Jiangsu Agricultural Sciences,2015,43(7):107-109.(in Chinese)

[3] Greenwood D J,Gastal F,Lemaire G,et al.Growth rate and %N of field grown crops: theory and experiments[J].Annals of Botany,1991,67(2):181-190.

[4] 王新,馬富裕,刁明,等.滴灌番茄臨界氮濃度、氮素吸收和氮營養(yǎng)指數(shù)模擬[J].農(nóng)業(yè)工程學(xué)報,2013,29(18):107-116.

Wang X,Ma F Y,Diao M,et al.Simulation of critical nitrogen concentration,nitrogen uptake and nitrogen nutrition index of processing tomato with drip irrigation[J].Transactions of the CSAE,2013,29(18):107-116.(in Chinese)

[5] Lemaire G,Gastal F.N uptake and distribution in plant canopies[M].Diagnosis of the Nitrogen Status in Crops, 1997,53(370):3-43.

[6] Justes E,Mary B,Meynard J M,et al.Determination of a critical nitrogen dilution curve for winter wheat crops[J]. Annals of Botany,1994,74(4):397-407.

[7] Sheehy J E,Dionora M J A,Mitchell P L,et al.Critical nitrogen concentrations: implications for high-yielding rice(L.) cultivars in the tropics[J].Field Crops Research, 1998,59(1):31-41.

[8] Colnenne C,Meynard J M,Reau R,et al.Determination of a critical nitrogen dilution curve for winter oilseed rape[J]. Annals of Botany,1998,81(2):311-317.

[9] Daniel P,Gilles L.Relationships between dynamics of nitrogen uptake and dry matter accumulation in maize crops: determination of critical N concentration[J].Plant and Soil, 1999,216(1-2):65-82.

[10] Debaeke P,Oosterom E J V,Justes E,et al.A species-specific critical nitrogen dilution curve for sunflower(L.)[J].Field Crops Research, 2012, 136: 76-84.

[11] Tei F,Benincasa P,Guiducci M.Critical nitrogen concentration in processing tomato[J].European Journal of Agronomy,2002, 18(1):45-55.

[12] Ekbladh G,Witter E.Determination of the critical nitrogen concentration of white cabbage[J].European Journal of Agronomy,2010,33(4):276-284.

[13] Giletto C M,Echeverría H E.Critical nitrogen dilution curve for processing potato in Argentinean humid pampas[J]. American Journal of Potato Research,2012,89(2):102-110.

[14] Oosterom E J V,Carberry P S,Muchow R C.Critical and minimum N contents for development and growth of grain sorghum[J].Field Crops Research,2001,70(1):55-73.

[15] Yasuor H,Ben-Gal A,Yermiyahu U,et al.Nitrogen management of greenhouse pepper production:agronomic nutritional, and environmental implications[J].Hort Science,2013,48(10): 1241-1249.

[16] 薛曉萍,沙奕卓,郭文琦,等.棉花蕾花鈴生物量、氮累積特征及臨界氮濃度稀釋模型[J].生態(tài)學(xué)報,2008,28(12): 6204-6211.

Xue X P,Sha Y Z,Guo W Q,et al.Accumulation characteristics of biomass and nitrogen and critical nitrogen concentration dilution model of cotton reproductive organ[J].Acta Ecologica Sinica,2008,28 (12):6204-6211.(in Chinese)

[17]梁效貴,張經(jīng)廷,周麗麗,等.華北地區(qū)夏玉米臨界氮稀釋曲線和氮營養(yǎng)指數(shù)研究[J].作物學(xué)報,2013,39(2):292-299.

Liang X G,Zhang J T,Zhou L L,et al.Critical nitrogen dilution curve and nitrogen nutrition index for summer maize in North China Plain[J].Acta Agronomica Sinica, 2013,39(2):292-299.(in Chinese)

[18] 強生才,張富倉,向友珍,等.關(guān)中平原不同降雨年型夏玉米臨界氮稀釋曲線模擬及驗證[J].農(nóng)業(yè)工程學(xué)報,2015,31 (17):168-175.

Qiang S C,Zhang F C,Xiang Y Z,et al.Simulation and verification of critical nitrogen dilution curve for summer maize in Guanzhong Plain during different rainfall years[J].Transactions of the CSAE,2015,31(17): 168-175.(in Chinese)

[19]向友珍,張富倉,范軍亮,等.基于臨界氮濃度模型的日光溫室甜椒氮營養(yǎng)診斷[J].農(nóng)業(yè)工程學(xué)報,2016,32(17):89-97.

Xiang Y Z,Zhang F C,Fan J L,et al.Nutrition diagnosis for N in bell pepper based on critical nitrogen model in solar greenhouse[J].Transactions of the CSAE,2016,32(17):89-97. (in Chinese )

[20]楊慧,曹紅霞,柳美玉,等.水氮耦合條件下番茄臨界氮濃度模型的建立及氮素營養(yǎng)診斷[J]. 植物營養(yǎng)與肥料學(xué)報,2015,21(5):1234-1242.

Yang H,Cao H X,Liu M Y,et al.Simulation of critical nitrogen concentration and nitrogen nutrition index of tomato under different water and nitrogen conditions[J]. Journal of Plant Nutrition and Fertilizer,2015,21(5):1234- 1242. (in Chinese )

[21] 劉明,呂愛鋒,武建軍,等.干旱對農(nóng)業(yè)生態(tài)系統(tǒng)影響研究進(jìn)展[J].中國農(nóng)學(xué)通報,2014,30(32):165-171.

Liu M,Lv A F,Wu J J,et al.A review of impacts of drought on agro-ecosystem[J].Chinese Agricultural Science Bulletin, 2014,30(32):165-171.(in Chinese)

[22] 王奎穩(wěn).設(shè)施葡萄秋季施肥技術(shù)[J].落葉果樹,2016,48(3): 30.

Wang K W.Autumn fertilization techniques for protected grapes[J].Deciduous Fruits,2016,48(3):30.(in Chinese)

[23]倪榮梅.不同水肥管理對設(shè)施葡萄產(chǎn)量和品質(zhì)的影響研究[J].新疆農(nóng)墾科技,2016,39(4):62-64.

Ni R M.Effects of different water and fertilizer management on yield and quality of protected grape[J].Xinjiang Farm Research of Science and Technology,2016,39(4):62-64.(in Chinese)

[24]Bao S D.Soil and agricultural chemistry analysis[M].Beijing:China Agriculture Press,1999.

[25]Gastal F,Lemaire G.N uptake and distribution in crops: an agronomical and ecophysiological perspective[J].Journal of Experimental Botany,2002,53(370):789-799.

[26]岳松華,劉春雨,黃玉芳,等.豫中地區(qū)冬小麥臨界氮稀釋曲線與氮營養(yǎng)指數(shù)模型的建立[J].作物學(xué)報,2016,42(6):909-916.

Yue S H,Liu C Y,Huang Y F,et al.Simulating critical nitrogen dilution curve and modeling nitrogen index in winter wheat in central Henan area[J].Acta Agronomica Sinica,2016,42(6):909-916.(in Chinese)

[27]Meynard J M,David G.Diagnostic de lelaboration du rendement des cultures[J].Cahiers Agriculture,1992,1(1):9-19.

[28]趙薇,張治平,繆旻珉.番茄臨界氮濃度模型的建立及氮素營養(yǎng)診斷[J].中國農(nóng)學(xué)通報,2018,34(22): 37-44.

Zhao W,Zhang Z P,Miao M M.Establishment of tomato critical nitrogen concentration model and diagnosis of nitrogen nutrition[J].Chinese Agricultural Science Bulletin,2018,34(22):37-44.(in Chinese)

Establishment of a Critical Nitrogen Dilution Model for Grapes and Nitrogen Nutrition Diagnosis under Different Water Conditions

LI Jia-shuai1, YANG Zai-qiang1,2, LI Yong-xiu1,2, JIANG Meng-yuan1, ZHAO He-li1, WEI Ting-ting1, ZHANG Xu-ran1

(1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China; 2.Jiangsu Provincial Key Laboratory of Agrometeorology, Nanjing 210044)

In this paper, a two-factor pot experiment was carried out in the Venlo-type greenhouse, using the annual grape of "Hongti" variety as experimental material. There were four irrigation levels and four nitrogen application levels in the two-factor pot experiment. Four irrigation levels were normal irrigation (W1, 70%?80% of the maximum water holding capacity in the field), slight drought stress (W2, 60%?70%), moderate drought stress (W3, 50%?60%) and extra severe drought tress (W4, 30%?40%). And four nitrogen application levels were 1.5 times normal recommended nitrogen application rate(N1, 25.5g×plant?1), normal recommended nitrogen application rate(N2, 17g×plant?1), 0.5 times normal recommended nitrogen application rate(N3, 8.5g×plant?1) and no nitrogen application(N4, 0g×plant?1). The nitrogen concentration in the plant and the aboveground biomass of the plant were observed every 10 days. The minimum nitrogen concentration corresponding to the maximum biomass obtained by grape plants under different water conditions in a certain period is called Nitrogen critical concentration (Nc), it can be used to construct the model of nitrogen critical concentration dilution curve of grape, which can further establish Nitrogen uptake model (Nupt) and Nitrogen Nutrition Index model (NNI) for quantitative diagnosis of nitrogen nutrition status of grapes under different water conditions.The results showed: there was a power function relationship between the nitrogen critical concentration and aboveground biomass of grape plants, which embodied as the nitrogen critical concentration, nitrogen uptake and aboveground biomass of grape plants all increased as the irrigation amount increasing. As nitrogen application increased, the biomass of grape plants under treatments of W1 and W2 increased, while that under W3 and W4, increased first and then decreased. Under the same water condition, the nitrogen concentration increased with the increase of nitrogen application and decreased with the progress of grapes growth. The established Nuptand NNI models can be used to effectively diagnose the deficiency of nitrogen nutrients in plants.

Grape; Nitrogen critical concentration; Nitrogen nutrition index; Nitrogen uptake; Dry matter accumulation

10.3969/j.issn.1000-6362.2019.08.005

李佳帥,楊再強,李永秀.不同水分條件下葡萄臨界氮稀釋曲線模型的建立及氮素營養(yǎng)診斷[J].中國農(nóng)業(yè)氣象,2019,40(8):523-533

2019?02?26

。E-mail：yzq@nuist.edu.cn

江蘇省科技支撐項目（BE2015693）

李佳帥（1993?），女，研究方向為設(shè)施農(nóng)業(yè)氣象。E-mail：596993884@qq.com