基于臨界氮濃度的滴灌棉花氮素營養診斷模型研究

馬露露 呂 新 張 澤 馬革新 海興巖

(1.石河子大學生命科學學院， 石河子 832003； 2.石河子大學農學院, 石河子 832003；3.新疆生產建設兵團綠洲生態農業重點實驗室， 石河子 832003)

0 引言

氮素是干旱區作物增產增效的重要限制因素，國內外關于滴灌作物養分運移與分配、作物生長與產量、施肥制度、氮肥優化管理等方面開展了大量研究[1]。但在作物生長過程中氮素需求和供應之間的同步性研究較少[2]，導致氮素利用效率偏低，同時造成溫室氣體排放加劇和大氣、土壤及水體的污染。因此，針對作物不同生育期氮素營養實際狀況優化供氮量，對于提高氮素利用效率和環境保護具有重要意義。

研究發現,不同作物、生態氣候區域及水肥管理方式之間的模型存在一定的差異，新疆屬于典型干旱區，滴灌水肥一體化模式對氮素吸收影響較大，因此需要進行臨界氮濃度稀釋模型參數的本地化校正。

本研究以我國重要經濟作物棉花為研究對象，在滴灌條件下開展不同施氮水平田間試驗，以確定滴灌棉花生長發育與氮素吸收利用規律為基礎，建立滴灌棉花臨界氮稀釋濃度模型，最終確定氮濃度營養指數，實現氮素準確、快速診斷，為進一步提高滴灌棉花的氮素利用效率、氮肥合理施用水平提供理論依據。

1 材料與方法

1.1 試驗材料與設計

試驗于2015—2016年在新疆瑪納斯縣六戶地鎮高產創田地(東經86°07′36″，北緯44°39′15″，圖1)進行，日照時間達2 800～3 000 h，年平均氣溫5℃，大于等于10℃有效積溫3 500～4 100℃，無霜期180 d。土壤質地為壤土，0～20 cm土層有機質質量比19.13 g/kg，堿解氮質量比50.8 mg/kg，速效磷質量比19.8 mg/kg，速效鉀質量比160.1 mg/kg，田間持水量29.6 g/m3，耕層土壤容重為1.43 g/cm3。

圖1 試驗地位置圖Fig.1 Experimental field location map

滴灌棉花的供試品種為新陸早45，試驗設6個氮素水平N0、N1、N2、N3、N4(0、120、240、360、480 kg/hm2)和N5(常規處理，施氮量為300 kg/hm2，為當地農戶施純氮量，作為模型的驗證)，氮肥施用均為30%作為基肥，剩余70%在滴灌棉花生育期隨水施入。磷肥(P2O5)150 kg/hm2和鉀肥(K2O)150 kg/hm2播前一次性施入。試驗設3次重復，隨機區組排列(圖2)。采用膜下滴灌的方式，氮肥(尿素，含氮量46%)施入時間與當地農戶同步。試驗小區面積為25 m2(2.5 m×10 m)，種植模式為“1膜3管6行”，膜寬2.05 m，株距10 cm，行距配置為10 cm+66 cm，種植密度為22.5萬株/hm2。試驗地為棉花連作，東西走向，四周設置保護行。2015年4月11日播種，并滴施出苗水，其他田間管理措施均按高產栽培要求進行。

圖2 試驗小區分布圖Fig.2 Field experiment district distribution map

1.2 測定內容與方法

2015和2016年試驗均在滴灌棉花出苗后開始計算采樣日期，直至棉花始絮期。第2次滴灌施肥5 d以后開始采樣，其中施氮和滴灌日期為6月17日、6月23日、7月1日、7月10日、7月22日、8月5日、8月20日，施氮比例分別為10%、10%、20%、20%、20%、10%、10%，采樣時間分別為出苗后第60、68、77、90、105、125天，2年試驗施氮比例、施氮日期、采樣時間相同。每個試驗小區選取具有代表性的棉花3～6株，按莖、葉和雷鈴3類不同器官分樣，稱量其各部位鮮質量，在105℃下殺青30 min，80℃干燥至恒質量，然后稱量其干物質量。粉碎后用凱氏定氮法測定各分器官的全氮含量，據此計算滴灌棉花單株氮素累積量。成熟后分小區收獲、測產。

1.3 滴灌棉花臨界氮濃度曲線模型的建立方法

1.3.1滴灌棉花臨界氮濃度曲線模型

根據GREENWOOD等[8]對臨界氮濃度的定義，對于滴灌棉花，地上干物質氮濃度值在其生長發育不受氮素制約的情況下達到最低值時為臨界氮濃度。

圖3 不同施氮水平下滴灌棉花地上部各器官氮素吸收分配規律Fig.3 Nitrogen absorption distribution rule of aboveground organs of drip-irrigation cotton under different nitrogen levels

臨界氮濃度稀釋曲線的模型為

(1)

式中Nc——臨界氮濃度，g/(100 g)

a——滴灌棉花地上部生物量為1 t/hm2時棉株的臨界氮濃度

Wmax——滴灌棉花地上部最大生物量

b——控制臨界氮濃度稀釋曲線斜率的參數

用均方根誤差(Root mean square error，RMSE)和標準化均方根誤差(n-RMSE)對所得模型進行驗證，然后建立模擬值和實測值1∶1直方圖顯示模型的擬合度和可靠程度。其中RMSE值越小，表示模擬值與實測值的偏差越小；n-RMSE小于10%，模型模擬性能極好；n-RMSE為[10%，20%]，模型模擬性能較好；n-RMSE為[20%，30%]，模型模擬性能一般；n-RMSE大于30%，模擬性能較差[23]。

1.3.2滴灌棉花的異速生長模型

滴灌棉花在臨界氮濃度下達到累積的最大地上部生物量時所需的氮吸收量Nuptc(kg/hm2)為

Nuptc=10NcWmax

(2)

將式(1)代入式(2)，得到滴灌棉花臨界氮吸收模型，即滴灌棉花氮累積量與地上部生物量之間的異速生長模型

(3)

1.3.3滴灌棉花臨界氮營養指數

氮營養指數(Nitrogen nutrition index，NNI)是指作物地上部植株實際的氮濃度與臨界氮濃度的比值，用來判定作物體內氮營養狀況[24-25]，是LEMAIRE等[26]基于臨界氮濃度稀釋模型提出了氮營養指數(NNI)的概念，計算式為

(4)

式中Ni——地上部實測氮濃度，g/(100 g)

若NNI=1，表明作物氮素營養水平處于最佳狀態，高于1為氮營養過剩，低于1則氮營養不足。

1.4 數據處理與分析

采用2015年的試驗數據進行回歸擬合，建立滴灌棉花臨界氮濃度稀釋曲線模型，利用2016年的常規試驗數據對模型進行驗證。采用Excel 2013進行數據的基礎整理，顯著性分析在SPSS 19.0 軟件中的LSD和Duncan進行，圖形通過Origin 8.0進行繪制。

2 結果與分析

2.1 施氮水平對滴灌棉花氮素吸收分配規律的影響

由圖3可以看出，滴灌棉花在不同施氮水平下，地上部各器官的氮素含量存在顯著性差異，其中葉片含量最高，蕾鈴次之，莖的含氮量最低。隨著滴灌棉花的生長發育，葉和莖的含氮量呈下降趨勢，其中葉片的含氮量從47.05 g/kg降到23.33 g/kg，降幅為50.41%，莖的含氮量從16.84 g/kg降到7.00 g/kg，降幅為58.43%。由于土壤基礎肥力的影響，N0、N1、N2和N3處理對葉和莖的含氮量沒有顯著影響，在N4水平下有顯著差異；葉和莖含氮量的降低與氮素向蕾鈴轉移有關，在棉花生長期內，蕾鈴的含氮量先下降后上升，這是由于后期幼鈴增多且成鈴數也在增加，蕾鈴的生物量大幅增加，單鈴的含氮量呈逐漸下降的趨勢；出苗后105 d以后即盛鈴期至始絮期，隨著棉花整個生育期施氮的完成，葉片和莖的含氮量顯著下降，同時由于棉花種子及纖維的發育，棉鈴的含氮量出現回升。

2.2 施氮水平對滴灌棉花地上生物量和氮濃度動態累積的影響

以W0、W1、W2、W3、W4 分別表示N0、N1、N2、N3、N4施氮水平下滴灌棉株的地上部生物量(kg/hm2)，表1試驗結果表明，通過顯著性分析，在不同施氮水平下的生物量由大到小依次為W3和W4和W2、W1、W0，滴灌棉花的地上部生物量隨著生物進程而增加，N0、N1、N2、N3、N4水平下地上部生物量從第60天到第125天增幅分別為243.49%、423.48%、519.93%、573.53%、561.79%，其中在N3水平下增幅最大，與N0、N1有顯著差異，與N2、N4無顯著差異；說明在其他條件平等的情況下，過多的施入氮肥并不能使滴灌棉花的生物量顯著增加。

表1 施氮水平對滴灌棉花生物量動態累積變化的影響Tab.1 Effect of nitrogen levels on dynamic change of nitrogen accumulation in aboveground biomass of drip-irrigation cotton

注：同一列的不同字母代表差異顯著(P<0.05)，下同。

如圖4所示，其氮濃度隨生物進程的進行而降低，出苗后68 d的氮濃度下降得最快。通過顯著性分析(表2)，其氮濃度在不同施氮水平下從大到小依次為N3和N4和N2、N1、N0；在N2、N3和N4條件下，棉株地上部氮濃度差異不顯著，說明在棉田土壤基礎肥力的影響下，棉花的生長發育不受N2、N3和N4施氮水平的顯著影響。

圖4 不同施氮水平下滴灌棉花地上部氮濃度動態累積變化Fig.4 Change of aboveground nitrogen concentration accumulation dynamic of drip-irrigation cotton under different nitrogen levels

結合JUSTES等[27]和薛曉萍等[21]提出的臨界氮濃度稀釋曲線的計算方法，本研究中滴灌棉花臨界氮濃度模型的建立步驟為：①分析不同施氮水平下，滴灌棉花地上部干物質與其氮濃度，通過對施氮水平進行方差分析，氮素是否制約棉花的正常生長發育。②對于不受氮素限制的施氮水平，采用其地上部干物質的平均值作為該采樣日期的最大干物質。③對于不能滿足棉花生長發育氮需求量的施氮水平，對其地上干物質與氮濃度值進行線性或者曲線擬合。④采樣日的臨界氮濃度為由步驟③的線性或曲線與橫坐標上最大干物質的垂線交點的縱坐標。

表2 不同施氮水平下滴灌棉花地上部氮濃度顯著性檢驗Tab.2 Significance test of aboveground nitrogen concentration of drip-irrigation cotton under different nitrogen levels

2.3 滴灌棉花臨界氮濃度稀釋曲線模型建立

由于新疆地區特殊的地理氣候條件，滴灌棉田的種植密度比內地大，因此單位面積棉花的生物量較高。在出苗后第60～125天的采樣時間里，滴灌棉花的生物量范圍在3～22 t/hm2，全部大于1 t/hm2；氮濃度在17～35 g/kg范圍內。通過分析得出，滴灌棉花在N0和N1條件下為限氮水平，由臨界氮濃度的確定條件和式(1)得出滴灌棉花臨界氮濃度稀釋曲線(圖5)及模型，表3為模型參數。分析實測數據得出在不同施氮水平下，同樣的地上生物量值其所對應氮濃度值不同，取每個采樣日氮濃度的最大、最小實測值模擬得到其最高(Nmax，%)、最低(Nmin，%)氮濃度稀釋模型，結果同樣也符合模型(1)，這與薛曉萍等[21,28]的結果相符，模型參數見表3。由式(2)和(3)得到異速生長模型，2個模型的參數見表3。其中10a值與薛曉萍等[21]和王子勝等[22]的結果存在差異，這與新疆土壤類型、農田生態和水肥供應方式有關。

圖5 滴灌棉花氮濃度稀釋曲線Fig.5 Nitrogen concentration dilution curves for drip-irrigation cotton

類型稀釋曲線abR2類型異速生長曲線10a1-bR2Nmax4.450.220.973Nuptcmax44.510.780.954Nmin3.420.220.882Nuptcmin34.210.780.897Nc3.910.240.906Nuptc39.140.760.933

2.4 滴灌棉花臨界氮濃度稀釋曲線模型的驗證

采用獨立試驗數據對模型進行驗證，研究選取2016年常規試驗即施純氮量為300 kg/hm2的大田試驗的全生育期的6次采樣的單株滴灌棉花地上部生物量來驗證模型。將獨立試驗的滴灌棉花地上部生物量(其中6月28日取樣3株、7月6日取樣3株、7月15日取樣3株、7月28日取樣3株、8月12日取樣3株、9月3日取樣4株，共計n=19，)代入所得滴灌棉花臨界氮濃度稀釋曲線模型中得到模擬值，然后做出模擬值和實測值的1∶1直方圖(圖6)，通過RMSE和n-RMSE來評價模型的精準度。其中RMSE為3.71 g/kg，n-RMSE為13.60%，n-RMSE為[10%，20%]，說明模型的模擬性能較好，可以進一步用于滴灌棉花的營養診斷。

圖6 滴灌棉花臨界氮濃稀釋曲線模型的驗證Fig.6 Validation of critical nitrogen concentration dilution curve model for drip-irrigation cotton

2.5 基于滴灌棉花臨界氮濃度稀釋曲線的氮營養診斷

本研究采用氮素營養指數作為判斷滴灌棉花臨界氮濃度稀釋模型的依據。由式(4)，通過氮營養指數得出滴灌棉花的營養狀況。如圖7所示，隨著施氮水平的增加，氮營養指數不斷提高。不同氮素處理的滴灌棉花的NNI值分別為：N0(0.74～0.96)、N1(0.84～0.97)和N2(0.92～0.99)水平在出苗后第60～125天內NNI值均小于1，說明N0、N1和N2水平下，滴灌棉花的生長受到了氮素的限制，N3(1.005～1.03)和N4(1.06～1.12)水平下NNI值均大于1，說明在N3、N4水平下滴灌棉花屬于過量氮素的施入，因此其適宜施氮量在N2和N3之間。

圖7 滴灌棉花氮營養指數動態變化Fig.7 Dynamic changes of nitrogen nutrient index of drip-irrigation cotton

3 討論

3.1 滴灌棉花臨界氮濃度稀釋曲線的模型特征

本研究建立了新疆干旱棉區滴灌棉花除苗期之外全部生育期臨界氮濃度稀釋曲線模型，其模型在形式上與薛曉萍等[21]、王子勝等[22]針對長江流域、黃河流域和東北特早熟棉區建立的棉花臨界氮濃度稀釋模型一致，其采樣時間(苗期后期)與薛曉萍等[21](初花期)的采樣時間大致一樣，這說明臨界氮濃度稀釋模型的建立與棉花的生育期無明顯的關系；所得模型的b值(0.24)與薛曉萍等[21]淹灌的模型b值(南京，0.131；安陽，0.131)差異較大，與王子勝等[22]所建立的模型b值(0.252)極為接近，但a值(3.91)與薛曉萍等[21](黃淮棉區，2.858；長江中下游棉區，3.387)和王子勝等[22](東北早熟棉區，4.377)的研究結果均有差異，說明新疆滴灌棉的氮吸收能力比安陽和南京的棉花氮吸收能力強，但是低于東北早熟棉區的氮吸收能力。其原因：①新疆干旱棉區采用滴灌方式，肥隨水入，每個生育期都會施入其所需的氮肥，這與薛曉萍等[21]分30%、70%施入和王子勝等[22]將氮肥一次性施入不同，說明施肥方式對棉花的氮素吸收能力有一定的影響。②新疆由于土壤類型和氣候原因，采用“密、早、矮、膜”的栽培方式，其密度(19.6萬株/hm2)為王子勝等[22]研究的最大密度(12萬株/hm2)的1.6倍，使得其在同一生育期的單位面積內生物量(表1)是東北早熟棉區2倍之多，但是在所得單位生物量下，所需的臨界氮濃度(表2)比之較小，說明新疆特有的高密度種植模式符合新疆滴灌棉花提高經濟效益的需求。

3.2 基于臨界氮濃度稀釋曲線和NNI滴灌棉花氮營養診斷

基于臨界氮濃度稀釋模型得出的氮營養指數能診斷滴灌棉花不同生長階段的氮營養狀況，還可以量化滴灌棉花氮素脅迫的程度。本研究通過滴灌棉花臨界氮濃度稀釋模型的建立得出，新疆滴灌棉區施氮量應不低于240 kg/hm2，不高于360 kg/hm2。

通過對新疆滴灌棉區不同施氮水平的臨界氮濃度稀釋模型和氮營養指數的研究，結合前人的研究，結果表明不同棉區的氮肥施入量、施肥方式和種植模式都會影響棉花的氮素吸收水平，進而對棉花的生物量、氮濃度和產量產生影響。新疆干旱棉區其特殊的滴灌模式和根據生育時期隨水施入的施肥模式能滿足干旱區棉花對于氮素養分的需要，更有效地提高氮肥利用率，降低過量施肥給土壤帶來的壓力，進一步改善土壤生態環境。

4 結論

(1)以滴灌棉花、設置5個施氮水平進行試驗，探討了施氮量對滴灌棉花氮素吸收分配規律的影響。滴灌棉花各器官含氮量的關系由大到小為葉片、蕾鈴、莖，葉片和莖的含氮量隨生物進程而降低，低于360 kg/hm2的施氮量對于葉片和莖的含氮量沒有顯著影響。

(2)在一定施氮水平下，滴灌棉花的地上部生物量隨施氮量的增加而增加，氮濃度隨植株生長發育的進行而降低；建立和驗證了滴灌棉花臨界氮濃度稀釋曲線模型，結果表明，滴灌棉花各生育期的最大生物量與氮濃度之間符合冪指數關系。

(3)基于臨界氮濃度稀釋曲線模型，通過氮營養指數對新疆干旱區滴灌棉花進行了氮營養診斷，結果顯示，在本試驗條件下，240～360 kg/hm2為滴灌棉花的適宜施氮量。

1 山立, 鄒宇鋒. 我國旱區農業的地位和發展潛力及政策建議[J]. 農業現代化研究, 2013, 34(4): 425-430.

SHAN Li, ZOU Yufeng. Status and potential for development of China’s dryland agriculture and suggestions on its policy making[J]. Research of Agricultural Modernization, 2013, 34(4): 425-430.(in Chinese)

2 BRYE K R, NORMAN J M, GOWER S T, et al. Methodological limitations and N-budget differences among a restored tallgrass prairie and maize agroecosystems[J]. Agriculture, Ecosystems & Environment, 2003, 97(1): 181-198.

3 REN T, CHRISTIE P, WANG J, et al. Root zone soil nitrogen management to maintain high tomato yields and minimum nitrogen losses to the environment[J]. Scientia Horticulturae, 2010, 125(1): 25-33.

4 ZHENG H, LIU Y, QIN Y, et al. Establishing dynamic thresholds for potato nitrogen status diagnosis with the SPAD chlorophyll meter[J]. Journal of Integrative Agriculture, 2015, 14(1): 190-195.

5 HANSEN P M, SCHJOERRING J K. Reflectance measurement of canopy biomass and nitrogen status in wheat crops using normalized difference vegetation indices and partial least squares regression[J]. Remote Sensing of Environment, 2003, 86(4): 542-553.

6 ULRICH A. Physiological bases for assessing the nutritional requirements of plants[J]. Annual Review of Plant Physiology, 1952, 3(1): 207-228.

7 ZIADI N, BRASSARD M, BéLANGER G, et al. Critical nitrogen curve and nitrogen nutrition index for corn in eastern Canada[J]. Agronomy Journal, 2008, 100(2): 271-276.

8 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.

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

10 DEBAEKE P, VAN OOSTEROM E J, JUSTES E, et al. A species-specific critical nitrogen dilution curve for sunflower (HelianthusannuusL.)[J]. Field Crops Research, 2012, 136: 76-84.

11 梁效貴, 張經廷, 周麗麗, 等. 華北地區夏玉米臨界氮稀釋曲線和氮營養指數研究[J]. 作物學報, 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)

12 BEN Zhao. Determining of a critical dilution curve for plant nitrogen concentration in winter barley[J]. Field Crops Research, 2014,160:64-72.

13 ATA-Ul-KARIM S T, YAO X, LIU X, et al. Determination of critical nitrogen dilution curve based on stem dry matter in rice[J]. PLOS ONE, 2014, 9(8): e104540.

14 陳玲玲, 毛培勝, 王明亞, 等. 臨界氮濃度稀釋模型在牧草種子生產中的應用前景[J]. 草地學報, 2015, 23(2): 231-236.

CHEN Lingling, MAO Peisheng, WANG Mingya,et al. Application prospect of critical nitrogen concentration dilution model in grass seed prodution[J]. Acta Agrestia Sinical, 2015, 23(2): 231-236.(in Chinese)

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

YANG Hui, CAO Hongxia, LIU Meiyu, et al. Simulation of critical nitrogen concentration and nitrogen nutrition index of tomato under different water and nitrogen conditions[J]. Plant Nutrition and Fertilizer Science ,2015, 21(5): 1234-1242.(in Chinese)

16 向友珍, 張富倉, 范軍亮, 等. 基于臨界氮濃度模型的日光溫室甜椒氮營養診斷[J]. 農業工程學報, 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 CSAE, 2016, 32(17):89-97.(in Chinese)

17 杜婭丹, 曹紅霞, 谷曉博, 等. 基質栽培番茄臨界氮濃度和氮營養指數研究[J]. 節水灌溉, 2016(9): 1-7.

DU Yadan, CAO Hongxia, GU Xiaobo, et al. Critical nitrogen dilution curve and nitrogen nutrition idex for soilless tomato[J]. Water Saving Irrigation,2016(9): 1-7.(in Chinese)

18 YUE S, MENG Q, ZHAO R, et al. Critical nitrogen dilution curve for optimizing nitrogen management of winter wheat production in the North China Plain[J]. Agronomy Journal, 2012, 104(2): 523-529.

19 ZIADI N, BéLANGER G, CLAESSENS A, et al. Determination of a critical nitrogen dilution curve for spring wheat[J]. Agronomy Journal, 2010, 102(1): 241-250.

20 HU D, SUN Z, LI T, et al. Nitrogen nutrition index and its relationship with N use efficiency, tuber yield, radiation use efficiency, and leaf parameters in potatoes[J]. Journal of Integrative Agriculture, 2014, 13(5): 1008-1016.

21 薛曉萍, 周治國, 張麗娟, 等. 棉花花后臨界氮濃度稀釋模型的建立及在施氮量調控中的應用[J].生態學報, 2006, 26(6): 1781-1791.

XUE Xiaoping, ZHOU Zhiguo, ZHANG Lijuan, et al. Development and application of critical nitrogen concentration dilution model for cotton after flowering[J]. Acta Ecologica Sinica,2006, 26(6): 1781-1791. (in Chinese)

22 王子勝, 金路路, 趙文青, 等. 東北特早熟棉區不同群體棉花氮臨界濃度稀釋模型的建立初探[J].棉花學報, 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)

23 CHEN X, ZHANG F, CUI Z, et al. Critical grain and stover nitrogen concentrations at harvest for summer maize production in China[J]. Agronomy Journal, 2010, 102(1): 289-295.

24 Ata-Ul-Karim S T, ZHU Y, YAO X, et al. Determination of critical nitrogen dilution curve based on leaf area index in rice[J]. Field Crops Research, 2014, 167: 76-85.

25 FIZGERALD G, RODRIGUEZ D, O’LEARY G. Measuring and predicting canopy nitrogen nutrition in wheat using a spectral index—the canopy chlorophyll content index (CCCI)[J]. Field Crops Research, 2010, 116(3): 318-324.

26 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.

27 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.

28 強生才, 張富倉, 向友珍, 等. 關中平原不同降雨年型夏玉米臨界氮稀釋曲線模擬及驗證[J]. 農業工程學報, 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 CSAE, 2015, 31(17): 168-175. (in Chinese)