單作套作大豆葉片氮素積累與光譜特征

諶俊旭 黃 山 范元芳 王 銳 劉沁林 楊文鈺 楊 峰

?

單作套作大豆葉片氮素積累與光譜特征

諶俊旭 黃 山 范元芳 王 銳 劉沁林 楊文鈺*楊 峰*

四川農(nóng)業(yè)大學(xué)農(nóng)學(xué)院 / 農(nóng)業(yè)部西南作物生理生態(tài)與耕作重點(diǎn)實(shí)驗(yàn)室 / 四川省作物帶狀復(fù)合種植工程技術(shù)研究中心, 四川成都 611130

種植模式和氮肥水平直接影響作物的生長(zhǎng)和氮素的吸收, 無(wú)損、即時(shí)監(jiān)測(cè)大豆葉片氮素水平對(duì)大豆生產(chǎn)中的氮肥精確管理十分重要。本研究設(shè)置4個(gè)氮肥水平, 分析單作套作下大豆在不同生育時(shí)期葉片氮素動(dòng)態(tài)和光譜特征, 明確對(duì)葉片氮素敏感的光譜特征參數(shù), 構(gòu)建單作套作大豆通用的葉片氮素積累量估測(cè)模型。結(jié)果表明, 隨大豆生育時(shí)期的推進(jìn), 單作套作種植模式下的大豆冠層葉片氮素積累量均呈現(xiàn)單峰變化趨勢(shì), 最大值出現(xiàn)在N3處理下的結(jié)莢期, 兩種模式兩年最大值平均分別為8.70 g m–2和8.38 g m–2; 不同生育時(shí)期和種植模式的大豆冠層原始反射光譜的變化規(guī)律與冠層葉片氮素變化規(guī)律均為先增加后降低, 原始反射光譜在700~1000 nm波段的反射率以結(jié)莢期為拐點(diǎn)先增大后減小, 最大反射率達(dá)到60%~70%左右; 通過(guò)對(duì)單作套作大豆冠層光譜一階導(dǎo)數(shù)變換, 紅邊幅值呈現(xiàn)先增加后降低的趨勢(shì), 同時(shí)紅邊位置隨葉片氮積累量的增加和減小出現(xiàn)“紅移”與“藍(lán)移”現(xiàn)象。經(jīng)波段自由組合和回歸分析表明, 以DSI (771、755)構(gòu)建的線性(= –1.249+3.209,2= 0.847)和乘冪(= –1.4701.676,2= 0.872)模型能較精確地估測(cè)不同生育時(shí)期大豆冠層葉片氮素狀況。

大豆; 氮素積累; 光譜反射率; 模型

高光譜遙感自20世紀(jì)70年代以來(lái), 就在農(nóng)作物監(jiān)測(cè)方面扮演著重要的角色。農(nóng)業(yè)遙感主要以作物和土壤為對(duì)象, 利用地物的光譜特性來(lái)監(jiān)測(cè)作物長(zhǎng)勢(shì)[1]、作物營(yíng)養(yǎng)[2]及作物病蟲害[3]等。氮素作為作物的重要營(yíng)養(yǎng)元素, 不僅對(duì)作物正常的生長(zhǎng)發(fā)育有著十分顯著的影響, 也對(duì)作物品質(zhì)起著十分重要的作用[4-7]。胡珍珠等[8]研究表明, 在核桃的不同生育時(shí)期, 與葉片氮素含量相關(guān)的特征光譜各有不同。同時(shí), Walburg等[9]研究也發(fā)現(xiàn), 缺氮時(shí)玉米冠層葉片光譜的反射率在紅光范圍內(nèi)增加, 但在近紅外波段降低。岳延濱等[10]研究發(fā)現(xiàn)在760~1350 nm范圍內(nèi), 用缺氮條件下冠層光譜反射率能夠比較準(zhǔn)確地估測(cè)辣椒的葉片全氮含量。田永超等[11]研究發(fā)現(xiàn), 綠光560 nm和紅光705 nm波段附近的反射率與近紅外波段組合而成的光譜比值指數(shù)能夠很好地監(jiān)測(cè)水稻葉層全氮含量。Wei等[12]利用優(yōu)化的指數(shù)構(gòu)建的估測(cè)模型在估測(cè)冬小麥冠層葉片氮含量時(shí)能夠部分消除品種差異、葉片水分含量、土壤性質(zhì)及生育時(shí)期的影響且估測(cè)精度高。同時(shí)氮素的含量對(duì)大豆品質(zhì)[13]、產(chǎn)量[14]都有明顯的影響。但是影響大豆光譜特性的因素較多, 且不同氮素營(yíng)養(yǎng)水平和種植模式下大豆光譜特性變化機(jī)理尚不完全明確, 還需要進(jìn)一步深入研究。

大豆-玉米套作模式作為近年來(lái)大豆增產(chǎn)的重要措施之一, 較以往大豆的單作模式具有提高土地利用率、減少肥料浪費(fèi)及比較效益高等特點(diǎn)[15]。但是由于高位作物玉米的蔭蔽作用, 套作低位作物大豆的冠層光環(huán)境較單作明顯不同, 對(duì)大豆產(chǎn)量的影響顯著[16]。高志英等[17]研究表明, 光照環(huán)境和施氮水平之間的互作對(duì)玉米氣孔特征的影響也是顯著的。同時(shí), 田間高光譜的測(cè)量由于受到高位作物玉米的蔭蔽, 對(duì)套作模式下大豆冠層光譜產(chǎn)生影響[18]。近年來(lái)還未見(jiàn)對(duì)套作和不同施氮水平大豆氮素遙感監(jiān)測(cè)研究的報(bào)道。

本研究利用單一波段、導(dǎo)數(shù)變換和兩波段自由組合等方法篩選大豆葉片氮素敏感波段, 構(gòu)建適合于單作套作大豆的氮素估測(cè)模型, 為不同種植模式下大豆氮素狀況監(jiān)測(cè)及氮肥精確管理提供理論支持。

1 材料與方法

1.1 試驗(yàn)地點(diǎn)與材料

2013年和2014年在仁壽縣四川現(xiàn)代糧食生產(chǎn)示范基地進(jìn)行試驗(yàn)。試驗(yàn)地土壤含有機(jī)質(zhì)12.96 g kg–1、全氮1.10 g kg–1、全磷0.68 g kg–1、全鉀14.66 g kg–1、堿解氮66.73 mg kg–1、速效磷3.26 mg kg–1、速效鉀178.74 mg kg–1, pH 6.85。供試大豆品種為南豆12, 玉米品種登海605。

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

采用二因素裂區(qū)試驗(yàn)設(shè)計(jì), 主因素為種植模式即套作和單作; 副因素為施純氮水平, 設(shè)N0為0 kg hm–2, N1為60 kg hm–2, N2為120 kg hm–2, N3為180 kg hm–2, 重復(fù)3次。

為保證土壤肥力的一致性, 前茬作物種植小麥勻地。寬窄行種植玉米, 寬行160 cm, 窄行40 cm, 3月下旬育苗移栽, 7月下旬收獲, 密度為6萬(wàn)株hm–2, 每穴單株。玉米底肥為氮肥(純N) 120 kg hm–2、磷肥(P2O5) 165 kg hm–2、鉀肥(K2O) 95 kg hm–2。玉米大喇叭口期追施氮肥(純氮) 120 kg hm–2。

大豆6月初播于玉米寬行內(nèi), 每幅種植2行, 穴距20 cm, 每穴2株, 行距40 cm。單作模式下, 大豆行距為50 cm, 穴距為40 cm, 每穴為2株。單作套作大豆密度一致, 均為10萬(wàn)株hm–2。兩種種植模式的小區(qū)面積均為6 m×6 m。大豆底肥為大豆底肥為磷肥(P2O5) 80 kg hm–2、鉀肥(K2O) 70 kg hm–2，大豆氮肥(純氮)按底肥∶始花期追肥=1∶1進(jìn)行配施。除草、噴藥等同大田操作。

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

單作和套作種植模式下在大豆苗期、分枝期、始花期、盛花期、始莢期和鼓粒期同步測(cè)定大豆冠層光譜和氮素積累量。

1.3.1 光譜測(cè)定 使用荷蘭AvaField-3型便攜式高光譜地物波譜儀測(cè)量光譜, 光譜范圍為350~1050 nm, 采樣間隔為0.6 nm, 視場(chǎng)角為25°。在6個(gè)生育期, 選擇晴朗無(wú)云, 11:00—14:00之間使用便攜式觀測(cè)架以距離大豆冠層1 m高度沿大豆條帶測(cè)定光譜, 測(cè)定前后及時(shí)參考白板的校正。每小區(qū)取值均為3次測(cè)量平均值。

1.3.2 氮素測(cè)定 對(duì)應(yīng)光譜測(cè)定位置(與光譜測(cè)定同步), 選取視場(chǎng)角范圍內(nèi)4株大豆葉片, 105℃殺青后在80℃烘干至恒重, 粉碎并過(guò)60目篩。用CE-440元素分析儀(美國(guó)加聯(lián)儀器有限公司生產(chǎn))測(cè)定葉片氮素。葉片氮素積累量(LNA)(g N m–2) = 干葉氮含量LNC (%)′干葉干物重LDW (g m–2)。每小區(qū)取值均為3次測(cè)量平均值。

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

采用Microsoft Excel 2013整理數(shù)據(jù), 使用Origin 9.0軟件繪制; 利用DPS 7.05及Matlab程序?qū)?shù)據(jù)進(jìn)行相關(guān)性分析, 模型構(gòu)建(2013年= 88)及驗(yàn)證(2014年= 40)。在大豆氮素積累量及光譜反射率變化規(guī)律的基礎(chǔ)上進(jìn)行分析, 探究氮素積累量與光譜指數(shù)間的相關(guān)性。本文引用的光譜指數(shù)如表1所示。

表1 本文引用的光譜指數(shù)

是相應(yīng)波長(zhǎng)的光譜反射率。

is the spectral reflectance rate at the corresponding wavelength.

2 結(jié)果與分析

2.1 不同種植模式和氮肥水平下大豆葉片氮素積累量變化規(guī)律

由表2可知, 對(duì)于同一施氮水平, 隨著生育時(shí)期的推進(jìn), 單作套作種植模式下的大豆葉片氮素積累量均呈現(xiàn)先增加后減少的趨勢(shì), 在結(jié)莢期達(dá)到最大, 且兩年表現(xiàn)出相同規(guī)律。而在不同的施肥水平和同一生育時(shí)期下, 隨著施氮量的增加, 葉片氮素積累量總體上均呈現(xiàn)增加的趨勢(shì); N3處理葉片氮素積累量最大, 2013年和2014年套作大豆冠層葉片氮素積累量最大值分別為9.36 g m–2和7.40 g m–2; 而單作大豆分別為8.79 g m–2和8.61 g m–2。通過(guò)對(duì)不同生育時(shí)期單作套作模式下大豆葉片氮素積累量的比較, 表明單作下氮素積累量的最快增長(zhǎng)時(shí)期為始花期到盛花期, 而套作條件下為盛花期到結(jié)莢期, 這可能歸因于本試驗(yàn)的套作模式下玉米收獲時(shí)大豆處于始花期, 始花期之前玉米對(duì)大豆的蔭蔽使得大豆生長(zhǎng)受到一定程度抑制, 導(dǎo)致葉片氮素積累量相對(duì)單作偏低, 而玉米收獲后蔭蔽解除, 受光影響的生理活動(dòng)隨冠層光環(huán)境的恢復(fù)在一定程度上逐步得到補(bǔ)償, 大豆在始花期到盛花期開始彌補(bǔ)前期蔭蔽所帶來(lái)的影響, 進(jìn)而可能使得根系對(duì)氮素的吸收增加, 但這個(gè)恢復(fù)過(guò)程是相對(duì)緩慢的, 因此造成單作套作下氮素最快增長(zhǎng)時(shí)期的差異。

表2 不同施氮水平下單作和套作大豆葉片氮素積累量變化

2.2 大豆冠層光譜變化特征

隨著生育時(shí)期的推進(jìn), 大豆冠層光譜反射率發(fā)生了明顯變化。如圖1所示(各處理平均值), 700~1000 nm內(nèi), 處理之間光譜反射率的差異變化十分明顯, 兩年的數(shù)據(jù)表明, 除2013年鼓粒期, 單作種植下大豆不同生育時(shí)期的冠層光譜反射率均高于套作; 而隨著生育時(shí)期的推進(jìn), 反射率均呈現(xiàn)先升高后降低的規(guī)律, 與大豆氮素積累量隨生育時(shí)期變化規(guī)律基本一致,

對(duì)于一階導(dǎo)數(shù)光譜, 光譜中紅邊位置(700 nm左右)是反映作物長(zhǎng)勢(shì)的敏感特征。圖2所示, 隨生育時(shí)期的推進(jìn), 光譜紅邊位置先后出現(xiàn)“紅移”及“藍(lán)移”, 且紅邊幅值也表現(xiàn)出先增大后降低之勢(shì), 同時(shí), 兩種變化之拐點(diǎn)在同一生育時(shí)期; 而以上規(guī)律均與相應(yīng)種植模式下大豆氮素積累量隨生育時(shí)期變化規(guī)律一致。套作模式下, 紅邊位置及紅邊幅值的拐點(diǎn)出現(xiàn)在結(jié)莢期, 單作下, 拐點(diǎn)出現(xiàn)的時(shí)期變化與冠層原始光譜反射率變化一致.

因此基于光譜原始反射率及其一階導(dǎo)數(shù)的變化規(guī)律與相應(yīng)的大豆葉片氮素積累量變化比較不難看出, 在700~1000 nm范圍內(nèi), 大豆原始光譜反射率及其一階導(dǎo)數(shù)與氮素的變化密切相關(guān), 但大豆氮素積累與葉片內(nèi)含氮有機(jī)物及含氮無(wú)機(jī)物的積累量又密切相關(guān), 而各物質(zhì)的光譜吸收波段不盡相同, 因此導(dǎo)致葉片氮素積累敏感波段不同于單一物質(zhì)那么簡(jiǎn)單。

圖1 單作和套作大豆葉片原始光譜變化規(guī)律

圖2 單作和套作大豆葉片一階導(dǎo)數(shù)光譜變化

2.3 大豆葉片氮素積累量與光譜指數(shù)的相關(guān)性分析

2.3.1 原始及一階導(dǎo)數(shù)光譜與氮素積累量相關(guān)性

圖3表明, 當(dāng)波長(zhǎng)小于720 nm時(shí), 葉片氮素積累量與原始光譜呈現(xiàn)負(fù)相關(guān), 其中在456~710 nm波段內(nèi)表現(xiàn)為顯著負(fù)相關(guān), 最小相關(guān)系數(shù)在692 nm處(= –0.497); 波長(zhǎng)處于720~1000 nm時(shí), 葉片氮素積累量與原始光譜呈現(xiàn)正相關(guān), 其中在730~1000 nm波段內(nèi)達(dá)到0.05顯著性水平, 最大相關(guān)系數(shù)在892 nm處(= 0.55)。相比于原始光譜, 一階導(dǎo)數(shù)光譜與葉片氮素積累量的相關(guān)性復(fù)雜得多, 其中在710~763 nm波段內(nèi), 一階導(dǎo)數(shù)光譜與葉片氮素積累量呈現(xiàn)0.05水平的顯著正相關(guān), 最大值在751 nm處(= 0.89); 在680~696 nm波段內(nèi), 一階導(dǎo)數(shù)光譜與葉片氮素積累量呈顯著負(fù)相關(guān), 最小值在688 nm處(= –0.61)。相關(guān)系數(shù)的最大和最小值都在紅邊范圍內(nèi), 處在冠層葉片原始光譜變化最大的波段內(nèi), 突出和強(qiáng)化了冠層原始光譜的信息。

圖3 大豆葉片氮素積累量與原始和一階導(dǎo)數(shù)光譜相關(guān)分析

2.3.2 光譜指數(shù)與氮素積累量相關(guān)性 通過(guò)對(duì)表1中列舉的光譜指數(shù)與本試驗(yàn)中大豆冠層葉片氮素積累量的相關(guān)性分析, 結(jié)果如表3所示, 除黃邊面積外, 其他光譜指數(shù)與氮素的相關(guān)性都達(dá)0.05水平上的顯著相關(guān), 且其中紅邊與藍(lán)邊面積之比達(dá)0.01水平的極顯著相關(guān)(0.871)。從圖3可以發(fā)現(xiàn), 與氮素敏感的波段很多, 因此基于波段自由組合, 構(gòu)建基于RSI、DSI和NDSI光譜指數(shù)與氮素相關(guān)性最優(yōu)的波段組合(圖4)。

基于原始光譜的相關(guān)系數(shù)最高的各光譜指數(shù)分別為RSI (767, 745)(= 0.901)、DSI (771, 755)(= 0.92)和NDSI (767, 745)(= 0.901); 基于一階導(dǎo)數(shù)光譜的相關(guān)系數(shù)最高的各光譜指數(shù)分別為RSI (D735, D713) (= 0.897)、DSI (D748, D508)(= 0.917)和NDSI (D748, D696)(= 0.890)。通過(guò)對(duì)基于原始光譜和一階導(dǎo)數(shù)光譜自由組合的光譜指數(shù)與大豆葉片氮素積累量的相關(guān)系數(shù)分析發(fā)現(xiàn), 大豆葉片氮素積累量的敏感區(qū)域集中于原始光譜中的700~800 nm的近紅外波段, 導(dǎo)數(shù)光譜集中于可見(jiàn)光和近紅外波段。其中在整個(gè)400~1000 nm波段范圍內(nèi), 基于原始光譜自由組合的光譜指數(shù)好于基于一階導(dǎo)數(shù)光譜自由組合的光譜指數(shù), 但兩者的最大相關(guān)系數(shù)值差別不大。

圖4 大豆葉片氮素積累量與光譜指數(shù)的相關(guān)性分析

2.4 大豆葉片氮素估測(cè)模型的構(gòu)建與檢驗(yàn)

2.4.1 大豆葉片氮素估測(cè)模型的構(gòu)建 基于前文篩選出相關(guān)系數(shù)達(dá)到0.900極顯著水平(<0.01)的光譜變量, 采用線性和非線性回歸方法建立估算大豆葉片氮素積累量的模型, 本試驗(yàn)采用以下5種擬合模型構(gòu)建估算模型。

線性函數(shù):=+

對(duì)數(shù)函數(shù):=+′ln ()

表3 大豆葉片氮素積累量與光譜指數(shù)的相關(guān)分析

*＜0.05 ;**＜0.01.

二次函數(shù):=++2

乘冪函數(shù):=′x

指數(shù)函數(shù):=′e

其中,代表大豆葉片的氮素積累量值,代表光譜指數(shù),、、均為常數(shù)。結(jié)果如表4所示, 對(duì)光譜指數(shù)最佳的擬合模型是基于DSI (771, 755)構(gòu)建的冪函數(shù)和線性函數(shù), 決定系數(shù)分別為0.872和0.847, 其中冪函數(shù)的值達(dá)到最大, 為588.223, 其次是線性函數(shù), 為475.534。

2.4.2 大豆葉片氮素估測(cè)模型的測(cè)試與檢驗(yàn) 將2014年的數(shù)據(jù)(= 40)帶入到上述所建模型, 結(jié)果如圖5顯示, 以DSI (771, 755)構(gòu)建的線性和乘冪兩模型來(lái)檢驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)精度, 選擇模型擬合2和預(yù)測(cè)2都較大。基于DSI (771, 755)的線性模型的2、RMSE和RPD分別為0.827、1.281、2.416; 乘冪模型的2、RMSE和RPD分別為0.800、1.389、2.227?；贒SI (771, 755)構(gòu)建的線性和乘冪模型的RPD均大于2, 說(shuō)明這2個(gè)模型的預(yù)測(cè)精度較好。

3 討論

氮素是作物生長(zhǎng)發(fā)育及產(chǎn)量形成的重要元素之一, 相關(guān)研究表明, 施氮水平的不同對(duì)油菜[19]、小麥[20]、玉米[21]等作物產(chǎn)量都有著重要的影響。同時(shí), 閆艷紅等[22]和劉小明等[23]的研究表明, 套作模式下, 施氮水平的不同對(duì)于大豆生育后期干物質(zhì)積累及產(chǎn)量形成都有著十分重要的影響。此外, 田艷洪等[24]的研究表明, 不同時(shí)期施氮對(duì)大豆的產(chǎn)量也有影響, 其中以鼓粒初期追肥對(duì)大豆產(chǎn)量的增加最為顯著。大豆的整個(gè)生育時(shí)期, 氮素主要來(lái)自根瘤固氮及人工使用氮肥, 因此, 根據(jù)大豆的氮素營(yíng)養(yǎng)水平, 及時(shí)、適量追施氮肥對(duì)于提高大豆產(chǎn)量及減少肥料的濫用是十分重要的。

表4 葉片氮素積累量與光譜參數(shù)的線性和非線性回歸分析

圖5 高光譜估測(cè)葉片氮素積累量的回歸分析

籽粒形成過(guò)程中物質(zhì)的積累主要來(lái)自葉片中物質(zhì)的轉(zhuǎn)運(yùn), 因此監(jiān)測(cè)大豆葉片的氮素水平有助大豆氮素的精準(zhǔn)調(diào)控。冠層葉片的光譜反射率與葉片含氮量有關(guān)[8,25], 而以紅邊特征參數(shù)能夠測(cè)定氮素含量[26]。本實(shí)驗(yàn)所構(gòu)建的預(yù)測(cè)模型與梁留鎖等[27]研究結(jié)果相比, 雖冪函數(shù)及線性函數(shù)模型的決定系數(shù)都略低, 但大田光譜測(cè)量受多種因素影響[28]。同時(shí), 玉-豆套作模式下, 大豆冠層有效光合輻射及透光率均顯著低于單作, 且受行間距不同配置及玉-豆套作不同組合的影響[16,29], 而陰影對(duì)光譜的反射也有影響。本試驗(yàn)所構(gòu)建模型很好估測(cè)單作套作模式下大豆冠層葉片氮素積累量, 增加了模型在不同環(huán)境下的適用性, 對(duì)高效指導(dǎo)農(nóng)業(yè)生產(chǎn)中大豆氮素監(jiān)測(cè)及施肥指導(dǎo)有著十分重要的理論支持。

本次試驗(yàn)中, 研究了單作與套作模式下大豆葉片的氮素估測(cè), 但尚未分析研究套作模式不同行間距和不同套作組合下的大豆葉片氮素估測(cè)。因此, 在后續(xù)的實(shí)驗(yàn)研究中, 需進(jìn)一步深入研究不同行距配置和組合下的廣適性遙感估測(cè)模型。

4 結(jié)論

充分利用高光譜的信號(hào)豐富、觀測(cè)能力強(qiáng)和即時(shí)性的特點(diǎn), 獲取了大豆在不同種植模式和氮肥差異下的不同生育時(shí)期光譜特征參數(shù)和葉片氮素積累量。基于波段優(yōu)選的DSI (771, 755)構(gòu)建的線性(= –1.249+3.209)和乘冪(= –1.4701.676)模型能較好預(yù)測(cè)單作套作兩種模式下大豆冠層葉片不同生育時(shí)期的氮素狀況。

[1] 李衛(wèi)國(guó), 李秉柏, 王志明, 張婭香, 黃曉軍. 作物長(zhǎng)勢(shì)遙感監(jiān)測(cè)應(yīng)用研究現(xiàn)狀和展望. 江蘇農(nóng)業(yè)科學(xué), 2006, (3): 12–15Li W G, Li B B, Wang Z M, Zhang Y C, Huang X J. Status and prospect on research and application of crop condition monitoring using remote sensing., 2006, (3): 12–15 (in Chinese with English abstract)

[2] 吳華兵, 朱艷, 田永超, 姚霞, 劉曉軍, 周治國(guó), 曹衛(wèi)星. 棉花冠層高光譜指數(shù)與葉片氮積累量的定量關(guān)系. 作物學(xué)報(bào), 2007, 33: 518–522 Wu H B, Zhu Y, Tian Y C, Yao X, Liu X J, Zhou Z G, Cao W X. Relationship between canopy hyperspectral index and leaf nitrogen accumulation in cotton., 2007, 33: 518–522 (in Chinese with English abstract)

[3] 馮偉, 王曉宇, 宋曉, 賀利, 王永華, 郭天財(cái). 基于冠層反射光譜的小麥白粉病嚴(yán)重度估測(cè). 作物學(xué)報(bào), 2013, 39: 1469–1477 Feng W, Wang X Y, Song X, He L, Wang Y H, Guo T C. Estimation of severity level of wheat powdery mildew based on ca-nopy spectral reflectance., 2013, 39: 1469–1477 (in Chinese with English abstract)

[4] 胡軍林.主要農(nóng)作物缺氮癥狀及防治措施. 農(nóng)技服務(wù), 2008, 25(1): 54, 108 Hu J L.Causes of nitrogen deficiency in main crops and preventive measures., 2008, 25(1): 54, 108 (in Chinese with English abstract)

[5] 郝小雨, 周寶庫(kù), 馬星竹, 高中超.長(zhǎng)期不同施肥措施下黑土作物產(chǎn)量與養(yǎng)分平衡特征. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(16): 178–185 Hao X Y, Zhou B K, Ma X Z, Gao Z C. Characteristics of crop yield and nutrient balance under different long-term fertilization practices in black soil., 2015, 31(16): 178–185 (in Chinese with English abstract)

[6] Bekele W. Effect of soybean varieties and nitrogen fertilizer rates on yield, yield components and productivity of associated crops under maize/soybean intercropping at Mechara, Eastern Ethiopia., 2016, 5(1): 1

[7] 王勁松, 焦曉燕, 丁玉川, 董二偉, 白文斌, 王立革, 武愛(ài)蓮.粒用高粱養(yǎng)分吸收、產(chǎn)量及品質(zhì)對(duì)氮磷鉀營(yíng)養(yǎng)的響應(yīng). 作物學(xué)報(bào), 2015, 41: 1269–1278Wang J S, Jiao X Y, Ding Y C, Dong E W, Bai W B, Wang L G, Wu A L. Response of nutrient uptake, yield and quality of grain sorghum to nutrition of nitrogen, phosphorus and potassium., 2015, 41: 1269–1278 (in Chinese with English abstract)

[8] 胡珍珠, 潘存德, 肖冰, 潘鑫.基于光譜特征參量的核桃葉片氮素含量估測(cè)模型. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(9): 180–186 Hu Z Z, Pan C D, Xiao B, Pan X. Spectral characteristic parameter-based models for foliar nitrogen concentration estimation of., 2015, 31(9): 180–186 (in Chinese with English abstract)

[9] Walburg G, Bauer M E, Daughtry C S T. Effects of nitrogen nutrition on the growth, yield, and reflectance characteristics of corn canopies., 1982, DOI: 10.2134/agronj1982.0002196200 7400040020x

[10] 岳延濱, 聶克艷, 黎瑞君, 李莉婕, 孫長(zhǎng)青, 彭志良, 趙澤英.不同施氮水平辣椒地上部全氮含量與冠層光譜反射率的相關(guān)性. 貴州農(nóng)業(yè)科學(xué), 2014, (11): 244–247 Yue Y B, Nie K Y, Li Y J, Li L J, Sun C Q, Peng Z L, Zhao Z Y. Correlation of total nitrogen content in overground part of pepper at different nitrogen application level with canopy spectral reflectance., 2014, (11): 244–247 (in Chinese with English abstract)

[11] 田永超, 楊杰, 姚霞, 曹衛(wèi)星, 朱艷.利用葉片高光譜指數(shù)預(yù)測(cè)水稻群體葉層全氮含量. 作物學(xué)報(bào), 2010, 36: 1529–1537 Tian Y C, Yang Jie, Yao X, Cao W X, Zhu Y. Monitoring canopy leaf nitrogen concentration based on leaf hyperspectral indices in rice., 2010, 36: 1529–1537 (in Chinese with English abstract)

[12] Wei F, Zhang H Y, Zhang Y S, Qi S L, Heng Y R, Guo B B, Ma D Y, Guo T C. Remote detection of canopy leaf nitrogen concentration in winter wheat by using water resistance vegetation indices fromhyperspectral data., 2016, 198: 238–246

[13] 宋英博.不同施氮量對(duì)大豆蛋白質(zhì)和脂肪含量的影響. 黑龍江農(nóng)業(yè)科學(xué), 2010, (7): 52–53 Song Y B. Effect of different nitrogen application on protein and fat content in soybean., 2010, (7): 52–53 (in Chinese with English abstract)

[14] 董守坤, 龔振平, 祖?zhèn)? 氮素營(yíng)養(yǎng)水平對(duì)大豆氮素積累及產(chǎn)量的影響. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2010, 16: 65–70 Dong S K, Gong Z P, Zu W. Effects of nitrogen nutrition levels on N-accumulation and yields of soybean., 2010, 16: 65–70 (in Chinese with English abstract)

[15] 楊文鈺, 雍太文, 任萬(wàn)軍, 樊高瓊, 牟錦毅, 盧學(xué)蘭. 發(fā)展套作大豆, 振興大豆產(chǎn)業(yè).大豆科學(xué), 2008, 27: 1–7 Yang W Y, Yong T W, Ren W J, Fan G Q, Mu J Y, Lu X L.Development of soybeans, the revitalization of soybean industry., 2008, 27: 1–7 (in Chinese with English abstract)

[16] 崔亮, 蘇本營(yíng), 楊峰, 楊文鈺.帶狀套作大豆群體冠層光能截獲與利用特征. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48: 43–54 Cui L, Su B Y, Yang F, Yang W Y. Relationship between light interception and light utilization of soybean canopy in relay strip intercropping system., 2015, 48: 43–54 (in Chinese with English abstract)

[17] 高志英, 丁圣彥, 谷艷芳, 邢倩.不同光環(huán)境與氮肥互作對(duì)玉米氣孔特征的影響. 河南農(nóng)業(yè)科學(xué), 2008, 37(9): 15–19 Gao Z Y, Ding S Y, Gu Y F, Xing Q. The interactive effects of light condition and nitrogen supply on stomatal characteristics of maize., 2008, 37(9): 15–19 (in Chinese with English abstract)

[18] 何挺, 劉榮, 王靜.野外波譜測(cè)量的影響因素研究.地理與地理信息科學(xué), 2003, 19(5): 6–10 He T, Liu R, Wang J. The influences factors on field spectrometry., 2003, 19(5): 6–10 (in Chinese with English abstract)

[19] 楊德金, 彭守華. 秦優(yōu)12不同播期、密度和施氮水平對(duì)產(chǎn)量的影響. 安徽農(nóng)學(xué)通報(bào), 2010, 16(13): 100–102 Yang D J, Peng S H.Effects of different sowing date, density and nitrogen application on yield., 2010, 16(13): 100–102 (in Chinese)

[20] 郭天財(cái), 宋曉, 馬冬云, 馮偉, 王永華.施氮水平對(duì)2種穗型冬小麥品種產(chǎn)量及氮素吸收利用的影響. 西北植物學(xué)報(bào), 2008, 28: 554–558 Guo T C, Song X, Ma D Y, Feng W, Wang Y H. Effects of nitrogen application on grain yield, n absorption and utilization rate in winter wheat with two spike-types., 2008, 28: 554–558 (in Chinese with English abstract)

[21] 寧堂原, 焦念元, 李增嘉, 張民, 趙春, 韓賓, 邵國(guó)慶.施氮水平對(duì)不同種植制度下玉米氮利用及產(chǎn)量和品質(zhì)的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2006, 17: 2332–2336 Ning T Y, Jiao N Y, Li Z J, Zhang M, Zhao C, Han B, Shao G Q. Effects of N application rate on N utilization, yield and quality of maize under different cropping systems., 2006, 17: 2332–2336 (in Chinese with English abstract)

[22] 閆艷紅, 楊文鈺, 張新全, 陳小林, 陳忠群.施氮量對(duì)套作大豆花后光合特性、干物質(zhì)積累及產(chǎn)量的影響. 草業(yè)學(xué)報(bào), 2011, 20: 233–238 Yan Y H, Yang W Y, Zhang X Q, Chen X L, Chen Z Q. Effects of different nitrogen levels on photosynthetic characteristics, dry matter accumulation and yield of relay strip intercroppingafter blooming., 2011, 20: 233–238 (in Chinese with English abstract)

[23] 劉小明, 雍太文, 蘇本營(yíng), 劉文鈺, 周麗, 宋春, 楊峰, 王小春, 楊文鈺.減量施氮對(duì)玉米-大豆套作系統(tǒng)中作物產(chǎn)量的影響. 作物學(xué)報(bào), 2014, 40: 1629–1638 Liu X M, Yong T W, Su B Y, Liu W Y, Zhou L, Song C, Yang F, Wang X C, Yang W Y. Effect of reduced N application on crop yield in maize-soybean intercropping system., 2014, 40: 1629–1638 (in Chinese with English abstract)

[24] 田艷洪, 劉元英, 張文釗, 羅盛國(guó).不同時(shí)期施用氮肥對(duì)大豆根瘤固氮酶活性及產(chǎn)量的影響. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2008, 39(5): 15–19 Tian Y H, Liu Y Y, Zhang W Z, Luo S G.Effect of N fertilization at different stage on nitrogenase activity and yield of soybean., 2008, 39(5): 15–19 (in Chinese with English abstract)

[25] 周麗麗, 馮漢宇, 閻忠敏, 劉克, 周順利.玉米葉片氮含量的高光譜估算及其品種差異. 農(nóng)業(yè)工程學(xué)報(bào), 2010, 26(8): 195–199 Zhou L L, Feng H Y, Yan Z M, Liu K, Zhou S L. Hyperspectral diagnosis of leaf N concentration of maize and variety difference., 2010, 26(8): 195–199 (in Chinese with English abstract)

[26] 馮偉, 朱艷, 姚霞, 田永超, 郭天財(cái), 曹衛(wèi)星.利用紅邊特征參數(shù)監(jiān)測(cè)小麥葉片氮素積累狀況. 農(nóng)業(yè)工程學(xué)報(bào), 2009, 25(11): 194–201 Feng W, Zhu Y, Yao X, Tian Y C, Guo T C, Cao W X. Monitoring nitrogen accumulation in wheat leaf with red edge characteristics parameters., 2009, 25(11): 194–201 (in Chinese with English abstract)

[27] 梁留鎖, 馬旭, 喬欣.大豆葉片氮素含量檢測(cè)技術(shù)的試驗(yàn)研究. 農(nóng)機(jī)化研究, 2010, 32(11): 121–124Liang L S, Ma X, Qiao X. The total leaf nitrogen determination of the soybean’s canopy based on spectral reflection., 2010, 32(11): 121–124 (in Chinese with English abstract)

[28] 何挺, 程燁, 王靜.野外地物光譜測(cè)量技術(shù)及方法. 中國(guó)土地科學(xué), 2002, 16(5): 30–36He T, Cheng Y, Wang J. The technology and method of field spectrometry., 2002, 16(5): 30–36 (in Chinese with English abstract)

[29] 崔亮, 蘇本營(yíng), 楊峰, 楊文鈺. 不同玉米-大豆帶狀套作組合條件下光合有效輻射強(qiáng)度分布特征對(duì)大豆光合特性和產(chǎn)量的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47: 1489–1501 Cui L, Su B Y, Yang F, Yang W Y. Effects of photosynthetically active radiation on photosynthetic characteristics and yield of soybean in different maize/soybean relay strip intercropping systems., 2014, 47: 1489–1501 (in Chinese with English abstract)

Remote Detection of Canopy Leaf Nitrogen Status in Soybean by Hyperspectral Data under Monoculture and Intercropping Systems

CHEN Jun-Xu, HUANG Shan, FAN Yuan-Fang, WANG Rui, LIU Qin-Lin, YANG Wen-Yu*, and YANG Feng*

Key Laboratory of Crop Eco-physiology and Farming System in Southwest, the Ministry of Agriculture / College of Agronomy, Sichuan Agricultural University / Sichuan Engineering Research Center for Crop Strip Intercropping System, Chengdu 611130, China

Non-destructive monitoring of soybean nitrogen status is important for precise N management in soybean production. In this study, the quantitative correlation between soybean leaf nitrogen status and canopy hyperspectral reflectance was investigated. Field experiments were conducted. With four nitrogen application rates for two years under monoculture and intercropping systems. The nitrogen accumulation of canopy leaves showed a single-peak changing trend in the process of soybean growth. The maximum value in monoculture and intercropping was 8.70 g m–2and 8.38 g m–2, respectively, at pod stage under N3 treatment. The raw hyperspectral reflectance and the leaf nitrogen accumulation had the same changing trend at different growth stages with different planting patterns. The peak value of the raw hyperspectral reflectance in the 700–1000 nm occurred at pod stage. In the first-order derivative spectrum, the red edge amplitude values increased first and then decreased. The position of the red edge changed as “Red shift” and “Blue shift” with the increase or decrease of leaf nitrogen accumulation. The results of the correlation analysis showed that the linear model and the power model by using the Difference Spectral Index (DSI: 771, 755) based on the best spectral band combination (BSBC) had the greatest accuracy to estimate the leaf nitrogen status of soybean.

Soybean; Nitrogen accumulation; Hyperspectral reflectance; Model

10.3724/SP.J.1006.2017.01835

本研究由國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFD0300602)資助。

The work was supported by the National Key Research and Development Program of China (2016YFD0300602).

楊峰, E-mail: f.yang@sicau.edu.cn; 楊文鈺, E-mail: mssiyangwy@sicau.edu.cn

E-mail: m17740901718@163.com

2017-03-09; Accepted(接受日期): 2017-05-10; Published online(網(wǎng)絡(luò)出版日期): 2017-06-08.

http://kns.cnki.net/kcms/detail/11.1809.S.20170608.1000.002.html