全 文 :犇犗犐:10.11686/犮狔狓犫2015199 犺狋狋狆://犮狔狓犫.犾狕狌.犲犱狌.犮狀
罗黎鸣,武建双,余成群,潘影,苗彦军,武俊喜,明升平,郭应杰.拉萨河谷山地灌丛草地植物多样性监测方法的比较研究.草业学报,2016,25(3):
2231.
LUOLiMing,WUJianShuang,YUChengQun,PANYing,MIAOYanJun,WUJunXi,MINGShengPing,GUOYingJie.Acomparative
studyofmethodsformonitoringplantdiversityinmontaneshrubgrasslandintheLhasaRiverBasin.ActaPrataculturaeSinica,2016,25(3):2231.
拉萨河谷山地灌丛草地植物多样性
监测方法的比较研究
罗黎鸣1,2,武建双1,余成群1,潘影1,苗彦军3,武俊喜1,明升平4,郭应杰4
(1.中国科学院地理科学与资源研究所,生态系统网络观测与模拟重点实验室,北京100101;2.西藏农牧科学院农业研究所,
西藏 拉萨850002;3.西藏大学农牧学院,西藏 林芝860000;4.西藏大学理学院,西藏 拉萨850001)
摘要:生物多样性是西藏生态安全监测与评估的重要内容之一。在气候变暖和放牧干扰背景下沿海拔梯度探讨西
藏拉萨河谷山地灌丛草地植物群落物种多样性的变化特征,对于该地区草地的保护与可持续利用具有重要意义。
以拉萨市林周县白朗村沟域的山地灌丛草地为研究对象,采用常规样方法、重复样线法和 Modified-Whittaker样
地法3种方法进行系统调查,研究草地群落物种组成及其沿海拔梯度的变化特征。结果表明,1)3种不同取样方
法α多样性的物种丰富度均随海拔的升高呈先增加后降低的趋势,其中50m样线法效果最好。2)表征β多样性
的Sorensen相似性指数和BrayCurtis指数均随样地间海拔高差的增大而减小。3)种-面积关系的截距犮、斜率狕
随海拔的升高呈先增加后降低的趋势,与α、β多样性随海拔梯度的变化特征吻合一致。4)依据种-面积关系研究
结果,山地灌丛草地设置100m2 的样地可覆盖80%左右的植物物种。山地灌丛草地的α和β多样性指数沿海拔
梯度均呈先增后减的变化特征,未来研究应注重植物多样性分布格局对山体气候梯度以及放牧干扰变化的响应和
适应机制。
关键词:西藏;海拔梯度;物种多样性;监测方法;山地灌丛草地
犃犮狅犿狆犪狉犪狋犻狏犲狊狋狌犱狔狅犳犿犲狋犺狅犱狊犳狅狉犿狅狀犻狋狅狉犻狀犵狆犾犪狀狋犱犻狏犲狉狊犻狋狔犻狀犿狅狀狋犪狀犲狊犺狉狌犫
犵狉犪狊狊犾犪狀犱犻狀狋犺犲犔犺犪狊犪犚犻狏犲狉犅犪狊犻狀
LUOLiMing1,2,WUJianShuang1,YUChengQun1,PANYing1,MIAOYanJun3,WUJunXi1,MING
ShengPing4,GUOYingJie4
1.犓犲狔犔犪犫狅狉犪狋狅狉狔狅犳犈犮狅狊狔狊狋犲犿犖犲狋狑狅狉犽犗犫狊犲狉狏犪狋犻狅狀犪狀犱犕狅犱犲犾犻狀犵,犐狀狊狋犻狋狌狋犲狅犳犌犲狅犵狉犪狆犺犻犮犛犮犻犲狀犮犲犪狀犱犖犪狋狌狉犪犾犚犲狊狅狌狉犮犲狊犚犲
狊犲犪狉犮犺,犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊,犅犲犻犼犻狀犵100101,犆犺犻狀犪;2.犃犵狉犻犮狌犾狋狌狉犪犾犚犲狊犲犪狉犮犺犐狀狊狋犻狋狌狋犲,犜犻犫犲狋犃犮犪犱犲犿狔狅犳犃犵狉犻犮狌犾狋狌狉犪犾
犪狀犱犃狀犻犿犪犾犛犮犻犲狀犮犲狊,犔犺犪狊犪850002,犆犺犻狀犪;3.犃犵狉犻犮狌犾狋狌狉犪犾犪狀犱犃狀犻犿犪犾犎狌狊犫犪狀犱狉狔犆狅犾犾犲犵犲狅犳犜犻犫犲狋犝狀犻狏犲狉狊犻狋狔,犔犻狀狕犺犻860000,
犆犺犻狀犪;4.犆狅犾犾犲犵犲狅犳犛犮犻犲狀犮犲,犜犻犫犲狋犝狀犻狏犲狉狊犻狋狔,犔犺犪狊犪850001,犆犺犻狀犪
犃犫狊狋狉犪犮狋:Biodiversityisanimportantconsiderationinthemonitoringandassessmentoftheecologicalstability
ofrangelandsinTibet.ClimatewarmingandgrazingdisturbanceboththreatenthestabilityoftheTibetan
rangelands;hence,researchintospeciessuccession,transitionpatternsandspeciesdiversityisimportantfor
theprotectionandutilizationofgrasslandsintheregion.Afieldsurveywascarriedouttomeasurethecommu
22-31
2016年3月
草 业 学 报
ACTAPRATACULTURAESINICA
第25卷 第3期
Vol.25,No.3
收稿日期:20150416;改回日期:20150930
基金项目:国家自然青年基金项目(41401070),国家科技支撑计划项目(2011BAC09B03,2011BAD17B05)和“西部之光”人才培养计划项目
(20100309)资助。
作者简介:罗黎鸣(1989),男,湖南邵阳人,硕士。Email:luolm1989tibet@126.com
通信作者Correspondingauthor.Email:wujx@igsnrr.ac.cn
nityspeciescompositionandthespecieschangealonganaltitudinalgradientintheBailangvilagevaley,Linz
houCountyintheLhasaRiverBasin.Threeconventionalsamplingmethods,randomsquarequadrats,repeat
samplinglinesandtheModified-Whittakerplots,werecomparedinthisstudy.Keyresultswere:(1)Al
threesamplingmethodsindicatedthatαdiversityinitialyincreasedandthendecreasedwithincreasingaltitude,
thoughthelinetransectmethodusingthree50mtransectsdetectedmorespeciesthantheothertwomethods.
(2)TheSorensensimilarityindexandtheBrayCurtisindexofβdiversitydecreasedwithincreasingaltitude.
(3)Theintercept犮andslope狕oftherelationshipbetweenspeciesandsamplingareashowedunimodalpatterns
withincreasingaltitudeasalsodidspeciesrichness.(4)Inordertocover80%oflocalplantspecies,themini
mumsamplingareashouldbeupto100m2.Futureresearchshouldfocusonthefactorsdeterminingbiodiversi
tydistributionandcommunityassemblyinresponsetoclimatechangeandgrazingdisturbance.
犓犲狔狑狅狉犱狊:Tibet;altitudinalgradient;speciesdiversity;monitoringmethods;montaneshrubgrassland
物种多样性通常从不同层次、尺度分α多样性、β多样性、γ多样性进行度量
[1]。其中,β多样性联系群落局
部α多样性和区域γ多样性,反映区域内的物种组成随时间、空间的变化,对于理解生态系统功能、加强生物多样
性保护和优化生态系统管理具有重要意义[14]。物种多样性野外调查方法很多,如样方法和样线法,各有优缺点。
样方法[5]有覆盖式的正方形镶嵌式取样方法,如巢式取样法;多尺度长方形镶嵌式取样方法,如 Whittaker样地
法和 Modified-Whittaker样地法。巢式样方因其简便实用在实际调查中经常应用,但其并不足以反映整个地
段上的群落特征[6]。而 Modified-Whittaker样地法特别适合空间异质性较大的生态系统,能更好地分析物种
丰富度[5],广泛应用于全球不同类型生态系统的物种调查中[78],也适合研究天然草地物种多样性的空间异质
性[9]。样线法是美国监测草地物种丰富度的主要方法[1011]。样线的长度,重复数多少,间隔距离多少,都会影响
样线法监测物种丰富度的精度。在物种多样性研究中取样方法和面积大小在很大程度上会影响研究结果[12],种
-面积关系是常用来确定特定群落取样最小的取样面积[13]。国内已有不少研究[1416]探讨了不同区域、不同草地
类型的最优取样法和最小取样面积。研究方法的一致性和可比性是保证数据质量的关键环节[15],同一地区采用
不同调查方法取得的物种多样性各指标存在差异,故此有必要比较多种调查监测方法以确定最优取样法和最小
取样面积。
青藏高原因其独特的高寒气候和复杂多样的自然条件,造就了丰富多样、独具特色的草地生态系统,是我国
宝贵的高原生物基因库和可贵的景观资源[1718]。高原地区的生物多样性是生态安全监测与评估的重要内容之
一[18]。孟凡栋等[19]已对西藏高原3种不同类型的高寒草地植物多样性和地上生物量监测进行了比较研究。董
世魁等[20]对青藏高原不同状况(健康和退化)的高寒草地植物多样性的最小样地面积进行研究。而针对青藏高
原河谷地区山地灌丛草地植物多样性监测方法的比较研究甚少。
本文采用常规样方法、重复样线法和 Modified-Whittaker样地法3种方法对拉萨山地灌丛草地生物多样
性进行系统调查。研究植物群落α、β多样性随海拔升高的变化特征,验证种-面积关系参数同α、β多样性的关
系,探究山地灌丛草地群落植物多样性的最优取样法和取样面积,旨在揭示西藏“一江两河”地区山地灌丛草地群
落物种组成及其沿海拔梯度的变化特征,为该地区天然草地的保护与可持续利用提供科学依据和理论基础。
1 材料与方法
1.1 研究区域概况
“一江两河”地区属高原温带季风半干旱气候区,年平均气温在4.7~8.3℃,年降水量251.7~580.0mm之
间[21]。土地面积6.65×106hm2,占西藏总面积的5.52%;天然草地4.75×106hm2,占该地区土地面积的
71.44%。白朗村沟域位于 “一江两河”的拉萨河流域支流澎波藏曲,地处于念青唐古拉山脉中段南翼,行政隶属
于西藏自治区拉萨市林周县卡孜乡白朗村。林周县域海拔3900~3950m,年日照时数>3000h,年均温7.5~
7.6℃,年均降水量440mm左右且主要集中在6-9月。
32第25卷第3期 草业学报2016年
试验地位于白朗村沟域的阴坡上,土壤类型由坡脚至山顶分别发育是亚高山草甸草原土、亚高山灌丛草甸
土、高山草甸土;对应的植被类型由下至上分别是亚高山草甸草原、亚高山灌丛草甸、高山草甸。其中,灌丛层片
以雪层杜鹃(犚犺狅犱狅犱犲狀犱狉狅狀狀犻狏犪犾犲)、微毛樱草杜鹃(犚犺狅犱狅犱犲狀犱狉狅狀狆狉犻犿狌犾犻犳犾狅狉狌犿)、金露梅(犘狅狋犲狀狋犻犾犾犪犳狉狌狋犻犮狅
狊犪)为主。草本层建群种主要为高山嵩草(犓狅犫狉犲狊犻犪狆狔犵犿犪犲犪)。在整个山体调查中,共记录植物种121个,分属
32科,70属。其中灌木植物4种,草本植物117种。样地基本信息见表1。
表1 样地基本信息
犜犪犫犾犲1 犘犾狅狋犮犺犪狉犪犮狋犲狉犻狊狋犻犮狊
海拔
Altitude(m)
经度
Longitude
纬度
Latitude
坡度
Slopegradient(°)
坡向
Exposure(°)
科
Family
属
Genus
种
Species
4000 E91°07′07″ N29°52′30″ 8 245 16 31 39
4300 E91°07′20″ N29°50′59″ 17 300 23 45 55
4600 E91°07′07″ N29°50′29″ 18 295 22 39 58
4900 E91°06′59″ N29°50′12″ 55 325 21 36 62
5100 E91°07′18″ N29°49′51″ 18 250 14 21 39
1.2 试验方法
图1 犕狅犱犻犳犻犲犱-犠犺犻狋狋犪犽犲狉取样法示意图
犉犻犵.1 犜犺犲狊犽犲狋犮犺犿犪狆狅犳犕狅犱犻犳犻犲犱-犠犺犻狋狋犪犽犲狉狊犪犿狆犾犲
于2014年8-9月,对林周县白朗村沟域的山地
灌丛草地沿海拔梯度进行常规样方法、重复样线法和
Modified-Whittaker样地法系统调查。从山脚4000
m至山顶5100m沿坡面按植被类型差异,每隔300m
设置一个样地,共设5个样地,样地大小约50m×50
m;在每样地随机选取6个面积为0.5m×0.5m的样
方,做常规草地样方调查,记录每个样方各物种的高度
和盖度;在每个样地间隔20m 处设置3条50m 样
线,样线方向垂直于样地海拔等高线,每隔5cm针刺记录所出现的物种和次数,用以测定各物种频度和物种丰富
度;在每样地采用 Modified-Whittaker样地法[22]进行种-面积关系调查及最小面积的确定,在20m×50m的
大样方里分别取如下小样方:10个0.5m×1m、2个2m×5m、1个5m×20m、1个20m×50m(图1)。
1.3 数据处理
1.3.1 各物种重要值[23] 物种重要值(犐)=(犆狉+犎狉+犉狉)/3,其中物种相对盖度(relativecoverage,犆狉)和
相对高度(relationheight,犎狉)为常规样方法所得,物种相对频度(relativefrequency,犉狉)为重复样线法统计所得。
1.3.2 物种多样性 本研究中物种多样性分α多样性和β多样性。α多样性的测定采用物种丰富度指数犛。
该指数是表示物种多样性最客观、最真实的指标[24]。β多样性的测定
[25]采用Sorensen指数和BrayCurtis指数。
Sorensen相似性指数表示不同样地间共有物种占每个样地群落物种总数的比例,能体现群落物种组成出现的重
叠效应。BrayCurtis指数是数量化Sorensen指数,更能体现群落物种组成的变化[26]。Sorensen指数,犆狊=2犼/
(犪+犫),式中,犼为两个群落或样地共有物种数;犪和犫分别为样地A和样地B的物种数;BrayCurtis指数,犆犐=
2犼犐/(犪犐+犫犐),式中,犼犐为样地A(犼犐犪)和样地B(犼犐犫)共有种物种重要值较小者之和,犼犐=∑min(犼犐犪+犼犐犫),式
中,犪犐、犫犐分别为样地A和B物种重要值。本研究中α多样性和β多样性Sorensen指数的测算分别采用常规样
方法、重复样线法和 Modified-Whittaker样地法中1000m2 样地中的物种。BrayCurtis指数仅采用常规样方
法获取的相关数据。
1.3.3 种-面积关系 本研究中,种-面积关系的公式拟合采用目前公认最简单且最有效的幂函数的对数形
式[2728],即log犛=log犮+狕log犃;公式描述的种-面积关系在双对数坐标轴上是一条直线。一般认为,其犮值可以
42 ACTAPRATACULTURAESINICA(2016) Vol.25,No.3
用于表征单位面积的物种多样性;狕值则可表征随采样面积增大物种数的积累速率,是物种多样性面积效应的直
接体现。因此,种-面积曲线中的犮值与狕值可以近似理解为α多样性和β多样性
[29]。
1.3.4 最小面积的确定 本研究中,选用非饱和曲线犛=犫+犪ln犃作为种-面积曲线的拟合模型[30]。要得到
群落总物种数一定比例ρ(0<ρ<1)的物种所需的最小面积(或临界抽样面积)
[14]为:犃min=exp[(ρ犛狋-犫)/犪],式
中,犃为样地面积,犛为面积犃 中出现的物种数,犪、犫都是模型中待估参数,犛狋是调查群落物种总数。
1.4 数据分析
将各海拔梯度上所采的样方数据计算均值,通过SPSS19.0软件,利用单因子方差分析(OnewayANOVA)
的最小显著差数法(LSD)进行分析,作图采用Origin8.0软件。
2 结果与分析
2.1 随海拔升高山地灌丛草地植物群落α多样性比较
在3种不同取样方法下,山地灌丛草地植物群落α多样性(物种丰富度)均随海拔的升高呈“单峰曲线”变化
特征(图2)。常规样方法(图2a)在海拔4300m处物种数最多(34),在5100m处最少(25);重复样线法(图2b)
图2 3种取样方法下α多样性随海拔升高的变化
犉犻犵.2 犞犪狉犻犪狋犻狅狀狊狅犳α犱犻狏犲狉狊犻狋狔犪犾狅狀犵犪狀犪犾狋犻狋狌犱犲犵狉犪犱犻犲狀狋狑犻狋犺狋犺狉犲犲狊犪犿狆犾犻狀犵犿犲狋犺狅犱狊
在海拔4600m 处物种数最多(54),在4000m处最少
(34);Modified-Whittaker样地法(图2c)在海拔
4900m处物种数最多(62),在4000和5100m处最少
(39)。从物种丰富度最大值来看,Modified-Whit
taker样地法>重复样线法>常规样方法。从拟合效
果(犚2)来看,重复样线法>常规样方法>Modified-
Whittaker样地法。综合考虑,本研究推荐50m样线
重复3次可较好地反映这一地区山地灌丛草地物种的
α多样性。
2.2 随海拔升高山地灌丛草地植物群落β多样性比
较
从3种不同取样方法下Sorensen指数测度结果
矩阵(表2)来看,β多样性整体上呈与α多样性相同的
变化趋势。即Sorensen指数均随样地间海拔高差的
增大而减小,群落物种组成差异变大。常规样方法
表2 各海拔样地β多样性犛狅狉犲狀狊犲狀指数测度结果矩阵
犜犪犫犾犲2 犕犪狋狉犻狓狅犳β犱犻狏犲狉狊犻狋狔狅犳犛狅狉犲狀狊犲狀犻狀犱犲狓
犪犿狅狀犵狆犾狅狋狊狅犳犪犾狋犻狋狌犱犲犵狉犪犱犻犲狀狋狊
方法
Method
海拔
Altitude(m)
海拔Altitude(m)
4000 4300 4600 4900
常规样方法
Randomsquare
plots
4300 0.5806
4600 0.2807 0.4127
4900 0.1818 0.2951 0.7143
5100 0.0755 0.2373 0.5926 0.6923
重复样线法
Repeatsampling
lines
4300 0.5185
4600 0.2619 0.4848
4900 0.1282 0.2796 0.6667
5100 0.0909 0.2469 0.5714 0.7436
样地法
Modified- Whit
takersampling
4300 0.5745
4600 0.3505 0.4956
4900 0.1584 0.3419 0.7167
5100 0.1795 0.3191 0.5567 0.6535
52第25卷第3期 草业学报2016年
Sorensen指数在海拔 4600 和 4900 m 之间最大
(0.7143),在海拔4000和5100m之间最小(0.0755);
重复样线法Sorensen指数在海拔4900和5100m之
间最大(0.7436),在海拔4000和5100m 之间最小
(0.0909);Modified-Whittaker样地法Sorensen指
数在海拔4600和4900m之间最大(0.7167),在海拔
4000和4900m之间最小(0.1584)。
从整体上看,β多样性BrayCurtis指数测度结果
矩阵(表3)整体趋势和各Sorensen指数(表2)随海拔
的变化趋势基本一致。这说明二者均可用于群落物种
组成差异性测度且无明显差异。本研究中BrayCur
tis指数随样地间海拔高差的增大而减小,表明随海拔
梯度差的增大,物种组成差异越大。在海拔4600和
4900m之间最大(2.7729),在海拔4000和5100m之
间最小(1.3630)。从BrayCurtis指数在各海拔梯度
之间变化速率来看,中海拔较低海拔、高海拔变化慢,
说明中海拔生境差异性较低,物种组成差异较小。
2.3 随海拔升高山地灌丛草地植物群落种-面积关
系的变化
在同一海拔,随着取样面积的增大物种数呈显著
增加。在取样面积相同的情况下,物种数随海拔的升
高呈先增加后减少的趋势。在海拔5100m处,取样
表3 各海拔样地β多样性犅狉犪狔犆狌狉狋犻狊指数测度结果矩阵
犜犪犫犾犲3 犕犪狋狉犻狓狅犳β犱犻狏犲狉狊犻狋狔狅犳犅狉犪狔犆狌狉狋犻狊犻狀犱犲狓
犪犿狅狀犵狆犾狅狋狊狅犳犪犾狋犻狋狌犱犲犵狉犪犱犻犲狀狋狊
海拔
Altitude(m)
海拔Altitude(m)
4000 4300 4600 4900
4300 2.7450
4600 1.9330 2.1434
4900 1.6506 2.1556 2.7729
5100 1.3630 1.7721 2.3843 2.4455
表4 种-面积关系随海拔升高的变化
犜犪犫犾犲4 犞犪狉犻犪狋犻狅狀狊狅犳狊狆犲犮犻犲狊犪狉犲犪狉犲犾犪狋犻狅狀狊犺犻狆
犪犾狅狀犵犪狀犪犾狋犻狋狌犱犲犵狉犪犱犻犲狀狋
海拔
Altitude(m)
面积 Area(m2)
1 10 100 1000
4000 19.4±4.9a 27.5±3.5 33 39
4300 18.9±3.9ab 36.5±0.7 42 55
4600 21.3±3.4a 39.0±5.7 52 58
4900 25.0±2.3a 34.5±3.5 54 62
5100 15.9±2.8b 19.5±6.4 28 39
注:同列不同小写字母表示差异显著(犘<0.05)。
Note:Differentsmallettersinthesamecolumnindicatesignificant
differencesat犘<0.05.
面积为1m2 时物种数最少(15.9±2.8);在海拔4900m处,取样面积为1000m2 时物种数最多(62)(表4)。
种-面积关系截距犮随海拔的升高呈先升高后降低的趋势,海拔5100m处犮值最小(15.601),海拔4900m
处犮值最大(25.401)。这表明截距犮表征的α多样性亦呈“单峰曲线”分布特征,与物种丰富度的分布特征相吻
合。随海拔的升高种-面积关系斜率狕亦呈先增加后降低的趋势,表明β多样性随海拔的升高呈先增加后降低
的趋势。海拔4000m处狕值最小(0.103),海拔4600m处狕值最大(0.157),表明斜率狕表征的β多样性在低海
拔处最小,中海拔处最大。这可能是因为低海拔处生境差异性较高,物种组成差异较大;中海拔处生境较为均一,
物种组成差异较小(图3)。
图3 不同海拔种-面积关系比较
犉犻犵.3 犆狅犿狆犪狉犻狊狅狀狅犳狊狆犲犮犻犲狊-犪狉犲犪狉犲犾犪狋犻狅狀狊犺犻狆狊犫犲狋狑犲犲狀犱犻犳犳犲狉犲狀狋犪犾狋犻狋狌犱犲犵狉犪犱犻犲狀狋狊
62 ACTAPRATACULTURAESINICA(2016) Vol.25,No.3
2.4 随海拔升高山地灌丛草地植物群落最小取样面积的确定
拉萨河谷山地灌丛草地种-面积关系拟合结果如表5所示。模型均通过犉检验(犘<0.0001),但各海拔样
地存在明显差异,拟合曲线犚2 整体随海拔的升高呈先增加后降低的趋势,说明种-面积关系受海拔影响明显。
模型参数斜率(犪)、截距(犫)随海拔的升高呈先增加后降低的趋势。同一海拔最小样地面积随精度ρ的增加而增
大。同一精度ρ最小样地面积随海拔的升高呈先增加后降低再升高的趋势。为满足研究精度ρ=60%的最低研
究,最小面积至少约为10m2;为满足研究精度ρ=80%的精确研究,最小面积则至少约为100m
2。
表5 拉萨河谷山地灌丛草地不同海拔样地群落最小面积
犜犪犫犾犲5 犕犻狀犻犿狌犿犮狅犿犿狌狀犻狋狔犪狉犲犪狅犳犱犻犳犳犲狉犲狀狋犪犾狋犻狋狌犱犻狀犪犾狆犾狅狋狊狅犳犿狅狀狋犪狀犲狊犺狉狌犫犵狉犪狊狊犾犪狀犱犻狀犪犔犺犪狊犪犚犻狏犲狉犅犪狊犻狀
样地海拔
Plotaltitude
(m)
模型参数
Modelparameters
斜率Slope 截距Intercept
模型检验
Modeltest
犘 犚2
不同精度下的最小样地面积
Minimumplotareasunderdifferentaccuracy(m2)
ρ=60% ρ=70% ρ=80% ρ=90%
4000 2.927 19.559 <0.0001 0.668 4 14 53 202
4300 5.042 19.410 <0.0001 0.852 15 44 131 391
4600 6.344 21.696 <0.0001 0.934 8 20 49 123
4900 5.400 24.783 <0.0001 0.945 10 31 99 312
5100 3.002 15.472 <0.0001 0.778 14 51 189 691
山体 Mountain 10.2 32.0 104.2 343.8
3 讨论
3.1 α多样性随海拔升高的变化特征
α多样性指某一确定采样面积内的物种多样性,反映群落物种组成和结构的多样化程度,是物种多样性研究
的基础[31]。玉龙雪山的研究[24]表明木本、草本、全部物种的α多样性(物种丰富度)随海拔的升高均呈“单峰曲
线”。青藏高原草甸区的研究[3233]表明,α多样性随海拔的升高呈“单峰曲线”,即物种丰富度、Shannon-Wiener
多样性指数、Pielous均匀度指数随海拔的升高均呈先增加后降低的趋势。拉萨河谷地区山地灌丛草地的已有研
究[34]也表明,物种α多样性随海拔的升高呈“单峰曲线”。本研究结果验证了拉萨河谷地区山地灌丛α多样性随
海拔的升高呈“单峰曲线”。该地区山地灌丛草地群落α多样性的这种变化特征可能是低海拔受水分制约,高海
拔受温度制约,中海拔常有云雾覆盖、水热条件相对均衡等气候差异造成。通过常规样方法、重复样线法和
Modified-Whittaker样地法3种取样方法的比较,我们认为常规样方法受草地空间异质性的影响太大,人为选
择样地会对结果造成较大偏差,研究中应尽量避免样地选择的主观性;Modified-Whittaker样地法样地面积大,
覆盖的物种更多,主要是涵盖了少量稀有种,但稳定性相对差些。重复样线法在评估样地物种丰富度时最全面、
最可靠,这可能是因为样线法覆盖了更多的微小生境,且增加样线重复次数,即可较为全面地代表当地物种库。
这与孟凡栋等[19]在西藏高寒草地的研究结果不一致,也可能是区域草地类型和数据处理方法差异的原因。
3.2 β多样性随海拔升高的变化特征
β多样性指沿着环境梯度不同生境群落之间物种组成的相异性或物种沿着环境梯度的更替速率,反映了生
境间的多样性[35]。秦岭太白山群落间的β多样性因坡向而有差异,北坡2800m以下变化不大,2800m以上随
海拔的升高而降低,南坡整体随海拔的升高而降低[36],说明不同植被类型及地形都会影响生境间物种的多样性。
玉龙雪山的木本和灌木的β多样性随海拔的升高呈线性下降趋势,而草本的β多样性随海拔的升高总体呈“单峰
曲线”[24],这可能是木本、灌木和草本植物对山地生境的适应策略差异造成的。青藏高原东部甘肃尕海湖草地群
落β多样性随海拔的升高而降低
[37],这可能是人类活动强烈影响的结果。在本研究中,拉萨河谷山地灌丛草地
随海拔的升高呈“单峰曲线”,说明该地区群落生境间在中海拔处相似性最大(差异最小),物种组成在中海拔处差
异较小,而在低海拔和高海拔处差异较大。这一特征与物种丰富度的海拔分布特征一致[34]。可能原因是中海拔
72第25卷第3期 草业学报2016年
处气候条件相比高海拔处稳定,且人类活动干扰相比低海拔处弱。
3.3 种-面积关系截距犮斜率狕
国外诸多研究[2729]表明种-面积关系中截距犮、斜率狕表征的生态学意义明确,截距犮表征了单位面积的物
种多样性;斜率狕表征随采样面积增大物种积累的速率。乔秀娟等[38]在天山的研究表明草甸和草原的截距犮显
著高于荒漠,草甸和草原的斜率狕显著低于荒漠,说明草甸和草原的α多样性高于荒漠,草甸和草原的β多样性
低于荒漠。本文研究中,种-面积关系截距犮、斜率狕随海拔的升高均呈先增加后降低的趋势,与α、β多样性随海
拔升高的变化特征吻合一致。说明种-面积关系中截距犮、斜率狕也可用来表征拉萨河谷山地灌丛草地生物多
样性的变化特征。
3.4 群落最小面积
种-面积关系法是确定最小样地面积测定中最常用的方法[20]。关于“面积”的取样方法有很多,国内大多采
用常规的巢式取样法[15,39]、改进的巢式样方[40]或以巢式为基础的组合取样法[16,41]。青藏高原针对不同类型[19]
和不同健康状况[20]高寒草地物种多样性的研究亦采用巢式样方。巢式样方是经典种-面积曲线的取样法,适用
性广,其尺度可从群落到区域乃至世界各大洲[42],其取样严格从总体样本抽取小样本,将大样方套小样方,区域
布点均匀,取样方便;但巢式样方只代表了群落的某一局部地段,其确定的最小面积并不足以反映整个地段上的
群落特征。巢式样方在一定程度上受到物种集群分布和环境异质性的影响[30]。而 Modified-Whittaker样地法
特别对空间异质性较大的生态系统具有较高的实用价值,该方法有助于降低空间尺度对物种丰富度的影响,使之
更为有效地量度和预测不同生态系统的物种丰富度[5]。为满足山地灌丛草地研究精度ρ=80%的精确研究,本
研究建议最小面积至少约为100m2。这与董世魁等[20]研究的高寒草甸最小面积为108m2 相近,也符合高山嵩
草草原草甸所采用的100m2 经验值[43]。
4 结论
拉萨河谷山地灌丛草地α多样性随海拔的升高先增加后降低,β多样性随海拔的升高亦呈“单峰曲线”,种-
面积关系的参数截距犮、斜率狕也印证了这种变化特征。说明这是该地区物种多样性呈现的一般规律,是植物对
山地气候因子沿海拔梯度变化的适应特征。未来研究应注重植物多样性分布格局对山体气候梯度以及放牧干扰
的响应和适应机制。
监测方法是影响野外调查的重要手段。本文中利用50m样线3次重复能较好地反映拉萨河谷山地灌丛草
地的α多样性。Modified-Whittaker样地法的结果很好地反映了拉萨河谷山地灌丛草地的β多样性。为满足
研究精度ρ=60%的最低研究,群落最小面积至少约为10m
2;为满足研究精度ρ=80%的精确研究,最小面积至
少约为100m2。我们认为应探究 Modified-Whittaker样地法在青藏高原其他草地类型的实用性和可靠性。
致谢:植物鉴定方面得到西藏大学农牧学院郑维列教授和罗建副教授的大力帮助,在此深表感谢!
犚犲犳犲狉犲狀犮犲狊:
[1] JacquemynH,HonnayO,PailerT.Rangesizevariation,nestednessandspeciesturnoveroforchidspeciesalonganaltitudi
nalgradientonRe’unionIsland:Implicationsforconservation.BiologicalConservation,2007,136:388397.
[2] LegendreP,BorcardD,PeresnetoPR.Analyzingorexplainingbetadiversity?Comment.Ecology,2008,89(11):3238
3244.
[3] ChenSB,OuYangZY,XuW H,犲狋犪犾.Areviewofbetadiversitystudies.BiodiversityScience,2010,18(4):323335.
[4] JostL,DevriesP,WalaT,犲狋犪犾.Partitioningdiversityforconservationanalyses.DiversityandDistributions,2010,16(1):
6576.
[5] LiangJY,YuCT,YangXH.ModifiedWhittaker:amultiscalenestedvegetationsamplingmethod.ResearchofSoiland
WaterConservation,2007,14:226230.
[6] ZhangDY,WangG,ZhaoSL.Studyonthenumberofplantcommunitysuccessionintheabandonedfarmlandofthealpine
meadowinSouthGansuⅠ.Thecharacteristicanalysisofthesuccessionofpioneercommunity.ChineseGrassland,1988,6:
82 ACTAPRATACULTURAESINICA(2016) Vol.25,No.3
1418.
[7] BashkinM,StohlgrenTJ,OtsukiY.SoilcharacteristicsandplantexoticspeciesinvasionsintheGrandStaircase-Escalante
NationalMonument,Utah,USA.AppliedSoilEcology,2002,22:6777.
[8] AdlerPB.NeutralmodelsfailtoreproduceobservedspeciestimeandspeciesarearelationshipsinKansasgrasslands.Ecolo
gy,2004,85:12651272.
[9] WangL,WangDL,BaiYG,犲狋犪犾.Spatialdistributionsofmultipleplantspeciesaffectherbivoreforagingselectivity.Oikos,
2010,119:401408.
[10] StohlgrenTJ,Bul KA,OtsukiY.ComparisonofrangelandvegetationsamplingtechniquesintheCentralGrasslands.
JournalofRangeManage,1998,(51):164172.
[11] WangSP,LiYH,WangYF,犲狋犪犾.Influenceofdifferentstockingratesonplantdiversityof犃狉狋犲犿犻狊犻犪犳狉犻犵犻犱犪community
inInnerMongoliaSteppe.ActaBiotanicaSinica,2001,43(1):8996.
[12] DrakareS,LennonJJ,HilebrandH.Theimprintofgeographical,evolutionaryandecologicalcontextonspeciesarearela
tionships.EcologyLetters,2006,9(2):215227.
[13] TangZY,QiaoXJ,FangJY.Speciesarearelationshipinbiologicalcommunities.BiodiversityScience,2009,17(6):549
559.
[14] LiuCR,MaKP,YuSL,犲狋犪犾.PlantcommunitydiversityinDonglingshanMountain,Beijing,China:Ⅶ.Thedetermina
tionofcriticalsamplingareasforseveraltypesofplantcommunities.ActaEcologySinica,1998,18(1):1523.
[15] WangGP,WangSP.Comparativestudyonthesamplingmethodsof犔犲狔犿狌狊犮犺犻狀犲狀狊犻狊grasslandofInnerMongoliaAuton
omousRegion.ActaAgrestiaSinica,2003,11(4):283288.
[16] ZhangR,ChenJQ,HouYC,犲狋犪犾.Relationshipsbetweenplantspeciesnumberandsamplingareaforsubalpinemeadow
plantcommunities.ChineseJournalofEcology,2013,32(9):22682274.
[17] ZhaoHX.Degradationstatus,causesandimprovingmeasuresofTibetgrassland.Tibet’sScienceandTechnology,2007,
2:4851.
[18] SunHL,ZhengD,YaoTD,犲狋犪犾.ProtectionandconstructionofthenationalecologicalsecurityshelterzoneonTibetan
Plateau.ActaGeographicaSinica,2012,67(1):312.
[19] MengFD,WangCS,ZhangZH,犲狋犪犾.Comparativestudyofmonitoringmethodsaboutplantdiversityandaboveground
biomassofalpinegrasslandsintheTibetanPlateau.ActaBotanicaBorealiOccidentaliaSinica,2013,33(9):19231929.
[20] DongSK,TangL,WangXX,犲狋犪犾.Minimumplotsizeforestimatingplantbiodiversityofthealpinegrasslandsonthe
QinghaiTibetanPlateau.BiodiversityScience,2013,21(6):651657.
[21] ZhangGL,OuYangH,ZhouCP,犲狋犪犾.Responseofagriculturalthermalresourcestoclimatechangeintheregionofthe
brahmaputrariveranditstwotributariesinTibetduringpast50years.ResourcesScience,2010,32(10):19431954.
[22] StohlgrenTJ,FalknerMB,SchelLD.A Modified-Whittakernestedvegetationsamplingmethod.Vegetation,1995,
117:113121.
[23] WuJS,LiXJ,ShenZX,犲狋犪犾.Speciesdiversitydistributionpatternofalpinegrasslandscommunitiesalongaprecipitation
gradientacrossNorthernTibetanPlateau.ActaPrataculturaeSinica,2012,21(3):1725.
[24] FengJM,WangXP,XuCD,犲狋犪犾.AltitudinalpatternsofplantspeciesdiversityandcommunitystructureonYulong
Mountains,Yunnan,China.JournalofMountainScience,2006,24(1):110116.
[25] GaoXM,MaKP,HuangJH,犲狋犪犾.StudiesonplantcommunitydiversityinDonglingshanMountain,Beijing,China:Ⅺ.
Theβdiversityofmountainmeadow.ActaEcologySinica,1998,18(1):2432.
[26] ZhangLM,GaoX,DongK,犲狋犪犾.Studyoncommunityβdiversityquantificationstructuresanditsmeasurements.Journal
ofYunnanAgriculturalUniversity,2014,29(4):578585.
[27] TjrveE.Shapesandfunctionsofspeciesareacurves:areviewofpossiblemodels.JournalofBiogeography,2003,30(6):
827835.
[28] ScheinerSM.Sixtypesofspeciesareacurves.GlobalEcologyandBiogeography,2003,12(6):441447.
[29] ConnorEF,McCoyED.Statisticsandbiologyofthespeciesarearelationship.AmericanNaturalist,1979,113:791833.
[30] GleasonHA.Speciesandarea.Ecology,1925,6:6674.
[31] HustonMA.BiologicalDiversity:TheCoexistenceofSpeciesonChanging[M].Cambridge:CambridgeUniversityPress,
1994:142143.
[32] WangCT,WangQJ,LongRJ,犲狋犪犾.Changesinplantspeciesdiversityandproductivityalonganelevationgradientinan
alpinemeadow.ActaPhytoecologicaSinica,2004,28(2):240245.
[33] DuanMJ,GaoQZ,GuoYQ,犲狋犪犾.Speciesdiversitydistributionpatternofalpinegrasslandcommunitiesalonganaltitudi
nalgradientintheNorthernTibet.PrataculturalScience,2011,28(10):18451850.
92第25卷第3期 草业学报2016年
[34] LuoLM,MiaoYJ,WuJS,犲狋犪犾.Variationinbiodiversityofmontaneshrubgrasslandcommunitiesalonganaltitudinal
gradientinaLhasaRiverbasinvaley.ActaPrataculturaeSinica,2014,23(6):320326.
[35] MaKM,YeW H,SangWG,犲狋犪犾.StudyonplantcommunitydiversityinDonglingshanMoutain,Beijing,China:Ⅹ.β
diversityandfractalanalysisontransectondifferentscales.ActaEcologySinica,1997,7(6):626634.
[36] TangZY,FangJY,ZhangL.PatternsofwoodyplantspeciesdiversityalongenvironmentalgradientsonMt.Taibai,Qin
lingMountains.BiodiversityScience,2004,12(1):115122.
[37] XuPB,DengJM,ZhaoCM,犲狋犪犾.Studyongrasslandcommunitycharacteristicsandspeciesdiversityalongaltitudinal
gradientsintheGahaiwetland,Gansuprovince.ActaPrataculturaeSinica,2012,21(2):219226.
[38] QiaoXJ,TangZY,MmatANWE,犲狋犪犾.Studyonspeciesarearelationshipsinthesouthernandnorthernslopesofthe
TianshanMountains.AridZoneResearch,2011,28(1):5459.
[39] JiangJ,ZhangCY,ZhaoXH.Plantspeciesarearelationshipina42hm2researchplotofconiferousandboardleavedmixed
forestinJiaohe,JilinProvince,China.ChineseJournalofPlantEcology,2012,36(1):3038.
[40] BaiYF,XuZX,LiDX.Studyonαdiversityoffour犛狋犻狆犪communitiesinInnerMongoliaPlateau.BiodiversityScience,
2000,8(4):353360.
[41] ChenH,LiYQ,ZhengSW,犲狋犪犾.Determinationofspeciesarearelationshipsandminimumsamplingareafortheshrub
communitiesinaridvaleyintheupperreachoftheMinjiangRiver.ActaEcologySinica,2007,27(5):18181825.
[42] WeiherE.Thecombinedeffectsofscaleandproductivityonspeciesrichness.JournalofEcology,1999,87:10051011.
[43] WangWY,WangQJ,WangHC.Theeffectoflandmanagementonplantcommunitycomposition,speciesdiversity,and
productivityofalpine犓狅犫犲狉狊犻犪steppemeadow.EcologicalResearch,2006,21(2):181187.
参考文献:
[3] 陈圣宾,欧阳志云,徐卫华,等.Beta多样性研究进展.生物多样性,2010,18(4):323335.
[5] 梁继业,于春堂,杨晓晖.改进的 ModifiedWhittaker:一种多尺度的嵌套式植物多样性取样方法.水土保持研究,2007,14:
226230.
[6] 张大勇,王刚,赵松岭.甘南亚高山草甸弃耕地植物群落演替的数量研究Ⅰ.演替先锋群落的特征分析.中国草地,1988,6:
1418.
[11] 汪诗平,李永宏,王艳芬,等.不同放牧率对内蒙古冷蒿草原植物多样性的影响.植物学报,2001,43(1):8996.
[13] 唐志尧,乔秀娟,方精云.生物群落种-面积关系.生物多样性,2009,17(6):549559.
[14] 刘灿然,马克平,于顺利,等.北京东灵山地区植物群落多样性研究Ⅶ.几种类型植物群落临界抽样面积的确定.生态学
报,1998,18(1):1523.
[15] 王国平,汪诗平.羊草草原植被监测方法的比较研究.草地学报,2003,11(4):283288.
[16] 张蕊,陈军强,侯尧宸,等.亚高山草甸植物物种数与取样面积的关系.生态学杂志,2013,32(9):22682274.
[17] 赵好信.西藏草地退化现状、成因及改良对策.西藏科技,2007,2:4851.
[18] 孙鸿烈,郑度,姚檀栋,等.青藏高原国家生态安全屏障保护与建设.地理学报,2012,67(1):312.
[19] 孟凡栋,王常顺,张振华,等.西藏高原高寒草地群落植物多样性和地上生物量监测方法的比较研究.西北植物学报,
2013,33(9):19231929.
[20] 董世魁,汤琳,王学霞,等.青藏高原高寒草地植物多样性测定的最小样地面积.生物多样性,2013,21(6):651657.
[21] 张戈丽,欧阳华,周才平,等.近50年来气候变化对西藏“一江两河”地区农业气候热量资源的影响.资源科学,2010,
32(10):19431954.
[23] 武建双,李晓佳,沈振西,等.藏北高寒草地样带物种多样性沿降水梯度的分布格局.草业学报,2012,21(3):1725.
[24] 冯建孟,王襄平,徐成东,等.玉龙雪山植物物种多样性和群落结构沿海拔梯度的分布格局.山地学报,2006,24(1):110
116.
[25] 高贤明,马克平,黄建辉,等.北京东灵山地区植物群落多样性的研究Ⅺ.山地草甸U多样性.生态学报,1998,18(1):24
32.
[26] 张立敏,高鑫,董坤,等.生物群落β多样性量化水平及其评价方法.云南农业大学学报,2014,29(4):578585.
[32] 王长庭,王启基,龙瑞军,等.高寒草甸群落植物多样性和初级生产力沿海拔梯度变化的研究.植物生态学报,2004,
28(2):240245.
[33] 段敏杰,高清竹,郭亚奇,等.藏北高寒草地植物群落物种多样性沿海拔梯度的分布格局.草业科学,2011,28(10):1845
1850.
[34] 罗黎鸣,苗彦军,武建双,等.拉萨河谷山地灌丛草地物种多样性随海拔升高的变化特征.草业学报,2014,23(6):320
326.
[35] 马克明,叶万辉,桑卫国,等.北京东灵山地区植物群落多样性研究Ⅹ.不同尺度下群落样带的β多样性及分形分析.生态
03 ACTAPRATACULTURAESINICA(2016) Vol.25,No.3
学报,1997,7(6):626634.
[36] 唐志尧,方精云,张玲.秦岭太白山木本植物物种多样性的梯度格局及环境解释.生物多样性,2004,12(1):115122.
[37] 徐鹏彬,邓建明,赵长明,等.甘肃尕海湿地不同海拔草地群落组分及物种多样性研究.草业学报,2012,21(2):219226.
[38] 乔秀娟,唐志尧,安尼瓦尔·买买提,等.天山南北坡植物种-面积关系.干旱区研究,2011,28(1):5459.
[39] 姜俊,张春雨,赵秀海.吉林蛟河42hm2 针阔混交林样地植物种-面积关系.植物生态学报,2012,36(1):3038.
[40] 白永飞,许志信,李德新.内蒙古高原针茅草原群落α多样性研究.生物多样性,2000,8(4):353360.
[41] 陈泓,黎燕琼,郑绍伟,等.岷江上游干旱河谷灌丛群落种-面积曲线的拟合及最小面积确定.生态学报,2007,27(5):
18181825.
13第25卷第3期 草业学报2016年