免费文献传递   相关文献

β diversity characteristic of vegetation community on south part of Gurbantunggut Desert and its interpretation

古尔班通古特沙漠南部植物群落β多样性及其解释



全 文 :书古尔班通古特沙漠南部植物群落
β多样性及其解释
赵怀宝1,2,刘彤3,雷加强2,桂东伟2,赵新俊4
(1.中国科学院研究生院,北京100048;2.中国科学院新疆生态与地理研究所,新疆 乌鲁木齐830000;
3.石河子大学生命科学学院,新疆 石河子832003;4.石河子大学师范学院,新疆 石河子832003)
摘要:β多样性是表示植物群落物种组成变化的一个重要指示因子。本研究以古尔班通古特沙漠南部为研究区域,
在67个样地植被学调查基础上,首先对研究区植物群落β多样性经向和纬向的变化进行了分析,其次利用方差分
解的方法分析环境变量,土壤种子库变量及空间变量对β多样性分异的贡献率,探讨沙漠植被β多样性变化机制。
结果表明,β多样性在经向上由西向东物种相似性增加,物种替代速率减弱;在纬向上自北向南物种相似性降低,物
种替代速率略有升高。对β多样性分异解释方面,环境变量中的30年月平均降水量等变量在沙漠植被群落分布中
起着决定性的作用,土壤变量与空间变量的解释比率则相对较低,而土壤种子库变化较土壤异质性对β多样性影
响更大。主环境变量、土壤变量的效应是尺度依赖的。空间结构决定的随机过程也发挥着作用。研究还发现除环
境变量(含土壤变量),土壤种子库变量及空间变量外,仍有近50%的β多样性分异未能解释,表明其他环境因子和
随机过程对植被的影响也较大。
关键词:β多样性;方差分解;主轴邻矩法;古尔班通古特沙漠
中图分类号:Q948.15  文献标识码:A  文章编号:10045759(2010)03002909
  物种多样性是生态系统的基本特征,对生态系统的功能和稳定起决定性作用[13]。β多样性是指沿着环境梯
度变化物种替代的程度,也称物种周转率,用来反映群落间物种组成的变化[46],是决定地区尺度物种多样性的重
要因素[7],是理解生态系统功能、生物多样性保护和生态系统管理的一个重要指标[8,9]。
目前,在研究β多样性分异机理方面,Borcard和Legendre
[10]提出了在多尺度下解开β多样性的新方法———
方差分解方法。该方法就是将群落组成数据表的总方差无偏分解成由各个解释变量所决定的子方差,调整的决
定系数提供了针对多元回归和典范冗余分析的无偏估计,以此来评估各变量对β多样性空间格局的贡献
[911]。
该方法已得到普遍认可[1216]。
古尔班通古特沙漠(以下简称“古漠”)位于准噶尔盆地,东西长412km,南北最宽处225km,是中国第二大
沙漠。古漠的气候较中亚相同纬度沙漠干燥,但是,由于均匀的季节性降水,特别是冬、春季的雨雪,植物多样性
丰富,群落类型多样,沙漠呈固定和半固定状态[17],尤其是古漠南部植被对防止沙漠移动,阻挡沙漠向绿洲扩张
起重要作用[18]。那么该地区物种多样性组成,空间分异特征研究对了解物种分布状况,提出生态保护和恢复措
施至关重要。
为此,本研究以古漠南部植被为研究对象,采用分层随机取样的方法,对研究区植被进行全面调查,利用主轴
邻矩法(principalcoordinatesofneighbormatricesanalysis,PCNM)确定空间变量[19],结合土壤理化性质、海拔、
降水等环境变量,以及土壤种子库变量,通过方差分解方法,重点研究:1)古漠南部植物群落β多样性特征;2)各
变量对β多样性分异的贡献率,对β多样性进行解释。
1 材料与方法
1.1 研究区概况
研究区位于古漠N44.40°~45.56°,E84.91°~88.83°区域(图1),东西约280km,南北约70km,样点覆盖
第19卷 第3期
Vol.19,No.3
草 业 学 报
ACTAPRATACULTURAESINICA   
29-37
2010年6月
 收稿日期:20100203;改回日期:20100322
基金项目:国家十一五科技支撑计划项目 (2007BAC17B03)资助。
作者简介:赵怀宝(1969),男,河南永城人,在读博士。Email:ktszhb@sina.com
通讯作者。Email:liutong1968@yahoo.com.cn
了整个沙漠近1/2。主要沙丘类型为梁窝状和树枝状沙垄,近南北走向,沙垄高度5~50m。沙垄西坡缓而长,东
坡短且陡,垄顶部呈穹状,并有一个活化流动带。古漠属典型的内陆荒漠气候。年均温6~10℃,年蒸发量2000
mm以上,年均降水量普遍不超过150mm,沙漠腹地仅有70~100mm,主要集中在早春到初夏时节,不同年份和
区域差异很大。边缘的潜水位在5m以下,腹地则深逾23m。起沙风集中于4-6月[20]。冬季积雪和春季降水
改善了土壤墒情,良好的水热配置为植物的生长创造了条件。主要植被有白梭梭(犎犪犾狅狓狔犾狅狀狆犲狉狊犻犮狌犿),沙拐枣
(犆犪犾犾犻犵狅狀狌犿犿狅狀犵狅犾犻犮狌犿)等,短命和类短命植物普遍发育,丘间低地有结皮分布[21]。
图1 研究区及样地
犉犻犵.1 犛狋狌犱狔犪狉犲犪犪狀犱狊犪犿狆犾犻狀犵狆犾狅狋狊
图中阴影部分为古漠的轮廓,黑点及右上角数字为样地及代号,样地非连续编号,共有67处样地
TheshadypartistheprofileofGurbantunggutdesert,thesoliddotsindicatethesampled
plotsanditcodenearby,theywerenamednonsequenced
1.2 植被调查
根据生境特点,采取分层随机取样的方法,设置永久性样地67处(图1)。2009年5月,在短命植物种子成熟
时,采用典型样带法进行植物多样性调查。在每样地随机选择典型沙垄,以垂直于沙垄走向的方向设置100~
200m长,10m宽的样带,覆盖整个沙垄,包括阳坡底部,阳坡中部,垄顶,阴坡中部,阴坡底部5个部位。将样带
划分成10m×10m的小样方,采用连续样方法记录每个样方中全部植物的种类、株数、株高、冠幅等数据。同时
记录样方位置和坡位。物种鉴定到种,物种命名依据中国科学院植物研究所(IBAS)[22]。
1.3 β多样性分析
1.3.1 β多样性测度方法 采用样地组间物种组成相似性(Srenson指数)和替代性(Cody指数)来反映β多样
性空间分异。
Srenson指数:犐犃犮=2犮/(犪+犫) (1)
Cody指数:β犆=(犪+犫-2犮)/2 (2)
03 ACTAPRATACULTURAESINICA(2010) Vol.19,No.3
式中,犪和犫分别为2个样地群落各自的物种数,犮为2个样地群落的共有物种数。
1.3.2 β多样性分析方法 在研究区内,经向上以15′为间隔,自西向东将样地划分成17组,同组内的各样地物
种进行合并,比较相邻2组之间β多样性指数的变化。纬向上以6′为间隔,自北向南将样地划分成13组,同组内
的各样地物种进行合并,比较相邻2组之间β多样性指数的变化,以此来研究β多样性经向和纬向的分异特征。
1.4 解释变量选取
1.4.1 环境变量 以各样地海拔、降水、土壤理化性质作为环境变量。海拔记录于样地。各样地降水数据是以
德国气象局全球降水气候中心(globalprecipitationclimatologycentre,GPCC)[23]30年降水数据为基础,经过差
值模型计算出各样地降水数据。
GPCC采用SPHEREMSAP方法(球坐标距离和角加权方案)将气象站获得的常规观测资料内插到规则的
0.5deg(°)×0.5deg(°)经、纬度网格上,获得不同区域的降水资料[24],并已在气候学、水文学得到广泛应
用[24,25]。具体计算方法如下:以GPCC提供的1978-2007年30年的0.5deg×0.5deg分辨率原始降水数据为
基础,采用二维三次样条插值算法对0.5deg×0.5deg分辨率数据进行插值处理,得到0.1deg×0.1deg分辨
率的GPCC基准数据。以此基准数据为基础,对比观测样点附近10个团场气象站30年对应空间位置的实测降
水记录,求得差异数据。对差异数据在二维网格上进行二维散乱数据点的三次样条插值(interpolationof
2dimensionalscatteringdata),将得到的结果与 GPCC基准数据合成,得出经过实测数据校正的样点1978-
2007年的各月降水数据。包括30年间各月的月平均降水量(12列数据)及30年的年平均降水量(1列数据),连
同海拔(1列数据),共构成67×14矩阵。
土壤理化性质是对样地沙垄5个坡位0~10cm土样进行理化指标测定,测定方法见参考文献[26]。指标包
括有机质(SOM)、全氮(TN)、全磷(TP)、速效磷(AP)、速效钾(AK)、电导率(EC)、pH、Na+、Cl-、SO42-、Ca2+、
Mg2+、CO32-、HCO3-。
为了更好地独立分析土壤理化性质对植被组成的影响,将土壤理化性质构成的变量单独定义为土壤变量,分
析中各指标采用取样地各坡位的平均值;而海拔与降水变量被定义为主环境变量。
1.4.2 土壤种子库变量 种源变化对植被群落物种组成变化具有一定影响[27,28],因此,对样地土壤种子库物种
组成进行研究。2008年11月,在白梭梭等植被种子散落后取样。由于沙漠中土壤种子约有80%~90%分布在
上层2~5cm处[29,30],在5个坡位上分为0~3,3~6,6~10cm三层,分别在植冠下与裸地取样,样方大小为0.25
m×0.25m。在沙漠边缘沙地进行野外萌发试验。2008年12月,将土样均匀铺于准备好的萌发格中,自然越
冬,2009年春、夏季滴灌补水,抚育管理,不间断统计其萌发情况,并进行幼苗鉴定,得到土壤种子库物种多度数
据表。数据分析采用整个样地的萌发物种多度数据。
1.4.3 PCNM变量确定 PCNM,即主轴邻矩法,是一种较新的、有效地被用于分析空间位置关系对植被组成
影响的方法。PCNM变量是基于样地空间关联矩阵的特征函数,代表了样地之间空间关系。具体算法实现需要
输入各样地间距离参数(通过样地经、纬度在ArcGIS地理信息系统软件中可以获得),最终得到的PCNM 变量
是一个向量集合,由于所得到的各特征向量彼此正交,通过前选后,对应于大、中、小尺度的显著性PCNM特征向
量可被直接整合到回归模型中,为分析空间位置关系对植被组成的贡献率做准备。Dray等[19]在其文献中给出
了详细的PCNM步骤。
1.5 多变量方差分解
为探讨主环境变量、土壤变量、土壤种子库及空间位置关系对β多样性(应变量)分异的影响,应用典范冗余
分析(canonicalredundancyanalysis,RDA)对群落组成数据进行方差分解,进而对β多样性的机制作出解
释[15,16,31,32]。Legendre等[16]建议的β多样性指数SS(Y)正好是不同样地群落物种组成差异的直接测度,对应于
β多样性的概念。
犛犛(犢)=∑

犻=1


犼=1
(狔犻犼-珔狔犼)2 (3)
式中,犻和犼分别代表样地和物种,狀和狆分别是样地和物种的数量,狔犻犼为物种犼在样地犻的多度,珔狔犼 为物种犼在
13第19卷第3期 草业学报2010年
所有样地的平均多度。方差分解中应变量直接采用样地-物种相对多度矩阵。群落物种组成相对多度数据包括
100个物种,为了降低稀有种对分析的影响,剔除了相对频度小于15/67的物种,并对物种多度进行了helinger
转换[16]。
1.6 统计、计算方法
所有统计分析均通过R软件实现,其中PCNM 特征函数的获得和前选分别通过Spacemaker、Packfor包来
实现;方差分解中典范冗余分析、显著性检验均通过R软件的Vegan等数据包实现[3335]。
2 结果与分析
2.1 β多样性的空间分异
随着经度的增加,自西向东Srenson指数总体上呈上升趋势(图2A),Cody指数则呈下降趋势(图2B),说
明随着经度的增加,群落组成更加相似,而且群落物种替代速率也在减弱。特别是在西部85.00°~85.75°范围
内,指数变化十分明显,而此区域正处于农业开垦、采油等人类干扰频繁的低矮沙丘和地形较平坦区域向大沙丘
过渡的区域。85.2°附近样地4为沙漠中一块湿地边缘,生境较为独特,指数也呈急剧变化。在88.25°以东区域,
靠近北部(样地22,23,24)沙丘活化明显,物种替代明显加强,出现了大量刺沙蓬(犃犵狉犻狅狆犺狔犾犾狌犿狊狇狌犪狉狉狅狊狌犿)、
蛇麻黄(犈狆犺犲犱狉犪犱犻狊狋犪犮犺狔犪)群落;该区域南部(样地19,20)则出现了大量矮刺苏(犆犺犪犿犪犲狊狆犺犪犮狅狊犻犾犻犮犻犳狅犾犻狌狊)、
心叶驼绒藜(犆犲狉犪狋狅犻犱犲狊犲狑犲狉狊犿犪狀狀犻犪狀犪),准噶尔无叶豆(犈狉犲犿狅狊狆犪狉狋狅狀狊狅狀犵狅狉犻犮狌犿)等物种组成的群落。
在纬向方面,自北向南Srenson指数总体上呈下降趋势(图2C),说明随着纬度的降低,群落组成相似性降
低,而群落物种替代速率则略有增加。在南缘44.5°~44.8°区域范围内,Srenson指数增加明显,但Cody指数
变化则不明显(图2D),这表明从和绿洲交接的沙漠边缘地带向沙漠腹地延伸,生境相似性增加,群落物种组成相
似性也在增加。在45.2°~45.5°区域,指数出现急剧变化,物种替代显著下降,相似性急剧提高,可能是由于距离
绿洲较远,生境较为均一所致,但在最北部又发生了相反的变化。综合经向、纬向来看,研究区南缘和西部β多样
性指数变化较为剧烈,而中间有所波动,但变化不大。
图2 各样地组犛狉犲狀狊狅狀指数和犆狅犱狔指数沿纬向、经向的变化
犉犻犵.2 犆犺犪狀犵犲狅犳犛狉犲狀狊狅狀犻狀犱犲狓犪狀犱犆狅犱狔犻狀犱犲狓犪犾狅狀犵狋犺犲犾犪狋犻狋狌犱犲,犾狅狀犵犻狋狌犱犲犱犻狉犲犮狋犻狅狀犪犿狅狀犵狆犾狅狋犵狉狅狌狆狊
23 ACTAPRATACULTURAESINICA(2010) Vol.19,No.3
2.2 β多样性的变量解释
在选定的4个解释变量中,通过方差分解分析各变量对β多样性的独立贡献率。为了更清晰地反映各变量
的贡献率,将主环境变量和PCNM变量分别与土壤变量,土壤种子库变量进行分析。
2.2.1 应用主环境、PCNM、土壤变量的β多样性方差分解 首先对主环境及土壤变量进行前选,前选后的显著
性主环境变量包括海拔和4,5,6,10,11月30年的月平均降水(α=0.05);土壤变量包括:Na+,SOM,Ca2+,
Mg2+,AP(α=0.1)。
PCNM分析得到45个特征值。其中,正特征值25个(犚犪2=0.2571),前选后得到10个显著的PCNM特征
函数来模拟群落的物种组成(犚2=0.3694,犚犪2=0.2568),它们分别是PCNM1,PCNM9,PCNM3,PCNM10,
PCNM18,PCNM8,PCNM22,PCNM20,PCNM2,分别代表了大尺度(PCNM1,2,3)、中尺度(PCNM8,9,10)、小
尺度(PCNM18,20,22)。
图3 方差分解图
犉犻犵.3 犞犪狉犻犪狋犻狅狀狆犪狉狋犻狋犻狅狀犻狀犵狉犲狊狌犾狋狊
解释变量A,B,C,D分别是主环境变量、土壤变量、PCNM变量、土壤种子库变量 Astandsforenvironmentalvariables;
Bforsoilvariables;CforPCNMeigenfunctions;Dstandsforsoilseedbankvariables
方差分解的结果如图3中左图所示。主环境变量、土壤变量及PCNM 变量三者共同解释了空间分异的
49.0%。它们分别独立解释了分异的15.0%,1.0%,3.0%,而和其他变量配合后分别解释了分异的45.0%,
18.0%,25.0%,这说明主环境变量影响最大,而PCNM变量的影响大于土壤变量,并协同作用。其中,主环境变
量和PCNM变量共同解释了47.8%,土壤变量和PCNM 变量共同解释了34.2%。在控制环境变量时,土壤变
量和PCNM变量分别只解释了1.0%和3.0%。在控制PCNM 变量时,主环境变量和土壤变量分别只解释了
22.4%和9.0%,说明它们的影响是尺度依赖的。同时,也说明群落组成数据能够被PCNM特征函数所解释。
2.2.2 应用主环境、PCNM、土壤种子库变量的β多样性方差分解 对土壤种子库变量前选后得到影响群落物
种组成变化的物种包括:簇花芹(犛狅狉犪狀狋犺狌狊犿犲狔犲狉犻)、齿稃草(犛犮犺犻狊犿狌狊犪狉犪犫犻犮狌狊)、倒披针叶虫实(犆狅狉犻狊狆犲狉犿狌犿
犾犲犺犿犪狀狀犻犪狀狌犿)、茧荚黄芪(犃狊狋狉犪犵犪犾狌狊犾犲犺犿犪狀狀犻犪狀狌狊)、蓝刺头(犈犮犺犻狀狅狆狊狊狆犺犪犲狉狅犮犲狆犺犪犾狌狊)、琉苞菊(犎狔犪犾犲犪
狆狌犾犮犺犲犾犾犪)、尖喙?牛儿苗(犈狉狅犱犻狌犿狊狋犲狆犺犪狀犻犪狀狌犿)。方差分解的结果如图3中右图所示。3个变量共同解释了
空间分异的52.0%,分别独立解释了分异的16.5%,4.0%,3.0%,而和其他变量配合后分别解释了分异的
45.0%,21.0%和25.0%。
土壤种子库变量较土壤变量有更高的解释比例(图3),说明土壤种子库变量(即种源)较土壤的异质性对沙
漠植被群落组成影响更大。在控制土壤种子库变量时,主环境变量和PCNM 变量分别只解释了27.7%和
14.4%。大环境和土壤种子库变量共同解释了分异的48.9%,主环境变量和PCNM变量共同解释了47.8%,而
在控制主环境变量时,土壤种子库变量和PCNM变量分别只解释了4.0%和3.0%,说明主环境变量对群落β多
33第19卷第3期 草业学报2010年
样性变化起决定性的作用。
3 讨论
β多样性空间分布格局可以决定保护区应该设置在哪里,β多样性高的区域需要给予更多的关注和更大的保
护面积[35]。对研究区的β多样性分异分析表明,在经向上,自西向东物种相似性增加,替代速率下降。人类扰动
频繁的西部区域,经度85.00°~85.75°范围内替代速率变化剧烈;自西向东随着海拔的增加,降水增加明显,适生
物种增加较快,替代性则有下降。88.25°以东靠北区域,指数变化也较为剧烈,调查也发现风蚀、沙丘活化严重,
物种单一。在纬向上,自北向南物种相似性降低,替代速率略有增加,这说明绿洲环境对距离较近的沙漠区植被
组成有影响。根据β多样性空间分异特征,研究区西部和东北部区域在生态监测和管理中急需密切关注。
本研究对β多样性方差分解的结果表明,主环境变量具有较高的解释比例,这说明降水等主环境变量在沙漠
植被群落分布中起决定性的作用。气候因子常被作为群落组成的决定因素[36]。降水对植被分布的影响多见报
道[37,38]。降水也被作为干旱、半干旱区生态系统的限制性因子。海拔则提供了一个生态因子在不同空间尺度变
化的复合梯度[39],其通过影响降水量,相对湿度和温度对植物产生影响[40]。本研究中海拔及降水组成的主环境
变量对β多样性较大的贡献率也支持此观点。
土壤变量解释比例较低,说明在较大尺度上,沙漠土壤异质性对植被组成变化影响不大。一些研究也表明,
除了气候因子之外,土壤资源的变化在小尺度上是植被特征重要的决定因子[41]。但是,本研究显示土壤异质性
在解释大尺度沙漠植被群落组成变化方面比率很低,影响植被组成的土壤因子主要有 Na+,有机质,Ca2+,
Mg2+,速效磷。
种子定居的成功不能单纯归因于生境的适宜性,种源是植物多度和分布的决定因素[4244],物种的分布是由生
境质量和种子的扩张能力共同决定的。土壤种子库变量与土壤变量相比,解释比例较大,这从另一侧面反映出种
源在沙漠植被组成上较生境土壤异质性影响更大。
PCNM分析是获得样本间空间关系的分解向量,显著的PCNM 变量能直接地解释空间尺度。在研究的尺
度范围内,PCNM 分析通过对空间距离建立截尾矩阵,运用主轴分析法进行处理,得到PCNM1,PCNM2,…,
PCNMn,PCNM1代表整个研究尺度内的空间信息,属大尺度,PCNM2、…、PCNMn则分别描述越来越小的尺度
信息[11]。本研究中获得的PCNM变量在解释β多样性方面有一定贡献,说明空间结构决定的(PCNM 变量)随
机过程也发挥着作用,但比率较小。这也显示与降水等气候条件相比,空间关系对沙漠植被群落组成影响贡献有
限。
即使将土壤种子库因素也考虑其中,仍有48.0%的β多样性分异未能解释。这可能是由于其他未能考虑进
的环境因子(温度及降水的波动)和其他的随机过程(如人为干扰)等因素。
致谢:感谢加拿大 Montreál大学Pierrelegendre教授,中国科学院植物所米湘成研究员,刘国芳、黄继红博士提
供的指导。感谢新疆农业大学王兵教授、石河子大学阎平教授、翟伟硕士对物种鉴定方面提供的帮助。感谢李
磊、唐诚、张荣、陈辉煌等硕士生在多样性调查、土壤分析测定中的无私帮助。
参考文献:
[1] ChapinSF,SchulzeED,MooneyHA.Biodiversityandecosystemprocess[J].TrendsinEcologyandEvolution,1992,7:
107108.
[2] HectorA,BagchiR.Biodiversityandecosystemmultifunctionality[J].Nature,2007,448:188190.
[3] HooperDU,ChapinFS,EwelJJ,犲狋犪犾.Effectsofbiodiversityonecosystemfunctioning:Aconsensusofcurrentknowledge[J].
EcologicalMonographs,2005,75:335.
[4] WhittakerRH.Evolutionandmeasurementofspeciesdiversity[J].Taxon,1972,21:213251.
[5] MagurranA.Ecologicaldiversityanditsmeasurement[M].NewJersey:PrincetonUniversityPress,1988.
[6] 左小安,赵哈林,赵学勇,等.科尔沁沙地不同恢复年限退化植被的物种多样性[J].草业学报,2009,18(4):916.
43 ACTAPRATACULTURAESINICA(2010) Vol.19,No.3
[7] HarrisonS,RossSJ,LawtonJH.BetadiversityongeographicalgradientsinBritain[J].JournalofAnimalEcology,1992,
61:151158.
[8] LegendreP,BorcardD,PereNetoPR.Analyzingbetadiversity:Partitioningthespatialvariationofcommunitycomposition
data[J].EcologicalMonographs,2005,75:435450.
[9] LegendreP.Studyingbetadiversity:Ecologicalvariationpartitioningbymultipleregressionandcanonicalanalysis[J].Journal
ofPlantEcology,2007,31:976981.
[10] BorcardD,LegendreP.Alscalespatialanalysisofecologicaldatabymeansofprincipalcoordinatesofneighbourmatrices[J].
EcologyModeling,2002,153:5168.
[11] BorcardD,LegendreP,AvoisJacquetC,犲狋犪犾.Dissectingthespatialstructureofecologicaldataatmultiplescales[J].
Ecology,2004,85:18261832.
[12] GilbertB,LechowiczMJ.Neutrality,niches,anddispersalinatemperateforestunderstory[J].ProceedingNationalAca
demicScienceUSA,2004,101:76517656.
[13] JonesM M,TuomistoH,BorcardD,犲狋犪犾.Explainingvariationintropicalplantcommunitycomposition:Influenceofenvi
ronmentalandspatialdataquality[J].Oecologia,2008,155:593604.
[14] KarstJ,GilbertB,LechowiczMJ.Ferncommunityassembly:therolesofchanceandtheenvironmentatlocalandinterme
diatescales[J].Ecology,2005,86:24732486.
[15] LalibertéE,PaquetteA,LegendreP,犲狋犪犾.Assessingthescale———specificimportanceofnichesandotherspatialprocesses
onbetadiversity:Acasestudyfromatemperateforest[J].Oecologia,2009,159:377388.
[16] LegendreP,MiXC,RenHB,犲狋犪犾.PartitioningbetadiversityinasubtropicalbroadleavedforestofChina[J].Ecology,
2009,90(3):663674.
[17] 张立运,陈昌笃.论古尔班通古特沙漠植物多样性的一般特点[J].生态学报,2002,22(11):19231932.
[18] 潘晓玲,王学才,雷加强.关于中国西部干旱区生态环境演变与控制的一些思考[J].地球科学进展,2001,16(1):2427.
[19] DrayS,LegendreP,PeresNetoPR.Spatialmodeling:Acomprehensiveframeworkforprincipalcoordinateanalysisof
neighbourmatrices(PCNM)[J].EcologicalModeling,2006,196:483493.
[20] 王雪芹,王涛,蒋进,等.古尔班通古特沙漠南部沙面稳定性研究[J].中国科学,2004,34(8):763768.
[21] ZhangYM,ChenJ,WangL.ThespatialdistributionpatternsofbiologicalsoilcrustsintheGurbantunggutDesert,North
ernXinjiang,China[J].JournaloftheAridEnvironment,2007,68:599610.
[22] IBAS(InstituteofBotany,AcademiaSinica).IconographiaCormophytorumSonicorum[A].Tomus15,Supplementum12
[M].Beijing:SciencePress,1983.
[23] Globalprecipitationclimatologycentre.Globalprecipitationanalysisproducts[EB/OL].(20081230)[20090415].http://
gpcc.dwd.de.
[24] AdlerRF,HuffmanGJ,ChangAlfred,犲狋犪犾.Theversion2globalprecipitationclimatologyproject(GPCP)monthlypre
cipitationanalysis(1979Present)[J].JournalofHydrometeorology,2003,4:11471167.
[25] NezlinNP,KostianoyAG,LiBL.Interannualvariabilityandinteractionofremotesensedvegetationindexandatmospher
icprecipitationintheAralSearegion[J].JournalofAridEnvironments,2005,62:677700.
[26] 中国土壤学会农业化学专业委员会.土壤农业化学常规分析方法[M].北京:科学出版社,1983.
[27] NathanR,MulerLandauHC.Spatialpatternsofseeddispersal,theirdeterminantsandconsequencesforrecruitment[J].
TrendsinEcologyandEvolution,2000,15:278285.
[28] 尚占环,徐鹏彬,任国华,等.土壤种子库研究综述———植被系统中的作用及功能[J].草业学报,2009,18(2):175183.
[29] GuoQF,RundelPW,GoodalDW.Structureofdesertseedbanks:ComparisonsacrossfourNorthAmericandesertsites[J].
JournalofAridEnvironments,1999,42:114.
[30] 尚占环,任国华,龙瑞军.土壤种子库研究综述———规模、格局及影响因素[J].草业学报,2009,18(1):144154.
[31] RaoCD.Theuseandinterpretationofprincipalcomponentsanalysisinappliedresearch[J].Sankhyaá,1964,26:329358.
[32] LegendreP,LegendreL.NumericalEcology[M].SecondEnglishEdition.Amsterdam:ElsevierScienceBV,1998.
[33] RDevelopmentCoreTeam.R:Alanguageandenvironmentforstatisticalcomputing[EB/OL].(20090824)[20091015].
53第19卷第3期 草业学报2010年
http://www.Rproject.org.
[34] OksanenJ,KindtR,LegendreP,犲狋犪犾.Vegan:Communityecologypackage[J].RPackageVersion,2007,1:925.
[35] WiersmaYF,UrbanDL.BetadiversityandnaturereservesystemdesignintheYukon,Canada[J].ConservationBiology,
2005,19(4):12621272.
[36] AguadoSantacruzA,GarcíaMoyaE.EnvironmentalfactorsandcommunitydynamicsatthesouthernpartoftheNorth
AmericanGraminetumI.Onthecontributionofclimaticfactorstotemporalvariationinspeciescomposition[J].PlantEcolo
gy,1998,135:1329.
[37] WhittakerRH,NieringW A.VegetationoftheSantaCatalinaMountains,ArizonaV.Biomass,productionanddiversity
alongtheelevationgradient[J].Ecology,1975,56:771790.
[38] NobelPS.RemarkableAgavesandCacti[M].NewYork:OxfordUniversityPress,1994.
[39] BarreraMD,FrangiJL,RichterLL,犲狋犪犾.Structuralandfunctionalchangesin犖狅狋犺狅犳犪犵狌狊狆狌犿犻犾犻狅forestsalonganalti
tudinalgradientinTierradelFuego,Argentina[J].JournalVegetationScience,2000,11:179188.
[40] HustonMA.BiologicalDiversity:TheCoexistenceofSpeciesonChangingLandscapes[M].Cambridge:CambridgeUniver
sityPress,1994.
[41] MikiT,KondohM.Feedbacksbetweennutrientcyclingandvegetationpredictplantspeciescoexistenceandinvasion[J].E
cologyLetter,2002,5:624633.
[42] OztasT,KocA,ComakliB.Changesinvegetationandsoilpropertiesalongaslopeonovergrazedanderodedrangelands[J].
JournalofAridEnvironments,2003,55:93100.
[43] EhrlenJ,MunzbergovaZ,DiekmannM,犲狋犪犾.Longtermassessmentofseedlimitationinplants:Resultsfroman11year
experiment[J].JournalofEcology,2006,94:12241232.
[44] ZobelM,OtsusM,LiraJ,犲狋犪犾.Issmalscalespeciesrichnesslimitedbyseedavailabilityormicrositeavailability?[J].
Ecology,2000,81(12):32743282.
[45] 李有志,张灿明,林鹏.土壤种子库评述[J].草业科学,2009,26(3):8390.
63 ACTAPRATACULTURAESINICA(2010) Vol.19,No.3
β犱犻狏犲狉狊犻狋狔犮犺犪狉犪犮狋犲狉犻狊狋犻犮狅犳狏犲犵犲狋犪狋犻狅狀犮狅犿犿狌狀犻狋狔狅狀狊狅狌狋犺狆犪狉狋狅犳
犌狌狉犫犪狀狋狌狀犵犵狌狋犇犲狊犲狉狋犪狀犱犻狋狊犻狀狋犲狉狆狉犲狋犪狋犻狅狀
ZHAOHuaibao1,2,LIUTong3,LEIJiaqiang2,GUIDongwei2,ZHAOXinjun4
(1.GraduateUniversityofCAS,Beijing100048,China;2.XinjiangInstituteofEcologyandGeography
CAS,Uruqmi830000,China;3.ColegeofLifeScience,ShiheziUniversity,Shihezi
832003,China;4.ColegeofNormal,ShiheziUniversity,Shihezi832003,China)
犃犫狊狋狉犪犮狋:βdiversityisanimportantindexcharacterizingcompositionvariationofvegetationcommunity.The
southernpartofGurbantunggutDesertistakenasthestudyareainthispaper,onthebasisofvegetationsur
veyin67plots,βdiversityvariationofvegetationcommunityinthestudyareaisanalyzedfirstinlongitudinal
andlatitudinaldirections,andthenthemethodofvariancepartitioningisusedforanalyzingcontributionratio
ofenvironmentvariable,soilseedbankvariableandspatialvariabletoβdiversitydifferentiation,discussingthe
mechanismofβdiversitydifferentiationofdesertvegetation.Itisshownwiththeresultsthatβdiversityisan
increaseofspeciessimilarityfromwesttoeastinlongitudinaldirection,andthespeciessuccessionratede
clines;adecreaseofspeciessimilarityfromnorthtosouthinlatitudinaldirection,andthespeciessuccession
rateincreasesslightly.Regardinginterpretationofβdiversitydifferentiation,theenvironmentvariablessuchas
30year’saveragemonthlyprecipitationplayadecisiveroleindistributionofdesertvegetationcommunity,
whiletheinterpretationratioofsoilvariableandspatialvariablearerelativelow.Comparingwithsoilheteroge
neity,thevariationofsoilseedbankhasabiggerimpactonβdiversitydifferentiation.Themainenvironmental
variables,soilvariableseffectwerescaledependent.Neutralprocessdeterminedbythespatialstructurealso
playsarole.Itisfoundoutthatalthoughenvironmentvariables(includingsoilvariable),soilseedbankvaria
bleandspatialvariablearetakenintoconsideration,nearly50%ofβdiversitydifferentiationcannotbeinterpre
ted,whichindicatesthatsomeotherenvironmentvariablesandrandomprocesshaveimpactonthevariationof
vegetationcommunity.
犓犲狔狑狅狉犱狊:βdiversity;variancepartitioning;principalcoordinatesofneighbormatricesmethod(PCNM);
GurbantunggutDesert
73第19卷第3期 草业学报2010年