海拔梯度下元謀干熱河谷植物群落特征

余嬌娥，司宏敏，吳雪濤，張彧，蘇文華，周睿*

1. 云南大學生態學與地植物學研究所，云南 昆明 650500；2. 吉林大學環境與資源學院，吉林 長春 130023

植物群落是植物在長期環境變化中自我調節和適應而形成的集合體（Lieth，1975），物種多樣性是群落最基本的特征，是形成群落結構的基礎。研究群落的物種組成及多樣性是認識植物群落結構、生態、動態的基礎（王國明等，2017；許傳陽等，2014）。一般認為，物種多樣性是由資源可利用程度決定的（Morin，2000），會隨著生境的干旱的程度而降低（Bello et al.，2006；劉小丹等，2014）。α多樣性以及β多樣性構成物種的總體多樣性。在植物群落多樣性研究中，群落α多樣性是反映群落內物種組成狀況的一個重要指標，即物種豐富度和均勻度，β多樣性表示物種沿環境梯度發生替代的程度或物種變化的速率，表示生境的變異程度或指示生境被物種分隔的程度（馬克平等，1995）。植物群落α、β多樣性是當前群落生態學乃至整個生態學研究中十分重要的內容（David，1994；馬斌等，2008；陳杰等，2012；張曉龍等，2017）。研究植物群落多樣性特征對揭示群落、生態系統結構和功能演變具有極為重要的意義（Pandey et al.，2003）。

植被蓋度（Vegetation coverage）通常是指植被冠層垂直投影面積占基準地表單位面積的比例或百分比（Purevdorj et al.，1998），是表示陸地表面植被數量及反映生態系統變化的重要基本參數，也是地面蒸散、光合作用等多種地表過程研究的控制性因子，也是反映森林資源的一個重要參數（Price，1993；張云霞等，2003）。植被覆蓋及其變化不僅是區域生態系統環境變化的重要指標，而且對水文、生態、全球變化等都具有重要的意義（陳曉光等，2007）。

干熱河谷是對分布于中國西南地區，具備干、熱兩個基本屬性河谷帶狀區域的總稱，該區域由于地形深陷封閉而產生的雨影作用以及海拔下降絕熱增溫和谷地焚風效應而形成了西南地區特殊的地形氣候類型（金振洲等，1987）。所謂“干熱”就是極端水分條件與熱量條件的結合。其中，“干”方面，全年的蒸發量大于降水量的3～6倍；“熱”方面，年均溫＞20 ℃，＞10 ℃年積溫＞7000 ℃。全年干濕季分明，水分和光熱條件失衡，熱量高、降水量少以及蒸發強是干熱河谷氣候的典型特征（何永彬等，2000；張建平等，2000）。這一特殊氣候特點導致該區域谷底（海拔約1100 m）植被為稀樹灌草叢組成的“河谷型”薩瓦納植被（Savanna of valley type，金振洲，1999），而河谷頂部（海拔約2100 m）則分布有常綠闊葉林。從植物區系來看，1600 m以下河谷中植物有明顯的熱帶性質，而河谷頂部亞熱帶成分增加（歐曉昆，1988），表現出明顯的物種更替現象。劇烈的氣候和群落變化被壓縮在1000 m的海拔落差中，其間物種適應特征、物種組成和生物多樣性也必然發生重大改變。該地區早期就有人類活動，是人類文明的發源地之一，是中國西部人口密集、開發較早的山地農業區，是云南農業經濟發展熱區（朱鐘麟等，2006）。社會資源的高度開發和利用以及其他人為因素（如放牧、開墾耕地等）直接或間接造成了生態環境的進一步退化和破壞（張建利等，2010）。且從 20世紀 50年代至今，該地區大肆開展植物造林工程，主要引進種植了以赤桉（Eucalyptus camaldulensis）、檸檬桉（Eucalyptus citriodora）、大葉相思（Acacia auriculaeformis）以及豆類等物種（李昆等，1999；李昆等，2004），迅速而直接地改變了干熱河谷的景觀類型、組分結構及空間格局。該區域干濕季分明的氣候特征也造就了植被外貌明顯的季節性變化，旱季很多物種采用避旱對策呈現假死狀態，只有很少完全適應當地極端干熱環境的物種能在旱季保持常綠，這些物種是組成當地植物群落的重要物種，是維持生物多樣性的基石。此外，雨季有人工種植物種逸出，這些物種未必適宜當地氣候，會對生物多樣性調查造成一定干擾。

本研究從干熱河谷底部到兩側山地沿海拔梯度設置樣方，探索沿海拔梯度群落結構的變化。試圖回答以下科學問題：該地區的群落結構變化是怎樣的？物種替代現象主要發生在哪個海拔區間？如何劃分干熱河谷與常綠闊葉林植被的分布區域？研究結果不僅可為干熱河谷地區物種多樣性的深入了解提供科學依據，還可為干熱河谷的植被恢復提供依據。

1 材料與方法

1.1 研究區概況

元謀干熱河谷位于金沙江流域干熱亞區（何毓榮 等 ， 1999 ） ， 地 理 位 置 25°31′～26°07′N ，101°36′～102°07′E，是金沙江干熱河谷的典型區域。屬南亞熱帶干熱季風氣候。隨海拔高度升高及順河而上氣溫下降，濕度和降水量增加以及蒸發量減小，輻射增強。谷底干熱及沿海拔梯度的水熱急劇變化是干旱河谷生態環境的基本特點之一（周志瓊等，2008）。年均降雨量為642.2 mm，年平均蒸發量為 3911.2 mm，蒸發量大約是降雨量的6.2倍，表現了干熱河谷干旱的特征；多年平均氣溫為21.9 ℃，≥10 ℃年積溫為8003 ℃，年日照時數約為 2670.4 h（崔永忠等，2006）。其氣候與稀樹草原（Savanna）氣候相似（Ratnam et al.，2016）?；鶐寥罏樵锛t土，另有變性土、薄層土及紫色土。自然植被以坡柳（Dodonaea viscosa）+扭黃茅（Heteropogon contortus）群系為主（金振洲等，1987）。主要喬木樹種有云南松（Pinus yunnanensis）、余甘子（Phyllanthus emblica）、滇刺棗（Ziziphus mauritiana）；主要灌木有坡柳、疏序黃荊（Vitex negundo）、小石積（Osteomeles anthyllidifolia）、鐵仔（Myrsine africana）、金合歡（Acacia farnesiana）、清香木（Pistacia weinmannifolia）；主要草本物種有扭黃茅、孔穎草（Bothriochloa pertusa）。

1.2 樣地設置

本研究采用樣方調查法，于元謀干熱河谷的2017年旱季時期（2017年4月）共設置了10個海拔梯度。從海拔1100 m左右開始，每隔100 m海拔根據地形條件設置4個10 m×10 m的樣方。在10 m×10 m的大樣方內再設置3個1 m×1 m的小樣方內調查草本層。共設置40個10 m×10 m的樣方，120個 1 m×1 m的小樣方。喬木層記錄所有胸徑（DBH）≥1 cm的物種及株數、高度、胸徑和蓋度，灌木層及草本層記錄所有種的株數、高度和蓋度。胸徑用圍尺測量，高度和蓋度采用目測法。其中，總蓋度為10 m×10 m樣方內所有植物垂直投影面積占樣方面積的比例，灌木層蓋度為10 m×10 m樣方內灌木植物垂直投影所占的比例，草本層蓋度為10 m×10 m樣方內草本植物垂直投影所占的比例，分層蓋度加和值大于或等于總蓋度。同一海拔4個樣地的蓋度為4個樣方的平均值。

1.3 數據分析

1.3.1 數據處理

運用Office WPS（北京金山辦公軟件股份有限公司和珠海金山辦公軟件有限公司）對實驗數據進行初步整理與計算。運用SPSS 13.0（Chicago，IL，USA）對海拔與物種蓋度、物種多樣性進行Pearson相關性分析。

1.3.2 α多樣性計算

本研究采用了4種常用的物種多樣性指數（姜雋，2009；劉雨芳，2003；馬克平等，1995；王飛等，2017）：豐富度指數（S）、Shannon-Wiener物種多樣性指數（H′）、Pielou群落均勻度指數（J）以及Simpson生態優勢度指數（D）進行群落特征分析，其計算公式如下：

式中，S為物種數目；N為群落中的物種數；pi為屬于種i的個體占全部個體的比例。

1.3.3 β多樣性計算

該指標是反映不同環境梯下群落之間物種組成的相異性或物種沿環境梯度的更替速率，表示生境的變異程度或指示生境被物種分隔的程度。不同環境梯度中不同樣地之間共有種越少，β多樣性就越高。一般采用相異性指數和Cody指數進行表征，其計算公式如下（孫榮等，2011）：

相異性指數（SI）：

Cody 指數（βc）：

式中，a和b分別為兩樣方的物種數；c為兩樣方的共有物種數；g(H)是沿生境梯度H增加的物種數目；l(H)是沿生境梯度 H失去的物種數目，即在上一個梯度中存在而在下一個梯度中沒有的物種數目；相異性指數（SI）反映群落或樣方間物種組成的相似性；Cody指數反映物種沿梯度的替代速率。

2 結果與分析

2.1 植被蓋度的變化趨勢

圖1 蓋度隨海拔的變化趨勢Fig. 1 Tendency of coverage to change with elevation

由圖1可知，隨著海拔的升高，群落總蓋度與灌木層蓋度整體上均呈現出了隨之增高的趨勢。植被總蓋度整體波動相比灌木層的蓋度波動大。在海拔1400～1500 m處植被總蓋度出現一個急劇下降的趨勢，海拔1500 m之后，蓋度又有一個回升的趨勢，且在海拔2000 m時，植被總蓋度達到最大值。灌木層蓋度隨海拔高度逐漸上升，在海拔1100 m時，其蓋度最低，海拔2000 m時，灌木層蓋度同樣達到最高值。草本層蓋度則呈緩慢下降的趨勢，但未有較大的變化。表1顯示，海拔高度與群落總蓋度在P＜0.05水平上顯著正相關；海拔高度和灌木層蓋度在P＜0.01水平上呈現極顯著正相關；同時群落總蓋度和灌木層蓋度以及群落總蓋度與草本層蓋度在P＜0.01水平上都存在極顯著正相關。由相關性分析可以得知，海拔與總蓋度存在顯著正相關關系，證明海拔的變化影響了植物蓋度的變化，海拔升高，植物總蓋度增加。分析其原因，河谷海拔越高，其水熱環境越好，越適宜植物生長故植物蓋度越高。而灌木及喬木層蓋度增加，導致了光熱條件的重新配置，使得草本層蓋度有所降低。海拔變化引起群落蓋度的變化，群落總蓋度受灌木層、草本層以及喬木層蓋度的共同作用。

表1 海拔與蓋度相關性分析Table 1 Analysis of the correlation between elevation and coverage

2.2 α物種多樣性的變化趨勢

表2顯示了元謀干熱河谷不同海拔梯度灌木群落的α物種多樣性測定結果。分析可知，Gleason指數（D1）和Marglef指數（D2）變化最大，說明隨著海拔升高，灌木群落物種多樣性變化較大，且總體呈現出隨著海拔升高，物種多樣性增強的趨勢。

表2 元謀干熱河谷不同海拔梯度灌木群落的的a多樣性指數Table 2 A diversity index of shrub communities at different elevation in dry-hot valley of Yuanmou

在海拔 2000 m時，物種豐富度達到最大值（15）；海拔1100～1500 m時，物種豐富度比較穩定，未發生較大變化。自海拔1500 m之后，物種豐富度迅速上升，且隨著海拔的升高呈現出物種豐富度增強的趨勢。Shannon-Wiener指數、Pielou指數隨海拔變化的趨勢基本與豐富度指數一致，在海拔2000 m處物種最為豐富，分布均勻，但Simpson指數的變化趨勢與之相反，整體上隨著海拔的升高呈現出降低的趨勢。

圖2顯示了元謀干熱河谷不同海拔梯度α物種多樣性中的Simpson指數和Shannon-Wiener指數的測定結果。從圖2中可以看出，Shannon-Wiener指數隨海拔梯度整體上呈現出上升的趨勢，Simpson指數的變化趨勢與 Shannon-Wiener指數的變化趨勢大致相反，整體上呈現出隨海拔梯度下降的趨勢。兩個指數隨海拔的變化均不呈現出線性的趨勢，且均在海拔1500 m處有較大的波動。

圖2 元謀干熱河谷不同海拔梯度α多樣性變化趨勢Fig. 2 Variation trend of alpha diversity at different Elevation gradients in dry-hot valley of Yuanmou

2.3 β物種多樣性的變化趨勢

元謀干熱河谷不同海拔梯度的β多樣性如圖3所示。在沿海拔梯度上，植物群落SI指數隨海拔的增加總體上呈現出“M”型變化趨勢。在海拔1400～1500 m之間的相似系數最大，為0.8。而在海拔 1100～1200、1200～1300、1300～1400 m之間的相似系數最小，僅為0.4。Cody指數與SI指數隨海拔的變化趨勢大致相反，在海拔1400～1500 m之間Cody指數達到最小值（0.5），且隨海拔的增加總體上呈現出先下降后上升的趨勢?？梢?，隨著海拔的增加，植物群落間的物種更替速率呈先降低再加快，最終達到最高點的特點。

圖3 元謀干熱河谷不同海拔梯度β多樣性變化趨勢Fig. 3 Variation trend of beta diversity at different elevation gradients in dry-hot valley of Yuanmou

3 討論

本研究聚焦于元謀海拔梯度上干熱河谷植物的物種多樣性?？傮w上，干熱河谷植物群落結構簡單，物種組成稀少，10個海拔梯度下的河谷群落植物只有4個群落的物種組成大于5種，這與該地區干旱惡劣的生境條件有關。本研究中，植被總蓋度與灌木層蓋度總體上呈現出隨著海拔的升高，蓋度增強的趨勢，但兩者的波動范圍不一致，與王飛（2018）關于白龍江干旱河谷不同坡向主要灌叢群落隨海拔梯度變化的物種多樣性研究研究結果不一致，造成這種差異的可能原因是元謀干熱河谷地區隨著海拔的升高，水分越充足（周志瓊等，2008），且灌木層物種趨于穩定，以坡柳為主，而木本植物與草本植物兩者的變化相對較大。

植物群落多樣性指數可以客觀反映群落內物種組成特征，在元謀所有海拔梯度下，植物群落 α多樣性指數偏低，物種豐富度指數、Shannon-Weiner多樣性指數和 Simpson優勢度指數最大值分別為15、1.06、0.09，植物群落物種多樣性與干熱河谷區植被物種多樣性水平相當（王飛等，2018；張建利等，2010；楊振寅等，2008）。Shannon-Wiener指數與Simpson指數之間的關系沿海拔呈現出相反的趨勢，Shannon-Wiener指數與海拔之間表現出較強的正相關關系，雖然受水熱條件以及人為干擾等因素的影響出現波動性，但總體呈現出隨著海拔升高，物種多樣性增強的趨勢（表 2，圖 2），與楊陽（2016）關于三江并流地區干旱河谷植物物種多樣性海拔梯度格局的比較研究結果，張大才等（2018）關于高寒草甸濕潤與干旱化生境中植物多樣性分布與植被的關系中的研究結果基本一致。

β多樣性用于反映環境梯度下植物群落間的差異性，當一個群落向另一個群落過渡時，環境梯度變化較大，不同群落或樣方間共有種越少，β多樣性越大（王慶鎖等，2000；張曉龍等，2017；韓路等，2017）。本研究中，Cody指數與SI指數之間的變化趨勢基本相反。隨著海拔的增加，SI相似性系數在海拔 1400～1500 m 的區域達到最大值（0.8），而在 1100～1200、1200～1300、1300～1400 m區域的相似系數最小，僅為0.4。同時，Cody指數在海拔1400～1500 m的區域值最小，為0.5，說明在該海拔區域的 β多樣性最小，即在海拔1400～1500 m區域的植物群落結構最穩定，物種更替速率最小。

1400～1500 m總蓋度最低，植被稀疏，物種豐富度偏低，推測該海拔是區域內環境條件最為嚴苛的地段，其穩定性亦說明，只有最適于當地環境的物種得以生存，其他物種難以進入。所以，這個地段的適生物種可以用于當地植被難恢復地區的生態建設。與楊銳等（2015）及宋璐璐等（2011）關于在海拔梯度上各種環境因子會表現出連續的梯度性變化的研究結果不一致，造成這些變化的主要原因可能與該地區的物種替代有關，即在海拔1400～1500 m處具有邊緣效應，是干熱河谷植被與常綠闊葉林植被的過渡帶。同時該研究結果與歐曉昆（1988）關于海拔1600 m是干熱河谷植被與常綠闊葉林植被的劃分帶的結論不一致，這可能是由于人為活動（植樹造林、放牧、開墾耕地）的影響而導致過渡帶下移。

4 結論

研究表明，元謀干熱河谷植物多樣性指數、植被蓋度與海拔之間均并不呈現連續的線性的變化趨勢，且都在海拔 1400～1500 m處出現明顯的波動，表明這個區段是生態恢復關鍵區段，這里的適生物種最適宜在干熱河谷環境中生存，是干熱河谷中的建群種。