青藏鐵路沿線植物群落β多樣性研究

馬濤 周金星

摘要 [目的]分析青藏鐵路沿線植物群落β多樣性，了解該地植物群落物種多樣性沿環境梯度的變化格局。[方法]應用β多樣性測度公式對青藏鐵路沿線植物群落多樣性特征進行定量分析。[結果]鐵路沿線各樣帶βW指數呈緩慢下降趨勢。最大值出現在第9樣帶，為4.135；最小值出現在24樣帶，為1.000。1號樣帶和其他各樣帶間βC和βT指數變化趨勢基本一致。βC、βT指數最大值分別出現在1號-12號樣帶間和1號-9號樣帶間，分別為44.5、4.118；最小值分別出現在1號-2號樣帶間和1號-24號樣帶間，分別為7.000、1.100。βCJ和βCS指數在1號樣帶與2、5、17號樣帶間的指數值為0.875、0.978、0.984以及0.778、0.935、0.968，與其他各樣帶間的βCI、βCS指數均為1.000。各相鄰樣帶間βC和βT指數變化趨勢基本一致。βC和βT指數最大值出現在24、25號樣帶間和2、3號樣帶間，分別為43.000、3.542；最小值出現在1、2號樣帶間和20、21樣帶間，分別為7.000、0.678。βCJ、βCS指數最大值均出現在2、3樣帶間，為1.000；最小值均出現在20、21號樣帶間，分別為0.440、0.333。[結論]研究結果為該區植被管理和生物多樣性保護提供了理論依據。

關鍵詞β多樣性；植物群落；青藏鐵路

中圖分類號Q948文獻標識碼A文章編號0517-6611（2017）16-0016-05

βdiversity of Plant Communities along the Line of QinghaiTibetan Railway

MA Tao1，2，ZHOU Jinxing3，4*，LIU Yuguo2（1.Soil and Water Conservation Institute of Gansu，Lanzhou，Gansu 730020；2.Institute of Desertification Studies，Chinese Academy of Forestry，Beijing 100091；3.Jianshui National Field Station，Key Laboratory of State Forestry Administration on Soil and Water Conservation，Beijing Forestry University，Beijing 100083；4.Engineering Research Centre of Forestry Ecological Engineering，Ministry of Education，Beijing Forestry University，Beijing 100083）

Abstract[Objective] The aim was to analyze βdiversity of plant communities along the line of QinghaiTibetan railway and understand change pattern of plant communities along the line of QinghaiTibetan railway. [Method] The plant communities diversity along the line of QinghaiTibetan railway had been quantitatively analyzed by using βdiversity index formulas. [Result] The βW index of all transect was descend along the railway. The max. was 4.135 at 9 transect and the min. was 1.000 at 24 transect. The βC index and βT index between the 1 transect and other transects had the same diversification trend. The max. of βC and βT was 44.500 and 4.188 between 1 and 12 transect and between 1 and 9 transect. The min. was 7.000 and 1.100 between 1 and 2 transect and between 1 and 24 transect. The value of βCJ and βCS was 0.875， 0.978， 0.984 as well as 0.778， 0.935， 0.968 between 1 and 2， 5， 17 transects， and the other values of βCJ and βCS was 1.000. The βC and βT index at any every two transects had similar diversification trend. The max. of βC and βT was 43.000 and 3.542 between 24 and 25 transect and between 2 and 3 transect. The min. was 7.000 and 0.678 between 1 and 2 transect and between 20 and 21 transect. The max. of βCJ and βCS was 1.000 all between 2 and 3 transect. The min. was 0.440 and 0.333 all between 20 and 21 transect. [Conclusion] The results provide theoretical basis for vegetation management and biodiversity conservation of QinghaiTibetan railway.

Key wordsβdiversity；Plant communities；QinghaiTibetan railway

β多樣性是群落多樣性研究的重要內容，可表征物種沿某一環境梯度的替代程度或速率、物種周轉速率和生物變化速率[1-2]，反映群落內或群落間環境異質性的大小及其對物種多樣性的影響[3]，也可為分析某一擾動因素及其強度對群落多樣性的干擾程度提供有價值的信息[4]。測度群落β多樣性的意義在于：它可以反映生境變化的程度或指示生境被物種分隔的程度；β多樣性的高低可以用來比較不同地點的生境多樣性；β多樣性與α多樣性一起構成了群落或生態系統總體多樣性或一定地段的生物異質性[5-7]。

青藏高原地區植物多樣性研究多以特定區域或者特定植被類型為主[8-11]，而青藏鐵路僅就沿線某一區段內的植物多樣性有研究報道[12-13]，但目前關于鐵路沿線植物群落β多樣性的研究鮮見報道。鑒于此，筆者應用β多樣性測度公式對青藏鐵路沿線植物群落多樣性特征進行了定量分析，以期為該區植被管理、生物多樣性保護乃至地區生態系統及景觀多樣性的持續提供理論依據。

1研究區概況與方法

1.1研究區概況

青藏鐵路（格爾木—拉薩）位于青藏高原腹地，地處90°30′～95°56′ E，29°48′～36°32′ N。青藏鐵路北起青海省格爾木市，經青海省西大灘、昆侖山口、五道梁、沱沱河，翻越唐古拉山，到西藏自治區安多、那曲、當雄、羊八井，終點位于拉薩市，全長1 142 km。其中，青海省境內562 km，西藏自治區境內548 km，海拔≥4 000 m地段960 km以上，最高點海拔5 072 m，位于唐古拉山口，多年連續凍土地段約550 km。自昆侖山至唐古拉山約450 km的路段處于長江源地區，橫跨我國三江源自然保護區和可可西里自然保護區；唐古拉山以南進入我國一江二河自然保護區，分布有大片濕地。因此，保護好青藏鐵路沿線生態環境是一項重要的任務，也是鐵路環境保護工作面臨的巨大挑戰。

根據地形地貌、植被類型等指標將青藏鐵路沿線劃分為5個生態類型區[14]，基本概況見表1。

1.2研究方法

1.2.1樣地設計與調查。

根據植物分布類型復雜程度將鐵路沿線分為格爾木—當雄和當雄—拉薩2段進行調查。各樣帶以鐵路線為中心，垂直于鐵路線布置。格爾木—當雄段長約980 km，每隔50 km設1條樣帶，共設置21條樣帶；當雄—拉薩段長約160 km，每隔25 km設1條樣帶，共設置6條樣帶。每條樣帶在鐵路兩側各設置樣地5個，樣地面積400 m2，樣地中心距鐵路線距離分別為5、25、75、150和300 m。5 m處樣地規格為10 m×40 m，其余為20 m×20 m。每個樣地按4頂點、對角線上4點及對角線交點設置9個樣方調查植被，樣方面積1 m×1 m，記錄植物種類、高度、蓋度、

株叢數及物候期，同時測量各樣地海拔、經緯度等因子。調

查工作于2005年6—8月完成。

1.2.2β多樣性的測度。

關于β多樣性測度，馬克平等[7]進行了比較系統的評述。該研究根據調查數據，并借鑒α多樣性測度的經驗，選用以下指數。

1.2.2.1Whittaker提出的β多樣性指數[15]。

βW=S/

式中，S為研究系統記載的所有物種數目；為各樣方或樣本的平均物種數。

1.2.2.2Cody β多樣性指數[16]。

βC=[g（H）+l（H）]/2

式中，g（H）是沿生境梯度H增加的物種數目；l（H）是沿生境梯度H失去的物種數目，即在上一個梯度中存在而在下一個梯度中沒有的物種數目。

1.2.2.3

Wilson-Schmida指數βT[17] 。該指數是把Cody指數（βC）與Whittaker指數（βW）結合形成的。

βT=[g（H）+l（H）]/2

式中，變量的含義與βC、βW指數計算公式中相應變量的含義相同。

1.2.2.4群落相似性系數及其變形。運用相似性系數測度群落或生境間的相似程度是植物群落研究常用的手段，應用較廣、效果較好的相似性系數是早期提出的Jaccard指數和Sorenson指數[18-19]。

Jaccard指數CJ=j/（a+b-j）

Sorenson指數CS=2j/（a+b）

式中，j為2個群落或樣地共有種數；a和b分別為樣地A和樣地B的物種數。

βCJ=1-CJ

βCS=1-CS

2結果與分析

2.1青藏鐵路沿線各樣帶植被調查結果

鐵路沿線植被分布由小灌木向草本轉化，各樣帶分布植物總數呈先上升后趨于平穩之勢（表2）。最大值出現在12號樣帶，為81種。從1號樣帶南山口到12號樣帶唐古拉山北，植物種以矮火絨草（Leontopodium nanum）、嵩草（Kobresia spp.）、針茅（Stipa spp.）、羊茅（Festuca spp.）及風毛菊（Sanguisorbo spp.）為主；13號樣帶扎加藏布到27號樣帶東嘎植物種以矮火絨草、苔草（Carex spp.）、嵩草、委陵菜（Potentilla spp.）以及墊狀植物為主。在20號樣帶當雄車站出現狼毒（Stellera chamaejasme），說明此處植物開始有退化趨勢。不論是哪個樣帶，優勢種蓋度均在30%以上，最高達80%，出現在2號樣帶納赤臺。

青藏高原屬高海拔、低氣壓、氣候寒冷的高原氣候，其腹地上的鐵路沿線植物分布格局是多個生態過程的產物，這些生態過程主要受物種進化、地理差異以及環境因子的控制。氣候條件決定了生物所能獲得的光、熱、水等生長必需條件，對物種空間分布起重要作用。相比年均氣溫、年降水量，鐵路沿線海拔高度沿各生態分區的變化最明顯，因此，海拔高度決定了物種的分布格局。

2.2各樣帶植物群落β多樣性分析

β多樣性指標測度是多樣性測度中一個非常重要的手段，其結果能夠反映群落沿某一環境梯度其物種組成的差異。為了比較不同樣帶間β多樣性的變化速率，該研究選用5個指標計算鐵路沿線27個樣帶的植物群落β多樣性。其中，Cody指數βC、Wilson和Shmida指數βT、群落相似性指數βCJ及βCS選擇1號樣地與其他各樣地之間，以及每2個相鄰樣地間進行計算。將各指數的計算結果繪于圖1～3。鐵路沿線海拔高度變化繪于圖4。

Whittacker指數βW能夠直觀地反映β多樣性與物種豐富度之間的關系。由圖1可知，βW指數呈緩慢下降趨勢，最大值出現在9樣帶，為4.135；最小值出現在24樣帶，為1.000。其中在2、9、24號樣帶值出現較大波動，分別為2.390、4.135、1.000。原因是：與1號樣帶相比，2號樣帶物種豐富度S雖有所增加，但平均物種數在2號樣帶比1號樣帶大，所以其βW值比1號樣地低。9號樣帶物種豐富度（55）比相鄰的8號樣帶（34）、10號樣帶（49）大，并且各樣帶差別較大，所以在9號樣帶βW指數出現較大波動。由于24號樣帶所在位置特殊，位于羊八井2號隧道入口處，坡度較大，并且僅在鐵路一側分布有零星植被，故該樣帶與S值相等，βW指數為1。另外，將βW曲線與海拔曲線比較可知，其變化趨勢與各樣帶海拔高度變化一致，只在第1樣帶有所不同，原因是在1號樣帶分布植物種數少（8種），值為2，βW指數較大（4.000）。

隨海拔梯度的變化，1號樣帶和其他各樣帶間Cody指數βC及由βC與βW結合形成的Wilson-Schmida指數βT有較為相同的變化趨勢（圖2）。βC和βT指數的最大值分別出現在1號-12號樣帶間和1號-9號樣帶間，分別為44.500和4.118；最小值分別出現在1號-2號樣帶間和1號-24號樣帶間，分別為7.000和1.100。各指數值的大小反映群落間物種替代速率陡慢變化懸殊的程度。值越大，共有種越少，即替代速率越大。基于群落相似系數CJ和CS導出的βCJ和βCS指數在1號樣帶與2、5、17號樣帶間的指數值分別為0.875、0.978、0.984以及0.778、0.935、0.968，而與其他各樣帶間的βCJ、βCS指數均為1.000，原因為從1號樣帶開始，植物種從小灌木向草本植物轉化，共有種數j為0，故βCJ、βCS指數均為1.000。

由圖3可知，各相鄰樣帶間βC和βT指數有非常相似的變化趨勢，反映了這2個多樣性指數具有較好的可加性。βC和βT指數最大值出現在24、25號樣帶間和2、3號樣帶間，分別為43.000、3.542；最小值出現在1、2號樣帶間和20、21號樣帶間，分別為7.000、0.678。βC、βT指數最大值均出現在2、3號樣帶間，為1.000，原因為植物種從小灌木向草本植物轉化導致；最小值均出現在20、21號樣帶間，分別為0.440、0.333。其值越大，表明兩群落間相似性越小，共有種越少，物種周轉速率越快。

以上各β多樣性指數的變化情況較好地反映了青藏鐵路沿線植物群落多樣性特征。當然，人類活動（放牧、旅游等）也會對這一地區的植被分布產生較大的影響，從而引起各樣帶間多樣性指數的變化；同時，生境梯度的變化也并非是均勻的，加之抽樣調查的誤差等都會在一定程度上對β多樣性的測度結果產生一定的影響。

3結論

青藏鐵路沿線各樣帶植被分布由小灌木向草本轉化。1號樣帶南山口到12號樣帶唐古拉山北，植物種以矮火絨草、嵩草、針茅、羊茅和風毛菊（Sanguisorbo spp.）為主；13號樣帶扎加藏布到27號樣帶東嘎植物種以矮火絨草、苔草、嵩草、委陵菜和墊狀植物為主。在20號樣帶當雄車站出現植被退化指示性植物狼毒。

鐵路沿線各樣帶指數βW呈緩慢下降趨勢，與海拔高度變化情況較為一致。最大值出現在第9號樣帶，為4.135；最小值出現在24號樣帶，為1.000。其中，在2、9、24號樣帶βW值出現較大變幅，其值分別為2.390、4.135、1.000。

1號樣帶和其他各樣帶間βC及βT指數變化趨勢基本一致。βC、βT指數最大值分別出現在1號-12號樣帶間和1號-9號樣帶間，指數值分別為44.500、4.118；最小值分別出現在1號-2號樣帶間和1號-24號樣帶間，指數值分別為7.000、1.100。βCJ和βCS指數在1號樣帶與2、5、17號樣帶間分別為0.875、0.978、0.984以及0.778、0.935、0.968，與其他各樣帶間的βCJ、βCS指數均為1.000。

各相鄰樣帶間βC和βT指數變化趨勢相似，反映了這2種指數具有較好的可加性。βC和βT指數最大值出現在24、25號樣帶間和2、3號樣帶間，分別為43.000、3.542；最小值出現在1、2號樣帶間和20、21號樣帶間，分別為7.000、0.678。βCJ、βCS指數最大值均出現在2、3號樣帶間，為1.000；最小值均出現在20、21號樣帶間，分別為0.440、0.333。

植物群落β多樣性指數的大小反映物種周轉速率的快慢。生境梯度、海拔、人類放牧及旅游活動等均是其影響因子。

參考文獻

[1] 郝占慶，于德永，吳鋼，等.長白山北坡植物群落β多樣性分析[J].生態學報，2001，21（12）：2018-2022.

[2] 白永飛，邢雪榮，許志信，等.內蒙古高原針茅草原群落β多樣性研究[J].應用生態學報，2000，11（3）：408-412.

[3] 陳世品，游水生，陳子英，等.運用β多樣性探討不同更新模式對米櫧群落高度級和物種多樣性的影響[J].植物資源與環境學報，2013，22（4）：62-69.

[4] 郝智如，侯瑞萍，張克斌，等.寧夏鹽池縣人工封育區植物群落的β多樣性[J].水土保持通報，2014，34（3）：257-261，299.

[5] 高賢明，馬克平，黃建輝，等.北京東靈山地區植物群落多樣性的研究Ⅺ.山地草甸β多樣性[J].生態學報，1998，18（1）：24-32.

[6] 林開敏，黃寶龍.杉木人工林林下植物物種β多樣性的研究[J].生物多樣性，2001，9（2）：157-161.

[7] 馬克平，劉燦然，劉玉明.生物群落多樣性的測度方法Ⅱ β多樣性的測度方法[J].生物多樣性，1995，3（1）：38-43.

[8] 羅亞勇，孟慶濤，張靜輝，等.青藏高原東緣高寒草甸退化過程中植物群落物種多樣性、生產力與土壤特性的關系[J].冰川凍土，2014，36（5）：1298-1305.

[9] 董世魁，湯琳，王學霞，等.青藏高原高寒草地植物多樣性測定的最小樣地面積[J].生物多樣性，2013，21（6）：651-657.

[10] 許岳飛，益西措姆，付娟娟，等.青藏高原高山嵩草草甸植物多樣性和土壤養分對放牧的響應機制[J].草地學報，2012，20（6）：1026-1032.

[11] 李中強，張萌，徐軍.青藏高原納木錯水生植物多樣性及群落生態學研究[J].水生生物學報，2011，32（1）：1-13.

[12]

殷代英，韓慶杰，屈建軍，等.青藏鐵路格爾木-錯那湖段沿線物種多樣性分析[J].干旱區資源與環境，2013，27（12）：152-160.

[13] 周金星，YANG J，董林水，等.青藏鐵路唐古拉山南段沿線植被多樣性及蓋度特征分析[J].北京林業大學學報，2008，30（3）：24-30.

[14] 馬濤，周金星.青藏鐵路沿線植物群落物種多樣性研究[J].水土保持通報，2012，32（2）：261-266.

[15] WHITTAKER R H.Vegetation of the Siskiyou Mountains[J].Ecol Monogr，1960，30：279-338.

[16] CODY M L.Towards a theory of continental species diversity bird distributions over Mediterranean habitat gradients[M]//CODY M L，DIAMOND J M.Ecology and Evolution of Communities.Cambridge：Harvard University Press，1975：214-257.

[17] WILSON M V，SCHMIDA A.Measuring beta diversity with presenceabsence data[J].J Ecol，1984，72：1055-1064.

[18] WHITAKER R H.Evolution and measurement of species diversity[J].Taxon，1972，21（2/3）：213-251.

[19] SOUTHWOOD T R.生態學研究方法：適用于昆蟲種群的研究[M].羅河清，等，譯.北京：科學出版社，1984.