王 建，王光鵬，斯貴才，沈妙根，張更新

(1.中國科學院青藏高原研究所，北京 100101； 2.中國科學院大學，北京100049；3.中國科學院地質與地球物理研究所蘭州油氣資源研究中心，甘肅 蘭州 730000)

生物多樣性與全球變化是目前各國面臨的主要生態(tài)環(huán)境問題[1]。生物多樣性的中心是物種多樣性[2]，物種多樣性是一個群落結構和功能復雜程度的量度[3]，它包括物種豐富程度、單位面積物種數(shù)目(物種密度)、特有物種比例等[4]，是當前群落生態(tài)學中重要研究內(nèi)容之一。研究物種多樣性對認識及保護生物多樣性具有重要意義。物種多樣性的研究主要集中在多樣性隨環(huán)境梯度的變化情況[5-7]等方面。大尺度水平上主要表現(xiàn)為緯度梯度變異，小尺度水平上主要表現(xiàn)為海拔梯度的變異。多數(shù)理論研究認為，隨緯度的增加物種多樣性下降[8]，形成這一格局的原因主要歸結于生產(chǎn)力或其它與能量和環(huán)境有關的因素[9]。海拔梯度由于涉及溫度、降水和太陽輻射等諸多環(huán)境因子而成為生物多樣性梯度格局研究的又一個重要方面[10]。此外，近年來，小尺度水平上有關單個環(huán)境因子對物種多樣性影響的研究也不斷增多，在青藏高原主要表現(xiàn)為模擬控制試驗尤其是增溫和放牧試驗[11-15]。但是，有關該區(qū)物種多樣性中起決定性因素的環(huán)境因子的研究并不多見。

草地生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的重要組成部分。目前，在生物多樣性的生態(tài)系統(tǒng)功能研究中，草地生態(tài)系統(tǒng)備受關注[16]。青藏公路地處世界屋脊，具有獨特的自然地理環(huán)境和草原草甸群落特征。研究沿線植物物種多樣性分布格局及其與環(huán)境因子的關系是人們關注的問題。目前，有關青藏公路沿線植被特征調查僅就特定的植被類型區(qū)做過一些植被多樣性的研究[17-20]，而對于完整的青藏公路沿線植物群落多樣性的研究還未有系統(tǒng)報道。為了分析青藏公路沿線植物物種多樣性分布格局及其與環(huán)境因子的關系，本研究對青藏公路沿線兩側植物進行樣方調查分析，以期為研究該區(qū)生物多樣性的保護與可持續(xù)利用提供科學參考。

1 研究方法

1.1研究區(qū)概況 青藏公路位于號稱地球“第三極”的青藏高原，全長1 917 km，全線平均海拔在4 000 m以上，是世界上海拔最高的公路。沿線降水量較少，年均溫較低，屬于干旱寒冷的典型生態(tài)敏感區(qū)。該區(qū)干濕兩季分明，冬長夏短，季節(jié)變化緩慢，太陽輻射強烈，日照時間長，年溫差小，日溫差大，受局部地形及小環(huán)境影響顯著，不同地域溫度差異很大[21]。青藏公路所跨范圍可以分為6個生態(tài)分區(qū)[18,22-23]，即西寧-天俊縣山地草原區(qū)，天俊縣-南山口山地荒漠區(qū)，南山口-昆侖山山地荒漠、山地草原區(qū)，昆侖山-唐古拉山高寒草甸、草原、荒漠區(qū)，唐古拉山-當雄高寒草甸區(qū)，當雄-拉薩灌叢草原區(qū)。草原草甸優(yōu)勢種主要有高山嵩草、藏嵩草、矮嵩草、紫花針茅和青藏苔草等，亞優(yōu)勢種主要有早熟禾、黃芪和委陵菜等，伴生種主要有冷地早熟禾、弱小火絨草、風毛菊和棘豆等[24-27]。

表1 樣點概況Table 1 Generalization of transect in different latitude gradients of Tibetan Plateau

1.2樣地設置 在植物生長最旺盛的8月，從西寧市出發(fā)，沿青藏公路進行植物群落調查，大約每隔50 km選擇一處樣地，一直到拉薩市共調查29處樣地(圖1)。在每處樣地選取一個有代表性的10 m×10 m的大樣方，在大樣方中隨機設置3個1 m×1 m的小樣方，進行樣方調查。樣方調查采用針刺法，調查記錄每個小樣方內(nèi)出現(xiàn)的物種數(shù)及每個物種的蓋度、高度、密度和頻度，同時，采集植物標本并記錄每個樣地所在位置的經(jīng)緯度、海拔等信息，然后在室內(nèi)準確鑒定植物標本，統(tǒng)計科屬種及其組成。各樣點經(jīng)緯度、海拔高度及主要植物種詳見表1。

圖1 樣地分布圖Fig.1 Location of the sample sites

1.3樣地多樣性計算 有關物種多樣性的計算方法很多，本研究根據(jù)野外蓋度、物種數(shù)、株數(shù)和頻數(shù)等基礎樣方數(shù)據(jù)，運用Excel進行數(shù)據(jù)的初步處理，同時采用Sigma Plot 10.0、ArcGIS 9.3軟件制圖及SPSS 19.0統(tǒng)計分析軟件進行相關分析和逐步回歸等統(tǒng)計分析。物種多樣性用豐富度指數(shù)、多樣性指數(shù)和均勻度指數(shù)3類指數(shù)表征[28-32]，其計算方法如下：

重要值(IV)=(相對頻度+相對蓋度)/2。

豐富度指數(shù)：R=S.

多樣性指數(shù)：

Simpson指數(shù)：D=1-∑Pi2；Pi=Ni/N.

Shannon-Wiener指數(shù)：

H=-∑Pi·lnPi.

均勻度指數(shù)：

Pielou指數(shù)：P=H/lnS.

Alatalo指數(shù)：A=[1/(D-1)]/(eH-1).

式中，S為調查樣方的物種數(shù)，由于個體數(shù)對于草本植物統(tǒng)計較為困難，本文采用各物種的重要值(IV)代替?zhèn)€體數(shù)進行計算[33-34]，Ni為樣方中第i種植物的重要值，N為樣方中所有物種的重要值之和。

1.4氣候數(shù)據(jù)來源和處理 對全國700多個標準氣象臺站1979-2010年間每月平均溫度和降水數(shù)據(jù)進行克里格插值分析，然后根據(jù)樣地的經(jīng)緯度從插值結果中提取樣地每月平均溫度、降水量、空氣比濕、短波輻射、長波輻射等數(shù)據(jù)，以各樣點32年的氣象數(shù)據(jù)平均值作為樣地的溫度和降水等數(shù)據(jù)。數(shù)據(jù)來源詳見Yang等[35]和Chen等[36]。

2 結果與分析

2.1青藏公路沿線各樣帶植物群落物種多樣性分析 由于青藏公路沿線氣候、地形等條件的差異，沿線29個樣點的物種多樣性之間存在一定的差異。

2.1.1各樣點總蓋度和豐富度特征 沿線各樣點1 m2樣方平均植被總蓋度和物種豐富度沿青藏公路呈先降低后升高然后趨于平穩(wěn)的趨勢(圖2)，而植被總蓋度這種變化趨勢更加明顯。

青海湖南岸1～4號樣點植被蓋度為90%左右，5號樣點由于靠近柴達木盆地東南緣，植被總蓋度降至青海湖的1/3，而物種豐富度較1～4號樣點有所減少。10～21號樣點植被總蓋度開始呈增加趨勢，10～12號3個樣點增加緩慢，14～16號3個樣點增加迅速，17～21號5個樣點增加緩慢，而物種豐富度最低點出現(xiàn)在海拔較高的風火山隧道(14號)。22號樣點植被總蓋度較17～21號5個樣點降低較明顯，物種豐富度略低于17～21號5個樣點。那曲地區(qū)3個樣點(23～25號)植被總蓋度較17～21號5個樣點降低較明顯，較22號樣點沒有明顯差異，物種豐富度沒有明顯變化趨勢。進入當雄灌叢草原區(qū)，樣點(26～29號)物種豐富度呈下降趨勢，而草原植被總蓋度并沒有受到灌叢的影響，反而較那曲地區(qū)有所增加。其中那賽村樣點(28號)植被總蓋度偏低，可能由于該樣點海拔(4 593 m)較其它3個樣點高所致(圖2)。

2.1.2各樣點多樣性指數(shù)特征 除個別樣點波動幅度較大外，物種多樣性辛普森指數(shù)和香農(nóng)指數(shù)總體從1號樣點到29號樣點呈降低趨勢(圖3)。

圖2 青藏公路沿線植被總蓋度和物種豐富度的變化Fig.2 Vegetation coverage and species richness along Qinghai-Tibet Highway

圖3 青藏公路沿線物種多樣性指數(shù)的變化Fig.3 Species diversity index along Qinghai-Tibet Highway

沿線樣點辛普森指數(shù)介于0.38～0.81，總體呈降低趨勢。1號樣點辛普森指數(shù)為0.75，一直到12號樣點辛普森指數(shù)緩慢降低，中間樣點辛普森指數(shù)略有波動。在13號和14號樣點辛普森指數(shù)急速下降至0.40左右，之后在15號樣點約增加至14號樣點的兩倍，然后一直到22號樣點呈現(xiàn)出遞減的趨勢。其中，15號到21號樣點辛普森指數(shù)緩慢降低，22號樣點降低幅度最大，約降至21號樣點的1/2。到25號樣點辛普森指數(shù)迅速增加至0.70左右，26號樣點緩慢降低，到27號樣點又迅速降低，最后兩個樣點又增加至26號樣點水平。辛普森指數(shù)對于物種豐富度的敏感性較低，而對于物種均勻度的敏感性較高[37]?？傮w來說，13號、14號、22號和27號4個樣點辛普森指數(shù)偏低，主要是由這些樣點物種分布不均勻所致[38]。

香農(nóng)指數(shù)與物種豐富度和分配均勻性有關，物種數(shù)目越多，多樣性越大，物種間個體分配的均勻性越大，多樣性越大。沿線樣點香農(nóng)指數(shù)所呈現(xiàn)的規(guī)律和辛普森指數(shù)基本一致，不同的是13號、14號、22號、23號和27號5個樣點波動幅度較辛普森指數(shù)小，而16號和17號兩個樣點波動幅度較辛普森指數(shù)大。

2.1.3各樣點均勻度指數(shù)特征 Pielou指數(shù)越大，物種分布越均勻，優(yōu)勢種越不明顯，Alatalo指數(shù)則相反。Pielou指數(shù)在1～12號樣點先保持穩(wěn)定趨勢，維持在0.73～0.88，快到風火山山頂時(13號)急速下降至0.60左右，之后增至風火山以前水平。過17號樣點后開始下降，直到安多山(22號)降至最低，之后呈增加趨勢，到25號樣點增至最大。翻越念青唐古拉山(27號)時又出現(xiàn)明顯降低趨勢，之后兩個樣點又呈快速增加趨勢，29號樣點約增至27號的兩倍。Alatalo指數(shù)變化趨勢與Pielou指數(shù)相反，不同的是這種相反趨勢較Pielou指數(shù)更加明顯，主要體現(xiàn)在個別Pielou指數(shù)偏低的樣點其Alatalo指數(shù)偏高，且波動很大，能達到其它樣點的3～4倍(圖4)。

圖4 青藏公路沿線物種均勻度指數(shù)的變化Fig.4 Species evenness index along Qinghai-Tibet Highway

青海湖南岸到柴達木盆地東南緣再到昆侖山區(qū)六十九道班(12號)，均勻度指數(shù)保持穩(wěn)定，波動幅度不大?？赡茉蚴乔嗪：貐^(qū)水熱條件好，植被生長茂盛，而柴達木盆地東南緣地區(qū)氣候干旱，只有少數(shù)耐旱植物生長，表現(xiàn)為植被覆蓋度低，長勢差等特點，兩個地區(qū)植被優(yōu)勢種均不明顯，植被分布較均勻。此外，高海拔地區(qū)物種均勻度指數(shù)波動較大。在風火山調查的14號樣點、安多山調查的22號樣點和念青唐古拉山調查的27號樣點，都表現(xiàn)出植被Pielou指數(shù)隨海拔升高而降低、Alatalo指數(shù)隨海拔升高而升高的趨勢。

2.2植被物種多樣性與環(huán)境因子的關系

2.2.1物種蓋度和豐富度與環(huán)境因子的關系 馬濤等[23]研究青藏鐵路格爾木至拉薩段植被樣方時發(fā)現(xiàn)，物種豐富度與海拔高度之間呈現(xiàn)極顯著線性正相關關系，說明海拔高度對物種多樣性有著重要影響[39]。植被多樣性不僅受海拔的影響，而且受其它環(huán)境因子的影響。青藏高原屬高海拔、空氣稀薄、氣壓低、氣候寒冷的高原氣候，其植被多樣性受高寒氣候條件下多種環(huán)境因子的影響。本研究分析了海拔高度、年降水量、生長季降水量、年均溫、生長季均溫、空氣比濕、生長季空氣比濕、短波輻射、生長季短波輻射、長波輻射、生長季長波輻射、壓強和風速這些環(huán)境因子與青藏公路沿線物種多樣性變化特征之間的關系及影響物種多樣性變化特征的關鍵性環(huán)境因子(表2)。物種總蓋度和物種豐富度受水分影響較大(圖5)。相關性分析表明，物種蓋度與年降水量、生長季降水量分別呈極顯著正相關關系(P<0.01)，與空氣比濕呈顯著正相關關系(P<0.05)。物種蓋度與年降水量、生長季降水量分別呈顯著正相關關系(P<0.05)，而溫度和其它環(huán)境因子與物種蓋度和物種豐富度之間均未達到顯著相關水平。

表2 青藏公路沿線環(huán)境因子與物種多樣性指數(shù)間的相關系數(shù)Table 2 Coefficients of correlation between species diversity and environmental factors

圖5 物種多樣性與部分環(huán)境因子之間的相關關系Fig.5 Relationships between species diversity and some environmental factors

對物種豐富度、蓋度分別與海拔高度、年降水量、生長季降水量、年均溫、生長季均溫、空氣比濕、生長季空氣比濕、短波輻射、生長季短波輻射、長波輻射、生長季長波輻射、壓強和風速進行逐步回歸分析，回歸方程為：

C=349.19+0.25gp-1.60s；

其中gpr=0.89，sr=-0.48，P<0.05；

R=7.55+0.09p-0.08gp-1.073gh；

其中pr=5.22，gpr=-4.41，ghr= -0.58，P<0.05。

式中，C、gp、S、gpr、sr、R、p、gh和ghr分別為總蓋度、生長季降水量、短波副射、生長季降水量標準偏回歸系數(shù)、短波輻射準偏回歸系數(shù)、物種豐富度、降水量、生長季比濕和生長季比濕標準偏回歸系數(shù)。

逐步回歸方程說明植被總蓋度隨生長季降水量增加而增加，隨短波輻射增加而減少。標準偏回歸系數(shù)絕對值大小表明生長季降水量對植被總蓋度的貢獻大于短波輻射。物種豐富度隨降水量的增加而增加，隨生長季降水量和生長季比濕的增加而減少，標準偏回歸系數(shù)絕對值大小表明降水量對植被總蓋度的貢獻大于生長季比濕。

2.2.2物種多樣性指數(shù)與環(huán)境因子的關系 物種多樣性指數(shù)與短波輻射關系密切(圖5)。香農(nóng)指數(shù)和辛普森指數(shù)均與短波輻射呈顯著負相關關系(P<0.05)。物種均勻度指數(shù)受多種環(huán)境因子共同作用。Pielou指數(shù)與海拔、降水量、短波輻射呈極顯著負相關關系(P<0.01)，與壓強呈極顯著正相關關系(P<0.05)。Alatalo指數(shù)與各環(huán)境因子相關性均沒有達到顯著水平。

此外，溫度、長波輻射、風速和物種多樣性指標之間關系均未達到顯著水平。該結果表明在青藏高原地區(qū)，物種多樣性受溫度、長波輻射、風速這3個環(huán)境因子影響較小。

對香農(nóng)指數(shù)、辛普森指數(shù)、Pielou均勻度指數(shù)和Alatalo均勻度指數(shù)分別和海拔高度、年降水量、生長季降水量、年均溫、生長季均溫、空氣比濕、生長季空氣比濕、短波輻射、生長季短波輻射、長波輻射、生長季長波輻射、壓強和風速進行逐步回歸分析，回歸方程為：

H=4.33-0.01s；

其中sr=-0.41，P<0.05；

D=1.99-0.006s；

其中sr=-0.45，P<0.05；

E=0.98-0.001p；

其中pr=-0.67，P<0.05；

A=-2.81+0.01s；

其中sr=-0.35，P<0.10。

式中，H、s、sr、D、E、p、pr和A分別為香農(nóng)指數(shù)、短波輻射、短波輻射標準偏回歸系數(shù)、辛普森指數(shù)、Pielou均勻度指數(shù)、降水量、降水量標準偏回歸系數(shù)和Alatalo均勻度指數(shù)。

回歸方程說明植被物種多樣性香農(nóng)和辛普森指數(shù)隨短波輻射增加而減少，物種Pielou均勻度指數(shù)隨降水量的增加而減少，物種Alatalo均勻度指數(shù)隨短波輻射的增加而增加。

3 討論

從物種多樣性和環(huán)境因子的相關性分析及逐步回歸分析發(fā)現(xiàn)，水分因子是影響該區(qū)物種總蓋度、物種豐富度和物種Pielou均勻度指數(shù)最主要的因子，而短波輻射是影響該區(qū)物種多樣性指數(shù)香農(nóng)指數(shù)和辛普森指數(shù)最主要的因子。這主要是由于青藏高原地區(qū)特殊的高海拔環(huán)境特征導致該區(qū)干旱缺水和太陽短波輻射強烈，其對植物的影響大于包括溫度在內(nèi)的其它環(huán)境因子，成為影響植物生長最主要的環(huán)境因子，進而決定物種多樣性沿青藏公路的變化趨勢。

沈振西等[39]通過野外控制實驗，探討矮嵩草草甸的群落特征對模擬降水的響應，發(fā)現(xiàn)增加降水可以顯著增加物種多樣性。楊元合等[40]通過對青海海北地區(qū)、唐古拉地區(qū)以及藏北地區(qū)的草原草甸物種豐富度研究和環(huán)境因子的相關性和逐步回歸分析也表明，降水是影響物種豐富度最主要的環(huán)境因子。這些研究結果與本研究結果一致，能較好地說明水分對干旱區(qū)物種多樣性的重要影響。有關太陽短波輻射對植物的影響的研究主要集中在對其形態(tài)結構以及生理特征的影響方面[41]，而對物種多樣性的影響方面還較少。但是，研究區(qū)強烈的太陽短波輻射使植物形態(tài)、結構和生理特征發(fā)生變化的同時，必定會改變物種的重要值，進而引起物種多樣性的改變。

4 小結

除個別樣點波動較大外，總體上看，隨緯度的增加，青藏高原植被蓋度呈先降低后升高的趨勢，物種豐富度呈先降低后升高再降低的趨勢，香農(nóng)指數(shù)和辛普森指數(shù)隨樣點呈降低趨勢，Pielou指數(shù)呈先降低后升高的趨勢，Alatalo指數(shù)變化趨勢與Pielou指數(shù)相反。

從物種多樣性和環(huán)境因子的相關性分析及逐步回歸分析發(fā)現(xiàn)，物種多樣性的變化趨勢主要受水分因子和太陽短波輻射因子控制。此外，一個地區(qū)的植被物種多樣性是多種環(huán)境因子綜合作用的產(chǎn)物，除了降水和短波輻射外，其它環(huán)境因子也對研究區(qū)物種多樣性也有一定程度的影響。

[1] 常學向,趙文智,趙愛芬.祁連山區(qū)不同海拔草地群落的物種多樣性[J].應用生態(tài)學報,2004,15(9):1599-1603.

[2] Tracy C R,Brussard P F.Preserving biodiversity-species in landscapes[J].Ecological Applications,1994,4(2):206-207.

[3] 黃忠良,孔國輝,何道泉.鼎湖山植物群落多樣性的研究[J].生態(tài)學報,2000,20(2):193-198.

[4] 周紅章.物種與物種多樣性[J].生物多樣性,2000,8(2):215-226.

[5] 高賢明,陳靈芝,黃建輝,等.秦嶺太白山棄耕地植物群落演替的生態(tài)學研究.I.演替系列的群落基本特征及其生活型譜分析[J].生態(tài)學報,1997,27(3):253-259.

[6] 沈澤昊,方精云,劉增力,等.貢嘎山東坡植被垂直帶譜的物種多樣性格局分析[J].植物生態(tài)學報,2001,25(6):721-732.

[7] 溫遠光,元昌安,李信賢,等.大明山中山植被恢復過程植物物種多樣性的變化[J].植物生態(tài)學報,1998,22(1):33-40.

[8] Brown J H,Lomolino M V.Biogeography,second edition[M].Sunderland,MA:Sinauer,1998.

[9] Currie D J.Energy and large-scale patterns of animal- and plant-species richness[J].American Naturalist,1991,137(1):27-49.

[10] Gaston K J.Global patterns in biodiversity[J].Nature,2000,405:220-226.

[11] 李英年,趙亮,趙新全,等.5年模擬增溫后矮嵩草草甸群落結構及生物量的變化[J].草地學報,2004，12(3):236-239.

[12] 汪詩平,李永宏,王艷芬,等.不同放牧率對內(nèi)蒙古冷嵩草原植物多樣性的影響[J].植物學報,2001,43(1):89-96.

[13] Klein J A,Harte J,Zhao X Q.Experimental warming causes large and rapid species loss,dampened by simulated grazing,on the Tibetan Plateau[J].Ecology Letters,2004,7(12):1170-1179.

[14] Kudo G,Suzuki S.Warming effects on growth,production,and vegetation structure of alpine shrub:A five-year experiment in northern Japan[J].Oecologia,2003,135:280-287.

[15] 周華坤,趙新全,唐艷鴻,等.長期放牧對青藏高原高寒灌叢植被的影響[J].中國草地,2004,26(6):1-11.

[16] Schl?pfer F,Schmid B.Ecosystem effects of biodiversity:A classification of hypotheses and exploration of empirical results[J].Ecological Applications,1999,9(3):893-912.

[17] 張鐿鋰,閻建忠.青藏公路對區(qū)域土地利用和景觀格局的影響——以格爾木至唐古拉山段為例[J].地理學報,2002,57(3):253-266.

[18] 陳輝,李雙成,鄭度.青藏公路鐵路沿線生態(tài)系統(tǒng)特征及道路修建對其影響[J].山地學報,2003,21(5):559-567.

[19] 鄢燕,張錦華,張建國.青藏公路沿線高寒草甸次生群落特征及生態(tài)修復[J].草地學報,2006,14(2):156-159.

[20] 丁明軍,沈振西,張鐿鋰,等.青藏公路與鐵路沿1981年-2001年植被覆蓋變化[J].資源科學,2005,27(5):128-133.

[21] 周金星,Yang J,董林水,等.青藏鐵路唐古拉山南段沿線植被多樣性及蓋度特征分析[J].北京林業(yè)大學學報,2008,30(3):24-30.

[22] 孫士云.青藏鐵路沿線的生態(tài)環(huán)境特點及保護對策[J].鐵道勞動安全衛(wèi)生與環(huán)保,2002,29(3):121-124.

[23] 馬濤,周金星.青藏鐵路沿線植物群落物種多樣性研究[J].水土保持學報,2012,32(2):261-266.

[24] 陳輝,李雙成,鄭度.基于人工神經(jīng)網(wǎng)絡的青藏公路鐵路沿線生態(tài)系統(tǒng)風險研究[J].北京大學學報(自然科學版),2005,41(4):586-593.

[25] 徐憲立,張科利,龐玲,等.青藏公路路堤邊坡產(chǎn)流產(chǎn)沙規(guī)律及影響因素分析[J].地理科學,2006,26(2):211-216.

[26] 汪雙杰,周文錦.青藏公路沿線環(huán)境演化及環(huán)境保護對策[J].公路,2005(1):37-41.

[27] 杜小冰,陳建業(yè),周劍.青藏公路建設對高寒草原生態(tài)的影響及治理對策[J].交通標準化,2011(10):91-94.

[28] 王文穎,王啟基.高寒嵩草草甸退化生態(tài)系統(tǒng)植物群落結構特征及物種多樣性分析[J].草業(yè)學報,2001,10(3):8-14.

[29] 孫儒泳,李慶芬,牛翠娟.基礎生態(tài)學[M].北京:高等教育出版社,2002:136-137.

[30] Margalef R.Perspectives in Ecological Theory[M].Chicago:University of Chicago,1968.

[31] Hill M,Bunce R,Shaw M.Indicator species analysis,a divisive polythetic method of classification,and its application to a survey of native pinewoods in Scotland[J].The Journal of Ecology,1975,63(2):597-613.

[32] Wang C T,Long R J,Wang Q J,etal.Effects of altitude on plant species diversity and productivity in an alpine meadow,Qinghai-Tibetan plateau[J].Australian Journal of Botany,2007,55(2):110-117.

[33] Peet R K.The measurement of species diversity[J].Annual Review of Ecology and Systematics,1974,5:285-307.

[34] 馬克平.生物多樣性的原理和方法[M].北京:中國科學技術出版社,1994:141-165.

[35] Yang K,He J,Tang W J,etal.On downward shortwave and longwave radiations over high altitude regions——Observation and modeling in the Tibetan Plateau[J].Agricultural and Forest Meteorology,2010,150:38-46.

[36] Chen Y Y,Yang K,He J,etal.Improving land surface temperature modeling for dry land of China[J].Journal of Geophysical Research,2011,116,D20104,doi-10.1029/2011JD015921.

[37] 董林水,張旭東,周金星,等.青藏鐵路沿線北段植被物種豐富度及蓋度的動態(tài)變化[J].長江流域資源與環(huán)境,2008,17(4):551-556.

[38] 許晴,張放,許中旗,等.Simpson指數(shù)和Shannon-Wiener 指數(shù)若干特征的分析及“稀釋效應”[J].草業(yè)科學,2011,28(4):527-531.

[39] 沈振西,周興民,陳佐忠,等.高寒矮嵩草草甸植物類群對模擬降水和施氮的響應[J].植物生態(tài)學報,2002,26(3):288-294.

[40] 楊元合,饒勝,胡會峰,等.青藏高原高寒草地植物物種豐富度及其與環(huán)境因子和生物量的關系[J].生物多樣性,2004,12(1):200-205.

[41] 趙新全.高寒草甸生態(tài)系統(tǒng)與全球變化[M].北京：科學出版社,2009:53.