基于快速葉綠素?zé)晒饧夹g(shù)的油菜冠層生化參數(shù)垂直異質(zhì)性分析

張佳菲 萬(wàn)亮 何勇 岑海燕

摘要： 準(zhǔn)確獲取作物冠層生化信息對(duì)監(jiān)測(cè)作物生長(zhǎng)和指導(dǎo)精準(zhǔn)施肥具有重要意義?，F(xiàn)有的作物生化參數(shù)的垂直分布研究以高光譜遙感反演為主，缺乏與光合生理的聯(lián)系。本研究主要探究了不同氮素處理水平下油菜苗期冠層內(nèi)的葉綠素、類胡蘿卜素、干物質(zhì)和水分等生化參數(shù)的垂直分布變化特性，同時(shí)利用快速葉綠素?zé)晒饧夹g(shù)測(cè)定了葉片的光合性能，并通過(guò)線性回歸分析和主成分分析進(jìn)一步剖析了熒光響應(yīng)與生化參數(shù)的內(nèi)在聯(lián)系。試驗(yàn)結(jié)果表明：（1）苗期中期油菜冠層的葉綠素含量、類胡蘿卜素含量、干物質(zhì)和水分含量均呈拋物線型的非均勻垂直分布，而葉綠素與類胡蘿卜素的比值具有與其他生化參數(shù)不同的垂直分布模式，其隨著葉位升高和施氮量的增加逐漸下降，與推動(dòng)力DFTotal、電子鏈末端量子產(chǎn)額φRo等熒光參數(shù)的垂直分布模式相同;（2）熒光參數(shù)，特別是DFTotal，對(duì)油菜葉片葉綠素與類胡蘿卜的比值、葉綠素和干物質(zhì)含量具有較強(qiáng)的評(píng)估能力;（3）缺氮會(huì)降低苗期油菜葉片的光系統(tǒng)I和II（PSI和PSII）性能，通過(guò)最大光化學(xué)效率φPo等熒光參數(shù)可對(duì)氮素脅迫進(jìn)行診斷;而不同葉位葉片在PSI性能即電子末端傳遞效率上具有顯著差異，通過(guò)DFTotal可有效表征冠層生化參數(shù)的垂直異質(zhì)性。上述結(jié)果表明，應(yīng)用快速葉綠素?zé)晒饧夹g(shù)對(duì)作物進(jìn)行生化信息的垂直異質(zhì)性監(jiān)測(cè)具有可行性，可為指導(dǎo)精準(zhǔn)施肥和提高優(yōu)質(zhì)優(yōu)產(chǎn)提供新思路和技術(shù)支撐。

關(guān)鍵詞： 快速葉綠素?zé)晒?JIP-測(cè)定;垂直異質(zhì)性;葉綠素;類胡蘿卜素;氮素

中圖分類號(hào)： S363??? ??? 文獻(xiàn)標(biāo)志碼： A??? 文章編號(hào)： 202103-SA005

引用格式：張佳菲， 萬(wàn)亮， 何勇， 岑海燕. 基于快速葉綠素?zé)晒饧夹g(shù)的油菜冠層生化參數(shù)垂直異質(zhì)性分析[J]. 智慧農(nóng)業(yè)（中英文）， 2021， 3 （1）： 40-50.

Citation：ZHANG Jiafei， WAN Liang， HE Yong， CEN Haiyan. Vertical heterogeneity analysis of biochemical parameters in oilseed rape canopy based on fast chlorophyll fluorescence technology[J]. Smart Agriculture， 2021， 3 （1）： 40-50. （

1? 引? 言

油菜是世界上僅次于大豆的第二大油料作物，在全球石油資源替代、動(dòng)物飼料和生物燃料供應(yīng)中發(fā)揮著重要作用[1]。油菜生化表型信息（例如色素、干物質(zhì)、水分等）的快速獲取對(duì)于優(yōu)良品質(zhì)選育、提高作物品質(zhì)和產(chǎn)量具有重要意義。而葉片生化特性在冠層內(nèi)往往呈現(xiàn)顯著的非均勻垂直分布，例如，光合作用較為活躍的上層葉片中的氮濃度往往較高[2，3]。作物冠層內(nèi)的這種空間異質(zhì)性可以被視為一種在有限的營(yíng)養(yǎng)資源和變化環(huán)境條件中不斷進(jìn)行葉片之間的調(diào)整分配，以最大限度地提高冠層光合速率的適應(yīng)性策略[2，4]。因此，有必要對(duì)作物冠層內(nèi)生化組分信息的垂直異質(zhì)性進(jìn)行定量研究與分析，為實(shí)施油菜生長(zhǎng)信息動(dòng)態(tài)監(jiān)測(cè)和精準(zhǔn)科學(xué)施肥提供有效的理論方法和技術(shù)支撐。

目前，對(duì)于作物生化組分垂直分布的研究以高光譜遙感評(píng)估作物的氮素垂直分布為主，通過(guò)統(tǒng)計(jì)模型建立冠層光譜與各個(gè)垂直層氮素的聯(lián)系[5-8]。然而，基于反射率的監(jiān)測(cè)方法往往包含了來(lái)自植物和土壤的混合信號(hào)，降低了反演的準(zhǔn)確性和穩(wěn)健性，且與作物光合作用等生理機(jī)制的聯(lián)系并不是很清楚[9，10]。此外，對(duì)于光合色素等與植物光合生理聯(lián)系十分緊密的生化組分的垂直分布研究相對(duì)較少，無(wú)法滿足田間作物優(yōu)良品質(zhì)選育和關(guān)鍵生育期施氮管理等依賴高通量表型信息獲取的要求。

在植物光合作用中，葉綠素分子吸收的能量一般有三種轉(zhuǎn)化形式：光化學(xué)、熱耗散和葉綠素?zé)晒鈁11]。也就是說(shuō)，未被光化學(xué)反應(yīng)消耗或轉(zhuǎn)化為熱量的剩余光能則以熒光的形式從植物中輻射出來(lái)。而植物衰老或逆境脅迫、高溫低溫、鹽脅迫以及干旱等都會(huì)引起這三種能量發(fā)生變化[12，13]。因此，通過(guò)葉綠素?zé)晒獾淖兓梢蕴綔y(cè)到植物葉片生理狀態(tài)的變化，如光能的吸收與轉(zhuǎn)換、反應(yīng)中心的狀態(tài)、過(guò)剩光能及其耗散、光系統(tǒng)II（Photosystem II，PSII）供體側(cè)和受體的活性、光合作用光抑制與光破壞等[14，15]。

與傳統(tǒng)的氣體交換測(cè)量方法相比，葉綠素?zé)晒饧夹g(shù)已經(jīng)成為一種了解光合功能、生理特性的時(shí)空動(dòng)態(tài)變化更有效、便捷且高通量的方法。其優(yōu)勢(shì)主要在于熒光信號(hào)只來(lái)自于植物本身，因此可以有效避免來(lái)自于土壤等其他信號(hào)的干擾，從而直接反映植物最真實(shí)的生理狀態(tài)。除了植物逆境生理響應(yīng)的檢測(cè)外，已有許多研究將葉綠素?zé)晒鈶?yīng)用于分析作物葉綠素含量、氮營(yíng)養(yǎng)指數(shù)（Nitrogen Nutrient Index，NNI）等生理生態(tài)指標(biāo)來(lái)指導(dǎo)作物生長(zhǎng)過(guò)程中的施肥。熒光參數(shù)會(huì)受氮素等營(yíng)養(yǎng)元素供應(yīng)的顯著影響，隨著施氮量增加，水稻中非光化學(xué)淬滅（Non-Photochemical Quenching，NPQ）也會(huì)增加[16]，而在棉花中過(guò)量施氮?jiǎng)t使作物最大光化學(xué)效率（Fv /Fm）和電子傳遞速率（Electron Transport Rate，ETR）等顯著下降[17]。此外，一些研究也利用熒光參數(shù)在不同葉層觀測(cè)到PSII活性和能量利用效率的顯著差異[18，19]，說(shuō)明了葉綠素?zé)晒忭憫?yīng)在對(duì)作物光合生理表征，特別是垂直異質(zhì)性表征的敏感性和巨大潛能。而目前對(duì)于作物冠層內(nèi)熒光參數(shù)的垂直異質(zhì)性的研究還較少。因此，利用葉綠素?zé)晒饧夹g(shù)探索作物冠層內(nèi)熒光參數(shù)的垂直異質(zhì)性以及與生化信息的關(guān)系，并作為一種評(píng)價(jià)作物生化參數(shù)的替代方法具有非常大的前景。

本研究主要用分光光度計(jì)測(cè)定并分析了油菜苗期冠層葉綠素、類胡蘿卜素、干物質(zhì)和水分等生化參數(shù)的垂直分布變化特性，同時(shí)利用快速葉綠素?zé)晒饧夹g(shù)探究了冠層內(nèi)光合生理特性的垂直異質(zhì)性，剖析了熒光響應(yīng)與生化參數(shù)的內(nèi)在聯(lián)系與生理機(jī)制，最終確定了最適宜表征油菜生化信息的熒光參數(shù)。

2? 材料與方法

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

本研究的試驗(yàn)對(duì)象為甘藍(lán)型油菜品種浙雙758（Brassica napus L.），于浙江大學(xué)紫金港校區(qū)農(nóng)業(yè)試驗(yàn)站試驗(yàn)田供試土壤設(shè)置3個(gè)氮素水平，分別為不施氮（N0，0 kg N/ha）、正常施氮（N1，215 kg N/ha）和過(guò)量施氮（N2，430 kg N/ha），即N0：N1：N2= 0：1：2。氮肥以3：1：2的比例分別施于油菜移栽前、苗期初期和蕾苔期，磷肥和鉀肥則均作基肥施用。試驗(yàn)田采用隨機(jī)區(qū)組排列，每個(gè)氮素水平3個(gè)重復(fù)，共9個(gè)小區(qū)。小區(qū)面積為6.6 m2（5.5 m×1.2 m），每個(gè)小區(qū)南北和東西方向各設(shè)置寬0.4 m的保護(hù)道。油菜于2019年11月初移栽，試驗(yàn)于2020年1月6日（苗期中期）進(jìn)行。

2.2 數(shù)據(jù)獲取

2.2.1 快速葉綠素?zé)晒鈹?shù)據(jù)采集

在每個(gè)氮素水平的3個(gè)小區(qū)內(nèi)，隨機(jī)選取3個(gè)長(zhǎng)勢(shì)相近且良好的油菜樣本，沿主莖自下而上的底1葉、底2葉、底3葉、底4葉、底5葉、底6葉和底7葉分別標(biāo)記為L(zhǎng)1、L2、L3、L4、L5、L6和L7，并根據(jù)所在垂直空間的相對(duì)位置劃分為基葉（L1、L2）、中葉（L3～L5）和頂葉（L6、L7）。試驗(yàn)時(shí)選取每個(gè)葉片葉脈右側(cè)中上部區(qū)域的點(diǎn)作為采樣點(diǎn)。在田間采用葉片夾（直徑4 mm）對(duì)葉片相應(yīng)的采樣點(diǎn)進(jìn)行20 min的暗適應(yīng)處理后，利用植物效率分析儀Handy-PEA（Handy Plant Efficiency Analyzer，Hansatech Instruments Ltd.，UK）進(jìn)行快速葉綠素?zé)晒庹T導(dǎo)動(dòng)力學(xué)曲線的測(cè)量。

2.2.2 生化參數(shù)測(cè)定

完成田間數(shù)據(jù)測(cè)量后，采集油菜各個(gè)葉位的葉片樣本，在實(shí)驗(yàn)室對(duì)每個(gè)葉位的采樣點(diǎn)用直徑0.85 cm的打孔器打孔稱量鮮重，并用1.8 mL 95%的乙醇在黑暗環(huán)境中浸泡24 h，用分光光度計(jì)（Epoch，BioTek Instruments，Winooski，USA）測(cè)定色素含量（葉綠素a、葉綠素b和類胡蘿卜素含量）。隨后將完成色素測(cè)定的樣本置于60 ℃的烘箱中烘干至恒重并稱量其干重，并用鮮重減去干重計(jì)算得到水分含量。

2.3 研究方法

2.3.1 OJIP曲線及JIP-測(cè)定參數(shù)

綠色植物或含有葉綠素的部分組織進(jìn)行一段時(shí)間的暗適應(yīng)處理，被突然暴露在可見(jiàn)光下便會(huì)發(fā)出一種強(qiáng)度不斷變化的暗紅色熒光信號(hào)，也稱為Kautsky效應(yīng)[20]。這期間，熒光信號(hào)先上升，后下降。將其剛暴露在光下時(shí)的最低熒光定義為O點(diǎn)，最高峰定義為P點(diǎn)，中間分別有J點(diǎn)（2 ms）和I點(diǎn)（30 ms）兩次階躍，快速葉綠素?zé)晒庹T導(dǎo)動(dòng)力學(xué)（OJIP）曲線指的就是熒光從O點(diǎn)到P點(diǎn)的變化過(guò)程[21，22]。

OJIP曲線形態(tài)的變化反映了植物葉片PSII的光合原初反應(yīng)電子傳遞鏈以及光合機(jī)構(gòu)的結(jié)構(gòu)和狀態(tài)的變化，通過(guò)JIP-測(cè)定的方法可以對(duì)該變化進(jìn)行定量分析和生理解釋。其理論基礎(chǔ)是生物膜能量流動(dòng)模型，天線色素將小部分吸收（Absorption，ABS）能轉(zhuǎn)化為熱和熒光，其余部分則被反應(yīng)中心（Reaction Center，RC）所捕獲（Trapping，TR），并將處于氧化態(tài)的初級(jí)電子受體QA激發(fā)還原為QA-，繼續(xù)往下進(jìn)行電子傳遞（Electron Transfer，ET），從而形成電子傳遞鏈[23]。以此過(guò)程為基礎(chǔ)推算出一系列葉綠素?zé)晒鈪?shù)的數(shù)據(jù)處理方法即為JIP-測(cè)定方法，具體的計(jì)算公式和參數(shù)含義如表1所示[21]。

2.3.2 主成分分析

主成分分析（Principal Component Analysis，PCA）以多變量數(shù)據(jù)為處理對(duì)象，旨在利用“降維”的思想，把多指標(biāo)轉(zhuǎn)化為少數(shù)幾個(gè)綜合指標(biāo)[24]?；贠JIP曲線的JIP-測(cè)定參數(shù)數(shù)量較多，且有些參數(shù)之間存在高度的相關(guān)性，從而增加了問(wèn)題分析的復(fù)雜性。因此利用PCA將數(shù)據(jù)從高維空間轉(zhuǎn)向低維空間，便于進(jìn)一步對(duì)不同氮素處理和不同葉位葉片熒光參數(shù)進(jìn)行統(tǒng)計(jì)分析。

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

對(duì)不同氮素處理水平下不同葉位的葉片生化參數(shù)，采用最小顯著差異法（Least Significant Difference，LSD）檢驗(yàn)顯著性差異（P < 0.05）。JIP-測(cè)定參數(shù)的PCA、生化參數(shù)和熒光參數(shù)之間的相關(guān)性分析以及顯著性分析均采用SPSS Statistics 19軟件（IBM Corporation，USA）。

3? 結(jié)果與分析

3.1 油菜冠層生化參數(shù)的垂直異質(zhì)性

油菜苗期中期各個(gè)葉位下葉片的葉綠素含量（Cab）、類胡蘿卜素含量（Cxc）、葉綠素與類胡蘿卜素含量的比值（Cab /Cxc）、葉綠素a和葉綠素b含量的比值（Chl a/b）、干物質(zhì)含量（Cm）以及水分含量（Cw）的變化結(jié)果如圖1所示。可以看到，上述生化參數(shù)值在不同葉位的葉片上均發(fā)生了顯著的變化，即垂直異質(zhì)性，且不同氮素處理水平下它們垂直空間上的變化基本一致。具體地，油菜葉片葉綠素（圖1（a））和類胡蘿卜素（圖1（b））含量從底部到頂部呈現(xiàn)開(kāi)口向下拋物線型的變化，干物質(zhì)（圖1（e））和水分（圖1（f））含量則呈現(xiàn)開(kāi)口向上拋物線型的變化，說(shuō)明在冠層垂直空間水平上中葉（L3～L5）的葉綠素、類胡蘿卜素含量最高，干物質(zhì)和水分含量則最低，而基葉（L1、L2）相較于頂葉（L6、L7）具有更高的葉綠素、水分含量和較低的類胡蘿卜素、干物質(zhì)含量。由于基葉一定程度上會(huì)受到遮蔽，其接受的光強(qiáng)相對(duì)頂葉較少，因此其所需要用來(lái)進(jìn)行光保護(hù)的類胡蘿卜素相對(duì)就少，且其干物質(zhì)需要向上部新葉轉(zhuǎn)運(yùn)而較低。特別的是，正常氮組（N1）的基葉葉綠素和類胡蘿卜素含量均高于缺氮組（N0）和過(guò)氮組（N2），而其中葉的葉綠素、類胡蘿卜素、干物質(zhì)、水分含量則均為最低。此外，頂葉的葉綠素含量隨著施氮量增加有下降的變化趨勢(shì)。這一結(jié)果表明氮素脅迫雖然對(duì)油菜冠層內(nèi)生化參數(shù)的垂直分布模式?jīng)]有顯著的影響，但在缺氮的情況下，油菜吸收的氮素會(huì)優(yōu)先分配給頂葉以滿足其生長(zhǎng)需求[25]，而氮素充足的情況下則會(huì)不斷供應(yīng)更多的頂部新葉生長(zhǎng)，從而出現(xiàn)了缺氮組頂葉葉綠素含量較高的情況。

葉綠素與類胡蘿卜素含量的比值（圖1（c））從底部到頂部呈現(xiàn)較為顯著的下降趨勢(shì)，進(jìn)一步反映了頂葉的葉綠素含量較少而類胡蘿卜素含量較高的特點(diǎn)，基葉則反之。葉綠素含量與類胡蘿卜素含量的比值有效反映了光保護(hù)葉黃素循環(huán)響應(yīng)[26，27]，說(shuō)明在吸收光能較多的頂葉中具有更活躍的葉黃素循環(huán)進(jìn)行光保護(hù)。此外，隨著施氮量的增加，比值也逐漸減小，這一結(jié)果也表明了葉綠素與類胡蘿卜素含量的比值與葉位（葉齡）、施氮水平均具有潛在的良好線性關(guān)系。葉綠素a/b反映了天線尺寸的大小，植物通過(guò)調(diào)整葉綠素a/b來(lái)適應(yīng)不同光強(qiáng)或氮素供應(yīng)的環(huán)境。在基葉中葉綠素a和葉綠素b含量的比值（圖1（d））顯著低于中葉和頂葉，且在各個(gè)葉位中基本都是隨著施氮量的增加而增加，這反映了基葉中的光合系統(tǒng)具有更大的天線尺寸（捕光色素葉綠素b所占的比例大）以適應(yīng)較低的外界光強(qiáng)[28]。

3.2 基于JIP-測(cè)定參數(shù)的油菜冠層光合生理垂直異質(zhì)性

以L1底葉的各個(gè)JIP-測(cè)定參數(shù)值為基準(zhǔn)，賦值為1，并將其他葉位的參數(shù)值與其分別做商，可有效觀察到油菜各個(gè)葉位葉片的JIP-測(cè)定參數(shù)差異，不同氮素處理水平下的結(jié)果如圖2所示。3個(gè)氮素水平中均出現(xiàn)較為顯著變化的參數(shù)是DFTotal、φRo、RE0 /RC、Sm和N，它們?cè)陧斎~中出現(xiàn)了極小值，說(shuō)明在頂葉中的電子鏈末端性能最弱。缺氮組（N0）和正常氮組（N1）中頂葉單位反應(yīng)中心或單位受光截面耗散的熱量（DI0 /RC和DI0 /CSm）均為最大，因?yàn)楣趯又懈鱾€(gè)葉位葉片所吸收的光能隨著葉位的升高是逐漸增加的，頂葉吸收較多的光能而使得其熱耗散也較高。但這個(gè)值在過(guò)量施氮組（N2）中卻并不是最大的，這可能與過(guò)量施氮組頂葉和中葉綠素a和葉綠素b含量的比值較高有關(guān)，天線尺寸相對(duì)較?。ú豆馍厝~綠素b較少）使其捕獲的光能并沒(méi)有比其他葉位的葉片多。

3.3 油菜葉片生化參數(shù)與JIP-測(cè)定參數(shù)的相關(guān)性分析

為進(jìn)一步研究油菜葉片生化參數(shù)與JIP-測(cè)定參數(shù)之間的潛在聯(lián)系，對(duì)兩者進(jìn)行皮爾遜相關(guān)性分析，所得結(jié)果如圖3所示。葉綠素和類胡蘿卜素比值的對(duì)數(shù)值與JIP-測(cè)定參數(shù)的相關(guān)性較比值本身有所提高。Cab、Cab /Cxc、log（Cab /Cxc）、Cm與多個(gè)JIP-參數(shù)之間存在較高的相關(guān)性，而Cxc、Chl a/b和Cw等與JIP-測(cè)定參數(shù)之間的相關(guān)性較低。其中，δRo、φRo和DFTotal與葉綠素含量呈顯著的高度正相關(guān)關(guān)系（r > 0.5，P < 0.01），DFTotal、φRo、PIABS、Sm和RE0 /RC與葉綠素和類胡蘿卜素的比值和其比值的對(duì)數(shù)均呈顯著的高度正相關(guān)關(guān)系，ABS/RC與干物質(zhì)含量呈高度正相關(guān)關(guān)系，而δRo和DFTotal則與其呈顯著的高度負(fù)相關(guān)關(guān)系（r < -0.5，P < 0.01）。這說(shuō)明一定程度上提高葉綠素含量和其與類胡蘿卜素的比值將顯著提高原初光化學(xué)反應(yīng)中電子傳遞到PSI傳遞鏈末端的量子產(chǎn)額和推動(dòng)力，而干物質(zhì)含量過(guò)高和水分含量過(guò)少則會(huì)阻礙電子在傳遞鏈中的傳遞過(guò)程，因?yàn)樽鳛樽畛蹼娮庸w的水分含量減少將使光合作用中水電解的過(guò)程受到影響。

從與葉綠素與類胡蘿卜比值的對(duì)數(shù)（圖4（a）～圖4（c））、葉綠素含量（圖4（d）～圖4（f））和干物質(zhì)（圖4（g）～圖4（i））相關(guān)性較高的3個(gè)JIP-測(cè)定參數(shù)進(jìn)行線性回歸的結(jié)果可以看到，3個(gè)生化參數(shù)均與DFTotal具有較顯著的線性關(guān)系（P < 0.0001），決定系數(shù)R2分別為0.45、0.35和0.27。結(jié)合不同葉位的散點(diǎn)分布情況可以看出，相比于葉綠素含量和干物質(zhì)含量，葉綠素含量與類胡蘿卜素含量比值的對(duì)數(shù)值能有效區(qū)分出基葉和頂葉，進(jìn)一步說(shuō)明了DFTotal參數(shù)具有反映光保護(hù)葉黃素循環(huán)等光合生理垂直異質(zhì)性的潛能。

3.4 基于PCA的JIP-測(cè)定參數(shù)聚類分析

通過(guò)油菜苗期中期葉片JIP-測(cè)定參數(shù)的主成分分析，進(jìn)一步探究在不同葉位和不同氮素處理下葉片光合機(jī)構(gòu)結(jié)構(gòu)和功能的變化規(guī)律，結(jié)果如圖5所示?？梢钥吹?，JIP-測(cè)定參數(shù)的分布基本形成了三個(gè)較好分離的簇。其中兩個(gè)（Cluster 1和Cluster 3）位于第一個(gè)主成分（PC1）上，一個(gè)（Cluster 2）位于第二個(gè)主成分（PC2）上。結(jié)合每一簇中JIP-測(cè)定參數(shù)的生理意義，可以發(fā)現(xiàn)它們分別對(duì)應(yīng)了不同的光合生理過(guò)程：光能吸收和捕獲階段（Cluster 1）、電子傳遞和熱耗散階段（Cluster 2）和電子鏈末端傳遞階段（Cluster 3）。進(jìn)一步可以推得，PC1可對(duì)應(yīng)PSI的活性，較高的值表示較高的PSI性能（高電子鏈末端傳輸效率），代表參數(shù)有φRo、δRo、RE0 /RC等;PC2則對(duì)應(yīng)PSⅠI活性，較高的值表示較高的PSⅠI性能（高光吸收、光化學(xué)和電子傳輸效率），代表參數(shù)包括ABS/CSm、φPo、ET0 /CSm等。

從不同氮素處理實(shí)驗(yàn)組的聚類結(jié)果可以發(fā)現(xiàn)，隨著施氮水平的提高，聚類的位置也在發(fā)生改變：缺氮組（N0）主要位于第三、四象限，正常施氮組（N1）主要分布在第二、三和四象限的Cluster1和Cluster3之間，而過(guò)量施氮組（N2）則主要分布于第一象限的Cluster2。結(jié)合JIP-測(cè)定參數(shù)分析可知，相較于正常施氮，過(guò)量施氮組油菜葉片具有更高的PSII和PSI性能，而缺氮組油菜葉片的PSII和PSI性能均較弱，特別是在光捕獲和能量傳遞階段;同一施氮水平下，隨著葉位的變化，聚類的位置也發(fā)生了變化：顏色較深的數(shù)據(jù)點(diǎn)所代表的底部葉片（L1）主要分布在第一、四象限，而顏色較淺的數(shù)據(jù)點(diǎn)所代表的頂部新葉（L7）則位于第二、三象限，表明油菜頂部新葉較底部成熟的葉片在PSI性能上有顯著的差異，即電子鏈末端傳輸效率比較低。

結(jié)合圖3、圖4中的結(jié)果，可以看到與葉片生化參數(shù)具有較高相關(guān)性的JIP-測(cè)定參數(shù)均分布在反映電子傳遞鏈末端性能的Cluster 3內(nèi)，其中DFTotal所在向量的方向正好也指向葉位升高的方向，這也與圖1（c）中的結(jié)果一致。DFTotal是電子能傳遞到電子鏈末端的推動(dòng)力，通過(guò)計(jì)算公式（30）可以看出該參數(shù)由PIABS、δRo兩個(gè)參數(shù)計(jì)算所得，這同時(shí)也解釋了PCA結(jié)果中這3個(gè)參數(shù)之間的夾角較小以及均與葉片生化參數(shù)有較高的相關(guān)性的原因。進(jìn)一步推測(cè)利用DFTotal這個(gè)綜合性能參數(shù)，可以有效地用于油菜冠層生化參數(shù)垂直異質(zhì)性的表征。

4? 結(jié) 論

本研究主要利用快速葉綠素?zé)晒饧夹g(shù)實(shí)現(xiàn)了對(duì)油菜冠層生化參數(shù)垂直異質(zhì)性的檢測(cè)與評(píng)估，得到以下研究結(jié)果。

（1）苗期中期油菜冠層的葉綠素含量、類胡蘿卜素含量、干物質(zhì)和水分含量均呈拋物線型的垂直分布，而葉綠素/類胡蘿卜素具有與其他生化參數(shù)不同的垂直分布模式，其隨著葉位升高和施氮量的增加逐漸下降，且與熒光參數(shù)DFTotal、φRo和RE0 /RC的垂直分布模式相同。

（2）熒光參數(shù)對(duì)油菜葉片生化參數(shù)中的葉綠素/類胡蘿卜、葉綠素和干物質(zhì)含量具有相對(duì)較強(qiáng)的評(píng)估能力。其中，電子能傳遞到電子鏈末端的推動(dòng)力性能參數(shù)DFTotal與葉綠素與類胡蘿卜素的比值、葉綠素含量均具有較高的正相關(guān)性（r > 0.5，P < 0.01），與干物質(zhì)含量則呈顯著的負(fù)相關(guān)（r < -0.5，P < 0.01）。

（3）缺氮一定程度上會(huì)降低苗期油菜葉片的PSII和PSI性能，通過(guò)φPo和RE0 /RC等熒光參數(shù)可以有效評(píng)估氮素脅迫。而不同葉位則在PSI性能即電子末端傳遞效率上具有顯著的差異，熒光參數(shù)DFTotal可以有效表征冠層生化參數(shù)的垂直異質(zhì)性。

上述結(jié)果表明，葉綠素?zé)晒饧夹g(shù)可作為一種評(píng)估作物生化參數(shù)垂直異質(zhì)性有效的替代方法，為實(shí)現(xiàn)基于高通量表型技術(shù)的作物生長(zhǎng)信息監(jiān)測(cè)提供了思路和技術(shù)支撐。擬進(jìn)一步對(duì)不同生長(zhǎng)期的油菜冠層內(nèi)生化參數(shù)進(jìn)行表征，實(shí)現(xiàn)油菜全生育期的生化參數(shù)垂直分布監(jiān)測(cè)。

參考文獻(xiàn)：

[1]??? ABDALLAH M， DUBOUSSET L， MEURIOT F， et al. Effect of mineral sulphur availability on nitrogen and sulphur uptake and remobilization during the vegetative growth of Brassica napus L[J]. Journal of Experimental Botany， 2010， 61（10）： 2635-2646.

[2]??? PAO Y， CHEN T， MOUALEU-NGANGUE D， et al. Environmental triggers for photosynthetic protein turnover determine the optimal nitrogen distribution and partitioning in the canopy[J]. Journal of Experimental Botany， 2019， 70（9）： 2419-2434.

[3]??? HIKOSAKA K. Optimality of nitrogen distribution among leaves in plant canopies[J]. Journal of Plant Research， 2016， 129（3）： 299-311.

[4]??? ADACHI S， YOSHIKAWA K， YAMANOUCHI U， et al. Fine mapping of carbon assimilation rate 8， a quantitative trait locus for flag leaf nitrogen content， stomatal conductance and photosynthesis in rice[J]. Frontiers in Plant Science， 2017， 8： ID 60.

[5]??? 王紀(jì)華， 王之杰， 黃文江， 等. 冬小麥冠層氮素的垂直分布及光譜響應(yīng)[J]. 遙感學(xué)報(bào)， 2004， 8（4）： 309-316.WANG J， WANG Z， HUANG W， et al. The vertical distribution characteristic and spectral response of canopy nitrogen in different layer of winter wheat[J]. Journal of Remote Sensing， 2004， 8（4）： 309-316.

[6]??? 葉曉青， 鄒勇， 余志虹， 等. 烤煙冠層光譜參數(shù)與氮素垂直分布相關(guān)性研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2013， 44（5）： 219-225.YE X， ZOU Y， YU Z， et al. Correlation between nitrogen vertical distribution and spectral characteristics of flue-cured tobacco[J]. Transactions of the CSAM， 2013， 44（5）： 219-225.

[7]??? FENG W， GUO B， WANG Z， et al. Measuring leaf nitrogen concentration in winter wheat using double-peak spectral reflection remote sensing data[J]. Field Crops Research， 2014， 159： 43-52.

[8]??? LI F， MISTELE B， HU Y， et al. Remotely estimating aerial N status of phenologically differing winter wheat cultivars grown in contrasting climatic and geographic zones in China and Germany[J]. Field Crops Research， 2012， 138： 21-32.

[9]??? STAGAKIS S， MARKOS N， SYKIOTI O， et al. Tracking seasonal changes of leaf and canopy light use efficiency in a Phlomis fruticosa Mediterranean ecosystem using field measurements and multi-angular satellite hyperspectral imagery[J]. ISPRS Journal of Photogrammetry and Remote Sensing， 2014， 97： 138-151.

[10]? ZHANG Q， CHEN J M， JU W， et al. Improving the ability of the photochemical reflectance index to track canopy light use efficiency through differentiating sunlit and shaded leaves[J]. Remote Sensing of Environment， 2017， 194：1-15.

[11]? STIRBET A， LAZ?R D， GUO Y， et al. Photosynthesis： Basics， history and modelling[J]. Annals of Botany， 2020， 126（4）： 511-537.

[12]? GUO Y， TAN J. Recent advances in the application of chlorophyll a fluorescence from photosystem II[J]. Photochemistry and Photobiology， 2015， 91（1）： 1-14.

[13]? GUANTER L， ZHANG Y， JUNG M， et al. Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence[J]. Proceedings of the National Academy of Sciences， 2014， 111（14）： 1327-1333.

[14]? PAPAGEORGIOU G. Chlorophyll a fluorescence： A signature of photosynthesis[M]. Dordrecht： Springer， 2004.

[15]? MAGYAR M， SIPKA G， KOV?CS L， et al.? Rate-limiting steps in the dark-to-light transition of Photosystem II-revealed by chlorophyll-a fluorescence induction[J]. Scientific Reports， 2018， 8（1）： 1-9.

[16]? SHRESTHA S， BRUECK H and ASCH F.?? Chlorophyll index， photochemical reflectance index and chlorophyll fluorescence measurements of rice leaves supplied with different N levels[J]. Journal of Photochemistry and Photobiology B： Biology， 2012， 113： 7-13.

[17]? LIU R， WANG Y， CHEN B， et al.? Effects of nitrogen levels on photosynthesis and chlorophyll fluorescence characteristics under drought stress in cotton flowering and boll-forming stage[J]. Acta Agronomica Sinica， 2008， 34（4）： 675-683.

[18]? FENG W， HE L， ZHANG H-Y， et al.? Assessment of plant nitrogen status using chlorophyll fluorescence parameters of the upper leaves in winter wheat[J]. European Journal of Agronomy， 2015， 64： 78-87.

[19]? LARBI A， V?ZQUEZ S， EL-JENDOUBI H， et al. Canopy light heterogeneity drives leaf anatomical， eco-physiological， and photosynthetic changes in olive trees grown in a high-density plantation[J]. Photosynthesis Research， 2015， 123（2）： 141-155.

[20]? KAUTSKY H， HIRSCH A. Neue versuche zur kohlens?ureassimilation[J]. Naturwissenschaften， 1931， 19（48）： 964-964.

[21]? STRASSER R J， TSIMILLI-MICHAEL M， SRIVASTAVA A. Analysis of the chlorophyll a fluorescence transient[M]. Dordrecht： Springer， 2004.

[22]? STIRBET A. On the relation between the Kautsky effect （chlorophyll a fluorescence induction） and photosystem II： Basics and applications of the OJIP fluorescence transient[J]. Journal of Photochemistry and Photobiology B： Biology， 2011， 104（1-2）： 236-257.

[23]? TSIMILLI-MICHAEL M， EGGENBERG P， BIRO B， et al. Synergistic and antagonistic effects of arbuscular mycorrhizal fungi and Azospirillum and Rhizobium nitrogen-fixers on the photosynthetic activity of alfalfa， probed by the polyphasic chlorophyll a fluorescence transient OJIP[J]. Applied Soil Ecology， 2000， 15（2）： 169-182.

[24]? ABDI H， WILLIAMS L J. Principal component analysis[J]. Wiley Interdisciplinary Reviews： Computational Statistics， 2010， 2（4）： 433-459.

[25]? HUANG W， YANG Q， PU R， et al. Estimation of nitrogen vertical distribution by bi-directional canopy reflectance in winter wheat[J]. Sensors， 2014， 14（11）： 20347-20359.

[26]? MERZLYAK M N， GITELSON A A， CHIVKUNOVA O B， et al. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening[J]. Physiologia Plantarum， 1999， 106（1）： 135-141.

[27]? ZHOU X， HUANG W， ZHANG J， et al. A novel combined spectral index for estimating the ratio of carotenoid to chlorophyll content to monitor crop physiological and phenological status[J]. International Journal of Applied Earth Observation and Geoinformation， 2019， 76： 128-142.

[28]? DIN? E， CEPPI M G， T?TH S Z， et al. The chl a fluorescence intensity is remarkably insensitive to changes in the chlorophyll content of the leaf as long as the chl a/b ratio remains unaffected[J]. Biochimica et Biophysica Acta （BBA）-Bioenergetics， 2012， 1817（5）： 770-779.

Vertical Heterogeneity Analysis of Biochemical Parameters in Oilseed Rape Canopy Based on Fast Chlorophyll Fluorescence Technology

ZHANG Jiafei1，2， WAN Liang1，2， HE Yong1，2，3， CEN Haiyan1，2，3*

（1.College of Biosystems Engineering and Food Science， Zhejiang University， Hangzhou 310058， China;

2.Key Laboratory of Spectroscopy Sensing， Ministry of Agriculture and Rural Affairs， Hangzhou 310058， China;

3.State Key Laboratory of Modern Spectroscopic Instruments， Zhejiang University， Hangzhou 310027， China）

Abstract： Accurate acquisition of crop canopy biochemical information is of great significance for monitoring crop growth and guiding precise fertilization. Previous vertical distribution researches of crop biochemical information were mainly based on hyperspectral inversion， which was lack of the association of plant photosynthesis physiology. This study mainly investigated the vertical distribution characteristics of biochemical parameters such as chlorophyll， carotenoid， dry matter， and water content in the oilseed rape canopy under different nitrogen treatments at the mid-seedling stage. The photosynthetic performance of leaves was measured by using fast chlorophyll fluorescence technology， and linear regression and principal component analysis were further implemented to explore the internal relationship between fluorescence response and biochemical parameters. The results showed that： （1） The chlorophyll content， carotenoid content， dry matter and water content of the rape canopy at the mid-seedling stage all showed a parabolic vertical distribution， while the ratio of chlorophyll to carotenoids content gradually decreases with the leaf position and nitrogen treatments， which was the same as the vertical distribution pattern of fluorescence parameters such as driving force comprehensive performance （DFTotal） and end electron chain quantum yield （φRo） and other fluorescence parameters could be used to diagnose nitrogen stress; （2） JIP-test parameters， especially DFTotal， had a good performance to evaluate the chlorophyll/carotenoids， chlorophyll and dry matter content of oilseed rape leaves; （3） Nitrogen deficiency would weaken the PSII and PSI performance of oilseed rape leaves at the mid-seedling stage， and the maximum photochemical efficiency （φPo） could be used to diagnose nitrogen stress. There was a significant difference in the PSI performance， namely electron transfer efficiency at the end acceptors of leaves in the different leaf position， hence the comprehensive performance parameter DFTotal could be an effective characterization of the vertical heterogeneity of canopy biochemical parameters. These findings indicated the feasibility of applying the rapid chlorophyll fluorescence technology to crop biochemical information heterogeneity monitoring and provided new ideas and technical support for guiding precise fertilization and achieving high-quality and high-yield.

Key words： fast chlorophyll fluorescence transient; JIP-test; vertical heterogeneity; chlorophyll; carotenoid; nitrogen