鐘培閣，孫洪仁，張吉萍，呂玉才，王應(yīng)海

（1.中國農(nóng)業(yè)大學草業(yè)科學與技術(shù)學院，北京 100193；2.北京六凱農(nóng)業(yè)科技有限公司北京 100095；3.凱風新農(nóng)（北京）科技有限公司，北京 100095；4.東方智感（浙江）科技股份有限公司，杭州 311200）

0 引言

土壤養(yǎng)分豐缺指標法是測土施肥的經(jīng)典方法。許多學者對我國諸多區(qū)域主要作物的土壤養(yǎng)分豐缺指標開展了研究，但針對甘蔗的相關(guān)研究較少。王獻華[1]采用盆栽試驗法對廣西甘蔗土壤有效磷豐缺指標進行了初步探討；唐栓虎等[2]研究了華南坡崗地甘蔗土壤有效磷和速效鉀豐缺指標；黃振瑞[3]建立了粵西甘蔗土壤堿解氮、有效磷和速效鉀土壤養(yǎng)分豐缺指標推薦施肥體系。絕大部分地區(qū)甘蔗測土施肥尚無據(jù)可依。孫洪仁研究團隊[4-11]采用零散實驗數(shù)據(jù)整合法，成功創(chuàng)建了中國北方紫花苜蓿、中國農(nóng)牧交錯帶燕麥和中國北方甜菜土壤養(yǎng)分豐缺指標推薦施肥系統(tǒng)。本研究擬采用零散實驗數(shù)據(jù)整合法進行中國南方甘蔗土壤全氮豐缺指標與適宜施氮量研究，為我國南方甘蔗測土施氮提供科學依據(jù)。

1 材料與方法

1.1 中國南方甘蔗土壤全氮養(yǎng)分豐缺指標

利用網(wǎng)絡(luò)數(shù)據(jù)庫資源，收集中國南方開展的甘蔗施肥試驗文獻。選擇含有施肥前試驗地土壤全氮含量、缺氮處理產(chǎn)量及全肥處理產(chǎn)量3項指標的文獻，并提取以上數(shù)據(jù)。利用式（1）計算甘蔗缺氮處理相對產(chǎn)量。

式（1）中R為“缺氮處理相對產(chǎn)量”，Y-N為“缺氮處理產(chǎn)量”，Y為“全肥處理產(chǎn)量”。

理論上，施肥前土壤全氮含量越高，獲得的缺氮處理相對產(chǎn)量越高；反之，施肥前土壤全氮含量越低，獲得的缺氮處理相對產(chǎn)量越低。利用Excel軟件，選擇適當模型建立甘蔗缺氮處理相對產(chǎn)量與土壤全氮含量回歸方程。為了適當提高推薦施肥精準度，參考“測土施肥土壤養(yǎng)分豐缺分級改良方案”[12]，確定第1～11級缺氮處理相對產(chǎn)量下限依次為100%、90%、80%……20%、10%和0%。將各級別缺氮處理相對產(chǎn)量起始點數(shù)值0%、10%、20%……80%、90%和100%，作為因變量值（y）分別代入本研究建立的甘蔗缺氮處理相對產(chǎn)量與土壤全氮含量回歸方程，反向計算得到的土壤全氮含量（x）即為各豐缺級別的起點指標。超出缺氮處理相對產(chǎn)量和土壤全氮含量試驗范圍的外推數(shù)據(jù)，低端和高端分別保留1個。

1.2 中國南方甘蔗推薦施肥量

采用“養(yǎng)分平衡—地力差減法”確定適宜施肥量的新應(yīng)用公式（2）[13]和“土壤養(yǎng)分豐缺指標法之不同豐缺級別土壤適宜施肥量檢索表”[14]，確定各豐缺級別甘蔗土壤的適宜施氮量。

式（2）中F為“適宜施氮量”，A為“目標產(chǎn)量甘蔗氮素移出量”，R為“缺氮處理相對產(chǎn)量”，E為“氮肥當季利用率”。

單位經(jīng)濟產(chǎn)量甘蔗氮素（N）移出量確定為1.9 kg/t[15]。根據(jù)生產(chǎn)現(xiàn)狀，確定甘蔗目標產(chǎn)量10 個，分別為45、60、75、90、105、120、135、150、165 和180 t/hm2。選取各豐缺級別的缺氮處理相對產(chǎn)量下限作為該級別之甘蔗缺氮處理相對產(chǎn)量，第1 級為100%，第2 級為90%，第3級為80%，……，第11級為0%。設(shè)置3個甘蔗氮肥當季利用率，分別為30%、40%和50%。

2 結(jié)果與分析

2.1 中國南方甘蔗施氮試驗文獻及與土壤全氮豐缺指標研究相關(guān)信息

搜集到在中國南方開展的同時含有“土壤全氮含量”、“缺氮處理產(chǎn)量”和“全肥處理產(chǎn)量”的甘蔗施氮試驗文獻28篇（2000—2009年8篇，2010—2019年20篇）[16-43]，提取出配套數(shù)據(jù)60組，進而獲得“土壤全氮含量”和根據(jù)公式（1）計算“甘蔗缺氮處理相對產(chǎn)量”配套數(shù)據(jù)60對（表1）。

表1 甘蔗施氮試驗文獻中與土壤全氮豐缺指標研究相關(guān)信息Table 1 Relevant information about the research of soil total N abundance-deficiency index in literature of sugarcane fertilization experiment

續(xù)表1 Continued Table 1

2.2 中國南方甘蔗土壤全氮含量與缺氮處理相對產(chǎn)量回歸方程

剔除5 對明顯不合理數(shù)據(jù)，利用55對配套數(shù)據(jù)建立“土壤全氮含量”（x）與“缺氮處理相對產(chǎn)量”（y）回歸方程，分別采用直線和自然對數(shù)兩種模型進行模擬，結(jié)果如圖1 和表2 所示。土壤全氮含量試驗數(shù)據(jù)范圍為0.15～1.61 g/kg，甘蔗缺氮處理相對產(chǎn)量試驗數(shù)據(jù)范圍為51.76%～100.78%。

圖1 中國南方甘蔗土壤全氮含量與缺氮處理相對產(chǎn)量回歸曲線Fig.1 The regression curve between soil total N and relative yield of no N fertilizer treatment for sugarcane in south China

表2 中國南方甘蔗土壤豐缺級別第1～7 級起訖點土壤全氮含量(g/kg)Table 2 The starting and ending point values of soil total N contents for the first to seventh grade of the abundance-deficiency level in sugarcane in the south of China(g/kg)

2.3 中國南方甘蔗土壤全氮豐缺指標

利用缺氮處理相對產(chǎn)量和土壤全氮含量回歸方程，計算中國南方甘蔗土壤豐缺級別第1～7 級起訖點土壤全氮含量，如表2 所示。100%～50%的缺氮處理相對產(chǎn)量值對應(yīng)計算得到的是土壤全氮豐缺級別第1～7級的起點指標。采用直線和自然對數(shù)兩種類型回歸方程的計算結(jié)果差別明顯。利用直線型回歸方程計算的結(jié)果中，缺氮處理相對產(chǎn)量50%對應(yīng)的土壤全氮含量為負值。利用自然對數(shù)型回歸方程計算的結(jié)果中，缺氮處理相對產(chǎn)量90%～100%對應(yīng)的土壤全氮含量為外推數(shù)據(jù)。利用直線型回歸方程計算得到的缺氮處理相對產(chǎn)量50%～60%對應(yīng)的土壤全氮含量皆明顯低于自然對數(shù)型。利用直線型回歸方程計算得到的缺氮處理相對產(chǎn)量80%對應(yīng)的土壤全氮含量略高于自然對數(shù)型，但二者頗為接近。利用直線型回歸方程計算得到的缺氮處理相對產(chǎn)量95%和100%對應(yīng)的土壤全氮含量皆明顯低于自然對數(shù)型。舍棄負值和超出試驗范圍的外推數(shù)據(jù)后，再根據(jù)公式（2）計算得到中國南方甘蔗土壤全氮1～7 級豐缺指標的適宜施氮量，如表3所示。

表3 中國南方甘蔗土壤全氮豐缺指標及不同目標產(chǎn)量對應(yīng)的適宜施氮量（N，kg/hm2）Table 3 Abundance-deficiency index of soil total N content and appropriate fertilizer application rates corresponding to target cane yield for sugarcane in the south of China

2.4 中國南方不同豐缺級別甘蔗土壤適宜施氮量

表3列出了中國南方甘蔗不同目標產(chǎn)量和不同肥料利用率下不同豐缺級別土壤的適宜施氮量。適宜施氮量與土壤全氮豐缺級別及氮肥當季利用率呈直線負相關(guān)，與甘蔗目標產(chǎn)量則呈現(xiàn)直線正相關(guān)。當?shù)世寐?0%、甘蔗目標產(chǎn)量45～180 t/hm2時，土壤全氮豐缺級別第1～7級的適宜施氮量范圍依次為0～0、21～86、43～171、64～257、86～342、107～428 和128～513 kg/hm2；當?shù)世寐侍岣叩?0%，適宜施氮量范圍依次為0～0、17～68、34～137、51～205、68～274、86～342、103～410 kg/hm2，比40%氮肥利用率分別節(jié)省施氮量0～0、4～18、9～34、13～52、18～68、21～86、25～105 kg/hm2。假如南方某地測得土壤全氮含量為1.64 g/kg，當?shù)世寐?0%時，要達到甘蔗目標產(chǎn)量90 t/hm2，那么由此推算適宜施氮量應(yīng)為86 kg/hm2。

3 討論

3.1 直線和自然對數(shù)模型的比較

自然對數(shù)模型在作物缺素處理相對產(chǎn)量與土壤養(yǎng)分含量回歸方程建立中應(yīng)用最為廣泛，而使用直線模型的研究則相對較少[41-46]。本研究結(jié)果表明，直線模型存在豐缺級別低端指標數(shù)值為負的不合理問題，而自然對數(shù)模型存在豐缺級別高端指標外推且偏高的問題。采用自然對數(shù)模型計算土壤養(yǎng)分豐缺級別之中、低端指標較為理想，而對于高端指標，直線模型計算結(jié)果可能更為符合實際情況。

3.2 自然對數(shù)模型最高豐缺級別土壤缺素處理相對產(chǎn)量下限的選擇

從理論上講，土壤缺素處理相對產(chǎn)量下限應(yīng)為100%，此時無需施肥也可達到目標產(chǎn)量。但在本研究自然對數(shù)模型得到的結(jié)果中，對應(yīng)于缺素處理相對產(chǎn)量100%之最高豐缺級別土壤養(yǎng)分含量指標過高，與正常認知有些不相符，并且在針對其他作物的研究中也出現(xiàn)了類似的現(xiàn)象[44-52]。自然對數(shù)模型缺素處理相對產(chǎn)量95%對應(yīng)的土壤養(yǎng)分含量與直線模型缺素處理相對產(chǎn)量100%對應(yīng)的土壤養(yǎng)分含量較為接近。因此，可選用95%來代替100%作為自然對數(shù)模型最高豐缺級別土壤缺素處理相對產(chǎn)量下限。

3.3 中國南方甘蔗土壤全氮豐缺指標

孫洪仁等[7，44-46]的研究表明，中國燕麥、玉米、小麥和馬鈴薯缺氮處理相對產(chǎn)量70%和80%對應(yīng)的土壤全氮含量范圍分別為0.6～1.4和0.9～2.0 g/kg。本中國南方甘蔗相應(yīng)研究結(jié)果皆介于其范圍之間。

3.4 中國南方甘蔗適宜施氮量

黃振瑞[3]針對粵西甘蔗的氮素推薦施用量為0～390 kg/hm2，本研究之甘蔗適宜施氮量范圍較其為大。本研究涉及文獻[16-43]針對甘蔗目標產(chǎn)量為40～190 t/hm2的氮素推薦施用量為69～656 kg/hm2。本研究之氮素適宜施用量與其接近。

4 結(jié)論

本研究建立了中國南方甘蔗土壤全氮豐缺指標，確定了不同目標產(chǎn)量和不同氮肥當季利用率情形下不同豐缺級別的甘蔗適宜施氮量，為我國南方甘蔗測土施氮提供了科學依據(jù)。