煙株葉片葉綠素的時空分布特征

肖慶禮 ，代先強，楊紅柯，秦平偉，張世渠，陳昆燕 ，汪長國 ，戴亞 *，丁偉

1 西南大學植物保護學院，重慶市北碚區天生路2號 400715；2 重慶中煙工業有限責任公司，重慶市南岸區南坪東路2號 400060；3 煙葉資源科學利用重慶市重點實驗室，重慶市南岸區南坪東路2號 400060；4 重慶市煙草公司，重慶市江北區五江路20號 400020

烤煙收獲的是葉片，在烤煙產量的形成過程中，葉綠素參與的光合作用貢獻了96%的干物質積累[1]。與此同時，在烤煙采收環節葉綠素也是判斷葉片成熟采收標準的主要指標之一[2-5]。在烤煙生長發育過程中，葉綠素不僅隨時間呈現明顯的變化，而且受到葉位、品種、栽培措施、地形地貌等各種主客觀因素的影響，因此烤煙葉綠素呈現較強的時空變異性特征[6-11]。關于時空變異性，國內外學者開展了大量的研究工作[12-15]。部分學者在研究土壤水分時空變異性時發現，土壤水分的空間分布格局具有隨時間持續不變的特性，Vachaud 等[16]最早將這種特性定義為時間穩定性。對時間穩定性最好的應用方式之一是尋找研究區的代表性樣點，其含水量可以代替研究區平均含水量[17-21]。

鑒于土壤水分時間穩定性應用的巨大潛力，一些學者也把時間穩定性的概念擴展應用到了的其它方面。徐國策等[22]研究了洛惠渠灌區地下水的電導率，發現地下水電導率具有很強的時間穩定性和周期性，可以利用地下水電導率代表性位置點來監測研究區地下水平均電導率的變化。邢旭光等[23]分析了膜下滴灌棉田主根區土壤的鹽分，發現土壤鹽分的時間穩定性隨土層深度增加呈現先增強后略微減弱趨勢。Douaik[24]分析了土壤鹽分的時間穩定性，發現低鹽點的時間穩定性較強，而高鹽點的時間穩定性則較弱。

目前，關于時間穩定性的擴展應用研究主要集中在土壤鹽分、水分電導率等方面，關于煙葉葉綠素時間穩定性的相關研究鮮見報道。由于煙葉葉綠素含量與SPAD 值之間存在極顯著線性關系[3,26-27]，本文采用Spad-502 plus 便攜式葉綠素儀測定大田煙株葉片的SPAD 值，進而換算出葉片的葉綠素含量。通過系統研究摩圍山小流域旱地煙株葉片葉綠素的時空分布特征，確定監測樣帶煙株葉片葉綠素的代表性樣點，從而為本區煙田的水肥管理及成熟采收提供理論支撐。

1 材料與方法

1.1 研究區概況

研究區位于重慶市彭水縣南部的摩圍山小流域。該區屬中亞熱帶溫潤季風氣候區，具有春早多夜雨，夏熱多伏旱，秋涼多綿雨，冬冷少霜雪，總體雨量充沛集中，光照偏少云霧多等氣候特點。該區年均氣溫17.5℃，年平均降雨量1104 mm，年平均蒸發量950 mm，降雨量略大于蒸發量[25]。該區屬于典型的巖溶地貌，地形起伏，溝壑縱橫，低階地、坡麓、巖溶洼地及小型山間盆地犬牙交錯。主要土壤類型為旱地黃壤，土層淺薄，水土流失較為嚴重。

1.2 試驗布設

在研究區選取典型的小型山間盆地作為試驗區，盆地長約520 m，寬50～150 m不等，為東北-西南走向，且東北至西南有一個約5°的下斜緩坡。農田土壤類型為旱地黃壤，前作為烤煙?？緹熞圃郧捌饓挪⒏采w地膜，壟體高度約為30 cm。施純氮量控制在105 kg/hm2，N : P2O5: K2O 為1 : 0.93 : 3.08。移栽時間為5月1 日-7 日，品種為云煙87。株行距為55 cm×110 cm，初花期打頂，有效留葉數為19-21 片。其它大田管理措施按優質煙葉生產技術要求進行。農田無灌溉條件，煙株耗水全部來自于天然降水和土壤水。在盆地東北-西南方向依次布設20 個監測樣點，樣點間隔距離為25 m（圖1b），在每個監測樣點分別選擇1 棵煙株，該煙株的長勢長相與該樣點附近的煙株長勢長相基本一致，然后掛牌進行長期定位監測。

圖1 （a）研究區位置、（b）樣帶布設示意圖Fig. 1 (a) Location of the study area, (b) The layout of the experimental plot

1.3 數據采集和處理

煙葉SPAD 值的測定使用Spad-502 plus 便攜式葉綠素儀進行。移栽結束后39 d 開始測定，測定時對樣點掛牌煙株按照自下而上逐葉位的方式進行。煙株最低的葉片為第一葉位，葉片自下而上葉位依次增加。烤煙生長期間各葉片的葉位固定不變。

烤煙生長期間共測定了18 次，每次間隔為5 d左右。測定時沿葉片主脈對稱的葉基、葉中和葉尖部兩邊各選3 個部位等距布局，每片葉共測定6 個部位，其平均值即為該片葉SPAD 值。

大量研究表明，煙葉葉綠素含量與SPAD 值存在極顯著線性關系，并建立了數學模型[3,26-27]。本研究借鑒汪強等[27]數學模型對煙葉葉綠素進行標定：

式中，y 為煙葉葉綠素含量，mg/g Fw；x 為煙葉SPAD 值。

煙株葉片按照下、中、上分為3 個部位，其中下部葉、中部葉和上部葉分別占煙株總葉數的30%、40%和30%。由于監測煙株葉片數為19～21 片，煙株葉數為19 片時，其下、中、上部葉分別為1～6 葉位、7～13 葉位和14～19 葉位；煙株葉數為20 片時，其下、中、上部葉分別為1～6 葉位、7～14 葉位和15～20 葉位；煙株葉數為21 片時，其下、中、上部葉分別為1～6葉位、7～15 葉位和16～21 葉位。煙株下、中、上部葉片葉綠素含量分別為對應部位葉片葉綠素含量的平均值。

1.4 分析方法

根據Vachaud 等[16]的研究，時間穩定性可以使用相對偏差法描述。給定時間特定樣點特定煙株部位葉綠素含量的相對偏差δikj 可以由下式計算得出：

式中，i 代表樣點位置（i=1, 2, 3, ……20），j 代表煙株部位（j=1, 2, 3），k 代表測定次數（k=1, 2, 3, ……M），為時間k 時，j 部位的平均葉綠素含量，計算公式為：

時間穩定性指數（ITS）可以作為單一指標來識別研究區域平均葉綠素含量的最佳樣點，同時該樣點具有較好的時間穩定性，公式如下：

變異系數Cv反映的是相對變異，即隨機變量的離散程度：

式中σ 為標準差，μ 為樣本平均值。根據變異系數將葉綠素含量的變異程度進行分類，當Cv＜10%為弱變異性，10%≤Cv≤100%為中等變異性，Cv＞100%為強變異性。

利用回歸方程的決定系數（R2）、均方根誤差（RMSE）和平均偏差（MBE）評價代表性樣點和對應葉片葉綠素的差異性:

RMSE 越小，MBE 越接近于0，說明代表性樣點煙株葉片葉綠素值越接近于樣帶對應部位煙株葉片葉綠素的平均值。

采用Excel 2013 和Surfer11.0 軟件對數據進行處理和作圖，采用SPSS Statistics 20.0 軟件對數據進行相關性分析。

2 結果與分析

2.1 煙株葉片葉綠素的時空變化特征

煙株不同部位葉片葉綠素均值及變異系數見圖2。由圖2a 可知，下部、中部和上部葉在生育前期均具有較高的葉綠素含量，隨著葉片的生長，葉綠素含量均明顯增加，且都達到最大值，其最大值分別為0.92、1.02 和1.18 mg/g Fw，隨后均進入下降通道，至生育末期下降至最低值，其最低值約為最高值的55%。就在整個生育期而言，下部葉葉綠素平均值（0.76 mg/g Fw）明顯小于中部葉（0.82 mg/g Fw），中部葉（0.82 mg/g Fw）又明顯小于上部葉（0.97 mg/g Fw）。

由圖2b 可知，在生育前期，下部葉葉綠素變異系數較低，其值僅為5.3%，隨著葉片的發育，變異系數明顯增大，最后達到最大值26.67%；中部葉在生育前期葉綠素變異系數較低，隨后快速增大，然后又明顯下降，最后進入上升通道，其最大值為19.74%；上部葉在生育前期葉綠素變異系數稍大，隨后快速降低至最低值（5.28%），最后進入上升通道，其最大值為23.28%。下部、中部和上部葉葉綠素變異系數雖然變化規律不同，但進入中后期其數值均隨葉片的生長發育而增大，后期均達到最大值?？傮w而言，在生育前期和中期，各部位葉綠素的變異系數基本低于10%，屬于弱變異；在生育后期，各部位葉綠素的變異系數均在10%～100%之間，屬于中等變異。

圖2 煙株葉綠素含量及變異系數的時間序列特征Fig. 2 Time series character of chlorophyll content and coefficient of variation in tobacco plan

為了更直觀的呈現煙株在生長期間葉綠素的動態變化，使用Surfer11.0 軟件繪制監測樣帶煙株葉綠素的等值分布圖。由圖3 可知，在樣帶尺度，下部葉生育前期，葉綠素等值線較為稀疏平滑；隨著葉片的發育，葉綠素值緩慢增加，其最大值出現在樣帶中部區域；隨后葉綠素等值線由疏變密，葉綠素值快速降低，其最低值出現在樣帶的中上部和中下部區域。中部葉生育前期，葉綠素等值線較為密集，樣帶上葉綠素值呈現斑塊狀分布的高值區域；隨著葉片發育，樣帶上葉綠素值呈斷崖式下降，生長后期葉綠素等值線較為碎片化。上部葉生育前期，樣帶上葉綠素值呈現斑塊狀分布的高值區域；隨著葉片的發育，樣帶上葉片葉綠素值進入緩慢下降通道，生長后期葉綠素等值線較為碎片化?？傮w而言，樣帶上煙株不同部位葉片生長發育期間，葉綠素值的變化趨勢具有相似性。

圖3 生長期間不同部位葉片葉綠素等值分布圖Fig. 3 Isoline distribution of chlorophyll content in different parts of tobaco leaves during growth period

2.2 煙株葉片葉綠素的時間穩定性分析

由圖4 可知，煙株下部、中部和上部葉葉綠素平均相對偏差分別介于-17.02%～20.81%、-12.59%～22.06%和-11.48%～18.15%，平均相對偏差的標準差分別介于2.82%～23.74%、4.04%～15.02%和2.89%～17.33%。時間穩定性指數（ITS）越小，說明樣點煙株葉綠素的時間穩定性越好，其最小值即可作為代表性樣點[28]。因此，下部、中部和上部葉葉綠素的代表性樣點分別為樣株15、樣株14 和樣株5。

圖4 煙株葉綠素的平均相對偏差；誤差棒表示平均相對偏差的標準差；空心標示的樣點為代表性樣點；刻度線為時間穩定性指數Fig. 4 Ranked mean relative differences (MRD) of chlorophyll content. Error bars represent the standard deviation of relative difference (SDRD). The most time-stable locations are indicated by the empty symbols. The curve indicates the index of time stability (ITS)

圖5a 顯示的是20 個樣點煙株下部、中部和上部葉葉綠素的平均相對偏差值。下部、中部和上部葉葉綠素的平均相對偏差的極差值分別為37.83%、34.65%和29.63%，變化分別介于-17.02%～20.81%、-12.59%～22.06%和-11.48%～18.15%。葉綠素平均相對偏差的極差值隨煙株部位的上升而逐漸減小。

圖5b 顯示的是所有樣點煙株葉綠素時間穩定性指數的累積值。從圖5b 可知，葉綠素累積的時間穩定性指數最大值（樣點7）為79.05%，最小值（樣點1）為19.90%，平均值為32.83%。就各部位而言，下部葉葉綠素的時間穩定性指數值最大，上部葉次之，中部葉最小。

圖5 所有樣點葉綠素含量的平均相對偏差（a）和時間穩定性指數（b）的累積值Fig. 5 Histogram of cumulated mean relative difference (MRD, a) and index of temporal stability (ITS, b) of each sampling point

2.3 代表性樣點的合理性檢驗

圖6 代表性樣點葉綠素含量與平均值的比較Fig. 6 Comparison between measured chlorophyll content at the most time-stable location (MTSL) and the estimated values

圖6 顯示的是代表性樣點葉綠素含量與對應部位葉綠素的平均含量。下部、中部和上部葉片葉綠素回歸方程的決定系數（R2）分別為0.9032、0.9509和0.9693，均方根誤差（RMSE）分別為0.0471、0.0359 和0.0476mg/g Fw，平均偏差（MBE）分別為-0.0038、0.0118 和-0.0366mg/g Fw。相關分析表明，代表性樣點的葉綠素含量均與對應部位平均葉綠素含量之間存在顯著性相關（P＜0.01）。因此可以用代表性樣點煙株葉綠素值估算對應部位葉綠素的平均值。

3 討論

葉綠素是煙草進行光合作用的重要色素，葉綠素含量不僅反映煙草的光合強度和營養水平，而且還與煙葉的成熟采收密切相關[1-5]。本研究結果顯示，在葉片整個生育期，葉綠素的動態變化表現為前期具有較高值，且隨時間快速增加，中后期則進入長期下降通道，至采收降至最低值。馮曉英等[6]通過對不同烤煙品種葉綠素研究，童淑媛等[29]通過對玉米葉綠素研究和黃瑞冬等[30]在玉米和高粱葉片葉綠素的研究均得出類似結論。原因可能是煙草葉片發育的前中期，較高的葉綠素含量有利于葉片干物質的積累[31]；葉片發育的后期，葉綠素的大量降解則有利于內含物質的代謝和煙葉香氣物質形成，并促使葉片呈現成熟特征[32-33]。本研究顯示，隨著葉片部位的上升，葉綠素含量明顯增大，其原因可能是植物的頂端優勢以及更好的光照條件導致煙草部位越高葉片光合作用越強，葉片厚度越厚，葉綠素含量也隨之越大[34]，李佛琳等[3]研究也得出類似的結論。

山區土層厚度不均[35]、土壤肥力多變[36]、坡度坡向各異[37]，這導致煙株長勢不均，進而造成葉片葉綠素具有較強的空間變異性。煙農在水肥管理和田間管理時則會采取針對性措施以盡量減緩這種不利影響，如對土層較薄、肥力較低、煙株長勢較弱的地塊采取施偏心肥以保障煙株的正常發育；在打頂留葉時也會根據煙株個體的營養水平分類打頂和合理留葉：對長勢稍弱的煙株扣心打頂并適當少留一兩片葉，長勢正常的采取初花打頂并合理留葉，而對長勢較好的煙株則采取盛花打頂并適當多留一兩片葉。這種差異化的施肥[36]和針對性的田間管理[38]等人為措施大幅提升了山區煙田煙株的整齊度，進而顯著降低煙株葉綠素空間的變異性。本研究中煙株葉綠素雖然具有變異性，但其變異性表現為弱變異和中等變異，樣帶尺度葉綠素呈明顯碎片化分布，但依然具有一定的規律性，這是自然因素和人為因素彼此相互影響的結果。

葉綠素的時間穩定性受到平均相對偏差和平均相對偏差標準差的雙重影響[23]。平均相對偏差的波動越小，葉綠素隨時間的穩定性越好；平均相對偏差標準差值越小，葉綠素的時間穩定性越強[39]。本研究顯示，中部葉片葉綠素的時間穩定性最強，上部葉次之，下部葉最弱。其原因可能是煙株在發育過程，中部葉片光照條件適宜，薄厚適中，受到水肥及田間管理等措施的影響相對較小，而下部葉光照條件較差，葉片較薄，受到的影響相對較大。市場上中部葉片具有最好的品質和價格也從側面印證了本研究結論的合理性[40]。

本研究針對摩圍山小流域典型旱地煙田煙株葉綠素含量選取的代表性樣點并對樣點煙株的葉片葉綠素進行估算，回歸分析和相關性分析表明，代表性樣點煙株葉片的葉綠素可以較好反映研究區煙田煙株對應部位葉片葉綠素含量，這在本區煙田水肥管理和成熟采收過程中具有重要作用?；跁r間穩定性選取的代表性樣點在不同年際間同樣具有較好的穩定性[22-23]。由于本區煙田常年由固定的煙農進行種植，田間管理在年際間也具有較好的一致性，因此選取的代表性樣點煙株葉片葉綠素含量在不同年際間也具有較好的穩定性。因此，在以后煙田的水肥管理和成熟采收過程中，只需測定代表性樣點煙株的葉綠素含量，再結合該煙株的長勢長相以及葉片的成熟特征即可判斷本區煙田的整體營養水平和待采葉片的田間成熟度。這為有效降低山區煙田水肥管理的難度，進一步提升煙葉的田間成熟度提供了一條便捷途徑。

4 結論

該文在摩圍山小流域，通過對旱地煙田煙株葉綠素進行了時空分布特征研究，得出以下結論：

（1）煙株葉片葉綠素的動態變化呈“倒勺子”型，且葉綠素值隨葉片部位的增加而增大；葉片發育前中期葉綠素值呈弱變異，后期則呈中等變異。

（2）煙株下、中、上部葉片葉綠素值均具有較好 的 時 間 穩 定 性（R2為0.9032～0.9693，RMSE 為0.0359～0.0476mg/g Fw，MBE 為-0.0366～0.0118mg/g Fw），利用代表性樣點煙株葉綠素值可以表征對應部位葉綠素的平均值。

（3）在煙田的水肥管理和成熟采收中，可利用代表性樣點煙株的葉綠素含量，以及該煙株的長勢長相和葉片的成熟特征即可判斷本區煙田的整體營養水平和待采葉片的田間成熟度。