油茶葉片氮磷鉀含量與經(jīng)濟(jì)性狀的關(guān)聯(lián)分析

曹永慶，任華東，王開(kāi)良，陳新建，汪舍平，俞春蓮，姚小華＊

(1. 中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所，浙江 杭州 311400；2. 浙江省常山縣林業(yè)局，浙江 常山 324200)

油茶（Camellia oleiferaAbel.）是我國(guó)南方重要的木本油料樹(shù)種，具有綜合利用價(jià)值高、生態(tài)效益好的優(yōu)點(diǎn)[1-2]，目前全國(guó)種植面積已達(dá)420 多萬(wàn)hm2，在脫貧攻堅(jiān)、鄉(xiāng)村振興等國(guó)家戰(zhàn)略中發(fā)揮著重要作用。油茶適生分布區(qū)多為pH 值4.5～6.5 的紅壤、黃壤或黃棕壤，在固定、淋溶等作用下，土壤礦質(zhì)養(yǎng)分虧缺的現(xiàn)象普遍存在[3-5]，開(kāi)展?fàn)I養(yǎng)管理研究對(duì)推動(dòng)產(chǎn)業(yè)提質(zhì)增效具有重要意義。

前期已對(duì)油茶主要礦質(zhì)的需求和積累特征進(jìn)行了分析[6-7]，并在配方施肥方面作了大量的試驗(yàn)研究[8-13]，提出了不同的施肥策略，但要實(shí)現(xiàn)科學(xué)精準(zhǔn)營(yíng)養(yǎng)管理，還應(yīng)根據(jù)葉片具體營(yíng)養(yǎng)狀態(tài)進(jìn)行配方施肥。袁軍[14]、周裕新[15]基于綜合營(yíng)養(yǎng)診斷施肥法（DRIS）或矢量分析法提出了高產(chǎn)林分的主要礦質(zhì)營(yíng)養(yǎng)參考標(biāo)準(zhǔn)；然而，不同區(qū)域油茶林地土壤養(yǎng)分狀態(tài)存在一定的差異[3]，且不同品系的營(yíng)養(yǎng)需求特征也不同[16]，這給營(yíng)養(yǎng)評(píng)估診斷和施肥管理帶來(lái)極大的不確定性。需要探索建立不依賴(lài)區(qū)域立地條件和品種、既穩(wěn)定又可靠的營(yíng)養(yǎng)診斷和評(píng)估技術(shù)。

基于離子平衡的等角對(duì)數(shù)比（ILR）數(shù)據(jù)處理分析方法，充分考慮了礦質(zhì)離子間的相互作用關(guān)系，被認(rèn)為是目前最適用于植物離子組學(xué)數(shù)據(jù)的描述方法[17]，與基于單比率的DRIS 分析方法相比，在數(shù)據(jù)去冗降噪等方面具有顯著優(yōu)勢(shì)，如Parent等[18]利用ILR 離子組特征分析方法，對(duì)175 個(gè)芒果（Mangifera indicaL.）園的營(yíng)養(yǎng)狀態(tài)進(jìn)行評(píng)估，并提出導(dǎo)致其低產(chǎn)的潛在影響因子。本文以廣泛栽培的油茶‘長(zhǎng)林4 號(hào)’和‘長(zhǎng)林53 號(hào)’為研究對(duì)象，利用ILR 分析方法對(duì)礦質(zhì)離子數(shù)據(jù)進(jìn)行處理，并研究分析油茶葉片氮、磷、鉀礦質(zhì)元素含量與產(chǎn)量、種仁含油率等經(jīng)濟(jì)性狀的相關(guān)性，以期為油茶樹(shù)體營(yíng)養(yǎng)狀態(tài)科學(xué)評(píng)估和營(yíng)養(yǎng)管理提供參考。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)地位于浙江省常山縣東案鄉(xiāng)油茶種植基地，屬低山丘陵區(qū)，具有典型亞熱帶季風(fēng)氣候，四季分明，光照充足，年平均氣溫17.4℃，年平均降水量1 725 mm，年平均日照時(shí)數(shù)1 975 h，年平均有效積溫5 514℃，年平均無(wú)霜期238 d；土壤以紅黃壤為主，土層深厚，有機(jī)質(zhì)含量13.9 g·kg?1，水解性氮75.4 mg·kg?1，有效磷10.9 mg·kg?1，速效鉀110 mg·kg?1?；刂饕筒枇挤N為 ‘長(zhǎng)林4 號(hào)’、‘長(zhǎng)林40 號(hào)’、‘長(zhǎng)林53 號(hào)’、‘長(zhǎng)林18 號(hào)’、‘長(zhǎng)林3 號(hào)’、‘長(zhǎng)林55 號(hào)’等“長(zhǎng)林”系列國(guó)家審定油茶良種，2013 年種植，株行距2 m × 3 m，正常管理。

1.2 試驗(yàn)方法

2018 年10 月果實(shí)成熟時(shí)，隨機(jī)選取‘長(zhǎng)林4 號(hào)’和‘長(zhǎng)林53 號(hào)’植株各30 株，分單株進(jìn)行取樣和測(cè)定。由東西、南北兩處測(cè)量樹(shù)冠寬度，記錄單株鮮果產(chǎn)量，隨機(jī)抽取20 個(gè)果實(shí)，進(jìn)行種仁含油率分析，并在樹(shù)冠中部東、西、南、北4 個(gè)方向，采取結(jié)果枝組生長(zhǎng)正常的當(dāng)年生春梢中部葉片50 片，進(jìn)行氮、磷、鉀礦質(zhì)元素含量分析；記錄葉片數(shù)和花芽數(shù)，并于2019 年1 月跟蹤調(diào)查座果率，每株調(diào)查20 個(gè)枝條。

1.3 測(cè)定方法

將葉片樣品剪碎，100～105℃下殺酶15 min，然后70～80℃下烘干至恒質(zhì)量并粉碎，用濃H2SO4-H2O2消解，以5020 型流動(dòng)注射分析儀（瑞典Tecafor 公司）測(cè)定氮元素含量，鉬銻抗比色法測(cè)定磷元素含量，火焰分光光度計(jì)法測(cè)定鉀元素含量。

采用正己烷索氏回流方法測(cè)定種仁含油率。將種子于60℃烘干，剝?nèi)シN皮，將種仁磨碎至粉末，稱(chēng)取1～2 g 油茶種仁粉末，用烘干的濾紙包好，105℃徹底烘干水分（多次稱(chēng)量，直至質(zhì)量不變，一般需2 h），最后上索氏回流管，回流6 h。用油脂占種仁總質(zhì)量的百分比計(jì)算種仁含油率。

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

按照《油茶主要性狀調(diào)查測(cè)定規(guī)范LY/T2955-2018》[19]計(jì)算冠幅面積，由產(chǎn)量和冠幅面積的比值計(jì)算單位面積冠幅產(chǎn)量；以座果數(shù)與花芽數(shù)的比值計(jì)算座果率；采用K 均值聚類(lèi)分析的方法，對(duì)產(chǎn)量、種仁含油率、花芽/葉指標(biāo)進(jìn)行分組。

參照Parent 等[20]的方法，通過(guò)ILR[N,P|K]=計(jì)算ILR[N,P|K]值，其中，g(CNCp)為N 和P元素含量的幾何平均值，g(CK)表示K 元素含量。

數(shù)據(jù)采用SPSS18.0 和Microsoft Excel 2003 軟件進(jìn)行統(tǒng)計(jì)分析并作圖。

2 結(jié)果與分析

2.1 油茶葉片氮、磷、鉀元素含量的差異和相關(guān)性分析

油茶葉片氮元素含量最高，其次為鉀元素，磷元素含量最低，且不同油茶品種的葉片氮、磷、鉀元素含量差異顯著（圖1）?！L(zhǎng)林53 號(hào)’葉片氮、磷、鉀礦質(zhì)元素的平均含量分別為12.72、0.89、3.65 g·kg?1，顯著高于‘長(zhǎng)林4 號(hào)’葉片氮、磷、鉀的含量，后者3 種元素含量分別為11.40、0.80、3.11 g·kg?1。

圖1 ‘長(zhǎng)林4 號(hào)’和‘長(zhǎng)林53 號(hào)’油茶葉片氮、磷、鉀元素的含量Fig.1 Nitrogen,phosphorus and potassium content in leaves of ‘Changlin 4’ and ‘Changlin 53’

油茶葉片氮、磷、鉀元素含量的相關(guān)性分析（圖2）表明：3 種礦質(zhì)元素含量間呈顯著正相關(guān)（P< 0.01），氮元素和磷元素、磷元素和鉀元素、氮元素和鉀元素含量的相關(guān)系數(shù)分別為0.55、0.52、0.60，單一元素吸收量的變化會(huì)影響其它2 個(gè)元素的吸收和含量。

圖2 油茶葉片氮、磷、鉀元素含量的相關(guān)性Fig.2 Correlations among nitrogen,phosphorus and potassium content in leaves of oil-tea camellia

2.2 油茶葉片氮、磷、鉀元素含量與產(chǎn)量的相關(guān)性分析

由圖3 可以看出，‘長(zhǎng)林4 號(hào)’單位面積冠幅產(chǎn)量與葉片氮元素含量和ILR[N,P|K]值呈顯著正相關(guān)（P< 0.05），R2值分別為0.470 3 和0.405 6，其相關(guān)性高于磷元素（R2= 0.133 4）和與鉀元素（R2= 0.001 5）。根據(jù)擬合曲線函數(shù)計(jì)算，當(dāng)葉片氮、磷元素含量以及ILR[N,P|K]值分別大于11.01、0.77 g·kg?1和0.40 時(shí)，‘長(zhǎng)林4 號(hào)’單位面積冠幅產(chǎn)量將達(dá)1.50 kg 以上。

圖3 ‘長(zhǎng)林4 號(hào)’葉片氮、磷、鉀元素含量與產(chǎn)量的相關(guān)性Fig.3 Correlations between nitrogen,phosphorus,potassium content and yield of oil-tea camellia ‘Changlin 4’

與‘長(zhǎng)林4 號(hào)’類(lèi)似，‘長(zhǎng)林53 號(hào)’單位面積冠幅產(chǎn)量與葉片氮、磷元素含量以及ILR[N,P|K]值也呈顯著正相關(guān)（P< 0.05），R2值分別為0.417 1、0.679 3、0.713 6，與鉀元素（R2= 0.026 2）含量無(wú)相關(guān)性（圖4）。不同的是，‘長(zhǎng)林53 號(hào)’單位面積冠幅產(chǎn)量與磷元素含量的相關(guān)性高于氮元素。當(dāng)葉片氮、磷元素含量以及ILR[N,P|K]值分別大于13.02、0.89 g·kg?1和0.40 時(shí)，‘長(zhǎng)林53 號(hào)’單位面積冠幅產(chǎn)量將達(dá)1.50 kg 以上。

圖4 ‘長(zhǎng)林53 號(hào)’葉片氮、磷、鉀元素含量與產(chǎn)量的相關(guān)性Fig.4 Correlations between nitrogen,phosphorus,potassium content and yield of oil-tea camellia ‘Changlin 53’

2.3 油茶葉片氮、磷、鉀元素含量與種仁含油率的相關(guān)性分析

圖5 表明：‘長(zhǎng)林4 號(hào)’果實(shí)種仁含油率與葉片氮、磷、鉀元素含量無(wú)顯著相關(guān)性（R2< 0.1），而與ILR[N,P|K]值呈顯著負(fù)相關(guān)（R2= 0.807 4，P<0.01）。根據(jù)擬合曲線函數(shù)計(jì)算，當(dāng)ILR[N,P|K]值小于0.45 時(shí)，‘長(zhǎng)林4 號(hào)’果實(shí)種仁含油率達(dá)平均值（45%）以上。

與‘長(zhǎng)林4 號(hào)’不同，‘長(zhǎng)林53 號(hào)’果實(shí)種仁含油率與葉片氮、磷元素含量以及ILR[N,P|K]值呈顯著負(fù)相關(guān)關(guān)系（圖6），且與磷元素（R2= 0.774 3）和ILR[N,P|K]值（R2= 0.857 3）的相關(guān)性高于氮元素（R2= 0.4663），與鉀元素（R2= 0.002 5）含量則無(wú)相關(guān)性。當(dāng)葉片氮、磷元素含量以及ILR[N,P|K]值分別低于13.36、0.90 g·kg?1和0.43 時(shí)，‘長(zhǎng)林53 號(hào)’果實(shí)種仁含油率達(dá)平均值（40%）以上。

2.4 油茶葉片ILR 值與花芽的相關(guān)性分析

圖7 表明：‘長(zhǎng)林4 號(hào)’、‘長(zhǎng)林53 號(hào)’油茶春梢花芽數(shù)和葉片數(shù)的比值與ILR[N,P|K]值均呈顯著正相關(guān)（P< 0.01），相關(guān)系數(shù)分別為0.943 7 和0.988 7。此外，由花芽量和座果率的關(guān)系（圖8）可知：隨著花芽/葉的增加，座果率表現(xiàn)出先升高后下降的趨勢(shì)，當(dāng)0.29 < 花芽/葉 < 0.96，即‘長(zhǎng)林4 號(hào)’0.39

3 討論

不同油茶品種因其適應(yīng)性的差異，植株體內(nèi)礦質(zhì)元素的含量也不同[21]，如具有磷高效吸收特性的‘長(zhǎng)林166 號(hào)’，其葉片磷元素的含量顯著高于‘長(zhǎng)林4 號(hào)’[22]。本研究中，‘長(zhǎng)林53 號(hào)’油茶葉片的氮、磷、鉀元素含量顯著高于‘長(zhǎng)林4 號(hào)’，這與品種特性有關(guān)。此外，氮、磷、鉀3 種礦質(zhì)元素含量間存在明顯的正相關(guān)關(guān)系，這也與作者前期的研究結(jié)果相一致[23]。

圖5 ‘長(zhǎng)林4 號(hào)’葉片氮、磷、鉀元素含量與種仁含油率的相關(guān)性Fig.5 Correlations between nitrogen,phosphorus,potassium content and oil content in kernel of oil-tea camellia ‘Changlin 4’

圖6 ‘長(zhǎng)林53 號(hào)’葉片氮、磷、鉀元素含量與種仁含油率的相關(guān)性Fig.6 Correlations between nitrogen,phosphorus,potassium content and oil content in kernel of oil-tea camellia ‘Changlin 53’

圖7 油茶葉片氮、磷、鉀元素平衡對(duì)花芽形成的影響Fig.7 Effect of nitrogen,phosphorus and potassium balance in leaves on flower bud formation

圖8 花芽/葉比值和座果率的關(guān)系Fig.8 Relations between flower bud/leaves ratio and fruit set ratio

鑒于3 種礦質(zhì)元素間相互平衡關(guān)系，利用等角對(duì)數(shù)比（ILR）轉(zhuǎn)換的數(shù)據(jù)處理分析方法對(duì)數(shù)據(jù)進(jìn)行去冗降噪，構(gòu)建ILR[N,P|K]參數(shù)[20]并進(jìn)一步分析了其與產(chǎn)量、種仁含油率等經(jīng)濟(jì)性狀的關(guān)系，發(fā)現(xiàn)ILR[N,P|K]值與產(chǎn)量和花芽/葉比值呈顯著正相關(guān)，與種仁含油率呈顯著負(fù)相關(guān)；隨著ILR[N,P|K]值增大，單位冠幅面積產(chǎn)量和花芽數(shù)量表現(xiàn)出升高的趨勢(shì)，種仁含油率則呈下降趨勢(shì)?？梢?jiàn)葉片氮、磷、鉀元素含量不僅能影響油茶產(chǎn)量和種仁油脂的合成積累[11-13]，還能顯著影響花芽的分化形成[24]。利用ILR[N,P|K]在評(píng)估單位面積冠幅產(chǎn)量和種仁含油率時(shí)，相較單一礦質(zhì)元素含量指標(biāo)更穩(wěn)定可靠，如以單位冠幅產(chǎn)量高于1.50 kg 為標(biāo)準(zhǔn)時(shí)，‘長(zhǎng)林4 號(hào)’和‘長(zhǎng)林53 號(hào)’ILR[N,P|K]值均為0.40 以上，而對(duì)應(yīng)的葉片氮元素含量閾值則差異較大，分別為11.01 g·kg?1和13.02 g·kg?1以上；在種仁含油率方面，以長(zhǎng)林4 號(hào)’種仁含油率45%、‘長(zhǎng)林53 號(hào)’種仁含油率40%為參考，ILR[N,P|K]值則宜控制在0.43～0.45 以下。

根據(jù)國(guó)家現(xiàn)行的油茶豐產(chǎn)栽培技術(shù)規(guī)程，一般要求5 個(gè)品種或以上進(jìn)行混合栽植，這就導(dǎo)致單以葉片礦質(zhì)元素含量為參考開(kāi)展田間營(yíng)養(yǎng)評(píng)估和管理時(shí)，會(huì)因品種營(yíng)養(yǎng)利用差異而導(dǎo)致評(píng)判標(biāo)準(zhǔn)難統(tǒng)一、不準(zhǔn)確的情況，如對(duì)‘長(zhǎng)林4 號(hào)’而言，葉片氮元素含量11.40 g·kg?1屬于正常營(yíng)養(yǎng)狀態(tài)，而對(duì)‘長(zhǎng)林53 號(hào)’而言，則屬于氮虧缺的營(yíng)養(yǎng)狀態(tài)。鑒于ILR[N,P|K]在評(píng)估單位面積冠幅產(chǎn)量和種仁含油率時(shí)品種間具有較好的穩(wěn)定性和可靠性，這為解決油茶營(yíng)養(yǎng)科學(xué)評(píng)估問(wèn)題提供了新的方法。以本研究為例，在充分平衡產(chǎn)量（1.50 kg）、種仁含油率和開(kāi)花數(shù)量指標(biāo)，以保障花芽和座果率的平衡，促進(jìn)高產(chǎn)穩(wěn)產(chǎn)為目標(biāo)，‘長(zhǎng)林4 號(hào)’和‘長(zhǎng)林53 號(hào)’混栽油茶林分的葉片ILR[N,P|K]值宜控制在0.40～0.43 之間，當(dāng)ILR[N,P|K]值低于0.40 時(shí)，應(yīng)適時(shí)補(bǔ)充氮磷肥，ILR[N,P|K]值高于0.43 時(shí)，適時(shí)補(bǔ)充鉀肥。

本研究?jī)H基于同一立地條件和管理水平，探索了礦質(zhì)元素含量多元數(shù)據(jù)等角對(duì)數(shù)比處理方法在油茶礦質(zhì)營(yíng)養(yǎng)分析及營(yíng)養(yǎng)評(píng)估應(yīng)用中的可行性，旨在為油茶營(yíng)養(yǎng)狀態(tài)的診斷和評(píng)估探索一種新的思路，在相關(guān)研究中未涉及其它的栽培品種（如長(zhǎng)林40 號(hào)、長(zhǎng)林18 號(hào)等）以及極端條件下的個(gè)體（如負(fù)載過(guò)度、營(yíng)養(yǎng)極度失衡的植株等），其穩(wěn)定性和普適性如何，仍有待于進(jìn)一步研究。此外，我國(guó)油茶栽培品種多，分布范圍廣，立地土壤類(lèi)型多樣復(fù)雜，不同栽培區(qū)劃的同一品種、同一栽培區(qū)劃內(nèi)的不同品種在生物學(xué)性狀上均表現(xiàn)出較大的差異[25]，不同土壤類(lèi)型也影響油茶樹(shù)體礦質(zhì)元素的含量[26]，在今后的研究中，有必要進(jìn)一步研究解析葉片礦質(zhì)離子組與土壤類(lèi)型、基因型以及重要性狀的關(guān)聯(lián)特征，為建立系統(tǒng)、科學(xué)、有效的油茶營(yíng)養(yǎng)狀態(tài)評(píng)估技術(shù)奠定基礎(chǔ)。

4 結(jié)論

不同油茶品種的葉片氮、磷、鉀礦質(zhì)元素含量差異顯著且存在顯著正相關(guān)關(guān)系，以ILR[N,P|K]值為參考指標(biāo)在評(píng)估不同品種的單位面積冠幅產(chǎn)量和種仁含油率時(shí)較單一元素指標(biāo)更穩(wěn)定可靠，油茶林分的葉片ILR[N,P|K]值宜控制在0.40～0.43 之間。本研究為油茶營(yíng)養(yǎng)狀態(tài)的診斷和評(píng)估提供了一種新的思路和方法。