石小虎，蔡煥杰

?

基于葉片SPAD估算不同水氮處理下溫室番茄氮營(yíng)養(yǎng)指數(shù)

石小虎1,2，蔡煥杰2※

（1.太原理工大學(xué)水利科學(xué)與工程學(xué)院，太原 030024； 2. 西北農(nóng)林科技大學(xué)中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院，楊凌 712100）

為了探討臨界氮稀釋曲線模型在西北地區(qū)溫室番茄不同水分處理下的適用性以及采用SPAD儀快速準(zhǔn)確診斷氮營(yíng)養(yǎng)狀況，該研究以“麗娜”番茄為材料，2013－2015年在陜西省楊凌區(qū)溫室內(nèi)進(jìn)行水分和氮素處理試驗(yàn)，水分處理設(shè)置4個(gè)水平，分別為全生育期充分灌水處理、僅苗期虧水50%、苗期開(kāi)花期連續(xù)虧水50%和全生育期虧水50%；氮素處理設(shè)置3個(gè)水平，施氮量分別為0、150和300 kg/hm2，通過(guò)2013－2015年試驗(yàn)數(shù)據(jù)對(duì)臨界氮濃度稀釋曲線模型進(jìn)行率定和驗(yàn)證，并將該模型參數(shù)與番茄全生育期平均日耗水量建立相關(guān)關(guān)系，提高了臨界氮濃度稀釋模型在不同水分條件下的適用性。結(jié)果表明通過(guò)番茄全生育期平均日耗水量和臨界氮濃度稀釋曲線模型估算得到的臨界氮濃度估算值和實(shí)際計(jì)算值有較好的一致性，其絕對(duì)誤差為0.13～0.34 g/(100 g)，標(biāo)準(zhǔn)誤差為0.14～0.39 g/(100 g)，決定系數(shù)為0.94～0.99，因此采用該方法可以對(duì)西北地區(qū)溫室番茄不同水分處理下臨界氮濃度稀釋進(jìn)行準(zhǔn)確估算。通過(guò)2013－2015年試驗(yàn)數(shù)據(jù)分析番茄不同葉位葉片SPAD值和氮營(yíng)養(yǎng)指數(shù)（nitrogen nutrition index，NNI）之間相關(guān)性，結(jié)果表明番茄中位葉片SPAD值與氮營(yíng)養(yǎng)指數(shù)（NNI）有良好的線性相關(guān)性（決定系數(shù)為0.77～0.98），且該相關(guān)系數(shù)值與番茄日耗水量呈極顯著相關(guān)關(guān)系，因此通過(guò)番茄日耗水量可以估算出NNI與中位葉片SPAD值之間的線性關(guān)系，估算出NNI=1時(shí)的中位葉片SPAD值，并以此SPAD值進(jìn)行氮營(yíng)養(yǎng)診斷。該研究可為西北地區(qū)溫室番茄實(shí)時(shí)氮營(yíng)養(yǎng)診斷和優(yōu)化氮素管理提供了較好的理論參考。

葉綠素；水分；溫室；番茄；水分處理；臨界氮濃度；氮營(yíng)養(yǎng)指數(shù)

0 引 言

溫室蔬菜栽培在中國(guó)西北地區(qū)種植規(guī)模發(fā)展迅速，番茄作為溫室蔬菜的主要品種，生育期短，在反季節(jié)蔬菜中占有很大比例。傳統(tǒng)的肥大水勤仍然是當(dāng)?shù)刂饕乃使芾矸绞剑蔬^(guò)量施用的現(xiàn)象比較普遍，過(guò)量施氮不僅能使蔬菜產(chǎn)量降低，還能導(dǎo)致土壤硝態(tài)氮積累，引發(fā)土壤次生鹽漬化[1-2]。臨界氮濃度是一定的生長(zhǎng)時(shí)期內(nèi)獲得最大生物量時(shí)的最小氮濃度值[3]。因此明確番茄干物質(zhì)形成過(guò)程中不同生育階段的臨界氮濃度是科學(xué)診斷植株氮營(yíng)養(yǎng)狀況，是實(shí)現(xiàn)番茄各生育階段氮肥合理施用的基礎(chǔ)。為了準(zhǔn)確地計(jì)算出作物不同生育期階段的臨界氮濃度，Greenwood等[4]在1990年提出了關(guān)于C3、C4作物臨界氮濃度與地上部生物量關(guān)系的通用模型，后經(jīng)Lemaire和Gastal等[5]大量試驗(yàn)，修正了其中的參數(shù)[6]和參數(shù)[7]，提出了關(guān)于C3、C4作物的新模型，由于模型是基于多個(gè)試驗(yàn)平均得到的結(jié)果，且供試作物不能代表所有作物。近些年國(guó)內(nèi)外學(xué)者對(duì)不同作物進(jìn)行臨界氮濃度稀釋模型的研究，主要集中在棉花[8]、小麥[9-10]、番茄[11-13]、高粱[14]、玉米[15-16]等作物，研究均表明臨界氮濃度稀釋曲線模型可較好的描述地上部生物量與氮濃度的關(guān)系，但由于試驗(yàn)地、供試作物及試驗(yàn)處理等因素不同，導(dǎo)致了氮濃度稀釋模型參數(shù)和也有較大的差異，因此需要根據(jù)實(shí)際情況對(duì)模型參數(shù)進(jìn)行校正。根據(jù)臨界氮濃度稀釋曲線，Lemaire[17]等定義了氮素營(yíng)養(yǎng)指數(shù)（nitrogen nutrition index，NNI），即地上部實(shí)測(cè)氮濃度和臨界氮濃度的比值。當(dāng)NNI=1時(shí)，表明作物體內(nèi)氮素營(yíng)養(yǎng)合適，NNI>1表明氮營(yíng)養(yǎng)過(guò)剩，NNI<1表明氮營(yíng)養(yǎng)不足。強(qiáng)生才等[16]研究表明在不同降雨年型下，所構(gòu)建的夏玉米臨界氮稀釋曲線模型有所不同，降雨量不同導(dǎo)致了模型參數(shù)的不同。楊慧等[12]在不同水分條件下構(gòu)建了溫室盆栽番茄臨界氮稀釋曲線模型，并表現(xiàn)出不同灌水量處理氮濃度稀釋模型參數(shù)不同。因此不同水分處理下臨界氮稀釋曲線模型參數(shù)有較大的差別，需要根據(jù)實(shí)際的種植和處理情況進(jìn)行氮素營(yíng)養(yǎng)診斷。

氮素營(yíng)養(yǎng)指數(shù)（NNI）診斷法作為傳統(tǒng)氮素診斷技術(shù)，這種方法判定準(zhǔn)確、可信度高，但是需要的樣本量較大、并且需要破壞植物樣本，且操作復(fù)雜、費(fèi)時(shí)、不能實(shí)時(shí)實(shí)地得出結(jié)果，限制了這種方法在實(shí)際操作中的推廣應(yīng)用[18]，而利用現(xiàn)代儀器氮素診斷技術(shù)則可以快速、便捷和非破壞性的進(jìn)行作物氮素診斷，如圖像及計(jì)算機(jī)視覺(jué)診斷、高光譜遙感診斷和SPAD儀診斷等。視覺(jué)技術(shù)是通過(guò)對(duì)人的視覺(jué)進(jìn)行計(jì)算機(jī)計(jì)算模擬，通過(guò)圖片采集工具對(duì)作物冠層葉片的圖像綠光值來(lái)判斷植株的氮營(yíng)養(yǎng)情況[19]，但是存在攜帶不方便，不能進(jìn)行實(shí)時(shí)實(shí)地的跟蹤觀察測(cè)定等缺陷；高光遙感技術(shù)是某一光譜區(qū)域?qū)庾V進(jìn)行詳細(xì)劃分，進(jìn)一步得到更多波段的光譜信息，與多波段遙感相比，高光譜遙感光譜分辨率較高，Maderia研究表明水稻葉片中的葉綠素含量高低與光譜特征之間兩者存在相關(guān)性，可以通過(guò)光譜檢測(cè)植株的氮營(yíng)養(yǎng)狀況[20]，但該方法存在數(shù)據(jù)獲得需經(jīng)過(guò)反復(fù)的轉(zhuǎn)換，影響因素較多，儀器昂貴等缺陷；葉綠素計(jì)（SPAD計(jì)）是一種便攜式光譜儀，具有攜帶方便、可實(shí)時(shí)實(shí)地進(jìn)行觀察測(cè)定等優(yōu)點(diǎn)，它是以葉綠素對(duì)對(duì)紅光和近紅外光的差異吸收，根據(jù)透過(guò)葉片的透射光特征來(lái)表達(dá)葉片的葉綠素相對(duì)含量，進(jìn)而可以根據(jù)葉片葉綠素含量與葉片氮含量的關(guān)系來(lái)推算植物的氮營(yíng)養(yǎng)狀況[21]，Errecart等[22]也通過(guò)建立葉片SPAD值與植株NNI值之間的相關(guān)性，利用葉片SPAD值來(lái)快速模擬作物NNI值，Yang等[23]比較水稻不同葉位的SPAD值和NNI之間的回歸關(guān)系，表明該擬合關(guān)系在年際間表現(xiàn)穩(wěn)定，推測(cè)利用SPAD值可判斷NNI值，用于水稻的氮素營(yíng)養(yǎng)診斷。本研究為了提高臨界氮濃度稀釋曲線模型在西北地區(qū)溫室番茄不同水分處理的適用性以及確定適當(dāng)葉位葉片SPAD值估算出番茄NNI值進(jìn)行氮素營(yíng)養(yǎng)診斷。因此本研究通過(guò)2013－2015年2 a不同水分處理處理下3個(gè)氮素處理的溫室番茄試驗(yàn)，率定和驗(yàn)證了臨界氮稀釋曲線模型在西北地區(qū)溫室番茄不同水分狀況的適用性，篩選出適宜的測(cè)定葉位以及根據(jù)葉片SPAD值確定氮營(yíng)養(yǎng)指數(shù)，旨在為不同水分條件下番茄氮素合理利用、氮營(yíng)養(yǎng)狀況的診斷及氮素優(yōu)化管理提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)區(qū)概況

試驗(yàn)于2013－2015年在陜西省楊凌區(qū)溫室內(nèi)進(jìn)行（34°18′N、108°4′E，海拔521 m）。試驗(yàn)溫室為非加熱型自然通風(fēng)溫室，主體為鋼架結(jié)構(gòu)，用塑料薄膜覆蓋，南北走向（長(zhǎng)×寬×高，50 m×6.5 m×4.6 m），溫室頂部和底部各設(shè) 1 m 寬通風(fēng)口，并配置手動(dòng)啟閉裝置，當(dāng)溫室內(nèi)溫度＞35 ℃或＜10 ℃時(shí)，通過(guò)開(kāi)啟或關(guān)閉通風(fēng)口來(lái)調(diào)節(jié)溫室內(nèi)溫度。溫室內(nèi)部番茄為南北方向種植。該地溫室0～60 cm土壤內(nèi)，黏粒（＜2m）為22%，粉粒（2～20m）為56%，砂粒（≥20～2 000m）為22%，有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)為3.44%，容重為1.41 g/cm3，飽和含水率θ為0.41 cm3/cm3，田間持水量θ為0.34 cm3/cm3，凋萎含水量θ為 0.14 cm3/cm3。

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

本研究設(shè)計(jì)2個(gè)因素：水分和氮素。4個(gè)灌水水平：全生育期充分灌水處理（W1）、僅苗期虧水50%（W2）、苗期開(kāi)花期連續(xù)虧水50%（W3）、全生育期虧水50%（W4）；3個(gè)施氮水平，分別為無(wú)氮（N0：0）、中氮（N150：150 kg/hm2）和高氮（N300：300 kg/hm2），水分和氮素進(jìn)行完全組合，如表1所示。各試驗(yàn)處理進(jìn)行3次重復(fù)共36個(gè)小區(qū)，各小區(qū)面積為6.5 m×2.4 m=15.6 m2，完全隨機(jī)布置，小區(qū)之間用埋深60 cm的塑料薄膜隔離，試驗(yàn)小區(qū)布置如圖1所示。

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

注：W1～W4分別為全生育期充分灌水處理、僅苗期虧水50%處理、苗期開(kāi)花期連續(xù)虧水50%處理和全部虧水50%處理；N0、N150和N300為施氮量為0、150和300 kg·hm-2處理；為灌水量，mm。

Note: W1-W4are full irrigation treatment, deficit of 50% only in seedling stage, deficit of 50% both in seedling stage and flowering stage, deficit of 50% in the whole growth stage; N0, N150and N300are no nitrogen application, nitrogen application rate of 150 kg·hm-2and nitrogen application rate of 300 kg·hm-2, respectively;is irrigation amount, mm.

圖1 溫室番茄試驗(yàn)小區(qū)處理布置

本研究采用膜下滴灌灌溉方式，灌水從定植后15d開(kāi)始，各水分處理灌水周期為7 d左右，設(shè)置充分灌溉處理（W1）灌水上限為田間持水率的90%[24]。石小虎[25]研究表明番茄根系在0～60 cm土層內(nèi)分布，因此本研究設(shè)置計(jì)劃濕潤(rùn)層深度為60 cm。根據(jù)王峰等[24]研究，計(jì)算出充分灌水處理（W1）各生育期灌水量，其他水分處理僅在灌水量上進(jìn)行減少，灌水時(shí)間和次數(shù)均與W1處理相同；根據(jù)石小虎[25]研究將高氮水平（N300）定為300 Kg/hm2（以N計(jì)），N0和N150處理只在施氮量上進(jìn)行減少。氮肥選用尿素（含氮質(zhì)量分?jǐn)?shù) 46%），定植之前基施 40%，剩余60%分別在定植后（days after transferring，DAT）70、90和110 d溶化到水中隨滴灌平均施入。

試驗(yàn)用品種為番茄“麗娜”，番茄生育期劃分為緩苗期（2013-08-10－08-31、2014-08-15－09-08），苗期（2013-09-01－10-05、2014-09-09－10-12），開(kāi)花期（2013-10-06－12-05、2014-10-13－12-12）和成熟期（2013-12-06－2014-01-10、2014-12-13－2015-01-21）。種植方式為當(dāng)?shù)氐湫偷钠饓鸥材ぴ耘嗄Ｊ?，壟?0 cm、壟寬80 cm，番茄幼苗按單穴單株定植在壟的兩側(cè)，其寬行距為80 cm，跨溝窄行距為40 cm，株距為40 cm，種植密度為4.2株/m2。定植前在溫室內(nèi)均勻施入等量的磷肥200 kg/hm2（以P計(jì)）和鉀肥300 kg/hm2（以K計(jì)）。定植時(shí)灌定植水20 mm，定植后14 d內(nèi)不灌水，以利于蹲苗，待番茄幼苗成活后再進(jìn)行試驗(yàn)處理。定植當(dāng)天沿溫室南北走向鋪設(shè)寬1.2 m，厚0.005 mm地膜，番茄開(kāi)花后用細(xì)繩將番茄懸吊在溫室上方的鐵絲上，并每3 d人工授粉1次，同時(shí)進(jìn)行噴藥等農(nóng)作管理。全生育期內(nèi)，每株番茄留3穗果后摘心，每穗留4～5個(gè)番茄。番茄成熟后每2 d進(jìn)行采摘1次，其他農(nóng)作管理按當(dāng)?shù)爻Ｒ?guī)進(jìn)行。

1.3 觀測(cè)項(xiàng)目與方法

1）土壤含水率

將Trime系列土壤水分測(cè)量?jī)x（IMKO Corp., Germany）埋設(shè)在距離植株 20 cm 位置處，分別測(cè)定各小區(qū)寬行、窄行和株間的土壤含水率，灌水前后各測(cè)1次，從表層到 60 cm 深每隔15 cm測(cè)量1次，計(jì)算時(shí)取其平均值。

2）充分灌水處理灌水量

灌水從定植后15 d開(kāi)始，充分灌水處理灌水上限為田間持水量（FC）的90%，其灌水量（mm）為

=10(0.9FC?i)Z（1）

式中θ為灌水前的土壤含水量，cm3/cm3；Z為計(jì)劃濕潤(rùn)層深度，cm，本文取60 cm。

3）番茄日耗水量

由于溫室內(nèi)地勢(shì)平坦，沒(méi)有地表徑流；當(dāng)?shù)氐叵滤裨O(shè)較深，地下水對(duì)番茄的補(bǔ)給可以忽略；溫室可以阻隔降雨的進(jìn)入，因此降雨量可以忽略；根據(jù)溫室內(nèi)自制的滲透儀可以得到60 cm土層處沒(méi)有深層滲漏。因此水量平衡公式可簡(jiǎn)化為

ET=(I+Δ)/（2）

式中ET為天內(nèi)平均日耗水量，mm/d；I為天內(nèi)的灌水量，mm；Δ為天內(nèi)土壤含水量的變化量，mm。

4）番茄地上部生物量

定苗后每隔20 d左右進(jìn)行破壞性取樣，每次均取3株。每次取樣稱量番茄地上部莖、葉、果鮮質(zhì)量，在105 ℃下烘15 min殺青，72 ℃下烘至恒質(zhì)量，計(jì)算其生物量。

5）植株各器官含氮量測(cè)定

將各處理的干物質(zhì)分器官粉碎后過(guò)篩，用H2SO4-H2O2消煮法和凱氏定氮儀（FOSS 2300 型）測(cè)定各器官全氮含量，并計(jì)算植株全氮含量。各器官氮累積量(kg/hm2)=器官含氮量(%)×器官生物量(kg/hm2)，所有器官氮累積量相加得到地上部植株氮累積量。植株含氮量(%)=植株氮累積量(kg/hm2)/植株生物量(kg/hm2)。

6）葉片SPAD值

各處理取整齊一致的3株植株，分別對(duì)不同節(jié)位葉片（由下向上依次為5、6、7、8、9、10、11、12節(jié)位）SPAD值采用SPAD-502儀測(cè)定，將第5、6和7節(jié)位所測(cè)數(shù)據(jù)平均后作為下位葉SPAD值，第8、9和10節(jié)位所測(cè)數(shù)據(jù)平均后作為中位葉SPAD值，第11和12節(jié)位所測(cè)數(shù)據(jù)平均后作為上位葉SPAD值。

1.4 模型描述

1.4.1 臨界氮稀釋曲線模型

根據(jù) Justes等[26]1994年提出的臨界氮濃度的定義及計(jì)算方法，綜合薛曉萍[27]、梁效貴等[28]、王新等[11]和楊慧等[12]關(guān)于棉花、夏玉米和番茄臨界氮濃度稀釋曲線模型的建模思路。本文臨界氮濃度稀釋曲線模型的構(gòu)建方法如下：

1）對(duì)比分析不同氮水平下每次取樣地上部生物量及相對(duì)應(yīng)的氮濃度值，通過(guò)方差分析對(duì)作物生長(zhǎng)受氮素營(yíng)養(yǎng)限制與否的氮素水平進(jìn)行分類；2）對(duì)于施氮量不能滿足作物生長(zhǎng)需求的試驗(yàn)資料，其地上部生物量與氮濃度值間的關(guān)系以線性曲線擬合；3）對(duì)于作物生長(zhǎng)不受氮素影響的施氮水平，其地上部生物量的平均值用以代表生物量的最大值；4）每次取樣日的理論臨界氮濃度由上述線性曲線與以最大生物量為橫坐標(biāo)的垂線的交點(diǎn)的縱坐標(biāo)決定。

依據(jù)Lemaire和Salette等[29]提出的臨界氮濃度與地上部生物量關(guān)系的方程式，臨界氮濃度稀釋曲線模型為：

N=·DW-b（3）

式中N為臨界氮濃度值，g/(100 g)；參數(shù)為作物地上部生物量為103kg/hm2時(shí)植株的臨界氮濃度；DW為作物地上部生物量的最大值，103kg/hm2；為決定臨界氮濃度稀釋曲線斜率的統(tǒng)計(jì)學(xué)參數(shù)。

1.4.2 氮營(yíng)養(yǎng)指數(shù)

為了進(jìn)一步明確作物的氮素營(yíng)養(yǎng)狀況，Lemaire等[29]提出了氮營(yíng)養(yǎng)指數(shù)（nitrogen nutrition index，NNI）的概念，可用式（4）來(lái)表示。

NNI=N/N（4）

式中NNI為氮營(yíng)養(yǎng)指數(shù)；N為地上部生物量氮濃度的實(shí)測(cè)值，g/(100 g)。

1.5 數(shù)據(jù)分析

采用 Microsoft Excel 2007和DPS軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行整理和分析，使用 Origin Pro 8.5軟件作圖。為了評(píng)價(jià)模型的精度，根據(jù)文獻(xiàn)[30-32]計(jì)算模型模擬值和計(jì)算值之間的決定系數(shù)（2）、絕對(duì)誤差指標(biāo)（mean absolute error，MAE）和標(biāo)準(zhǔn)誤差（root mean square error，RMSE）。

2 結(jié)果與分析

2.1 番茄地上部生物量與氮濃度值分析

不同水分和氮肥處理下番茄地上部生物量和植株氮素含量隨定植后天數(shù)變化過(guò)程如表2所示。由表2可以得知，番茄地上部生物量隨時(shí)間的推進(jìn)呈逐漸增大的趨勢(shì)，而番茄植株氮素含量隨時(shí)間推進(jìn)呈逐漸減小的趨勢(shì)，在番茄拉秧時(shí)（DAT=150 d）番茄地上部干物質(zhì)量達(dá)到最大值，為7.38×103～11.12×103kg/hm2，而植株氮素含量達(dá)到最小值，為1.43～1.96 g/(100 g)。充分灌水（W1）時(shí)，番茄地上部生物量和植株氮素含量均隨施氮量的增加呈顯著增加的趨勢(shì)；與充分灌水處理相比，僅苗期虧水（W2）時(shí)不會(huì)顯著影響拉秧時(shí)番茄地上部生物量和植株氮素含量，隨著虧水天數(shù)的增加，番茄地上部生物量和植株氮素含量減少量逐漸增加，全生育期虧水（W4）時(shí)，地上部生物量和植株氮素含量減少到最低。表明充分灌水時(shí)，增加施氮量有顯著的增加地上部干物質(zhì)和植株氮素含量的效應(yīng)，而非充分灌水時(shí)，施氮量增加到一定量（150 kg/hm2）繼續(xù)增加施氮量時(shí)不會(huì)顯著增加番茄地上部生物量。

表2 2013－2015年不同水分和氮素處理對(duì)番茄地上部生物量的影響

注：W1～W4分別為全生育期充分灌水處理、僅苗期虧水50%處理、苗期開(kāi)花期連續(xù)虧水50%處理和全部虧水50%處理；N0、N150和N300為施氮量為0、150和300 kg·hm-2處理；表中數(shù)據(jù)為各重復(fù)平均值，同一行數(shù)據(jù)具有相同字母的表示與同一年處理間沒(méi)有達(dá)到顯著性檢驗(yàn)（<0.05），下同。

Note: W1-W4are full irrigation treatment, deficit of 50% only in seedling stage, deficit of 50% both in seedling stage and flowering stage, deficit of 50% in the whole growth stage; N0, N150and N300are no nitrogen application, nitrogen application rate of 150 kg·hm-2and nitrogen application rate of 300 kg·hm-2, respectively; Data of the table represent average value and those with the same letters are not significantly different (<0.05), same as below.

2.2 臨界氮濃度稀釋模型構(gòu)建及驗(yàn)證

由于2013－2014年試驗(yàn)中在定植30 d時(shí)的各處理地上部生物量均小于1×103kg /hm2，生物量較小，粉碎后的樣品無(wú)法測(cè)定氮含量，故舍去此數(shù)據(jù)。利用Justes等[26]所描述的方法，分別將2013－2014年不同灌水條件下不同氮肥處理番茄地上部生物量和對(duì)應(yīng)的氮濃度進(jìn)行計(jì)算，得到各取樣日的臨界氮濃度。根據(jù)地上部生物量及對(duì)應(yīng)的臨界氮濃度，建立番茄臨界氮稀釋曲線，如圖2a所示。由圖2a可以看出不同水分處理臨界氮稀釋曲線的決定系數(shù)（2）為別為0.92～0.98，其擬合度均達(dá)到極顯著水平（＜0.01），表明該模型在不同水分處理下均可以較好的反映番茄臨界氮濃度和地上部生物量之間的關(guān)系。

圖2 不同水分處理下番茄氮濃度稀釋曲線構(gòu)建和驗(yàn)證

不同水分處理時(shí)臨界氮稀釋曲線參數(shù)、均不同，因此采用2013－2014年不同水分處理番茄全生育期平均日耗水量與對(duì)應(yīng)處理臨界氮稀釋曲線參數(shù)、進(jìn)行擬合，得到臨界氮稀釋曲線系數(shù)、與番茄全生育期平均日耗水量的相關(guān)關(guān)系，如圖3所示，番茄全生育期平均日耗水量如表3所示。

注：a和b為不同水分處理經(jīng)驗(yàn)系數(shù)；ET為不同水分處理番茄全生育期日耗水量。

表3 2013－2015年不同水氮處理番茄耗水量

運(yùn)用2014－2015年試驗(yàn)數(shù)據(jù)對(duì)擬合得到臨界氮濃度稀釋曲線和圖3中的經(jīng)驗(yàn)公式的可行性進(jìn)行驗(yàn)證，其步驟為：通過(guò)2014－2015年不同水分處理下番茄全生育期平均日耗水量、全生育期地上部生物量、式（3）和圖3中的經(jīng)驗(yàn)公式擬合得到的臨界氮濃度稀釋計(jì)算臨界氮含量模擬值，將臨界氮濃度模擬值與計(jì)算值進(jìn)行比較，結(jié)果如圖2b所示。由圖2b可以看出通過(guò)番茄生育期平均日耗水量、全生育期地上部生物量、式（3）和圖3中的經(jīng)驗(yàn)公式模擬得到的臨界氮濃度與計(jì)算值之間有較好的一致性，其MAE為0.13～0.34 g/(100 g)，RMSE為0.14～0.39 g/(100 g)，2為0.94～0.99，因此可以運(yùn)用不同灌水處理氮濃度稀釋曲線和番茄全生育期平均日耗水量模擬番茄臨界氮濃度，提高了氮濃度稀釋曲線模型在溫室番茄的適用性。

2.3 不同水分處理下番茄氮素營(yíng)養(yǎng)指數(shù)變化分析

氮素營(yíng)養(yǎng)指數(shù)（NNI）作為實(shí)際植株含氮量與臨界含氮量的比值，可以直觀的反映作物體內(nèi)氮素的營(yíng)養(yǎng)狀況，當(dāng)NNI=1，氮素營(yíng)養(yǎng)狀況最為適宜；當(dāng)NNI>1，表現(xiàn)為氮素營(yíng)養(yǎng)過(guò)剩；當(dāng)NNI<1，表現(xiàn)為氮素營(yíng)養(yǎng)虧缺。圖4為根據(jù)公式（4）計(jì)算所得的不同水分和氮素處理下番茄氮營(yíng)養(yǎng)指數(shù)（NNI）的動(dòng)態(tài)變化。由圖4可以看出不同灌水量時(shí)，NNI隨施氮量的增加而增加，其值范圍為0.66～1.19。定植后70 d后番茄開(kāi)始開(kāi)花，此時(shí)營(yíng)養(yǎng)生長(zhǎng)和生殖生長(zhǎng)旺盛，植株對(duì)氮素的需求量較大，不同處理NNI的差距增大，此時(shí)開(kāi)始追施氮肥來(lái)滿足植株氮素的需求，充分灌水處理（W1）和僅苗期虧水（W2）時(shí)，N0和N150處理全生育期NNI始終小于1，不能滿足植株對(duì)氮素的需求，N300處理NNI在1附近波動(dòng)，因此充分灌水和僅苗期虧水適宜的施氮量為N300；苗期開(kāi)花期連續(xù)虧水（W3）和全生育期虧水（W4）時(shí)，N0處理全生育期NNI均小于1，隨著追施氮素的增加，N150和N300在定植90 d后NNI均大于1，且N150處理的NNI較為接近1，表明苗期開(kāi)花期連續(xù)虧水和全生育期虧水處理時(shí)，N0和N300處理因?yàn)榈夭蛔慊蛘叩靥鄷?huì)抑制植株的生長(zhǎng)，最適宜的施氮量為N150。

圖4 2013－2015年不同水氮處理番茄氮素營(yíng)養(yǎng)指標(biāo)（NNI）動(dòng)態(tài)變化

2.4 不同水氮處理番茄葉片SPAD值動(dòng)態(tài)變化

不同水氮處理下番茄不同葉位SPAD值隨定植后天數(shù)變化如表4所示。由表4可以得知，番茄不同葉位SPAD均值隨定植天數(shù)表現(xiàn)出先增加后減小的變化趨勢(shì)，在DAT=110 d時(shí)各處理不同葉位SPAD值均達(dá)到最大值，分別為42.8～70.3，隨后番茄葉片開(kāi)始老化，不同葉位SPAD值開(kāi)始下降。

表4 2013－2015年不同水氮處理番茄不同葉位SPAD值

充分灌水（W1）時(shí)，番茄不同葉位葉SPAD值變化范圍為21.0～70.3，番茄不同葉位SPAD值隨施氮量的增加呈顯著增加的趨勢(shì)；苗期虧水（W2）、苗期開(kāi)花期虧水（W3）和全生育期虧水（W4）時(shí)，虧水時(shí)間內(nèi)番茄不同葉位SPAD值隨施氮量的增加呈顯著增加的趨勢(shì)，且N150和N300處理之間不同葉位SPAD值無(wú)顯著性差異；與充分灌水處理相比，僅苗期虧水（W2）時(shí)不會(huì)顯著影響復(fù)水后番茄不同葉位SPAD值，隨著虧水天數(shù)的增加，苗期開(kāi)花期連續(xù)虧水成熟期復(fù)水后（W3）不同葉位SPAD值仍顯著低于充分灌水處理（W1）；全生育期虧水（W4）時(shí)，不同生育期葉片不同葉位SPAD值均顯著低于充分灌水處理。表明充分灌水時(shí)，增加施氮量有顯著的增加不同葉位SPAD值的效應(yīng)，而非充分灌水時(shí)，施氮量增加到一定量（150 kg/hm2）繼續(xù)增加施氮時(shí)不會(huì)顯著增加番茄不同葉位SPAD值。施氮量相同時(shí)，充分灌水處理番茄不同葉位SPAD值顯著高于非充分灌水處理番茄不同葉位SPAD值，表明增加灌水量可以增加番茄SPAD值，且各處理中位葉SPAD值高于上位葉和下位葉SPAD值。

2.5 不同葉位SPAD和氮營(yíng)養(yǎng)指數(shù)（NNI）之間的關(guān)系

圖5為2013－2015年不同水氮處理不同葉位（上、中、下葉位）SPAD值與對(duì)應(yīng)處理的NNI之間的相關(guān)關(guān)系。由圖5可以看出，不同水氮處理上位葉和下位葉SPAD值和NNI之間決定系數(shù)（2）為0.05～0.93，上位葉和下位葉SPAD值和NNI之間的回歸關(guān)系受不同處理影響較大，除W3處理DAT=90d外，SPAD值與之間無(wú)顯著線性關(guān)系；而中位葉SPAD值和NNI之間有良好的相關(guān)性，決定系數(shù)（2）為0.77～0.98，除W3和W4處理DAT=30 d外，其他均達(dá)到顯著水平以上，表明番茄中位葉SPAD值和NNI之間呈顯著正相關(guān)關(guān)系，NNI隨著SPAD的增加而增加，且擬合關(guān)系表現(xiàn)比較穩(wěn)定、密切，因此在不同水分處理時(shí)可以利用中位葉SPAD值對(duì)NNI進(jìn)行估算，這樣結(jié)合SPAD值監(jiān)測(cè)方便、快捷和NNI預(yù)測(cè)精度高的優(yōu)點(diǎn)可以對(duì)番茄植株氮含量進(jìn)行更加快捷精準(zhǔn)的估算。

圖5 2013－2015年番茄不同葉位SPAD值與NNI之間的關(guān)系

將2013－2015年不同水分處理不同定植后天數(shù)番茄NNI和中位葉SPAD值進(jìn)行擬合（圖5），由圖5可以看出不同水分處理番茄NNI和中位葉SPAD值的相關(guān)關(guān)系不同，其擬合的形式為

NNI=·SPAD+（5）

式中NNI為不同處理氮素營(yíng)養(yǎng)指數(shù)；SPAD為不同處理中位葉SPAD值；、為擬合公式參數(shù)，不同水分處理擬合得到的參數(shù)、如圖5所示。

不同水分處理番茄中位葉SPAD與NNI之間的關(guān)系不同，因此采用2013－2015年試驗(yàn)數(shù)據(jù)對(duì)不同水分處理不同定植后天數(shù)番茄平均日耗水量與參數(shù)、進(jìn)行擬合，得到參數(shù)、與番茄日耗水量（表3）的擬合關(guān)系，如圖6所示。由圖6可以看出番茄中位葉SPAD與NNI之間擬合關(guān)系參數(shù)、與對(duì)應(yīng)番茄日耗水量有顯著的關(guān)系，因此可以采用不同處理日耗水量得到參數(shù)、，進(jìn)而得到番茄中位葉SPAD與NNI的線性關(guān)系，進(jìn)而根據(jù)番茄中位葉SPAD值估算得到番茄NNI進(jìn)行氮營(yíng)養(yǎng)診斷。將NNI=1作為適宜施氮量的標(biāo)準(zhǔn)，NNI＞1或NNI＜1，表明氮素施加過(guò)量或不足營(yíng)養(yǎng)，因此運(yùn)用該方法模擬得到不同水分處理的NNI，進(jìn)一步為氮營(yíng)養(yǎng)診斷提供指導(dǎo)。

注：m和k為不同水分處理經(jīng)驗(yàn)系數(shù)；ETi為不同水分處理番茄i天內(nèi)平均日耗水量。

3 討 論

3.1 臨界氮濃度稀釋曲線模型在西北地區(qū)溫室番茄的適用性

水分和氮素作為影響植株生長(zhǎng)的重要因素，不同水分和氮素處理影響作物發(fā)育和干物質(zhì)累計(jì)[27]，進(jìn)而影響植株對(duì)氮素的吸收[28]。楊慧等[12]在研究盆栽番茄表明不同水分處理下臨界施氮量有所不同，高水處理下植株的臨界氮濃度較大，灌水可以促進(jìn)植株對(duì)氮素的吸收。本研究也表明增加灌水量和施氮量可以增加番茄植株的臨界氮濃度，通過(guò)擬合臨界氮稀釋曲線得到不同水分處理番茄臨界氮稀釋曲線參數(shù)和參數(shù)。不同因素對(duì)臨界氮稀釋曲線參數(shù)和參數(shù)的影響有所不同，趙犇等[33]研究表明由于品種不同植株吸收和同化氮的能力不同，參數(shù)與品種蛋白質(zhì)含量呈正相關(guān)關(guān)系；向友珍等[34]研究氮濃度稀釋模型在溫室辣椒的適用性，得到臨界氮稀釋曲線參數(shù)與灌水量呈先增加后降低的趨勢(shì)，參數(shù)不隨灌水量的變化而顯著變化；強(qiáng)生才等[16]在研究大田作物時(shí)表明，參數(shù)具有較好的穩(wěn)定性，不會(huì)隨降雨年型的改變而改變，且降水量較多時(shí)可以增加植株對(duì)土壤氮素的利用，從而減緩了植株氮含量的稀釋過(guò)程，最終導(dǎo)致參數(shù)明顯偏小。不同作物品種對(duì)水分和氮素的敏感性不同，溫室番茄和辣椒對(duì)水分和氮素的敏感性遠(yuǎn)遠(yuǎn)高于大田玉米，水分和氮素影響著植株對(duì)氮素的吸收和植株地上部的生長(zhǎng)，進(jìn)而影響植株氮素的吸收和稀釋。本研究表明臨界氮稀釋曲線參數(shù)和參數(shù)均隨番茄全生育期平均日耗水量增加呈現(xiàn)出先降低后增加，番茄全生育期耗水量較低時(shí)植株吸收氮素的速率高于植株生長(zhǎng)的速率，減緩了植株氮含量的稀釋速率，降低了植株臨界氮濃度值；當(dāng)耗水量高于某一值時(shí)，繼續(xù)增加耗水量時(shí)，植株吸收氮素的速率低于植株生長(zhǎng)的速率，加快了植株氮含量的稀釋速率，增大植株臨界氮濃度值，提高番茄植株對(duì)氮素營(yíng)養(yǎng)的容納能力，促進(jìn)植株生長(zhǎng)。本文建立的基于全生育期平均日耗水量的不同水分處理臨界氮稀釋模型可以較為有效的估算出不同水分處理下植株臨界氮濃度，與計(jì)算值相比，其絕對(duì)誤差（MAE）為0.13～0.34 g/(100 g)，標(biāo)準(zhǔn)誤差（RMSE）為0.14～0.39 g/(100 g)，絕對(duì)系數(shù)（2）為0.94～0.99，因此可以通過(guò)該方法準(zhǔn)確估算不同水分處理的植株臨界氮濃度，提高了臨界氮稀釋曲線在溫室番茄不同水分處理的適用性。

3.2 基于不同葉位SPAD值的氮營(yíng)養(yǎng)指數(shù)估算方法

SPAD儀具有攜帶方便、可實(shí)時(shí)實(shí)地進(jìn)行測(cè)定等優(yōu)點(diǎn)，被普遍用于監(jiān)測(cè)棉花[35]、小麥[36]和玉米[37]等多種作物的施氮水平。葉片SPAD值與NNI之間的相關(guān)關(guān)系在玉米[38]和小麥[39]等植物上被證實(shí)，且Ziadi等[38]和 Prost and Jeuffroy[39]研究表明玉米和小麥葉片相對(duì)SPAD值與NNI之間的非線性關(guān)系，Debaeke等[36]研究表明小麥葉片相對(duì)SPAD值與NNI之間的非線性關(guān)系受年份、品種和生育時(shí)期影響不顯著。Yang 等[23]研究表明作物不同葉位SPAD值和NNI表現(xiàn)出不同的線性關(guān)系，且不同葉位的穩(wěn)定性不同，其中某一葉位葉片SPAD值和NNI擬合度較高，穩(wěn)定性較好，可以作為診斷氮素的理想葉位。本研究也表明番茄不同水分處理不同葉位葉片SPAD和NNI之間呈線性相關(guān)關(guān)系，其中上位葉片和下位葉片SPAD值與NNI之間擬合度較差，擬合度受年份和處理的影響較大，穩(wěn)定性較差，而中位葉片SPAD值和NNI之間呈顯著線性正相關(guān)關(guān)系，穩(wěn)定性較好，因此可以將番茄中位葉片作為適宜的診斷葉片進(jìn)行氮素診斷。且番茄中位葉片SPAD值和NNI之間擬合參數(shù)與番茄平均耗水量呈顯著非線性相關(guān)關(guān)系，因此通過(guò)不同水分處理日耗水量得到番茄中位葉片SPAD值與NNI之間的線性關(guān)系，并以NNI=1作為理想施氮營(yíng)養(yǎng)狀態(tài)指標(biāo)，根據(jù)番茄不同耗水量得到NNI=1時(shí)的SPAD值，并以此SPAD值作為適宜值進(jìn)行氮素營(yíng)養(yǎng)診斷。當(dāng)試驗(yàn)處理中位葉片SPAD值大于適宜SPAD值時(shí)，表明該處理氮素施加過(guò)量應(yīng)該適當(dāng)減少，當(dāng)試驗(yàn)處理中位葉片SPAD值小于適宜SPAD值時(shí)，表明該處理氮素施加較少應(yīng)該適當(dāng)增加施氮量。因此可以通過(guò)比較不同處理SPAD值與適宜SPAD值來(lái)準(zhǔn)確判斷氮營(yíng)養(yǎng)狀況，來(lái)及時(shí)調(diào)整施氮量。

4 結(jié) 論

本研究依據(jù)2 a溫室番茄不同水分處理下3個(gè)氮素水平的試驗(yàn)數(shù)據(jù)，率定和驗(yàn)證臨界氮濃度稀釋曲線模型在西北地區(qū)溫室番茄不同水分處理下的適用性，并利用不同葉位葉片SPAD值對(duì)溫室番茄氮營(yíng)養(yǎng)指數(shù)進(jìn)行估算，并進(jìn)一步進(jìn)行氮營(yíng)養(yǎng)診斷，結(jié)果表明：

1）番茄臨界氮濃度與采樣日地上部最大生物量之間符合冪指數(shù)關(guān)系（N=DW-b），其中參數(shù)為2.84～3.44，參數(shù)為0.22～0.35，參數(shù)、與番茄全生育期平均日耗水量之間存在顯著的非線性相關(guān)關(guān)系，因此可以通過(guò)不同水分處理全生育期平均日耗水量估算參數(shù)和參數(shù)，進(jìn)而準(zhǔn)確估算不同水分處理的植株臨界氮濃度，提高了臨界氮稀釋曲線在西北地區(qū)溫室番茄不同水分處理的適用性。

2）利用不同葉片SPAD值進(jìn)行番茄氮營(yíng)養(yǎng)診斷表明，氮營(yíng)養(yǎng)診斷理想葉位是中位葉片。番茄中位葉片SPAD值與氮素營(yíng)養(yǎng)指數(shù)NNI有良好的線性關(guān)系（NNI=·SPAD+），且參數(shù)、與番茄日耗水量有顯著的非線性相關(guān)關(guān)系，因此可以通過(guò)不同水分處理日耗水量估算出參數(shù)、，進(jìn)而模擬出不同定植后天數(shù)NNI=1時(shí)的SPAD值，并以此SPAD值作為判斷氮素營(yíng)養(yǎng)狀況的標(biāo)準(zhǔn)來(lái)及時(shí)精準(zhǔn)診斷植株氮營(yíng)養(yǎng)狀態(tài)。

[1] 高兵，任濤，李俊良，等. 灌溉策略及氮肥施用對(duì)設(shè)施番茄產(chǎn)量及氮素利用的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào)，2008，14(6)：1104－1109.

Gao Bing, Ren Tao, Li Junliang, et al. Effects of irrigation strategies and N sidedressing on the yield and N utilization of greenhouse tomato[J]. Plant Nutrition and Fertilizer Science, 2008, 14(6): 1104－1109. (in Chinese with English abstract)

[2] 郭文忠，劉聲峰，李丁仁，等. 設(shè)施蔬菜土壤次生鹽漬化發(fā)生機(jī)理的研究現(xiàn)狀與展望[J]. 土壤，2004，36(1)：25－29.

Guo Wenzhong, Liu Shengfeng, Li Dingren, et al. Mechanism of soil salinization in protected cultivation[J]. Soil, 2004, 36(1): 25－29. (in Chinese with English abstract)

[3] Ziadi N, Brassard M, Bélanger G, et al. Critical nitrogen curve and nitrogen nutrition index for corn in eastern canada[J]. Agron J, 2008, 100(2): 271－276.

[4] Greenwood D J, Lemaire G. Decline in percentage N of C3and C4crops with increasing plant mass[J]. Annals of Botany, 1990, 66(4): 425－436.

[5] Lemaire G, Gastal F. Relationships between plant-N, plant mass and relative growth rate for C3and C4crops[C]// Proceedings first ESA Congress, Paris, 1990, 1.

[6] Lemaire G, van Oosterom E, Sheehy J, et al. Is crop N demand more closely related to dry matter accumulation or leaf area expansion during vegetative growth[J]. Field Crops Research, 2007, 100(1): 91－106.

[7] 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.

[8] 王子勝，金路路，趙文青，等. 東北特早熟棉區(qū)不同群體棉花氮臨界濃度稀釋模型的建立初探[J]. 棉花學(xué)報(bào)，2012，24(5)：427－434.

Wang Zisheng, Jin Lulu, Zhao Wenqing, et al. Preliminary development and application of a critical nitrogen concentration dilution model for cotton grown in northeast China[J]. Cotton Science, 2012, 24(5): 427－434. (in Chinese with English abstract)

[9] J?rgen E. Olesen, J?rgen Berntsen. Crop nitrogen demand and canopy area expansion in winter wheat during vegetative growth[J]. European Journal of Agronomy, 2002, 16(4): 279－294.

[10] 強(qiáng)生才，張富倉(cāng)，田建柯，等. 基于葉片干物質(zhì)的冬小麥臨界氮稀釋曲線模擬研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2015，46(11)：121－128.

Qiang Shengcai, Zhang Fucang, Tian Jianke, et al. development of critical nitrogen dilution curve in winter wheat based on leaf dry matter[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(11): 121－128. (in Chinese with English abstract).

[11] 王新，馬富裕，刁明，等. 滴灌番茄臨界氮濃度、氮素吸收和氮營(yíng)養(yǎng)指數(shù)模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(18)：99－108.

Wang Xin, Ma Fuyu, Diao Ming, et al. Simulation of critical nitrogen concentration, nitrogen uptake and nitrogen nutrition index of processing tomato with drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(18):99－108. (in Chinese with English abstract).

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

Yang Hui, Cao Hongxia, Liu Meiyu, et al. Simulation of critical nitrogen concentration and nutrition index of tomato under different water and nitrogen conditions[J]. Journal of Plant Nutrition and Fertilizer. 2015, 21(5): 1234－1242. (in Chinese with English abstract).

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

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

[15] 李正鵬，宋明丹，馮浩. 關(guān)中地區(qū)玉米臨界氮濃度稀釋曲線的建立和驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(13)：135－141.

Li Zhengpeng, Song Mingdan, Feng Hao. Development and validation of critical nitrogen content curve for maize in Guanzhong area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015. 31(13): 135－141. (in Chinese with English abstract).

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

Qiang Shengcai，Zhang Fucang，Xiang Youzhen，et al. Simulation and validation of critical nitrogen dilution curve for summer maize in Guanzhong Plain during different rainfall years[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(17): 168－175. (in Chinese with English abstract).

[17] 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[J]. European Journal of Agronomy, 2008, 28(4): 614－624.

[18] Lemaire G, Gastal F, Salette J. N uptake and distributionin plant canopies: Diagnosis of the nitrogen status in crops[J]. Berlin: Spriger. 1997.

[19] 賈良良，范明生，張福鎖，等. 應(yīng)用數(shù)碼相機(jī)進(jìn)行水稻氮營(yíng)養(yǎng)診斷[J]. 光譜學(xué)與光譜分析，2009，29(8)：2176－2179.

Jia Liangliang, Fan Mingsheng, Zhang Fusuo, et al. Nitrogen status diagnosis of rice by using a digital camera[J]. Spectroscopy and Spectral Analysis, 2009, 29(8): 2176－2179. (in Chinese with English abstract).

[20] Madeire A C, Mendonca A, Ferreire M E, et al. Relationship between spectroradiometric and chlorophyll measurements in green beans[J]. Communications in Soil Science & Plant Analysis, 2000, 7(2/3): 695－704.

[21] 張金恒，王柯，王人潮. 葉綠素計(jì)SPAD-502在水稻氮素營(yíng)養(yǎng)診斷中的應(yīng)用[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2003，31(2)：177－180.

Zhang Jinheng, Wang Ke,Wang Renchao. Application o f chlorophyll meter SPAD-502 in diagnosis of nitrogen status and nitrogenous fertilizer in rice[J]. Jour. of Northwest Sci-Tech Univ. of Agri. and For.: Nat. Sci. Ed. 2003, 31(2): 177－180. (in Chinese with English abstract).

[22] Errecart P M, Agnusdei M G, Lattanzi F A, et al. Leaf nitrogen concentration and chlorophyll meter readings as predictors of tall fescue nitrogen nutrition status[J]. Field Crops Research, 2012,129(1): 46－58.

[23] Yang H, Yang J, Lv Y, et al. SPAD values and nitrogen nutrition index for the evaluation of rice nitrogen status[J]. Plant Production Science, 2014, 17(1): 81－92.

[24] 王峰，杜太生，邱讓建，等. 虧缺灌溉對(duì)溫室番茄產(chǎn)量與水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(9)：46－52.

Wang Feng, Du Taisheng, Qiu Rangjian, et al. Effects of deficit irrigation on yield and water use efficiency of tomato in solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(9): 46－52. (in Chinese with English abstract)

[25] 石小虎. 溫室膜下滴灌番茄對(duì)水氮耦合的響應(yīng)研究[D]. 楊陵：西北農(nóng)林科技大學(xué)，2013.

Shi Xiaohu. Response of Greenhouse Tomato with Drip Irrigation under Plastic Mulch to Water-nitrogen Coupling[D]. Yangling: Northwest A＆F University, 2013. (in Chinese with English abstract)

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

[27] 薛曉萍，陳兵林，郭文琦. 棉花臨界需氮量動(dòng)態(tài)定量模型[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2006，17(12)：2363－2370.

Xue Xiaoping, Chen Binglin, Guo Wenqi. Dynamic quantitative model of critical nitrogen demand of cotton[J]. Chinese Journal of Applied Ecology, 2006, 17(12): 2363－2370. (in Chinese with English abstract)

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

Liang Xiaogui, Zhang Jingting, Zhou Lili, 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 with English abstract).

[29] Lemaire G, Salette J. Croissance estivale en matiere seche de peuplements de fetuque elevee (Festuca arundinacea Schreb) et de dactyle (Dactylis glomerata L.) danslouest de la France: I. Etude en conditions de nutrition azotee et dalimentation hydrique non limitantes[J]. Agronomie, 1987, 7(6): 373－380.

[30] Nash J E, Sutcliffe J V. River flow forecasting through conceptual models, part I-A discussion of principles[J]. Journal of Hydrology, 1970, 10: 282－290.

[31] Willmott C J, Matsuura K. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance[J]. Climate Research, 2005, 30(1): 79－82.

[32] Wu Fuzhong, Bao Weikai, Li Fanglan, et al. Effects of drought stress and N supply on the growth, biomass partitioning and water-use efficiency of Sophora davidii seedlings[J]. Environ Exp Bot, 2008, 63(1/2/3): 248－255.

[33] 趙犇，姚霞，田永超，等. 基于臨界氮濃度的小麥地上部氮虧缺模型[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2012，23(11)：3141－3148.

Zhao Ben, Yao Xia, Tian Yongchao, et al. Accumulative nitrogen deficit models of wheat aboveground part based on critical nitrogen concentration[J]. Chinese Journal of Applied Ecology, 2012, 23(11): 3141－3148. (in Chinese with English abstract).

[34] 向友珍，張富倉(cāng)，范軍亮，等. 基于臨界氮濃度模型的日光溫室甜椒氮營(yíng)養(yǎng)診斷[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32（17）：89－97.

Xiang Youzhen, Zhang Fucang, Fan Junliang, et al. Nutrition diagnosis for N in bell pepper based on critical nitrogen model in solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(17): 89－97. (in Chinese with English abstract)

[35] Wu F B, Wu L H, Xu F H. Chlorophyll meter to predict nitrogen sidedress requirements for short-season cotton (Gossypium hirsutum L.)[J]. Field Crops Res. 1998, 56:309－314.

[36] Debaeke P, Rouet P, Justes E. Relationship between the normalized SPAD index and the nitrogen nutrition index: Application to durum wheat[J]. Journal of Plant Nutrition. 2006, 29: 75－92.

[37] Singh V, Singh Y, Singh B, et al. Calibrating the leaf colour chart for need based fertilizer nitrogen management in different maize (Zea mays L.) genotypes[J]. Field Crops Res. 2011, 120: 276－282.

[38] Ziadi N, Brassard M G, Belanger A N, et al. Critical nitrogen curve and nitrogen nutrition index for corn in eastern Canada[J]. Agron. J.2008, 100: 271－276.

[39] Prost L, Jeuffroy M H. 2007. Replacing the nitrogen nutrition index by the chlorophyll meter to assess wheat Nstatus[J]. Agronomy Sustainable Development. 2007, 27: 321－330.

Estimation of nitrogen nutrition index of greenhouse tomato under different water and nitrogen fertilizer treatments based on leaf SPAD

Shi Xiaohu1,2, Cai Huanjie2※

(1.,030024,; 2.,,712100,)

Water and nitrogen have great influences on the production of tomato in greenhouses. This study aims to determine the applicability of the critical nitrogen concentration dilution curve model in different water conditions, and to suggest a rapid and accurate method of identifying the nitrogen condition of tomato. Greenhouse experiments with different water-nitrogen combinations were conducted in 2013-2015 in Yangling, Shaanxi. Four water treatments were investigated, i.e., full irrigation through the whole life span of tomato, deficit of 50% at seedling stage, deficit of 50% at flowering stage, and deficit of 50% through the whole life span. Three nitrogen treatments were investigated, i.e., 0, 150, and 300 kg/hm2. We validated and verified the critical nitrogen concentration dilution curve model with the data in 2013-2015. We verified the applicability of the model under different water treatments of greenhouse tomato in Northwest China. We built relationship between the model parameters and the average daily water consumption of tomato during the whole growth stage, and improved the applicability of the critical nitrogen concentration dilution curve model under different water treatments. Both the parameter a and the parameter b of the critical nitrogen concentration dilution curve model decreased first and then increased with the increase of the average daily water consumption during the whole growth period of the tomato.We analyzed the relationship between SPAD (soil and plant analyzer development) values and nitrogen nutrition index (NNI) at different leaf positions with the test data in 2013-2015. There was a linear correlation between SPAD and NNI at different leaf positions of tomato under different water treatments. The fitting relationship between the SPAD values and the NNI of the upper and lower leaves was poor, and the degree of fitting was greatly affected by the year and the treatment. The fitting relationship between the median leaf’s SPAD value and NNI was a significant linear positive correlation relationship and had good stability. Therefore, the tomato median leaf can be used as a suitable diagnostic leaf for nitrogen diagnosis. The results showed that: 1) The model could well estimate the critical nitrogen concentration of tomato with the absolute error of 0.13-0.34 g/ 100 g, the standard error of 0.14-0.39 g/ 100 g, and the decisive factor (2) of 0.94-0.99. 2) At the median position of canopy, there was a good linear correlation between the SPAD value and NNI, and the SAPD values of the middle leaves can be used to determine the nitrogen conditions of tomato. 3) The SPAD values of middle leaves at the NNI of 1 can be the indicators for appropriate nitrogen fertilizing. When the SPAD value of the median leaf in the different treatments is greater than the appropriate SPAD value, it indicates that the application of nitrogen in the treatment is excessive and should be appropriately reduced. When the SPAD value of the medial leaf in the test treatment is less than the appropriate SPAD value, it indicates that the application of nitrogen should be appropriately increased. Therefore, we have determined and validated the applicability of the critical nitrogen concentration dilution curve model under different water treatments of greenhouse tomato in Northwest China and constructed a nitrogen nutrition diagnosis method based on the daily water consumption of tomato, which can support nitrogen diagnosis and nitrogen management of tomato in greenhouses in this region.

chlorophyll; water; greenhouse; tomato; water treatments; critical nitrogen concentration; nitrogen nutrition index

10.11975/j.issn.1002-6819.2018.17.016

S152.7

A

1002-6819(2018)-17-0116-11

2018-03-24

2018-07-27

國(guó)家科技支撐計(jì)劃（2011BAD29B01）；國(guó)家自然基金（51179162）；2011年度高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金資助課題（20110204130004）

石小虎，博士，主要從事農(nóng)業(yè)節(jié)水理論研究。 Email：shixiaohu2006@126.com

蔡煥杰，教授，博士生導(dǎo)師，主要從事農(nóng)業(yè)節(jié)水與水資源高效利用研究。Email：caih j@nwsuaf.edu.cn

石小虎，蔡煥杰. 基于葉片SPAD估算不同水氮處理下溫室番茄氮營(yíng)養(yǎng)指數(shù)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(17)：116－126. doi：10.11975/j.issn.1002-6819.2018.17.016 http://www.tcsae.org

Shi Xiaohu, Cai Huanjie. Estimation of nitrogen nutrition index of greenhouse tomato under different water and nitrogen fertilizer treatments based on leaf SPAD[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(17): 116－126. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.17.016 http://www.tcsae.org