肖蕊

（長治職業(yè)技術(shù)學(xué)院園藝系，山西長治 046011）

辣椒（Capsicum annumL.）在我國各地均有不同程度的分布和栽培，是人們?nèi)粘Ｉ钪凶钍軞g迎的蔬菜作物之一（王繼榜，2013）。氮素是影響蔬菜作物產(chǎn)量和品質(zhì)的關(guān)鍵性元素（張妮 等，2021；袁巧麗 等，2022），在農(nóng)業(yè)生產(chǎn)中，大多數(shù)種植戶為了追求經(jīng)濟(jì)效益最大化，不是按照作物需肥特性和土壤供肥性能來精準(zhǔn)施肥，而是普遍采用過量施入氮肥來提質(zhì)增產(chǎn)（聶金 等，2021；王赫 等，2021），從而導(dǎo)致土壤肥力不平衡、氮肥利用效率低，威脅農(nóng)田生態(tài)環(huán)境，制約農(nóng)業(yè)可持續(xù)發(fā)展（吳玥 等，2021）。因此，尋找蔬菜作物產(chǎn)量最大化下的臨界氮濃度值，對實現(xiàn)科學(xué)合理施肥具有重要意義。

優(yōu)化作物生長階段所需的氮濃度對于評價作物氮素營養(yǎng)狀況至關(guān)重要（王濤 等，2013）。前人針對測土配方施肥（Hansen &Schjoerring，2003）、SPAD 計快速診斷（Zheng et al.，2015）、光譜遙感（Ren et al.，2010）和機(jī)器視覺（賈彪和馬富裕，2016）等方法對植株氮素營養(yǎng)精確診斷做了大量研究，但是這些作物氮素營養(yǎng)診斷技術(shù)受成本、技術(shù)普及度的影響，存在地域和年際間的不穩(wěn)定性，推廣存在一定困難。Greenwood 等（1991）總結(jié)了作物生長和氮素吸收的規(guī)律，提出了臨界氮濃度（Nc）的概念，即作物最大生長所需的最低氮濃度。Nc因其在作物氮素營養(yǎng)診斷中的準(zhǔn)確性和穩(wěn)定性而受到研究者的廣泛關(guān)注，因此可以作為評價作物氮素營養(yǎng)是否足夠的優(yōu)選技術(shù)?；谧魑镎晟锪浚▍稳銤?等，2018；石小虎和蔡煥杰，2018；劉秋霞等，2019；付江鵬 等，2020）和器官生物量（Zhao et al.，2017）的臨界氮濃度稀釋曲線已在多種農(nóng)作物上構(gòu)建和應(yīng)用。向友珍等（2016）、Hoogmoed 和Sadras（2016）研究表明，由于氣候環(huán)境差異、生育期長短和品種差異的影響，模型參數(shù)存在變異性。因此，對模型的本地化研究是非常重要的。目前眾多學(xué)者對辣椒開展了氮肥試驗研究，唐恒朋（2016）研究了不同氮素用量對辣椒的形態(tài)指標(biāo)、光合特性、產(chǎn)量和品質(zhì)的影響，結(jié)果表明在貴州黔南地區(qū)長辣6 號和火焰山氮素施用量以207 kg·hm-2較優(yōu)，實際生產(chǎn)中可在此基礎(chǔ)上適當(dāng)調(diào)整；葉潔（2017）研究了控釋氮肥用量對辣椒生長生理和養(yǎng)分利用的影響，結(jié)果表明常規(guī)控釋氮肥用量（172.2 kg·hm-2）可作為辣椒栽培中控釋氮肥用量的參考依據(jù)。辣椒作為貴州省的特色優(yōu)勢經(jīng)濟(jì)作物，在全國辣椒產(chǎn)業(yè)中具有舉足輕重的地位。由于不同基因型和環(huán)境的變化對Nc稀釋曲線參數(shù)的響應(yīng)不同，本試驗以當(dāng)?shù)刂髟缘? 個品種辣研102、蘇椒5 號為研究對象開展辣椒氮肥試驗，建立當(dāng)?shù)刂髟岳苯菲贩NNc稀釋曲線，將其與現(xiàn)有不同作物品種的Nc曲線進(jìn)行比較，并對該曲線的可靠性進(jìn)行評價，以期為區(qū)域辣椒氮肥精準(zhǔn)管理及高產(chǎn)栽培提供科學(xué)參考。

1 材料與方法

1.1 試驗地概況

于2019—2020 年連續(xù)兩年在貴州省農(nóng)業(yè)科學(xué)院試驗基地（E 107°05′17″，N 26°16′33″）進(jìn)行4 個田間地塊定位試驗。試驗地海拔高度為1 000 m，年均溫度16 ℃，年均降雨量1 270 mm。土壤質(zhì)地為紅壤土，在播種前采集耕層0～30 cm土壤樣品，風(fēng)干后測定基本理化性質(zhì)（表1）。

表1 試驗地土壤基礎(chǔ)肥力

1.2 試驗設(shè)計

供試?yán)苯菲贩N為當(dāng)?shù)胤N植面積較廣的辣研102、蘇椒5 號。試驗設(shè)5 個氮肥施用水平，分別為0、70、140、210、280 kg·hm-2，以N0、N70、N140、N210、N280 表示；小區(qū)面積為10 m2，3 次重復(fù)，隨機(jī)區(qū)組排列；定植株距30 cm，行距50 cm，每小區(qū)種植60 株。氮肥為尿素（總N≥46.4%），磷肥為磷酸二氫鉀（138 kg ·hm-2），鉀肥為硫酸鉀（120 kg·hm-2），均作基肥一次性施入土壤。采用穴盤育苗法，待幼苗長至6～7片真葉、株高15 cm 左右時，帶土定植。

1.3 測定指標(biāo)及方法

在定植后30、50、70、90、110、130 d，每小區(qū)隨機(jī)選取長勢一致的植株5 株，將地上部分分離成莖、葉、果實3 部分，采用干燥法對各器官干質(zhì)量進(jìn)行測定，最后折算成辣椒植株地上部干質(zhì)量。分別將各處理的干樣粉碎、研磨、過篩，采用微量凱氏定氮法測定各器官氮含量，折算成辣椒植株地上部氮濃度。于第一次采收果實到最后一次采收果實，分別統(tǒng)計各批次小區(qū)產(chǎn)量，最后計算辣椒總產(chǎn)量。

1.4 臨界氮濃度稀釋曲線模型描述

根據(jù)Justes 等（1994）提出的臨界氮濃度稀釋曲線計算方法，建模步驟如下：①對每次取樣的植株地上部干質(zhì)量及對應(yīng)的氮濃度進(jìn)行方差分析，將其分為限氮營養(yǎng)和非限氮營養(yǎng)2 類；②對受氮素影響的氮素水平，將其植株地上部干質(zhì)量與氮濃度進(jìn)行曲線擬合；③對不受氮素影響的氮素水平，取其植株地上部干質(zhì)量的均值代表最大干質(zhì)量；④采樣日的臨界氮濃度值由上述線性曲線與以最大植株地上部干質(zhì)量為橫坐標(biāo)的垂線交點的縱坐標(biāo)決定?；谥仓甑厣喜扛少|(zhì)量的臨界氮濃度稀釋曲線模型為：

Nc＝aPDM－b

式中，Nc為臨界氮濃度值（%），PDM為植株地上部干質(zhì)量的最大值（t·hm-2），a和b均為模型的參數(shù)。

1.5 氮素營養(yǎng)指數(shù)測定

氮素營養(yǎng)指數(shù)（nitrogen nutrition index，NNI）根據(jù)Yue 等（2014）描述的方法計算。

式中，Nc為臨界氮濃度值（%），PNC為植株地上部氮濃度（%）。若NNI＜1，表明氮素不足；NNI＝1，表明氮素恰好適量；NNI＞1，表明氮素過量。

1.6 數(shù)據(jù)處理與分析

采用Excel 2010 軟件進(jìn)行數(shù)據(jù)整理，采用IBM SPSS Statistics 22.0 軟件對PDM和PNC進(jìn)行單因素方差分析（one-way analysis of variance，ANOVA），預(yù)設(shè)顯著性水平P＜0.05，采用最小顯著性差異法檢驗PDM和PNC的差異顯著性。采用2019 年的數(shù)據(jù)，構(gòu)建Nc稀釋曲線；利用2020 年的數(shù)據(jù)，采用均方根誤差RMSE（root mean square error）和標(biāo)準(zhǔn)化均方根誤差（n-RMSE）進(jìn)行模型驗證（Yao et al.，2014）。

2 結(jié)果與分析

2.1 辣椒植株地上部干質(zhì)量及氮濃度動態(tài)變化

施氮水平對辣椒植株地上部干物質(zhì)積累有顯著影響，植株地上部干質(zhì)量大多呈隨著施氮量的增加逐漸增大的變化趨勢，且隨著植株生育進(jìn)程的推進(jìn)逐漸增加，而除2019 年定植50 d 的辣研102外，N210 與N280 處理的植株地上部干質(zhì)量均無顯著差異（表2）。辣研102、蘇椒5 號2019 年的植株地上部干質(zhì)量變化范圍分別為0.24～15.01 t·hm-2和0.64～14.22 t·hm-2，2020 年的變化范圍分別為0.34～15.14 t·hm-2和0.61～14.98 t·hm-2。2個品種植株地上部干質(zhì)量均滿足不等式：N0（P＝0.263）＜N70（P＝0.672）＜N140（P＝0.564）＜N210（P＝0.041）≈N280（P＝0.033）。

表2 辣椒植株地上部干質(zhì)量動態(tài)變化

同一取樣時期，2 個辣椒品種植株地上部氮濃度均隨著施氮量的提高而增加。從整個生育期來看，植株地上部氮濃度均隨著定植天數(shù)的增加呈下降趨勢（圖1）。2019 年2 個品種的植株地上部氮濃度變化范圍分別為1.38%～4.11%和1.39%～4.21%，2020 年的變化范圍分別為1.21%～4.39%和1.30%～4.31%。

圖1 辣椒植株地上部氮濃度動態(tài)變化

2.2 臨界氮濃度稀釋曲線模型構(gòu)建

利用2019 年試驗獲取的數(shù)據(jù)資料，計算2 個辣椒品種定植后30～130 d 的臨界氮濃度值。由于辣椒定植后30 d 各施氮水平植株地上部干質(zhì)量之間大多沒有顯著性差異（表2），導(dǎo)致植株地上部氮濃度變化較為穩(wěn)定，因此該時期的臨界氮濃度值由不受氮素限制處理的最小氮濃度與受氮素限制處理的最大氮濃度的平均值代表（向友珍 等，2016），經(jīng)計算辣研102、蘇椒5 號的臨界氮濃度常數(shù)為4.45%、4.38%。定植后50、70、90、110、130 d辣研102 的臨界氮濃度值分別為3.79%、2.88%、2.56%、2.04%和1.97%，蘇椒5 號分別為3.82%、3.03%、2.49%、2.01%和1.92%。對上述臨界氮濃度值與其對應(yīng)的最大植株地上部干質(zhì)量進(jìn)行擬合，得到了辣研102、蘇椒5 號2 個辣椒品種的臨界氮濃度稀釋曲線模型，方程的決定系數(shù)分別為0.94、0.95，均達(dá)到顯著水平（圖2）。

基于植株地上部干質(zhì)量構(gòu)建的2 個辣椒品種臨界氮濃度稀釋曲線模型中，參數(shù)a分別為4.67和4.23，參數(shù)b分別為-0.29 和-0.27（圖2）。為了進(jìn)一步分析這2 個品種之間的顯著性差異，首先將冪函數(shù)模型進(jìn)行直線化處理，即lnNc＝lna+blnPDM，辣研102、蘇椒5 號的直線化模型分別為lnNc＝1.58 -0.29 lnPDM和lnNc＝1.54-0.27 lnPDM。采用協(xié)方差分析方法，分別分析2 個辣椒品種臨界氮濃度稀釋曲線模型之間的斜率與截距間差異，結(jié)果顯示辣研102 和蘇椒5 號斜率與截距的P值分別為0.923 和0.535，都大于0.05，表明2 個辣椒品種之間沒有顯著性差異。因此，將這2 個品種的曲線并置擬合，形成統(tǒng)一的辣椒臨界氮濃度稀釋曲線模型：Nc＝4.39PDM-0.28，R2為0.93（圖3）。

圖2 不同品種辣椒臨界氮濃度稀釋曲線模型構(gòu)建

圖3 統(tǒng)一的辣椒臨界氮濃度稀釋曲線模型

2.3 臨界氮濃度稀釋曲線驗證

利用2020 年的獨(dú)立試驗數(shù)據(jù)資料對基于植株地上部干質(zhì)量的Nc模型進(jìn)行檢驗。結(jié)果表明（圖4），利用方程擬合的臨界氮濃度實測值與模擬值之間呈顯著的線性關(guān)系，決定系數(shù)R2為0.87。將獨(dú)立數(shù)據(jù)中的最大植株地上部干質(zhì)量代入模型后，計算臨界氮濃度模擬值，通過1∶1 直方圖來對比臨界氮濃度實測值和模擬值之間的擬合度。經(jīng)計算RMSE 為0.26，n-RMSE 為10.08%，表明模型具有較好的穩(wěn)定性，可用于辣椒氮素營養(yǎng)的評估與診斷。

圖4 辣椒臨界氮濃度稀釋曲線模型驗證結(jié)果

2.4 氮素營養(yǎng)指數(shù)動態(tài)變化

從圖5 可以看出，不同的施氮處理、生長時期、生育時期以及品種間的氮素營養(yǎng)指數(shù)（NNI）均有明顯差異。NNI 值均隨著施氮水平的提高而不斷上升，2019 年辣研102、蘇椒5 號的NNI 值變化范圍為0.61～1.30 和0.62～1.18，2020 年NNI 值變化范圍為0.60～1.19 和0.55～1.14。在辣椒整個生育期內(nèi)，N0、N70 和N140 處理下2 個品種的NNI 值均小于1，表明植株地上部氮濃度偏低，氮肥施用量不足；在N280 處理下NNI 值大于1，表明植株地上部氮濃度過高，氮肥施用過量；在N210 處理下NNI 值在1 上下波動，表明在當(dāng)前土壤肥力水平下氮肥施入量為210 kg·hm-2較為適宜。

圖5 辣椒氮素營養(yǎng)指數(shù)動態(tài)變化

此外，采用線性加平臺模型模擬NNI 與相對產(chǎn)量（RY）之間的關(guān)系，回歸方程決定系數(shù)為0.80，達(dá)到極顯著水平（圖6）。當(dāng)NNI≥0.95 時，RY獲得最大值（0.97），當(dāng)NNI＜0.95 時，RY 值隨NNI 的降低而減小。

圖6 辣椒氮營養(yǎng)指數(shù)與相對產(chǎn)量的關(guān)系

3 討論

臨界氮濃度稀釋曲線可快速有效診斷及評價植株氮素營養(yǎng)狀況（呂茹潔 等，2018；石小虎和蔡煥杰，2018；劉秋霞 等，2019；付江鵬 等，2020）。從數(shù)學(xué)的角度來講，參數(shù)a代表植株地上部干質(zhì)量為1 t·hm-2時的植株氮濃度，表示作物生育初期內(nèi)在的需氮特性；參數(shù)b表示植株氮濃度隨植株地上部干質(zhì)量的變化情況，其值取決于植株氮素吸收量與干質(zhì)量之間的比例關(guān)系（Gastal &Lemaire，2002）。本試驗建立的辣椒臨界氮濃度稀釋曲線模型（Nc＝4.39PDM-0.28）與向友珍等（2016）針對西北日光溫室甜椒在最佳灌水量下建立的臨界氮濃度稀釋曲線模型（Nc＝4.71PDM-0.320）參數(shù)相比明顯偏低，造成參數(shù)a、b值偏低的原因一方面可能是西北日光溫室甜椒供試土壤為重壤土，其氮素礦化率比貴州地區(qū)黃壤土高，而且在溫室中溫濕度可控制的前提下，適宜的土壤溫度和濕度也使微生物的活性更強(qiáng)，提高了土壤的供氮能力（Justes et al.，1994）；另一方面可能是由于不同辣椒品種對氮素敏感性程度不同，在氮素吸收效率和利用效率等方面存在差異，導(dǎo)致模型參數(shù)不同。

辣椒屬于C3作物中的茄科類蔬菜，與同科類的番茄生長習(xí)性具有較大的相似性，結(jié)合王新等（2013）針對滴灌番茄所建立的模型（Nc＝4.352PDM-0.274）參數(shù)a值和b值與本試驗所構(gòu)建的模型參數(shù)極為接近，表明兩個生態(tài)區(qū)域氣候、土壤等環(huán)境條件對臨界氮濃度稀釋曲線斜率影響較小，即其存在較好的穩(wěn)定性，提高了模型的精度和普適性。此外，本試驗利用2020 年的數(shù)據(jù)對辣椒臨界氮濃度稀釋曲線模型進(jìn)行了驗證，RMSE 為0.26。從驗證結(jié)果可知，基于植株地上部干質(zhì)量構(gòu)建的辣椒臨界氮濃度稀釋曲線模型模擬效果較好。與向友珍等（2016）針對西北日光溫室甜椒建立的臨界氮濃度稀釋曲線模型模擬效果相比，本試驗建立的模型RMSE 相對較小，表明模型模擬的效果相對較好，這可能與施肥方式、不同生長環(huán)境下作物的不同氮濃度變化特性相關(guān)。此外，辣椒臨界氮濃度稀釋曲線模型n-RMSE 值為10.08%，表明模型的穩(wěn)定性較好，可以作為辣椒氮素營養(yǎng)狀況判斷的工具之一。

NNI 是衡量作物氮素狀態(tài)的理想指標(biāo)（向友珍等，2016；Hoogmoed &Sadras，2016；劉秋霞 等，2019）。本試驗中，2 個辣椒品種不同施氮處理下NNI 值在0.55～1.30 之間變化。在現(xiàn)蕾期（定植后30 d），不同氮素處理下的NNI 值均最小，而開花期（定植后50 d），除2020 年種植的蘇椒5 號外，NNI 值均隨著氮肥用量的增加呈現(xiàn)不同程度的上升趨勢。這一方面可能是現(xiàn)蕾期是辣椒生長發(fā)育最旺盛的階段，此時期辣椒需要養(yǎng)分的絕對數(shù)量和相對數(shù)量都最大且吸收速度也最快；另一方面可能是由于該時期辣椒開始由營養(yǎng)生長向生殖生長轉(zhuǎn)變，其對氮肥的需求不如現(xiàn)蕾期迫切，這表明施肥量和施肥時期均能影響作物的營養(yǎng)狀況，且NNI 能夠很好地對其作出反應(yīng)，與前人研究結(jié)果一致（Lemaire et al.，2008）。

關(guān)于辣椒最佳施氮量前人做了大量研究。唐恒朋（2016）研究表明，從辣椒優(yōu)質(zhì)高產(chǎn)的角度考慮，推薦氮素施用量以207 kg·hm-2較優(yōu)；葉潔（2017）研究表明，從氮肥減量的角度考慮，施氮量為172.2 kg·hm-2可作為辣椒栽培中控釋氮肥用量的參考依據(jù)。根據(jù)辣椒氮素營養(yǎng)指數(shù)的動態(tài)變化，本試驗認(rèn)為較適宜的辣椒施氮量為210 kg·hm-2，這與唐恒朋（2016）的研究結(jié)果基本一致。當(dāng)辣椒NNI ＜0.95 時，NNI 與RY 呈極顯著正相關(guān)，可以解釋不同施氮水平下辣椒相對產(chǎn)量的變化。Ziadi 等（2008）也驗證了玉米NNI ＜0.93 時，NNI 與RY 呈顯著正相關(guān)。因此，利用NNI 進(jìn)行作物氮營養(yǎng)狀況和產(chǎn)量的估計具有較大的潛力（Ata-Ul-Karim et al.，2016，2017）。

4 結(jié)論

隨生育進(jìn)程的推進(jìn)，辣椒植株地上部氮濃度逐漸降低，干質(zhì)量呈升高趨勢。根據(jù)植株地上部干質(zhì)量和氮濃度的關(guān)系，建立了辣椒臨界氮濃度稀釋曲線模型：Nc＝4.39PDM-0.28。模型的RMSE 為0.26，n-RMSE 為10.08%，表明模型有較好的穩(wěn)定性。辣椒氮素營養(yǎng)指數(shù)隨施氮量的增加而升高，施氮量為210 kg·hm-2時NNI 基本接近于1。根據(jù)辣椒相對產(chǎn)量與氮素營養(yǎng)指數(shù)的關(guān)系及NNI 的綜合表現(xiàn)，推薦貴州地區(qū)辣椒較適宜的施氮量為210 kg·hm-2。