甘肅省扎尕梁北坡頭花杜鵑枝葉性狀特征

楊克彤，陳國鵬*，鮮駿仁，俞筱押，張金武，王立

（1.甘肅農業大學林學院，甘肅 蘭州 730070；2.四川農業大學環境學院，四川 成都 611130；3.黔南民族師范學院旅游與資源環境學院，貴州 都勻 558000）

植物與生境之間的相互關系是當今植物生態學研究的熱點問題［1-3］。植物功能性狀是植物長期進化形成的能夠反映物種對外界環境適應策略的核心屬性，不僅可以揭示植物對環境的響應和適應機制，還可以反映植物對資源的分配和利用規律［3-5］。闡明不同器官功能性狀對異質環境的響應及其相互關系是理解植物進化與適應的基礎［6-10］。加強對環境變化響應敏感器官的研究，可以及時反映植物的生長和適應動態，在揭示物種性狀的進化表達、資源獲取與物種分布等方面具有重要意義。

葉片和小枝作為植物與外界環境接觸最密切的部分，是對環境變化感知最敏銳的構件之一［8-11］，其性狀可塑性能夠較為全面地展示植物對異質環境的適應與生存策略。葉片是光合產出的最重要器官，其功能性狀不僅直接影響植物碳收益，還可以平衡、轉換、彌補植物與外界環境的差異從而決定植物的生存狀態。小枝作為葉片生長的平臺，主要為葉片輸送養分并調整葉片受光面的夾角。此外，小枝的更新發育進一步決定了植物的外部結構和空間特征［11-13］，不僅維持了植物生存空間的穩定，其向外擴展的特性還劇烈影響著物種間的競爭格局［14］。葉片和小枝的自我更新能力明顯優于其他構件，能夠迅速對環境變化做出響應并保持相對穩定［15-16］，為深入探究植物對環境的適應提供了基礎，因其性狀特征具有普適性與易測性，也因此成了研究植物在異質生境中性狀表達的首選指標。

海拔梯度包含了溫度、氣壓、水分和輻射等一系列具有連續分布特征的環境因子，是研究植物適應異質生境變化的最佳場所［1，3，17-21］。不同植物在適應海拔異質生境時的生存策略各不相同，多數植物主要通過功能性狀的權衡與適應來尋求最佳的生存狀態［21-24］，具體表現在：隨海拔上升個體大小差異、器官間生物量分配、異速生長與權衡、枝葉形態異化或內含物含量的波動以及超微結構變化等多個方面［22-25］。由海拔變化引起的環境差異會導致植物產生生理及形態上的多重響應，并逐步經過“瀑布式上升效應”可能會引起整個生態系統的轉變［1，26-27］，因此海拔梯度細微的變化有時也會帶來巨大的結果差異。雖然近年來部分研究在群落功能性狀隨海拔梯度變化方面已有所成果，但部分結論存在差異甚至相反，并且群落水平的性狀變化掩蓋了物種的個體特征，導致其對生態系統的預測精度不高［28-29］，而在單物種尺度上開展枝葉性狀與環境關系的研究可以深入了解關鍵植物適應機制并凸顯環境對其功能性狀的影響。

山地生態系統是對全球氣候變化最敏感的區域，是研究植被與氣候關系的天然實驗室［1，30］。頭花杜鵑（Rhododendron capitatum）作為典型的高山植物，是維持青藏高原東北緣山地生態系統穩定的重要類群。本研究以甘肅扎尕梁北坡的優勢種頭花杜鵑為研究對象，測定其23組枝葉功能性狀，探究海拔梯度下頭花杜鵑枝葉性狀的適應與響應模式，以進一步解析高山植物對海拔梯度的應對策略。基于此，本研究試圖回答以下問題：1）隨海拔梯度上升，頭花杜鵑枝、葉性狀的變化有何特征與規律？2）性狀間的相互關系如何？3）不同性狀的表型分化有何特征？

1 材料與方法

1.1 研究區概況

研究區位于甘肅省舟曲縣插崗梁自然保護區（103°57′05″-104°42′05″E，33°14′32″-33°53′52″N，海拔1173～4536 m），保護區處在青藏高原東北緣與秦嶺西翼及岷山山脈交匯地帶，屬西秦嶺地質構造帶南部隴南山地［7］。該區域處于我國溫帶季風、亞熱帶季風和高原山地氣候帶的匯合區域及半濕潤和半干旱區的分界線上。近30年來，林區年平均降水量為951.0 mm，年平均蒸發量918.8 mm，相對濕度為82.0%，干燥度0.5，年平均氣溫4.3℃，最冷月（1月）平均氣溫為-13.3℃，最熱月（7月）平均氣溫為20.8℃，全年無霜期96.7 d（最多143 d，最少53 d）；日平均氣溫≥5℃的日數為178 d，年積溫1994.6℃；日平均氣溫≥10℃的日數為103 d，年積溫1387.9℃。年日照時數為1398.4 h，日照百分率為32.0%。植被分布的垂直帶譜明顯，從低到高依次為干旱河谷灌叢、亞高山次生林、亞高山暗針葉林和高山杜鵑灌叢［25］。頭花杜鵑植株主要分布在該區域海拔3100～3700 m范圍內，在3400 m以下主要為單株零星分布或可形成小規模種群，3400 m以上與隴蜀杜鵑（Rhododendron przew alskii）混生形成灌叢群落。

1.2 數據采集與測定

于2019年8月在研究區進行全面勘查后，在插崗梁保護區的扎尕梁北坡的林線到山頂處自下而上沿等高線按海拔高度3400、3500、3600和3700 m設置調查樣線。在樣線中隨機選取株高50 cm左右的標準木，以保證植株年齡基本一致，并進行小枝與葉片的采樣。每個海拔梯度共采集30余個小枝，并放于自封袋內，葉片夾放在濕潤濾紙中，并置于便攜式低溫箱中帶回。本研究所述小枝定義為植株頂端無分支的1年生莖。

小枝與葉片帶回后即刻測定各項指標，每個海拔隨機測量15個樣本。對于小枝，首先剪下小枝上的葉片，然后測量小枝基徑（stem diameter，SD）、小枝長（stem length，SL），再用掃描儀（EPSON V 700）對所有葉片進行掃描并求出總葉面積（total leaf area，TLA）和葉片數（leaf number，LNO.），用電子天平測量小枝鮮重（stem fresh weight，SFW）、小枝干重（stem dry weight，SDW）、總葉鮮重（total leaf fresh weight，TLFW）和總葉干重（total leaf dry weight，TLDW）；對于單葉，先用天平稱量葉片鮮重（leaf fresh weight，LFW），再用掃描儀進行掃描并求出單葉面積（leaf area，LA），隨后測量葉厚（leaf thickness，LT）與葉干重（leaf dry weight，LDW），測量葉厚時靠近葉片中心但避開主脈，測量3次求取平均值后作為葉片厚度。測量小枝與葉片干重時，將其置于玻璃皿放入烘箱，120℃殺青30 min后轉為75℃烘至恒重后稱取干質量。葉長（leaf length，LL）、葉寬（leaf width，LW）、單葉面積（LA）與總葉面積（TLA）均使用Image J面積測量軟件進行。鮮、干質量的測定均使用電子天平（精確至0.0001 g），長、寬、厚、徑均使用數顯游標卡尺（精確至0.01 mm）測量。由于頭花杜鵑葉柄極短，難以測量，所以用解剖刀將葉柄除去而不納入測定指標，因此本研究所述葉片性狀均不包括葉柄。

相關枝葉性狀的推算如下［7，10］：葉形指數（leaf shape index，LSI）=葉長/葉寬，比葉面積（specific leaf area，SLA）=葉面積/葉干質量，葉干物質含量（leaf dry matter content，LDMC）=葉干重/葉鮮重，莖干物質含量（stem dry matter content，SDMC）=莖干重/莖鮮重，葉體積（leaf volume，LV）=葉面積×葉厚，莖體積（stem volume，SV）=π×莖長×（莖徑/2）2×1/3，葉密度（leaf density，LDE）=葉干重/葉體積，莖密度（stem density，SDE）=莖干重/莖體積，纖細率（stem slender ratio，SSR）=莖長/莖徑，總比葉面積(total specific leaf area,TSLA)=總葉面積/總葉干質量，總葉干物質含量(total leaf dry matter content,TLDMC)=總葉干重/總葉鮮重，總葉莖重比（total leaf and stem weight ratio，TLSWR）=總葉干重/總莖干重。

1.3 數據分析

采用單因素方差分析（one-way ANOVA）檢驗海拔高度對頭花杜鵑枝葉功能性狀的影響，并采用Tukey’s HSD多重比較檢驗性狀在不同海拔間的差異顯著性；利用方差分析提取海拔間與海拔內及隨機誤差的方差分量并計算各性狀的表型分化系數（Vst），表型分化系數Vst=δ2t/s/（δ2t/s+δ2s），其中δ2t/s為海拔間方差值，δ2s為海拔內方差值［31-33］；采用Pearson相關系數、雙尾檢驗進行海拔與性狀間的相關分析，最后用單因素方差分析檢驗海拔內與海拔間方差分量的顯著性。利用Excel 2010和SPSS 20.0進行數據統計、分析，用Origin 9.1繪圖。

2 結果與分析

2.1 頭花杜鵑枝葉性狀隨海拔的變化

頭花杜鵑葉性狀隨海拔梯度變化差異顯著，除LSI無顯著變化外，其他性狀均有顯著變化（圖1），LL、LW、LA、LV和LDW在3600 m時顯著小于其他海拔；LDMC和LDE隨海拔上升逐漸增大，并在海拔3700 m時顯著高于其他海拔；而SLA逐漸減小，海拔3700 m時顯著小于其他海拔；LT整體呈增大趨勢。小枝性狀除SD、SDMC、LNO.和SDE隨海拔梯度無顯著變化外，其他性狀均變化顯著（圖2），其中SL、SDW、TLDW、TLA、SSR和SV在3600 m顯著小于其他海拔，TSLA呈下降趨勢，TLDMC呈上升趨勢，TLSWR則先增大后減小。

圖1 頭花杜鵑葉性狀隨海拔的變化Fig.1 Variation of R.capitatum leaf traits with altitude

圖2 頭花杜鵑小枝性狀隨海拔的變化Fig.2 Variation of R.capitatum stem tr aits with altitude

2.2 枝葉性狀間的相互關系

頭花杜鵑大部分枝葉性狀間均為顯著正相關（圖3），LL、LW、LA、LV和LDW均與SD、SL、SDW、TLDW、TLSWR和SV顯著正相關，LT、LSI和SLA與多數性狀顯著負相關，小枝水平的SDMC、TLSWR、LNO.、SSR和SDE與部分性狀的相關性較弱，總體而言小枝水平不及葉片水平的相關性強。

圖3 頭花杜鵑枝葉性狀間的相關關系Fig.3 Relationship between stem and leaf traits of R.capitatum

2.3 枝葉性狀的表型分化

頭花杜鵑枝葉性狀在海拔間和海拔內的方差分量百分比均值分別為17.30%和70.76%，隨機誤差為11.94%（表1），表型分化系數整體介于3.19%～57.73%，均值為20.65%，僅8組性狀表型分化系數高于均值，占所有性狀（23組）比例的34.78%。頭花杜鵑在同一海拔內的方差分量顯著高于海拔間，其中LSI、SD、SDMC、LNO.和SDE的表型分化系數較低（小于11.11%，下四分位數），LT、LDW、LDMC、SLA和LDE的表型分化系數較高（大于25.41%，上四分位數），即在應對海拔梯度引起的異質生境時，這10組性狀分別表現出較低的變異性與較高的可塑性，表明頭花杜鵑枝葉性狀變化范圍大、表型變異豐富。

表1 頭花杜鵑枝葉性狀的方差分量與表型分化系數Table 1 Var iance component and phenotype differentiation coefficient of stem-leaf tr aits of R.capitatum

3 討論

海拔梯度引起的生境變化會顯著影響植物功能性狀，并且部分性狀也會隨之表現出梯度變化，這取決于性狀對環境的敏感程度與自身功能特性［32-34］。本研究中，頭花杜鵑葉長、葉寬、葉面積、葉體積、葉干重、莖長、莖干重、總葉干重、總葉面積、纖細率和莖體積均在3600 m處顯著減小，整體呈“V”形變化（圖1和圖2）。可能是由于山地植被逐步從林線過渡到樹線、再過渡到高寒杜鵑灌叢后，杜鵑灌叢抗環境干擾能力強，土壤理化性質相對穩定，頭花杜鵑產生適應性策略，因此3600 m可能是生境適應連續變化梯度的轉折點［25］。葉厚、葉干物質含量、葉片密度和總葉干物質含量均隨著海拔梯度上升而增大，與以往研究結果一致［34-35］。隨著海拔升高，溫度降低、光照輻射增強，葉片將資源集中于構建防御結構［36］，一方面可以降低輻射傷害，另一方面增加葉厚有助于抵抗低溫凍害。莖干物質含量降低則可能與單葉面積減小有關，Corner［28］法則認為：大葉片需要較大的枝條支撐并輸送養分，隨著海拔升高，單葉面積減小，由此導致植株對莖的投資減少，這種策略實現了有限資源的最大利用效率。通常情況下，分布在資源有限生境及生長緩慢的物種比葉面積都較低，不論草本、灌木與喬木，還是極地、溫帶與熱帶，植物的比葉面積變化都隨著海拔上升而降低［2-3，35］，導致這種情況的原因則主要與溫度下降有關，低溫環境下，葉肉細胞增大、表皮細胞壁增厚，導致葉片厚度增加，因此比葉面積下降［14］。

植物在應對異質環境時，并非某一性狀采取適變策略，而是多個性狀間的組合（即功能性狀組）產生共變，這種共變機制更有利于應對環境引起的激變［22，30］。本研究中，僅7組性狀與海拔梯度為極顯著（或顯著）相關，其余16組性狀均不顯著，即這7組性狀可能是枝葉功能性狀的組合中心。由這7組性狀產生較為劇烈的應對策略并組合其他性狀產生微小變化便足以應對海拔梯度引起的生境異質，這種應對策略減少了能量損耗與成本投資，具有高效率、低能耗的特點［22］。此外，葉長、葉寬、葉面積、葉體積、葉干重均與莖徑、莖長、莖干重、莖體積、總葉干重極顯著正相關，說明頭花杜鵑枝葉性狀支持Corner［28］法則，即較大的葉片需要大的枝條來承擔；另一方面可能隨著海拔升高，溫度降低，土壤水分難以利用，導致頭花杜鵑產生干旱脅迫，加之較小枝條在低溫環境運送水分與養分時容易產生栓塞，較大且健壯的枝條有助于緩解栓塞并抵御高山地區大風引起的損傷［29-30］。

枝葉性狀的變異與表型分化共同決定了其變異程度與范圍，較強的變異與較大的可塑性有助于應對高山地區極端氣候。本研究發現，頭花杜鵑枝葉性狀表型分化系數均值為20.65%（表1），低于紅椿（Toona ciliata，80.73%）、砂生槐（Sophora moorcroftiana，39.25%）、珙桐（Davidia involucrata，37.00%）、尾葉櫻桃（Cerasus dielsiana，30.78%），而高于云南紅豆杉（Taxusyunnanensis，14.66%）、紅松（Pinuskoraiensis，12.39%）、馬尾松（Pinus massoniana，6.44%），說明其表型變異程度低于一般的闊葉樹種卻高于針葉樹種。有研究認為生存環境會極大地影響表型分化程度［31-35］，頭花杜鵑長期生活于低溫與強光輻射等不利環境中，環境因子的多變與復雜可能是造成其表型分化相對保守的主要原因。另外，枝葉性狀中僅8組（葉寬、葉厚、葉面積、葉體積、葉干重、葉干物質含量、比葉面積和葉密度）性狀的表型分化系數高于平均值（20.65%），說明這8組性狀具有相對豐富的表型變異，尤其比葉面積與葉密度表型分化系數高達55.14%和57.73%，可能是應對高山環境海拔異質的核心性狀［32-36］。葉形指數、莖徑、莖干物質含量、葉片數和莖密度表型分化系數較低（<11.11%），說明這些性狀具有穩定特質，并不輕易隨海拔梯度做出改變，其穩定特質可能對于維持植物枝葉生長發育有特殊意義，如：極端環境中維持生存的下限等［13-17］。另外，頭花杜鵑枝葉性狀海拔內變異普遍高于海拔間變異，即同一海拔內的居群是產生豐富表型變異的主要來源，可能是應對生境異質的主要貢獻者，也可能由于在研究區范圍內頭花杜鵑分布集中，并未產生地理隔離，從而導致海拔間居群表型變異較弱。

頭花杜鵑枝葉功能性狀對環境的適應有助于其應對高寒生境，可在一定程度上改善其存活、生長和繁殖狀態。頭花杜鵑73.91%（17組/23組）的枝葉性狀均在3700 m顯著增大，這種策略可能是為了避免較低海拔（3400 m）的天敵或中海拔（3500～3600 m）的物種競爭；也可能因為頭花杜鵑屬于典型的高山植物，在長期的生物演化與自然選擇過程中已經適應了高山特殊環境，其固有的生態特性可能更有助于頭花杜鵑在更高海拔而非低海拔生存［7-12］。

4 結論

頭花杜鵑主要以枝葉功能性狀組合與協變共同應對高山生境的不利因素；部分性狀在海拔3600 m顯著減小，大部分性狀在海拔3700 m達到最大；海拔內表型變異普遍大于海拔間，葉片的表型分化系數顯著大于小枝。頭花杜鵑枝葉功能性狀對高山生境的積極響應與適應在提高頭花杜鵑生境適合度、占據穩定生態位、維持高山物種競爭格局與生物多樣性等方面具有重要意義，在未來的研究中應繼續加強對高山物種的關注。