氮磷葉面噴施對云南松苗木葉綠素含量及其異速生長關系的影響

黃 鍵王德新楊 松鄭超凡陸莊躍許玉蘭

（1. 云南省林木種苗工作總站，云南 昆明 650215；2. 西南林業大學林學院，云南 昆明 650233；3. 西南林業大學云南生物多樣性研究院，云南 昆明 650233）

光合作用是植株進行的重要生理活動，葉綠素是參與光合作用非常重要的一種光合色素[1]，葉綠素含量與葉片光合速率密切相關[2]，葉綠素含量可以反映植株的光合作用能力及營養健康狀態[3]，成為衡量光合產物的重要指標[4]。葉綠素是一個復雜的數量性狀，極易受外界環境，比如光、溫、水、肥等因素的影響[5-7]。研究表明，增施氮（N）對葉綠素合成和光系統功能有明顯的促進作用[8-10]。但不同用量下的響應不一致，外源性肥料添加促進作用與用量、植物自身元素虧缺狀態有關[9-11]。N、磷（P）作為植物的基本營養元素，在植物生長和各種生理調節機制中發揮著重要作用[12-15]，其中控制好N、P 配比尤為重要。尹意婷等[12]在N、P 配施對青岡（Quercus glauca）、苦櫧（Castanopsis sclerophylla）幼苗生物量分配及葉片養分含量的影響試驗發現，少量的N 和P 對苗木生物量的積累影響不大，而合理的施N 和施P 能夠顯著的增加苗木生物量。因此合理施肥對促進苗木正常生長，提高生產力、改善植物體內元素含量具有十分重要的作用。施肥的方式也較多，其中葉面施肥相對于傳統施肥方式具有利用效率高、用量少、養分吸收快、便于噴施等優勢，目前葉面施肥的研究越來越廣泛[16-17] 。

云南松（Pinus yunnanensis）為松科松屬的常綠高大針葉喬木，具有喜光、耐干旱，而且對土壤條件的要求較低，不僅能在酸性紅壤、紅黃壤和棕色森林土或微石灰性土壤上生長，也能在其他樹種不能生長的貧瘠地或水土流失嚴重的荒山生長[18]。云南松苗期生長緩慢，對其生長促成培育非常有必要，其中施肥是較為常見簡易的途徑，前期的研究表明，N、P 添加可促進云南松苗木的生長[17,19]，但不同的N、P 用量如何影響植株的葉綠素含量不清。本研究闡述了不同N、P 葉面噴施下葉綠素a、葉綠素b 指標的變化特征，以探索N、P 添加對云南松苗木葉綠素含量、比值及其異速生長關系的影響，進而篩選云南松生長的適宜N、P 噴施量，以期為云南松苗木的促成培育提供基礎數據，為云南松苗木合理施肥和高產高效栽培提供理論基礎。

1 材料與方法

1.1 研究材料

試驗地設置在西南林業大學的苗圃內，地處云南省昆明市，位于北緯25°04'00″、東經102°45'41″，海拔約1 945 m，屬北亞熱帶半濕潤高原季風氣候。本研究以2 年生云南松幼苗為材料，施肥前選取株長勢基本一致，健康無病蟲害的幼苗組成不同的試驗小組，插牌標記并進行N、P 葉面噴施試驗。試驗用土以西南林業大學校內紅壤土與育苗基質公司的腐殖土按體積比為1∶2 混合均勻作為育苗基質。

1.2 研究方法

試驗采用N、P 各2 因素3 水平試驗設計，外源N、P 按高、中、低濃度設置3 個梯度[19]，N采用尿素（CH4N2O）、P 采用磷酸二氫鉀（KH2PO4），N水平N1、N2、N3分別為0%、0.25%、0.5%（以N 含量計算，質量分數），P 肥水平P1、P2、P3分別為0%、0.2%、0.4%（以P 含量計算，質量分數），具體方案如表1 所示。試驗共9 個處理，每個處理40 株，重復3 次，合計1 080 株云南松苗木。試驗于2020 年6 月底至8 月底采用等量施肥法每隔7 d 噴施1 次，共噴施8 次，噴施周期為2 個月，噴施時盡可能使每個處理每個單株的施肥量相近。并對云南松苗木進行定期澆水、除草，松土。

表 1 云南松苗木N、P 葉面噴施用量Table 1 Experimental design of foliar spraying of nitrogen and phosphorus on P. yunnanensisseedlings

1.3 數據測定

本試驗數據測量時間為2020 年8 月、10 月、12 月底各測量1 次，共測3 次，各處理各個重復內隨機選8 株長勢一致的苗木，用混合采樣法采摘當年生成熟針葉，用丙酮提取法對葉綠素進行提取，用分光光度法對葉綠素含量進行測定，重復3 次。然后計算葉綠素a 含量、葉綠素b 含量、葉綠素（a + b）含量以及葉綠素a/b。

1.4 數據分析

對所獲得的數據進行統計整理，使用Excel 2016、SPSS 等軟件對其N、P 以及交互作用對葉綠素影響進行雙因素方差分析，利用Duncan 多重比較（α=0.05）檢驗各處理間的差異顯著性[20-21]，用平均值 ± 標準誤差進行繪圖[22]。利用SPSS 構建葉綠素指標與N、P 用量間的二元二次回歸方程，進一步采用降維法對N、P 配施效果進行單因素效應分析，確定最佳N、P 用量及其最高理論產量值和葉綠素含量的適宜N、P 濃度配比[19]。

2 結果與分析

2.1 N、P 葉面噴施對云南松苗木葉綠素含量的變異來源分析

在測定的8 月、10 月、12 月，N、P 交互作用對云南松苗木的葉綠素a、葉綠素b、葉綠素（a + b）含量的影響達極顯著水平（P＜0.01），而8 月和12 月N、P 交互作用對葉綠素a/b 無顯著相互作用。從N 葉面噴施單個因素來看，對8 月、10 月、12月葉綠素a 含量的影響達極顯著水平、對8 月葉綠素b 含量的影響達極顯著水平、對8 月和10 月葉綠素（a + b）含量的影響達極顯著水平、對10 月和12 月葉綠素a/b 的影響達極顯著水平（P＜0.01），而對8 月和12 月葉綠素b、12 月葉綠素（a + b）含量均無顯著影響。從P 葉面噴施單個因素來看，對各月的葉綠素b 含量均無顯著影響，對葉綠素a 含量影響在8 月和10 月達到極顯著水平（P＜0.01），在12月達顯著水平（P＜0.05）?？傮w來看，N、P 葉面噴施對云南松苗木葉綠素含量影響明顯，且隨著時間推移發生改變。

2.2 N、P 葉面噴施對云南松苗木葉綠素含量的影響

2.2.1 N、P 葉面噴施對云南松苗木葉綠素a 含量的影響

由圖1 可知，N、P 葉面噴施對云南松苗木葉綠素a 含量有不同程度的影響。8 月，處理3、4、5 的葉綠素a 含量均高于處理1，處理3 最高，處理9 最低。10 月，所有施肥處理葉綠素a 含量均大于處理1，處理9 中最大為13.407 mg/g，處理1 中最小。12 月，處理4、7、8、9 的葉綠素a 含量均大于處理1，處理4 最大，處理6 最?。?.122 mg/g）。綜合來看，在測定時期內，處理4 葉綠素a 含量均大于處理1，即處理4（中N低P）對云南松苗木葉綠素a 含量產生明顯的促進作用。

圖 1 各處理間葉綠素a 含量的動態變化圖Fig. 1 The dynamic change of chlorophyll a content among different treatments

2.2.2 N、P 葉面噴施對云南松苗木葉綠素b 含量的影響

由圖2 可知，N、P 葉面噴施對云南松苗木葉綠素b 的含量也有一定的影響。其中，8 月各處理的葉綠素b 波動于3.220～4.600 mg/g，以處理5 最高、處理6 最低，處理2、6、8、9 的葉綠素b 含量均小于處理1。10 月，所有施肥處理中葉綠素b 含量均大于處理1，且在處理2 中觀測到最大值5.237 mg/g，處理1 中觀測到最小值3.559 mg/g。12 月，處理2、6 的葉綠素b 含量小于處理1，以處理4 最高、處理6 最低。綜合來看，在整個試驗周期中，處理3、4、5、7 的葉綠素b 含量均大于處理1，以處理4 最大，即處理4（中N低P）對云南松苗木葉綠素b 含量促進作用明顯。

圖 2 各處理間葉綠素b 含量的動態變化圖Fig. 2 The dynamic change of chlorophyll b content among different treatments

2.2.3 N、P 葉面噴施對云南松苗木葉綠素a/b 含量的影響

由圖3 可知，8 月，處理2、6、7、8、9 的葉綠素（a + b）含量均小于處理1，以處理3 最大，處理9 最小。10 月，所有試驗處理中葉綠素（a + b）含量均大于處理1，處理9 中觀測到最大值18.605 mg/g，處理1 中觀測到最小值11.684 mg/g。12 月，處理4、7、8、9 均大于處理1，以處理4 最高、處理6 最小。綜上來看，處理4（中N 低P）最佳，對云南松苗木葉綠素（a + b）含量產生明顯的促進作用。

圖 3 各處理間葉綠素（a + b）含量動態變化圖Fig. 3 The dynamic change of chlorophyll (a + b) content among different treatments

2.2.4 N、P 葉面噴施對云南松苗木葉綠素a/b 的影響

由圖4 可知，不同N、P 葉面噴施對云南松苗木葉綠素a/b 有影響。其中在8 月，除處理2 中葉綠素a/b 含量大于處理1，其他處理均小于處理1，處理2 中觀測到最大值3.146，處理7 中觀測到最小值2.338。10 月，所有試驗處理中葉綠素a/b 均大于處理1，以處理5 最大，處理1 最小。12 月，除了處理9 中葉綠素a/b 大于處理1，其他處理均小于處理1，且以處理9 最大，處理5 最小。綜上所述，8 月和10 月的處理2 以及10 月和12 月的處理9 中葉綠素a/b 大于處理1，其中，6—8 月處理2 中的葉綠素a/b 含量存在最大值，處理2、3、4、5、6、7、8 的葉綠素a/b 除了在10 月大于處理1 之外，其余均小于處理1。

圖 4 各處理間葉綠素a/b 的動態變化圖Fig. 4 The dynamic change of chlorophyll a/b content among different treatments

2.3 N、P 葉面噴施對云南松苗木葉綠素a/b 含量間的相對生長關系分析

由表2 可知，云南松苗木葉綠素a 含量與葉綠素b 含量間的相對生長關系隨觀測時間發生變化，其中處理4 在8 月、10 月、12 月均表現為異速生長關系，處理7 則表現為等速生長關系。處理2 和處理3 在8 月、10 月均為等速生長關系，12 月為異速生長關系。處理5 和處理9 在8 月、10 月均為異速生長關系，在12 月為等速生長關系。處理1 在8 月為等速生長關系，在10 月和12 月為異速生長關系。

各處理在不同測定時間的斜率也發生變化，除處理7、8、9 外，其他處理的相對積累速率隨測定時間發生變化，其中處理8 在8 月、10 月、12 月均表現為大于1，即葉綠素a 含量的積累速率大于葉綠素b 含量；處理7、9 在8 月、10 月、12 月均表現為小于1，即葉綠素a 含量的積累速率小于葉綠素b 含量。其余的處理均隨苗木生長發生變化。處理1、2、3 隨測定時間表現為：8 月斜率小于1、10 月斜率大于1、12 月斜率小于1。處理4、6 在8 月的斜率大于1 而在10 月和12 月均小于1，處理5 在8 月和10 月的斜率小于1 而在12 月大于1。綜合來看，施肥改變葉綠素a 含量和葉綠素b 含量的積累速率，且這種相對生長關系隨時間會發生改變，在單施P 的情況下（處理1、2、3）葉綠素a 含量的積累速率相對于葉綠素b 含量呈現先升高后降低的趨勢，而在單施N 或N、P 配施情況下變化較為復雜。

表 2 云南松苗木葉綠素a/b 含量間的異速生長關系Table 2 The allometric relationship between chlorophyll a and chlorophyll b content of P. yunnanensis seedlings

2.4 N、P 葉面噴施對云南松苗木葉綠素含量的主導因子分析

由表3 可知，N、P 都對云南松苗木葉綠素a 含量起到了不同程度影響，其中，8 月、10 月、12 月的主導因子均為N，促進云南松葉綠素a 含量增加的理論最佳水平組合分別為N3P1、N3P3和N1P3。從葉綠素b 含量來看，8 月、10 月的主導因子均為N，12 月的主導因子為P，促進云南松葉綠素b 含量增加的最優水平組合分別為N3P2、N1P3和N3P3。從葉綠素（a + b）含量來看，8 月、10 月、12 月的主導因子均為N，理論最優水平組合分別為N3P1、N1P3、N3P3。從葉綠素a/b 來看，8 月和10 月主導因子為P、12 月的主導因子為N，總體各水平間差異較少。綜合來看，對葉綠素影響的主導因素以N 為主。

表 3 N、P 葉面噴施對云南松苗木葉綠素的極差分析Table 3 Range analysis of chlorophyll of P. yunnanensisseedlings by foliar spraying of nitrogen and phosphorus

2.5 云南松苗木N、P 葉面噴施肥料效應與最佳施肥量的確定

由表4 可知，所構建的苗木葉綠素含量與肥料用量的二元二次回歸方程，葉綠素a 含量的方程為Y=11.68 + 4.837N-4.154P-12.657N2+ 11.785P2-13.399N × P，葉綠素b 含量的方程為Y=3.781 +3.579N-0.875P-4.191N2+ 4.525P2-9.13N × P。進一步對方程進行降維可知，2 年生云南松適宜的N、P 濃度施用范圍分別為N：0.20%～0.43%，P：0.10～0.17%，單施下的理論產值，對于葉綠素a 而言，適宜的N、P 配比（質量比）為0.20∶0.17；對于葉綠素b 而言，適宜的N、P 質量比為0.43∶0.10。

表 4 N、P 用量與葉綠素含量間的回歸方程Table 4 Regression equation between nitrogen, phosphorus consumption and chlorophyll content

3 結論與討論

N、P 是植物的基本營養元素，其中N 元素是遺傳物質的基礎和蛋白質構成的重要成分，P 是遺傳物質、生物膜、核糖體和能量載體等的組成成分[23]，與植物的光合作用密切相關，而光合作用與葉綠素含量有著密切的聯系。本研究表明，不同N、P 葉面噴施對云南松苗木葉綠素含量具有明顯的影響，其中處理4 可提高葉綠素a 含量、葉綠素b 含量和葉綠素（a + b）含量，這在其他植物的葉面噴施中也有類似的報道[6]，N、P 添加對葉綠素含量有一定的促進作用。但N、P 噴施用量不同，其表現的作用有所差異，前期云南松及相關研究中也有報道[8-11,17,19]，一定N、P 用量可提高葉綠素含量，而高濃度下作用不明顯[10-11,19]，適宜外源N 可以提高葉綠素含量，進而增強植物的光合能力，但高N 水平光合速率反而下降，可能是N 導致植物營養失衡或自遮蔽產生了負面效應[10,24]。但隨著苗木的生長，不同N、P 添加處理對葉綠素含量的影響也發生一定的變化，這可能是因為植物營養需求具有明顯的階段性[16]。當然，葉綠素受環境影響較大，是一個復雜的數量性狀[5-7]，有必要繼續跟蹤調查，并針對苗木的光合特性，進一步開展苗木生長、葉綠素含量等之間的遺傳剖析。

環境中有限的N、P 往往會成為限制生態系統生產力的因素[25]，其中N 含量的高低決定著植物生長速率的快慢[26]。本研究分析表明，對云南松苗木葉綠素影響的主導因子為N，這可能與N 素添加提高葉綠素、保留葉綠體結構等有關[24,27-28]。N 是植物光合作用及生長的限制因子，梁文斌等[24]研究表明，N 添加可以增加植物的可利用N，在一定范圍內，植物的光合速率隨葉片N 含量的增加呈線性增加，進而有利于植物碳的積累和生產力的增加[9]。胡玉玲等[16]利用N、P、鉀（K）等多種元素的肥料對油茶（Camelliaoleifera）進行葉面噴施結果表明，不同生長階段施肥對油茶葉片葉綠素含量的影響最大，各試驗因素對油茶葉片葉綠素含量都存在顯著影響，其中N 肥的作用較為明顯，施N 肥可以明顯提高植物葉片葉綠素含量[29-30]。在云南松的研究中也表明，N 受限比較明顯[31-32]，土壤中N 素營養較缺乏，這也體現了N 素是植物重要的營養元素[16,33]，同時也是植物體內蛋白質、核酸、葉綠素等的重要組成部分[33]，N 參與影響葉綠素的生物合成，一定量的N 條件下，葉綠素合成途徑暢通，葉綠體結構完整，葉綠素含量增加[27,34-35]。

葉綠素含量間的相對生長關系分析表明，N、P 葉面噴施對云南松苗木葉綠素a 含量與葉綠素b 含量間相對生長關系的影響隨時間發生變化，各處理在不同測定時間的斜率也發生變化，這種異速生長關系反映了各屬性之間的表型可塑性[36]，從前期的變異來源分析可以看出，葉綠素a 含量在各測定階段N、P、N × P 下均表現為各處理間差異顯著或極顯著，而葉綠素b 含量在N × P 下表現為各處理間差異極顯著，而在P 下無顯著差異，在N 下僅10 月為差異極顯著。由此表明，葉綠素a 含量與葉綠素b 含量對N、P 葉面噴施的響應不一樣，葉綠素a 的變化較為敏感，這種可塑性體現在葉綠素a 與葉綠素b 間的相對生長關系發生改變，且這種變化隨苗木生長也有所不同[37]。因此，苗木各階段的葉綠素a/b 也發生改變，造成這種差異性的原因有待進一步探討。

綜上所述，N、P 葉面噴施對云南松苗木葉綠素含量有影響，且不同N、P 施用量對葉綠素含量的影響不同，適宜的N、P 用量可提高云南松苗木葉綠素含量。從N、P 兩因素影響效應上來看，N 對綠素a 含量、葉綠素b 含量及其葉綠素（a + b）含量的影響大于P，N 起是影響葉綠素含量的主導因子。不同N、P 葉面噴施也影響了葉綠素a 含量與葉綠素b 含量間的相對生長關系，且隨苗木的生長發生變化。構建了N、P 用量與葉綠素含量間的回歸方程，2 年生云南松適宜的N、P 質量濃度施用范圍分別為N：0.20%～0.43%，P：0.10%～0.17%。